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

reworked restarting for compatibility with abaqus (not yet fully working) 

added new orientation feature for direct simulation:
  component type (random) asigns random orientation to an entire grain
This commit is contained in:
Franz Roters 2011-05-28 09:44:43 +00:00
parent 2418dfe96d
commit 08d39342e4
7 changed files with 78 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
code/CPFEM.f90
View File

@ -116,7 +116,7 @@ subroutine CPFEM_init()
use FEsolving, only: parallelExecution, & use FEsolving, only: parallelExecution, &
symmetricSolver, & symmetricSolver, &
restartRead, & restartRead, &
restartJob FEmodelGeometry
use mesh, only: mesh_NcpElems, & use mesh, only: mesh_NcpElems, &
mesh_maxNips mesh_maxNips
use material, only: homogenization_maxNgrains, & use material, only: homogenization_maxNgrains, &
@ -149,31 +149,31 @@ subroutine CPFEM_init()
write(6,'(a)') '<< CPFEM >> Restored state variables of last converged step from binary files' write(6,'(a)') '<< CPFEM >> Restored state variables of last converged step from binary files'
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
endif endif
if (IO_read_jobBinaryFile(777,'recordedPhase',restartJob,size(material_phase))) then if (IO_read_jobBinaryFile(777,'recordedPhase',FEmodelGeometry,size(material_phase))) then
read (777,rec=1) material_phase read (777,rec=1) material_phase
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedF',restartJob,size(crystallite_F0))) then if (IO_read_jobBinaryFile(777,'convergedF',FEmodelGeometry,size(crystallite_F0))) then
read (777,rec=1) crystallite_F0 read (777,rec=1) crystallite_F0
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedFp',restartJob,size(crystallite_Fp0))) then if (IO_read_jobBinaryFile(777,'convergedFp',FEmodelGeometry,size(crystallite_Fp0))) then
read (777,rec=1) crystallite_Fp0 read (777,rec=1) crystallite_Fp0
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedLp',restartJob,size(crystallite_Lp0))) then if (IO_read_jobBinaryFile(777,'convergedLp',FEmodelGeometry,size(crystallite_Lp0))) then
read (777,rec=1) crystallite_Lp0 read (777,rec=1) crystallite_Lp0
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergeddPdF',restartJob,size(crystallite_dPdF0))) then if (IO_read_jobBinaryFile(777,'convergeddPdF',FEmodelGeometry,size(crystallite_dPdF0))) then
read (777,rec=1) crystallite_dPdF0 read (777,rec=1) crystallite_dPdF0
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedTstar',restartJob,size(crystallite_Tstar0_v))) then if (IO_read_jobBinaryFile(777,'convergedTstar',FEmodelGeometry,size(crystallite_Tstar0_v))) then
read (777,rec=1) crystallite_Tstar0_v read (777,rec=1) crystallite_Tstar0_v
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedStateConst',restartJob)) then if (IO_read_jobBinaryFile(777,'convergedStateConst',FEmodelGeometry)) then
m = 0_pInt m = 0_pInt
do i = 1,homogenization_maxNgrains; do j = 1,mesh_maxNips; do k = 1,mesh_NcpElems do i = 1,homogenization_maxNgrains; do j = 1,mesh_maxNips; do k = 1,mesh_NcpElems
do l = 1,size(constitutive_state0(i,j,k)%p) do l = 1,size(constitutive_state0(i,j,k)%p)
@ -183,7 +183,7 @@ subroutine CPFEM_init()
enddo; enddo; enddo enddo; enddo; enddo
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergedStateHomog',restartJob)) then if (IO_read_jobBinaryFile(777,'convergedStateHomog',FEmodelGeometry)) then
m = 0_pInt m = 0_pInt
do k = 1,mesh_NcpElems; do j = 1,mesh_maxNips do k = 1,mesh_NcpElems; do j = 1,mesh_maxNips
do l = 1,homogenization_sizeState(j,k) do l = 1,homogenization_sizeState(j,k)
@ -193,7 +193,7 @@ subroutine CPFEM_init()
enddo; enddo enddo; enddo
close (777) close (777)
endif endif
if (IO_read_jobBinaryFile(777,'convergeddcsdE',restartJob,size(CPFEM_dcsdE))) then if (IO_read_jobBinaryFile(777,'convergeddcsdE',FEmodelGeometry,size(CPFEM_dcsdE))) then
read (777,rec=1) CPFEM_dcsdE read (777,rec=1) CPFEM_dcsdE
close (777) close (777)
endif endif

