Effect of Load Path on Parameter Identification for Plasticity Models using Bayesian Methods

by   Ehsan Adeli, et al.

To evaluate the cyclic behavior under different loading conditions using the kinematic and isotropic hardening theory of steel, a Chaboche viscoplastic material model is employed. The parameters of a constitutive model are usually identified by minimization of the distance between model response and experimental data. However, measurement errors and differences in the specimens lead to deviations in the determined parameters. In this article, the Choboche model is used and a stochastic simulation technique is applied to generate artificial data which exhibit the same stochastic behavior as experimental data. Then the model parameters are identified by applying an estimation using Bayes's theorem. The Gauss-Markov-Kalman filter using functional approximation is introduced and employed to estimate the model parameters in the Bayesian setting. Identified parameters are compared with the true parameters in the simulation, and the efficiency of the identification method is discussed. In the end, the effect of the load path on parameter identification is investigated.



There are no comments yet.


page 1

page 2

page 3

page 4


Parameter Identification in Viscoplasticity using Transitional Markov Chain Monte Carlo Method

To evaluate the cyclic behavior under different loading conditions using...

Parameter identification in elasto-plasticity: distance between parameters and impact of measurement errors

A special aspect of parameter identification in finite-strain elasto-pla...

Inspection of ratcheting models for pathological error sensitivity and overparametrization

Accurate analysis of plastic strain accumulation under stress-controlled...

Real-time parameter inference in reduced-order flame models with heteroscedastic Bayesian neural network ensembles

The estimation of model parameters with uncertainties from observed data...

Bayesian stochastic multi-scale analysis via energy considerations

In this paper physical multi-scale processes governed by their own princ...

Parameters identification of photovoltaic cell models using enhanced exploratory salp chains-based approach

The integration of photovoltaic systems (PVSs) in future power systems g...

Parameter Identification for Multirobot Systems Using Optimization Based Controllers (Extended Version)

This paper considers the problem of parameter identification for a multi...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

In order to predict the behavior of loaded metallic materials, constitutive models are applied, which present a mathematical frame for the description of elastic and inelastic deformation. The models by Miller, Krempl, Korhonen, Aubertin, Chan, and Bodner are well-known constitutive models for isotropic materials [1, 2, 3, 4, 5]. In 1983, Chaboche [6, 7] put forward what has become known as the unified Chaboche viscoplasticity constitutive model, which has been widely accepted.

All inelastic constitutive models contain parameters which have to be identified for a given material from experiments. In the literature only few investigations can be found dealing with identification problems using stochastic approaches. Klosowski and Mleczek have applied the least-squares method in the Marquardt-Levenberg variant to estimate the parameters of an inelastic model [8]. Gong et al. have also used some modification of the least-squares method to identify the parameters [9]. Harth and Lehn identified the model parameters of a model by employing some generated artificial data instead of experimental data using a stochastic technique [10]. A similar study by Harth and Lehn has been done for other constitutive models like Lindholm and Chan [11].

In this paper, a viscoplastic model of Chaboche is studied. The model contains five material parameters which have to be determined from experimental data. It should be noted that here virtual data are employed instead of real experimental data. A cyclic tension-compression test is applied in order to extract the virtual data.

The model is described in Section 2, whereas Section 3 explains how to propagate the uncertainty in the model and how to perform the update. The probabilistic model is reformulated from the deterministic model, and once the forward model is provided, the model parameters are updated using a Bayesian approach.

In Section 4 the desired parameters are identified from the measured data. In fact, the parameters which have been considered as uncertain parameters are updated and the uncertainties of the them are reduced while the random variables representing the uncertain parameters are updated during the process. The results are thoroughly studied and the identified parameters as well as the corresponding model responses are analyzed. Finally the prediction of the models is compared with the measured data for different applied load paths. It is also explained why different load paths cause different identification of model parameters.

2 Model problem

The mathematical description of metals under cyclic loading beyond the yield limit that includes viscoplastic material behavior as well as the characterization of compulsory isotropic-kinematic hardening is here given in terms of a modified Chaboche model introduced in [12]. As we consider classical infinitesimal strains, we assume an additive strain decomposition. The material behavior is described for the elastic part by isotropic homogeneous elasticity, and for viscoplasticity the dissipation potential is given by


with and and as the material parameters. Here is the over-stress, defined via the equivalent stress () which reads


where denotes the deviatoric part and is the back-stress of kinematic hardening. The over-stress is given by


