# Parameter identification of two-time-scale nonlinear transient models

Guillaume Puel, Denis Aubry

► To cite this version:

Guillaume Puel, Denis Aubry. Parameter identification of two-time-scale nonlinear transient models. Advanced Modeling and Simulation in Engineering Sciences, 2015, 2, 10.1186/s40323-015-0030-z. hal-01832555

HAL Id: hal-01832555

https://hal.science/hal-01832555

Submitted on 13 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Parameter identification of two-time-scale nonlinear transient models

Guillaume Puel* and Denis Aubry

*Correspondence:
guillaume.puel@centralesupelec.fr
Laboratoire MSSMat
(CentraleSupélec-CNRS
UMR8579-Université
Paris-Saclay), Grande
voie des vignes,
92290 Châtenay-Malabry,
France

**Abstract**

The aim of this paper is to study two-time-scale nonlinear transient models and their associated parameter identification. When it is possible to consider two well-separated time scales, and when the fast component of the applied loading is periodic, a periodic time homogenization scheme, similar to what exists in space homogenization, can be used to derive a homogenized model. A parameter identification process for this latter is then proposed, and applied to an academic example, which allows to show the benefits of such a strategy.

**Keywords:** Inverse problem, Identification, Adjoint state, Time multiscale, Periodic homogenization

**Background**

**Context**

Material fatigue is one of the key phenomena whose study is essential to design and develop modern industrial parts: many of the classical methods used in the industry, such as those described in Ref. [1], allow to estimate fatigue life through cumulative formulas. However, these latter do not take into account the history associated with the loading, and tend to give poor results as soon as the loading frequencies are high and the inertia effects should be considered. One particularly difficult case of study corresponds to combined cycle fatigue (CCF), where two different periodic loads (generally, one being 'slow', the other 'fast') are involved. In the meantime, more and more complex material laws have been proposed to describe cyclic behaviors, and the ideal way to use these models consists in solving them in the time domain.

As far as the numerical calculation of a time-dependent model is concerned, the question of the computational cost can be of utmost relevance, especially when the considered model deals with fast phenomena, which require the use of very small time steps, when compared with the length of the time interval of study. In order to drastically reduce the computational cost, a periodic time homogenization method can be used when two well-separated time scales exist, and when the components of the applied loading are periodic with respect to the fast time scale [2]. The resulting homogenized model can then be solved with a drastically reduced computational cost, using time steps related to the slow time scale only, whereas the fast time scale is taken into account in an average way in the homogenization scheme.

<figure><img src="image_4.png" /></figure>
In order to give accurate predictions, such time-homogenized models have to be compared with experimental data. The key point is to define an identification strategy able to deal with such models and using a process that remains cheap and efficient. The aim of this paper is to analyze on a specific academic example how a relevant identification strategy can be proposed, and to comment on the different choices that can be made throughout the whole identification process. First, the basic elements of the periodic time homogenization method are described, as well as one particular case of study showing the efficiency of the method; then, the general parameter identification strategy that will be adapted to the case of time-homogenized models is introduced. After this introduction, the main part of the paper focuses on a simple academic example allowing to give the theoretical developments associated with a truly two-time-scale identification strategy, and to analyze the associated numerical results upon which conclusions can be drawn.

## Periodic time homogenization method

### Main ingredients

Periodic time homogenization, as it was initially proposed in Ref. [2], can be seen as a transposition to time of the classical periodic space homogenization methods, such as the techniques described in Ref. [3] or in Ref. [4]: it is typically based on the introduction of a two-scale asymptotic expansion for every time-dependent variable sought for in the equations of the reference problem:

$$ \alpha(t, \tau) = \alpha_0(t, \tau) + \xi\alpha_1(t, \tau) + \xi^2\alpha_2(t, \tau) + O(\xi^3) \quad (1) $$

where $t$ stands for the so-called slow time scale and $\tau = t/\xi$ for the fast time scale. Indeed, when the ratio $\xi$ is very small, it is possible to separate the two time scales by considering that they are independent one from the other, implying that any derivative with respect to time has to use the partial derivatives with respect to the two time scales:

$$ d_t \alpha = \partial_t \alpha + \frac{1}{\xi} \partial_{\tau} \alpha = \dot{\alpha} + \frac{1}{\xi} \alpha' \quad (2) $$

where $d_t \bullet, \partial_t \bullet = \dot{\bullet}$ and $\partial_{\tau} \bullet = \bullet'$ stand for the total time derivative, the partial derivative with respect to the slow time and the partial derivative with respect to the fast time respectively. Moreover, if the applied loading has a component, which is periodic with respect to the fast time $\tau$, it is possible to assume that each variable is quasiperiodic, meaning that it is periodic with respect to $\tau$. Figure 1 gives an illustration of such a quantity. This quasiperiodicity assumption is all the more justified as the ratio $\xi$ is small, and can then be written as:

$$ \alpha(t, \tau) = \alpha \left( t, \tau + \frac{T_f}{\xi} \right) \quad \forall t, \tau \quad (3) $$

where $T_f$ is the period associated with the fast periodic component of the loading.

Using the asymptotic expansion (1) for each quantity in the reference equations, and balancing the different orders of $\xi$, the time-homogenized equations are determined by averaging over a fast period the different quantities:

$$ \langle \alpha \rangle(t) = \frac{\xi}{T_f} \int_0^{T_f/\xi} \alpha(t, \tau) d\tau \quad (4) $$
<figure><img src="image_5.png" /><figcaption>Figure 1 Periodic time homogenization. Illustration of the two main assumptions required in the method.</figcaption></figure>

This fast-time average actually allows to separate slow-evolving phenomena and fast-time periodic components by using the fact that the quasiperiodicity assumption can actually be rewritten for any time-dependent variable as:

$$ \langle \alpha' \rangle = 0 \qquad (5) $$

In addition to the different time-homogenized quantities $\langle \alpha \rangle$, the residuals associated with this fast-time average are then denoted as $\alpha^* = \alpha - \langle \alpha \rangle$, and depend on both time scales $t$ and $\tau$ a priori. They usually have to verify a very simple problem (typically a linear, elastic one), which will be solved for a time interval corresponding to one fast period $T_f$. Eventually, the obtained time-homogenized equations are solved relatively to the slow time scale only, by introducing the averaged influence of the fast cycles corresponding to the solution of the fast problem at each slow time step, hence allowing to solve for all the zeroth-order time-homogenized variables with a drastic reduction in the computational cost when compared with what would be required to solve the reference problem (for example by using the classical rule of thumb of 20 time steps per fast cycle).

