Combining Set Propagation with Finite Element Methods for Time Integration in Transient Solid Mechanics Problems

The Finite Element Method (FEM) is the gold standard for spatial discretization in numerical simulations for a wide spectrum of real-world engineering problems. Prototypical areas of interest include linear heat transfer and linear structural dynamics problems modeled with partial differential equations (PDEs). While different algorithms for direct integration of the equations of motion exist, exploring all feasible behaviors for varying loads, initial states and fluxes in models with large numbers of degrees of freedom remains a challenging task. In this article we propose a novel approach, based in set propagation methods and motivated by recent advances in the field of Reachability Analysis. Assuming a set of initial states and inputs, the proposed method consists in the construction of a union of sets (flowpipe) that enclose the infinite number of solutions of the spatially discretized PDE. We present the numerical results obtained in five examples to illustrate the capabilities of our approach, and compare its performance against reference numerical integration methods. We conclude that, for problems with single known initial conditions, the proposed method is accurate. For problems with uncertain initial conditions included in sets, the proposed method can compute all the solutions of the system more efficiently than numerical integration methods.

[1]  A. Girard,et al.  Reachability analysis of linear systems using support functions , 2010 .

[2]  K. Bathe,et al.  The Bathe time integration method revisited for prescribing desired numerical dissipation , 2019, Computers & Structures.

[3]  Xin Chen,et al.  Flow*: An Analyzer for Non-linear Hybrid Systems , 2013, CAV.

[4]  Matthias Althoff,et al.  An Introduction to CORA 2015 , 2015, ARCH@CPSWeek.

[5]  Wooram Kim,et al.  An improved implicit time integration algorithm: The generalized composite time integration algorithm , 2018 .

[6]  David Moens,et al.  Robust uncertainty quantification in structural dynamics under scarse experimental modal data: A Bayesian-interval approach , 2020 .

[7]  Eric Li,et al.  Large time steps in the explicit formulation of transient heat transfer , 2017 .

[8]  R. Fariña,et al.  Finite element and morphological analysis in extant mammals’ claws and quaternary sloths’ ungual phalanges , 2019, Historical Biology.

[9]  Pierre Ladevèze,et al.  Error Estimate Procedure in the Finite Element Method and Applications , 1983 .

[10]  Alan Edelman,et al.  Julia: A Fresh Approach to Numerical Computing , 2014, SIAM Rev..

[11]  P. Forooghi,et al.  Quantifying structural uncertainties in Reynolds-averaged Navier–Stokes simulations of wind turbine wakes , 2021 .

[12]  Sergiy Bogomolov,et al.  Reach Set Approximation through Decomposition with Low-dimensional Sets and High-dimensional Matrices , 2018, HSCC.

[13]  Benjamin J. Ellis,et al.  FEBio: finite elements for biomechanics. , 2012, Journal of biomechanical engineering.

[14]  Guijun Wang,et al.  The Mass Concrete Temperature Simulation Analysis with Fuzzy Factors , 2010, Annual Conference of Fuzzy Information and Engineering.

[15]  A. Al-Jumaily,et al.  Prediction of upper airway dryness and optimal continuous positive airway pressure conditions , 2020, Journal of Biomechanics.

[16]  Pablo J. Blanco,et al.  Multidimensional modelling for the carotid artery blood flow , 2006 .

[17]  Federico Perotti,et al.  Numerical analysis of the non-linear dynamic behaviour of suspended cables under turbulent wind excitation , 2001 .

[18]  Jean-Marc Battini Analysis of Dampers for Stay Cables Using Non Linear Beam Elements , 2018 .

[19]  JuliaReach , 2019, Proceedings of the 22nd ACM International Conference on Hybrid Systems: Computation and Control.

[20]  K. Bathe,et al.  Insight into an implicit time integration scheme for structural dynamics , 2012 .

[21]  Edward L. Wilson,et al.  Finite element analysis of linear and nonlinear heat transfer , 1974 .

[22]  M. Géradin,et al.  Mechanical Vibrations: Theory and Application to Structural Dynamics , 1994 .

[23]  T. Henzinger,et al.  Algorithmic analysis of nonlinear hybrid systems , 1998, IEEE Trans. Autom. Control..

[24]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[25]  Taylor T. Johnson,et al.  Numerical Verification of Affine Systems with up to a Billion Dimensions , 2019 .

