From 50e874f8873e412fb82304dab085c9b979deb2e2 Mon Sep 17 00:00:00 2001
From: Pratheek Shanthraj
Date: Tue, 22 Jan 2013 15:48:47 +0000
Subject: [PATCH] fixed bug in elasticity matrix calculation
---
code/constitutive_dislotwin.f90 | 2 +-
code/constitutive_j2.f90 | 2 +-
code/constitutive_none.f90 | 2 +-
code/constitutive_nonlocal.f90 | 2 +-
code/constitutive_phenopowerlaw.f90 | 2 +-
code/constitutive_titanmod.f90 | 2 +-
6 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90
index daefcf26d..48e089f46 100644
--- a/code/constitutive_dislotwin.f90
+++ b/code/constitutive_dislotwin.f90
@@ -543,7 +543,7 @@ do i = 1_pInt,maxNinstance
!* Elasticity matrix and shear modulus according to material.config
constitutive_dislotwin_Cslip_66(1:6,1:6,i) = lattice_symmetrizeC66(constitutive_dislotwin_structureName(i),&
- constitutive_dislotwin_Cslip_66)
+ constitutive_dislotwin_Cslip_66(:,:,i))
constitutive_dislotwin_Gmod(i) = &
0.2_pReal*(constitutive_dislotwin_Cslip_66(1,1,i)-constitutive_dislotwin_Cslip_66(1,2,i)) &
+0.3_pReal*constitutive_dislotwin_Cslip_66(4,4,i)
diff --git a/code/constitutive_j2.f90 b/code/constitutive_j2.f90
index 938a17199..b4de1bbc7 100644
--- a/code/constitutive_j2.f90
+++ b/code/constitutive_j2.f90
@@ -305,7 +305,7 @@ subroutine constitutive_j2_init(myFile)
constitutive_j2_sizeState(i) = 1_pInt
constitutive_j2_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_j2_structureName(i),&
- constitutive_j2_Cslip_66)
+ constitutive_j2_Cslip_66(:,:,i))
constitutive_j2_Cslip_66(1:6,1:6,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_j2_Cslip_66(1:6,1:6,i)))
diff --git a/code/constitutive_none.f90 b/code/constitutive_none.f90
index c1eb16bef..dddb99448 100644
--- a/code/constitutive_none.f90
+++ b/code/constitutive_none.f90
@@ -182,7 +182,7 @@ subroutine constitutive_none_init(myFile)
constitutive_none_sizeState(i) = 1_pInt
constitutive_none_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_none_structureName(i),&
- constitutive_none_Cslip_66)
+ constitutive_none_Cslip_66(:,:,i))
constitutive_none_Cslip_66(:,:,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_none_Cslip_66(:,:,i)))
diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90
index 53146e950..192ee8513 100644
--- a/code/constitutive_nonlocal.f90
+++ b/code/constitutive_nonlocal.f90
@@ -855,7 +855,7 @@ do i = 1,maxNinstance
!*** elasticity matrix and shear modulus according to material.config
constitutive_nonlocal_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_nonlocal_structureName(i),&
- constitutive_nonlocal_Cslip_66)
+ constitutive_nonlocal_Cslip_66(:,:,i))
constitutive_nonlocal_Gmod(i) = 0.2_pReal * ( constitutive_nonlocal_Cslip_66(1,1,i) - constitutive_nonlocal_Cslip_66(1,2,i) &
+ 3.0_pReal*constitutive_nonlocal_Cslip_66(4,4,i) ) ! (C11iso-C12iso)/2 with C11iso=(3*C11+2*C12+4*C44)/5 and C12iso=(C11+4*C12-2*C44)/5
constitutive_nonlocal_nu(i) = ( constitutive_nonlocal_Cslip_66(1,1,i) + 4.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
diff --git a/code/constitutive_phenopowerlaw.f90 b/code/constitutive_phenopowerlaw.f90
index b71331995..7cbbdb816 100644
--- a/code/constitutive_phenopowerlaw.f90
+++ b/code/constitutive_phenopowerlaw.f90
@@ -451,7 +451,7 @@ subroutine constitutive_phenopowerlaw_init(myFile)
myStructure = constitutive_phenopowerlaw_structure(i)
constitutive_phenopowerlaw_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_phenopowerlaw_structureName(i),&
- constitutive_phenopowerlaw_Cslip_66)
+ constitutive_phenopowerlaw_Cslip_66(:,:,i))
! assign elasticity tensor
constitutive_phenopowerlaw_Cslip_66(:,:,i) = &
math_Mandel3333to66(math_Voigt66to3333(constitutive_phenopowerlaw_Cslip_66(:,:,i)))
diff --git a/code/constitutive_titanmod.f90 b/code/constitutive_titanmod.f90
index 25432d067..b67738b59 100644
--- a/code/constitutive_titanmod.f90
+++ b/code/constitutive_titanmod.f90
@@ -868,7 +868,7 @@ do i = 1_pInt,maxNinstance
!* Elasticity matrix and shear modulus according to material.config
constitutive_titanmod_Cslip_66(:,:,i) = lattice_symmetrizeC66(constitutive_titanmod_structureName(i),&
- constitutive_titanmod_Cslip_66)
+ constitutive_titanmod_Cslip_66(:,:,i))
constitutive_titanmod_Gmod(i) = &
0.2_pReal*(constitutive_titanmod_Cslip_66(1,1,i)-constitutive_titanmod_Cslip_66(1,2,i))&
+ 0.3_pReal*constitutive_titanmod_Cslip_66(4,4,i)