## Cases of study

References on the time homogenization method still tend to be quite scarce, as showed a recent tentative review in Ref. [5]. Nevertheless, the possible applications are numerous, as we can see with a significant number of references dealing with topics as diverse as the behavior of viscoelastic materials in Ref. [6], the vibration of preloaded beams in Ref. [7], the ultrasonic welding of composites in Ref. [8], the ultrasonic imaging of tissues in Ref. [9], the cyclic loading of normally consolidated clay in Ref. [10] or the viscoelastic-visco-plastic behavior of polymers in Ref. [11].

Recently, we focused on validating the method for different cases of simulations of structures withstanding fatigue loads with two periodic components:

* material fatigue with a viscoplastic law defining two hardening variables in Ref. [12];

* material fatigue with an isotropic damage law in Ref. [13];

* extension of the method to three different time scales in Refs. [5] and [14];

* extension of the method to the case of a resonant excitation in Ref. [5].
In this latter case, a real case example, coming from the European Project PREMECCY on fatigue, related to gas turbine blades, and briefly described in Ref. [15], was presented: indeed, these latter, during operation, withstand a loading with two periodic components whose time evolution is sketched in Figure 2. The PREMECCY project's principle was to experimentally study CCF with several blade-shaped specimens, such as the one depicted in Figure 3. Numerical simulations using the periodic time homogenization method have been run on such a geometry, with a loading implying frequencies

<figure><img src="image_6.png" /><figcaption>**Figure 2** PREMECCY project. Schematic representation of the loading, after [15].</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 3** Numerical results for a specific specimen. Bending II mode at 1,400 Hz and final local longitudinal plastic strain after 100 slow cycles and 1,000,000 fast cycles.</figcaption></figure>
$F = 0.14$ Hz and $F/\xi = 1,400$ Hz, this latter value corresponding to the specimen's second bending eigenfrequency. A specific result is shown in Figure 3, where it is easily seen that the highest values of the longitudinal plastic strain are located at the specimen's edges, and that these values nearly vanish everywhere else. More details can be found in Ref. [5].

### Gradient-based parameter identification using adjoint state formulations

To improve the quality of the model's predictions, especially the fatigue life estimate, it is mandatory to compare computed quantities with experimental data and to define an identification scheme allowing to update the model's parameters judiciously. In this section, a general description of the proposed identification strategy is given, which, as we will see, is relevant for both reference and time-homogenized problems: it will thus be illustrated in this latter case in the “Methods”.

The forward problem is considered as an implicit formulation with a vector function $\mathcal{F}$ over a time interval $[0, T]$:

$$ 
\begin{aligned}
\mathcal{F}(\mathbf{u}(t), d_t\mathbf{u}(t), d_t^2\mathbf{u}(t), \mathbf{p}, t) &= \mathbf{0} \\
\mathbf{u}(0) &= \mathbf{U}_0 \\
d_t\mathbf{u}(0) &= \mathbf{V}_0
\end{aligned}
\qquad (6) $$

where $d_t$ and $d_t^2$ are the first- and second-order time derivatives respectively. $\mathbf{U}_0$ and $\mathbf{V}_0$ stand for the initial conditions of the dynamic problem. Whereas $\mathbf{u}$ is the state vector of size $N$, composed of all the time-dependent degrees of freedom (DOFs) describing the studied problem, $\mathbf{p}$ stands for the vector containing the $P$ scalar parameters associated with the differential equation (6).

### Formulation of the identification problem

The identification problem consists in finding the parameter vector $\mathbf{p}_{opt}$ such that the solution $\mathbf{u}(t; \mathbf{p}_{opt})$ of (6) obtained with the parameters $\mathbf{p}_{opt}$ is as close to the available experimental data as possible. These latter are indeed compared with the corresponding quantities $\mathbf{A}\mathbf{u}(t; \mathbf{p})$, where $\mathbf{A}$ is a projection operator allowing to select, for each quantity, the closest DOF to the experimental measurement point. In order to use consistent notations, the corresponding experimental quantity is denoted $\mathbf{A}\mathbf{u}_{exp}(t)$; however, it does not mean that such a vector $\mathbf{u}_{exp}(t)$ actually exists.

The following misfit function is then introduced: it consists of a norm measuring the discrepancy between the quantities predicted with the forward model (6) and experimental data:

$$ \mathcal{J}(\mathbf{p}) = \frac{1}{2} \int_{0}^{T} \left| \mathbf{A}(\mathbf{u}(t; \mathbf{p}) - \mathbf{u}_{exp}(t)) \right|^{2} dt + \frac{1}{2} \left| \mathbf{R}(\mathbf{p} - \mathbf{p}_{0}) \right|^{2} \qquad (7) $$

where $\mathbf{u}(t; \mathbf{p})$ satisfies Eq. (6). The $L^2$-norm proposed here is completed with a Tikhonov regularization term, allowing to deal with the ill-posedness of the identification problem, by bounding the magnitude of the parameter vector $\mathbf{p}$ to be identified: this regularization term uses a vector $\mathbf{p}_0$ containing nominal values corresponding to a priori experience,
and a diagonal weighing matrix $\mathbf{R}$. Eventually, the solution of the identification problem can be sought as the parameter vector $\mathbf{p}_{opt}$ minimizing the misfit function $\mathcal{J}(\mathbf{p})$:

$$ \mathbf{p}_{opt} = \arg \min_{\mathbf{p}} \mathcal{J}(\mathbf{p}) \qquad (8) $$

### Adjoint state formulation

The determination of this minimum is achieved using gradient-based minimization methods, therefore the question of avoiding local minima by means of an appropriate regularization process should be carefully addressed. In some cases, rather than using the classical Tikhonov regularization term, the a priori experience can be introduced in some specific ways, as in Ref. [16]. Similarly, the fact of using a homogenized model in the parameter identification process can introduce a regularizing effect, just as explained in Ref. [17]. However, we will not address here this specific question, but rather focus on the identification process itself.

To estimate the gradient of the misfit function, we solve here an adjoint state problem. A typical example in mechanical engineering is given in Ref. [18], where the parameters of an elastoplastic material law are identified with indentation tests. In the strategy proposed here, the generic form of the adjoint state problem, which can be obtained by expressing the stationarity of a Lagrangian as in Ref. [19], is as follows:

