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

hybrid IA sampling now respects sample symmetry

This commit is contained in:
Philip Eisenlohr 2007-10-24 15:13:10 +00:00
parent bf35fa4868
commit 20693b9d38
2 changed files with 24 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
trunk/constitutive_dislo.f90
View File

@ -584,23 +584,22 @@ ODFmap = 0_pInt ! blank mapping
sampleCount = 0_pInt ! count orientations assigned per texture sampleCount = 0_pInt ! count orientations assigned per texture
do texID=1,texture_maxN do texID=1,texture_maxN
if (texture_ODFfile(texID)=='') then select case (texture_symmetry(texID)) ! set symmetry factor
case ('orthotropic')
Nsym(texID) = 4_pInt
case ('monoclinic')
Nsym(texID) = 2_pInt
case default
Nsym(texID) = 1_pInt
end select
if (texture_ODFfile(texID)=='') then ! texture components
sumVolfrac(texID) = sum(texture_gauss(5,:,texID))+sum(texture_fiber(6,:,texID)) sumVolfrac(texID) = sum(texture_gauss(5,:,texID))+sum(texture_fiber(6,:,texID))
if (sumVolfrac(texID)<1.0_pReal) texture_NRandom(texID) = 1_pInt ! check whether random component missing if (sumVolfrac(texID)<1.0_pReal) texture_NRandom(texID) = 1_pInt ! check whether random component missing
select case (texture_symmetry(texID)) ! set symmetry factor
case ('orthotropic')
Nsym(texID) = 4_pInt
case ('monoclinic')
Nsym(texID) = 2_pInt
case default
Nsym(texID) = 1_pInt
end select
Ncomponents(texID) = texture_NGauss(texID)+texture_NFiber(texID)+texture_NRandom(texID) Ncomponents(texID) = texture_NGauss(texID)+texture_NFiber(texID)+texture_NRandom(texID)
else ! hybrid IA else ! hybrid IA
o = o+1 o = o+1
ODFmap(texID) = o ! remember mapping ODFmap(texID) = o ! remember mapping
Ncomponents(texID) = 1_pInt ! single "component" Ncomponents(texID) = 1_pInt ! single "component"
Nsym(texID) = 1_pInt ! no symmetry (use full ODF instead)
endif endif
! adjust multiplicity and number of grains per IP of components ! adjust multiplicity and number of grains per IP of components
multiplicity(texID) = max(1_pInt,texture_Ngrains(texID)/Ncomponents(texID)/Nsym(texID)) multiplicity(texID) = max(1_pInt,texture_Ngrains(texID)/Ncomponents(texID)/Nsym(texID))
@ -624,10 +623,10 @@ do i=1,mesh_NcpElems
enddo enddo
! generate hybridIA samplings for ODFfile textures to later draw from these populations ! generate hybridIA samplings for ODFfile textures to later draw from these populations
allocate(hybridIA_population(3,maxval(texture_totalNgrains,ODFmap /= 0),o)) allocate(hybridIA_population(3,maxval(texture_totalNgrains/Nsym,ODFmap /= 0),o))
do texID = 1,texture_maxN do texID = 1,texture_maxN
if (ODFmap(texID) > 0) & if (ODFmap(texID) > 0) &
hybridIA_population(:,:,ODFmap(texID)) = IO_hybridIA(texture_totalNgrains(texID),texture_ODFfile(texID)) hybridIA_population(:,:,ODFmap(texID)) = IO_hybridIA(texture_totalNgrains(texID)/Nsym(texID),texture_ODFfile(texID))
enddo enddo
!* Array allocation !* Array allocation
@ -695,7 +694,6 @@ do e=1,mesh_NcpElems
endif endif
do s = 1,Nsym(texID)*o ! loop over orientations to be assigned to ip (ex multiplicity) do s = 1,Nsym(texID)*o ! loop over orientations to be assigned to ip (ex multiplicity)
g = g+1 ! next "grain" g = g+1 ! next "grain"
sampleCount(texID) = sampleCount(texID)+1 ! next member of population
constitutive_matID(g,i,e) = matID ! copy matID of element constitutive_matID(g,i,e) = matID ! copy matID of element
constitutive_texID(g,i,e) = texID ! copy texID of element constitutive_texID(g,i,e) = texID ! copy texID of element
constitutive_MatVolFrac(g,i,e) = 1.0_pReal ! singular material (so far) constitutive_MatVolFrac(g,i,e) = 1.0_pReal ! singular material (so far)
@ -708,6 +706,7 @@ do e=1,mesh_NcpElems
constitutive_state_new(l,g,i,e) = material_rho0(matID) constitutive_state_new(l,g,i,e) = material_rho0(matID)
end forall end forall
enddo ! components enddo ! components
sampleCount(texID) = sampleCount(texID)+1 ! next member of hybrid IA population
enddo ! multiplicity enddo ! multiplicity
enddo ! ip enddo ! ip
enddo ! cp_element enddo ! cp_element