where is the yield stress and models the isotropic hardening which is introduced in the following. The partial derivative of the dissipation potential with respect to leads to the equation for the inelastic strain rate


The viscoplastic model allows for isotropic and kinematic hardening, which is considered in order to describe different specifications. Assuming and with and to describe isotropic and kinematic hardening respectively, the evaluation equations for these two are




respectively. In the evaluation equations of the both hardening, is the viscoplastic multiplier rate given as:


which describes the rate of accumulated plastic strains. The parameter indicates the speed of stabilization, whereas the value of the parameter is an asymptotic value according to the evolution of the isotropic hardening. Similarly, the parameter denotes the speed of saturation and the parameter is the asymptotic value of the kinematic hardening variables. The complete model is stated in Table 1. Note that

represents the elasticity tensor.

Strain Hooke’s Law Flow Rule Hardening

Initial Conditions
,  ,  
 (Yield Stress) ,  (Flow Rule) , , ,  (Hardening)
Table 1: The constitutive model of Chaboche

By gathering all the desired material parameters to identify into the vector

, where and are bulk and shear modulus, respectively, which determine the isotropic elasticity tensor, the goal is to estimate given measurement displacement data, i.e.


in which represents the measurement operator and the measurement (also possibly the model) error. Being an ill-posed problem, the estimation of given is not an easy task and usually requires regularization. This can be achieved either in a deterministic or a probabilistic setting. Here, the latter one is taken into consideration as further described in the text.

3 Bayesian identification

By using additional (prior) knowledge on the parameter set next to the observation data, the probabilistic approach regularizes the problem of estimating with the help of Bayes’s theorem


in which the likelihood describes how likely the measurement data are given prior knowledge . This in turn requires the reformulation of the deterministic model into a probabilistic one, and hence the propagation of material uncertainties through the model —the so-called forward problem— in order to obtain the likelihood [13, 14].

The main difficulty in using Equation (9

) lies in the computation of the likelihood. Various numerical algorithms can be applied, the most popular example of which are the Markov chain Monte Carlo methods. Being constructed on the fundamentals of ergodic Markov theory, these methods are characterized by very slow convergence. To avoid this, an approximate method based on Kolmogorov’s definition of conditional expectation as already presented in

[15] is considered here.

Let the material parameters

be modeled as random variables on a probability space

. Here, denotes the space of elementary events , is the -algebra and stands for the probability measure. This alternative formulation of Bayes’s rule can be achieved by expressing the conditional probabilities in Equation (9) in terms of conditional expectation. Following the mathematical derivation in [15, 17, 16, 18], this approach boils down to a quadratic minimization problem by considering the forecast random variable and the update of the forecast random variable :


where is the orthogonal projection operator of onto the space of the new information in which the space is the space of random variables generated by the measurement . Due to the Doob-Dynkin lemma, is a function of the observation, where is the forecast, and the assimilated value is .

Constraining the space of all functions to the subspace of linear maps, the minimization problem in Equation (10) leads to a unique solution . Note that the projection is performed over a smaller space than . An implication of this is that available information is not completely used in the process of updating, introducing an approximation error. This gives an affine approximation of Equation (10)


also known as a linear Bayesian posterior estimate or the so-called Gauss-Markov-Kalman filter (GMKF). Here, represents the prior random variable, is the posterior approximation, is the predicted measurement and represents the very well-known Kalman gain


which can be easily evaluated if the appropriate covariance matrices , and are known.

An advantage of Equation (11) compared to Equation (9) is that the inference in Equation (11) is given in terms of random variables instead of conditional densities. Namely, , , and denote the random variables used to model the posterior, prior, observation, and predicted observation, respectively.

In this light the linear Bayesian procedure can be reduced to a simple algebraic method. Starting from the functional representation of the prior


where are multivariate Hermite polynomials, and by considering the proxy in Equation (13), one may discretize Equation (11) as:


where , etc. are the PCE coefficients. As the measurement is a deterministic value, has only a zero-th order tensor. The covariances for the Kalman gain Equation 12 are easily computed, e.g.


4 Numerical results

The identification of the material constants in the Chaboche unified viscoplasticity model is a reverse process, here based on virtual data. In case of the Chaboche model the best way of parameter identification is using the results of the cyclic tests, since more information can be obtained from cyclic test rather than creep and relaxation tests, specifically information regarding hardening parameters. The aim of the parameter identification is to find a parameter vector introduced in the previous section. The bulk modulus (), the shear modulus (), the isotropic hardening coefficient (), the kinematic hardening coefficient () and the yield stress () are considered as the uncertain parameters of the constitutive model.

