Homogenised C per grain can now depend on microstructure (important for deformation twinning)

This commit is contained in:
Luc Hantcherli 2007-12-07 10:34:21 +00:00
parent 3c0d51542e
commit 975c113ae8
1 changed files with 2 additions and 1 deletions

View File

@ -403,7 +403,7 @@ if(mod(CPFEM_cn,2)==0) then
return return
endif endif
C_66 = constitutive_HomogenizedC(grain, CPFEM_in, cp_en)
A = matmul(Fg_new,invFp_old) ! actually Fe A = matmul(Fg_new,invFp_old) ! actually Fe
A = matmul(transpose(A), A) A = matmul(transpose(A), A)
@ -422,6 +422,7 @@ state: do ! outer iteration: state
return return
endif endif
call constitutive_Microstructure(state_new,CPFEM_Temperature(CPFEM_in,cp_en),grain,CPFEM_in,cp_en) call constitutive_Microstructure(state_new,CPFEM_Temperature(CPFEM_in,cp_en),grain,CPFEM_in,cp_en)
C_66 = constitutive_HomogenizedC(state_new, grain, CPFEM_in, cp_en)
iStress = 0_pInt iStress = 0_pInt
stress: do ! inner iteration: stress stress: do ! inner iteration: stress
iStress = iStress+1 iStress = iStress+1