La loi d'élastohystérésis 3D permet de représenter le comportement élastoplastique.
Le comportement est décomposé en 2 parties: une sphérique qui est élastique et une déviatorique qui correspond à la composante d'hystérésis.
Un exemple de loi:
#==========================================================================
LOI_acier_elasto_plastique LOI_ADDITIVE_EN_SIGMA
#------------ hyper part --------
ISOHYPER3DFAVIER3
#----------------------------------------
- K | Qor | muzero | mu_inf |
#----------------------------------------
5.25e5 3. 1000. 800.
#------------ non reversible part -----------
HYSTERESIS_3D
#-------------------------------------------------------
- para de prager mu limite maxi de plasticité |
#-------------------------------------------------------
np= 2. mu= 80000. Qzero= 430 avec_parametres_de_reglage_
type_de_resolution_ 2 cas_kutta_ 5
erreurAbsolue_ 1.e-4 erreurRelative_ 1.e-4
tolerance_coincidence_ 1.e-7 mini_rayon_ 1.e-12
nb_maxInvCoinSurUnPas_ 10 min_angle_trajet_neutre_ 0.008
fin_parametres_reglage_Hysteresis_
fin_liste_lois_elementaires
#=========================================================================
Le K = 5.25e5 de la partie hyper représente le comportement volumique élastique de la loi
Qor = 3 : de la partie hyper indique un seuil quasiment nul
muzero = 1000. : de la partie hyper indique un écrouissage cinématique
Dans le cas où on veut une plasticité parfaite, il faut prendre une loi bulk uniquement: ISOHYPERBULK3 par exemple.
Pour la partie hystérésis: mu= 80000. indique la pente à l'origine en cisaillement tau/gamma
ce qui correspond au G de la loi de Hooke G = E/(2*(1+nu))
Qzero= 430 indique le seuil de saturation en cisaillement: tau maxi = Qzero/sqrt(2)
np contrôle la courbure du passage entre la pente initiale et la saturation.
Pour avoir une idée de la loi, il est intéressant de faire un petit test de traction simple avec un cube fixé isostatiquement et soumis à un déplacement imposé.