SUMMARYA rigorous computational framework for the dimensional reduction of discrete, high-fidelity, nonlinear, finite element structural dynamics models is presented. It is based on the pre-computation of solution snapshots, their compression into a reduced-order basis, and the Galerkin projection of the given discrete high-dimensional model onto this basis. To this effect, this framework distinguishes between vector-valued displacements and manifold-valued finite rotations. To minimize computational complexity, it also differentiates between the cases of constant and configuration-dependent mass matrices. Like most projection-based nonlinear model reduction methods, however, its computational efficiency hinges not only on the ability of the constructed reduced-order basis to capture the dominant features of the solution of interest but also on the ability of this framework to compute fast and accurate approximations of the projection onto a subspace of tangent matrices and/or force vectors. The computation of the latter approximations is often referred to in the literature as hyper reduction. Hence, this paper also presents the energy-conserving sampling and weighting (ECSW) hyper reduction method for discrete (or semi-discrete), nonlinear, finite element structural dynamics models. Based on mesh sampling and the principle of virtual work, ECSW is natural for finite element computations and preserves an important energetic aspect of the high-dimensional finite element model to be reduced. Equipped with this hyper reduction procedure, the aforementioned Galerkin projection framework is first demonstrated for several academic but challenging problems. Then, its potential for the effective solution of real problems is highlighted with the realistic simulation of the transient response of a vehicle to an underbody blast event. For this problem, the proposed nonlinear model reduction framework reduces the CPU time required by a typical high-dimensional model by up to four orders of magnitude while maintaining a good level of accuracy.
The computational efficiency of a typical, projection-based, nonlinear model reduction method hinges on the efficient approximation, for explicit computations, of the scalar projections onto a subspace of a residual vector. For implicit computations, it also hinges on the additional efficient approximation of similar projections of the Jacobian of this residual with respect to the solution. The computation of both approximations is often referred to in the literature as hyper reduction. To this effect, this paper focuses on the analysis and comparative performance study for nonlinear finite element reduced-order models of solids and structures of the recently developed energy-conserving mesh sampling and weighting (ECSW) hyper reduction method. Unlike most alternative approaches, this method approximates the scalar projections of residuals and/or Jacobians directly, instead of approximating first these vectors and matrices then projecting the resulting approximations onto the subspaces of interest. In this paper, it is shown that ECSW distinguishes itself furthermore from other hyper reduction methods through its preservation of the Lagrangian structure associated with Hamilton's principle. For second-order dynamical systems, this enables it to also preserve the numerical stability properties of the discrete system to which it is applied. It is also shown that for a fixed set of parameter values, the approximation error committed online by ECSW is bounded by its counterpart error committed off-line during the training of this method. Therefore, this error can be estimated in this case a priori and controlled. The performance of ECSW is demonstrated first for two academic but nevertheless interesting nonlinear dynamic response problems. For both of them, ECSW is shown to preserve numerical stability and deliver the desired level of accuracy, whereas the discrete empirical interpolation method and its recently introduced unassembled variant are shown to be susceptible to failure because of numerical instability. The potential of ECSW for enabling the effective reduction of nonlinear finite element dynamic models of solids and structures is also highlighted with the realistic simulation of the nonlinear transient dynamic response of a complete car engine to thermal and combustion pressure loads using an implicit scheme. For this simulation, ECSW is reported to enable the reduction of the CPU time required by the high-dimensional nonlinear finite element dynamic analysis by more than four orders of magnitude, while achieving a very good level of accuracy.
A dual-primal variant of the FETI-H domain decomposition method is designed for the fast, parallel, iterative solution of large-scale systems of complex equations arising from the discretization of acoustic scattering problems formulated in bounded computational domains. The convergence of this iterative solution method, named here FETI-DPH, is shown to scale with the problem size, the number of subdomains, and the wave number. Its solution time is also shown to scale with the problem size. CPU performance results obtained for the acoustic signature analysis in the mid-frequency regime of mockup submarines reveal that the proposed FETI-DPH solver is significantly faster than the previous generation FETI-H solution algorithm.
SUMMARYProblems of the form Z( )u( ) = f( ), where Z is a given matrix, f is a given vector, and is a circular frequency or circular frequency-related parameter arise in many applications including computational structural and fluid dynamics, and computational acoustics and electromagnetics. The straightforward solution of such problems for fine increments of is computationally prohibitive, particularly when Z is a large-scale matrix. This paper discusses an alternative solution approach based on the efficient computation of u and its successive derivatives with respect to at a few sample values of this parameter, and the reconstruction of the solution u( ) in the frequency band of interest using multi-point Padé approximants. This computational methodology is illustrated with applications from structural dynamics and underwater acoustic scattering. In each case, it is shown to reduce the CPU time required by the straightforward approach to frequency sweep computations by two orders of magnitude.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.