$$ \begin{gathered} \nabla_{\mathbf{u}} \mathcal{F}^{\mathrm{T}} \mathbf{z} - d_t (\nabla_{d_t \mathbf{u}} \mathcal{F}^{\mathrm{T}} \mathbf{z}) + d_t^2 (\nabla_{d_t^2 \mathbf{u}} \mathcal{F}^{\mathrm{T}} \mathbf{z}) = \mathbf{A}^{\mathrm{T}} \mathbf{A} (\mathbf{u} - \mathbf{u}_{\mathrm{exp}}) \\ (\nabla_{d_t^2 \mathbf{u}} \mathcal{F}^{\mathrm{T}} \mathbf{z})|_{t=T} = \mathbf{0} \\ (\nabla_{d_t \mathbf{u}} \mathcal{F}^{\mathrm{T}} \mathbf{z})|_{t=T} = \mathbf{0} \end{gathered} \qquad (9) $$

where $\nabla_{\mathbf{u}} \mathcal{F}$, $\nabla_{d_t \mathbf{u}} \mathcal{F}$ and $\nabla_{d_t^2 \mathbf{u}} \mathcal{F}$ stand for the directional derivatives of $\mathcal{F}$ with respect to $\mathbf{u}$, $d_t \mathbf{u}$ and $d_t^2 \mathbf{u}$ respectively. The adjoint state problem is then a time-backward differential equation with two final conditions, and where the first-order sensitivities of the forward problem are concerned.

Once the adjoint state problem (9) is solved, it can be shown that the misfit function's gradient with respect to the parameter vector $\mathbf{p}$ can be expressed as:

$$ \nabla_{\mathbf{p}} \mathcal{J}(\mathbf{p}) = \mathbf{R}^{\mathrm{T}} \mathbf{R}(\mathbf{p} - \mathbf{p}_0) - \int_{0}^{T} \nabla_{\mathbf{p}} \mathcal{F}^{\mathrm{T}} \mathbf{z}(t) \mathrm{d}t \qquad (10) $$

This specific way of estimating the misfit function's gradient can be compared with a classical finite difference formula, such as the central finite difference scheme: in this latter case, when the parameter vector is of size $P$, the gradient calculation is obtained by evaluating the misfit function in $2P$ additional 'points', each couple of points corresponding to two symmetrical perturbations of the misfit function associated with each parameter in the vector $\mathbf{p}$. The resulting computational cost for each gradient evaluation consists of $2P$ solutions of the forward problem (6) and $2P$ time integral evaluations. By contrast, when the adjoint state solution is used, only two differential equation solutions are required: one for the forward problem (6) and one for the adjoint problem (9). The associated computational cost for each gradient evaluation is then two differential equation solutions, and $P$ time integral evaluations. The resulting gain is as high as the number of parameters to be identified. Moreover, it is easier to control the accuracy of the
gradient estimate with the adjoint state method than with finite difference formulas, for which the choice of the discretization steps has a strong influence on the final estimate.

## Methods

We propose to study an academic example in order to discuss the different steps of the parameter identification problem associated with a time-homogenized model. To the best of our knowledge, it is the first time that such a two-time-scale identification strategy using an adjoint state formulation is proposed: only one alternative has been briefly studied in Ref. [20] with the use of a genetic method. First, we will recall how the identification process is formulated on a classical reference problem, i.e. with no time homogenization. Then, we will apply the identification strategy on the time-homogenized version of the problem, and see what connections exist with the identification process associated with the reference problem.

The case of study detailed here consists of a straight bar of length L, clamped at one end, and withstanding at the other end a normal force with two periodic components. The measured displacement at this end is then used to determine the parameter values of the material elastic viscoplastic law.

## Reference problem

### Forward problem

The dynamic equilibrium of the bar is a scalar equation, with the normal (nonviscous) stress $\sigma(x, t)$ and the longitudinal displacement $u(x, t)$ defined at each point $x \in [0, L]$ and for any $t \in [0, T]$:

$$ \partial_x \sigma + c_K d_t \partial_x \sigma = \rho d_t^2 u \quad (11) $$

where $\partial_x$ is the partial space derivative, $\rho$ is the mass density, and $c_K$ is the damping ratio where damping is assumed as proportional to stiffness. Homogeneous initial conditions for the displacement and the velocity are assumed. No load is applied along the bar, and, whereas this latter is clamped at $x = 0$, a surface force $f_s(t)$ is applied at $x = L$ for any $t \in [0, T]$:

$$ u|_{x=0} = 0 \quad (12) $$

$$ (\sigma + c_K d_t \sigma)|_{x=L} = f_s \quad (13) $$

This surface force is biperiodic, with the ratio of the associated low frequency over the fast one equal to $10^{-4}$.

The constitutive relation links the normal (nonviscous) stress to the longitudinal strain along the bar for any $(x, t) \in [0, L] \times [0, T]$:

$$ \sigma = E(\partial_x u - \varepsilon^P) \quad (14) $$

where $E$ is the Young's modulus, and $\varepsilon^P(x, t)$ stands for the longitudinal plastic strain. One assumes that this latter verifies a Norton's viscoplastic evolution law for any $(x, t) \in [0, L] \times [0, T]$, which depends on the nonviscous stress only:

$$ d_t \varepsilon^P = \left( \frac{|\sigma|}{K} \right)^n \text{sign} \sigma \quad (15) $$
with $\epsilon^P|_{t=0} = 0 \forall x \in [0, L]$. K and n are the two constant parameters to be identified.

The numerical implementation consists of quadratic 1D finite elements under MAT-LAB for the spatial discretization of the displacement, whereas the longitudinal plastic strain is linearly interpolated using the 'external' nodes of the different elements. The time integration uses the 'ode45' procedure, based on a fourth-order embedded Runge-Kutta formula according to Ref. [21]: depending on the number of finite elements used, the time step to be chosen can verify the classical rule of thumb of 20 time steps per fast cycle, or must be smaller than the maximal time step associated with the Courant-Friedrichs-Lewy condition, making in both cases the computational cost become unaffordable as soon as a high number of cycles is calculated.

### Parameter identification process

It is assumed that experimental data consist of the knowledge of the displacement $u_{exp}(t)$ measured at $x = L$. The proposed misfit function is then:

$$ \mathcal{J}(K, n) = \frac{1}{2} \int_{0}^{T} (u(L, t; K, n) - u_{exp}(t))^2 dt \\ + \frac{\alpha_{K}}{2}(K - K_{0})^2 + \frac{\alpha_{n}}{2}(n - n_{0})^2 \quad (16) $$

where $u(x, t; K, n)$ is the solution of the forward problem (11)–(15) with parameter values K and n.

