Inverse mass matrix for isogeometric explicit transient analysis via the method of localized Lagrange multipliers

Abstract

A variational framework is employed to generate inverse mass matrices for isogeometric analysis (IGA). As different dual bases impact not only accuracy but also computational overhead, several dual bases are extensively investigated.Specifically, locally discontinuous biorthogonal basis functions are evaluated in detail for B-splines of high continuity and Bézier elements with a standard C 0 continuous finite element structure. The boundary conditions are enforced by the method of localized Lagrangian multipliers after generating the inverse mass matrix for completely free body. Thus, unlike inverse mass matrix methods without employing the method of Lagrange multipliers, no modifications in the reciprocal basis functions are needed to account for the boundary conditions. Hence, the present method does not require internal modifications of existing IGA software structures. Numerical examples show that globally continuous dual basis functions yield better accuracy than locally discontinuous biorthogonal functions, but with much higher computational overhead. Locally discontinuous dual basis functions are found to be an economical alternative to lumped mass matrices when combined with mass parameterization. The resulting inverse mass matrices are tested in several vibration problems and applied to explicit transient analysis of structures. KEYWORDSBézier extraction, explicit transient analysis, free vibration, inverse mass matrix, isogeometric analysis, localized Lagrange multipliers Int J Numer Methods Eng. 2019;117:939-966. wileyonlinelibrary.com/journal/nme

Large‐step explicit time integration via mass matrix tailoring

“…Two mass matrix tailoring approaches, namely, the high-frequency filtering method Equation (26) and the low-frequency preservation method Equation (31), have been presented and their performance is assessed. We now offer the following observations:…”

“…To derive the partitioned equations of motion of a linear structural dynamical system, we use the variational formulation proposed by Park and Felippa where the problem is treated like if all bodies were entirely free. Then, the total virtual work of the complete system is obtained by summing up the contributions of each substructure, plus the contribution of the interface constraints via the method of localized Lagrange multipliers (LLM): terms corresponding to the virtual work of the free‐floating substructures and localized interface constraints, respectively.…”

“…where c is a user-specified cutoff frequency beyond which all the higher frequencies will be filtered out. Substituting Equation (26) in the result of Proposition 1, we observe that the modified frequencies take the form:…”

Reciprocal mass matrices and a feasible time step estimator for finite elements with Allman's rotations

“…Algebraic constructions of the reciprocal mass matrix are introduced in papers. [27][28][29][30] These formulations apply to solid finite elements with displacement degrees of freedom and plate finite elements with rotations. The starting point of the formulation is a diagonal mass matrix and a stiffness matrix.…”

Time Step Estimates for Reciprocal Mass Matrices Using Ostrowski's Bounds

“…Thirdly, nodal estimates are so far the best choice for reciprocal mass matrices, where element-based estimates may not be conservative [12]. Several constructions for the reciprocal mass matrices are given in [9,13,6,11,7] and they require efficient time step estimates. In this contribution, the focus is on further development of nodal time step estimates for reciprocal mass matrices and on understanding the factors that limit the sharpness of nodal time step estimate for reciprocal mass matrices.…”