I am beginner in crystal plasticity FEM .I need to update Huang's UMAT subroutine with Bassini and Wu hardening parameters. Can anyone guide in this basic query.
Thank you in advance.
Huang’s UMAT allows to input Bassani and Wu’s self-hardening parameters. The .inp files for the subroutine are defined with 8 properties per line. You need to write the corresponding modelling parameters from line 13 to line 18 (marked in bold).
*Material, name=MATERIAL-1
*Depvar
125,
*User Material, constants=160, unsymm
168400.,121400., 75400., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
1., 0., 0., 0., 0., 0., 0., 0.
1., 1., 1., 1., 1., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
-1., 0., 1., 0., 0., 1., 0., 0.
0., 1., 0., 0., 1., 0., 0., 0.
10., 0.001, 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
541.5, 109.5, 60.8, 0., 0., 0., 0., 0.
1., 1., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0.5, 1., 0., 0., 0., 0., 0., 0.
1., 10., 1e-05, 0., 0., 0., 0., 0.
Huang’s UMAT can generate up to three sets of slip systems for a cubic crystal. Lines 13 and 14 correspond to the hardening parameters of the first set of slip systems. In line 13 you must input Bassani and Wu’s parameters in the following order:
h0, τs, τ0, hs, γ0, γ1, f0, f1
In the example above, replace the values of line 13 (541.5, 109.5, 60.8, 0., 0., 0., 0., 0.) with the particular values you want to use.
Line 14 is used for the latent hardening parameters:
q, q1
If you have additional sets of slip systems, you can use the same logic to complete the remaining lines.
Mauro Francisco Arcidiacono Thank you so much for your reply.
In addition to this, the literature I am following used five (apart from h0, τs, τ0, hs, γ0, γ1, f0, f1) Bassani and Wu hardening law namely, a1 (no junction)=1.75, a2 (Hirth lock=1.75), a3 (coplanar junction=1.75), a4 (glissile junction=2) and a5 (sessile junction=2.25).
Can you please help me here, how to add these five parameters in the Huang script?
To add those five parameters you have to modify the functions HSELF and HLATNT. They are called in the LATENTHARDEN subroutine. You will need to pass the parameters through the PROPS array and adjust the code in the functions so it performs the calculations you require.
Mauro Francisco Arcidiacono Thank you very much. Being a beginner in this will be challenging for me.
- Badges
- Science topic
- Similar topics
- Filing