11,858 research outputs found
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
Subdivision Shell Elements with Anisotropic Growth
A thin shell finite element approach based on Loop's subdivision surfaces is
proposed, capable of dealing with large deformations and anisotropic growth. To
this end, the Kirchhoff-Love theory of thin shells is derived and extended to
allow for arbitrary in-plane growth. The simplicity and computational
efficiency of the subdivision thin shell elements is outstanding, which is
demonstrated on a few standard loading benchmarks. With this powerful tool at
hand, we demonstrate the broad range of possible applications by numerical
solution of several growth scenarios, ranging from the uniform growth of a
sphere, to boundary instabilities induced by large anisotropic growth. Finally,
it is shown that the problem of a slowly and uniformly growing sheet confined
in a fixed hollow sphere is equivalent to the inverse process where a sheet of
fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless,
quasi-static, elastic limit.Comment: 20 pages, 12 figures, 1 tabl
Approximation of tensor fields on surfaces of arbitrary topology based on local Monge parametrizations
We introduce a new method, the Local Monge Parametrizations (LMP) method, to
approximate tensor fields on general surfaces given by a collection of local
parametrizations, e.g.~as in finite element or NURBS surface representations.
Our goal is to use this method to solve numerically tensor-valued partial
differential equations (PDE) on surfaces. Previous methods use scalar
potentials to numerically describe vector fields on surfaces, at the expense of
requiring higher-order derivatives of the approximated fields and limited to
simply connected surfaces, or represent tangential tensor fields as tensor
fields in 3D subjected to constraints, thus increasing the essential number of
degrees of freedom. In contrast, the LMP method uses an optimal number of
degrees of freedom to represent a tensor, is general with regards to the
topology of the surface, and does not increase the order of the PDEs governing
the tensor fields. The main idea is to construct maps between the element
parametrizations and a local Monge parametrization around each node. We test
the LMP method by approximating in a least-squares sense different vector and
tensor fields on simply connected and genus-1 surfaces. Furthermore, we apply
the LMP method to two physical models on surfaces, involving a tension-driven
flow (vector-valued PDE) and nematic ordering (tensor-valued PDE). The LMP
method thus solves the long-standing problem of the interpolation of tensors on
general surfaces with an optimal number of degrees of freedom.Comment: 16 pages, 6 figure
H-matrix accelerated second moment analysis for potentials with rough correlation
We consider the efficient solution of partial differential equationsfor strongly elliptic operators with constant coefficients and stochastic Dirichlet data by the boundary integral equation method. The computation of the solution's two-point correlation is well understood if the two-point correlation of the Dirichlet data is known and sufficiently smooth.Unfortunately, the problem becomes much more involved in case of rough data. We will show that the concept of the H-matrix arithmetic provides a powerful tool to cope with this problem. By employing a parametric surface representation, we end up with an H-matrix arithmetic based on balanced cluster trees. This considerably simplifies the implementation and improves the performance of the H-matrix arithmetic. Numerical experiments are provided to validate and quantify the presented methods and algorithms
A global approach to the refinement of manifold data
A refinement of manifold data is a computational process, which produces a
denser set of discrete data from a given one. Such refinements are closely
related to multiresolution representations of manifold data by pyramid
transforms, and approximation of manifold-valued functions by repeated
refinements schemes. Most refinement methods compute each refined element
separately, independently of the computations of the other elements. Here we
propose a global method which computes all the refined elements simultaneously,
using geodesic averages. We analyse repeated refinements schemes based on this
global approach, and derive conditions guaranteeing strong convergence.Comment: arXiv admin note: text overlap with arXiv:1407.836
Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures
The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly
transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit
parametrization of boundary surfaces to impose traction and displacement boundary conditions, which
constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing
traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model
generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as
initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer
all boundary terms in the variational formulation into volumetric terms. We show that in the context of the
voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions
defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field,
the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method
by analyzing stresses in a human femur and a vertebral body
Multiscale Representations for Manifold-Valued Data
We describe multiscale representations for data observed on equispaced grids and taking values in manifolds such as the sphere , the special orthogonal group , the positive definite matrices , and the Grassmann manifolds . The representations are based on the deployment of Deslauriers--Dubuc and average-interpolating pyramids "in the tangent plane" of such manifolds, using the and maps of those manifolds. The representations provide "wavelet coefficients" which can be thresholded, quantized, and scaled in much the same way as traditional wavelet coefficients. Tasks such as compression, noise removal, contrast enhancement, and stochastic simulation are facilitated by this representation. The approach applies to general manifolds but is particularly suited to the manifolds we consider, i.e., Riemannian symmetric spaces, such as , , , where the and maps are effectively computable. Applications to manifold-valued data sources of a geometric nature (motion, orientation, diffusion) seem particularly immediate. A software toolbox, SymmLab, can reproduce the results discussed in this paper
- …