From 44d62357774975483889ed3f8e4c69c4256ceb73 Mon Sep 17 00:00:00 2001 From: Luc Hantcherli Date: Fri, 7 Dec 2007 14:17:42 +0000 Subject: [PATCH] Constitutive_dislo.f90 contains now a modified version of HomogenisedC that accounts for twin reoriented regions. ONLY constitutive_dislo was changed, NOT constitutive.pheno --- trunk/constitutive_dislo.f90 | 12 +++++++++++- 1 file changed, 11 insertions(+), 1 deletion(-) diff --git a/trunk/constitutive_dislo.f90 b/trunk/constitutive_dislo.f90 index b781acb0c..5bdbe8d06 100644 --- a/trunk/constitutive_dislo.f90 +++ b/trunk/constitutive_dislo.f90 @@ -750,6 +750,7 @@ do i=1,material_maxN enddo enddo enddo + !* Mapping back to 66-formulation of the matices material_Ctwin_66(:,:,j,i) = math_Mandel3333to66(material_Ctwin_3333(:,:,:,:,j,i)) enddo enddo @@ -790,11 +791,20 @@ implicit none !* Definition of variables integer(pInt) ipc,ip,el +integer(pInt) matID,i real(pReal), dimension(6,6) :: constitutive_homogenizedC 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 -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 end function