In order to derive the adjoint state problem, the following Lagrangian is introduced:

$$ \begin{aligned} \mathcal{L}(u, \epsilon^P, K, n, z, \lambda) = & \frac{1}{2} \int_0^T (u(L, t) - u_{exp}(t))^2 dt \\ & + \frac{\alpha_K}{2}(K - K_0)^2 + \frac{\alpha_n}{2}(n - n_0)^2 \\ & - \int_0^T \int_0^L \rho(d_t^2 u) z \, dx \, dt + \int_0^T f_s z|_{x=L} \, dt \\ & - \int_0^T \int_0^L [E(\partial_x u - \epsilon^P) + c_K E d_t(\partial_x u - \epsilon^P)] \partial_x z \, dx \, dt \\ & + \int_0^T \int_0^L \left[ d_t \epsilon^P - \left( \frac{E|\partial_x u - \epsilon^P|}{K} \right)^n \operatorname{sign}(\partial_x u - \epsilon^P) \right] \lambda \, dx \, dt \end{aligned} \quad (17) $$

where $(u, \epsilon^P, K, n, z, \lambda)$ are considered as independent quantities. The minimization of the misfit function (16) under the constraints (11)–(15) is then equivalent to express the stationarity of the Lagrangian (17) with no constraint [19].

The adjoint states $z(x, t)$ and $\lambda(x, t)$ are Lagrange multipliers verifying, for any $(x, t) \in (0, L) \times [0, T]$, equations coming from the stationarity of the Lagrangian with respect to u:

$$ \partial_x \left[ E \partial_x z - c_K E d_t \partial_x z + n \frac{E}{K} \left( \frac{E |\partial_x u - \epsilon^P|}{K} \right)^{n-1} \lambda \right] = \rho d_t^2 z \quad (18) $$

$$ \left( E \partial_x z - c_K E d_t \partial_x z + n \frac{E}{K} \left( \frac{E |\partial_x u - \epsilon^P|}{K} \right)^{n-1} \lambda \right) \Big|_{x=L} = u|_{x=L} - u_{exp} \quad (19) $$
$$z_{|t=T} = 0, \quad d_t z_{|t=T} = 0 \qquad (20)$$

and $\varepsilon^P$:

$$d_t \lambda = n \frac{E}{K} \left( \frac{E |\partial_x u - \varepsilon^P|}{K} \right)^{n-1} \lambda + E |\partial_x z - c_K E d_t \partial_x z \qquad (21)$$

$$\lambda_{|t=T} = 0 \qquad (22)$$

Formally, the latter equations correspond to a dynamic elasto-viscoplastic problem (with negative damping), when one considers the first Lagrange multiplier $z$ as a kind of displacement (even if it is not homogeneous to a length) and the second Lagrange multiplier $\lambda$ as a kind of plastic strain (even it is not dimensionless). The discrepancy between the forward model's predictions and experimental data is directly imposed as a Neumann boundary condition at $x = L$, i.e. where the measurements are available. One should notice that this problem is time-backward with final conditions (at $t = T$): once the adjoint state problem has been properly discretized with respect to space, the obtained time differential equations can be solved with classical time integration methods by applying a change in variables $\theta = T - t$ in order to deal with initial conditions rather than final ones.

Once solved, the adjoint state solution is used to evaluate the misfit function's gradient, whose components are equal to the partial derivatives of the Lagrangian (17) with respect to the parameters:

$$\begin{aligned} \frac{\partial \mathcal{J}}{\partial K}(K, n) &= \frac{\partial \mathcal{L}}{\partial K}(u(K, n), \varepsilon^P(K, n), K, n, z(u, \varepsilon^P), \lambda(u, \varepsilon^P)) \\ &= \int_0^T \int_0^L \frac{n}{K} \left( \frac{E |\partial_x u - \varepsilon^P|}{K} \right)^n \mathrm{sign} \left( \partial_x u - \varepsilon^P \right) \lambda \, dx \, dt \\ &\quad + \alpha_K (K - K_0) \end{aligned} \qquad (23)$$

$$\begin{aligned} \frac{\partial \mathcal{J}}{\partial n}(K, n) &= \frac{\partial \mathcal{L}}{\partial n}(u(K, n), \varepsilon^P(K, n), K, n, z(u, \varepsilon^P), \lambda(u, \varepsilon^P)) \\ &= - \int_0^T \int_0^L \left( \frac{E |\partial_x u - \varepsilon^P|}{K} \right)^n \mathrm{sign} \left( \partial_x u - \varepsilon^P \right) \log \left( \frac{E |\partial_x u - \varepsilon^P|}{K} \right) \lambda \, dx \, dt \\ &\quad + \alpha_n (n - n_0) \end{aligned} \qquad (24)$$

One should notice that, beside the forward state solutions, only the second Lagrange multiplier $\lambda$ appears in the previous expressions: this is logical here, because one is interested in the parameters $(K, n)$ of the evolution law (15) only. The first Lagrange multiplier $z$ would have appeared in these expressions if one had tried to identify parameters associated with the dynamic equilibrium equation, such as the Young's modulus $E$ or the damping ratio $c_K$.

These estimates of the misfit function's gradient components are used by the minimization algorithm to determine the best search direction and step size. However, even if the strategy using the adjoint state problem is the cheapest one in terms of computational cost, this latter can be awfully prohibitive if a high number of fast cycles has to
calculated, as it is the case here, since viscoplasticity is a slow-evolving phenomenon, implying that the change in the longitudinal viscoplastic strain is very small for each individual fast cycle.

### Time-homogenized problem

In order to obtain a drastic reduction in the computational cost, the periodic time homogenization method is applied. Following the different steps briefly described in the “Back-ground”, the time-homogenized equations to be solved are obtained.

### Dynamic equilibrium equation

The first equation to be dealt with is the dynamic equilibrium equation, which can be rewritten in a normalized way in order to compare the orders of magnitude of the different terms:

$$ \partial_{\hat{x}} \hat{\sigma} + F c_K d_t \partial_{\hat{x}} \hat{\sigma} = \frac{\rho L^2 F^2}{E} d_t^2 \hat{u} \quad (25) $$

where F is the frequency of the slow component of the loading, and the different normalized quantities are denoted as $\hat{\bullet}$. Different cases arise, depending on the orders of magnitude of $\rho L^2 F^2/E$ and $F c_K$. Here, we assume that:

$$ \frac{\rho L^2 F^2}{E} = \beta \xi^2 \quad (26) $$