A preliminary study is on a regular cube, modeled with one 8 node element, completely restrained on the back face, and with normal traction on the opposite (front) face. Two cases are considered in order to compare the effect of applied force on identified parameters. For both cases the magnitude of the normal traction and a stress in the plane of the front face are plotted in Fig. 1 and Fig. 2, respectively. Blue and red colors represent the stress value in normal and in plane directions, respectively. As it is seen, the magnitude of the applied force for the case 1 is constant all time but for the case 2 the magnitude of the applied force grows gradually by time.

Figure 1: Decomposed applied force on desired node according to time- Case 1
Figure 2: Decomposed applied force on desired node according to time- Case 2

Considering the parameters listed in Table 2, the related - hysteretic graph obtained for the applied force case 1 and 2 which can be seen in Fig. 3 and Fig. 4, respectively.

Table 2: The model parameters

Figure 3: - for node on the front surface in plane and normal directions- Case 1

Figure 4: - for node on the front surface in plane and normal directions- Case 2

The displacements of one of the nodes on the front surface in normal and in plane directions are observed as the virtual data in this study. Applying the Gauss-Markov-Kalman filter with functional approximation as explained in the previous chapter and introducing measurement error in such a way that 15 percent of mean values are equal to the coefficient of variation for the related parameter, the probability density function (PDF) of prior and posterior of the identified parameters can be seen in Fig. 

5 and Fig. 6 for the first and second case, respectively.

Figure 5: PDF of identified parameters- Case 1

Figure 6: PDF of identified parameters- case 2

Summarizing the results, the true values

and the mean and standard deviation of the estimated parameters,

and respectively, for both cases are compared in Table 3.

50 52.36 3.71 50.27 0.29
50 52.04 3.01 50.19 0.53
Table 3: The identified model parameters

4.1 Discussion of the results

From the sharpness of the posterior PDF of , and , it can be concluded that enough information from virtual data is received and updating the parameters considering their uncertainty is done very properly for the both cases, as the standard deviation of the residual uncertainty is below of the mean.

For the posterior PDF of and , it can be inferred that better updating is done for the second case compared to the first case. Not only are the more accurate estimations of the exact hardening parameters, and , predicted for the second case, but the uncertainty of the estimated hardening parameters are also reduced much more for the second case.

One reason that can be mentioned is that the process is not always in the states that hardening equations are involved like the elastic states. Therefore less information from the whole simulation can be analyzed for estimating the hardening parameters and updating their parameters’ uncertainties. Fig. 7 and Fig. 8 prove this fact that since more states are out of the von Mises yield criterion for the second case compared to the first case, in which the hardening equations are involved only in these states, the better identification can be done for the second case, where a gradually varying increasing applied force is considered, for hardening parameters in comparison with the first case where a constant magnitude applied force is employed. In fact, the cyclic applied force in the second case causes the more activation of the desired parameters in the studied set of equations comparing to the first case and accordingly a better determination of the parameters can be carried out for the second case. It should be pointed out that the von Mises yield criterion is illustrated by the green cylinder in the mentioned figures i.e. inside and outside of the cylinder respectively refer to the elasticity and plasticity states, respectively. Also the blue color represents the principal stresses in these figures.

Figure 7: Principal stresses of applied force in 3D considering the von Mises yield criterion- Case 1

Figure 8: Principal stresses of applied force in 3D considering the von Mises yield criterion- Case 2

5 Summary

Using the Gauss-Markov-Kalman Filter method explained in Section 3 to identify the model parameters of the Choboche model indicates that it is possible to identify the model parameters by employing this method using functional approximation. The parameters are well estimated and the uncertainty of the parameters is reduced while the random variables of the parameters are updated during the process [20]

. The other conclusion that can be made is that the more information we receive, the better parameter identification we can do using the Gauss-Markov-Kalman filter method. This fact is observed by comparing the posterior probability density functions of hardening parameters,

and , for case 1 and 2. Therefore in terms of mechanical models, it should be always considered that the applied force should be applied in such a way that all time all equations should be involved. In other words, the applied load path should lead to activation of all uncertain parameters in the set of equations, as here a cyclic gradually varying increasing applied force leads to a better determination of the parameters. Otherwise only the parameters which are in the involved equations are updated.

