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

697 lines
28 KiB
Fortran
Raw Normal View History

some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
! $Id$
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief CPFEM engine
!--------------------------------------------------------------------------------------------------
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
module CPFEM
some doxygen corrections
2013-03-06 20:11:15 +05:30
use prec, only: &
pReal, &
pInt
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
implicit none
introduced dummy heat calculation, available as an crystallite output temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-19 00:27:28 +05:30
private
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
#if defined(Marc4DAMASK) || defined(Abaqus)
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
real(pReal), parameter, private :: &
CPFEM_odd_stress = 1e15_pReal, & !< return value for stress in case of ping pong dummy cycle
CPFEM_odd_jacobian = 1e50_pReal !< return value for jacobian in case of ping pong dummy cycle
real(pReal), dimension (:,:,:), allocatable, private :: &
CPFEM_cs !< Cauchy stress
real(pReal), dimension (:,:,:,:), allocatable, private :: &
CPFEM_dcsdE !< Cauchy stress tangent
real(pReal), dimension (:,:,:,:), allocatable, private :: &
CPFEM_dcsdE_knownGood !< known good tangent
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
#endif
bugfix for last commit, too many variables declared as private
2013-10-23 15:59:38 +05:30
logical, public, protected :: &
CPFEM_init_done = .false., & !< remember whether init has been done already
CPFEM_calc_done = .false. !< remember whether first ip has already calced the results
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
!--------------------------------------------------------------------------------------------------
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
subroutine CPFEM_initAll(temperature,el,ip)
some doxygen corrections
2013-03-06 20:11:15 +05:30
use prec, only: &
prec_init
use numerics, only: &
numerics_init
use debug, only: &
debug_init
use FEsolving, only: &
FE_init
use math, only: &
math_init
use mesh, only: &
mesh_init
use lattice, only: &
lattice_init
use material, only: &
material_init
use constitutive, only: &
constitutive_init
use crystallite, only: &
crystallite_init
use homogenization, only: &
homogenization_init
use IO, only: &
IO_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use DAMASK_interface
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#ifdef FEM
use FEZoo, only: &
FEZoo_init
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
implicit none
comments
2015-03-15 21:00:14 +05:30
integer(pInt), intent(in) :: el, & !< FE el number
ip !< FE integration point number
real(pReal), intent(in) :: temperature !< temperature
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
!$OMP CRITICAL (init)
if (.not. CPFEM_init_done) then
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
#if defined(Spectral) || defined(FEM)
call DAMASK_interface_init ! Spectral and FEM interface to commandline
reordered initialization for spectral method, corrected bug of deformation BC parsing when prescribing velocity gradient resulting in wrong average deformation
2012-07-30 21:21:48 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call prec_init
call IO_init
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#ifdef FEM
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 FEZoo_init
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 06:08:29 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call numerics_init
call debug_init
call math_init
call FE_init
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
call mesh_init(ip, el) ! pass on coordinates to alter calcMode of first ip
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
call lattice_init
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
call material_init(temperature)
call constitutive_init(temperature)
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
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
call 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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
call DAMASK_interface_init ! Spectral solver and FEM init is already done
introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only. Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver. Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile. Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 21:40:21 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
CPFEM_init_done = .true.
endif
exchanged race condition prevention from waiting random amount of time to omp critical
2014-05-14 19:27:25 +05:30
!$OMP END CRITICAL (init)
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine CPFEM_initAll
Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Major changes: CPFEM.f90 => 1. Moving the initialization out of CPFEM_general into a separate subroutine, which is directly called by the hypela2 (Beware, the Abaqus version must also be modified in order to adapt with this change). 2. Restore primary state variables in CPFEM_init from binary files when requested (Marc flag: restart read). 3. Writing primary state variables into binary files (Marc flag: restart write). FEsolving.f90 => 1. Adding functions to recognize Marc restart flags: read and write and the corresponding restart file (parent job). 2. Change the initial value of cycleCounter = -1 in conjuction with the change made the ping-pong scheme homogenization_RGC.f90 => 1. Just syntax polishing. IO.f90 => 1. Adding functions/subroutines to open binary files for writing the primary state variables for restart purpose. mpie_cpfem_marc.f90 1. Modification of the general scheme for collection and calculation in order to accommodate the newly added restart feature.
2010-11-03 20:09:18 +05:30
some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief allocate the arrays defined in module CPFEM and initialize them
!--------------------------------------------------------------------------------------------------
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
subroutine CPFEM_init
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
use prec, only: &
pInt
use IO, only: &
introduced #EOF# as module wide parameter in IO.f90 renamed read/write binary file to read/write real/int removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 19:37:55 +05:30
IO_read_realFile,&
IO_read_intFile, &
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
IO_timeStamp, &
IO_error
use numerics, only: &
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: &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
symmetricSolver, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
restartRead, &
modelName
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
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, &
Finished applying newstate to Homogenization.
2014-09-19 23:29:06 +05:30
plasticState, &
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,'(a)') ' $Id$'
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"
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
#if defined(Marc4DAMASK) || defined(Abaqus)
! 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
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
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
! *** 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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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)
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
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
restartRead = .false.
endif
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
#if defined(Marc4DAMASK) || defined(Abaqus)
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
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
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
endif
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
#endif
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 17:18:29 +05:30
flush(6)
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine CPFEM_init
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
some doxygen corrections
2013-03-06 20:11:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief perform initialization at first call, update variables and call the actual material model
!--------------------------------------------------------------------------------------------------
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
#if defined(Marc4DAMASK) || defined(Abaqus)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
subroutine CPFEM_general(mode, parallelExecution, ffn, ffn1, temperature, dt, elFE, ip, cauchyStress, jacobian)
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
#else
subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip)
#endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
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, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
debug_stressMaxLocation, &
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, &
#endif
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: &
outdatedFFN1, &
terminallyIll, &
cycleCounter, &
theInc, &
theTime, &
theDelta, &
FEsolving_execElem, &
FEsolving_execIP, &
restartWrite
use math, only: &
math_identity2nd, &
math_mul33x33, &
math_det33, &
math_transpose33, &
math_I3, &
math_Mandel3333to66, &
math_Mandel66to3333, &
math_Mandel33to6, &
math_Mandel6to33
use mesh, only: &
mesh_FEasCP, &
mesh_NcpElems, &
mesh_maxNips, &
mesh_element
use material, only: &
microstructure_elemhomo, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState, &
damageState, &
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, &
age results and write restart data for vacancy concentration
2014-12-08 14:39:54 +05:30
vacancyState,&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
mappingConstitutive, &
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, &
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, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
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, &
materialpoint_dPdF, &
materialpoint_results, &
materialpoint_sizeResults, &
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 21:41:30 +05:30
#endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
materialpoint_stressAndItsTangent, &
temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use.
2014-10-13 23:38:54 +05:30
materialpoint_postResults, &
field_putFieldTemperature
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
real(pReal), intent(in) :: temperature !< temperature
simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional) commented and cleaned up the marc interface. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! For marc simulations, run ./code/setup/setup_code.sh !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 15:38:00 +05:30
real(pReal), intent(in) :: dt !< time increment
real(pReal), dimension (3,3), intent(in) :: ffn, & !< deformation gradient for t=t0
ffn1 !< deformation gradient for t=t1
integer(pInt), intent(in) :: mode !< computation mode 1: regular computation plus aging of results
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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
#else
logical, parameter :: parallelExecution = .true.
#endif
integer(pInt) elCP, & ! crystal plasticity element number
undid some marc/abaqus related changes
2014-09-23 13:49:15 +05:30
i, j, k, l, m, n, ph, homog
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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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)
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
#else
elCP = elFE
#endif
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
further improvements on implementing the plain mode
2013-04-16 22:37:27 +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
#if defined(Marc4DAMASK) || defined(Abaqus)
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
!*** backup jacobian
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
!*** restore jacobian
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_RESTOREJACOBIAN) /= 0_pInt) &
CPFEM_dcsde = CPFEM_dcsde_knownGood
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
#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
!*** 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
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
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
forall ( i = 1:size(damageState)) damageState(i)%state0 = damageState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array
forall ( i = 1:size(thermalState)) thermalState(i)%state0= thermalState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array
age results and write restart data for vacancy concentration
2014-12-08 14:39:54 +05:30
forall ( i = 1:size(vacancyState)) vacancyState(i)%state0= vacancyState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array
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, &
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
plasticState(mappingConstitutive(2,1,debug_i,debug_e))%state(:,mappingConstitutive(1,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
homogState(homog)%state0 = homogState(homog)%state
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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
#endif
corrected indentation level of parts of the code moved debugging output of stress and stiffness to different position in code
2013-08-01 21:40:56 +05:30
endif
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
temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use.
2014-10-13 23:38:54 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
call field_putFieldTemperature(ip,elCP,temperature)
#endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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
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)
temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use.
2014-10-13 23:38:54 +05:30
#endif
comments
2015-03-15 21:00:14 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use.
2014-10-13 23:38:54 +05:30
call field_putFieldTemperature(ip,elCP,temperature)
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
#endif
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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
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
#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
!*** 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)
change debug level for OUTDATED info from extensive to basic
2014-08-26 19:51:24 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
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
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
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
reporting of elFE (instead of elCP) when advertised as such.
2014-08-27 21:24:50 +05:30
#if defined(Marc4DAMASK) || defined(Abaqus)
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
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
#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
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
#if defined(Marc4DAMASK) || defined(Abaqus)
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
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
#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