Neural networks for the approximation of Euler's elastica

Elena Celledoni,Ergys Çokaj,Andrea Leone,Sigrid Leyendecker,Davide Murari,Brynjulf Owren,Rodrigo T. Sato Martín de Almagro,Martina Stavole
2024-06-04
Abstract:Euler's elastica is a classical model of flexible slender structures, relevant in many industrial applications. Static equilibrium equations can be derived via a variational principle. The accurate approximation of solutions of this problem can be challenging due to nonlinearity and constraints. We here present two neural network based approaches for the simulation of this Euler's elastica. Starting from a data set of solutions of the discretised static equilibria, we train the neural networks to produce solutions for unseen boundary conditions. We present a $\textit{discrete}$ approach learning discrete solutions from the discrete data. We then consider a $\textit{continuous}$ approach using the same training data set, but learning continuous solutions to the problem. We present numerical evidence that the proposed neural networks can effectively approximate configurations of the planar Euler's elastica for a range of different boundary conditions.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to effectively approximate the static equilibrium configurations of Euler's elastica using neural networks. Specifically, the authors propose two neural - network - based methods for this nonlinear and constrained mechanical problem, aiming to provide accurate approximate solutions for the configurations of Euler's elastica under different boundary conditions. ### Problem Background The Euler's elastica model is a classic model of flexible slender structures and is of great significance in many industrial applications. Its static equilibrium equations can be derived through the variational principle. However, due to the nonlinearity of the problem and the existence of constraint conditions, it is challenging to solve these equations exactly. ### Paper Objectives The main objective of the paper is to develop and evaluate two neural - network methods to simulate the static equilibrium configurations of Euler's elastica: 1. **Discrete Method**: Learn discrete solutions from discrete data sets. 2. **Continuous Method**: Use the same training data set but learn continuous solutions. The core of these two methods is to use the existing data set of discrete static equilibrium solutions to train the neural network so that it can generate solutions for unseen boundary conditions. ### Specific Problem Description - **Discrete Method**: Starting directly from the discretized numerical solutions, train the neural network to predict the displacements and tangent vectors at the internal nodes. - **Continuous Method**: Introduce the arc - length parameterization into the neural network so that the network can output a smooth curve and can ensure the unit - length constraint of the tangent vector through the differential operator. ### Mathematical Model The energy functional of Euler's elastica is defined as: \[ E=\int_{0}^{L} \kappa(s)^{2} \, ds \] where $\kappa(s)$ is the curvature of the curve $q(s)$. To ensure the inextensibility of the curve (i.e., the norm of the tangent vector is 1), the constraint condition is introduced: \[ \Phi(q, q')=\|q'\|^{2}-1 = 0 \] By introducing the Lagrange multiplier $\Lambda(s)$, the augmented Lagrangian function can be constructed: \[ \tilde{L}(q, q', q'', \Lambda)=L(q, q', q'')+\Lambda \Phi(q, q') \] ### Main Contributions The main contributions of the paper include: 1. A detailed analysis of the performance of the discrete neural network in approximating the configurations of Euler's elastica. 2. Pointing out the limitations of the discrete method. 3. Proposing a new parameterization strategy - the continuous neural network - to overcome these limitations. Through these methods, the paper demonstrates the potential of neural networks in dealing with complex mechanical problems, especially in efficiently generating the configurations of elastica under different boundary conditions.