This work is partially supported by the DFG through GRK 2075.


  • [1] A. Miller. An Inelastic Constitutive Model for Monotonic, Cyclic, and Creep Deformation: Part I–Equations Development and Analytical Procedures. J. Eng. Mater. Technol. 98(2), 97–105 (1976).
  • [2] E. Krempl, J. J. McMahon, and D. Yao. Viscoplasticity Based on Overstress with a Differential Growth Law for the Equilibrium Stress. Mechanics of Materials 5, 35–48 (1986).
  • [3] R. K. Korhonen, M. S. Laasanen, J. Toyras, R. Lappalainen, H. J. Helminen, and J. S. Jurvelin. Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal, proteoglycan depleted and collagen degraded articular cartilage. J Biomech 36, 1373–1379.
  • [4] Michel Aubertin, Denis E. Gill, and Branko Ladanyi. A unified viscoplastic model for the inelastic flow of alkali halides. Mechanics of Materials 11, 63–82 (1991).
  • [5] K. S. Chan, S. R. Bodner, A.F. Fossum, and D.E. Munson. A constitutive model for inelastic flow and damage evolution in solids under triaxial compression. Mechanics of Materials 14, 1–14 (1992).
  • [6] J. L. Chaboche, and G. Rousselier. On the plastic and viscoplastic constitutive equations - part 1: rules developed with internal variable concept. J. Press. Vessel Technol., 105, 153–158 (1983).
  • [7] J. L. Chaboche, and G. Rousselier. On the plastic and viscoplastic constitutive equations - part 2: application of internal variable concepts to the 316 stainless steel. J. Press. Vessel Technol., 105, 159–164 (1983).
  • [8] P. Kłosowski, and A. Mleczek. Parameters’ Identification of Perzyna and Chaboche Viscoplastic Models for Aluminum Alloy at Temperature of . Engng. Trans. 62, 3, 291–305 (2014).
  • [9] Y. Gong, C. Hyde, W. Sun, and T. Hyde. Determination of material properties in the Chaboche unified viscoplasticity model. Journal of Materials Design and Applications, 224(1), 19–29 (2010).
  • [10] T. Harth, and Jürgen Lehn. Identification of Material Parameters for Inelastic Constitutive Models Using Stochastic Methods. GAMM-Mitt. 30, No. 2, 409–429 (2007).
  • [11] K. S. Chan, S. R. Bodner, and U. S. Lindholm. Phenomenological Modelling of Hardening and Thermal Recovery in Metals. Journal of Engineering Materials and Technology 110, 1–8 (1988).
  • [12] J. Velde. 3D Nonlocal Damage Modeling for Steel Structures under Earthquake Loading. Department of Architecture, Civil Engineering and Environmental Sciences University of Braunschweig - Institute of Technology (2010).
  • [13]

    H. G. Matthies. Stochastic finite elements: Computational Approaches to Stochastic Partial Differential Equations. Journal of Journal of Applied Mathematics and Mechanics, volume 88, 849–873 (2008).

  • [14] H. G. Matthies. Uncertainty Quantification with Stochastic Finite Elements. Encyclopedia of Computational Mechanics, edited by E. Stein, R. de Borst, T. R. J. Hughes, John Wiley and Sons, Chichester, (2007).
  • [15] H. G. Matthies, E. Zander, B. V. Rosić, and A. Litvinenko. Parameter Estimation via Conditional Expectation: A Bayesian Inversion. Journal of Advanced Modeling and Simulation in Engineering Sciences, 3–24 (2016).
  • [16] B. V. Rosić, A. Litvinenko, O. Pajonk, H. G. Matthies. Sampling-free linear Bayesian update of polynomial chaos representations. Journal of Computational Physics 231(17):5761–5787, (2012).
  • [17] H. G. Matthies, E. Zander, B. V. Rosić, A. Litvinenko, and Oliver Pajonk. Inverse Problems in a Bayesian Setting. Journal of Computational Methods for Solids and Fluids, volume 41, 245–286 (2016).
  • [18] B. V. Rosić and H. G. Matthies. Identification of Properties of Stochastic Elastoplastic Systems. Computational Methods in Stochastic Dynamics, 237–253, (2013).
  • [19] S. Ghahramani. Fundamentals of probability, with stochastic processes. 3rd ed. New Jersey, USA: Pearson, Prentice Hall (2005).
  • [20] E. Adeli, B. V. Rosić, H. G. Matthies and S. Reinstädler. Bayesian parameter identification in plasticity. XIV International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIV E. Oñate, D.R.J. Owen, D. Peric and M. Chiumenti (Eds), 247–256, DOI: 10.13140/RG.2.2.22323.89124, (2017).