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

653 lines
26 KiB
Fortran
Raw Normal View History

some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
! $Id$
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief CPFEM engine
!--------------------------------------------------------------------------------------------------
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
module CPFEM
some doxygen corrections
2013-03-06 20:11:15 +05:30
use prec, only: &
pReal, &
pInt
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
implicit none
introduced dummy heat calculation, available as an crystallite output temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-19 00:27:28 +05:30
private
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), parameter, private :: &
CPFEM_odd_stress = 1e15_pReal, & !< return value for stress in case of ping pong dummy cycle
CPFEM_odd_jacobian = 1e50_pReal !< return value for jacobian in case of ping pong dummy cycle
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), dimension (:,:,:), allocatable, private :: &
CPFEM_cs !< Cauchy stress
real(pReal), dimension (:,:,:,:), allocatable, private :: &
CPFEM_dcsdE !< Cauchy stress tangent
real(pReal), dimension (:,:,:,:), allocatable, private :: &
CPFEM_dcsdE_knownGood !< known good tangent
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
logical, public, protected :: &
CPFEM_init_done = .false., & !< remember whether init has been done already
CPFEM_calc_done = .false. !< remember whether first ip has already calced the results
integer(pInt), parameter, public :: &
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
CPFEM_CALCRESULTS = 2_pInt**0_pInt, &
CPFEM_AGERESULTS = 2_pInt**1_pInt, &
CPFEM_BACKUPJACOBIAN = 2_pInt**2_pInt, &
CPFEM_RESTOREJACOBIAN = 2_pInt**3_pInt, &
parallelExecution flag now passed as input variable by CPFEM_general; flag is set on the solver level (DAMASK_abaqus_exp.f, DAMASK_marc.f90, etc.)
2013-08-02 16:50:11 +05:30
CPFEM_COLLECT = 2_pInt**4_pInt
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
public :: &
CPFEM_general, &
added "public" (and where possible "protected") statements to functions/subroutines/variables/parameters where needed added -fmodule-private to compile options for gfortran, i.e. gfortran won't compile anymore if an entity is not explicitly declared as public but used in another module to ensure stronger encapsulation
2013-03-28 13:10:30 +05:30
CPFEM_initAll
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
contains
forgot to add those two (rather crucial) files .. ;-)
2009-03-04 19:31:36 +05:30
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief call (thread safe) all module initializations
!--------------------------------------------------------------------------------------------------
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
subroutine CPFEM_initAll(temperature,el,ip)
some doxygen corrections
2013-03-06 20:11:15 +05:30
use prec, only: &
prec_init
use numerics, only: &
numerics_init
use debug, only: &
debug_init
use FEsolving, only: &
FE_init
use math, only: &
math_init
use mesh, only: &
mesh_init
use lattice, only: &
lattice_init
use material, only: &
material_init
use constitutive, only: &
constitutive_init
use crystallite, only: &
crystallite_init
use homogenization, only: &
homogenization_init
use IO, only: &
IO_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use DAMASK_interface
implicit none
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
integer(pInt), intent(in) :: el, & ! FE el number
ip ! FE integration point number
real(pReal), intent(in) :: temperature ! temperature
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
!$OMP CRITICAL (init)
if (.not. CPFEM_init_done) then
reordered initialization for spectral method, corrected bug of deformation BC parsing when prescribing velocity gradient resulting in wrong average deformation
2012-07-30 21:21:48 +05:30
#ifdef Spectral
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call DAMASK_interface_init() ! Spectral solver is interfacing to commandline
reordered initialization for spectral method, corrected bug of deformation BC parsing when prescribing velocity gradient resulting in wrong average deformation
2012-07-30 21:21:48 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call prec_init
call IO_init
call numerics_init
call debug_init
call math_init
call FE_init
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
call mesh_init(ip, el) ! pass on coordinates to alter calcMode of first ip
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call lattice_init
call material_init
call constitutive_init
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
call crystallite_init(temperature) ! (have to) use temperature of first ip for whole model
call homogenization_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call CPFEM_init
introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only. Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver. Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile. Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 21:40:21 +05:30
#ifndef Spectral
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call DAMASK_interface_init() ! Spectral solver init is already done
introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only. Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver. Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile. Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 21:40:21 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
CPFEM_init_done = .true.
endif
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
!$OMP END CRITICAL (init)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine CPFEM_initAll
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief allocate the arrays defined in module CPFEM and initialize them
!--------------------------------------------------------------------------------------------------
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
subroutine CPFEM_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
use prec, only: &
pInt
use IO, only: &
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
IO_read_realFile,&
IO_read_intFile, &
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
IO_timeStamp, &
IO_error
use numerics, only: &
parallelExecution flag now passed as input variable by CPFEM_general; flag is set on the solver level (DAMASK_abaqus_exp.f, DAMASK_marc.f90, etc.)
2013-08-02 16:50:11 +05:30
DAMASK_NumThreadsInt
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use debug, only: &
debug_level, &
debug_CPFEM, &
debug_levelBasic, &
debug_levelExtensive
use FEsolving, only: &
symmetricSolver, &
restartRead, &
modelName
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase
use constitutive, only: &
constitutive_state0
use crystallite, only: &
crystallite_F0, &
crystallite_Fp0, &
crystallite_Lp0, &
crystallite_dPdF0, &
crystallite_Tstar0_v, &
crystallite_localPlasticity
use homogenization, only: &
homogenization_sizeState, &
homogenization_state0
implicit none
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
integer(pInt) :: i,j,k,l,m
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
fixed bug in coordinate calculation when no average F was given (scaled wrong by ncp_elems**2) other changes: just polishing + some more comments
2013-03-27 17:58:55 +05:30
write(6,'(/,a)') ' <<<+- CPFEM init -+>>>'
write(6,'(a)') ' $Id$'
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
#include "compilation_info.f90"
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
! initialize stress and jacobian to zero
allocate(CPFEM_cs(6,mesh_maxNips,mesh_NcpElems)) ; CPFEM_cs = 0.0_pReal
allocate(CPFEM_dcsdE(6,6,mesh_maxNips,mesh_NcpElems)) ; CPFEM_dcsdE = 0.0_pReal
allocate(CPFEM_dcsdE_knownGood(6,6,mesh_maxNips,mesh_NcpElems)) ; CPFEM_dcsdE_knownGood = 0.0_pReal
! *** restore the last converged values of each essential variable from the binary file
if (restartRead) then
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
write(6,'(a)') '<< CPFEM >> restored state variables of last converged step from binary files'
!$OMP END CRITICAL (write2out)
endif
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_intFile(777,'recordedPhase',modelName,size(material_phase))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) material_phase
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedF',modelName,size(crystallite_F0))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) crystallite_F0
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedFp',modelName,size(crystallite_Fp0))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) crystallite_Fp0
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedLp',modelName,size(crystallite_Lp0))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) crystallite_Lp0
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergeddPdF',modelName,size(crystallite_dPdF0))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) crystallite_dPdF0
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedTstar',modelName,size(crystallite_Tstar0_v))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) crystallite_Tstar0_v
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedStateConst',modelName)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
m = 0_pInt
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)
m = m+1_pInt
read(777,rec=m) constitutive_state0(i,j,k)%p(l)
enddo
enddo; enddo; enddo
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergedStateHomog',modelName)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
m = 0_pInt
do k = 1,mesh_NcpElems; do j = 1,mesh_maxNips
do l = 1,homogenization_sizeState(j,k)
m = m+1_pInt
read(777,rec=m) homogenization_state0(j,k)%p(l)
enddo
enddo; enddo
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_read_realFile(777,'convergeddcsdE',modelName,size(CPFEM_dcsdE))
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
read (777,rec=1) CPFEM_dcsdE
close (777)
restartRead = .false.
endif
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a32,1x,6(i8,1x))') 'CPFEM_cs: ', shape(CPFEM_cs)
write(6,'(a32,1x,6(i8,1x))') 'CPFEM_dcsdE: ', shape(CPFEM_dcsdE)
write(6,'(a32,1x,6(i8,1x),/)') 'CPFEM_dcsdE_knownGood: ', shape(CPFEM_dcsdE_knownGood)
further improvements on implementing the plain mode
2013-04-16 22:37:27 +05:30
write(6,*) 'symmetricSolver: ', symmetricSolver
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
flush(6)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine CPFEM_init
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief perform initialization at first call, update variables and call the actual material model
!--------------------------------------------------------------------------------------------------
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
subroutine CPFEM_general(mode, parallelExecution, ffn, ffn1, temperature, dt, elFE, ip, cauchyStress, jacobian)
use numerics, only: &
defgradTolerance, &
iJacoStiffness
use debug, only: &
debug_level, &
debug_CPFEM, &
debug_levelBasic, &
debug_levelExtensive, &
debug_levelSelective, &
debug_e, &
debug_i, &
debug_stressMaxLocation, &
debug_stressMinLocation, &
debug_jacobianMaxLocation, &
debug_jacobianMinLocation, &
debug_stressMax, &
debug_stressMin, &
debug_jacobianMax, &
debug_jacobianMin
use FEsolving, only: &
outdatedFFN1, &
terminallyIll, &
cycleCounter, &
theInc, &
theTime, &
theDelta, &
FEsolving_execElem, &
FEsolving_execIP, &
restartWrite
use math, only: &
math_identity2nd, &
math_mul33x33, &
math_det33, &
math_transpose33, &
math_I3, &
math_Mandel3333to66, &
math_Mandel66to3333, &
math_Mandel33to6, &
math_Mandel6to33
use mesh, only: &
mesh_FEasCP, &
mesh_NcpElems, &
mesh_maxNips, &
mesh_element
use material, only: &
homogenization_maxNgrains, &
microstructure_elemhomo, &
some changes related to new state
2014-05-09 14:34:09 +05:30
#ifdef NEWSTATE
plasticState,&
#endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
material_phase
use constitutive, only: &
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
constitutive_state0, &
#ifdef NEWSTATE
mappingConstitutive,&
#endif
constitutive_state
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
use crystallite, only: &
crystallite_partionedF,&
crystallite_F0, &
crystallite_Fp0, &
crystallite_Fp, &
crystallite_Lp0, &
crystallite_Lp, &
crystallite_dPdF0, &
crystallite_dPdF, &
crystallite_Tstar0_v, &
crystallite_Tstar_v, &
crystallite_temperature
use homogenization, only: &
homogenization_sizeState, &
homogenization_state, &
homogenization_state0, &
materialpoint_F, &
materialpoint_F0, &
materialpoint_P, &
materialpoint_dPdF, &
materialpoint_results, &
materialpoint_sizeResults, &
materialpoint_stressAndItsTangent, &
materialpoint_postResults
use IO, only: &
IO_write_jobRealFile, &
IO_warning
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use DAMASK_interface
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
implicit none
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
integer(pInt), intent(in) :: elFE, & !< FE element number
ip !< integration point number
real(pReal), intent(in) :: temperature !< temperature
simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional) commented and cleaned up the marc interface. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! For marc simulations, run ./code/setup/setup_code.sh !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 15:38:00 +05:30
real(pReal), intent(in) :: dt !< time increment
real(pReal), dimension (3,3), intent(in) :: ffn, & !< deformation gradient for t=t0
ffn1 !< deformation gradient for t=t1
integer(pInt), intent(in) :: mode !< computation mode 1: regular computation plus aging of results
parallelExecution flag now passed as input variable by CPFEM_general; flag is set on the solver level (DAMASK_abaqus_exp.f, DAMASK_marc.f90, etc.)
2013-08-02 16:50:11 +05:30
logical, intent(in) :: parallelExecution !< flag indicating parallel computation of requested IPs
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional) commented and cleaned up the marc interface. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! For marc simulations, run ./code/setup/setup_code.sh !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 15:38:00 +05:30
real(pReal), dimension(6), intent(out), optional :: cauchyStress !< stress vector in Mandel notation
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
real(pReal), dimension(6,6), intent(out), optional :: jacobian !< jacobian in Mandel notation (Consistent tangent dcs/dE)
simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional) commented and cleaned up the marc interface. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! For marc simulations, run ./code/setup/setup_code.sh !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 15:38:00 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
real(pReal) J_inverse, & ! inverse of Jacobian
rnd
real(pReal), dimension (3,3) :: Kirchhoff, & ! Piola-Kirchhoff stress in Matrix notation
cauchyStress33 ! stress vector in Matrix notation
real(pReal), dimension (3,3,3,3) :: H_sym, &
H, &
jacobian3333 ! jacobian in Matrix notation
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
integer(pInt) elCP, & ! crystal plasticity element number
fixed bug in coordinate calculation when no average F was given (scaled wrong by ncp_elems**2) other changes: just polishing + some more comments
2013-03-27 17:58:55 +05:30
i, j, k, l, m, n
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
logical updateJaco ! flag indicating if JAcobian has to be updated
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
elCP = mesh_FEasCP('elem',elFE)
!if elCP = 0_pInt return ToDo: needed?
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
parallelExecution flag now passed as input variable by CPFEM_general; flag is set on the solver level (DAMASK_abaqus_exp.f, DAMASK_marc.f90, etc.)
2013-08-02 16:50:11 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt &
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
.and. elCP == debug_e .and. ip == debug_i) then
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
!$OMP CRITICAL (write2out)
write(6,'(/,a)') '#############################################'
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a1,a22,1x,i8,a13)') '#','element', elCP, '#'
write(6,'(a1,a22,1x,i8,a13)') '#','ip', ip, '#'
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
write(6,'(a1,a22,1x,f15.7,a6)') '#','theTime', theTime, '#'
write(6,'(a1,a22,1x,f15.7,a6)') '#','theDelta', theDelta, '#'
write(6,'(a1,a22,1x,i8,a13)') '#','theInc', theInc, '#'
write(6,'(a1,a22,1x,i8,a13)') '#','cycleCounter', cycleCounter, '#'
write(6,'(a1,a22,1x,i8,a13)') '#','computationMode',mode, '#'
write(6,'(a,/)') '#############################################'
flush (6)
!$OMP END CRITICAL (write2out)
endif
further improvements on implementing the plain mode
2013-04-16 22:37:27 +05:30
Changed order of blocks associated with different tasks: backup jacobian restore jacobian age results collect calc This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 1) backup jaco, 2) age results, 3) calc
2013-08-02 17:06:51 +05:30
!*** backup or restore jacobian
if (iand(mode, CPFEM_BACKUPJACOBIAN) /= 0_pInt) &
CPFEM_dcsde_knownGood = CPFEM_dcsde
if (iand(mode, CPFEM_RESTOREJACOBIAN) /= 0_pInt) &
CPFEM_dcsde = CPFEM_dcsde_knownGood
!*** age results and write restart data if requested
further improvements on implementing the plain mode
2013-04-16 22:37:27 +05:30
if (iand(mode, CPFEM_AGERESULTS) /= 0_pInt) then
crystallite_F0 = crystallite_partionedF ! crystallite deformation (_subF is perturbed...)
crystallite_Fp0 = crystallite_Fp ! crystallite plastic deformation
crystallite_Lp0 = crystallite_Lp ! crystallite plastic velocity
crystallite_dPdF0 = crystallite_dPdF ! crystallite stiffness
crystallite_Tstar0_v = crystallite_Tstar_v ! crystallite 2nd Piola Kirchhoff stress
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
forall ( i = 1:homogenization_maxNgrains, &
j = 1:mesh_maxNips, &
k = 1:mesh_NcpElems ) &
constitutive_state0(i,j,k)%p = constitutive_state(i,j,k)%p ! microstructure of crystallites
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
some changes related to new state
2014-05-09 14:34:09 +05:30
#ifdef NEWSTATE
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
forall ( i = 1:size(plasticState)) plasticState(i)%state0= plasticState(i)%state ! copy state in this lenghty way because A component cannot be an array if the encompassing structure is an array
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
write(6,'(a)') '<< CPFEM >> aging states'
if (debug_e <= mesh_NcpElems .and. debug_i <= mesh_maxNips) 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(mappingConstitutive(2,1,debug_i,debug_e))%state(:,mappingConstitutive(1,1,debug_i,debug_e))
endif
!$OMP END CRITICAL (write2out)
endif
some changes related to new state
2014-05-09 14:34:09 +05:30
#endif
fixed some errors in changes related to newstate
2014-05-12 06:14:44 +05:30
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
write(6,'(a)') '<< CPFEM >> aging states'
if (debug_e <= mesh_NcpElems .and. debug_i <= mesh_maxNips) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a,1x,i8,1x,i2,1x,i4,/,(12x,6(e20.8,1x)),/)') '<< CPFEM >> aged state of elFE ip grain',&
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
debug_e, debug_i, 1, constitutive_state(1,debug_i,debug_e)%p
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!$OMP END CRITICAL (write2out)
endif
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!$OMP PARALLEL DO
do k = 1,mesh_NcpElems
do j = 1,mesh_maxNips
if (homogenization_sizeState(j,k) > 0_pInt) &
homogenization_state0(j,k)%p = homogenization_state(j,k)%p ! internal state of homogenization scheme
enddo
enddo
!$OMP END PARALLEL DO
! * dump the last converged values of each essential variable to a binary file
if (restartWrite) then
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
write(6,'(a)') '<< CPFEM >> writing state variables of last converged step to binary files'
!$OMP END CRITICAL (write2out)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'recordedPhase',size(material_phase))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) material_phase
close (777)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedF',size(crystallite_F0))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) crystallite_F0
close (777)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedFp',size(crystallite_Fp0))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) crystallite_Fp0
close (777)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedLp',size(crystallite_Lp0))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) crystallite_Lp0
close (777)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergeddPdF',size(crystallite_dPdF0))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) crystallite_dPdF0
close (777)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedTstar',size(crystallite_Tstar0_v))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) crystallite_Tstar0_v
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedStateConst')
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
m = 0_pInt
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)
m = m+1_pInt
write(777,rec=m) constitutive_state0(i,j,k)%p(l)
enddo
enddo; enddo; enddo
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergedStateHomog')
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
m = 0_pInt
do k = 1,mesh_NcpElems; do j = 1,mesh_maxNips
do l = 1,homogenization_sizeState(j,k)
m = m+1_pInt
write(777,rec=m) homogenization_state0(j,k)%p(l)
enddo
enddo; enddo
close (777)
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
call IO_write_jobRealFile(777,'convergeddcsdE',size(CPFEM_dcsdE))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write (777,rec=1) CPFEM_dcsdE
close (777)
endif
endif
Changed order of blocks associated with different tasks: backup jacobian restore jacobian age results collect calc This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 1) backup jaco, 2) age results, 3) calc
2013-08-02 17:06:51 +05:30
!*** collection of FEM input with returning of randomize odd stress and jacobian
!* In case that no parallel execution is required, there is no need to collect FEM input
if (.not. parallelExecution) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
crystallite_temperature(ip,elCP) = temperature
materialpoint_F0(1:3,1:3,ip,elCP) = ffn
materialpoint_F(1:3,1:3,ip,elCP) = ffn1
Changed order of blocks associated with different tasks: backup jacobian restore jacobian age results collect calc This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 1) backup jaco, 2) age results, 3) calc
2013-08-02 17:06:51 +05:30
elseif (iand(mode, CPFEM_COLLECT) /= 0_pInt) then
call random_number(rnd)
if (rnd < 0.5_pReal) rnd = rnd - 1.0_pReal
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
crystallite_temperature(ip,elCP) = temperature
materialpoint_F0(1:3,1:3,ip,elCP) = ffn
materialpoint_F(1:3,1:3,ip,elCP) = ffn1
CPFEM_cs(1:6,ip,elCP) = rnd * CPFEM_odd_stress
CPFEM_dcsde(1:6,1:6,ip,elCP) = CPFEM_odd_jacobian * math_identity2nd(6)
Changed order of blocks associated with different tasks: backup jacobian restore jacobian age results collect calc This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 1) backup jaco, 2) age results, 3) calc
2013-08-02 17:06:51 +05:30
CPFEM_calc_done = .false.
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
endif
Changed order of blocks associated with different tasks: backup jacobian restore jacobian age results collect calc This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 1) backup jaco, 2) age results, 3) calc
2013-08-02 17:06:51 +05:30
!*** calculation of stress and jacobian
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
if (iand(mode, CPFEM_CALCRESULTS) /= 0_pInt) then
!*** deformation gradient outdated or any actual deformation gradient differs more than relevantStrain from the stored one
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
if (terminallyIll .or. outdatedFFN1 .or. any(abs(ffn1 - materialpoint_F(1:3,1:3,ip,elCP)) > defgradTolerance)) then
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
! if (.not. terminallyIll .and. .not. outdatedFFN1) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
if (any(abs(ffn1 - materialpoint_F(1:3,1:3,ip,elCP)) > defgradTolerance)) then
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a,1x,i8,1x,i2)') '<< CPFEM >> OUTDATED at elFE ip',elCP,ip
introduced dummy heat calculation, available as an crystallite output temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-19 00:27:28 +05:30
write(6,'(a,/,3(12x,3(f10.6,1x),/))') '<< CPFEM >> FFN1 old:',&
math_transpose33(materialpoint_F(1:3,1:3,ip,elCP))
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
write(6,'(a,/,3(12x,3(f10.6,1x),/))') '<< CPFEM >> FFN1 now:',math_transpose33(ffn1)
!$OMP END CRITICAL (write2out)
endif
outdatedFFN1 = .true.
endif
call random_number(rnd)
if (rnd < 0.5_pReal) rnd = rnd - 1.0_pReal
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
CPFEM_cs(1:6,ip,elCP) = rnd*CPFEM_odd_stress
CPFEM_dcsde(1:6,1:6,ip,elCP) = CPFEM_odd_jacobian*math_identity2nd(6)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!*** deformation gradient is not outdated
else
updateJaco = mod(cycleCounter,iJacoStiffness) == 0
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
!* no parallel computation, so we use just one single elFE and ip for computation
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
if (.not. parallelExecution) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
FEsolving_execElem(1) = elCP
FEsolving_execElem(2) = elCP
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
if (.not. microstructure_elemhomo(mesh_element(4,elCP)) .or. & ! calculate unless homogeneous
(microstructure_elemhomo(mesh_element(4,elCP)) .and. ip == 1_pInt)) then ! and then only first ip
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
FEsolving_execIP(1,elCP) = ip
FEsolving_execIP(2,elCP) = ip
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
!$OMP CRITICAL (write2out)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a,i8,1x,i2)') '<< CPFEM >> calculation for elFE ip ',elCP,ip
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
!$OMP END CRITICAL (write2out)
endif
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
call materialpoint_stressAndItsTangent(updateJaco, dt) ! calculate stress and its tangent
call materialpoint_postResults()
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!* parallel computation and calulation not yet done
elseif (.not. CPFEM_calc_done) then
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then
!$OMP CRITICAL (write2out)
introduced dummy heat calculation, available as an crystallite output temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-19 00:27:28 +05:30
write(6,'(a,i8,a,i8)') '<< CPFEM >> calculation for elements ',FEsolving_execElem(1),&
' to ',FEsolving_execElem(2)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!$OMP END CRITICAL (write2out)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
call materialpoint_stressAndItsTangent(updateJaco, dt) ! calculate stress and its tangent (parallel execution inside)
call materialpoint_postResults()
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
CPFEM_calc_done = .true.
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!* map stress and stiffness (or return odd values if terminally ill)
if ( terminallyIll ) then
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call random_number(rnd)
if (rnd < 0.5_pReal) rnd = rnd - 1.0_pReal
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
CPFEM_cs(1:6,ip,elCP) = rnd * CPFEM_odd_stress
CPFEM_dcsde(1:6,1:6,ip,elCP) = CPFEM_odd_jacobian * math_identity2nd(6)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
else
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
if (microstructure_elemhomo(mesh_element(4,elCP)) .and. ip > 1_pInt) then ! me homogenous? --> copy from first ip
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
materialpoint_P(1:3,1:3,ip,elCP) = materialpoint_P(1:3,1:3,1,elCP)
materialpoint_F(1:3,1:3,ip,elCP) = materialpoint_F(1:3,1:3,1,elCP)
materialpoint_dPdF(1:3,1:3,1:3,1:3,ip,elCP) = materialpoint_dPdF(1:3,1:3,1:3,1:3,1,elCP)
introduced dummy heat calculation, available as an crystallite output temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-19 00:27:28 +05:30
materialpoint_results(1:materialpoint_sizeResults,ip,elCP) = &
materialpoint_results(1:materialpoint_sizeResults,1,elCP)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
endif
! translate from P to CS
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
Kirchhoff = math_mul33x33(materialpoint_P(1:3,1:3,ip,elCP), math_transpose33(materialpoint_F(1:3,1:3,ip,elCP)))
J_inverse = 1.0_pReal / math_det33(materialpoint_F(1:3,1:3,ip,elCP))
CPFEM_cs(1:6,ip,elCP) = math_Mandel33to6(J_inverse * Kirchhoff)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
! translate from dP/dF to dCS/dE
H = 0.0_pReal
do i=1,3; do j=1,3; do k=1,3; do l=1,3; do m=1,3; do n=1,3
H(i,j,k,l) = H(i,j,k,l) + &
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
materialpoint_F(j,m,ip,elCP) * &
materialpoint_F(l,n,ip,elCP) * &
materialpoint_dPdF(i,m,k,n,ip,elCP) - &
math_I3(j,l) * materialpoint_F(i,m,ip,elCP) * materialpoint_P(k,m,ip,elCP) + &
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
0.5_pReal * (math_I3(i,k) * Kirchhoff(j,l) + math_I3(j,l) * Kirchhoff(i,k) + &
math_I3(i,l) * Kirchhoff(j,k) + math_I3(j,k) * Kirchhoff(i,l))
enddo; enddo; enddo; enddo; enddo; enddo
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
forall(i=1:3, j=1:3,k=1:3,l=1:3) &
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
H_sym(i,j,k,l) = 0.25_pReal * (H(i,j,k,l) + H(j,i,k,l) + H(i,j,l,k) + H(j,i,l,k))
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
CPFEM_dcsde(1:6,1:6,ip,elCP) = math_Mandel3333to66(J_inverse * H_sym)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
endif
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
!* remember extreme values of stress and jacobian
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
cauchyStress33 = math_Mandel6to33(CPFEM_cs(1:6,ip,elCP))
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
if (maxval(cauchyStress33) > debug_stressMax) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
debug_stressMaxLocation = [elCP, ip]
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
debug_stressMax = maxval(cauchyStress33)
endif
if (minval(cauchyStress33) < debug_stressMin) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
debug_stressMinLocation = [elCP, ip]
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
debug_stressMin = minval(cauchyStress33)
endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
jacobian3333 = math_Mandel66to3333(CPFEM_dcsdE(1:6,1:6,ip,elCP))
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
if (maxval(jacobian3333) > debug_jacobianMax) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
debug_jacobianMaxLocation = [elCP, ip]
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
debug_jacobianMax = maxval(jacobian3333)
endif
if (minval(jacobian3333) < debug_jacobianMin) then
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
debug_jacobianMinLocation = [elCP, ip]
moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location
2013-08-01 21:57:37 +05:30
debug_jacobianMin = minval(jacobian3333)
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!* report stress and stiffness
if ((iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) &
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
.and. ((debug_e == elCP .and. debug_i == ip) &
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
.or. .not. iand(debug_level(debug_CPFEM), debug_levelSelective) /= 0_pInt)) then
!$OMP CRITICAL (write2out)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
write(6,'(a,i8,1x,i2,/,12x,6(f10.3,1x)/)') '<< CPFEM >> stress/MPa at elFE ip ', &
elCP, ip, CPFEM_cs(1:6,ip,elCP)/1.0e6_pReal
write(6,'(a,i8,1x,i2,/,6(12x,6(f10.3,1x)/))') '<< CPFEM >> Jacobian/GPa at elFE ip ', &
elCP, ip, transpose(CPFEM_dcsdE(1:6,1:6,ip,elCP))/1.0e9_pReal
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
flush(6)
!$OMP END CRITICAL (write2out)
endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!*** warn if stiffness close to zero
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
if (all(abs(CPFEM_dcsdE(1:6,1:6,ip,elCP)) < 1e-10_pReal)) call IO_warning(601,elCP,ip)
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
!*** copy to output if required (FEM solver)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
if(present(cauchyStress)) cauchyStress = CPFEM_cs(1:6,ip,elCP)
if(present(jacobian)) jacobian = CPFEM_dcsdE(1:6,1:6,ip,elCP)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine CPFEM_general
forgot to add those two (rather crucial) files .. ;-)
2009-03-04 19:31:36 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
end module CPFEM