223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

simplified

This commit is contained in:
Martin Diehl 2020-11-18 14:01:52 +01:00
parent 89c748fdf1
commit 1c9028d0b1
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
python/damask/mechanics.py
View File

@ -60,8 +60,7 @@ def stress_Cauchy(P,F):
Cauchy stress. Cauchy stress.
""" """
sigma = _np.einsum('...,...ij,...kj',1.0/_np.linalg.det(F),P,F) return tensor.symmetric(_np.einsum('...,...ij,...kj',1.0/_np.linalg.det(F),P,F))
return tensor.symmetric(sigma)
def deviatoric_part(T): def deviatoric_part(T):
@ -79,7 +78,7 @@ def deviatoric_part(T):
Deviatoric part of T. Deviatoric part of T.
""" """
return T - _np.einsum('...ij,...',_np.eye(3),spherical_part(T)) return T - spherical_part(T,tensor=True)
def maximum_shear(T_sym): def maximum_shear(T_sym):
@ -157,8 +156,7 @@ def stress_second_Piola_Kirchhoff(P,F):
Second Piola-Kirchhoff stress. Second Piola-Kirchhoff stress.
""" """
S = _np.einsum('...jk,...kl',_np.linalg.inv(F),P) return tensor.symmetric(_np.einsum('...ij,...jk',_np.linalg.inv(F),P))
return tensor.symmetric(S)
def rotational_part(T): def rotational_part(T):