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

bugfix: 

wrong results when calculating multiple strain tensors. Introduced in cd6a4d1cfd
This commit is contained in:
Martin Diehl 2019-12-23 06:48:42 +01:00
parent 3ed9126845
commit 3c4c364107
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
processing/post/addStrainTensors.py
View File

@ -15,12 +15,12 @@ scriptID = ' '.join([scriptName,damask.version])
def operator(stretch,strain,eigenvalues): def operator(stretch,strain,eigenvalues):
"""Albrecht Bertram: Elasticity and Plasticity of Large Deformations An Introduction (3rd Edition, 2012), p. 102.""" """Albrecht Bertram: Elasticity and Plasticity of Large Deformations An Introduction (3rd Edition, 2012), p. 102."""
return { return {
'V#ln': np.log(eigenvalues) , 'V#ln': np.log(eigenvalues) ,
'U#ln': np.log(eigenvalues) , 'U#ln': np.log(eigenvalues) ,
'V#Biot': ( np.ones(3,'d') - 1.0/eigenvalues ) , 'V#Biot': ( np.ones(3,'d') - eigenvalues**-1.0) ,
'U#Biot': ( eigenvalues - np.ones(3,'d') ) , 'U#Biot': ( eigenvalues**1.0 - np.ones(3,'d')) ,
'V#Green': ( np.ones(3,'d') - 1.0/eigenvalues**2.0) *0.5, 'V#Green': ( np.ones(3,'d') - eigenvalues**-2.0)*0.5,
'U#Green': ( eigenvalues**2.0 - np.ones(3,'d')) *0.5, 'U#Green': ( eigenvalues**2.0 - np.ones(3,'d')) *0.5,
}[stretch+'#'+strain] }[stretch+'#'+strain]
@ -142,7 +142,7 @@ for name in filenames:
for theStretch in stretches: for theStretch in stretches:
(D,V) = np.linalg.eigh((stretch[theStretch]+stretch[theStretch].T)*0.5) # eigen decomposition (of symmetric(ed) matrix) (D,V) = np.linalg.eigh((stretch[theStretch]+stretch[theStretch].T)*0.5) # eigen decomposition (of symmetric(ed) matrix)
for theStrain in strains: for theStrain in strains:
d = operator(theStretch,theStrain,D) # operate on eigenvalues of U or V d = operator(theStretch,theStrain,D) # operate on eigenvalues of U or V
eps = np.dot(V,np.dot(np.diag(d),V.T)).reshape(9) # build tensor back from eigenvalue/vector basis eps = np.dot(V,np.dot(np.diag(d),V.T)).reshape(9) # build tensor back from eigenvalue/vector basis