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

665 lines
28 KiB
Fortran
Raw Normal View History

some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @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 :: &
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 18:59:42 +05:30
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
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), dimension (:,:,:), allocatable, private :: &
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 18:59:42 +05:30
CPFEM_cs !< Cauchy stress
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), dimension (:,:,:,:), allocatable, private :: &
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 18:59:42 +05:30
CPFEM_dcsdE !< Cauchy stress tangent
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), dimension (:,:,:,:), allocatable, private :: &
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 18:59:42 +05:30
CPFEM_dcsdE_knownGood !< known good tangent
integer(pInt), public :: &
cycleCounter = 0_pInt, & !< needs description
theInc = -1_pInt, & !< needs description
lastLovl = 0_pInt, & !< lovl in previous call to marc hypela2
lastStep = 0_pInt !< kstep in previous call to abaqus umat
real(pReal), public :: &
theTime = 0.0_pReal, & !< needs description
theDelta = 0.0_pReal
logical, public :: &
outdatedFFN1 = .false., & !< needs description
lastIncConverged = .false., & !< needs description
outdatedByNewInc = .false. !< needs description
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
logical, public, protected :: &
calc done is used only internally
2016-01-17 19:45:38 +05:30
CPFEM_init_done = .false. !< remember whether init has been done already
logical, private :: &
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 18:59:42 +05:30
CPFEM_calc_done = .false. !< remember whether first ip has already calced the results
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
integer(pInt), parameter, public :: &
changed CPFEM_COLLECT to be 2^0=1 such that computation modes with collection step are always an odd number. (used to be 2^4 with no particular intrinsic meaning...)
2014-08-26 23:09:52 +05:30
CPFEM_COLLECT = 2_pInt**0_pInt, &
CPFEM_CALCRESULTS = 2_pInt**1_pInt, &
CPFEM_AGERESULTS = 2_pInt**2_pInt, &
CPFEM_BACKUPJACOBIAN = 2_pInt**3_pInt, &
CPFEM_RESTOREJACOBIAN = 2_pInt**4_pInt
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
public :: &
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
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
!--------------------------------------------------------------------------------------------------
remove deprecated spectral load case definition of temperature
2015-07-24 20:27:29 +05:30
subroutine CPFEM_initAll(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
more general name (should include parsing of debug and numerics)
2018-06-14 10:09:49 +05:30
use config, only: &
config_init
some doxygen corrections
2013-03-06 20:11:15 +05:30
use FEsolving, only: &
FE_init
use math, only: &
math_init
use mesh, only: &
mesh_init
use material, only: &
material_init
initializing lattice after materials allows to use the pre-parsed material.config
2018-06-03 00:29:30 +05:30
use lattice, only: &
lattice_init
some doxygen corrections
2013-03-06 20:11:15 +05:30
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
comments
2015-03-15 21:00:14 +05:30
integer(pInt), intent(in) :: el, & !< FE el number
ip !< FE integration point number
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
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
call DAMASK_interface_init ! Spectral and FEM interface to commandline
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
more general name (should include parsing of debug and numerics)
2018-06-14 10:09:49 +05:30
call config_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call math_init
call FE_init
calcmode only needed for Abaqus and MSC.Marc
2018-09-23 21:36:18 +05:30
call mesh_init(ip, el)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call lattice_init
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
call material_init
call constitutive_init
removed redundant crystallite_temperature/heat and materialpoint_heat and using temperature from the thermal module
2014-10-10 17:58:57 +05:30
call crystallite_init
remove deprecated spectral load case definition of temperature
2015-07-24 20:27:29 +05:30
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
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
compiler_options and compiler_version supported by new Intel compiler
2018-02-02 17:06:09 +05:30
#if __INTEL_COMPILER >= 1800
suppress warning when compiling with gfortran
2017-10-05 20:05:34 +05:30
use, intrinsic :: iso_fortran_env, only: &
compiler_version, &
compiler_options
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
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: &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
worldrank
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: &
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
material_phase, &
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
homogState, &
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
phase_plasticity, &
material.config is read in centrally moving data from material to config_material. use statements need to change. All aspects of reading from file will be removed from the individual modules
2018-06-10 21:31:52 +05:30
plasticState
more general name (should include parsing of debug and numerics)
2018-06-14 10:09:49 +05:30
use config, only: &
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
material_Nhomogenization
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use crystallite, only: &
crystallite_F0, &
crystallite_Fp0, &
crystallite_Lp0, &
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
crystallite_Fi0, &
crystallite_Li0, &
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
crystallite_dPdF0, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
crystallite_Tstar0_v
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
implicit none
removed unused variables
2014-09-23 02:04:42 +05:30
integer(pInt) :: k,l,m,ph,homog
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
character(len=1024) :: rankStr
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
mainProcess: if (worldrank == 0) then
write(6,'(/,a)') ' <<<+- CPFEM init -+>>>'
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
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
#include "compilation_info.f90"
fixing spectral cutback hiccup and multiple cleanups flush(6) at better places, added dedicated CPFEM_age subroutine, cleaned up cutback logic, fixed broken assignment of old timeinc, continueCalculation is now a logical, rearrnaged interfaces for utilities_constitutiveResponse and utilities_calculateRate, handling of stressBC more understandable, added more comments and explanations
2018-01-18 21:44:06 +05:30
flush(6)
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
endif mainProcess
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
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
! initialize stress and jacobian to zero
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
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
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
write(6,'(a)') '<< CPFEM >> restored state variables of last converged step from binary files'
flush(6)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
write(rankStr,'(a1,i0)')'_',worldrank
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_intFile(777,'recordedPhase'//trim(rankStr),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)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedF'//trim(rankStr),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)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedFp'//trim(rankStr),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)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedFi'//trim(rankStr),modelName,size(crystallite_Fi0))
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
read (777,rec=1) crystallite_Fi0
close (777)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedLp'//trim(rankStr),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)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedLi'//trim(rankStr),modelName,size(crystallite_Li0))
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
read (777,rec=1) crystallite_Li0
close (777)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergeddPdF'//trim(rankStr),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)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedTstar'//trim(rankStr),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)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedStateConst'//trim(rankStr),modelName)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
m = 0_pInt
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
readPlasticityInstances: do ph = 1_pInt, size(phase_plasticity)
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
do k = 1_pInt, plasticState(ph)%sizeState
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 23:20:01 +05:30
do l = 1, size(plasticState(ph)%state0(1,:))
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
m = m+1_pInt
read(777,rec=m) plasticState(ph)%state0(k,l)
enddo; enddo
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
enddo readPlasticityInstances
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
close (777)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_read_realFile(777,'convergedStateHomog'//trim(rankStr),modelName)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
m = 0_pInt
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
readHomogInstances: do homog = 1_pInt, material_Nhomogenization
do k = 1_pInt, homogState(homog)%sizeState
do l = 1, size(homogState(homog)%state0(1,:))
m = m+1_pInt
read(777,rec=m) homogState(homog)%state0(k,l)
enddo; enddo
enddo readHomogInstances
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
close (777)
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +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_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)
corrected linking: use PETSC MPI (when available) rather then system wide version. little polishing (e.g. debug output for new state)
2014-07-23 18:56:05 +05:30
write(6,'(a32,l1)') 'symmetricSolver: ', symmetricSolver
fixing spectral cutback hiccup and multiple cleanups flush(6) at better places, added dedicated CPFEM_age subroutine, cleaned up cutback logic, fixed broken assignment of old timeinc, continueCalculation is now a logical, rearrnaged interfaces for utilities_constitutiveResponse and utilities_calculateRate, handling of stressBC more understandable, added more comments and explanations
2018-01-18 21:44:06 +05:30
flush(6)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
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
!--------------------------------------------------------------------------------------------------
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
subroutine CPFEM_general(mode, parallelExecution, ffn, ffn1, temperature_inp, dt, elFE, ip, cauchyStress, 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
use numerics, only: &
defgradTolerance, &
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
iJacoStiffness, &
worldrank
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 debug, only: &
debug_level, &
debug_CPFEM, &
debug_levelBasic, &
debug_levelExtensive, &
debug_levelSelective, &
debug_stressMaxLocation, &
debug_stressMinLocation, &
debug_jacobianMaxLocation, &
debug_jacobianMinLocation, &
debug_stressMax, &
debug_stressMin, &
debug_jacobianMax, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
debug_jacobianMin, &
debug_e, &
debug_i
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 FEsolving, only: &
terminallyIll, &
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: &
microstructure_elemhomo, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState, &
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
sourceState, &
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
homogState, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
thermalState, &
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
damageState, &
exchanged mappingConstitutive array for more descriptive "phaseAt and phasememberAt"
2016-01-15 05:49:44 +05:30
phaseAt, phasememberAt, &
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
material_phase, &
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
phase_plasticity, &
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
temperature, &
thermalMapping, &
update temperature to FEM provided temperature only in case of thermal_condutcion, i.e. thermo-mechanical FEM simulation
2017-04-24 12:40:37 +05:30
thermal_type, &
THERMAL_conduction_ID, &
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
phase_Nsources, &
clearer code homogenization is defined per element, not per IP hence, use material_homogenizationAt instead of deprecated material_homog which pretends a dependency on the integration point
2018-10-04 10:09:03 +05:30
material_homogenizationAt
more general name (should include parsing of debug and numerics)
2018-06-14 10:09:49 +05:30
use config, only: &
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
material_Nhomogenization
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, &
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
crystallite_Fi0, &
crystallite_Fi, &
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_Lp0, &
crystallite_Lp, &
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
crystallite_Li0, &
crystallite_Li, &
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_dPdF0, &
crystallite_dPdF, &
crystallite_Tstar0_v, &
removed redundant crystallite_temperature/heat and materialpoint_heat and using temperature from the thermal module
2014-10-10 17:58:57 +05:30
crystallite_Tstar_v
fixed wrong line continuation
2014-09-20 01:08:59 +05:30
use homogenization, only: &
materialpoint_F, &
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_F0, &
materialpoint_P, &
materialpoint_dPdF, &
materialpoint_results, &
materialpoint_sizeResults, &
materialpoint_stressAndItsTangent, &
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
materialpoint_postResults
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 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
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
remove deprecated spectral load case definition of temperature
2015-07-24 20:27:29 +05:30
real(pReal), intent(in) :: temperature_inp !< temperature
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
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
real(pReal), dimension(6), intent(out) :: cauchyStress !< stress vector in Mandel notation
real(pReal), dimension(6,6), intent(out) :: 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
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
real(pReal) J_inverse, & ! inverse of Jacobian
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
rnd
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
real(pReal), dimension (3,3) :: Kirchhoff, & ! Piola-Kirchhoff stress in Matrix notation
cauchyStress33 ! stress vector in Matrix notation
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
real(pReal), dimension (3,3,3,3) :: H_sym, &
H, &
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
jacobian3333 ! jacobian in Matrix notation
integer(pInt) elCP, & ! crystal plasticity element number
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
i, j, k, l, m, n, ph, homog, mySource
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
logical updateJaco ! flag indicating if JAcobian has to be updated
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
character(len=1024) :: rankStr
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
write(6,'(/,a)') '#############################################'
write(6,'(a1,a22,1x,i8,a13)') '#','element', elCP, '#'
write(6,'(a1,a22,1x,i8,a13)') '#','ip', ip, '#'
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, '#'
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
if (terminallyIll) &
write(6,'(a,/)') '# --- terminallyIll --- #'
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
write(6,'(a,/)') '#############################################'; flush (6)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +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
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
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
crystallite_Fi0 = crystallite_Fi ! crystallite intermediate deformation
crystallite_Li0 = crystallite_Li ! crystallite intermediate velocity
further improvements on implementing the plain mode
2013-04-16 22:37:27 +05:30
crystallite_dPdF0 = crystallite_dPdF ! crystallite stiffness
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
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
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +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
do i = 1, size(sourceState)
do mySource = 1,phase_Nsources(i)
sourceState(i)%p(mySource)%state0 = sourceState(i)%p(mySource)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array
enddo; enddo
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 23:20:01 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
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, &
exchanged mappingConstitutive array for more descriptive "phaseAt and phasememberAt"
2016-01-15 05:49:44 +05:30
plasticState(phaseAt(1,debug_i,debug_e))%state(:,phasememberAt(1,debug_i,debug_e))
fixed small allocation flaws for new state
2014-05-12 18:30:37 +05:30
endif
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
endif
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
do homog = 1_pInt, material_Nhomogenization
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 22:32:23 +05:30
homogState (homog)%state0 = homogState (homog)%state
thermalState (homog)%state0 = thermalState (homog)%state
damageState (homog)%state0 = damageState (homog)%state
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
enddo
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
! * dump the last converged values of each essential variable to a binary file
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
if (restartWrite) then
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 23:20:01 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) &
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)') '<< CPFEM >> writing state variables of last converged step to binary files'
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
write(rankStr,'(a1,i0)')'_',worldrank
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'recordedPhase'//trim(rankStr),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
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedF'//trim(rankStr),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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedFp'//trim(rankStr),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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedFi'//trim(rankStr),size(crystallite_Fi0))
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
write (777,rec=1) crystallite_Fi0
close (777)
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedLp'//trim(rankStr),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)
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedLi'//trim(rankStr),size(crystallite_Li0))
changes related to intermediate configuration kinematics: - switched Fi and Li from state variables to crystallite variables - Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration - T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi - Updated analytic jacobian to take into account tangents wrt Fi - Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 18:39:00 +05:30
write (777,rec=1) crystallite_Li0
close (777)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergeddPdF'//trim(rankStr),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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedTstar'//trim(rankStr),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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedStateConst'//trim(rankStr))
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
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
writePlasticityInstances: do ph = 1_pInt, size(phase_plasticity)
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
do k = 1_pInt, plasticState(ph)%sizeState
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 23:20:01 +05:30
do l = 1, size(plasticState(ph)%state0(1,:))
fixed reading in of state for restarting
2014-07-24 14:08:52 +05:30
m = m+1_pInt
write(777,rec=m) plasticState(ph)%state0(k,l)
enddo; enddo
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
enddo writePlasticityInstances
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
close (777)
corrected linking: use PETSC MPI (when available) rather then system wide version. little polishing (e.g. debug output for new state)
2014-07-23 18:56:05 +05:30
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 21:36:19 +05:30
call IO_write_jobRealFile(777,'convergedStateHomog'//trim(rankStr))
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
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
writeHomogInstances: do homog = 1_pInt, material_Nhomogenization
do k = 1_pInt, homogState(homog)%sizeState
do l = 1, size(homogState(homog)%state0(1,:))
m = m+1_pInt
write(777,rec=m) homogState(homog)%state0(k,l)
enddo; enddo
enddo writeHomogInstances
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
close (777)
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +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,'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)
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +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
endif
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
endif ! results aging
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
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
!* If no parallel execution is required, there is no need to collect FEM input
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
if (.not. parallelExecution) then
fixed syntax errors
2017-04-28 16:09:01 +05:30
chosenThermal1: select case (thermal_type(mesh_element(3,elCP)))
case (THERMAL_conduction_ID) chosenThermal1
clearer code homogenization is defined per element, not per IP hence, use material_homogenizationAt instead of deprecated material_homog which pretends a dependency on the integration point
2018-10-04 10:09:03 +05:30
temperature(material_homogenizationAt(elCP))%p(thermalMapping(material_homogenizationAt(elCP))%p(ip,elCP)) = &
update temperature to FEM provided temperature only in case of thermal_condutcion, i.e. thermo-mechanical FEM simulation
2017-04-24 12:40:37 +05:30
temperature_inp
fixed syntax errors
2017-04-28 16:09:01 +05:30
end select chosenThermal1
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_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
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +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)
fixed syntax errors
2017-04-28 16:09:01 +05:30
chosenThermal2: select case (thermal_type(mesh_element(3,elCP)))
case (THERMAL_conduction_ID) chosenThermal2
clearer code homogenization is defined per element, not per IP hence, use material_homogenizationAt instead of deprecated material_homog which pretends a dependency on the integration point
2018-10-04 10:09:03 +05:30
temperature(material_homogenizationAt(elCP))%p(thermalMapping(material_homogenizationAt(elCP))%p(ip,elCP)) = &
update temperature to FEM provided temperature only in case of thermal_condutcion, i.e. thermo-mechanical FEM simulation
2017-04-24 12:40:37 +05:30
temperature_inp
fixed syntax errors
2017-04-28 16:09:01 +05:30
end select chosenThermal2
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_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
CPFEM_calc_done = .false.
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
endif ! collection
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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
!*** 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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
!*** deformation gradient outdated or any actual deformation gradient differs more than relevantStrain from the stored one
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
validCalculation: if (terminallyIll &
.or. outdatedFFN1 &
.or. any(abs(ffn1 - materialpoint_F(1:3,1:3,ip,elCP)) > defgradTolerance)) then
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
if (any(abs(ffn1 - materialpoint_F(1:3,1:3,ip,elCP)) > defgradTolerance)) then
if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
write(6,'(a,1x,i8,1x,i2)') '<< CPFEM >> OUTDATED at elFE ip',elFE,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)
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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
!*** deformation gradient is not outdated
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
else validCalculation
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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 (.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
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) &
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
write(6,'(a,i8,1x,i2)') '<< CPFEM >> calculation for elFE ip ',elFE,ip
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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
!* parallel computation and calulation not yet done
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
elseif (.not. CPFEM_calc_done) then
remove OMP Critical for output in regions not parallel
2014-05-16 19:31:27 +05:30
if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) &
write(6,'(a,i8,a,i8)') '<< CPFEM >> calculation for elements ',FEsolving_execElem(1),&
' to ',FEsolving_execElem(2)
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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +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
!* map stress and stiffness (or return odd values if terminally ill)
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
terminalIllness: 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)
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
else terminalIllness
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) = &
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +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
! 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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
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))
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
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)
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
endif terminalIllness
endif validCalculation
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
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
!* report stress and stiffness
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) &
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
write(6,'(a,i8,1x,i2,/,12x,6(f10.3,1x)/)') &
'<< CPFEM >> stress/MPa at elFE ip ', elFE, ip, CPFEM_cs(1:6,ip,elCP)*1.0e-6_pReal
write(6,'(a,i8,1x,i2,/,6(12x,6(f10.3,1x)/))') &
'<< CPFEM >> Jacobian/GPa at elFE ip ', elFE, ip, transpose(CPFEM_dcsdE(1:6,1:6,ip,elCP))*1.0e-9_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)
endif
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
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)
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 23:20:01 +05:30
!*** copy to output if using commercial FEM solver
disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM: command: /usr/bin/time -v OLD DAMASK terminated on: Date: 24/07/2014 Time: 12:03:58 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 12230.57 System time (seconds): 45.98 Percent of CPU this job got: 353% Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 243392 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 13839 Voluntary context switches: 452541 Involuntary context switches: 2233168 Swaps: 0 File system inputs: 0 File system outputs: 556112 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0 NEW DAMASK terminated on: Date: 24/07/2014 Time: 14:12:16 STOP 0 Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom" User time (seconds): 11948.80 System time (seconds): 37.28 Percent of CPU this job got: 365% Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97 Average shared text size (kbytes): 0 Average unshared data size (kbytes): 0 Average stack size (kbytes): 0 Average total size (kbytes): 0 Maximum resident set size (kbytes): 235280 Average resident set size (kbytes): 0 Major (requiring I/O) page faults: 0 Minor (reclaiming a frame) page faults: 14200 Voluntary context switches: 362509 Involuntary context switches: 1808383 Swaps: 0 File system inputs: 0 File system outputs: 556192 Socket messages sent: 0 Socket messages received: 0 Signals delivered: 0 Page size (bytes): 4096 Exit status: 0
2014-07-24 17:49:15 +05:30
cauchyStress = CPFEM_cs (1:6, ip,elCP)
jacobian = CPFEM_dcsdE(1:6,1:6,ip,elCP)
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
!*** remember extreme values of stress ...
cauchyStress33 = math_Mandel6to33(CPFEM_cs(1:6,ip,elCP))
if (maxval(cauchyStress33) > debug_stressMax) then
debug_stressMaxLocation = [elCP, ip]
debug_stressMax = maxval(cauchyStress33)
endif
if (minval(cauchyStress33) < debug_stressMin) then
debug_stressMinLocation = [elCP, ip]
debug_stressMin = minval(cauchyStress33)
endif
!*** ... and Jacobian
jacobian3333 = math_Mandel66to3333(CPFEM_dcsdE(1:6,1:6,ip,elCP))
if (maxval(jacobian3333) > debug_jacobianMax) then
debug_jacobianMaxLocation = [elCP, ip]
debug_jacobianMax = maxval(jacobian3333)
endif
if (minval(jacobian3333) < debug_jacobianMin) then
debug_jacobianMinLocation = [elCP, ip]
debug_jacobianMin = minval(jacobian3333)
endif
fixed reading in of state for restarting
2014-07-24 14:08:52 +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