$$ F c_K = \gamma \xi \quad (27) $$

where $\beta \le O(1)$ and $\gamma \le O(1)$. Physically, the first assumption (26) is equivalent to having $L/\lambda = \sqrt{\beta}$, where $\lambda$ is the wavelength associated with the fast loading frequency: this means that the order of magnitude of $\lambda$ should be at least equal to the length of the bar.

Using these two assumptions and introducing an asymptotic expansion of the displacement $u = u_0 + \xi u_1 + \xi^2 u_2 + O(\xi^3)$ allows to express the acceleration as follows:

$$ \rho d_t^2 u = \frac{\beta E}{L^2 F^2} u_0'' + \xi \frac{\beta E}{L^2 F^2} (2\dot{u}_0' + u_1'') + \xi^2 \frac{\beta E}{L^2 F^2} (\ddot{u}_0 + 2\dot{u}_1' + u_2'') + O(\xi^3) \quad (28) $$

It can then be shown [5] that the two time scales are separable and that the homogenization method can be applied: the zeroth-order dynamic equilibrium equation becomes:

$$ \partial_x \sigma_0 + \frac{\gamma}{F} \partial_x \sigma_0' = \frac{\beta E}{L^2 F^2} u_0'' \quad (29) $$

and the associated zeroth-order boundary condition at $x = L$ reads:

$$ (\sigma_0 + \frac{\gamma}{F} \sigma_0')_{|x=L} = f_s \quad (30) $$

Eventually, by taking the fast-time average (4) of the two previous equations, and by using the quasiperiodicity assumption (5), one can obtain the time-homogenized zeroth-order dynamic equilibrium equation:

$$ \partial_x \langle \sigma_0 \rangle = 0 \quad (31) $$
for any $(x, t) \in (0, L) \times [0, T]$, and the time-homogenized zeroth-order boundary condition:

$$ \langle \sigma_0 \rangle_{x=L} = \langle f_s \rangle \qquad (32) $$

which corresponds to the solution of a quasistatic problem.

### Evolution law

In the same way, the evolution law (15) gives, up to order zero:

$$ \frac{1}{\xi} \epsilon_0^{p'} + (\dot{\epsilon}_0^p + \epsilon_1^{p'}) + O(\xi) = \left(\frac{|\sigma_0|}{K}\right)^n \text{sign} \sigma_0 + O(\xi) \qquad (33) $$

