2. Benchmark solution#
2.1. Calculation of the reference solution (cf. bib [1] and [3])#
For a tensile test in direction \(x\), the Kirchhoff \(\tau\) and Cauchy \(\sigma\) tensors are of the form:
and
with
The change in volume
is given by the resolution of
where
is thermal deformation. This applies to an austenitic — bainitic transformation:
Note:
The K coefficient is the compression modulus (not to be confused with the coefficients)
relating to the law of plasticity of transformation)
In plastic filler, for isotropic work hardening
linear, such as:
the cumulative plastic deformation
Worth
with
The gradient tensors of the transformation
and
and the plastic deformation tensor
are of the form:
The law of evolution of plastic deformation
is written:
For
, we have
. There is no transformation plasticity. We then obtain:
For
, we have
.
To integrate the law of evolution of plastic deformation, it must be assumed that the Kirchhoff stress \(\tau\) varies very little, that is to say that the change in volume
is very small. Under this hypothesis, we obtain
The component
of the transformation gradient is given by the resolution of:
Finally, the displacement field \(u\) (in the initial configuration) is of the form
. The components are given by:
2.2. note#
In test case HSNV101 (modeling B), the coefficients of the material were chosen in such a way as to have no classical plasticity
during the metallurgical transformation that takes place between moments \(60\) and \(\mathrm{122s}\). In fact, if we write the charge—discharge criterion in this time interval, we obtain
with
which is cancelled out for only one value of the cumulative plastic deformation
.
For the law of behavior written in large deformations, the charge—discharge criterion is written between these two moments
with
In this case, as long as the variable
remains less than the value obtained in time \(t\mathrm{=}\mathrm{60s}\), we will have
. However, the value of
is a function only of the thermal deformation value (stress \(\sigma\) is constant) and the coefficient
is independent of metallurgical phases and temperature).
In this time interval, thermal deformation
is given by the following equation:
The thermal deformation as well as the volume variation are plotted below.
, solution of the 3rd degree equation, as a function of time.
Thermal deformation as a function of time
Volume variation
depending on the time
It can be seen that the variable
decreases and increases in the same way as thermal deformation. In this case, to find out the time from which the variable
is greater than the value obtained at time \(\mathrm{60s}\), it is sufficient to know the instant for which the thermal deformation is identical to that obtained at time \(t\mathrm{=}\mathrm{60s}\). By solving the equation above, we find \(t\mathrm{=}\mathrm{84.46s}\).
2.3. Uncertainty about the solution#
The solution is analytical. Two mistakes are made with this solution. The first relates to the calculation of the proportion of bainitic phase created. The preliminary metallurgical calculation does not exactly restore the equation of [§1.2] giving
as a function of time, which is why the reference results presented below are calculated with the proportion of bainitic phase calculated by Code_Aster.
The second error is the assumption made on the Kirchhoff constraint \(\tau\) which is not constant over the time interval between \(60\) and \(\mathrm{176s}\). This will impact the displacement calculation.
and plastic deformation
.
2.4. Benchmark results#
As reference results, we will adopt the displacement in the direction of tensile loading, the Cauchy stress \(\sigma\), the Boolean plasticity indicator \(\chi\) and the cumulative plastic deformation.
. The different calculation times are \(t\mathrm{=}\mathrm{47,}\mathrm{48,}\mathrm{60,}\mathrm{83,}\mathrm{84,}85\) and \(\mathrm{176s}\). For the calculation of the displacement, the initial length of the bar in the loading direction is \(\mathrm{0.2m}\).
In any case, we have
\(\mathrm{3K}\mathrm{=}500000\mathit{MPa}\) (compression module) \(\mu \mathrm{=}76923.077\mathit{MPa}\)
At time \(t\mathrm{=}\mathrm{47s}\), we have
,
,
At time \(t\mathrm{=}\mathrm{48s}\), we have
,
,
At time \(t\mathrm{=}\mathrm{60s}\), we have
,
,
At time \(t\mathrm{=}\mathrm{83s}\), we have
,
,
At time \(t\mathrm{=}84\) s, we have
,
,
At time \(t\mathrm{=}\mathrm{85s}\), we have
,
,
At time \(t\mathrm{=}\mathrm{176s}\), we have
,
,
2.5. Bibliographical references#
Reference may be made to:
CANO, E. LORENTZ: Introduction to the Code_Aster of a model of behavior in large elastoplastic deformations with isotropic work hardening — Internal note EDF DER HI‑74/98/006/0
A.M. DONORE, F. WAECKEL: Influence of structural transformations in elasto—plastic behavior laws Note HI—74/93/024
WAECKEL, V. CANO: Law of behavior of large elasto (visco) plastic deformations with metallurgical transformations [R4.04.03]