5,349 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
Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization
In this paper, we propose a general framework for constructing IGA-suitable
planar B-spline parameterizations from given complex CAD boundaries consisting
of a set of B-spline curves. Instead of forming the computational domain by a
simple boundary, planar domains with high genus and more complex boundary
curves are considered. Firstly, some pre-processing operations including
B\'ezier extraction and subdivision are performed on each boundary curve in
order to generate a high-quality planar parameterization; then a robust planar
domain partition framework is proposed to construct high-quality patch-meshing
results with few singularities from the discrete boundary formed by connecting
the end points of the resulting boundary segments. After the topology
information generation of quadrilateral decomposition, the optimal placement of
interior B\'ezier curves corresponding to the interior edges of the
quadrangulation is constructed by a global optimization method to achieve a
patch-partition with high quality. Finally, after the imposition of
C1=G1-continuity constraints on the interface of neighboring B\'ezier patches
with respect to each quad in the quadrangulation, the high-quality B\'ezier
patch parameterization is obtained by a C1-constrained local optimization
method to achieve uniform and orthogonal iso-parametric structures while
keeping the continuity conditions between patches. The efficiency and
robustness of the proposed method are demonstrated by several examples which
are compared to results obtained by the skeleton-based parameterization
approach
Volumetric Untrimming: Precise decomposition of trimmed trivariates into tensor products
3D objects, modeled using Computer Aided Geometric Design tools, are
traditionally represented using a boundary representation (B-rep), and
typically use spline functions to parameterize these boundary surfaces.
However, recent development in physical analysis, in isogeometric analysis
(IGA) in specific, necessitates a volumetric parametrization of the interior of
the object. IGA is performed directly by integrating over the spline spaces of
the volumetric spline representation of the object. Typically, tensor-product
B-spline trivariates are used to parameterize the volumetric domain. A general
3D object, that can be modeled in contemporary B-rep CAD tools, is typically
represented using trimmed B-spline surfaces. In order to capture the generality
of the contemporary B-rep modeling space, while supporting IGA needs, Massarwi
and Elber (2016) proposed the use of trimmed trivariates volumetric elements.
However, the use of trimmed geometry makes the integration process more
difficult since integration over trimmed B-spline basis functions is a highly
challenging task. In this work, we propose an algorithm that precisely
decomposes a trimmed B-spline trivariate into a set of (singular only on the
boundary) tensor-product B-spline trivariates, that can be utilized to simplify
the integration process in IGA. The trimmed B-spline trivariate is first
subdivided into a set of trimmed B\'ezier trivariates, at all its internal
knots. Then, each trimmed B\'ezier trivariate, is decomposed into a set of
mutually exclusive tensor-product B-spline trivariates, that precisely cover
the entire trimmed domain. This process, denoted untrimming, can be performed
in either the Euclidean space or the parametric space of the trivariate. We
present examples on complex trimmed trivariates' based geometry, and we
demonstrate the effectiveness of the method by applying IGA over the
(untrimmed) results.Comment: 18 pages, 32 figures. Contribution accepted in International
Conference on Geometric Modeling and Processing (GMP 2019
A BEM based on the BĂ©zier/Bernstein polynomial for acoustic waveguide modelization
42nd International Conference on Boundary Elements and other Mesh Reduction Methods, BEM/MRM 2019; ITeCons-University of CoimbraCoimbra; Portugal; 2 July 2019 through 4 July 2019; Code 155806. Publicado en WIT Transactions on Engineering Sciences, Vol 126This paper proposes a novel boundary element approach formulated on the BĂ©zier–Bernstein basis to yield a geometry-independent field approximation. The proposed method is geometrically based on both computer aided design (CAD) and isogeometric analysis (IGA), but field variables are independently approximated from the geometry. This approach allows the appropriate approximation functions for the geometry and variable field to be chosen. We use the BĂ©zier–Bernstein form of a polynomial as an approximation basis to represent both geometry and field variables. The solution of the element interpolation problem in the BĂ©zier–Bernstein space defines generalised Lagrange interpolation functions that are used as element shape functions. The resulting Bernstein–Vandermonde matrix related to the BĂ©zier–Bernstein interpolation problem is inverted using the Newton–Bernstein algorithm. The applicability of the proposed method is demonstrated by solving the Helmholtz equation over an unbounded region in a two-and-a-half dimensional (2.5D) domain.Ministerio de EconomĂa y Competitividad BIA2016-75042-C2-1-
- …