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

Constitutive_dislo.f90 contains now a modified version of HomogenisedC that accounts for twin reoriented regions. 

ONLY constitutive_dislo was changed, NOT constitutive.pheno
This commit is contained in:
Luc Hantcherli 2007-12-07 14:17:42 +00:00
parent b69aa2b112
commit 44d6235777
1 changed files with 11 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
trunk/constitutive_dislo.f90
View File

@ -750,6 +750,7 @@ do i=1,material_maxN
enddo enddo
enddo enddo
enddo enddo
!* Mapping back to 66-formulation of the matices
material_Ctwin_66(:,:,j,i) = math_Mandel3333to66(material_Ctwin_3333(:,:,:,:,j,i)) material_Ctwin_66(:,:,j,i) = math_Mandel3333to66(material_Ctwin_3333(:,:,:,:,j,i))
enddo enddo
enddo enddo
@ -790,11 +791,20 @@ implicit none
!* Definition of variables !* Definition of variables
integer(pInt) ipc,ip,el integer(pInt) ipc,ip,el
integer(pInt) matID,i
real(pReal), dimension(6,6) :: constitutive_homogenizedC real(pReal), dimension(6,6) :: constitutive_homogenizedC
real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: state real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: state
!* Get the material-ID from the triplet(ipc,ip,el)
matID = constitutive_matID(ipc,ip,el)
!* Homogenization scheme !* Homogenization scheme
constitutive_homogenizedC=material_Cslip_66(:,:,constitutive_matID(ipc,ip,el)) constitutive_homogenizedC=(1-sum(state((material_Nslip(matID)+1):(material_Nslip(matID)+material_Ntwin(matID)))))*&
material_Cslip_66(:,:,constitutive_matID(ipc,ip,el))
do i=1,material_Ntwin(matID)
constitutive_homogenizedC=constitutive_homogenizedC+state((material_Nslip(matID)+i))*&
material_Ctwin_66(:,:,i,constitutive_matID(ipc,ip,el))
enddo
return return
end function end function