Fast assembly of Galerkin matrices for 3D solid laminated composites using finite element and isogeometric discretizations

This work presents a novel methodology for speeding up the assembly of stiffness matrices for laminate composite 3D structures in the context of isogeometric and finite element discretizations. By splitting the involved terms into their in-plane and out-of-plane contributions, this method computes the problems's 3D stiffness matrix as a combination of 2D (in-plane) and 1D (out-of-plane) integrals. Therefore, the assembly's computational complexity is reduced to the one of a 2D problem. Additionally, the number of 2D integrals to be computed becomes independent of the number of material layers that constitute the laminated composite, it only depends on the number of different materials used (or different orientations of the same anisotropic material). Hence, when a high number of layers is present, the proposed technique reduces by orders of magnitude the computational time required to create the stiffness matrix with standard methods, being the resulting matrices identical up to machine precision. The predicted performance is illustrated through numerical experiments.Comment: 40 pages, 13 figure

A priori error for unilateral contact problems with Lagrange multiplier and IsoGeometric Analysis

In this paper, we consider unilateral contact problem without friction between a rigid body and deformable one in the framework of isogeometric analysis. We present the theoretical analysis of the mixed problem using an active-set strategy and for a primal space of NURBS of degree $p$ and $p-2$ for a dual space of B-Spline. A inf-sup stability is proved to ensure a good property of the method. An optimal a priori error estimate is demonstrated without assumption on the unknown contact set. Several numerical examples in two- and three-dimensional and in small and large deformation demonstrate the accuracy of the proposed method

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

Analysis-aware defeaturing of complex geometries with Neumann features

Local modifications of a computational domain are often performed in order to simplify the meshing process and to reduce computational costs and memory requirements. However, removing geometrical features of a domain often introduces a non-negligible error in the solution of a differential problem in which it is defined. In this work, we extend the results from [1] by studying the case of domains containing an arbitrary number of distinct Neumann features, and by performing an analysis on Poisson's, linear elasticity, and Stokes' equations. We introduce a simple, computationally cheap, reliable, and efficient a posteriori estimator of the geometrical defeaturing error. Moreover, we also introduce a geometric refinement strategy that accounts for the defeaturing error: Starting from a fully defeatured geometry, the algorithm determines at each iteration step which features need to be added to the geometrical model to reduce the defeaturing error. These important features are then added to the (partially) defeatured geometrical model at the next iteration, until the solution attains a prescribed accuracy. A wide range of two- and three-dimensional numerical experiments are finally reported to illustrate this work.Comment: 38 page

Fast parametric analysis of trimmed multi-patch isogeometric Kirchhoff-Love shells using a local reduced basis method

This contribution presents a model order reduction framework for real-time efficient solution of trimmed, multi-patch isogeometric Kirchhoff-Love shells. In several scenarios, such as design and shape optimization, multiple simulations need to be performed for a given set of physical or geometrical parameters. This step can be computationally expensive in particular for real world, practical applications. We are interested in geometrical parameters and take advantage of the flexibility of splines in representing complex geometries. In this case, the operators are geometry-dependent and generally depend on the parameters in a non-affine way. Moreover, the solutions obtained from trimmed domains may vary highly with respect to different values of the parameters. Therefore, we employ a local reduced basis method based on clustering techniques and the Discrete Empirical Interpolation Method to construct affine approximations and efficient reduced order models. In addition, we discuss the application of the reduction strategy to parametric shape optimization. Finally, we demonstrate the performance of the proposed framework to parameterized Kirchhoff-Love shells through benchmark tests on trimmed, multi-patch meshes including a complex geometry. The proposed approach is accurate and achieves a significant reduction of the online computational cost in comparison to the standard reduced basis method.Comment: 43 pages, 21 figures, 3 table

Dynamics of high-speed railway bridges: review of design issues and new research for lateral dynamics

Revisión y puesta al día de la publicaciones relacionadas con el diseño de puentes de ferrocarril de alta velocidad y nuevas investigaciones sobre dinámica lateral

Modelos para la interacción dinámica lateral entre viaductos y vehículos ferroviarios de alta velocidad

En este trabajo se proponen modelos num\'ericos para el estudio de la interacción dinámica lateral entre viaductos y vehículos ferroviarios de Alta Velocidad. Los efectos de dinámica vertical, que aparecen en los viaductos al paso de trenes de Alta Velocidad, han sido profusamente estudiados y son bien conocidos hoy en día. Pero, en los últimos diez años, ha aparecido un inter\'es, cada vez mayor, por los fen\'omenos dinámicos laterales que, por lo general, no comprometen la seguridad de la estructura, pero sí la de los vehículos y el confort de los viajeros. Para el estudio de estos efectos laterales se han desarrollado modelos que representan el comportamiento de las estructuras, mediante el m\'etodo de los elementos finitos, y el de los vehículos, mediante sistemas multicuerpo. El acoplamiento entre ambos, pieza fundamental del problema, se realiza mediante elementos de interacción que han sido creados específicamente para este trabajo. Estos elementos tienen en cuenta los movimientos relativos que existen entre los trenes y los puentes, calculan e introducen en la dinámica global del sistema, las fuerzas de contacto que aparecen en las ruedas y los carriles, cuya importancia es crucial en la dinámica lateral de vehículos ferroviarios. Debido a la orografía de España, existen, y se están construyendo, viaductos para las líneas de Alta Velocidad con pilas muy altas, en los que aparecen vientos transversales muy fuertes, y que presentan frecuencias de vibración lateral muy bajas. Los trenes que circulen sobre estas estructuras pueden ser susceptibles de sufrir efectos dinámicos de caracter lateral, y serán el objeto de aplicación de la metodología desarrollada en este trabajo

Coupled models for the dynamics of bridges under high-speed rail traffic

The dynamic effects of high-speed trains on viaducts are important issues for the design of the structures, as well as for determining safe running conditions of trains. In this work we start by reviewing the relevance of some basic moving load models for the dynamic action of vertical traffic loads. The study of lateral dynamics of running trains on bridges is of importance mainly for the safety of the traffic, and may be relevant for laterally compliant bridges. These studies require 3D coupled vehicle-bridge models and consideration of wheel to rail contact. We describe here a fully nonlinear coupled model, formulated in absolute coordinates and incorporated into a commercial finite element framework. An application example is presented for a vehicle subject to a strong wind gust traversing a bridge, showing the relevance of the nonlinear wheel-rail contact model as well as the interaction between bridge and vehicle