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/constitutive.f90

1149 lines
58 KiB
Fortran
Raw Normal View History

added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
!--------------------------------------------------------------------------------------------------
some doxygen corrections
2013-03-06 20:11:15 +05:30
! $Id$
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
some doxygen corrections
2013-03-06 20:11:15 +05:30
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
!> @brief elasticity, plasticity, internal microstructure state
!--------------------------------------------------------------------------------------------------
module constitutive
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use prec, only: &
pInt, &
pReal, &
p_vec
implicit none
private
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifndef NEWSTATE
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
type(p_vec), public, dimension(:,:,:), allocatable :: &
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
constitutive_state0, & !< pointer array to microstructure at start of BVP inc
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
constitutive_partionedState0, & !< pointer array to microstructure at start of homogenization inc
constitutive_subState0, & !< pointer array to microstructure at start of crystallite inc
constitutive_state, & !< pointer array to current microstructure (end of converged time step)
constitutive_state_backup, & !< pointer array to backed up microstructure (end of converged time step)
constitutive_dotState, & !< pointer array to evolution of current microstructure
constitutive_deltaState, & !< pointer array to incremental change of current microstructure
constitutive_previousDotState,& !< pointer array to previous evolution of current microstructure
constitutive_previousDotState2,& !< pointer array to 2nd previous evolution of current microstructure
constitutive_dotState_backup, & !< pointer array to backed up evolution of current microstructure
constitutive_RK4dotState, & !< pointer array to evolution of microstructure defined by classical Runge-Kutta method
constitutive_aTolState !< pointer array to absolute state tolerance
type(p_vec), public, dimension(:,:,:,:), allocatable :: &
constitutive_RKCK45dotState !< pointer array to evolution of microstructure used by Cash-Karp Runge-Kutta method
integer(pInt), public, dimension(:,:,:), allocatable :: &
constitutive_sizeDotState, & !< size of dotState array
constitutive_sizeState, & !< size of state array per grain
constitutive_sizePostResults !< size of postResults array per grain
integer(pInt), private :: &
constitutive_maxSizeState
more changes related to new state
2014-05-22 20:54:12 +05:30
#else
integer(pInt), public, dimension(:,:,:), allocatable :: &
constitutive_sizePostResults !< size of postResults array per grain
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
#endif
more work on the new state
2014-05-27 20:16:03 +05:30
integer(pInt), public, protected :: &
constitutive_maxSizePostResults, &
constitutive_maxSizeDotState
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
public :: &
constitutive_init, &
constitutive_homogenizedC, &
constitutive_microstructure, &
constitutive_LpAndItsTangent, &
constitutive_TandItsTangent, &
constitutive_collectDotState, &
constitutive_collectDeltaState, &
constitutive_postResults
private :: &
constitutive_hooke_TandItsTangent
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
contains
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief allocates arrays pointing to array of the various constitutive modules
!--------------------------------------------------------------------------------------------------
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 constitutive_init
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
#ifdef HDF
use hdf5, only: &
HID_T
use IO, only : &
added a little bit more HDF5 functionality made abbreviations in configure options consistently capitals
2014-03-12 22:21:01 +05:30
HDF5_mappingConstitutive
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
#endif
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use debug, only: &
debug_level, &
debug_constitutive, &
debug_levelBasic
use numerics, only: &
numerics_integrator
use IO, only: &
IO_error, &
IO_open_file, &
IO_open_jobFile_stat, &
IO_write_jobFile, &
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_write_jobIntFile, &
reporting current time for most module inits
2013-02-25 22:04:59 +05:30
IO_timeStamp
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use mesh, only: &
mesh_maxNips, &
mesh_NcpElems, &
mesh_element, &
FE_Nips, &
FE_geomtype
use material, only: &
material_phase, &
material_Nphase, &
material_localFileExt, &
material_configFile, &
phase_name, &
phase_elasticity, &
phase_plasticity, &
phase_plasticityInstance, &
phase_Noutput, &
homogenization_Ngrains, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
homogenization_maxNgrains, &
ELASTICITY_HOOKE_ID, &
PLASTICITY_NONE_ID, &
PLASTICITY_J2_ID, &
PLASTICITY_PHENOPOWERLAW_ID, &
PLASTICITY_DISLOTWIN_ID, &
PLASTICITY_TITANMOD_ID, &
PLASTICITY_NONLOCAL_ID ,&
ELASTICITY_HOOKE_label, &
PLASTICITY_NONE_label, &
PLASTICITY_J2_label, &
PLASTICITY_PHENOPOWERLAW_label, &
PLASTICITY_DISLOTWIN_label, &
PLASTICITY_TITANMOD_label, &
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
#ifdef NEWSTATE
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
#endif
#if defined(HDF) || defined(NEWSTATE)
mappingConstitutive, &
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
PLASTICITY_NONLOCAL_label
new plasticity: none see material.config for necessary parameters (lattice_structure, elastic constants)
2012-07-03 16:46:38 +05:30
use constitutive_none
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
use constitutive_j2
use constitutive_phenopowerlaw
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
use constitutive_dislotwin
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
use constitutive_titanmod
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
use constitutive_nonlocal
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
integer(pInt), parameter :: FILEUNIT = 200_pInt
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
integer(pInt) :: &
g, & !< grain number
i, & !< integration point number
e, & !< element number
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
cMax, & !< maximum number of grains
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
iMax, & !< maximum number of integration points
eMax, & !< maximum number of elements
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
phase, &
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
s, &
more work on the new state
2014-05-27 20:16:03 +05:30
p, &
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
instance,&
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
myNgrains
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
integer(pInt), dimension(:,:), pointer :: thisSize
character(len=64), dimension(:,:), pointer :: thisOutput
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
character(len=32) :: outputName !< name of output, intermediate fix until HDF5 output is ready
fixed bug in nonlocal state init that was introduced in rev 2438 and used to reset states from other constitutions
2013-05-24 19:13:44 +05:30
logical :: knownPlasticity, nonlocalConstitutionPresent
nonlocalConstitutionPresent = .false.
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
! parse plasticities from config file
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
if (.not. IO_open_jobFile_stat(FILEUNIT,material_localFileExt)) & ! no local material configuration present...
call IO_open_file(FILEUNIT,material_configFile) ! ... open material.config file
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
if (any(phase_plasticity == PLASTICITY_NONE_ID)) call constitutive_none_init(FILEUNIT)
if (any(phase_plasticity == PLASTICITY_J2_ID)) call constitutive_j2_init(FILEUNIT)
if (any(phase_plasticity == PLASTICITY_PHENOPOWERLAW_ID)) call constitutive_phenopowerlaw_init(FILEUNIT)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
if (any(phase_plasticity == PLASTICITY_DISLOTWIN_ID)) call constitutive_dislotwin_init(FILEUNIT)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
if (any(phase_plasticity == PLASTICITY_TITANMOD_ID)) call constitutive_titanmod_init(FILEUNIT)
more changes related to new state
2014-05-22 20:54:12 +05:30
if (any(phase_plasticity == PLASTICITY_NONLOCAL_ID)) call constitutive_nonlocal_init(FILEUNIT)
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
close(FILEUNIT)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
init reporting of constitutive_*.f90 had less/to much spaces, renamed label to LABEL because it is a parameter. removed debug output of geom_fromEuclideanDistance.py
2013-05-28 23:01:55 +05:30
write(6,'(/,a)') ' <<<+- constitutive init -+>>>'
write(6,'(a)') ' $Id$'
fixed some warnings issued by gfortran (type conversion, not needed use statements)
2013-05-08 14:53:47 +05:30
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
fixed problem with gnu compiler: allocation of variables in constitutive_nonlocal_stateInit was not safe when not using any nonlocal constitution put header of debugging output in constitutive_init at the start of the function in order to make debugging easier
2012-10-09 18:04:57 +05:30
#include "compilation_info.f90"
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
! write description file for constitutive phase output
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
call IO_write_jobFile(FILEUNIT,'outputConstitutive')
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
do phase = 1_pInt,material_Nphase
instance = phase_plasticityInstance(phase) ! which instance of a plasticity is present phase
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
knownPlasticity = .true. ! assume valid
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
select case(phase_plasticity(phase)) ! split per constititution
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONE_ID)
outputName = PLASTICITY_NONE_label
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
thisOutput => null() ! constitutive_none_output
thisSize => null() ! constitutive_none_sizePostResult
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_J2_ID)
outputName = PLASTICITY_J2_label
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
thisOutput => constitutive_j2_output
thisSize => constitutive_j2_sizePostResult
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_PHENOPOWERLAW_ID)
homogenization: added enums and sourced allocation for RGC, some higher verbosity for isostrain
2013-12-18 12:58:01 +05:30
outputName = PLASTICITY_PHENOPOWERLAW_label
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
thisOutput => constitutive_phenopowerlaw_output
thisSize => constitutive_phenopowerlaw_sizePostResult
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
outputName = PLASTICITY_DISLOTWIN_label
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
thisOutput => constitutive_dislotwin_output
thisSize => constitutive_dislotwin_sizePostResult
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
outputName = PLASTICITY_TITANMOD_label
thisOutput => constitutive_titanmod_output
thisSize => constitutive_titanmod_sizePostResult
case (PLASTICITY_NONLOCAL_ID)
outputName = PLASTICITY_NONLOCAL_label
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
thisOutput => constitutive_nonlocal_output
thisSize => constitutive_nonlocal_sizePostResult
case default
knownPlasticity = .false.
end select
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
write(FILEUNIT,'(/,a,/)') '['//trim(phase_name(phase))//']'
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (knownPlasticity) then
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
write(FILEUNIT,'(a)') '(plasticity)'//char(9)//trim(outputName)
added thermal and damage effects to elastic hooke's law
2014-06-25 04:52:52 +05:30
if (phase_plasticity(phase) /= PLASTICITY_NONE_ID) then
do e = 1_pInt,phase_Noutput(phase)
write(FILEUNIT,'(a,i4)') trim(thisOutput(e,instance))//char(9),thisSize(e,instance)
enddo
endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
enddo
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
close(FILEUNIT)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
! allocation of states
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
cMax = homogenization_maxNgrains
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
iMax = mesh_maxNips
eMax = mesh_NcpElems
more changes related to new state
2014-05-22 20:54:12 +05:30
allocate(constitutive_sizePostResults(cMax,iMax,eMax), source=0_pInt)
#ifndef NEWSTATE
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
! lumped into new state
fixed previous incomplete commit and few other trivial changes in newstate
2014-06-11 22:22:18 +05:30
allocate(constitutive_state0(cMax,iMax,eMax))
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
allocate(constitutive_partionedState0(cMax,iMax,eMax))
allocate(constitutive_subState0(cMax,iMax,eMax))
allocate(constitutive_state(cMax,iMax,eMax))
allocate(constitutive_state_backup(cMax,iMax,eMax))
allocate(constitutive_dotState(cMax,iMax,eMax))
allocate(constitutive_deltaState(cMax,iMax,eMax))
allocate(constitutive_dotState_backup(cMax,iMax,eMax))
allocate(constitutive_aTolState(cMax,iMax,eMax))
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
! not needed anymore for new state
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
allocate(constitutive_sizeDotState(cMax,iMax,eMax), source=0_pInt)
allocate(constitutive_sizeState(cMax,iMax,eMax), source=0_pInt)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
allocate(constitutive_previousDotState(cMax,iMax,eMax))
allocate(constitutive_previousDotState2(cMax,iMax,eMax))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
allocate(constitutive_RK4dotState(cMax,iMax,eMax))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
allocate(constitutive_RKCK45dotState(6,cMax,iMax,eMax))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
ElemLoop:do e = 1_pInt,mesh_NcpElems ! loop over elements
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
myNgrains = homogenization_Ngrains(mesh_element(3,e))
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
IPloop:do i = 1_pInt,FE_Nips(FE_geomtype(mesh_element(2,e))) ! loop over IPs
GrainLoop:do g = 1_pInt,myNgrains ! loop over grains
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
select case(phase_elasticity(material_phase(g,i,e)))
case default ! so far no output for elasticity
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end select
some minor polishing on the way
2014-04-04 20:10:30 +05:30
phase = material_phase(g,i,e)
instance = phase_plasticityInstance(phase)
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
select case(phase_plasticity(material_phase(g,i,e)))
case (PLASTICITY_NONE_ID)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#ifndef NEWSTATE
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_state0(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_subState0(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_state(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_state_backup(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_aTolState(g,i,e)%p(constitutive_none_sizeState(instance)))
allocate(constitutive_dotState(g,i,e)%p(constitutive_none_sizeDotState(instance)))
allocate(constitutive_deltaState(g,i,e)%p(constitutive_none_sizeDotState(instance)))
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_none_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_none_sizeDotState(instance)))
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_none_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_none_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_none_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
enddo
endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_state0(g,i,e)%p = 0.0_pReal
constitutive_aTolState(g,i,e)%p = 1.0_pReal
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_sizeState(g,i,e) = 0_pInt
constitutive_sizeDotState(g,i,e) = 0_pInt
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_sizePostResults(g,i,e) = 0_pInt
more changes related to new state
2014-05-22 20:54:12 +05:30
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
case (PLASTICITY_J2_ID)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#ifndef NEWSTATE
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_state0(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_subState0(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_state(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_state_backup(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_aTolState(g,i,e)%p(constitutive_j2_sizeState(instance)))
allocate(constitutive_dotState(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
allocate(constitutive_deltaState(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_j2_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_j2_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
enddo
endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_state0(g,i,e)%p = constitutive_j2_stateInit(instance)
constitutive_aTolState(g,i,e)%p = constitutive_j2_aTolState(instance)
constitutive_sizeState(g,i,e) = constitutive_j2_sizeState(instance)
constitutive_sizeDotState(g,i,e) = constitutive_j2_sizeDotState(instance)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_sizePostResults(g,i,e) = constitutive_j2_sizePostResults(instance)
more changes related to new state
2014-05-22 20:54:12 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_PHENOPOWERLAW_ID)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#ifndef NEWSTATE
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_state0(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_subState0(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_state(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_state_backup(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_aTolState(g,i,e)%p(constitutive_phenopowerlaw_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_dotState(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_deltaState(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_phenopowerlaw_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
enddo
endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_state0(g,i,e)%p = constitutive_phenopowerlaw_stateInit(instance)
constitutive_aTolState(g,i,e)%p = constitutive_phenopowerlaw_aTolState(instance)
constitutive_sizeState(g,i,e) = constitutive_phenopowerlaw_sizeState(instance)
constitutive_sizeDotState(g,i,e) = constitutive_phenopowerlaw_sizeDotState(instance)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_sizePostResults(g,i,e) = constitutive_phenopowerlaw_sizePostResults(instance)
more changes related to new state
2014-05-22 20:54:12 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifndef NEWSTATE
allocate(constitutive_state0(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_subState0(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_state(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_state_backup(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_aTolState(g,i,e)%p(constitutive_dislotwin_sizeState(instance)),source=0.0_pReal)
allocate(constitutive_dotState(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_deltaState(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
new state related changes
2014-06-11 17:57:41 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_dislotwin_sizeDotState(instance)),source=0.0_pReal)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
enddo
endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_state0(g,i,e)%p = constitutive_dislotwin_stateInit(instance,material_phase(g,i,e))
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_aTolState(g,i,e)%p = constitutive_dislotwin_aTolState(instance)
constitutive_sizeState(g,i,e) = constitutive_dislotwin_sizeState(instance)
constitutive_sizeDotState(g,i,e) = constitutive_dislotwin_sizeDotState(instance)
new state related changes
2014-06-11 17:57:41 +05:30
#endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_sizePostResults(g,i,e) = constitutive_dislotwin_sizePostResults(instance)
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifndef NEWSTATE
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_state0(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_subState0(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_state(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_state_backup(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_aTolState(g,i,e)%p(constitutive_titanmod_sizeState(instance)))
allocate(constitutive_dotState(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
allocate(constitutive_deltaState(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
if (any(numerics_integrator == 1_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_titanmod_sizeDotState(instance)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
enddo
endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_state0(g,i,e)%p = constitutive_titanmod_stateInit(instance,material_phase(g,i,e))
constitutive_aTolState(g,i,e)%p = constitutive_titanmod_aTolState(instance)
constitutive_sizeState(g,i,e) = constitutive_titanmod_sizeState(instance)
constitutive_sizeDotState(g,i,e) = constitutive_titanmod_sizeDotState(instance)
new state related changes
2014-06-11 17:57:41 +05:30
#endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_sizePostResults(g,i,e) = constitutive_titanmod_sizePostResults(instance)
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
fixed bug in nonlocal state init that was introduced in rev 2438 and used to reset states from other constitutions
2013-05-24 19:13:44 +05:30
nonlocalConstitutionPresent = .true.
removed few bugs from nonlocal (newstate)
2014-06-24 14:54:59 +05:30
#ifdef NEWSTATE
plasticState(mappingConstitutive(2,g,i,e))%nonlocal = .true.
#endif
added check if nGrains==1 when non-local plasticity is used
2013-02-04 20:04:01 +05:30
if(myNgrains/=1_pInt) call IO_error(252_pInt, e,i,g)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifndef NEWSTATE
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_state0(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_partionedState0(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_subState0(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_state(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_state_backup(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_aTolState(g,i,e)%p(constitutive_nonlocal_sizeState(instance)))
allocate(constitutive_dotState(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
allocate(constitutive_deltaState(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
allocate(constitutive_dotState_backup(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (any(numerics_integrator == 1_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_previousDotState(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
allocate(constitutive_previousDotState2(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 4_pInt)) then
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RK4dotState(g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
if (any(numerics_integrator == 5_pInt)) then
do s = 1_pInt,6_pInt
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
allocate(constitutive_RKCK45dotState(s,g,i,e)%p(constitutive_nonlocal_sizeDotState(instance)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
enddo
endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_aTolState(g,i,e)%p = constitutive_nonlocal_aTolState(instance)
constitutive_sizeState(g,i,e) = constitutive_nonlocal_sizeState(instance)
constitutive_sizeDotState(g,i,e) = constitutive_nonlocal_sizeDotState(instance)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
constitutive_sizePostResults(g,i,e) = constitutive_nonlocal_sizePostResults(instance)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
end select
more changes related to new state
2014-05-22 20:54:12 +05:30
enddo GrainLoop
enddo IPloop
enddo ElemLoop
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
if (nonlocalConstitutionPresent) &
#ifdef NEWSTATE
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
call constitutive_nonlocal_stateInit()
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#else
fixed bug in nonlocal state init that was introduced in rev 2438 and used to reset states from other constitutions
2013-05-24 19:13:44 +05:30
call constitutive_nonlocal_stateInit(constitutive_state0(1,1:iMax,1:eMax))
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#endif
removed few bugs from nonlocal (newstate)
2014-06-24 14:54:59 +05:30
#ifdef NEWSTATE
do e = 1_pInt,mesh_NcpElems ! loop over elements
myNgrains = homogenization_Ngrains(mesh_element(3,e))
forall(i = 1_pInt:FE_Nips(FE_geomtype(mesh_element(2,e))), g = 1_pInt:myNgrains)
plasticState(mappingConstitutive(2,g,i,e))%partionedState0(:,mappingConstitutive(1,g,i,e)) = &
plasticState(mappingConstitutive(2,g,i,e))%State0(:,mappingConstitutive(1,g,i,e)) ! need to be defined for first call of constitutive_microstructure in crystallite_init
plasticState(mappingConstitutive(2,g,i,e))%State(:,mappingConstitutive(1,g,i,e)) = &
plasticState(mappingConstitutive(2,g,i,e))%State0(:,mappingConstitutive(1,g,i,e)) ! need to be defined for first call of constitutive_microstructure in crystallite_init
endforall
enddo
#else
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
do e = 1_pInt,mesh_NcpElems ! loop over elements
myNgrains = homogenization_Ngrains(mesh_element(3,e))
forall(i = 1_pInt:FE_Nips(FE_geomtype(mesh_element(2,e))), g = 1_pInt:myNgrains)
constitutive_partionedState0(g,i,e)%p = constitutive_state0(g,i,e)%p
constitutive_state(g,i,e)%p = constitutive_state0(g,i,e)%p ! need to be defined for first call of constitutive_microstructure in crystallite_init
endforall
enddo
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
removed few bugs from nonlocal (newstate)
2014-06-24 14:54:59 +05:30
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
#ifdef HDF
small changes regarding restart and HDF5
2014-03-26 14:11:45 +05:30
call HDF5_mappingConstitutive(mappingConstitutive)
modified reading in of values (realizes, if phase is not used in the model at all), started to include HDF5 output and new state
2014-04-15 15:39:20 +05:30
do phase = 1_pInt,material_Nphase
instance = phase_plasticityInstance(phase) ! which instance of a plasticity is present phase
select case(phase_plasticity(phase)) ! split per constititution
case (PLASTICITY_NONE_ID)
case (PLASTICITY_J2_ID)
end select
enddo
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
#endif
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifndef NEWSTATE
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
! write out state size file
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_jobIntFile(777,'sizeStateConst', size(constitutive_sizeState))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
write (777,rec=1) constitutive_sizeState
close(777)
changed line endings to unix style
2009-03-06 15:32:36 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
! report
constitutive_maxSizeState = maxval(constitutive_sizeState)
constitutive_maxSizeDotState = maxval(constitutive_sizeDotState)
constitutive_maxSizePostResults = maxval(constitutive_sizePostResults)
if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 0_pInt) then
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
write(6,'(a32,1x,7(i8,1x))') 'constitutive_state0: ', shape(constitutive_state0)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_partionedState0: ', shape(constitutive_partionedState0)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_subState0: ', shape(constitutive_subState0)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_state: ', shape(constitutive_state)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_aTolState: ', shape(constitutive_aTolState)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_dotState: ', shape(constitutive_dotState)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_deltaState: ', shape(constitutive_deltaState)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_sizeState: ', shape(constitutive_sizeState)
write(6,'(a32,1x,7(i8,1x))') 'constitutive_sizeDotState: ', shape(constitutive_sizeDotState)
write(6,'(a32,1x,7(i8,1x),/)') 'constitutive_sizePostResults: ', shape(constitutive_sizePostResults)
write(6,'(a32,1x,7(i8,1x))') 'maxSizeState: ', constitutive_maxSizeState
write(6,'(a32,1x,7(i8,1x))') 'maxSizeDotState: ', constitutive_maxSizeDotState
write(6,'(a32,1x,7(i8,1x))') 'maxSizePostResults: ', constitutive_maxSizePostResults
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
endif
flush(6)
fixed missing sizeState initialization causing problem with output (only for new state)
2014-05-27 17:40:16 +05:30
#else
constitutive_maxSizePostResults = maxval(constitutive_sizePostResults)
more work on the new state
2014-05-27 20:16:03 +05:30
constitutive_maxSizeDotState = 0_pInt
do p = 1, size(plasticState)
constitutive_maxSizeDotState = max(constitutive_maxSizeDotState, plasticState(p)%sizeDotState)
enddo
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end subroutine constitutive_init
changed line endings to unix style
2009-03-06 15:32:36 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns the homogenize elasticity matrix
!--------------------------------------------------------------------------------------------------
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
pure function constitutive_homogenizedC(ipc,ip,el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use material, only: &
phase_plasticity, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
material_phase, &
PLASTICITY_TITANMOD_ID, &
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
plasticState,&
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
PLASTICITY_DISLOTWIN_ID
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
use constitutive_titanmod, only: &
constitutive_titanmod_homogenizedC
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
use constitutive_dislotwin, only: &
constitutive_dislotwin_homogenizedC
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
use lattice, only: &
lattice_C66
forgot to add those two (rather crucial) files .. ;-)
2009-03-04 19:31:36 +05:30
implicit none
real(pReal), dimension(6,6) :: constitutive_homogenizedC
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
integer(pInt), intent(in) :: &
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
changed line endings to unix style
2009-03-06 15:32:36 +05:30
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
constitutive_homogenizedC = constitutive_dislotwin_homogenizedC &
(plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
ipc,ip,el)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#else
constitutive_homogenizedC = constitutive_dislotwin_homogenizedC &
(constitutive_state(ipc,ip,el),ipc,ip,el)
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifdef NEWSTATE
constitutive_homogenizedC = constitutive_titanmod_homogenizedC &
(plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_homogenizedC = constitutive_titanmod_homogenizedC(constitutive_state(ipc,ip,el), &
ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case default
constitutive_homogenizedC = lattice_C66(1:6,1:6,material_phase(ipc,ip,el))
now with first draft of nonlocal constitutive law debugging memory leak closed debugging counters corrected center of gravity stored in mesh state updated is now split into a collecting loop and an execution updateState and updateTemperature fill sequentially separate logicals and evaluate afterwards to converged added 3x3 transposition function, norm for 3x1 matrix and 33x3 matrix multiplication in math non-converged crystallite triggers materialpoint cutback (used to respond elastically) non-converged materialpoint raises terminal illness which in turn renders whole FE increment useless by means of odd stress/stiffness and thus waits for FE cutback
2009-08-11 22:01:57 +05:30
forgot to add those two (rather crucial) files .. ;-)
2009-03-04 19:31:36 +05:30
end select
changed line endings to unix style
2009-03-06 15:32:36 +05:30
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end function constitutive_homogenizedC
changed line endings to unix style
2009-03-06 15:32:36 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief calls microstructure function of the different constitutive models
!--------------------------------------------------------------------------------------------------
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 constitutive_microstructure(temperature, Fe, Fp, ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use material, only: &
phase_plasticity, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
material_phase, &
PLASTICITY_DISLOTWIN_ID, &
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
PLASTICITY_TITANMOD_ID, &
PLASTICITY_NONLOCAL_ID
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_titanmod, only: &
constitutive_titanmod_microstructure
use constitutive_nonlocal, only: &
constitutive_nonlocal_microstructure
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
use constitutive_dislotwin, only: &
constitutive_dislotwin_microstructure
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in) :: &
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
temperature
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(3,3) :: &
Fe, & !< elastic deformation gradient
Fp !< plastic deformation gradient
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
call constitutive_dislotwin_microstructure(temperature, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
ipc,ip,el)
#else
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
call constitutive_dislotwin_microstructure(temperature,constitutive_state(ipc,ip,el), &
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
ipc,ip,el)
#endif
new state related changes
2014-06-11 17:57:41 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifdef NEWSTATE
call constitutive_titanmod_microstructure(temperature, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_titanmod_microstructure(temperature,constitutive_state(ipc,ip,el), &
ipc,ip,el)
new state related changes
2014-06-11 17:57:41 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
call constitutive_nonlocal_microstructure(Fe,Fp,ipc,ip,el)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#else
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 constitutive_nonlocal_microstructure(constitutive_state,Fe,Fp,ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
end select
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end subroutine constitutive_microstructure
openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules: DON'T use implicit array subscripts: example: real, dimension(3,3) :: A,B A(:,2) = B(:,1) <--- DON'T USE A(1:3,2) = B(1:3,1) <--- BETTER USE In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks. Enclose all write statements with the following: !$OMP CRITICAL (write2out) <your write statement> !$OMP END CRITICAL (write2out) Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 16:16:17 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
!> @brief contains the constitutive equation for calculating the velocity gradient
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +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
subroutine constitutive_LpAndItsTangent(Lp, dLp_dTstar, Tstar_v, temperature, ipc, ip, el)
use math, only: &
math_identity2nd
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use material, only: &
phase_plasticity, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
material_phase, &
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
plasticState,&
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
PLASTICITY_NONE_ID, &
PLASTICITY_J2_ID, &
PLASTICITY_PHENOPOWERLAW_ID, &
PLASTICITY_DISLOTWIN_ID, &
PLASTICITY_TITANMOD_ID, &
PLASTICITY_NONLOCAL_ID
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_j2, only: &
constitutive_j2_LpAndItsTangent
use constitutive_phenopowerlaw, only: &
more changes related to new state
2014-05-22 20:54:12 +05:30
constitutive_phenopowerlaw_LpAndItsTangent
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_dislotwin, only: &
constitutive_dislotwin_LpAndItsTangent
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
use constitutive_titanmod, only: &
constitutive_titanmod_LpAndItsTangent
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_nonlocal, only: &
constitutive_nonlocal_LpAndItsTangent
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in) :: &
Temperature
real(pReal), intent(in), dimension(6) :: &
Tstar_v !< 2nd Piola-Kirchhoff stress
real(pReal), intent(out), dimension(3,3) :: &
Lp !< plastic velocity gradient
real(pReal), intent(out), dimension(9,9) :: &
dLp_dTstar !< derivative of Lp with respect to Tstar (4th-order tensor)
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONE_ID)
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
Lp = 0.0_pReal
dLp_dTstar = math_identity2nd(9)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_J2_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
call constitutive_j2_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_j2_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
constitutive_state(ipc,ip,el),ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_PHENOPOWERLAW_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
call constitutive_phenopowerlaw_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_phenopowerlaw_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
constitutive_state(ipc,ip,el),ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
call constitutive_dislotwin_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v,temperature, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_dislotwin_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
temperature,constitutive_state(ipc,ip,el),ipc,ip,el)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#endif
new state related changes
2014-06-11 17:57:41 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifdef NEWSTATE
call constitutive_titanmod_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v,temperature, &
plasticState(mappingConstitutive(2,ipc,ip,el))%state(:,mappingConstitutive(1,ipc,ip,el)), &
ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_titanmod_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v, &
temperature,constitutive_state(ipc,ip,el),ipc,ip,el)
new state related changes
2014-06-11 17:57:41 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
call constitutive_nonlocal_LpAndItsTangent(Lp, dLp_dTstar, Tstar_v, &
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
temperature, ipc,ip,el)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
call constitutive_nonlocal_LpAndItsTangent(Lp, dLp_dTstar, Tstar_v, &
temperature, constitutive_state(ipc,ip,el), ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end select
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end subroutine constitutive_LpAndItsTangent
changed line endings to unix style
2009-03-06 15:32:36 +05:30
openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules: DON'T use implicit array subscripts: example: real, dimension(3,3) :: A,B A(:,2) = B(:,1) <--- DON'T USE A(1:3,2) = B(1:3,1) <--- BETTER USE In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks. Enclose all write statements with the following: !$OMP CRITICAL (write2out) <your write statement> !$OMP END CRITICAL (write2out) Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 16:16:17 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns the 2nd Piola-Kirchhoff stress tensor and its tangent with respect to
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!> the elastic deformation gradient depending on the selected elastic law (so far no case switch
!! because only hooke is implemented
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
pure subroutine constitutive_TandItsTangent(T, dT_dFe, Fe, ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(3,3) :: &
Fe !< elastic deformation gradient
real(pReal), intent(out), dimension(3,3) :: &
T !< 2nd Piola-Kirchhoff stress tensor
real(pReal), intent(out), dimension(3,3,3,3) :: &
dT_dFe !< derivative of 2nd P-K stress with respect to elastic deformation gradient
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
call constitutive_hooke_TandItsTangent(T, dT_dFe, Fe, ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end subroutine constitutive_TandItsTangent
shifted location of elasticity part closer to plasticity... (constitutive has been shaken back and forth quite a bit today ;-] )
2012-03-15 15:21:33 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns the 2nd Piola-Kirchhoff stress tensor and its tangent with respect to
improved comments
2013-10-14 20:05:41 +05:30
!> the elastic deformation gradient using hookes law
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
pure subroutine constitutive_hooke_TandItsTangent(T, dT_dFe, Fe, ipc, ip, el)
removed unnecsessary omp statements, removed one transpose in constitutitve
2013-01-08 16:39:20 +05:30
use math, only : &
Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged. In contrast, the analytical jacobian will probably be affected by this change! @Pratheek: Can you check with me how this can be fixed?
2014-01-22 14:08:13 +05:30
math_mul3x3, &
removed unnecsessary omp statements, removed one transpose in constitutitve
2013-01-08 16:39:20 +05:30
math_mul33x33, &
minor changes on buggy latest commit (wrong use statement)
2013-01-09 03:24:25 +05:30
math_mul3333xx33, &
removed unnecsessary omp statements, removed one transpose in constitutitve
2013-01-08 16:39:20 +05:30
math_Mandel66to3333, &
minor changes on buggy latest commit (wrong use statement)
2013-01-09 03:24:25 +05:30
math_transpose33, &
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
MATH_I3
added thermal and damage effects to elastic hooke's law
2014-06-25 04:52:52 +05:30
#ifdef NEWSTATE
use material, only: &
mappingConstitutive, &
damageState, &
phase_damage, &
DAMAGE_gradient_ID, &
thermalState, &
phase_thermal, &
THERMAL_conduction_ID, &
THERMAL_adiabatic_ID
use lattice, only: &
lattice_referenceTemperature, &
lattice_thermalExpansion33
#endif
fixed bug in crystallite_FPI: stateDamper always has to be defined for each grain added some OMP FLUSH statements were necessary replaced openmp do by forall construct where possible; this is much safer and perhaps even as fast for small loops
2012-11-07 21:13:29 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(3,3) :: &
Fe !< elastic deformation gradient
real(pReal), intent(out), dimension(3,3) :: &
T !< 2nd Piola-Kirchhoff stress tensor
real(pReal), intent(out), dimension(3,3,3,3) :: &
dT_dFe !< dT/dFe
Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged. In contrast, the analytical jacobian will probably be affected by this change! @Pratheek: Can you check with me how this can be fixed?
2014-01-22 14:08:13 +05:30
integer(pInt) :: i, j, k, l
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), dimension(3,3) :: FeT
real(pReal), dimension(3,3,3,3) :: C
added thermal and damage effects to elastic hooke's law
2014-06-25 04:52:52 +05:30
real(pReal) :: damage
integer(pInt) :: phase, constituent
shifted location of elasticity part closer to plasticity... (constitutive has been shaken back and forth quite a bit today ;-] )
2012-03-15 15:21:33 +05:30
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
C = math_Mandel66to3333(constitutive_homogenizedC(ipc,ip,el))
shifted location of elasticity part closer to plasticity... (constitutive has been shaken back and forth quite a bit today ;-] )
2012-03-15 15:21:33 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
FeT = math_transpose33(Fe)
added thermal and damage effects to elastic hooke's law
2014-06-25 04:52:52 +05:30
T = 0.5_pReal*math_mul3333xx33(C,math_mul33x33(FeT,Fe)-MATH_I3)
Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged. In contrast, the analytical jacobian will probably be affected by this change! @Pratheek: Can you check with me how this can be fixed?
2014-01-22 14:08:13 +05:30
dT_dFe = 0.0_pReal
forall (i=1_pInt:3_pInt, j=1_pInt:3_pInt, k=1_pInt:3_pInt, l=1_pInt:3_pInt) &
dT_dFe(i,j,k,l) = math_mul3x3(C(i,j,l,1:3),Fe(k,1:3)) ! dT*_ij/dFe_kl
shifted location of elasticity part closer to plasticity... (constitutive has been shaken back and forth quite a bit today ;-] )
2012-03-15 15:21:33 +05:30
added thermal and damage effects to elastic hooke's law
2014-06-25 04:52:52 +05:30
#ifdef NEWSTATE
phase = mappingConstitutive(2,ipc,ip,el)
constituent = mappingConstitutive(1,ipc,ip,el)
select case (phase_damage(phase))
case (DAMAGE_gradient_ID)
damage = damageState(phase)%state(3,constituent) &
* damageState(phase)%state(3,constituent)
T = damage*T
dT_dFe = damage*dT_dFe
end select
select case (phase_thermal(phase))
case (THERMAL_conduction_ID)
T = T - math_mul3333xx33(C, (thermalState(phase)%state(2,constituent) - &
lattice_referenceTemperature(phase)) &
* lattice_thermalExpansion33(1:3,1:3,phase))
case (THERMAL_adiabatic_ID)
T = T - math_mul3333xx33(C, (thermalState(phase)%state(1,constituent) - &
lattice_referenceTemperature(phase)) &
* lattice_thermalExpansion33(1:3,1:3,phase))
end select
#endif
shifted location of elasticity part closer to plasticity... (constitutive has been shaken back and forth quite a bit today ;-] )
2012-03-15 15:21:33 +05:30
end subroutine constitutive_hooke_TandItsTangent
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief contains the constitutive equation for calculating the rate of change of microstructure
!--------------------------------------------------------------------------------------------------
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
subroutine constitutive_collectDotState(Tstar_v, FeArray, FpArray, Temperature, subdt, subfracArray,&
ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use prec, only: &
pLongInt
use debug, only: &
debug_cumDotStateCalls, &
debug_cumDotStateTicks, &
debug_level, &
debug_constitutive, &
debug_levelBasic
use mesh, only: &
mesh_NcpElems, &
fixed some warnings issued by gfortran (type conversion, not needed use statements)
2013-05-08 14:53:47 +05:30
mesh_maxNips
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use material, only: &
phase_plasticity, &
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
material_phase, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
homogenization_maxNgrains, &
PLASTICITY_NONE_ID, &
PLASTICITY_J2_ID, &
PLASTICITY_PHENOPOWERLAW_ID, &
PLASTICITY_DISLOTWIN_ID, &
PLASTICITY_TITANMOD_ID, &
PLASTICITY_NONLOCAL_ID
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_j2, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_j2_dotState
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_phenopowerlaw, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_phenopowerlaw_dotState
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_dislotwin, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_dislotwin_dotState
use constitutive_titanmod, only: &
constitutive_titanmod_dotState
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_nonlocal, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_nonlocal_dotState
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in) :: &
Temperature, &
subdt !< timestep
real(pReal), intent(in), dimension(homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: &
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
subfracArray !< subfraction of timestep
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(3,3,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: &
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
FeArray, & !< elastic deformation gradient
FpArray !< plastic deformation gradient
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(6) :: &
Tstar_v !< 2nd Piola Kirchhoff stress tensor (Mandel)
integer(pLongInt) :: &
tick, tock, &
tickrate, &
maxticks
new plasticity: none see material.config for necessary parameters (lattice_structure, elastic constants)
2012-07-03 16:46:38 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (iand(debug_level(debug_constitutive), debug_levelBasic) /= 0_pInt) &
call system_clock(count=tick,count_rate=tickrate,count_max=maxticks)
constitutive_nonlocal: - corrected flux term - multiplication is now aware of dislocation type - corrected change rate for "dipole size" dupper - corrected term for dipole dissociation by stress change - added transmissivity term in fluxes which accounts for misorientation between two neighboring grains (yet hardcoded transmissivity according to misorientation angle) - added more output variables constitutive: - 2 additional variables "previousDotState" and "previousDotState2", which are used to store the previous and second previous dotState (used in crystallite for acceleration/stabilization of state integration) - timer for dotState now measures the time for calls to constitutive_ collectState (used to reside in crystallite_updateState, which is not critical in terms of calculation time anymore) crystallite: - convergence check for nonlocal elments is now done at end of crystallite loop, not at the beginning; we simple set all elements to not converged if there is at least one nonlocal element that did not converge - need call to microstructure before first call to collect dotState for dependent states - stiffness calculation (jacobian): if there are nonlocal elements, we also have to consider changes in our neighborhood's states; so for every perturbed component in a single ip, we have to loop over all elements; since this is extremely time-consuming, we just perturb one component per cycle, starting with the one that changes the most during regular time step. - updateState gets a damping prefactor for our dotState that helps to improve convergence; prefactor is calculated according to change of dotState IO: - additional warning message for unknown crystal symmetry
2009-12-15 13:50:31 +05:30
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONE_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_J2_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
plasticState(mappingConstitutive(2,ipc,ip,el))%dotState(:,mappingConstitutive(1,ipc,ip,el)) &
= constitutive_j2_dotState(Tstar_v,plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)), ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_dotState(ipc,ip,el)%p = constitutive_j2_dotState(Tstar_v,&
constitutive_state(ipc,ip,el), ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_PHENOPOWERLAW_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
plasticState(mappingConstitutive(2,ipc,ip,el))%dotState(:,mappingConstitutive(1,ipc,ip,el)) &
= constitutive_phenopowerlaw_dotState(Tstar_v,plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)), ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_dotState(ipc,ip,el)%p = constitutive_phenopowerlaw_dotState(Tstar_v,&
constitutive_state(ipc,ip,el), ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
case (PLASTICITY_DISLOTWIN_ID)
fixed array out of bounds error due to misplaced #ifdef
2014-06-05 20:52:35 +05:30
#ifdef NEWSTATE
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
plasticState(mappingConstitutive(2,ipc,ip,el))%dotState(:,mappingConstitutive(1,ipc,ip,el)) &
= constitutive_dislotwin_dotState(Tstar_v,Temperature,&
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)), ipc,ip,el)
#else
constitutive_dotState(ipc,ip,el)%p = constitutive_dislotwin_dotState(Tstar_v,Temperature,&
constitutive_state(ipc,ip,el), ipc,ip,el)
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifdef NEWSTATE
plasticState(mappingConstitutive(2,ipc,ip,el))%dotState(:,mappingConstitutive(1,ipc,ip,el)) &
= constitutive_titanmod_dotState(Tstar_v,Temperature,&
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)), ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_dotState(ipc,ip,el)%p = constitutive_titanmod_dotState(Tstar_v,Temperature,&
constitutive_state(ipc,ip,el), ipc,ip,el)
new state related changes
2014-06-11 17:57:41 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
!*
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
plasticState(mappingConstitutive(2,ipc,ip,el))%dotState(:,mappingConstitutive(1,ipc,ip,el)) = &
constitutive_nonlocal_dotState(Tstar_v, FeArray, FpArray, &
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
Temperature, subdt, &
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
subfracArray, ipc, ip, el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_dotState(ipc,ip,el)%p = constitutive_nonlocal_dotState(Tstar_v, FeArray, FpArray, &
Temperature, constitutive_state, constitutive_state0, subdt, &
subfracArray, ipc, ip, el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end select
if (iand(debug_level(debug_constitutive), debug_levelBasic) /= 0_pInt) then
call system_clock(count=tock,count_rate=tickrate,count_max=maxticks)
!$OMP CRITICAL (debugTimingDotState)
debug_cumDotStateCalls = debug_cumDotStateCalls + 1_pInt
debug_cumDotStateTicks = debug_cumDotStateTicks + tock-tick
!$OMP FLUSH (debug_cumDotStateTicks)
if (tock < tick) debug_cumDotStateTicks = debug_cumDotStateTicks + maxticks
!$OMP END CRITICAL (debugTimingDotState)
endif
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end subroutine constitutive_collectDotState
changed line endings to unix style
2009-03-06 15:32:36 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
!> @brief for constitutive models having an instantaneous change of state (so far, only nonlocal)
!> will return false if delta state is not needed/supported by the constitutive model
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
logical function constitutive_collectDeltaState(Tstar_v, ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use prec, only: &
pLongInt
use debug, only: &
debug_cumDeltaStateCalls, &
debug_cumDeltaStateTicks, &
debug_level, &
debug_constitutive, &
debug_levelBasic
use material, only: &
phase_plasticity, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
material_phase, &
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#endif
PLASTICITY_NONLOCAL_ID
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_nonlocal, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_nonlocal_deltaState
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(6) :: &
Tstar_v !< 2nd Piola-Kirchhoff stress
integer(pLongInt) :: &
tick, tock, &
tickrate, &
maxticks
if (iand(debug_level(debug_constitutive), debug_levelBasic) /= 0_pInt) &
call system_clock(count=tick,count_rate=tickrate,count_max=maxticks)
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState).
2012-05-16 20:13:26 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
constitutive_collectDeltaState = .true.
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
call constitutive_nonlocal_deltaState(Tstar_v,ip,el)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#else
fixed floating point exception in nonlocal (calculation of dUpper at several places)
2014-06-27 01:31:38 +05:30
call constitutive_nonlocal_deltaState(constitutive_deltaState(ipc,ip,el)%p,&
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_state(ipc,ip,el), Tstar_v,ipc,ip,el)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
case default
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
constitutive_collectDeltaState = .false.
#ifndef NEWSTATE
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_deltaState(ipc,ip,el)%p = 0.0_pReal !ToDo: needed or will it remain zero anyway?
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end select
all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState).
2012-05-16 20:13:26 +05:30
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
if (iand(debug_level(debug_constitutive), debug_levelBasic) /= 0_pInt) then
call system_clock(count=tock,count_rate=tickrate,count_max=maxticks)
!$OMP CRITICAL (debugTimingDeltaState)
debug_cumDeltaStateCalls = debug_cumDeltaStateCalls + 1_pInt
debug_cumDeltaStateTicks = debug_cumDeltaStateTicks + tock-tick
!$OMP FLUSH (debug_cumDeltaStateTicks)
if (tock < tick) debug_cumDeltaStateTicks = debug_cumDeltaStateTicks + maxticks
!$OMP END CRITICAL (debugTimingDeltaState)
endif
all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState).
2012-05-16 20:13:26 +05:30
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 12:24:49 +05:30
end function constitutive_collectDeltaState
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns array of constitutive results
!--------------------------------------------------------------------------------------------------
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
function constitutive_postResults(Tstar_v, FeArray, temperature, ipc, ip, el)
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
mappingConstitutive, &
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
phase_plasticity, &
material_phase, &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
homogenization_maxNgrains, &
PLASTICITY_NONE_ID, &
PLASTICITY_J2_ID, &
PLASTICITY_PHENOPOWERLAW_ID, &
PLASTICITY_DISLOTWIN_ID, &
PLASTICITY_TITANMOD_ID, &
PLASTICITY_NONLOCAL_ID
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_j2, only: &
modified reading in of values (realizes, if phase is not used in the model at all), started to include HDF5 output and new state
2014-04-15 15:39:20 +05:30
#ifdef HDF
constitutive_j2_postResults2,&
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_j2_postResults
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_phenopowerlaw, only: &
more changes related to new state
2014-05-22 20:54:12 +05:30
constitutive_phenopowerlaw_postResults
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_dislotwin, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_dislotwin_postResults
use constitutive_titanmod, only: &
constitutive_titanmod_postResults
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
use constitutive_nonlocal, only: &
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_nonlocal_postResults
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
implicit none
integer(pInt), intent(in) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
ipc, & !< grain number
ip, & !< integration point number
el !< element number
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
real(pReal), dimension(constitutive_sizePostResults(ipc,ip,el)) :: &
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
constitutive_postResults
real(pReal), intent(in) :: &
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
temperature
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(3,3,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: &
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
FeArray !< elastic deformation gradient
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
real(pReal), intent(in), dimension(6) :: &
Tstar_v !< 2nd Piola Kirchhoff stress tensor (Mandel)
constitutive_postResults = 0.0_pReal
changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw
2013-09-19 13:16:01 +05:30
select case (phase_plasticity(material_phase(ipc,ip,el)))
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONE_ID)
case (PLASTICITY_TITANMOD_ID)
new state related changes
2014-06-11 17:57:41 +05:30
#ifdef NEWSTATE
constitutive_postResults = constitutive_titanmod_postResults(&
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_postResults = constitutive_titanmod_postResults(&
constitutive_state(ipc,ip,el),ipc,ip,el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_J2_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
constitutive_postResults= constitutive_j2_postResults(Tstar_v, &
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_postResults = constitutive_j2_postResults(Tstar_v,&
more changes related to new state
2014-05-22 20:54:12 +05:30
constitutive_state(ipc,ip,el),ipc,ip,el)
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_PHENOPOWERLAW_ID)
more changes related to new state
2014-05-22 20:54:12 +05:30
#ifdef NEWSTATE
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_postResults = constitutive_phenopowerlaw_postResults(Tstar_v,&
more changes related to new state
2014-05-22 20:54:12 +05:30
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
constitutive_postResults = constitutive_phenopowerlaw_postResults(Tstar_v,&
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
constitutive_state(ipc,ip,el),ipc,ip,el)
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_DISLOTWIN_ID)
prepared new state for dislotwin
2014-06-03 19:16:42 +05:30
#ifdef NEWSTATE
constitutive_postResults = constitutive_dislotwin_postResults(Tstar_v,Temperature,&
plasticState(mappingConstitutive(2,ipc,ip,el))% &
state(:,mappingConstitutive(1,ipc,ip,el)),ipc,ip,el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_postResults = constitutive_dislotwin_postResults(Tstar_v,Temperature,&
constitutive_state(ipc,ip,el),ipc,ip,el)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#endif
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
case (PLASTICITY_NONLOCAL_ID)
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-14 02:23:17 +05:30
#ifdef NEWSTATE
constitutive_postResults = constitutive_nonlocal_postResults(Tstar_v, FeArray, &
mappingConstitutive, ipc, ip, el)
#else
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
constitutive_postResults = constitutive_nonlocal_postResults(Tstar_v, FeArray, &
constitutive_state, constitutive_dotstate(ipc,ip,el), ipc, ip, el)
more changes related to new state
2014-05-22 20:54:12 +05:30
#endif
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 15:44:57 +05:30
end select
openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules: DON'T use implicit array subscripts: example: real, dimension(3,3) :: A,B A(:,2) = B(:,1) <--- DON'T USE A(1:3,2) = B(1:3,1) <--- BETTER USE In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks. Enclose all write statements with the following: !$OMP CRITICAL (write2out) <your write statement> !$OMP END CRITICAL (write2out) Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 16:16:17 +05:30
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
end function constitutive_postResults
changed line endings to unix style
2009-03-06 15:32:36 +05:30
introduced keyword elasticity in material.config added constitutive_TandItsTangent and special version constitutive_hooke_TandItsTangent those are however not yet used and need to be tested
2012-03-14 21:46:11 +05:30
removed unused variables, added pure statements, and declared external functions as external (using gfortrans debug option)
2013-02-11 16:13:45 +05:30
end module constitutive