48
code/FEsolving.f90
View File

@ -36,7 +36,7 @@
logical, dimension(:,:), allocatable :: calcMode logical, dimension(:,:), allocatable :: calcMode
integer(pInt), dimension(:,:), allocatable :: FEsolving_execIP integer(pInt), dimension(:,:), allocatable :: FEsolving_execIP
integer(pInt), dimension(2) :: FEsolving_execElem integer(pInt), dimension(2) :: FEsolving_execElem
character(len=1024) restartJob character(len=1024) FEmodelGeometry
CONTAINS CONTAINS
@ -48,17 +48,20 @@
use prec, only: pInt use prec, only: pInt
use debug, only: debug_verbosity use debug, only: debug_verbosity
use DAMASK_interface, only: getModelName, FEsolver
use IO use IO
implicit none implicit none
integer(pInt), parameter :: fileunit = 222 integer(pInt), parameter :: fileunit = 222
integer(pInt), parameter :: maxNchunks = 6 integer(pInt), parameter :: maxNchunks = 6
integer(pInt) i
integer(pInt), dimension(1+2*maxNchunks) :: positions integer(pInt), dimension(1+2*maxNchunks) :: positions
character(len=64) tag character(len=64) tag
character(len=1024) line character(len=1024) line
FEmodelGeometry = getModelName()
if (IO_open_inputFile(fileunit)) then if (IO_open_inputFile(fileunit,FEmodelGeometry)) then
rewind(fileunit) rewind(fileunit)
do do
@ -75,6 +78,11 @@
positions = IO_stringPos(line,maxNchunks) positions = IO_stringPos(line,maxNchunks)
restartWrite = iand(IO_intValue(line,positions,1),1_pInt) > 0_pInt restartWrite = iand(IO_intValue(line,positions,1),1_pInt) > 0_pInt
restartRead = iand(IO_intValue(line,positions,1),2_pInt) > 0_pInt restartRead = iand(IO_intValue(line,positions,1),2_pInt) > 0_pInt
case ('*restart')
do i=2,positions(1)
restartWrite = IO_lc(IO_StringValue(line,positions,i)) == 'write'
! restartRead = IO_lc(IO_StringValue(line,positions,i)) == 'read'
enddo
end select end select
enddo enddo
else else
@ -83,17 +91,29 @@
100 close(fileunit) 100 close(fileunit)
if (restartRead .and. IO_open_logFile(fileunit)) then if (restartRead) then
rewind(fileunit) if(FEsolver == 'Marc' .and. IO_open_logFile(fileunit)) then
do rewind(fileunit)
read (fileunit,'(a1024)',END=200) line do
positions = IO_stringPos(line,maxNchunks) read (fileunit,'(a1024)',END=200) line
if ( IO_lc(IO_stringValue(line,positions,1)) == 'restart' .and. & positions = IO_stringPos(line,maxNchunks)
IO_lc(IO_stringValue(line,positions,2)) == 'file' .and. & if ( IO_lc(IO_stringValue(line,positions,1)) == 'restart' .and. &
IO_lc(IO_stringValue(line,positions,3)) == 'job' .and. & IO_lc(IO_stringValue(line,positions,2)) == 'file' .and. &
IO_lc(IO_stringValue(line,positions,4)) == 'id' ) & IO_lc(IO_stringValue(line,positions,3)) == 'job' .and. &
restartJob = IO_StringValue(line,positions,6) IO_lc(IO_stringValue(line,positions,4)) == 'id' ) &
enddo FEmodelGeometry = IO_StringValue(line,positions,6)
enddo
elseif (FEsolver == 'Abaqus' .and. IO_open_jobFile(fileunit, 'com')) then
rewind(fileunit)
do
read (fileunit,'(a1024)',END=200) line
positions = IO_stringPos(line,maxNchunks)
! if ( IO_lc(IO_stringValue(line,positions,?)) == 'oldjob?') &
! FEmodelGeometry = IO_StringValue(line,positions,?)
enddo
else
call IO_error(106) ! cannot open file for old job info
endif
endif endif
200 close(fileunit) 200 close(fileunit)
@ -106,7 +126,7 @@
if (debug_verbosity > 0) then if (debug_verbosity > 0) then
write(6,*) 'restart writing: ', restartWrite write(6,*) 'restart writing: ', restartWrite
write(6,*) 'restart reading: ', restartRead write(6,*) 'restart reading: ', restartRead
if (restartRead) write(6,*) 'restart Job: ', trim(restartJob) if (restartRead) write(6,*) 'restart Job: ', trim(FEmodelGeometry)
write(6,*) write(6,*)
endif endif
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)

