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

first commit mechanics

This commit is contained in:
Samad Vakili 2020-05-24 18:30:44 +02:00
parent 45e608e71b
commit 987c4a9e8d
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
python/damask/mechanics.py
View File

@ -17,7 +17,8 @@ def Cauchy(P,F):
if _np.shape(F) == _np.shape(P) == (3,3):
sigma = 1.0/_np.linalg.det(F) * _np.dot(P,F.T)
else:
sigma = _np.einsum('i,ijk,ilk->ijl',1.0/_np.linalg.det(F),P,F)
#sigma = _np.einsum('i,ijk,ilk->ijl',1.0/_np.linalg.det(F),P,F)
sigma = _np.einsum('...,...ij,...kj->...ij',1.0/_np.linalg.det(F),P,F)
return symmetric(sigma)
@ -190,7 +191,7 @@ def spherical_part(T,tensor=False):
sph = _np.trace(T)/3.0
return sph if not tensor else _np.eye(3)*sph
else:
sph = _np.trace(T,axis1=1,axis2=2)/3.0
sph = _np.trace(T,axis2=-2,axis1=-1)/3.0
if not tensor:
return sph
else: