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 New Approach in the Design of Microstructured Ultralight Components to Achieve Maximum Functional Performance

In the energy and aeronautics industry, some components need to be very light but with high strength. For instance, turbine blades and structural components under rotational centrifugal forces, or internal supports, ask for low weight, and in general, all pieces in energy turbine devices will benefit from weight reductions. In space applications, a high ratio strength/weight is even more important. Light components imply new optimal design concepts, but to be able to be manufactured is the real key enable technology. Additive manufacturing can be an alternative, applying radical new approaches regarding part design and components’ internal structure. Here, a new approach is proposed using the replica of a small structure (cell) in two or three orders of magnitude. Laser Powder Bed Fusion (L-PBF) is one of the most well-known additive manufacturing methods of functional parts (and prototypes as well), for instance, starting from metal powders of heat-resistant alloys. The working conditions for such components demand high mechanical properties at high temperatures, Ni-Co superalloys are a choice. The work here presented proposes the use of “replicative” structures in different sizes and orders of magnitude, to manufacture parts with the minimum weight but achieving the required mechanical properties. Printing process parameters and mechanical performance are analyzed, along with several examples.Thanks are owed to H2020-FETOPEN-2018-2019-2020-01 ADAM2 PROJECT Analysis, Design, And Manufacturing using Microstructures and Authors are grateful to Basque government group IT IT1337-19 and the Ministry of Mineco REF DPI2016-74845-R, PID2019-109340RB-I00, KK-2020/00102, KK-2020/00042 and PID2019-104488RB-I00

IGA-based Multi-Index Stochastic Collocation for random PDEs on arbitrary domains

This paper proposes an extension of the Multi-Index Stochastic Collocation (MISC) method for forward uncertainty quantification (UQ) problems in computational domains of shape other than a square or cube, by exploiting isogeometric analysis (IGA) techniques. Introducing IGA solvers to the MISC algorithm is very natural since they are tensor-based PDE solvers, which are precisely what is required by the MISC machinery. Moreover, the combination-technique formulation of MISC allows the straight-forward reuse of existing implementations of IGA solvers. We present numerical results to showcase the effectiveness of the proposed approach.Comment: version 3, version after revisio

A minimal stabilization procedure for Isogeometric methods on trimmed geometries

Trimming is a common operation in CAD, and, in its simplest formulation, consists in removing superfluous parts from a geometric entity described via splines (a spline patch). After trimming the geometric description of the patch remains unchanged, but the underlying mesh is unfitted with the physical object. We discuss the main problems arising when solving elliptic PDEs on a trimmed domain. First we prove that, even when Dirichlet boundary conditions are weakly enforced using Nitsche's method, the resulting method suffers lack of stability. Then, we develop novel stabilization techniques based on a modification of the variational formulation, which allow us to recover well-posedness and guarantee accuracy. Optimal a priori error estimates are proven, and numerical examples confirming the theoretical results are provided

Heterogeneous Parametric Trivariate Fillets

Blending and filleting are well established operations in solid modeling and computer-aided geometric design. The creation of a transition surface which smoothly connects the boundary surfaces of two (or more) objects has been extensively investigated. In this work, we introduce several algorithms for the construction of, possibly heterogeneous, trivariate fillets, that support smooth filleting operations between pairs of, possibly heterogeneous, input trivariates. Several construction methods are introduced that employ functional composition algorithms as well as introduce a half Volumetric Boolean sum operation. A volumetric fillet, consisting of one or more tensor product trivariate(s), is fitted to the boundary surfaces of the input. The result smoothly blends between the two inputs, both geometrically and material-wise (properties of arbitrary dimension). The application of encoding heterogeneous material information into the constructed fillet is discussed and examples of all proposed algorithms are presented

Interpolatory Catmull-Clark volumetric subdivision over unstructured hexahedral meshes for modeling and simulation applications

International audienceVolumetric modeling is an important topic for material modeling and isogeometric simulation. In this paper, two kinds of interpolatory Catmull-Clark volumetric subdivision approaches over unstructured hexahedral meshes are proposed based on the limit point formula of Catmull-Clark subdivision volume. The basic idea of the first method is to construct a new control lattice, whose limit volume by the CatmullClark subdivision scheme interpolates vertices of the original hexahedral mesh. The new control lattice is derived by the local push-back operation from one CatmullClark subdivision step with modified geometric rules. This interpolating method is simple and efficient, and several shape parameters are involved in adjusting the shape of the limit volume. The second method is based on progressive-iterative approximation using limit point formula. At each iteration step, we progressively modify vertices of an original hexahedral mesh to generate a new control lattice whose limit volume interpolates all vertices in the original hexahedral mesh. The convergence proof of the iterative process is also given. The interpolatory subdivision volume has C 2-smoothness at the regular region except around extraordinary vertices and edges. Furthermore, the proposed interpolatory volumetric subdivision methods can be used not only for geometry interpolation, but also for material attribute interpolation in the field of volumetric material modeling. The application of the proposed volumetric subdivision approaches on isogeometric analysis is also given with several examples