25
trunk/constitutive_pheno.f90
View File

@ -562,23 +562,22 @@ ODFmap = 0_pInt ! blank mapping
sampleCount = 0_pInt ! count orientations assigned per texture sampleCount = 0_pInt ! count orientations assigned per texture
do texID=1,texture_maxN do texID=1,texture_maxN
if (texture_ODFfile(texID)=='') then select case (texture_symmetry(texID)) ! set symmetry factor
case ('orthotropic')
Nsym(texID) = 4_pInt
case ('monoclinic')
Nsym(texID) = 2_pInt
case default
Nsym(texID) = 1_pInt
end select
if (texture_ODFfile(texID)=='') then ! texture components
sumVolfrac(texID) = sum(texture_gauss(5,:,texID))+sum(texture_fiber(6,:,texID)) sumVolfrac(texID) = sum(texture_gauss(5,:,texID))+sum(texture_fiber(6,:,texID))
if (sumVolfrac(texID)<1.0_pReal) texture_NRandom(texID) = 1_pInt ! check whether random component missing if (sumVolfrac(texID)<1.0_pReal) texture_NRandom(texID) = 1_pInt ! check whether random component missing
select case (texture_symmetry(texID)) ! set symmetry factor
case ('orthotropic')
Nsym(texID) = 4_pInt
case ('monoclinic')
Nsym(texID) = 2_pInt
case default
Nsym(texID) = 1_pInt
end select
Ncomponents(texID) = texture_NGauss(texID)+texture_NFiber(texID)+texture_NRandom(texID) Ncomponents(texID) = texture_NGauss(texID)+texture_NFiber(texID)+texture_NRandom(texID)
else ! hybrid IA else ! hybrid IA
o = o+1 o = o+1
ODFmap(texID) = o ! remember mapping ODFmap(texID) = o ! remember mapping
Ncomponents(texID) = 1_pInt ! single "component" Ncomponents(texID) = 1_pInt ! single "component"
Nsym(texID) = 1_pInt ! no symmetry (use full ODF instead)
endif endif
! adjust multiplicity and number of grains per IP of components ! adjust multiplicity and number of grains per IP of components
multiplicity(texID) = max(1_pInt,texture_Ngrains(texID)/Ncomponents(texID)/Nsym(texID)) multiplicity(texID) = max(1_pInt,texture_Ngrains(texID)/Ncomponents(texID)/Nsym(texID))
@ -603,10 +602,10 @@ do i=1,mesh_NcpElems
enddo enddo
! generate hybridIA samplings for ODFfile textures to later draw from these populations ! generate hybridIA samplings for ODFfile textures to later draw from these populations
allocate(hybridIA_population(3,maxval(texture_totalNgrains,ODFmap /= 0),o)) allocate(hybridIA_population(3,maxval(texture_totalNgrains/Nsym,ODFmap /= 0),o))
do texID = 1,texture_maxN do texID = 1,texture_maxN
if (ODFmap(texID) > 0) & if (ODFmap(texID) > 0) &
hybridIA_population(:,:,ODFmap(texID)) = IO_hybridIA(texture_totalNgrains(texID),texture_ODFfile(texID)) hybridIA_population(:,:,ODFmap(texID)) = IO_hybridIA(texture_totalNgrains(texID)/Nsym(texID),texture_ODFfile(texID))
enddo enddo
!* Array allocation !* Array allocation
@ -665,7 +664,6 @@ do e=1,mesh_NcpElems
endif endif
do s = 1,Nsym(texID)*o ! loop over orientations to be assigned to ip (ex multiplicity) do s = 1,Nsym(texID)*o ! loop over orientations to be assigned to ip (ex multiplicity)
g = g+1 ! next "grain" g = g+1 ! next "grain"
sampleCount(texID) = sampleCount(texID)+1 ! next member of population
constitutive_matID(g,i,e) = matID ! copy matID of element constitutive_matID(g,i,e) = matID ! copy matID of element
constitutive_texID(g,i,e) = texID ! copy texID of element constitutive_texID(g,i,e) = texID ! copy texID of element
constitutive_MatVolFrac(g,i,e) = 1.0_pReal ! singular material (so far) constitutive_MatVolFrac(g,i,e) = 1.0_pReal ! singular material (so far)
@ -678,6 +676,7 @@ do e=1,mesh_NcpElems
constitutive_state_new(l,g,i,e) = material_s0_slip(matID) constitutive_state_new(l,g,i,e) = material_s0_slip(matID)
end forall end forall
enddo ! components enddo ! components
sampleCount(texID) = sampleCount(texID)+1 ! next member of hybrid IA population
enddo ! multiplicity enddo ! multiplicity
enddo ! ip enddo ! ip
enddo ! cp_element enddo ! cp_element