crystallite should be responsible of crystallite variables
This commit is contained in:
Martin Diehl
parent
7f3a1a0238
commit
7d7eff0d94
|
|
@ -94,6 +94,7 @@ subroutine CPFEM_initAll(el,ip)
|
||||||
call crystallite_init
|
call crystallite_init
|
||||||
call homogenization_init
|
call homogenization_init
|
||||||
call CPFEM_init
|
call CPFEM_init
|
||||||
|
call CPFEM_initX
|
||||||
CPFEM_init_done = .true.
|
CPFEM_init_done = .true.
|
||||||
endif
|
endif
|
||||||
!$OMP END CRITICAL (init)
|
!$OMP END CRITICAL (init)
|
||||||
|
|
@ -174,35 +175,7 @@ subroutine CPFEM_general(mode, parallelExecution, ffn, ffn1, temperature_inp, dt
|
||||||
CPFEM_dcsde = CPFEM_dcsde_knownGood
|
CPFEM_dcsde = CPFEM_dcsde_knownGood
|
||||||
|
|
||||||
!*** age results
|
!*** age results
|
||||||
if (iand(mode, CPFEM_AGERESULTS) /= 0_pInt) then
|
if (iand(mode, CPFEM_AGERESULTS) /= 0_pInt) call CPFEM_forward
|
||||||
crystallite_F0 = crystallite_partionedF ! crystallite deformation
|
|
||||||
crystallite_Fp0 = crystallite_Fp ! crystallite plastic deformation
|
|
||||||
crystallite_Lp0 = crystallite_Lp ! crystallite plastic velocity
|
|
||||||
crystallite_Fi0 = crystallite_Fi ! crystallite intermediate deformation
|
|
||||||
crystallite_Li0 = crystallite_Li ! crystallite intermediate velocity
|
|
||||||
crystallite_S0 = crystallite_S ! crystallite 2nd Piola Kirchhoff stress
|
|
||||||
|
|
||||||
forall (i = 1:size(plasticState)) plasticState(i)%state0 = plasticState(i)%state
|
|
||||||
do i = 1, size(sourceState)
|
|
||||||
do mySource = 1,phase_Nsources(i)
|
|
||||||
sourceState(i)%p(mySource)%state0 = sourceState(i)%p(mySource)%state
|
|
||||||
enddo; enddo
|
|
||||||
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then
|
|
||||||
write(6,'(a)') '<< CPFEM >> aging states'
|
|
||||||
if (debug_e <= discretization_nElem .and. debug_i <=discretization_nIP) then
|
|
||||||
write(6,'(a,1x,i8,1x,i2,1x,i4,/,(12x,6(e20.8,1x)),/)') &
|
|
||||||
'<< CPFEM >> aged state of elFE ip grain',debug_e, debug_i, 1, &
|
|
||||||
plasticState(material_phaseAt(1,debug_e))%state(:,material_phasememberAt(1,debug_i,debug_e))
|
|
||||||
endif
|
|
||||||
endif
|
|
||||||
|
|
||||||
do homog = 1_pInt, material_Nhomogenization
|
|
||||||
homogState (homog)%state0 = homogState (homog)%state
|
|
||||||
thermalState (homog)%state0 = thermalState (homog)%state
|
|
||||||
damageState (homog)%state0 = damageState (homog)%state
|
|
||||||
enddo
|
|
||||||
endif
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
!*** collection of FEM input with returning of randomize odd stress and jacobian
|
!*** collection of FEM input with returning of randomize odd stress and jacobian
|
||||||
|
|
|
||||||
107
src/CPFEM2.f90
107
src/CPFEM2.f90
|
|
@ -59,6 +59,7 @@ subroutine CPFEM_initAll
|
||||||
call crystallite_init
|
call crystallite_init
|
||||||
call homogenization_init
|
call homogenization_init
|
||||||
call CPFEM_init
|
call CPFEM_init
|
||||||
|
call CPFEM_initX
|
||||||
|
|
||||||
end subroutine CPFEM_initAll
|
end subroutine CPFEM_initAll
|
||||||
|
|
||||||
|
|
@ -68,117 +69,11 @@ end subroutine CPFEM_initAll
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
subroutine CPFEM_init
|
subroutine CPFEM_init
|
||||||
|
|
||||||
integer :: i
|
|
||||||
integer(HID_T) :: fileHandle, groupHandle
|
|
||||||
character(len=pStringLen) :: fileName, datasetName
|
|
||||||
|
|
||||||
write(6,'(/,a)') ' <<<+- CPFEM init -+>>>'; flush(6)
|
write(6,'(/,a)') ' <<<+- CPFEM init -+>>>'; flush(6)
|
||||||
|
|
||||||
if (interface_restartInc > 0) then
|
|
||||||
write(6,'(/,a,i0,a)') ' reading restart information of increment ', interface_restartInc, ' from file'
|
|
||||||
|
|
||||||
write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
|
|
||||||
fileHandle = HDF5_openFile(fileName)
|
|
||||||
|
|
||||||
call HDF5_read(fileHandle,crystallite_F0, 'F')
|
|
||||||
call HDF5_read(fileHandle,crystallite_Fp0,'Fp')
|
|
||||||
call HDF5_read(fileHandle,crystallite_Fi0,'Fi')
|
|
||||||
call HDF5_read(fileHandle,crystallite_Lp0,'Lp')
|
|
||||||
call HDF5_read(fileHandle,crystallite_Li0,'Li')
|
|
||||||
call HDF5_read(fileHandle,crystallite_S0, 'S')
|
|
||||||
|
|
||||||
groupHandle = HDF5_openGroup(fileHandle,'constituent')
|
|
||||||
do i = 1,size(phase_plasticity)
|
|
||||||
write(datasetName,'(i0,a)') i,'_omega_plastic'
|
|
||||||
call HDF5_read(groupHandle,plasticState(i)%state0,datasetName)
|
|
||||||
enddo
|
|
||||||
call HDF5_closeGroup(groupHandle)
|
|
||||||
|
|
||||||
groupHandle = HDF5_openGroup(fileHandle,'materialpoint')
|
|
||||||
do i = 1, material_Nhomogenization
|
|
||||||
write(datasetName,'(i0,a)') i,'_omega_homogenization'
|
|
||||||
call HDF5_read(groupHandle,homogState(i)%state0,datasetName)
|
|
||||||
enddo
|
|
||||||
call HDF5_closeGroup(groupHandle)
|
|
||||||
|
|
||||||
call HDF5_closeFile(fileHandle)
|
|
||||||
endif
|
|
||||||
|
|
||||||
end subroutine CPFEM_init
|
end subroutine CPFEM_init
|
||||||
|
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
|
||||||
!> @brief Forward data after successful increment.
|
|
||||||
! ToDo: Any guessing for the current states possible?
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
|
||||||
subroutine CPFEM_forward
|
|
||||||
|
|
||||||
integer :: i, j
|
|
||||||
|
|
||||||
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0) &
|
|
||||||
write(6,'(a)') '<< CPFEM >> aging states'
|
|
||||||
|
|
||||||
crystallite_F0 = crystallite_partionedF
|
|
||||||
crystallite_Fp0 = crystallite_Fp
|
|
||||||
crystallite_Lp0 = crystallite_Lp
|
|
||||||
crystallite_Fi0 = crystallite_Fi
|
|
||||||
crystallite_Li0 = crystallite_Li
|
|
||||||
crystallite_S0 = crystallite_S
|
|
||||||
|
|
||||||
do i = 1, size(plasticState)
|
|
||||||
plasticState(i)%state0 = plasticState(i)%state
|
|
||||||
enddo
|
|
||||||
do i = 1, size(sourceState)
|
|
||||||
do j = 1,phase_Nsources(i)
|
|
||||||
sourceState(i)%p(j)%state0 = sourceState(i)%p(j)%state
|
|
||||||
enddo; enddo
|
|
||||||
do i = 1, material_Nhomogenization
|
|
||||||
homogState (i)%state0 = homogState (i)%state
|
|
||||||
thermalState(i)%state0 = thermalState(i)%state
|
|
||||||
damageState (i)%state0 = damageState (i)%state
|
|
||||||
enddo
|
|
||||||
|
|
||||||
end subroutine CPFEM_forward
|
|
||||||
|
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
|
||||||
!> @brief Write current restart information (Field and constitutive data) to file.
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
|
||||||
subroutine CPFEM_restartWrite
|
|
||||||
|
|
||||||
integer :: i
|
|
||||||
integer(HID_T) :: fileHandle, groupHandle
|
|
||||||
character(len=pStringLen) :: fileName, datasetName
|
|
||||||
|
|
||||||
write(6,'(a)') ' writing field and constitutive data required for restart to file';flush(6)
|
|
||||||
|
|
||||||
write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
|
|
||||||
fileHandle = HDF5_openFile(fileName,'a')
|
|
||||||
|
|
||||||
call HDF5_write(fileHandle,crystallite_partionedF,'F')
|
|
||||||
call HDF5_write(fileHandle,crystallite_Fp, 'Fp')
|
|
||||||
call HDF5_write(fileHandle,crystallite_Fi, 'Fi')
|
|
||||||
call HDF5_write(fileHandle,crystallite_Lp, 'Lp')
|
|
||||||
call HDF5_write(fileHandle,crystallite_Li, 'Li')
|
|
||||||
call HDF5_write(fileHandle,crystallite_S, 'S')
|
|
||||||
|
|
||||||
groupHandle = HDF5_addGroup(fileHandle,'constituent')
|
|
||||||
do i = 1,size(phase_plasticity)
|
|
||||||
write(datasetName,'(i0,a)') i,'_omega_plastic'
|
|
||||||
call HDF5_write(groupHandle,plasticState(i)%state,datasetName)
|
|
||||||
enddo
|
|
||||||
call HDF5_closeGroup(groupHandle)
|
|
||||||
|
|
||||||
groupHandle = HDF5_addGroup(fileHandle,'materialpoint')
|
|
||||||
do i = 1, material_Nhomogenization
|
|
||||||
write(datasetName,'(i0,a)') i,'_omega_homogenization'
|
|
||||||
call HDF5_write(groupHandle,homogState(i)%state,datasetName)
|
|
||||||
enddo
|
|
||||||
call HDF5_closeGroup(groupHandle)
|
|
||||||
|
|
||||||
call HDF5_closeFile(fileHandle)
|
|
||||||
|
|
||||||
end subroutine CPFEM_restartWrite
|
|
||||||
|
|
||||||
|
|
||||||
!--------------------------------------------------------------------------------------------------
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
|
|
||||||
|
|
@ -43,6 +43,8 @@ module DAMASK_interface
|
||||||
logical, public :: symmetricSolver
|
logical, public :: symmetricSolver
|
||||||
character(len=*), parameter, public :: INPUTFILEEXTENSION = '.dat'
|
character(len=*), parameter, public :: INPUTFILEEXTENSION = '.dat'
|
||||||
|
|
||||||
|
integer(pInt), parameter, public :: interface_restartInc = 0
|
||||||
|
|
||||||
public :: &
|
public :: &
|
||||||
DAMASK_interface_init, &
|
DAMASK_interface_init, &
|
||||||
getSolverJobName
|
getSolverJobName
|
||||||
|
|
|
||||||
|
|
@ -11,6 +11,8 @@
|
||||||
module crystallite
|
module crystallite
|
||||||
use prec
|
use prec
|
||||||
use IO
|
use IO
|
||||||
|
use HDF5_utilities
|
||||||
|
use DAMASK_interface
|
||||||
use config
|
use config
|
||||||
use debug
|
use debug
|
||||||
use numerics
|
use numerics
|
||||||
|
|
@ -104,7 +106,10 @@ module crystallite
|
||||||
crystallite_stressTangent, &
|
crystallite_stressTangent, &
|
||||||
crystallite_orientations, &
|
crystallite_orientations, &
|
||||||
crystallite_push33ToRef, &
|
crystallite_push33ToRef, &
|
||||||
crystallite_results
|
crystallite_results, &
|
||||||
|
CPFEM_restartWrite, &
|
||||||
|
CPFEM_forward, &
|
||||||
|
CPFEM_initX
|
||||||
|
|
||||||
contains
|
contains
|
||||||
|
|
||||||
|
|
@ -1844,4 +1849,118 @@ logical function stateJump(ipc,ip,el)
|
||||||
|
|
||||||
end function stateJump
|
end function stateJump
|
||||||
|
|
||||||
|
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
!> @brief Write current restart information (Field and constitutive data) to file.
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
subroutine CPFEM_restartWrite
|
||||||
|
|
||||||
|
integer :: i
|
||||||
|
integer(HID_T) :: fileHandle, groupHandle
|
||||||
|
character(len=pStringLen) :: fileName, datasetName
|
||||||
|
|
||||||
|
write(6,'(a)') ' writing field and constitutive data required for restart to file';flush(6)
|
||||||
|
|
||||||
|
write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
|
||||||
|
fileHandle = HDF5_openFile(fileName,'a')
|
||||||
|
|
||||||
|
call HDF5_write(fileHandle,crystallite_partionedF,'F')
|
||||||
|
call HDF5_write(fileHandle,crystallite_Fp, 'Fp')
|
||||||
|
call HDF5_write(fileHandle,crystallite_Fi, 'Fi')
|
||||||
|
call HDF5_write(fileHandle,crystallite_Lp, 'Lp')
|
||||||
|
call HDF5_write(fileHandle,crystallite_Li, 'Li')
|
||||||
|
call HDF5_write(fileHandle,crystallite_S, 'S')
|
||||||
|
|
||||||
|
groupHandle = HDF5_addGroup(fileHandle,'constituent')
|
||||||
|
do i = 1,size(phase_plasticity)
|
||||||
|
write(datasetName,'(i0,a)') i,'_omega_plastic'
|
||||||
|
call HDF5_write(groupHandle,plasticState(i)%state,datasetName)
|
||||||
|
enddo
|
||||||
|
call HDF5_closeGroup(groupHandle)
|
||||||
|
|
||||||
|
groupHandle = HDF5_addGroup(fileHandle,'materialpoint')
|
||||||
|
do i = 1, material_Nhomogenization
|
||||||
|
write(datasetName,'(i0,a)') i,'_omega_homogenization'
|
||||||
|
call HDF5_write(groupHandle,homogState(i)%state,datasetName)
|
||||||
|
enddo
|
||||||
|
call HDF5_closeGroup(groupHandle)
|
||||||
|
|
||||||
|
call HDF5_closeFile(fileHandle)
|
||||||
|
|
||||||
|
end subroutine CPFEM_restartWrite
|
||||||
|
|
||||||
|
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
!> @brief Forward data after successful increment.
|
||||||
|
! ToDo: Any guessing for the current states possible?
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
subroutine CPFEM_forward
|
||||||
|
|
||||||
|
integer :: i, j
|
||||||
|
|
||||||
|
crystallite_F0 = crystallite_partionedF
|
||||||
|
crystallite_Fp0 = crystallite_Fp
|
||||||
|
crystallite_Lp0 = crystallite_Lp
|
||||||
|
crystallite_Fi0 = crystallite_Fi
|
||||||
|
crystallite_Li0 = crystallite_Li
|
||||||
|
crystallite_S0 = crystallite_S
|
||||||
|
|
||||||
|
do i = 1, size(plasticState)
|
||||||
|
plasticState(i)%state0 = plasticState(i)%state
|
||||||
|
enddo
|
||||||
|
do i = 1, size(sourceState)
|
||||||
|
do j = 1,phase_Nsources(i)
|
||||||
|
sourceState(i)%p(j)%state0 = sourceState(i)%p(j)%state
|
||||||
|
enddo; enddo
|
||||||
|
do i = 1, material_Nhomogenization
|
||||||
|
homogState (i)%state0 = homogState (i)%state
|
||||||
|
thermalState(i)%state0 = thermalState(i)%state
|
||||||
|
damageState (i)%state0 = damageState (i)%state
|
||||||
|
enddo
|
||||||
|
|
||||||
|
end subroutine CPFEM_forward
|
||||||
|
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
!> @brief allocate the arrays defined in module CPFEM and initialize them
|
||||||
|
!--------------------------------------------------------------------------------------------------
|
||||||
|
subroutine CPFEM_initX
|
||||||
|
|
||||||
|
integer :: i
|
||||||
|
integer(HID_T) :: fileHandle, groupHandle
|
||||||
|
character(len=pStringLen) :: fileName, datasetName
|
||||||
|
|
||||||
|
write(6,'(/,a)') ' <<<+- CPFEM init -+>>>'; flush(6)
|
||||||
|
|
||||||
|
if (interface_restartInc > 0) then
|
||||||
|
write(6,'(/,a,i0,a)') ' reading restart information of increment ', interface_restartInc, ' from file'
|
||||||
|
|
||||||
|
write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
|
||||||
|
fileHandle = HDF5_openFile(fileName)
|
||||||
|
|
||||||
|
call HDF5_read(fileHandle,crystallite_F0, 'F')
|
||||||
|
call HDF5_read(fileHandle,crystallite_Fp0,'Fp')
|
||||||
|
call HDF5_read(fileHandle,crystallite_Fi0,'Fi')
|
||||||
|
call HDF5_read(fileHandle,crystallite_Lp0,'Lp')
|
||||||
|
call HDF5_read(fileHandle,crystallite_Li0,'Li')
|
||||||
|
call HDF5_read(fileHandle,crystallite_S0, 'S')
|
||||||
|
|
||||||
|
groupHandle = HDF5_openGroup(fileHandle,'constituent')
|
||||||
|
do i = 1,size(phase_plasticity)
|
||||||
|
write(datasetName,'(i0,a)') i,'_omega_plastic'
|
||||||
|
call HDF5_read(groupHandle,plasticState(i)%state0,datasetName)
|
||||||
|
enddo
|
||||||
|
call HDF5_closeGroup(groupHandle)
|
||||||
|
|
||||||
|
groupHandle = HDF5_openGroup(fileHandle,'materialpoint')
|
||||||
|
do i = 1, material_Nhomogenization
|
||||||
|
write(datasetName,'(i0,a)') i,'_omega_homogenization'
|
||||||
|
call HDF5_read(groupHandle,homogState(i)%state0,datasetName)
|
||||||
|
enddo
|
||||||
|
call HDF5_closeGroup(groupHandle)
|
||||||
|
|
||||||
|
call HDF5_closeFile(fileHandle)
|
||||||
|
endif
|
||||||
|
|
||||||
|
end subroutine CPFEM_initX
|
||||||
|
|
||||||
end module crystallite
|
end module crystallite
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue