Accelerating Construction of Non-Intrusive Nonlinear Structural Dynamics Reduced Order Models through Hyperreduction

Alexander Saccani,Paolo Tiso
2024-11-22
Abstract:We present a novel technique to significantly reduce the offline cost associated to non-intrusive nonlinear tensors identification in reduced order models (ROMs) of geometrically nonlinear, finite elements (FE)-discretized structural dynamics problems. The ROM is obtained by Galerkin-projection of the governing equations on a reduction basis (RB) of Vibration Modes (VMs) and Static Modal Derivatives (SMDs), resulting in reduced internal forces that are cubic polynomial in the reduced coordinates. The unknown coefficients of the nonlinear tensors associated with this polynomial representation are identified using a modified version of Enhanced Enforced Displacement (EED) method which leverages Energy Conserving Sampling and Weighting (ECSW) as hyperreduction technique for efficiency improvement. Specifically, ECSW is employed to accelerate the evaluations of the nonlinear reduced tangent stiffness matrix that are required within EED. Simulation-free training sets of forces for ECSW are obtained from displacements corresponding to quasi-random samples of a nonlinear second order static displacement manifold. The proposed approach is beneficial for the investigation of the dynamic response of structures subjected to acoustic loading, where multiple VMs must be added in the RB, resulting in expensive nonlinear tensor identification. Superiority of the novel method over standard EED is demonstrated on FE models of a shallow curved clamped panel and of a nine-bay aeronautical reinforced panel modelled, using the commercial finite element program Abaqus.
Numerical Analysis,Computational Physics
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the problem of high off - line computational cost in the nonlinear tensor identification process of non - intrusive reduced - order models (ROMs) in geometrically nonlinear structural dynamics problems. Specifically, the author proposes a novel technique to significantly reduce the off - line computational cost associated with nonlinear tensor identification, especially when dealing with geometrically nonlinear, finite - element - discretized structural dynamics problems. #### Main problems: 1. **High off - line computational cost**: When using the Enhanced Enforced Displacement (EED) method for nonlinear tensor identification, as the size of the reduced - order models (ROMs) increases, the computational cost becomes extremely high. 2. **Need for efficiency improvement**: For scenarios that need to be implemented multiple times (such as geometric optimization or uncertainty quantification in the presence of geometric defects), this high computational cost can become a major problem. 3. **Limitations of existing methods**: Although the traditional EED method can reduce the computational cost, there are still bottlenecks in large - scale ROMs. #### Solutions: The author proposes an improved EED method combined with the Energy Conserving Sampling and Weighting (ECSW) technique to accelerate ROM tensor identification in the following ways: - Use ECSW as a hyper - reduction technique to improve computational efficiency. - Use the nonlinear static displacement manifold generated by quasi - random samples to obtain the training set, thereby avoiding relying on expensive full - fidelity model (HFM) simulations. - Construct a hyper - reduced grid for rapid evaluation of the tangent stiffness matrix for use in the EED scheme. This method can not only significantly reduce the computational cost but also make ROM construction based on commercial finite - element software packages more efficient and practical, suitable for applications in an industrial environment. #### Application background: This method is particularly suitable for studying the dynamic response of aerospace panels under acoustic loads. Such problems usually need to include multiple vibration modes (VMs), resulting in high costs for nonlinear tensor identification. By introducing the new method, the dynamic behavior of these complex structures can be analyzed more effectively, fatigue life can be predicted, and the structural integrity during service can be ensured. #### Formula representation: - The nonlinear internal force vector \(f(q)\) can be represented as a cubic polynomial of the nodal displacement \(q\): \[ f(q)=K^{(1)}\cdot q + K^{(2)}:(q\otimes q)+K^{(3)}:(q\otimes q\otimes q) \] where \(K^{(1)}\in\mathbb{R}^{n\times n}\) is the linear stiffness matrix, and \(K^{(2)}\in\mathbb{R}^{n\times n\times n}\) and \(K^{(3)}\in\mathbb{R}^{n\times n\times n\times n}\) are tensors related to quadratic and cubic forces respectively. - The reduced - order internal force vector \(\tilde{f}(\eta)\) can be represented as: \[ \tilde{f}(\eta)=\tilde{K}^{(1)}\cdot\eta+\tilde{K}^{(2)}:(\eta\otimes\eta)+\tilde{K}^{(3)}:(\eta\otimes\eta\otimes\eta) \] In this way, the author has successfully solved the problem of high computational cost in nonlinear tensor identification and improved the efficiency of ROM construction.