206 research outputs found
Subdivision surfaces with isogeometric analysis adapted refinement weights
Subdivision surfaces provide an elegant isogeometric analysis framework for
geometric design and analysis of partial differential equations defined on
surfaces. They are already a standard in high-end computer animation and
graphics and are becoming available in a number of geometric modelling systems
for engineering design. The subdivision refinement rules are usually adapted
from knot insertion rules for splines. The quadrilateral Catmull-Clark scheme
considered in this work is equivalent to cubic B-splines away from
extraordinary, or irregular, vertices with other than four adjacent elements.
Around extraordinary vertices the surface consists of a nested sequence of
smooth spline patches which join continuously at the point itself. As
known from geometric design literature, the subdivision weights can be
optimised so that the surface quality is improved by minimising
short-wavelength surface oscillations around extraordinary vertices. We use the
related techniques to determine weights that minimise finite element
discretisation errors as measured in the thin-shell energy norm. The
optimisation problem is formulated over a characteristic domain and the errors
in approximating cup- and saddle-like quadratic shapes obtained from
eigenanalysis of the subdivision matrix are minimised. In finite element
analysis the optimised subdivision weights for either cup- or saddle-like
shapes are chosen depending on the shape of the solution field around an
extraordinary vertex. As our computations confirm, the optimised subdivision
weights yield a reduction of and more in discretisation errors in the
energy and norms. Although, as to be expected, the convergence rates are
the same as for the classical Catmull-Clark weights, the convergence constants
are improved.Partial support through Trimble Inc and Cambridge Trust is gratefully acknowledged
Isogeometric FEM-BEM coupled structural-acoustic analysis of shells using subdivision surfaces
We introduce a coupled finite and boundary element formulation for acoustic
scattering analysis over thin shell structures. A triangular Loop subdivision
surface discretisation is used for both geometry and analysis fields. The
Kirchhoff-Love shell equation is discretised with the finite element method and
the Helmholtz equation for the acoustic field with the boundary element method.
The use of the boundary element formulation allows the elegant handling of
infinite domains and precludes the need for volumetric meshing. In the present
work the subdivision control meshes for the shell displacements and the
acoustic pressures have the same resolution. The corresponding smooth
subdivision basis functions have the continuity property required for the
Kirchhoff-Love formulation and are highly efficient for the acoustic field
computations. We validate the proposed isogeometric formulation through a
closed-form solution of acoustic scattering over a thin shell sphere.
Furthermore, we demonstrate the ability of the proposed approach to handle
complex geometries with arbitrary topology that provides an integrated
isogeometric design and analysis workflow for coupled structural-acoustic
analysis of shells
Fast Isogeometric Boundary Element Method based on Independent Field Approximation
An isogeometric boundary element method for problems in elasticity is
presented, which is based on an independent approximation for the geometry,
traction and displacement field. This enables a flexible choice of refinement
strategies, permits an efficient evaluation of geometry related information, a
mixed collocation scheme which deals with discontinuous tractions along
non-smooth boundaries and a significant reduction of the right hand side of the
system of equations for common boundary conditions. All these benefits are
achieved without any loss of accuracy compared to conventional isogeometric
formulations. The system matrices are approximated by means of hierarchical
matrices to reduce the computational complexity for large scale analysis. For
the required geometrical bisection of the domain, a strategy for the evaluation
of bounding boxes containing the supports of NURBS basis functions is
presented. The versatility and accuracy of the proposed methodology is
demonstrated by convergence studies showing optimal rates and real world
examples in two and three dimensions.Comment: 32 pages, 27 figure
Isogeometric Analysis on V-reps: first results
Inspired by the introduction of Volumetric Modeling via volumetric
representations (V-reps) by Massarwi and Elber in 2016, in this paper we
present a novel approach for the construction of isogeometric numerical methods
for elliptic PDEs on trimmed geometries, seen as a special class of more
general V-reps. We develop tools for approximation and local re-parametrization
of trimmed elements for three dimensional problems, and we provide a
theoretical framework that fully justify our algorithmic choices. We validate
our approach both on two and three dimensional problems, for diffusion and
linear elasticity.Comment: 36 pages, 44 figures. Reviewed versio
Isogeometric FEM-BEM coupled structural-acoustic analysis of shells using subdivision surfaces
We introduce a coupled finite and boundary element formulation for acoustic
scattering analysis over thin shell structures. A triangular Loop subdivision
surface discretisation is used for both geometry and analysis fields. The
Kirchhoff-Love shell equation is discretised with the finite element method and
the Helmholtz equation for the acoustic field with the boundary element method.
The use of the boundary element formulation allows the elegant handling of
infinite domains and precludes the need for volumetric meshing. In the present
work the subdivision control meshes for the shell displacements and the
acoustic pressures have the same resolution. The corresponding smooth
subdivision basis functions have the continuity property required for the
Kirchhoff-Love formulation and are highly efficient for the acoustic field
computations. We validate the proposed isogeometric formulation through a
closed-form solution of acoustic scattering over a thin shell sphere.
Furthermore, we demonstrate the ability of the proposed approach to handle
complex geometries with arbitrary topology that provides an integrated
isogeometric design and analysis workflow for coupled structural-acoustic
analysis of shells
Shape optimisation with multiresolution subdivision surfaces and immersed finite elements
We develop a new optimisation technique that combines multiresolution
subdivision surfaces for boundary description with immersed finite elements for
the discretisation of the primal and adjoint problems of optimisation. Similar
to wavelets multiresolution surfaces represent the domain boundary using a
coarse control mesh and a sequence of detail vectors. Based on the
multiresolution decomposition efficient and fast algorithms are available for
reconstructing control meshes of varying fineness. During shape optimisation
the vertex coordinates of control meshes are updated using the computed shape
gradient information. By virtue of the multiresolution editing semantics,
updating the coarse control mesh vertex coordinates leads to large-scale
geometry changes and, conversely, updating the fine control mesh coordinates
leads to small-scale geometry changes. In our computations we start by
optimising the coarsest control mesh and refine it each time the cost function
reaches a minimum. This approach effectively prevents the appearance of
non-physical boundary geometry oscillations and control mesh pathologies, like
inverted elements. Independent of the fineness of the control mesh used for
optimisation, on the immersed finite element grid the domain boundary is always
represented with a relatively fine control mesh of fixed resolution. With the
immersed finite element method there is no need to maintain an analysis
suitable domain mesh. In some of the presented two- and three-dimensional
elasticity examples the topology derivative is used for creating new holes
inside the domain.The partial support of the EPSRC through grant # EP/G008531/1 and EC through Marie Curie Actions (IAPP) program CASOPT project are gratefully acknowledged.This is the final version of the article. It was first available from Elsevier via http://dx.doi.org/10.1016/j.cma.2015.11.01
- …