[26]  Qing Nie,et al.  DifferentialEquations.jl – A Performant and Feature-Rich Ecosystem for Solving Differential Equations in Julia , 2017, Journal of Open Research Software.

[27]  Matthias Althoff,et al.  Combining zonotopes and support functions for efficient reachability analysis of linear systems , 2016, 2016 IEEE 55th Conference on Decision and Control (CDC).

[28]  Sergiy Bogomolov,et al.  JuliaReach: a toolbox for set-based reachability , 2019, HSCC.

[29]  R. Barretta,et al.  Variational Formulation of the First Principle of Continuum Thermodynamics , 2010 .

[30]  Eugenio Oñate Ibáñez de Navarra,et al.  Structural analysis with the Finite Element Method. , 2009 .

[31]  Antoine Girard,et al.  Set Propagation Techniques for Reachability Analysis , 2021, Annu. Rev. Control. Robotics Auton. Syst..

[32]  Klaus-Jürgen Bathe,et al.  Accurate solution of wave propagation problems in elasticity , 2021 .

[33]  C. C. Law,et al.  ParaView: An End-User Tool for Large-Data Visualization , 2005, The Visualization Handbook.

[34]  Indra Narayan Kar,et al.  Wind Turbine Structural Modeling Consideration for Dynamic Studies of DFIG Based System , 2017, IEEE Transactions on Sustainable Energy.

[35]  Masanobu Shinozuka,et al.  Monte Carlo solution of structural dynamics , 1972 .

[36]  L. G. Jaeger,et al.  Dynamics of structures , 1990 .

[37]  Matthias Althoff,et al.  Adaptive Parameter Tuning for Reachability Analysis of Linear Systems , 2020, 2020 59th IEEE Conference on Decision and Control (CDC).

[38]  Lu Zhen-zhou,et al.  Interval arithmetic and static interval finite element method , 2001 .

[39]  Mohammad Tahersima,et al.  Finite element modeling of hydration heat in a concrete slab-on-grade floor with limestone blended cement , 2017 .

[40]  Robert L. Mullen,et al.  Formulation of Fuzzy Finite‐Element Methods for Solid Mechanics Problems , 1999 .

[41]  Christian Soize,et al.  Uncertainty quantification in computational linear structural dynamics for viscoelastic composite structures , 2016 .

[42]  K. Bathe Finite Element Procedures , 1995 .

[43]  Antoine Girard,et al.  SpaceEx: Scalable Verification of Hybrid Systems , 2011, CAV.

[44]  Matthias Althoff,et al.  ARCH-COMP20 Category Report: Continuous and Hybrid Systems with Linear Continuous Dynamics , 2020 .

[45]  Long Chen FINITE ELEMENT METHOD , 2013 .

[46]  Matthias Althoff,et al.  Reachability analysis of nonlinear systems with uncertain parameters using conservative linearization , 2008, 2008 47th IEEE Conference on Decision and Control.

[47]  Antoine Girard,et al.  Reachability Analysis of Nonlinear Systems Using Conservative Approximation , 2003, HSCC.

[48]  A. Bayen,et al.  Assessment of Uncertainty Propagation in the Dynamic Response of Single-Degree-of-Freedom Structures Using Reachability Analysis , 2014 .

[50]  Robert L. Mullen,et al.  Interval Based Finite Elements for Uncertainty Quantification in Engineering Mechanics , 2011, WoCoUQ.

[51]  O. C. Zienkiewicz,et al.  Introductory Lectures on the Finite Element Method , 1967 .

[52]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[53]  Antoine Girard,et al.  Efficient Computation of Reachable Sets of Linear Time-Invariant Systems with Inputs , 2006, HSCC.

[54]  Matthias Althoff,et al.  Reachability Analysis of Large Linear Systems With Uncertain Inputs in the Krylov Subspace , 2017, IEEE Transactions on Automatic Control.

[55]  Colas Le Guernic Reachability Analysis of Hybrid Systems with Linear Continuous Dynamics. (Calcul d'Atteignabilité des Systèmes Hybrides à Partie Continue Linéaire) , 2009 .

[56]  Davide Bresolin,et al.  Reachability computation for hybrid systems with Ariadne , 2008 .

[57]  H. Saunders,et al.  Finite element procedures in engineering analysis , 1982 .