Optimal analysis-aware parameterization of computational domain in 3D isogeometric analysis

International audienceIn isogeometric analysis framework, computational domain is exactly described using the same representation as that employed in the CAD process. For a CAD object, we can construct various computational domain with same shape but with different parameterization. One basic requirement is that the resulting parameterization should have no self-intersections. In this paper, a linear and easy-to-check sufficient condition for injectivity of trivariate B-spline parameterization is proposed. By an example of 3D thermal conduction problem, we show that different parameterization of computational domain has different impact on the simulation result and efficiency in isogeometric analysis. For problems with exact solutions, we propose a shape optimization method to obtain optimal parameterization of computational domain. The proposed injective condition is used to check the injectivity of initial trivariate B-spline parameterization constructed by discrete Coons volume method, which is the generalization of discrete Coons patch method. Several examples and comparisons are presented to show the effectiveness of the proposed method. Compared with the initial parameterization during refinement, the optimal parameterization can achieve the same accuracy but with less degrees of freedom

Parameterization adaption for 3D shape optimization in aerodynamics

When solving a PDE problem numerically, a certain mesh-refinement process is always implicit, and very classically, mesh adaptivity is a very effective means to accelerate grid convergence. Similarly, when optimizing a shape by means of an explicit geometrical representation, it is natural to seek for an analogous concept of parameterization adaptivity. We propose here an adaptive parameterization for three-dimensional optimum design in aerodynamics by using the so-called "Free-Form Deformation" approach based on 3D tensorial B\'ezier parameterization. The proposed procedure leads to efficient numerical simulations with highly reduced computational costs

The modal analysis of disc resonator gyroscope using isogeometric analysis

This paper presents a free vibration analysis of the disc resonator gyroscope (DRG) based on the Isogeometric Analysis (IGA). The non-uniform rational B-splines (NURBS) basis functions are employed to represent both the exact geometry of DRG and displacement fields. It provides a unified process between the geometric design and analysis which significantly reduces the design cycle by eliminating the re-meshing process. Moreover, the modeling process of IGA and the numerical example of three-circle cascaded model with four different refinement schemes are presented and the computational accuracy and efficiency of the present numerical method are verified by comparing with the standard finite method (ABAQUS)

An integrated design-analysis framework for three dimensional composite panels

We present an integrated design-analysis framework for three dimensional composite panels. The main components of the proposed framework consist of (1) a new curve/surface offset algorithm and (2) the isogeometric concept recently emerged in the computational mechanics community. Using the presented approach, finite element analysis of composite panels can be performed with the only input is the geometry representation of the composite surface. In this paper, non-uniform rational B-splines (NURBS) are used to represent the panel surfaces. A stress analysis of curved composite panel with stiffeners is provided to demonstrate the proposed framework

Analysis-suitable volume parameterization of multi-block computational domain in isogeometric applications

International audienceParameterization of computational domain is a key step in isogeometric analysis just as mesh generation is in finite element analysis. In this paper, we study the volume parameterization problem of multi-block computational domain in isogeometric version, i.e, how to generate analysis-suitable parameterization of the multi-block computational domain bounded by B-spline surfaces. Firstly, we show how to find good volume parameterization of single-block computational domain by solving a constraint optimization problem, in which the constraint condition is the injectivity sufficient conditions of B-spline volume parametrization, and the optimization term is the minimization of quadratic energy functions related to the first and second derivatives of B-spline volume parameterization. By using this method, the resulted volume parameterization has no self-intersections, and the isoparametric structure has good uniformity and orthogonality. Then we extend this method to the multi-block case, in which the continuity condition between the neighbor B-spline volume should be added to the constraint term. The effectiveness of the proposed method is illustrated by several examples based on three-dimensional heat conduction problem

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