The LATIN-PGD methodology to nonlinear dynamics and quasi-brittle materials for future earthquake engineering applications

Sebastian Rodriguez,Pierre-Etienne Charbonnel,Pierre Ladevèze,David Néron

2024-08-09

Abstract:This paper presents a first implementation of the LArge Time INcrement (LATIN) method along with the model reduction technique called Proper Generalized Decomposition (PGD) for solving nonlinear low-frequency dynamics problems when dealing with a quasi-brittle isotropic damage constitutive relations. The present paper uses the Time-Discontinuous Galerkin Method (TDGM) for computing the temporal contributions of the space-time separate-variables solution of the LATIN-PGD approach, which offers several advantages when considering a high number of DOFs in time. The efficiency of the method is tested for the case of a 3D bending beam, where results and benchmarks comparing LATIN-PGD to classical time-incremental Newmark/Quasi-Newton nonlinear solver are presented. This work represents a first step towards taking into account uncertainties and carrying out more complex parametric studies imposed by seismic risk assessment.

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in earthquake engineering applications, to develop an efficient numerical calculation method for non - linear low - frequency dynamics problems and quasi - brittle materials (such as concrete) in order to reduce the associated computational costs. Specifically, the authors combine the LArge Time INcrement (LATIN) method with the Proper Generalized Decomposition (PGD) technique to solve non - linear low - frequency dynamics problems with a quasi - brittle isotropic damage constitutive relationship.

### Key points of the solution:

1. **Combining the LATIN and PGD methods**:
- **LATIN method**: It is used to handle non - linear problems and gradually approximates the solution of the problem through alternating non - linear and linear steps.
- **PGD technique**: It provides a numerically economical low - rank approximate solution by the method of separating variables, thus reducing the amount of calculation.

2. **Time - discrete method**:
- Use the **Time - Discontinuous Galerkin Method (TDGM)** to calculate the time contribution of the space - time separated variable solution. This method is especially effective when dealing with problems in long - time intervals, such as long - time excitation or fatigue problems in earthquake engineering.

3. **Application background**:
- The focus is on concrete materials because they play a key role in seismic - resistant structural design, but the proposed method is general and applicable to any material with non - linear behavior.

4. **Objective**:
- Improve computational efficiency so that a single earthquake event can be quickly simulated and complex parameter studies can be carried out to evaluate the variability of different parameters and their impact on structural performance.

### Formula representation:

- **Dynamic equilibrium equation**:
\[
\rho \ddot{u}=\text{div}(\sigma)+\rho f
\]
where \(\rho\) is the material density, \(\ddot{u}\) is the second - order time derivative of the displacement field, \(\sigma\) is the stress tensor, and \(f\) is the body force.

- **Constitutive equation**:
\[
\tilde{\sigma}_{me}=\frac{\sigma_{me}}{1 - d}
\]
\[
\tilde{\sigma}_{cr}=\frac{\sigma_{cr}}{d}
\]
\[
\sigma=\sigma_{me}+\sigma_{cr}
\]
where \(d\) is the damage variable, and \(\sigma_{me}\) and \(\sigma_{cr}\) are the effective stress tensor and the stress tensor of closed micro - cracks respectively.

Through the combination of these methods and techniques, the authors hope to significantly reduce the computational cost of non - linear dynamics problems while ensuring computational accuracy, especially in earthquake engineering applications.