AFT Neural Function Approximators for 1D Nonlinear Force Laws
Authors: Miriam Goldack, Johann Groß, Malte Krack, Merten Stender
Organizations: Technische Universität Berlin, Chair of Cyber-Physical Systems in Mechanical Engineering, Straße des 17. Juni 135, 10623 Berlin, Germany · University of Stuttgart, Institute of Aircraft Propulsion Systems, Pfaffenwaldring 6, 70569 Stuttgart, Germany
Nonlinear contacts and friction strongly influence the vibration response of assembled structures, but their accurate numerical treatment is computationally demanding. The harmonic balance method is widely used to compute periodic steady-state responses, yet the required alternating frequency-time scheme becomes costly for nonsmooth and hysteretic nonlinearities and must be repeated throughout the nonlinear solution process. Here we show that this procedure can be replaced by neural networks that directly map displacement Fourier coefficients to nonlinear force coefficients and provide the corresponding Jacobian through automatic differentiation. The surrounding solver and continuation algorithms remain unchanged for the computation of frequency response curves. The neural networks exclusively learn individual nonlinear elements rather than complete system responses. Physics-based nondimensionalization and phase normalization facilitate the learning process and enable a single trained network to cover a wide range of parameter combinations. Building on the cubic spring, unilateral spring, and Jenkins elements considered here, the approach points toward a reusable library of nonlinear-element surrogates that can be combined in arbitrary number and location within a mechanical system. By bypassing the iterative force evaluation in time domain, the method offers favorable computational scaling for high-resolution analyses and systems with many nonlinear elements.
Figures & tables
Figure 1: Comparison of the AFT- and SFN-based workflows within the HBM. Green and red indicate vectorizable and non-vectorizable operations, respectively; solid and dashed borders denote the SFN and AFT workflows.
Cubic spring
Unilateral spring
Jenkins element
Table 1: Schematics of a single-mass damped oscillator with three different nonlinearities, force-displacement diagrams, and nonlinear force history for different numbers of time samples N .
Cubic spring
Unilateral spring
Jenkins element
Force coefficients: global relative L2 -error εf^,2glob
0.00135
0.13635
0.00314
Force coefficients: component-wise normalized RMSE
Jacobian: mean pointwise relative Frobenius norm error εˉJ,F
0.00058
0.16013
0.04257
Table 2: Comparison of the force-coefficient and Jacobian predictions along the considered frequency-response continuation. All global error metrics are evaluated using the same nΩ points on the solution curve and identical inputs for the AFT and SFN evaluations.
Appendix figures & tables4 assets
Supplementary material from the paper’s appendix.
Appendix
Test case
Samples
Coefficient
Sampling rule
Cubic spring
105
a1′
U[0,5]
a3′
U[−5,5]
b3′
U[−5,5]
Unilateral spring
106
a0∗,′
U[−2,2]
a1∗,′
U[0,2]
a2∗,′
U[−2,2]
Appendix
Table 3: Sampling domains and dataset sizes used for the generation of the SFN training data.
Regime
Samples
Coeff.
Sampling rule
Open/ onset
2500
a0∗,′
0
a1∗,′
N[0.01,1.1](1.0,0.3)
a2∗,′
0
b2∗,′
0
Contact- active
7500
a0∗,′
N[−0.2,0.2](−0.05,0.04)−0.6(a1∗,′−1)
a1∗,′
N[1.0,2.0](1.5,0.3)
Appendix
Table 4: Application-specific sampling strategy for unilateral spring training data.
Cubic spring
Unilateral spring
Jenkins element
Network type
Feedforward neural network
Feedforward neural network
Feedforward neural network
Input dimension
3
4
3
Output dimension
4
5
4
Hidden layers
3
5
5
Neurons per hidden layer
128
128
128
Activation
GELU
GELU
GELU
Appendix
Table 5: Architecture and training strategy of the SFNs used for the three application cases.
Test case
Parameters
Method
ncont
∑nit
nit
∑nFC
Cubic spring
k3=0.1,F0=0.1
AFT
28
54
1.93
82
SFN
28
54
1.93
82
k3=0.1,F0=0.2
AFT
44
85
1.93
129
SFN
44
85
1.93
129
k3=0.1,F0=0.3
AFT
57
106
1.86
163
SFN
57
106
1.86
163
Appendix
Table 6: Solver statistics for the adaptive continuation runs. The AFT-based formulation serves as reference.
We use tools from nonlinear algebra to study the equilibria of small linear translational spring networks. Specifically we use the techniques of homotopy continuation, monodromy, and parameter homotopy (a.k.a. cheater homotopy) to solve all rigid linear translational spring networks up to 5 nodes in both 2 and 3 dimensions. We describe a method of implementing parameter homotopy that arises naturally from the physical structure of the system. We give precise total degree bounds on the maximum number of solutions for general planar spring networks. We discuss further efficiency gains obtained from polyhedral homotopy methods. We compare the computation efficiency of these techniques against a baseline of Newton's method.
Luke Oeding, Ethan Clayton, Jackson Elsea +2
Auburn University, Department of Mathematics and Statistics, Auburn, AL, USA · New York University, Courant Institute of Mathematical Sciences, New York, NY, USA · Air Force Research Laboratory, Eglin Air Force Base, Florida, USA
Extrapolative prediction of complex nonlinear dynamics remains a central challenge in engineering. This study proposes a one-shot learning method to identify global frequency-response curves from a single excitation time history by learning governing equations. We introduce MEv-SINDy (Multi-frequency Evolutionary Sparse Identification of Nonlinear Dynamics) to infer the governing equations of non-autonomous and multi-frequency systems. The methodology leverages the Generalized Harmonic Balance (GHB) method to decompose complex forced responses into a set of slow-varying evolution equations. We validated the capabilities of MEv-SINDy on two critical Micro-Electro-Mechanical Systems (MEMS). These applications include a nonlinear beam resonator and a MEMS micromirror. Our results show that the model trained on a single point accurately predicts softening/hardening effects and jump phenomena across a wide range of excitation levels. This approach significantly reduces the data acquisition burden for the characterization and design of nonlinear microsystems.
This paper proposes a novel data-driven algorithm to approximate the dominant eigenfunctions (aka.~modes) of the Koopman operator of nonlinear dynamical systems using neural networks. The relevance of learning the dominant Koopman modes is to approximate nonlinear dynamics by linear ones in a lifted space, thereby enabling simplified control and analysis. To fight the curse of dimensionality arising from using expressive templates (here neural networks) for the mode approximation, the proposed method leverages a power-iteration scheme that directly learns the dominant Koopman modes without explicitly constructing the projection of the Koopman operator on the template of functions. Our approach connects to other approaches in the literature that avoid the curse of dimensionality by learning small dictionaries of functions, but differs from them in that we do not require ``anti-collapse mechanisms'' to ensure that the learned dictionary is expressive enough to approximate the Koopman operator since our power-iteration scheme is designed to converge toward the dominant modes of the projected Koopman operator. The approach is fully data-driven, requiring only sampled state transitions. Theoretical guarantees are provided, showing convergence under increasing sample size and network width (in connection with the neural tangent kernel theorem). Numerical experiments demonstrate that the method achieves accurate and smooth approximations of dominant modes while avoiding the limitations of traditional techniques such as extended dynamic mode decomposition.