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

random distribution of initial dislocation density now really independent of meshsize

This commit is contained in:
Christoph Kords 2012-10-02 15:26:58 +00:00
parent d5ce49c471
commit bce8012c89
1 changed files with 4 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
code/constitutive_nonlocal.f90
View File

@ -957,8 +957,7 @@ real(pReal), dimension(2) :: noise
real(pReal), dimension(4) :: rnd real(pReal), dimension(4) :: rnd
real(pReal) meanDensity, & real(pReal) meanDensity, &
totalVolume, & totalVolume, &
linelength, & densityBinning, &
totalLinelength, &
minimumIpVolume minimumIpVolume
@ -983,12 +982,11 @@ do myInstance = 1_pInt,maxNinstance
endif endif
enddo enddo
enddo enddo
linelength = constitutive_nonlocal_rhoSglRandomBinning(myInstance) * minimumIpVolume ** (1.0_pReal / 3.0_pReal) densityBinning = constitutive_nonlocal_rhoSglRandomBinning(myInstance) / minimumIpVolume ** (2.0_pReal / 3.0_pReal)
! subsequently fill random ips with dislocation segments until we reach the desired overall density ! subsequently fill random ips with dislocation segments until we reach the desired overall density
meanDensity = 0.0_pReal meanDensity = 0.0_pReal
totalLinelength = 0.0_pReal
do while(meanDensity < constitutive_nonlocal_rhoSglRandom(myInstance)) do while(meanDensity < constitutive_nonlocal_rhoSglRandom(myInstance))
call random_number(rnd) call random_number(rnd)
el = nint(rnd(1)*real(mesh_NcpElems,pReal)+0.5_pReal,pInt) el = nint(rnd(1)*real(mesh_NcpElems,pReal)+0.5_pReal,pInt)
@ -997,9 +995,8 @@ do myInstance = 1_pInt,maxNinstance
.and. myInstance == phase_plasticityInstance(material_phase(1,ip,el))) then .and. myInstance == phase_plasticityInstance(material_phase(1,ip,el))) then
s = nint(rnd(3)*real(ns,pReal)+0.5_pReal,pInt) s = nint(rnd(3)*real(ns,pReal)+0.5_pReal,pInt)
t = nint(rnd(4)*4.0_pReal+0.5_pReal,pInt) t = nint(rnd(4)*4.0_pReal+0.5_pReal,pInt)
totalLinelength = totalLinelength + linelength meanDensity = meanDensity + densityBinning * mesh_ipVolume(ip,el) / totalVolume
meanDensity = totalLinelength / totalVolume state(1,ip,el)%p((t-1)*ns+s) = state(1,ip,el)%p((t-1)*ns+s) + densityBinning
state(1,ip,el)%p((t-1)*ns+s) = state(1,ip,el)%p((t-1)*ns+s) + linelength / mesh_ipVolume(ip,el)
endif endif
enddo enddo