First, the $1/\xi$-order has to be considered, giving: $\epsilon_0^{p'} = 0$. This is equivalent to say that the zeroth-order plastic strain only depends on the slow time variable:

$$ \epsilon_0^p(x, t, \tau) = \epsilon_0^p(x, t) \qquad (34) $$

The physical interpretation of this latter result is that viscoplasticity is a slow-evolving phenomenon, even if the material withstands a high-frequency loading. This slow variation is described by the zeroth-order evolution law, which can be expressed as:

$$ \dot{\epsilon}_0^p + \epsilon_1^{p'} = \left(\frac{|\sigma_0|}{K}\right)^n \text{sign} \sigma_0 \qquad (35) $$

In order to make the contribution of the first-order plastic strain disappear from this relation, the fast-time average (4) of this latter is evaluated. Since, on the one hand, Eq. (34) is equivalent to say that $\langle \dot{\epsilon}_0^p \rangle = \dot{\epsilon}_0^p$, and, on the other hand, the quasiperiodicity assumption (5) implies that $\langle \epsilon_1^{p'} \rangle = 0$, this average eventually gives that, for any $(x, t) \in [0, L] \times [0, T]$:

$$ \dot{\epsilon}_0^p = \left\langle \left( \frac{|\sigma_0|}{K} \right)^n \text{sign} \sigma_0 \right\rangle \qquad (36) $$

with a zero initial condition.

### Zeroth-order time-homogenized problem

The final system consists of the following equations to be solved at the slow time scale only, for any $(x, t) \in (0, L) \times [0, T]$:

$$ \partial_x \langle \sigma_0 \rangle = 0 \qquad (37) $$

$$ \langle \sigma_0 \rangle = E(\partial_x \langle u_0 \rangle - \epsilon_0^p) \qquad (38) $$

$$ \langle u_0 \rangle_{x=0} = 0 \qquad (39) $$

$$ \langle \sigma_0 \rangle_{x=L} = \langle f_s \rangle \qquad (40) $$

$$ \dot{\epsilon}_0^p = \left\langle \left( \frac{|\sigma_0|}{K} \right)^n \text{sign} \sigma_0 \right\rangle \qquad (41) $$
The resulting system then corresponds to a quasistatic, elastic viscoplastic problem, where the loading consists of the slow component of the surface force.

Of course, the fast component does have an effect on these homogenized equations, by means of the evolution law (41): since this latter is nonlinear, it is required to calculate the zeroth-order 'instantaneous' stress field $\sigma_0$ instead of directly using the corresponding homogenized quantity $\langle\sigma_0\rangle$. To do this, it is mandatory to estimate the residual quantity $\sigma_0^*(x, t, \tau)$, defined according to the following relation:

$$ \sigma_0(x, t, \tau) = \langle\sigma_0\rangle(x, t) + \sigma_0^*(x, t, \tau) \quad (42) $$

This residual is then the solution of the following system of residual equations:

$$ \partial_x \sigma_0^* + \frac{\gamma}{F} \partial_x \sigma_0^{*'} = \frac{\beta E}{L^2 F^2} u_0^{*''} \quad (43) $$

$$ \sigma_0^* = E \partial_x u_0^* \quad (44) $$

$$ u_0^{*|x=0} = 0 \quad (45) $$

$$ \left( \sigma_0^* + \frac{\gamma}{F} \sigma_0^{*'} \right)_{|x=L} = f_s^* \quad (46) $$

obtained by subtracting the equations of the time-homogenized problem (37)–(40) to the zeroth-order equations of the reference problem. This residual system consists in solving, for the fast time scale $\tau$, and at each slow time step $t_k$, a dynamic problem for a purely viscoelastic material, with a surface force $f_s^* = f_s - \langle f_s \rangle$. Since viscoplasticity (which is the only source of nonlinearity here) appears in the zeroth-order time-homogenized equations only, this residual problem is linear, and can be conveniently solved in the frequency domain. If needed, additional orders could be addressed, as detailed in Ref. [5].

The evaluation of the fast-time average appearing in (41) can be achieved by means of a numerical integration formula, such as the trapezoidal rule with $N + 1$ points, defined as:

$$ \langle\alpha\rangle(t) \approx \frac{1}{N} \left( \frac{\alpha(t, 0)}{2} + \sum_{j=1}^{N-1} \alpha\left(t, \frac{j}{N} \frac{1}{F}\right) + \frac{\alpha\left(t, \frac{1}{F}\right)}{2} \right) \quad (47) $$

This choice is not only related to its formal simplicity, but can actually be justified using the conclusions presented in Ref. [22]: the author shows that the trapezoidal rule converges extremely fast when integrating smooth periodic functions, as it is the case here, and that no substantial additional gain would be obtained with the use of more elaborated formulas, such as Simpson's for example.

### Associated parameter identification strategy

The identification strategy using adjoint state formulations described in the “Background” is applied now on the time-homogenized problem. First, a consistent misfit function has to
be selected, then the adjoint state problem required to estimate the misfit function's gradient is introduced.

### Choice of the misfit function

The first step consists in describing, through the misfit function, the discrepancy between the time-homogenized model's predictions and experimental data: indeed, on the one hand, the model has been solved on slow time steps only, whereas, on the other hand, the experimental data can be available on a much finer scale.

The most efficient choice in terms of computation cost is to use time-homogenized quantities in the misfit function (16), since it allows to address the time integral on slow time steps only:

$$ \mathcal{J}_0(K, n) = \frac{1}{2} \int_0^T |\langle u_0 \rangle(L, t; K, n) - \langle u_{exp} \rangle(t)|^2 dt \\ + \frac{\alpha_K}{2}(K - K_0)^2 + \frac{\alpha_n}{2}(n - n_0)^2 \quad (48) $$

Whereas $\langle u_0 \rangle(L, t; E, K, n)$ is the solution of the zeroth-order time-homogenized forward problem (37)–(41), $\langle u_{exp} \rangle(t)$ stands for the corresponding experimental quantity, which is obtained by fast-averaging the experimental data for each slow time step $t_k$ of the time-homogenized displacement:

$$ \langle u_{exp} \rangle(t_k) = \frac{F}{\xi} \int_{t_k}^{t_k + \frac{\xi}{F}} u_{exp}(t) dt \quad (49) $$

### Adjoint state problem

The misfit function's gradient is evaluated in the same way as in the section associated with the reference problem, using the solution of an adjoint state problem. Starting from the Lagrangian associated with the misfit function (48) and the time-homogenized forward problem (37)–(41), it is first obtained that the second Lagrange multiplier depends on the slow time scale only:

$$ \lambda_0(x, t, \tau) = \lambda_0(x, t) \quad (50) $$

because it is associated with the zeroth-order time-homogenized evolution law (41), which depends on $t$ only. On the contrary, the first Lagrange multiplier is decomposed as $\langle z_0 \rangle(x, t) + z_0^*(x, t, \tau)$ to deal with the two systems (37)–(40) and (43)–(46) respectively.

Then, the two following PDEs corresponding to the adjoint state problem are obtained, for any $(x, t) \in (0, L) \times [0, T]$:

$$ \partial_x \left[ E \partial_x \langle z_0 \rangle + n \frac{E}{K} \left\langle \left( \frac{E |\partial_x u_0 - \varepsilon_0^P|}{K} \right)^{n-1} \right\rangle \lambda_0 \right] = 0 \quad (51) $$

$$ \left( E \partial_x \langle z_0 \rangle + n \frac{E}{K} \left\langle \left( \frac{E |\partial_x u_0 - \varepsilon_0^P|}{K} \right)^{n-1} \right\rangle \lambda_0 \right) \Bigg|_{x=L} = \langle u_0 \rangle_{|x=L} - \langle u_{exp} \rangle \quad (52) $$
$$ \dot{\lambda}_0 = n \frac{E}{K} \left\langle \left( \frac{E |\partial_x u_0 - \varepsilon_0^P|}{K} \right)^{n-1} \right\rangle \lambda_0 + E \partial_x \langle z_0 \rangle \quad (53) $$

$$ \lambda_0|_{t=T} = 0 \quad (54) $$

and they actually correspond to the zeroth-order time-homogenized versions of Eqs. (18)–(22) giving the adjoint state solutions associated with the identification process for the reference problem. Figures 4 and 5 show the comparison between the two

<figure><img src="image_8.png" /><figcaption>**Figure 4** Adjoint state. Adjoint state solutions for reference (in red) and time-homogenized (in blue) problems, at $x = 0$ and $x = L$.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 5** Adjoint state. Zoom of Figure 4 (same color conventions).</figcaption></figure>
corresponding adjoint state solutions $\lambda$ and $\lambda_0$ for the identification process detailed in the “Results and discussion”: the adjoint state solution $\lambda_0$ corresponding to the time-homogenized problem is homogeneous along the bar, and is very close to the (homogeneous) fast-time average of the (heterogeneous) adjoint state solution $\lambda$ associated with the reference problem.

Indeed, in this case, the previous equations can be decoupled in order to get the PDE verified by the second Lagrange multiplier $\lambda_0$ for any $(x, t) \in [0, L] \times [0, T]$:

$$ \partial_x \dot{\lambda}_0 = 0 \qquad (55) $$

$$ \dot{\lambda}_0|_{x=L} = \langle u_0 \rangle_{x=L} - \langle u_{exp} \rangle \qquad (56) $$

$$ \lambda_0|_{t=T} = 0 \qquad (57) $$

finally implying that, for any $(x, t) \in [0, L] \times [0, T]$:

$$ \dot{\lambda}_0 = \langle u_0 \rangle_{x=L} - \langle u_{exp} \rangle \qquad (58) $$

$$ \lambda_0|_{t=T} = 0 \qquad (59) $$

which is a time-backward ODE with a final condition equal to zero. This equation can be solved using the slow time steps $t_k$ only, which allows to derive the solution in a way as efficient as for the time-homogenized forward solution.

### Misfit function's gradient

The misfit function's gradient then consists of the two following partial derivatives:

$$ \frac{\partial \mathcal{J}_0}{\partial K} = \int_0^T \int_0^L \left\langle \frac{n}{K} \left( \frac{E|\partial_x u_0 - \varepsilon_0^p|}{K} \right)^n \text{sign} \left( \partial_x u_0 - \varepsilon_0^p \right) \right\rangle \dot{\lambda}_0 \, dx \, dt \qquad (60) \\ + \alpha_K (K - K_0) $$

$$ \frac{\partial \mathcal{J}_0}{\partial n} = - \int_0^T \int_0^L \left\langle \left( \frac{E|\partial_x u_0 - \varepsilon_0^p|}{K} \right)^n \text{sign} \left( \partial_x u_0 - \varepsilon_0^p \right) \log \left( \frac{E|\partial_x u_0 - \varepsilon_0^p|}{K} \right) \right\rangle \dot{\lambda}_0 \, dx \, dt \qquad (61) \\ + \alpha_n (n - n_0) $$

Once again, these relations correspond to the zeroth-order time-homogenized estimates of the two misfit function's gradient components (23)–(24) obtained for the identification problem associated with the reference problem. Actually, this is a result that already occurs in periodic space homogenization, as shown for example in Ref. [23].

## Results and discussion

### Validation of the time homogenization method

The studied bar, which is 1 m long, is discretized in ten quadratic finite elements. The surface force at $x = L$ is a combined cycle load, defined as the sum of two sines of frequencies $F$ and $F/\xi$. For this case of study, the assumption (26) is verified, since the two loading frequencies are $F = 0.05$ Hz and $F/\xi = 500$ Hz (cyclic loads with respective
amplitudes 1 kN and 0.25 kN) and the material considered is a typical steel (E = 200 GPa and ρ = 7, 800 kg m⁻³). Moreover, the damping ratio cK is equal to 10⁻⁵, hence verifying assumption (27).

The reference calculations consist in solving the Eqs. (11)–(15) of the reference problem for a time interval corresponding to the first slow loading period [0, 20]s: the time step is chosen according to the classical rule of thumb of 20 time steps per fast loading period, i.e. 10⁻⁴ s. MATLAB's 'ode45' procedure is used.

The zeroth-order time-homogenized equations (37)–(41) are solved with the same algorithm, but the chosen time step is 0.01 s, which reduces by 100 the number of iterations required for the calculation. Since the loading is the sum of two periodic components (one slow and one fast), the residual problem, consisting in Eqs. (43)–(46), depends on the fast time scale only, and can be solved in the frequency domain once and for all.

The associated results are compared in Figure 6 presenting how evolves, for the first slow period, the zeroth-order longitudinal plastic strain at $x = 0$ and $x = L$: the estimates are in excellent agreement, as confirmed by Table 1, and tend to indicate that the predictions can remain good even if a high number of slow cycles is applied. Figure 7 shows the previous time evolutions (at $x = 0$) between two close time steps, and allows us to see that the zeroth-order longitudinal plastic strain evolves smoothly, while the reference strain increases step by step for each fast loading period. All these results allow to validate with this typical example the periodic time homogenization method, and to get an estimate of the computational gain it allows through the drastic reduction in the

<figure><img src="image_10.png" /><figcaption>**Figure 6** Validation of the periodic time homogenization method. Reference (in red) and zeroth-order time-homogenized (in blue) estimates, at $x = 0$ and $x = L$, of the longitudinal plastic strain.</figcaption></figure>

<table><caption>**Table 1** Validation of the periodic time homogenization method</caption><thead><tr><th></th><th>x = 0</th><th>x = L</th></tr></thead><tbody><tr><td>Reference</td><td>1.874891 × 10<sup>-6</sup></td><td>1.652784 × 10<sup>-6</sup></td></tr><tr><td>Time-homogenized</td><td>1.874898 × 10<sup>-6</sup></td><td>1.652791 × 10<sup>-6</sup></td></tr></tbody></table>

Reference and zeroth-order time-homogenized estimates, at $t = 20$ s, of the longitudinal plastic strain at $x = 0$ and $x = L$.
<figure><img src="image_11.png" /><figcaption>**Figure 7** Validation of the periodic time homogenization method. Zoom of Figure 6 (same color conventions).</figcaption></figure>

number of required time steps for the numerical simulation. This number is indeed the main criterion concerning the computational cost of the time homogenization method, since, when compared with the reference problem, the only additional calculations consist in solving the fast residual system (once and for all) and in estimating the fast-time averages with the trapezoidal rule (47).

## Parameter identification

In order to evaluate how the misfit function performs, synthetic data $u_{exp}(t)$ are created by solving the Eqs. (11)–(15) of the reference forward problem, using $K_{exp} = 100 \times 10^6$ SI Units and $n_{exp} = 10$ as parameter values. As previously, the bar withstands a combined cycle loading with $F = 0.05$ Hz and $F/\xi = 500$ Hz. $K_0 = 50 \times 10^6$ SI Units and $n_0 = 5$ are chosen as initial parameter values for the identification process, which is based on an interior-reflective Newton method [24] with a BFGS formula in order to minimize the misfit function $J_0$. Since we do not want to address specifically the question of regularization, no regularization term is added to the misfit function (48). Eventually, given the fact that the two parameters have very different orders of magnitude, a normalization is introduced, using the values of $K_0$ and $n_0$, in order to prevent numerical inaccuracies in the gradient estimate. The logarithm of the misfit function is depicted in Figure 8, showing a well-defined global minimum despite a crescent-shaped valley.

Subsequent results are listed in Table 2: the parameter values identified are very close to the ones used to create the synthetic data. Figure 9 shows the comparison between the identified model and the synthetic reference, more precisely the variations of the longitudinal plastic strain, which is not directly observable.

In addition, to demonstrate the robustness of the process, centered Gaussian noise is added to the synthetic data used in the identification strategy: different levels (corresponding to a standard deviation proportional to the mean of the measured
<figure><img src="image_12.png" /><figcaption>**Table 2 Parameter identification process**</figcaption></figure>

<table><thead><tr><td>Noise level</td><td>Number of iterations</td><td>K<sub>id</sub> (SIUnits)</td><td>n<sub>id</sub></td></tr></thead><tbody><tr><td>0</td><td>24</td><td>100.89 × 10<sup>-6</sup></td><td>9.94</td></tr><tr><td>10%</td><td>23</td><td>100.96 × 10<sup>-6</sup></td><td>9.95</td></tr><tr><td>20%</td><td>22</td><td>100.87 × 10<sup>-6</sup></td><td>9.94</td></tr><tr><td>50%</td><td>37</td><td>99.40 × 10<sup>-6</sup></td><td>9.78</td></tr><tr><td>100%</td><td>35</td><td>97.36 × 10<sup>-6</sup></td><td>9.56</td></tr></tbody></table>

Results obtained for synthetic data (with $K_{exp} = 100 \times 10^6$ SI Units and $n_{exp} = 10$) and different noise levels.

displacement) are proposed, and the associated identification results are listed in Table 2. They show that the identification works, even for high levels of noise: the choice of the misfit function (48) with the application of the fast-time average to the noisy data has allowed to filter significantly the introduced noise, hence strongly limiting its impact on the identification results. This is a particularly interesting property, whereas identifying using the reference model and the whole set of noisy data undoubtedly would have led to poor results: indeed, as said above, the fact of using a homogenized model in a parameter identification process generally introduces a regularizing effect, as explained in Ref. [17]; unfortunately, even in this academic case of study, the numerical cost for solving the reference identification problem is too prohibitive to rigorously check this property.

Eventually, the computational cost associated with the identification process is significantly reduced when compared with what is obtained when the inverse problem related
<figure><img src="image_13.png" /><figcaption>**Figure 9** Parameter identification of a time-homogenized model. Experimental (in red) and predicted (in blue) longitudinal plastic strain.</figcaption></figure>

to the reference problem (11)–(15) is considered, which allows to address the case of slow-evolving phenomena, requiring the simulation over a large time interval.

## Conclusions

Here we have proposed a first preliminary study of a two-time-scale parameter identification process, using time-homogenized models: the adaptation of a classical identification strategy based on an adjoint state formulation to estimate the misfit function's gradient can be used in this specific framework: this leads, on the proposed example, to the determination of the time-homogenized counterpart of the adjoint solution associated with the reference identification problem. Despite its simplicity, the academic example studied here showed the relevance of the strategy and its reduced computational cost: these results can be viewed as a first step before dealing with more complex cases of study.

To come back more specifically to CCF life prediction, a possible prospect could be to use the experimental information available (with classical high cycle fatigue tests for instance) to tune correctly the parameters of the considered material laws. In addition, to further improve the reduction in the computational cost, the extension to models with three different time scales, such as those described in Ref. [5], should be straightforward as well.

## Authors' contributions

DA started the project and the initial coding. GP developed the current version of the code, analysed the data and wrote the manuscript. Both collaborated to develop the analyses and expertise on the different aspects of the model. Both authors read and approved the final manuscript.

## Acknowledgements

The authors would like to acknowledge the support of the European Commission through the project PREMECCY, described in [15], which allowed to initiate this study.

## Compliance with ethical guidelines

### Competing interests

The authors declare that they have no competing interests.
## References

1. Fatemi A, Yang L (1998) Cumulative fatigue damage and life prediction theories: a survey of the state of the art for homogeneous materials. Int J Fatigue 20(1):9–34

2. Guennouni T, Aubry D (1986) Réponse homogénéisée en temps de structures sous chargements cycliques. Comptes rendus de l'Académie des sciences. Série II 303(20):1765–1768 (in French)

3. Bensoussan A, Lions JL, Papanicolaou G (1978) Asymptotic analysis for periodic structures. Elsevier, Burlington

4. Sanchez-Palencia E (1980) Non-homogeneous media and vibration theory. Springer, Berlin

5. Puel G, Aubry D (2014) Efficient fatigue simulation using periodic homogenization with multiple time scales. Int J Multiscale Comput Eng 12(4):291–318

6. Yu Q, Fish J (2002) Temporal homogenization of viscoelastic and viscoplastic solids subjected to locally periodic loading. Comput Mech 29(3):199–211

7. Shishkina EV, Blekhman II, Cartmell MP, Gavrilov SN (2008) Application of the method of direct separation of motions to the parametric stabilization of an elastic wire. Nonlinear Dyn 54(4):313–331

8. Lévy A, Le Corre S, Poitou A, Soccard E (2011) Ultrasonic welding of thermoplastic composites: modeling of the process using time homogenization. Int J Multiscale Comput Eng 9(1):53–72

9. Dontsov EV, Guzina BB (2012) Acoustic radiation force in tissue-like solids due to modulated sound field. J Mech Phys Solids 60(10):1791–1813

10. Papon A, Yin Z-Y, Riou Y, Hicher P-Y (2013) Time homogenization for clays subjected to large numbers of cycles. Int J Numer Anal Methods Geomech 37(11):1470–1491

11. Haouala S, Doghri I (2015) Modeling and algorithms for two-scale time homogenization of viscoelastic-viscoplastic solids under large numbers of cycles. Int J Plast 70:98–125

12. Aubry D, Puel G (2010) CCF modelling with use of a two-timescale homogenization model. In: Proceedings of the international fatigue conference (Fatigue 2010), vol 2. Procedia Engineering, pp 787–796

13. Devulder A, Aubry D, Puel G (2010) Two-time scale fatigue modelling: application to damage. Comput Mech 45(6):637–646

14. Puel G, Aubry D (2012) Material fatigue simulation using a periodic time homogenization method. Eur J Comput Mech 21(3–6):312–324

15. PREMECCY (2012) Predictive methods for combined cycle fatigue in gas turbine blades. Final Project Report AST5-CT-2006-030889, European Commission—6th RTD Framework Program

16. Puel G, Bourgeteau B, Aubry D (2013) Parameter identification of transient nonlinear models for the multibody simulation of a vehicle chassis. Mech Syst Signal Process 36(2):354–369

17. Chung C-B, Kravaris C (1988) Identification of spatially discontinuous parameters in second-order parabolic systems by piecewise regularisation. Inverse Probl 4(4):973–994

18. Constantinescu A, Tardieu N (2001) On the identification of elastoviscoplastic constitutive laws from indentation tests. Inverse Probl Eng 9(1):19–44

19. Plessix R-E (2006) A review of the adjoint-state method for computing the gradient of a functional with geophysical applications. Geophys J Int 167(2):495–503

20. Papon A (2010) Numerical modeling of soil behavior under very large numbers of cycles: time homogenization and parameter identification. Ph.D. thesis, Ecole Centrale de Nantes (in French)

21. Dormand JR, Prince PJ (1980) A family of embedded Runge-Kutta formulae. J Comput Appl Math 6(1):19–26

22. Weideman JAC (2002) Numerical integration of periodic functions: a few examples. Am Math Mon 109(1):21–36

23. Kesavan S, Saint Jean Paulin J (1997) Homogenization of an optimal control problem. SIAM J Control Optim 35(5):1557–1573

24. Coleman TF, Li Y (1994) On the convergence of interior-reflective Newton methods for nonlinear minimization subject to bounds. Math Progr 67(1–3):189–224

Submit your manuscript to a SpringerOpen journal and benefit from:

* Convenient online submission

* Rigorous peer review

* Immediate publication on acceptance

* Open access: articles freely available online

* High visibility within the field

* Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com