Approximating tensor product Bézier surfaces with tangent plane continuity

AbstractWe present a simple method for degree reduction of tensor product Bézier surfaces with tangent plane continuity in L2-norm. Continuity constraints at the four corners of surfaces are considered, so that the boundary curves preserve endpoints continuity of any order α. We obtain matrix representations for the control points of the degree reduced surfaces by the least-squares method. A simple optimization scheme that minimizes the perturbations of some related control points is proposed, and the surface patches after adjustment are C∞ continuous in the interior and G1 continuous at the common boundaries. We show that this scheme is applicable to surface patches defined on chessboard-like domains

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

The Construction of Optimized High-Order Surface Meshes by Energy-Minimization

Despite the increasing popularity of high-order methods in computational fluid dynamics, their application to practical problems still remains challenging. In order to exploit the advantages of high-order methods with geometrically complex computational domains, coarse curved meshes are necessary, i.e. high-order representations of the geometry. This dissertation presents a strategy for the generation of curved high-order surface meshes. The mesh generation method combines least-squares fitting with energy functionals, which approximate physical bending and stretching energies, in an incremental energy-minimizing fitting strategy. Since the energy weighting is reduced in each increment, the resulting surface representation features high accuracy. Nevertheless, the beneficial influence of the energy-minimization is retained. The presented method aims at enabling the utilization of the superior convergence properties of high-order methods by facilitating the construction of coarser meshes, while ensuring accuracy by allowing an arbitrary choice of geometric approximation order. Results show surface meshes of remarkable quality, even for very coarse meshes representing complex domains, e.g. blood vessels

IST Austria Thesis

Fabrication of curved shells plays an important role in modern design, industry, and science. Among their remarkable properties are, for example, aesthetics of organic shapes, ability to evenly distribute loads, or efficient flow separation. They find applications across vast length scales ranging from sky-scraper architecture to microscopic devices. But, at the same time, the design of curved shells and their manufacturing process pose a variety of challenges. In this thesis, they are addressed from several perspectives. In particular, this thesis presents approaches based on the transformation of initially flat sheets into the target curved surfaces. This involves problems of interactive design of shells with nontrivial mechanical constraints, inverse design of complex structural materials, and data-driven modeling of delicate and time-dependent physical properties. At the same time, two newly-developed self-morphing mechanisms targeting flat-to-curved transformation are presented. In architecture, doubly curved surfaces can be realized as cold bent glass panelizations. Originally flat glass panels are bent into frames and remain stressed. This is a cost-efficient fabrication approach compared to hot bending, when glass panels are shaped plastically. However such constructions are prone to breaking during bending, and it is highly nontrivial to navigate the design space, keeping the panels fabricable and aesthetically pleasing at the same time. We introduce an interactive design system for cold bent glass façades, while previously even offline optimization for such scenarios has not been sufficiently developed. Our method is based on a deep learning approach providing quick and high precision estimation of glass panel shape and stress while handling the shape multimodality. Fabrication of smaller objects of scales below 1 m, can also greatly benefit from shaping originally flat sheets. In this respect, we designed new self-morphing shell mechanisms transforming from an initial flat state to a doubly curved state with high precision and detail. Our so-called CurveUps demonstrate the encodement of the geometric information into the shell. Furthermore, we explored the frontiers of programmable materials and showed how temporal information can additionally be encoded into a flat shell. This allows prescribing deformation sequences for doubly curved surfaces and, thus, facilitates self-collision avoidance enabling complex shapes and functionalities otherwise impossible. Both of these methods include inverse design tools keeping the user in the design loop

Non-Uniform Rational B-Splines and Rational Bezier Triangles for Isogeometric Analysis of Structural Applications

Isogeometric Analysis (IGA) is a major advancement in computational analysis that bridges the gap between a computer-aided design (CAD) model, which is typically constructed using Non-Uniform Rational B-splines (NURBS), and a computational model that traditionally uses Lagrange polynomials to represent the geometry and solution variables. In IGA, the same shape functions that are used in CAD are employed for analysis. The direct manipulation of CAD data eliminates approximation errors that emanate from the process of converting the geometry from CAD to Finite Element Analysis (FEA). As a result, IGA allows the exact geometry to be represented at the coarsest level and maintained throughout the analysis process. While IGA was initially introduced to streamline the design and analysis process, this dissertation shows that IGA can also provide improved computational results for complex and highly nonlinear problems in structural mechanics. This dissertation addresses various problems in structural mechanics in the context of IGA, with the use of NURBS and rational Bézier triangles for the description of the parametric and physical spaces. The approaches considered here show that a number of important properties (e.g., high-order smoothness, geometric exactness, reduced number of degrees of freedom, and increased flexibility in discretization) can be achieved, leading to improved numerical solutions. Specifically, using B-splines and a layer-based discretization, a distributed plasticity isogeometric frame model is formulated to capture the spread of plasticity in large-deformation frames. The modeling approach includes an adaptive analysis where the structure of interest is initially modeled with coarse mesh and knots are inserted based on the yielding information at the quadrature points. It is demonstrated that improvement on efficiency and convergence rates is attained. With NURBS, an isogeometric rotation-free multi-layered plate formulation is developed based on a layerwise deformation theory. The derivation assumes a separate displacement field expansion within each layer, and considers transverse displacement component as C0-continuous at dissimilar material interfaces, which is enforced via knot repetition. The separate integration of the in-plane and through-thickness directions allows to capture the complete 3D stresses in a 2D setting. The proposed method is used to predict the behavior of advanced materials such as laminated composites, and the results show advantages in efficiency and accuracy. To increase the flexibility in discretizing complex geometries, rational Bézier triangles for domain triangulation is studied. They are further coupled with a Delaunay-based feature-preserving discretization algorithm for static bending and free vibration analysis of Kirchhoff plates. Lagrange multipliers are employed to explicitly impose high-order continuity constraints and the augmented system is solved iteratively without increasing the matrix size. The resulting discretization is geometrically exact, admits small geometric features, and constitutes C1-continuity. The feature-preserving rational Bézier triangles are further applied to smeared damage modeling of quasi-brittle materials. Due to the ability of Lagrange multipliers to raise global continuity to any desired order, the implicit fourth- and sixth-order gradient damage models are analyzed. The inclusion of higher-order terms in the nonlocal Taylor expansion improves solution accuracy. A local refinement algorithm that resolves marked regions with high resolution while keeping the resulting mesh conforming and well-conditioned is also utilized to improve efficiency. The outcome is a unified modeling framework where the feature-preserving discretization is able to capture the damage initiation and early-stage propagation, and the local refinement technique can then be applied to adaptively refine the mesh in the direction of damage propagation.PHDCivil EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/147668/1/ningliu_1.pd