2. Benchmark solution#
2.1. Calculation method used for the reference solution#
It is identical to that of test SSNP15.
In plan \((\sigma ,\tau \sqrt{3})\), the vonMises standard is translated into the classical distance, so that we can immediately predict the charging and discharging phases during the loading path, since they are respectively the phases in which the norm increases or decreases:
2.1.1. Resolution process#
Mechanically, it is a \(\mathrm{0D}\) stress-controlled test, the material being elastoplastic with VonMises criterion and linear isotropic work hardening. For stress-controlled loading, the cumulative plastic deformation is easily determined:
\(F(\sigma ,p)={\sigma }_{\mathrm{eq}}-{\sigma }^{Y}-R\text{'}p\le 0\Rightarrow p=\frac{{\sigma }_{\mathrm{eq}}-{\sigma }^{Y}}{R\text{'}}\) in charge eq 2.1.1-2.1.1-1
The integration of plastic deformation is of course more difficult. The flow equation is written as:
\(\dot{{\varepsilon }^{p}}=\frac{3}{2}\dot{p}\frac{\tilde{\sigma }}{{\sigma }_{\mathrm{eq}}}\Rightarrow \dot{{\varepsilon }^{p}}=\frac{3}{2R\text{'}}\frac{\dot{{\sigma }_{\mathrm{eq}}}}{{\sigma }_{\mathrm{eq}}}\tilde{\sigma }\) in charge eq 2.1.1-2.1.1-2
Finally, we will deduce the deformation via the state relationship:
\(\varepsilon ={\varepsilon }^{p}+{E}^{-1}:\sigma \Rightarrow {\varepsilon }_{\mathrm{xx}}={\varepsilon }_{\mathrm{xx}}^{p}+\frac{\sigma }{E}\) and \({\varepsilon }_{\mathrm{xy}}={\varepsilon }_{\mathrm{xy}}^{p}+\frac{\tau }{2\mu }\) eq 2.1.1-2.1.1-3
2.1.2. Radial loading phase treatment#
Note that during the radial loading phase, the flow law [éq 2.1.2-1] integrates directly:
\({\varepsilon }^{p}=\frac{3}{2}p\frac{\tilde{\sigma }}{{\sigma }_{\mathrm{eq}}}\) eq 2.1.2-1
The cumulative plastic deformation is then given by [éq 2.1.1-1], the plastic deformation by [éq2.1.2-1] and the total deformation by [éq 2.1.1-3]. With:
\(E\) |
|
|
|
|
|
|
We get: |
||||||
\(p(A)\) |
|
|
|
|
|
|
\({\varepsilon }_{\mathrm{xy}}^{p}(A)\) |
|
|
|
2.1.3. Treatment of the non-radial loading phase#
In the non-radial loading phase \(\mathrm{B0}-B\), the stress path can be parameterized by:
\(\sigma (q)={\sigma }^{{B}^{0}}+q\underset{\text{direction fixe}}{\underset{\underbrace{}}{({\sigma }^{B}-{\sigma }^{{B}^{0}})}}\) with \(0\le q\le 1\) eq 2.1.3-1
As the loading path remains confined in the traction-shear plane \((\sigma ,\tau )\), it will be useful to represent the stress state by a complex number:
\(\Sigma =\sigma +i\sqrt{(3)}\tau \Rightarrow {\sigma }_{\mathrm{eq}}=∣(\Sigma )∣\) and \(\Sigma (q)={\Sigma }^{{B}^{0}}+q\underset{\text{direction fixe}}{\underset{\underbrace{}}{({\Sigma }^{{B}^{0}}-{\Sigma }^{B})}}\) eq 2.1.3-2
The integration of the flow law [éq 2.1.1-2], followed by an integration by part, makes it possible to express the plastic deformation:
\(\frac{\mathrm{2R}\text{'}}{3}{\left[{\varepsilon }^{p}\right]}_{0}^{1}=\underset{0}{\overset{1}{\int }}\frac{\dot{{\sigma }_{\mathrm{eq}}}}{{\sigma }_{\mathrm{eq}}}\tilde{\sigma }\mathrm{dq}={\left[\mathrm{ln}({\sigma }_{\mathrm{eq}})\tilde{\sigma }\right]}_{0}^{1}-\frac{1}{2}\underset{{\tilde{\sigma }}^{B}-{\tilde{\sigma }}^{{B}^{0}}}{\underset{\underbrace{}}{\dot{\tilde{\sigma }}}}\underset{0}{\overset{1}{\int }}\mathrm{ln}({\sigma }_{\mathrm{eq}}^{2})\mathrm{dq}\)
The adoption of the complex plan allows an easy calculation of the last integral:
\(\underset{0}{\overset{1}{\int }}\mathrm{ln}({\sigma }_{\mathrm{eq}}^{2})\mathrm{dq}=\underset{0}{\overset{1}{\int }}\mathrm{ln}(\Sigma \overline{\Sigma })\mathrm{dq}=\underset{0}{\overset{1}{\int }}\mathrm{ln}(\Sigma )\mathrm{dq}+\underset{0}{\overset{1}{\int }}\mathrm{ln}(\overline{\Sigma })\mathrm{dq}=2\Re \left[\underset{0}{\overset{1}{\int }}\mathrm{ln}(\Sigma )\mathrm{dq}\right]=2\Re {\left[\frac{\Sigma \mathrm{ln}(\Sigma )-\Sigma }{{\Sigma }^{B}-{\Sigma }^{{B}^{0}}}\right]}_{0}^{1}\)
Finally, the plastic deformation increment on path \(\mathrm{B0}-B\) is equal to:
\({\left[{\varepsilon }^{p}\right]}_{{B}^{0}}^{B}=\frac{3}{2R\text{'}}{\left[\mathrm{ln}({\sigma }_{\mathrm{eq}})\tilde{\sigma }\right]}_{{B}^{0}}^{B}-\frac{3}{2R\text{'}}\Re {\left[\frac{\Sigma \mathrm{ln}(\Sigma )-\Sigma }{{\Sigma }^{B}-{\Sigma }^{{B}^{0}}}\right]}_{{B}^{0}}^{B}({\tilde{\sigma }}^{B}-{\tilde{\sigma }}^{{B}^{0}})\) eq 2.1.3-4
2.2. Benchmark results#
By calculating the cumulative plastic deformation by [éq 2.1.1-1], the plastic deformation by [éq 2.1.3-4], and the total deformation by [éq 2.1.1-3], we get:
\(p(B)\) |
|
|
|
|
|
|
We get: |
\({\varepsilon }_{\mathrm{xy}}^{p}(B)\) |
|
|
|
We will focus on the values of the stresses, deformations and the cumulative plastic deformation at points \(A\) and \(B\) of the loading path.
2.3. Bibliographical references#
French Society of Mechanics. Guide to the validation of structural calculation software packages (VPCS). AFNOR Technical, 1990.