14
code/IO.f90
View File

@ -26,7 +26,8 @@
!--------------------------- !---------------------------
! function IO_abaqus_assembleInputFile ! function IO_abaqus_assembleInputFile
! function IO_open_file(unit,relPath) ! function IO_open_file(unit,relPath)
! function IO_open_inputFile(unit) ! function IO_open_inputFile(unit, model)
! function IO_open_logFile(unit)
! function IO_hybridIA(Nast,ODFfileName) ! function IO_hybridIA(Nast,ODFfileName)
! private function hybridIA_reps(dV_V,steps,C) ! private function hybridIA_reps(dV_V,steps,C)
! function IO_stringPos(line,maxN) ! function IO_stringPos(line,maxN)
@ -175,27 +176,28 @@ end function
! : changed the function to open *.inp_assembly, which is basically ! : changed the function to open *.inp_assembly, which is basically
! the input file without comment lines and possibly assembled includes ! the input file without comment lines and possibly assembled includes
!******************************************************************** !********************************************************************
logical function IO_open_inputFile(unit) logical function IO_open_inputFile(unit,model)
use prec, only: pReal, pInt use prec, only: pReal, pInt
use DAMASK_interface use DAMASK_interface
implicit none implicit none
integer(pInt), intent(in) :: unit integer(pInt), intent(in) :: unit
character(len=*) model
IO_open_inputFile = .false. IO_open_inputFile = .false.
if (FEsolver == 'Abaqus') then if (FEsolver == 'Abaqus') then
open(unit+1,status='old',err=100,& open(unit+1,status='old',err=100,&
file=trim(getSolverWorkingDirectoryName())//& file=trim(getSolverWorkingDirectoryName())//&
trim(getModelName())//InputFileExtension) trim(model)//InputFileExtension)
open(unit,err=100,file=trim(getSolverWorkingDirectoryName())//& open(unit,err=100,file=trim(getSolverWorkingDirectoryName())//&
trim(getModelName())//InputFileExtension//'_assembly') trim(model)//InputFileExtension//'_assembly')
IO_open_inputFile = IO_abaqus_assembleInputFile(unit,unit+1) ! strip comments and concatenate any "include"s IO_open_inputFile = IO_abaqus_assembleInputFile(unit,unit+1) ! strip comments and concatenate any "include"s
close(unit+1) close(unit+1)
else else
open(unit,status='old',err=100,file=trim(getSolverWorkingDirectoryName())//& open(unit,status='old',err=100,file=trim(getSolverWorkingDirectoryName())//&
trim(getModelName())//InputFileExtension) trim(model)//InputFileExtension)
IO_open_inputFile = .true. IO_open_inputFile = .true.
endif endif
@ -1183,6 +1185,8 @@ endfunction
msg = 'FFTW init error' msg = 'FFTW init error'
case (105) case (105)
msg = 'Error reading from ODF file' msg = 'Error reading from ODF file'
case (106)
msg = 'Error reading info on old job'
case (110) case (110)
msg = 'No homogenization specified via State Variable 2' msg = 'No homogenization specified via State Variable 2'
case (120) case (120)

3
code/config/material.config
View File

@ -393,3 +393,6 @@ symmetry orthotropic # or monoclinic
[fiber example] [fiber example]
# fiber axis in spherical coordinates: alpha crystal system, beta sample system # fiber axis in spherical coordinates: alpha crystal system, beta sample system
(fiber) alpha1 123 alpha2 123 beta1 12 beta2 45 scatter 15 fraction 0.333 (fiber) alpha1 123 alpha2 123 beta1 12 beta2 45 scatter 15 fraction 0.333
[random single crystals]
(random) scatter 0.000 fraction 1.000

21
code/material.f90
View File

@ -445,7 +445,7 @@ subroutine material_parseTexture(file,myPart)
use prec, only: pInt, pReal use prec, only: pInt, pReal
use IO use IO
use math, only: inRad use math, only: inRad, math_sampleRandomOri
implicit none implicit none
character(len=*), intent(in) :: myPart character(len=*), intent(in) :: myPart
@ -467,8 +467,9 @@ subroutine material_parseTexture(file,myPart)
allocate(texture_Ngauss(Nsections)); texture_Ngauss = 0_pInt allocate(texture_Ngauss(Nsections)); texture_Ngauss = 0_pInt
allocate(texture_Nfiber(Nsections)); texture_Nfiber = 0_pInt allocate(texture_Nfiber(Nsections)); texture_Nfiber = 0_pInt
texture_Ngauss = IO_countTagInPart(file,myPart,'(gauss)',Nsections) texture_Ngauss = IO_countTagInPart(file,myPart,'(gauss)', Nsections) + &
texture_Nfiber = IO_countTagInPart(file,myPart,'(fiber)',Nsections) IO_countTagInPart(file,myPart,'(random)',Nsections)
texture_Nfiber = IO_countTagInPart(file,myPart,'(fiber)', Nsections)
texture_maxNgauss = maxval(texture_Ngauss) texture_maxNgauss = maxval(texture_Ngauss)
texture_maxNfiber = maxval(texture_Nfiber) texture_maxNfiber = maxval(texture_Nfiber)
allocate(texture_Gauss (5,texture_maxNgauss,Nsections)); texture_Gauss = 0.0_pReal allocate(texture_Gauss (5,texture_maxNgauss,Nsections)); texture_Gauss = 0.0_pReal
@ -511,6 +512,19 @@ subroutine material_parseTexture(file,myPart)
texture_symmetry(section) = 1 texture_symmetry(section) = 1
end select end select
case ('(random)')
gauss = gauss + 1
texture_Gauss(1:3,gauss,section) = math_sampleRandomOri()
do i = 2,4,2
tag = IO_lc(IO_stringValue(line,positions,i))
select case (tag)
case('scatter')
texture_Gauss(4,gauss,section) = IO_floatValue(line,positions,i+1)*inRad
case('fraction')
texture_Gauss(5,gauss,section) = IO_floatValue(line,positions,i+1)
end select
enddo
case ('(gauss)') case ('(gauss)')
gauss = gauss + 1 gauss = gauss + 1
do i = 2,10,2 do i = 2,10,2
@ -528,6 +542,7 @@ subroutine material_parseTexture(file,myPart)
texture_Gauss(5,gauss,section) = IO_floatValue(line,positions,i+1) texture_Gauss(5,gauss,section) = IO_floatValue(line,positions,i+1)
end select end select
enddo enddo
case ('(fiber)') case ('(fiber)')
fiber = fiber + 1 fiber = fiber + 1
do i = 2,12,2 do i = 2,12,2

1
code/math.f90
View File

@ -166,6 +166,7 @@ real(pReal), dimension(4,36), parameter :: math_symOperations = &
call random_seed(get=randInit) call random_seed(get=randInit)
if (debug_verbosity > 0) then if (debug_verbosity > 0) then
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
! this critical block did cause trouble at IWM
write(6,*) 'random seed: ',randInit(1) write(6,*) 'random seed: ',randInit(1)
write(6,*) write(6,*)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)

6
code/mesh.f90
View File

@ -265,7 +265,7 @@
use DAMASK_interface use DAMASK_interface
use prec, only: pInt use prec, only: pInt
use IO, only: IO_error,IO_open_InputFile,IO_abaqus_hasNoPart use IO, only: IO_error,IO_open_InputFile,IO_abaqus_hasNoPart
use FEsolving, only: parallelExecution, FEsolving_execElem, FEsolving_execIP, calcMode, lastMode use FEsolving, only: parallelExecution, FEsolving_execElem, FEsolving_execIP, calcMode, lastMode, FEmodelGeometry
implicit none implicit none
@ -279,9 +279,9 @@
write(6,*) write(6,*)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
call mesh_build_FEdata() ! --- get properties of the different types of elements call mesh_build_FEdata() ! --- get properties of the different types of elements
if (IO_open_inputFile(fileUnit)) then ! --- parse info from input file... if (IO_open_inputFile(fileUnit,FEmodelGeometry)) then ! --- parse info from input file...
select case (FEsolver) select case (FEsolver)
case ('Spectral') case ('Spectral')