fixed wrong? homogenizedC calculation
This commit is contained in:
parent
ea66511160
commit
1a0e0f4fa2
|
@ -1032,7 +1032,7 @@ function constitutive_dislokmc_homogenizedC(ipc,ip,el)
|
|||
constitutive_dislokmc_homogenizedC = (1.0_pReal-sumf)*lattice_C66(1:6,1:6,ph)
|
||||
do i=1_pInt,nt
|
||||
constitutive_dislokmc_homogenizedC = &
|
||||
constitutive_dislokmc_homogenizedC + plasticState(ph)%state(3_pInt*ns+i, of)*lattice_C66(1:6,1:6,ph)
|
||||
constitutive_dislokmc_homogenizedC + plasticState(ph)%state(3_pInt*ns+i, of)*lattice_Ctwin66(1:6,1:6,ph)
|
||||
enddo
|
||||
|
||||
end function constitutive_dislokmc_homogenizedC
|
||||
|
|
|
@ -1057,7 +1057,7 @@ function constitutive_dislotwin_homogenizedC(ipc,ip,el)
|
|||
constitutive_dislotwin_homogenizedC = (1.0_pReal-sumf)*lattice_C66(1:6,1:6,ph)
|
||||
do i=1_pInt,nt
|
||||
constitutive_dislotwin_homogenizedC = &
|
||||
constitutive_dislotwin_homogenizedC + plasticState(ph)%state(3_pInt*ns+i, of)*lattice_C66(1:6,1:6,ph)
|
||||
constitutive_dislotwin_homogenizedC + plasticState(ph)%state(3_pInt*ns+i, of)*lattice_Ctwin66(1:6,1:6,ph)
|
||||
enddo
|
||||
|
||||
end function constitutive_dislotwin_homogenizedC
|
||||
|
|
Loading…
Reference in New Issue