1,042 research outputs found
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
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
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
Algebraic level sets for CAD/CAE integration and moving boundary problems
Boundary representation (B-rep) of CAD models obtained from solid modeling kernels are commonly used in design, and analysis applications outside the CAD systems. Boolean operations between interacting B-rep CAD models as well as analysis of such multi-body systems are fundamental operations on B-rep geometries in CAD/CAE applications. However, the boundary representation of B-rep solids is, in general, not a suitable representation for analysis operations which lead to CAD/CAE integration challenges due to the need for conversion from B-rep to volumetric approximations. The major challenges include intermediate mesh generation step, capturing CAD features and associated behavior exactly and recurring point containment queries for point classification as inside/outside the solid. Thus, an ideal analysis technique for CAD/CAE integration that can enable direct analysis operations on B-rep CAD models while overcoming the associated challenges is desirable. ^ Further, numerical surface intersection operations are typically necessary for boolean operations on B-rep geometries during the CAD and CAE phases. However, for non-linear geometries, surface intersection operations are non-trivial and face the challenge of simultaneously satisfying the three goals of accuracy, efficiency and robustness. In the class of problems involving multi-body interactions, often an implicit knowledge of the boolean operation is sufficient and explicit intersection computation may not be needed. Such implicit boolean operations can be performed by point containment queries on B-rep CAD models. However, for complex non-linear B-rep geometries, the point containment queries may involve numerical iterative point projection operations which are expensive. Thus, there is a need for inexpensive, non-iterative techniques to enable such implicit boolean operations on B-rep geometries. ^ Moreover, in analysis problems with evolving boundaries (ormoving boundary problems), interfaces or cracks, blending functions are used to enrich the underlying domain with the known behavior on the enriching entity. The blending functions are typically dependent on the distance from the evolving boundaries. For boundaries defined by free form curves or surfaces, the distance fields have to be constructed numerically. This may require either a polytope approximation to the boundary and/or an iterative solution to determine the exact distance to the boundary. ^ In this work a purely algebraic, and computationally efficient technique is described for constructing signed distance measures from Non-Uniform Rational B-Splines (NURBS) boundaries that retain the geometric exactness of the boundaries while eliminating the need for iterative and non-robust distance calculation. The proposed technique exploits the NURBS geometry and algebraic tools of implicitization. Such a signed distance measure, also referred to as the Algebraic Level Sets, gives a volumetric representation of the B-rep geometry constructed by purely non-iterative algebraic operations on the geometry. This in turn enables both the implicit boolean operations and analysis operations on B-rep geometries in CAD/CAE applications. Algebraic level sets ensure exactness of geometry while eliminating iterative numerical computations. Further, a geometry-based analysis technique that relies on hierarchical partition of unity field compositions (HPFC) theory and its extension to enriched field modeling is presented. The proposed technique enables direct analysis of complex physical problems without meshing, thus, integrating CAD and CAE. The developed techniques are demonstrated by constructing algebraic level sets for complex geometries, geometry-based analysis of B-rep CAD models and a variety of fracture examples culminating in the analysis of steady state heat conduction in a solid with arbitrary shaped three-dimensional cracks. ^ The proposed techniques are lastly applied to investigate the risk of fracture in the ultra low-k (ULK) dies due to copper (Cu) wirebonding process. Maximum damage induced in the interlayer dielectric (ILD) stack during the process steps is proposed as an indicator of the reliability risk. Numerical techniques based on enriched isogeometric approximations are adopted to model damage in the ULK stacks using a cohesive damage description. A damage analysis procedure is proposed to conduct damage accumulation studies during Cu wirebonding process. Analysis is carried out to identify weak interfaces and potential sites for crack nucleation as well as damage nucleation patterns. Further, the critical process condition is identified by analyzing the damage induced during the impact and ultrasonic excitation stages. Also, representative ILD stack designs with varying Cu percentage are compared for risk of fracture
TiGL - An Open Source Computational Geometry Library for Parametric Aircraft Design
This paper introduces the software TiGL: TiGL is an open source high-fidelity
geometry modeler that is used in the conceptual and preliminary aircraft and
helicopter design phase. It creates full three-dimensional models of aircraft
from their parametric CPACS description. Due to its parametric nature, it is
typically used for aircraft design analysis and optimization. First, we present
the use-case and architecture of TiGL. Then, we discuss it's geometry module,
which is used to generate the B-spline based surfaces of the aircraft. The
backbone of TiGL is its surface generator for curve network interpolation,
based on Gordon surfaces. One major part of this paper explains the
mathematical foundation of Gordon surfaces on B-splines and how we achieve the
required curve network compatibility. Finally, TiGL's aircraft component module
is introduced, which is used to create the external and internal parts of
aircraft, such as wings, flaps, fuselages, engines or structural elements
- …