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

1453 lines
77 KiB
Fortran
Raw Normal View History

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
!--------------------------------------------------------------------------------------------------
removed copyright from *.f90 and *.f files, will be added using the deployMe script for release/nightly build. Now based on presence of svn:keywords property set (usually "Id"). Similar thing will be done for .py, .sh, and .config files but more involved because they might be executable
2014-03-29 13:50:36 +05:30
! $Id$
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
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables
2013-06-11 22:05:04 +05:30
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!> @brief Parses material config file, either solverJobName.materialConfig or material.config
!> @details reads the material configuration file, where solverJobName.materialConfig takes
!! precedence over material.config and parses the sections 'homogenization', 'crystallite',
!! 'phase', 'texture', and 'microstucture'
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
!--------------------------------------------------------------------------------------------------
module material
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use prec, only: &
pReal, &
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
pInt, &
introduced new state structure, compile with STATE=NEW
2014-04-15 14:50:38 +05:30
tState, &
polished field state related changes and brittle damage
2014-09-22 23:45:19 +05:30
tFieldData, &
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
p_intvec
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
implicit none
private
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=*), parameter, public :: &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
ELASTICITY_hooke_label = 'hooke', &
PLASTICITY_none_label = 'none', &
PLASTICITY_j2_label = 'j2', &
PLASTICITY_phenopowerlaw_label = 'phenopowerlaw', &
PLASTICITY_dislotwin_label = 'dislotwin', &
PLASTICITY_dislokmc_label = 'dislokmc', &
PLASTICITY_titanmod_label = 'titanmod', &
PLASTICITY_nonlocal_label = 'nonlocal', &
LOCAL_DAMAGE_none_label = 'none', &
LOCAL_DAMAGE_brittle_label = 'brittle', &
LOCAL_DAMAGE_ductile_label = 'ductile', &
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
LOCAL_DAMAGE_gurson_label = 'gurson', &
started implementing anisotropic damage
2014-10-15 20:32:29 +05:30
LOCAL_DAMAGE_anisotropic_label = 'anisotropic', &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
LOCAL_THERMAL_isothermal_label = 'isothermal', &
LOCAL_THERMAL_adiabatic_label = 'adiabatic', &
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
LOCAL_VACANCY_constant_label = 'constant', &
LOCAL_VACANCY_generation_label = 'generation', &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
FIELD_DAMAGE_local_label = 'local', &
FIELD_DAMAGE_nonlocal_label = 'nonlocal', &
FIELD_THERMAL_local_label = 'local', &
FIELD_THERMAL_nonlocal_label = 'nonlocal', &
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
FIELD_VACANCY_local_label = 'local', &
FIELD_VACANCY_nonlocal_label = 'nonlocal', &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
HOMOGENIZATION_none_label = 'none', &
HOMOGENIZATION_isostrain_label = 'isostrain', &
HOMOGENIZATION_rgc_label = 'rgc'
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
file was missing in last commit
2013-11-27 13:35:23 +05:30
enum, bind(c)
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
enumerator :: ELASTICITY_undefined_ID, &
ELASTICITY_hooke_ID
end enum
enum, bind(c)
enumerator :: PLASTICITY_undefined_ID, &
PLASTICITY_none_ID, &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
PLASTICITY_j2_ID, &
file was missing in last commit
2013-11-27 13:35:23 +05:30
PLASTICITY_phenopowerlaw_ID, &
PLASTICITY_dislotwin_ID, &
added dislokmc (dislotwin modification)
2014-08-08 16:34:40 +05:30
PLASTICITY_dislokmc_ID, &
file was missing in last commit
2013-11-27 13:35:23 +05:30
PLASTICITY_titanmod_ID, &
PLASTICITY_nonlocal_ID
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
end enum
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
enum, bind(c)
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
enumerator :: LOCAL_DAMAGE_none_ID, &
LOCAL_DAMAGE_brittle_ID, &
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
LOCAL_DAMAGE_ductile_ID, &
started implementing anisotropic damage
2014-10-15 20:32:29 +05:30
LOCAL_DAMAGE_gurson_ID, &
LOCAL_DAMAGE_anisotropic_ID
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
end enum
enum, bind(c)
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
enumerator :: LOCAL_THERMAL_isothermal_ID, &
LOCAL_THERMAL_adiabatic_ID
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
end enum
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
enum, bind(c)
enumerator :: LOCAL_VACANCY_constant_ID, &
LOCAL_VACANCY_generation_ID
end enum
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
enum, bind(c)
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
enumerator :: FIELD_DAMAGE_local_ID ,&
FIELD_DAMAGE_nonlocal_ID
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
end enum
enum, bind(c)
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
enumerator :: FIELD_THERMAL_local_ID, &
FIELD_THERMAL_nonlocal_ID
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
end enum
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
enum, bind(c)
enumerator :: FIELD_VACANCY_local_ID, &
FIELD_VACANCY_nonlocal_ID
end enum
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
enum, bind(c)
enumerator :: HOMOGENIZATION_undefined_ID, &
introduced homogenization_none to substitute isostrain with ngrains 1. cleaned up homogenization isostrain (has no state)
2014-03-14 04:50:50 +05:30
HOMOGENIZATION_none_ID, &
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
HOMOGENIZATION_isostrain_ID, &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
HOMOGENIZATION_rgc_ID
file was missing in last commit
2013-11-27 13:35:23 +05:30
end enum
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=*), parameter, public :: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
MATERIAL_configFile = 'material.config', & !< generic name for material configuration file
MATERIAL_localFileExt = 'materialConfig' !< extension of solver job name depending material configuration file
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=*), parameter, public :: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
MATERIAL_partHomogenization = 'homogenization', & !< keyword for homogenization part
MATERIAL_partCrystallite = 'crystallite', & !< keyword for crystallite part
MATERIAL_partPhase = 'phase' !< keyword for phase part
file was missing in last commit
2013-11-27 13:35:23 +05:30
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
integer(kind(ELASTICITY_undefined_ID)), dimension(:), allocatable, public, protected :: &
some small improvements, still something wrong, probably in crystallite.f90
2013-12-13 19:44:17 +05:30
phase_elasticity !< elasticity of each phase
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
integer(kind(PLASTICITY_undefined_ID)), dimension(:), allocatable, public, protected :: &
some small improvements, still something wrong, probably in crystallite.f90
2013-12-13 19:44:17 +05:30
phase_plasticity !< plasticity of each phase
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
integer(kind(LOCAL_DAMAGE_none_ID)), dimension(:), allocatable, public, protected :: &
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
phase_damage !< local damage of each phase
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
integer(kind(LOCAL_THERMAL_isothermal_ID)), dimension(:), allocatable, public, protected :: &
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
phase_thermal !< local thermal of each phase
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
integer(kind(LOCAL_VACANCY_constant_ID)), dimension(:), allocatable, public, protected :: &
phase_vacancy !< local vacancy model of each phase
integer(kind(FIELD_DAMAGE_local_ID)), dimension(:), allocatable, public, protected :: &
field_damage_type !< field damage of each phase
integer(kind(FIELD_THERMAL_local_ID)), dimension(:), allocatable, public, protected :: &
field_thermal_type !< field thermal of each phase
integer(kind(FIELD_VACANCY_local_ID)), dimension(:), allocatable, public, protected :: &
field_vacancy_type !< field vacancy of each phase
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
some small improvements, still something wrong, probably in crystallite.f90
2013-12-13 19:44:17 +05:30
integer(kind(HOMOGENIZATION_undefined_ID)), dimension(:), allocatable, public, protected :: &
file was missing in last commit
2013-11-27 13:35:23 +05:30
homogenization_type !< type of each homogenization
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=64), dimension(:), allocatable, public, protected :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
phase_name, & !< name of each phase
homogenization_name, & !< name of each homogenization
crystallite_name !< name of each crystallite setting
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
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
integer(pInt), public, protected :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
homogenization_maxNgrains, & !< max number of grains in any USED homogenization
material_Nphase, & !< number of phases
material_Nhomogenization, & !< number of homogenizations
material_Nmicrostructure, & !< number of microstructures
material_Ncrystallite !< number of crystallite settings
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
integer(pInt), dimension(:), allocatable, public, protected :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
homogenization_Ngrains, & !< number of grains in each homogenization
homogenization_Noutput, & !< number of '(output)' items per homogenization
phase_Noutput, & !< number of '(output)' items per phase
phase_elasticityInstance, & !< instance of particular elasticity of each phase
phase_plasticityInstance, & !< instance of particular plasticity of each phase
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_damageInstance, & !< instance of particular damage of each phase
phase_thermalInstance, & !< instance of particular thermal of each phase
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
phase_vacancyInstance, & !< instance of particular vacancy model of each phase
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
crystallite_Noutput, & !< number of '(output)' items per crystallite setting
homogenization_typeInstance, & !< instance of particular type of each homogenization
microstructure_crystallite !< crystallite setting ID of each microstructure
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
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 13:03:51 +05:30
integer(pInt), dimension(:,:,:), allocatable, public :: &
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
material_phase !< phase (index) of each grain,IP,element
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
integer(pInt), dimension(:,:), allocatable, public :: &
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
material_homog !< homogenization (index) of each IP,element
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
type(tState), allocatable, dimension(:), public :: &
introduced new state structure, compile with STATE=NEW
2014-04-15 14:50:38 +05:30
plasticState, &
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
damageState, &
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
thermalState,&
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
vacancyState,&
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
homogState
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
polished field state related changes and brittle damage
2014-09-22 23:45:19 +05:30
type(tFieldData), allocatable, dimension(:), public :: &
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
fieldDamage
polished field state related changes and brittle damage
2014-09-22 23:45:19 +05:30
type(tFieldData), allocatable, dimension(:), public :: &
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
fieldThermal
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
type(tFieldData), allocatable, dimension(:), public :: &
fieldVacancy
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
integer(pInt), dimension(:,:,:), allocatable, public, protected :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
material_texture !< texture (index) of each grain,IP,element
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
real(pReal), dimension(:,:,:,:), allocatable, public, protected :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
material_EulerAngles !< initial orientation of each grain,IP,element
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
logical, dimension(:), allocatable, public, protected :: &
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
microstructure_active, &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
microstructure_elemhomo, & !< flag to indicate homogeneous microstructure distribution over element's IPs
phase_localPlasticity !< flags phases with local constitutive law
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=*), parameter, private :: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
MATERIAL_partMicrostructure = 'microstructure', & !< keyword for microstructure part
MATERIAL_partTexture = 'texture' !< keyword for texture part
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=64), dimension(:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
microstructure_name, & !< name of each microstructure
texture_name !< name of each texture
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
character(len=256), dimension(:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
texture_ODFfile !< name of each ODF file
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
integer(pInt), private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
material_Ntexture, & !< number of textures
microstructure_maxNconstituents, & !< max number of constituents in any phase
texture_maxNgauss, & !< max number of Gauss components in any texture
texture_maxNfiber !< max number of Fiber components in any texture
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
integer(pInt), dimension(:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
microstructure_Nconstituents, & !< number of constituents in each microstructure
texture_symmetry, & !< number of symmetric orientations per texture
texture_Ngauss, & !< number of Gauss components per texture
texture_Nfiber !< number of Fiber components per texture
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
integer(pInt), dimension(:,:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
microstructure_phase, & !< phase IDs of each microstructure
microstructure_texture !< texture IDs of each microstructure
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
real(pReal), dimension(:,:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
microstructure_fraction !< vol fraction of each constituent in microstructure
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
real(pReal), dimension(:,:,:), allocatable, private :: &
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
material_volume, & !< volume of each grain,IP,element
texture_Gauss, & !< data of each Gauss component
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
texture_Fiber, & !< data of each Fiber component
renamed texture_rotation into texture_transformation changed keyword for material.config to axes in line with geom_fromAng.py former keyword rotation is still recognized for compatibility
2013-07-24 16:39:39 +05:30
texture_transformation !< transformation for each texture
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
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
logical, dimension(:), allocatable, private :: &
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
homogenization_active
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
integer(pInt), dimension(:,:,:,:), allocatable, public, protected :: mappingConstitutive
integer(pInt), dimension(:,:,:), allocatable, public, protected :: mappingCrystallite
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
integer(pInt), dimension(:,:,:), allocatable, public, protected :: mappingHomogenization
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +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
public :: &
file was missing in last commit
2013-11-27 13:35:23 +05:30
material_init, &
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
ELASTICITY_hooke_ID ,&
file was missing in last commit
2013-11-27 13:35:23 +05:30
PLASTICITY_none_ID, &
PLASTICITY_J2_ID, &
PLASTICITY_phenopowerlaw_ID, &
PLASTICITY_dislotwin_ID, &
added dislokmc (dislotwin modification)
2014-08-08 16:34:40 +05:30
PLASTICITY_dislokmc_ID, &
file was missing in last commit
2013-11-27 13:35:23 +05:30
PLASTICITY_titanmod_ID, &
PLASTICITY_nonlocal_ID, &
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
LOCAL_DAMAGE_none_ID, &
LOCAL_DAMAGE_brittle_ID, &
started adding a ductile damage model
2014-09-26 16:08:13 +05:30
LOCAL_DAMAGE_ductile_ID, &
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
LOCAL_DAMAGE_gurson_ID, &
started implementing anisotropic damage
2014-10-15 20:32:29 +05:30
LOCAL_DAMAGE_anisotropic_ID, &
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
LOCAL_THERMAL_isothermal_ID, &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
LOCAL_THERMAL_adiabatic_ID, &
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
LOCAL_VACANCY_constant_ID, &
LOCAL_VACANCY_generation_ID, &
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
FIELD_DAMAGE_local_ID, &
FIELD_DAMAGE_nonlocal_ID, &
FIELD_THERMAL_local_ID, &
FIELD_THERMAL_nonlocal_ID, &
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
FIELD_VACANCY_local_ID, &
FIELD_VACANCY_nonlocal_ID, &
introduced homogenization_none to substitute isostrain with ngrains 1. cleaned up homogenization isostrain (has no state)
2014-03-14 04:50:50 +05:30
HOMOGENIZATION_none_ID, &
file was missing in last commit
2013-11-27 13:35:23 +05:30
HOMOGENIZATION_isostrain_ID, &
small changes to HDF5, put link to homepage into README not to have double effort
2014-03-25 22:51:47 +05:30
#ifdef HDF
material_NconstituentsPhase, &
#endif
file was missing in last commit
2013-11-27 13:35:23 +05:30
HOMOGENIZATION_RGC_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
private :: &
material_parseHomogenization, &
material_parseMicrostructure, &
material_parseCrystallite, &
material_parsePhase, &
material_parseTexture, &
material_populateGrains
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses material configuration file
!> @details figures out if solverJobName.materialConfig is present, if not looks for
!> material.config
!--------------------------------------------------------------------------------------------------
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
subroutine material_init(temperature_init)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +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 documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use IO, only: &
IO_error, &
IO_open_file, &
reporting current time for most module inits
2013-02-25 22:04:59 +05:30
IO_open_jobFile_stat, &
IO_timeStamp
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use debug, only: &
debug_level, &
debug_material, &
debug_levelBasic, &
debug_levelExtensive
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
use mesh, only: &
mesh_maxNips, &
mesh_NcpElems, &
mesh_element, &
FE_Nips, &
FE_geomtype
only output from the root processor for parallel runs
2014-10-10 01:53:06 +05:30
use numerics, only: &
worldrank
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
real(pReal), intent(in) :: temperature_init !< initial field temperature
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
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
integer(pInt) :: m,c,h, myDebug, myHomog
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
integer(pInt) :: &
g, & !< grain number
i, & !< integration point number
e, & !< element number
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
phase, &
homog, &
NofMyField
integer(pInt), dimension(:), allocatable :: ConstitutivePosition
integer(pInt), dimension(:), allocatable :: CrystallitePosition
integer(pInt), dimension(:), allocatable :: HomogenizationPosition
changed variable name "debug_what" to "debug_level"
2012-07-05 15:24:50 +05:30
myDebug = debug_level(debug_material)
numerics.config is now optional; if this file does not exist in the working directory, standard values from numerics.f90 are used; either way the output file tells which values were taken. added some output to the init subroutines of debug.f90, FEsolving.f90, lattice.f90, material.f90 and mesh.f90.
2009-06-18 19:58:02 +05:30
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
mainProcess: if (worldrank == 0) then
write(6,'(/,a)') ' <<<+- material init -+>>>'
write(6,'(a)') ' $Id$'
write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
made code standard conform to Fortran 2008 (ignoring warning concerning comments beyond character 132). Basically, changing "x" format specifier to "Nx" ("1x") plus removing $ format specifier added compiler switches for gfortran and ifort to check for standard conformity old gnu compilers <4.4 are not longer supported because they don't provide the c binding for fftw
2012-02-01 00:48:55 +05:30
#include "compilation_info.f90"
updated MPI reporting in line with recent changes
2014-10-10 21:51:10 +05:30
endif mainProcess
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
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
call material_parseHomogenization(FILEUNIT,material_partHomogenization)
file was missing in last commit
2013-11-27 13:35:23 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) write(6,'(a)') ' Homogenization parsed'
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 material_parseMicrostructure(FILEUNIT,material_partMicrostructure)
file was missing in last commit
2013-11-27 13:35:23 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) write(6,'(a)') ' Microstructure parsed'
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 material_parseCrystallite(FILEUNIT,material_partCrystallite)
file was missing in last commit
2013-11-27 13:35:23 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) write(6,'(a)') ' Crystallite parsed'
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 material_parseTexture(FILEUNIT,material_partTexture)
file was missing in last commit
2013-11-27 13:35:23 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) write(6,'(a)') ' Texture parsed'
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 material_parsePhase(FILEUNIT,material_partPhase)
file was missing in last commit
2013-11-27 13:35:23 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) write(6,'(a)') ' Phase parsed'
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)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
introduced new state structure, compile with STATE=NEW
2014-04-15 14:50:38 +05:30
allocate(plasticState(material_Nphase))
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
allocate(damageState (material_Nphase))
allocate(thermalState(material_Nphase))
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
allocate(vacancyState(material_Nphase))
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
allocate(homogState (material_Nhomogenization))
allocate(fieldDamage (material_Nhomogenization))
allocate(fieldThermal(material_Nhomogenization))
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
allocate(fieldVacancy(material_Nhomogenization))
general polishing
2013-03-28 19:20:20 +05:30
do m = 1_pInt,material_Nmicrostructure
file was missing in last commit
2013-11-27 13:35:23 +05:30
if(microstructure_crystallite(m) < 1_pInt .or. &
microstructure_crystallite(m) > material_Ncrystallite) &
call IO_error(150_pInt,m,ext_msg='crystallite')
if(minval(microstructure_phase(1:microstructure_Nconstituents(m),m)) < 1_pInt .or. &
maxval(microstructure_phase(1:microstructure_Nconstituents(m),m)) > material_Nphase) &
call IO_error(150_pInt,m,ext_msg='phase')
if(minval(microstructure_texture(1:microstructure_Nconstituents(m),m)) < 1_pInt .or. &
maxval(microstructure_texture(1:microstructure_Nconstituents(m),m)) > material_Ntexture) &
call IO_error(150_pInt,m,ext_msg='texture')
if(microstructure_Nconstituents(m) < 1_pInt) &
call IO_error(151_pInt,m)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
general polishing
2013-03-28 19:20:20 +05:30
debugOut: if (iand(myDebug,debug_levelExtensive) /= 0_pInt) then
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
write(6,'(/,a,/)') ' MATERIAL configuration'
write(6,'(a32,1x,a16,1x,a6)') 'homogenization ','type ','grains'
general polishing
2013-03-28 19:20:20 +05:30
do h = 1_pInt,material_Nhomogenization
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
write(6,'(1x,a32,1x,a16,1x,i6)') homogenization_name(h),homogenization_type(h),homogenization_Ngrains(h)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
enddo
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
write(6,'(/,a14,18x,1x,a11,1x,a12,1x,a13)') 'microstructure','crystallite','constituents','homogeneous'
general polishing
2013-03-28 19:20:20 +05:30
do m = 1_pInt,material_Nmicrostructure
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
write(6,'(1x,a32,1x,i11,1x,i12,1x,l13)') microstructure_name(m), &
general polishing
2013-03-28 19:20:20 +05:30
microstructure_crystallite(m), &
microstructure_Nconstituents(m), &
microstructure_elemhomo(m)
if (microstructure_Nconstituents(m) > 0_pInt) then
do c = 1_pInt,microstructure_Nconstituents(m)
write(6,'(a1,1x,a32,1x,a32,1x,f7.4)') '>',phase_name(microstructure_phase(c,m)),&
texture_name(microstructure_texture(c,m)),&
microstructure_fraction(c,m)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
enddo
write(6,*)
endif
enddo
general polishing
2013-03-28 19:20:20 +05:30
endif debugOut
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 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
call material_populateGrains
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
allocate(mappingConstitutive(2,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems),source=0_pInt)
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
allocate(mappingHomogenization(2,mesh_maxNips,mesh_NcpElems),source=0_pInt)
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
allocate(mappingCrystallite (2,homogenization_maxNgrains,mesh_NcpElems),source=0_pInt)
allocate(ConstitutivePosition(material_Nphase),source=0_pInt)
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
allocate(HomogenizationPosition(material_Nhomogenization),source=0_pInt)
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
allocate(CrystallitePosition(material_Nphase),source=0_pInt)
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
ElemLoop:do e = 1_pInt,mesh_NcpElems ! loop over elements
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
myHomog = mesh_element(3,e)
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
IPloop:do i = 1_pInt,FE_Nips(FE_geomtype(mesh_element(2,e))) ! loop over IPs
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
HomogenizationPosition(myHomog) = HomogenizationPosition(myHomog) + 1_pInt
mappingHomogenization(1:2,i,e) = [HomogenizationPosition(myHomog),myHomog]
added damage and thermal state vectors and phase instances. only works with new state layout.
2014-06-24 19:21:17 +05:30
GrainLoop:do g = 1_pInt,homogenization_Ngrains(mesh_element(3,e)) ! loop over grains
phase = material_phase(g,i,e)
ConstitutivePosition(phase) = ConstitutivePosition(phase)+1_pInt ! not distinguishing between instances of same phase
mappingConstitutive(1:2,g,i,e) = [ConstitutivePosition(phase),phase]
enddo GrainLoop
enddo IPloop
enddo ElemLoop
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
do homog = 1,material_Nhomogenization
NofMyField=count(material_homog==homog)
select case(field_damage_type(homog))
case (FIELD_DAMAGE_LOCAL_ID)
fieldDamage(homog)%sizeField = 0_pInt
fieldDamage(homog)%sizePostResults = 0_pInt
allocate(fieldDamage(homog)%field(fieldDamage(homog)%sizeField,NofMyField), source = 1.0_pReal)
case (FIELD_DAMAGE_NONLOCAL_ID)
fieldDamage(homog)%sizeField = 1_pInt
fieldDamage(homog)%sizePostResults = 1_pInt
allocate(fieldDamage(homog)%field(fieldDamage(homog)%sizeField,NofMyField), source = 1.0_pReal)
end select
enddo
do homog = 1,material_Nhomogenization
NofMyField=count(material_homog==homog)
select case(field_thermal_type(homog))
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
case (FIELD_THERMAL_local_ID)
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
fieldThermal(homog)%sizeField = 0_pInt
fieldThermal(homog)%sizePostResults = 0_pInt
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
allocate(fieldThermal(homog)%field(fieldThermal(homog)%sizeField,NofMyField))
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
case (FIELD_THERMAL_nonlocal_ID)
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
fieldThermal(homog)%sizeField = 1_pInt
fieldThermal(homog)%sizePostResults = 1_pInt
allocate(fieldThermal(homog)%field(fieldThermal(homog)%sizeField,NofMyField), &
temperature initialised to value from load case file
2014-10-13 23:24:27 +05:30
source = temperature_init)
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
end select
enddo
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
do homog = 1,material_Nhomogenization
NofMyField=count(material_homog==homog)
select case(field_vacancy_type(homog))
case (FIELD_VACANCY_local_ID)
fieldVacancy(homog)%sizeField = 0_pInt
fieldVacancy(homog)%sizePostResults = 0_pInt
allocate(fieldVacancy(homog)%field(fieldVacancy(homog)%sizeField,NofMyField), &
source = 0.0_pReal)
case (FIELD_VACANCY_nonlocal_ID)
fieldVacancy(homog)%sizeField = 1_pInt
fieldVacancy(homog)%sizePostResults = 1_pInt
allocate(fieldVacancy(homog)%field(fieldVacancy(homog)%sizeField,NofMyField), &
source = 0.0_pReal)
end select
enddo
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine material_init
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses the homogenization part in the material configuration 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
subroutine material_parseHomogenization(fileUnit,myPart)
general polishing
2013-03-28 19:20:20 +05:30
use IO, only: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
IO_read, &
general polishing
2013-03-28 19:20:20 +05:30
IO_globalTagInPart, &
IO_countSections, &
IO_error, &
IO_countTagInPart, &
IO_lc, &
IO_getTag, &
IO_isBlank, &
IO_stringValue, &
IO_intValue, &
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
IO_stringPos, &
IO_EOF
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use mesh, only: &
mesh_element
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
Update in order to include the c/a-ratio as additional information about the structure (modules lattice, constitutive_ changed accordingly) Correction of the expression for the velocity of mobile dislocation densities in constitutive_dislobased.f90 Correction of some Fortran syntax, incompatible with Compaq Visual Fortran (e.g. integer(8))
2009-03-20 20:04:24 +05:30
character(len=*), intent(in) :: myPart
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), intent(in) :: fileUnit
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
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 :: MAXNCHUNKS = 2_pInt
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
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), dimension(1+2*MAXNCHUNKS) :: positions
file was missing in last commit
2013-11-27 13:35:23 +05:30
integer(pInt) :: Nsections, section, s
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
character(len=65536) :: &
tag, line
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
logical :: echo
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
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
echo = IO_globalTagInPart(fileUnit,myPart,'/echo/')
Nsections = IO_countSections(fileUnit,myPart)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
material_Nhomogenization = Nsections
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
if (Nsections < 1_pInt) call IO_error(160_pInt,ext_msg=myPart)
major restructuring of code. homogenization as well as constitutive are now free to choose. the runtime got somewhat longer (25% on simple tests) compared to a hardcoded isostrain homogenization. this might be a point of further optimization at a later stage... please use homogenization_isostrain.f90 as starting point / example for future developments of homog-schemes. the homogenization scheme now can additionally output certain results. hence, the userdata structure at each integration point now looks like this: - sizeHomogPostResults - block of that size containing homogPostResults then for each grain: - sizeGrainPostResults - block of that size containing crystallitePostResults, which consist of: + phaseID + volFrac + Eulers (3) + any constitutive post results requested
2009-05-07 21:57:36 +05:30
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
allocate(homogenization_name(Nsections)); homogenization_name = ''
introduced homogenization_none to substitute isostrain with ngrains 1. cleaned up homogenization isostrain (has no state)
2014-03-14 04:50:50 +05:30
allocate(homogenization_type(Nsections), source=HOMOGENIZATION_undefined_ID)
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
allocate(FIELD_DAMAGE_type(Nsections), source=FIELD_DAMAGE_local_ID)
allocate(FIELD_THERMAL_type(Nsections), source=FIELD_THERMAL_local_ID)
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
allocate(FIELD_VACANCY_type(Nsections), source=FIELD_VACANCY_local_ID)
introduced homogenization_none to substitute isostrain with ngrains 1. cleaned up homogenization isostrain (has no state)
2014-03-14 04:50:50 +05:30
allocate(homogenization_typeInstance(Nsections), source=0_pInt)
allocate(homogenization_Ngrains(Nsections), source=0_pInt)
allocate(homogenization_Noutput(Nsections), source=0_pInt)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
allocate(homogenization_active(Nsections), source=.false.) !!!!!!!!!!!!!!!
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
forall (s = 1_pInt:Nsections) homogenization_active(s) = any(mesh_element(3,:) == s) ! current homogenization used in model? Homogenization view, maximum operations depend on maximum number of homog schemes
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
homogenization_Noutput = IO_countTagInPart(fileUnit,myPart,'(output)',Nsections)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
rewind(fileUnit)
line = '' ! to have it initialized
section = 0_pInt ! - " -
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= myPart) ! wind forward to <homogenization>
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
line = IO_read(fileUnit)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(/,1x,a)') trim(line) ! echo part header
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
do while (trim(line) /= IO_EOF) ! read through sections of material part
line = IO_read(fileUnit)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
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 (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
exit
endif
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(2x,a)') trim(line) ! echo back read lines
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
section = section + 1_pInt
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
homogenization_name(section) = IO_getTag(line,'[',']')
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
if (section > 0_pInt) 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
positions = IO_stringPos(line,MAXNCHUNKS)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case(tag)
case ('type')
file was missing in last commit
2013-11-27 13:35:23 +05:30
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
introduced homogenization_none to substitute isostrain with ngrains 1. cleaned up homogenization isostrain (has no state)
2014-03-14 04:50:50 +05:30
case(HOMOGENIZATION_NONE_label)
homogenization_type(section) = HOMOGENIZATION_NONE_ID
homogenization_Ngrains(section) = 1_pInt
file was missing in last commit
2013-11-27 13:35:23 +05:30
case(HOMOGENIZATION_ISOSTRAIN_label)
homogenization_type(section) = HOMOGENIZATION_ISOSTRAIN_ID
case(HOMOGENIZATION_RGC_label)
homogenization_type(section) = HOMOGENIZATION_RGC_ID
case default
call IO_error(500_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
homogenization_typeInstance(section) = &
count(homogenization_type==homogenization_type(section)) ! count instances
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
case ('field_damage')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
case(FIELD_DAMAGE_LOCAL_label)
FIELD_DAMAGE_type(section) = FIELD_DAMAGE_LOCAL_ID
case(FIELD_DAMAGE_NONLOCAL_label)
FIELD_DAMAGE_type(section) = FIELD_DAMAGE_NONLOCAL_ID
case default
call IO_error(500_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
case ('field_thermal')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
case(FIELD_THERMAL_local_label)
FIELD_THERMAL_type(section) = FIELD_THERMAL_local_ID
case(FIELD_THERMAL_nonlocal_label)
FIELD_THERMAL_type(section) = FIELD_THERMAL_nonlocal_ID
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
case default
call IO_error(500_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
case ('field_vacancy')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
case(FIELD_VACANCY_local_label)
FIELD_VACANCY_type(section) = FIELD_VACANCY_local_ID
case(FIELD_VACANCY_nonlocal_label)
FIELD_VACANCY_type(section) = FIELD_VACANCY_nonlocal_ID
case default
call IO_error(500_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case ('nconstituents','ngrains')
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
homogenization_Ngrains(section) = IO_intValue(line,positions,2_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
endif
enddo
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
homogenization_maxNgrains = maxval(homogenization_Ngrains,homogenization_active)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
general polishing
2013-03-28 19:20:20 +05:30
end subroutine material_parseHomogenization
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses the microstructure part in the material configuration 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
subroutine material_parseMicrostructure(fileUnit,myPart)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
use IO
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use mesh, only: &
mesh_element, &
mesh_NcpElems
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
Update in order to include the c/a-ratio as additional information about the structure (modules lattice, constitutive_ changed accordingly) Correction of the expression for the velocity of mobile dislocation densities in constitutive_dislobased.f90 Correction of some Fortran syntax, incompatible with Compaq Visual Fortran (e.g. integer(8))
2009-03-20 20:04:24 +05:30
character(len=*), intent(in) :: myPart
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), intent(in) :: fileUnit
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
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 :: MAXNCHUNKS = 7_pInt
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
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), dimension(1_pInt+2_pInt*MAXNCHUNKS) :: positions
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
integer(pInt) :: Nsections, section, constituent, e, i
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
character(len=65536) :: &
tag, line
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
logical :: echo
parts of the material.config can now be set to echo their content back into the log-file (i.e. STDOUT / unit 6) use keyword "/echo/" before first section in the respective parts that you want to see echoed back into the log file.
2012-06-26 15:54:54 +05:30
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
echo = IO_globalTagInPart(fileUnit,myPart,'/echo/')
did a lot of polishing: - removed unnecessary "return" before end of subroutine or function: - changed undetermined array length (:) to (1:3) To prevent problems with some code analysing tools: - "3D oneliner loops" (with ";) only for "do" and "enddo" at the same time - removed line continuation in OMP statements made the makefile more flexible, removed heap-arrays switch
2011-09-13 21:24:06 +05:30
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
Nsections = IO_countSections(fileUnit,myPart)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
material_Nmicrostructure = Nsections
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
if (Nsections < 1_pInt) call IO_error(160_pInt,ext_msg=myPart)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
allocate(microstructure_name(Nsections)); microstructure_name = ''
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(microstructure_crystallite(Nsections), source=0_pInt)
allocate(microstructure_Nconstituents(Nsections), source=0_pInt)
allocate(microstructure_active(Nsections), source=.false.)
allocate(microstructure_elemhomo(Nsections), source=.false.)
quite some changes: # non-greedy memory allocation # generation of outputConstitutive to allow for script-based T16 extraction # exchange of phenomenological by more general phenopowerlaw # lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.) # nicer debugging output # changed some error/warning codes # plus potentially some minor additional brushes here and there
2009-07-22 21:37:19 +05:30
added consistency check in material.f90: Microstructure index in geometry must not exceed number of sections in material.config better error messages in case of recursive file input in IO.f90. also supports absolute path now
2014-05-15 15:10:43 +05:30
if(any(mesh_element(4,1:mesh_NcpElems) > Nsections)) &
call IO_error(155_pInt,ext_msg='Microstructure in geometry > Sections in material.config')
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
forall (e = 1_pInt:mesh_NcpElems) microstructure_active(mesh_element(4,e)) = .true. ! current microstructure used in model? Elementwise view, maximum N operations for N elements
added 'spectralPictureMode' for reading in geometries with 1to1 mapping (e.g. EBSD pattern) !incremental update, wait for commit of damask_spectral.f90 before checking out
2011-10-18 14:48:05 +05:30
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
microstructure_Nconstituents = IO_countTagInPart(fileUnit,myPart,'(constituent)',Nsections)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
microstructure_maxNconstituents = maxval(microstructure_Nconstituents)
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
microstructure_elemhomo = IO_spotTagInPart(fileUnit,myPart,'/elementhomogeneous/',Nsections)
introduction of a possibility to have homogeneous element (all ips in one element are identical, sort of reduced integration) and bugs-fixing in crystallite.f90, homogenization_RGC.f90, numerics.f90.
2009-11-24 20:30:25 +05:30
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(microstructure_phase (microstructure_maxNconstituents,Nsections),source=0_pInt)
allocate(microstructure_texture (microstructure_maxNconstituents,Nsections),source=0_pInt)
allocate(microstructure_fraction(microstructure_maxNconstituents,Nsections),source=0.0_pReal)
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
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
rewind(fileUnit)
general polishing
2013-03-28 19:20:20 +05:30
line = '' ! to have it initialized
section = 0_pInt ! - " -
constituent = 0_pInt ! - " -
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= myPart) ! wind forward to <microstructure>
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
line = IO_read(fileUnit)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(/,1x,a)') trim(line) ! echo part header
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
do while (trim(line) /= IO_EOF) ! read through sections of material part
line = IO_read(fileUnit)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
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 (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
exit
endif
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(2x,a)') trim(line) ! echo back read lines
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
section = section + 1_pInt
constituent = 0_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
microstructure_name(section) = IO_getTag(line,'[',']')
endif
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
if (section > 0_pInt) 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
positions = IO_stringPos(line,MAXNCHUNKS)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case(tag)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
case ('crystallite')
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
microstructure_crystallite(section) = IO_intValue(line,positions,2_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case ('(constituent)')
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
constituent = constituent + 1_pInt
do i=2_pInt,6_pInt,2_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
tag = IO_lc(IO_stringValue(line,positions,i))
select case (tag)
case('phase')
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
microstructure_phase(constituent,section) = IO_intValue(line,positions,i+1_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('texture')
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
microstructure_texture(constituent,section) = IO_intValue(line,positions,i+1_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('fraction')
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
microstructure_fraction(constituent,section) = IO_floatValue(line,positions,i+1_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
enddo
end select
endif
enddo
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
end subroutine material_parseMicrostructure
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses the crystallite part in the material configuration 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
subroutine material_parseCrystallite(fileUnit,myPart)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use IO, only: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
IO_read, &
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
IO_countSections, &
IO_error, &
IO_countTagInPart, &
IO_globalTagInPart, &
IO_getTag, &
IO_lc, &
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
IO_isBlank, &
IO_EOF
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
implicit none
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
character(len=*), intent(in) :: myPart
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), intent(in) :: fileUnit
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
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
integer(pInt) :: Nsections, &
section
character(len=65536) :: line
logical :: echo
parts of the material.config can now be set to echo their content back into the log-file (i.e. STDOUT / unit 6) use keyword "/echo/" before first section in the respective parts that you want to see echoed back into the log file.
2012-06-26 15:54:54 +05:30
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
echo = IO_globalTagInPart(fileUnit,myPart,'/echo/')
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
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
Nsections = IO_countSections(fileUnit,myPart)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
material_Ncrystallite = Nsections
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
if (Nsections < 1_pInt) call IO_error(160_pInt,ext_msg=myPart)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
allocate(crystallite_name(Nsections)); crystallite_name = ''
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(crystallite_Noutput(Nsections), source=0_pInt)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
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
crystallite_Noutput = IO_countTagInPart(fileUnit,myPart,'(output)',Nsections)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
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
rewind(fileUnit)
line = '' ! to have it initialized
section = 0_pInt ! - " -
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= myPart) ! wind forward to <Crystallite>
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
line = IO_read(fileUnit)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
enddo
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(/,1x,a)') trim(line) ! echo part header
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
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
do while (trim(line) /= IO_EOF) ! read through sections of material part
line = IO_read(fileUnit)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
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 (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
exit
endif
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(2x,a)') trim(line) ! echo back read lines
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
section = section + 1_pInt
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
crystallite_name(section) = IO_getTag(line,'[',']')
endif
enddo
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
end subroutine material_parseCrystallite
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses the phase part in the material configuration 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
subroutine material_parsePhase(fileUnit,myPart)
general polishing
2013-03-28 19:20:20 +05:30
use IO, only: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
IO_read, &
general polishing
2013-03-28 19:20:20 +05:30
IO_globalTagInPart, &
IO_countSections, &
IO_error, &
IO_countTagInPart, &
IO_getTag, &
IO_spotTagInPart, &
IO_lc, &
IO_isBlank, &
IO_stringValue, &
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
IO_stringPos, &
IO_EOF
general polishing
2013-03-28 19:20:20 +05:30
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
Update in order to include the c/a-ratio as additional information about the structure (modules lattice, constitutive_ changed accordingly) Correction of the expression for the velocity of mobile dislocation densities in constitutive_dislobased.f90 Correction of some Fortran syntax, incompatible with Compaq Visual Fortran (e.g. integer(8))
2009-03-20 20:04:24 +05:30
character(len=*), intent(in) :: myPart
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), intent(in) :: fileUnit
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
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 :: MAXNCHUNKS = 2_pInt
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
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), dimension(1+2*MAXNCHUNKS) :: positions
fixed bug in initializtion of phase instance, added svn properties for FEM code
2014-10-15 17:52:09 +05:30
integer(pInt) :: Nsections, section, p
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
character(len=65536) :: &
tag,line
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
logical :: echo
parts of the material.config can now be set to echo their content back into the log-file (i.e. STDOUT / unit 6) use keyword "/echo/" before first section in the respective parts that you want to see echoed back into the log file.
2012-06-26 15:54:54 +05:30
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
echo = IO_globalTagInPart(fileUnit,myPart,'/echo/')
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
Nsections = IO_countSections(fileUnit,myPart)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
material_Nphase = Nsections
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
if (Nsections < 1_pInt) call IO_error(160_pInt,ext_msg=myPart)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
general polishing
2013-03-28 19:20:20 +05:30
allocate(phase_name(Nsections)); phase_name = ''
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(phase_elasticity(Nsections), source=ELASTICITY_undefined_ID)
allocate(phase_elasticityInstance(Nsections), source=0_pInt)
allocate(phase_plasticity(Nsections) , source=PLASTICITY_undefined_ID)
allocate(phase_plasticityInstance(Nsections), source=0_pInt)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
allocate(phase_damage(Nsections) , source=LOCAL_DAMAGE_none_ID)
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
allocate(phase_damageInstance(Nsections), source=0_pInt)
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
allocate(phase_thermal(Nsections) , source=LOCAL_THERMAL_isothermal_ID)
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
allocate(phase_thermalInstance(Nsections), source=0_pInt)
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
allocate(phase_vacancy(Nsections) , source=LOCAL_VACANCY_constant_ID)
allocate(phase_vacancyInstance(Nsections), source=0_pInt)
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(phase_Noutput(Nsections), source=0_pInt)
allocate(phase_localPlasticity(Nsections), source=.false.)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
phase_Noutput = IO_countTagInPart(fileUnit,myPart,'(output)',Nsections)
phase_localPlasticity = .not. IO_spotTagInPart(fileUnit,myPart,'/nonlocal/',Nsections)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
rewind(fileUnit)
line = '' ! to have it initialized
section = 0_pInt ! - " -
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= myPart) ! wind forward to <Phase>
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
line = IO_read(fileUnit)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(/,1x,a)') trim(line) ! echo part header
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
do while (trim(line) /= IO_EOF) ! read through sections of material part
line = IO_read(fileUnit)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
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 (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
exit
endif
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(2x,a)') trim(line) ! echo back read lines
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
section = section + 1_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
phase_name(section) = IO_getTag(line,'[',']')
endif
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
if (section > 0_pInt) 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
positions = IO_stringPos(line,MAXNCHUNKS)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case(tag)
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
case ('elasticity')
file was missing in last commit
2013-11-27 13:35:23 +05:30
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
case (ELASTICITY_HOOKE_label)
phase_elasticity(section) = ELASTICITY_HOOKE_ID
case default
call IO_error(200_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
replaced keyword constitution with plasticity in preparation of upcomming generalized elasticity
2012-03-12 19:39:37 +05:30
case ('plasticity')
file was missing in last commit
2013-11-27 13:35:23 +05:30
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
case (PLASTICITY_NONE_label)
phase_plasticity(section) = PLASTICITY_NONE_ID
case (PLASTICITY_J2_label)
phase_plasticity(section) = PLASTICITY_J2_ID
case (PLASTICITY_PHENOPOWERLAW_label)
phase_plasticity(section) = PLASTICITY_PHENOPOWERLAW_ID
case (PLASTICITY_DISLOTWIN_label)
phase_plasticity(section) = PLASTICITY_DISLOTWIN_ID
added dislokmc (dislotwin modification)
2014-08-08 16:34:40 +05:30
case (PLASTICITY_DISLOKMC_label)
phase_plasticity(section) = PLASTICITY_DISLOKMC_ID
file was missing in last commit
2013-11-27 13:35:23 +05:30
case (PLASTICITY_TITANMOD_label)
phase_plasticity(section) = PLASTICITY_TITANMOD_ID
case (PLASTICITY_NONLOCAL_label)
phase_plasticity(section) = PLASTICITY_NONLOCAL_ID
case default
call IO_error(201_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
case ('damage')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
case (LOCAL_DAMAGE_none_label)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
phase_damage(section) = LOCAL_DAMAGE_none_ID
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
case (LOCAL_DAMAGE_brittle_label)
more work on homogenisation new state and introduced field state
2014-09-03 22:42:06 +05:30
phase_damage(section) = LOCAL_DAMAGE_BRITTLE_ID
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
case (LOCAL_DAMAGE_ductile_label)
started adding a ductile damage model
2014-09-26 16:08:13 +05:30
phase_damage(section) = LOCAL_DAMAGE_DUCTILE_ID
started implementing a void fraction based ductile damage model.
2014-10-10 18:12:12 +05:30
case (LOCAL_DAMAGE_gurson_label)
phase_damage(section) = LOCAL_DAMAGE_gurson_ID
started implementing anisotropic damage
2014-10-15 20:32:29 +05:30
case (LOCAL_DAMAGE_anisotropic_label)
phase_damage(section) = LOCAL_DAMAGE_anisotropic_ID
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
case default
call IO_error(200_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
case ('thermal')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal
2014-09-26 21:37:26 +05:30
case (LOCAL_THERMAL_ISOTHERMAL_label)
phase_thermal(section) = LOCAL_THERMAL_isothermal_ID
changed thermal heatGen to thermal adiabatic
2014-10-09 19:38:32 +05:30
case (LOCAL_THERMAL_ADIABATIC_label)
phase_thermal(section) = LOCAL_THERMAL_adiabatic_ID
moved mappings for new state layout to material
2014-06-23 00:28:29 +05:30
case default
call IO_error(200_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage
2014-10-11 02:25:09 +05:30
case ('vacancy')
select case (IO_lc(IO_stringValue(line,positions,2_pInt)))
case (LOCAL_VACANCY_CONSTANT_label)
phase_vacancy(section) = LOCAL_VACANCY_constant_ID
case (LOCAL_VACANCY_GENERATION_label)
phase_vacancy(section) = LOCAL_VACANCY_generation_ID
case default
call IO_error(200_pInt,ext_msg=trim(IO_stringValue(line,positions,2_pInt)))
end select
fixed bug in initializtion of phase instance, added svn properties for FEM code
2014-10-15 17:52:09 +05:30
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
endif
enddo
fixed bug in initializtion of phase instance, added svn properties for FEM code
2014-10-15 17:52:09 +05:30
do p=1_pInt, Nsections
phase_elasticityInstance(p) = count(phase_elasticity(1:p) == phase_elasticity(p))
phase_plasticityInstance(p) = count(phase_plasticity(1:p) == phase_plasticity(p))
phase_damageInstance(p) = count(phase_damage(1:p) == phase_damage(p))
phase_thermalInstance(p) = count(phase_thermal(1:p) == phase_thermal(p))
phase_vacancyInstance(p) = count(phase_vacancy(1:p) == phase_vacancy(p))
enddo
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
fixed bug in initializtion of phase instance, added svn properties for FEM code
2014-10-15 17:52:09 +05:30
end subroutine material_parsePhase
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief parses the texture part in the material configuration 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
subroutine material_parseTexture(fileUnit,myPart)
general polishing
2013-03-28 19:20:20 +05:30
use IO, only: &
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
IO_read, &
general polishing
2013-03-28 19:20:20 +05:30
IO_globalTagInPart, &
IO_countSections, &
IO_error, &
IO_countTagInPart, &
IO_getTag, &
IO_spotTagInPart, &
IO_lc, &
IO_isBlank, &
IO_floatValue, &
IO_stringValue, &
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
IO_stringPos, &
IO_EOF
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use math, only: &
inRad, &
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
math_sampleRandomOri, &
math_I3, &
math_inv33
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
Update in order to include the c/a-ratio as additional information about the structure (modules lattice, constitutive_ changed accordingly) Correction of the expression for the velocity of mobile dislocation densities in constitutive_dislobased.f90 Correction of some Fortran syntax, incompatible with Compaq Visual Fortran (e.g. integer(8))
2009-03-20 20:04:24 +05:30
character(len=*), intent(in) :: myPart
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), intent(in) :: fileUnit
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
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 :: MAXNCHUNKS = 13_pInt
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
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), dimension(1+2*MAXNCHUNKS) :: positions
general polishing
2013-03-28 19:20:20 +05:30
integer(pInt) :: Nsections, section, gauss, fiber, j
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
character(len=65536) :: tag
character(len=65536) :: line
logical :: echo
parts of the material.config can now be set to echo their content back into the log-file (i.e. STDOUT / unit 6) use keyword "/echo/" before first section in the respective parts that you want to see echoed back into the log file.
2012-06-26 15:54:54 +05:30
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
echo = IO_globalTagInPart(fileUnit,myPart,'/echo/')
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
Nsections = IO_countSections(fileUnit,myPart)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
material_Ntexture = Nsections
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
if (Nsections < 1_pInt) call IO_error(160_pInt,ext_msg=myPart)
reworked crystallite part to allow for flexible user output --> new "crystallite" part in config file --> new "crystallite" option for microstructures --> new output file "...job.outputCrystallite" to be used in conjunction with marc_addUserOutput for meaningful naming of User Defined Vars.
2010-02-25 23:09:11 +05:30
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(texture_name(Nsections)); texture_name=''
allocate(texture_ODFfile(Nsections)); texture_ODFfile=''
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
allocate(texture_symmetry(Nsections), source=1_pInt)
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(texture_Ngauss(Nsections), source=0_pInt)
allocate(texture_Nfiber(Nsections), source=0_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
texture_Ngauss = IO_countTagInPart(fileUnit,myPart,'(gauss)', Nsections) + &
IO_countTagInPart(fileUnit,myPart,'(random)',Nsections)
texture_Nfiber = IO_countTagInPart(fileUnit,myPart,'(fiber)', Nsections)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
texture_maxNgauss = maxval(texture_Ngauss)
texture_maxNfiber = maxval(texture_Nfiber)
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(texture_Gauss (5,texture_maxNgauss,Nsections), source=0.0_pReal)
allocate(texture_Fiber (6,texture_maxNfiber,Nsections), source=0.0_pReal)
allocate(texture_transformation(3,3,Nsections), source=0.0_pReal)
texture_transformation = spread(math_I3,3,Nsections)
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
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
rewind(fileUnit)
general polishing
2013-03-28 19:20:20 +05:30
line = '' ! to have in initialized
section = 0_pInt ! - " -
gauss = 0_pInt ! - " -
fiber = 0_pInt ! - " -
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= myPart) ! wind forward to <texture>
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
line = IO_read(fileUnit)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(/,1x,a)') trim(line) ! echo part header
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
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
do while (trim(line) /= IO_EOF)
line = IO_read(fileUnit)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
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 (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
exit
endif
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 21:18:13 +05:30
if (echo) write(6,'(2x,a)') trim(line) ! echo back read lines
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
section = section + 1_pInt
gauss = 0_pInt
fiber = 0_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
texture_name(section) = IO_getTag(line,'[',']')
endif
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
if (section > 0_pInt) 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
positions = IO_stringPos(line,MAXNCHUNKS)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
general polishing
2013-03-28 19:20:20 +05:30
textureType: select case(tag)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
renamed texture_rotation into texture_transformation changed keyword for material.config to axes in line with geom_fromAng.py former keyword rotation is still recognized for compatibility
2013-07-24 16:39:39 +05:30
case ('axes', 'rotation') textureType
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
do j = 1_pInt, 3_pInt ! look for "x", "y", and "z" entries
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
tag = IO_lc(IO_stringValue(line,positions,j+1_pInt))
select case (tag)
case('x', '+x')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [ 1.0_pReal, 0.0_pReal, 0.0_pReal] ! original axis is now +x-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case('-x')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [-1.0_pReal, 0.0_pReal, 0.0_pReal] ! original axis is now -x-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case('y', '+y')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [ 0.0_pReal, 1.0_pReal, 0.0_pReal] ! original axis is now +y-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case('-y')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [ 0.0_pReal,-1.0_pReal, 0.0_pReal] ! original axis is now -y-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case('z', '+z')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [ 0.0_pReal, 0.0_pReal, 1.0_pReal] ! original axis is now +z-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case('-z')
fixed wrong initialized symmetry value
2013-12-16 16:26:56 +05:30
texture_transformation(j,1:3,section) = [ 0.0_pReal, 0.0_pReal,-1.0_pReal] ! original axis is now -z-axis
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
case default
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
call IO_error(157_pInt,section)
introduced new keyword 'rotation' to specify texture rotation ATTENTION: so far it is read but NOT yet applied!!! waiting for updated code for orientation assignement
2013-05-02 14:05:37 +05:30
end select
enddo
general polishing
2013-03-28 19:20:20 +05:30
case ('hybridia') textureType
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
texture_ODFfile(section) = IO_stringValue(line,positions,2_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
general polishing
2013-03-28 19:20:20 +05:30
case ('symmetry') textureType
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
tag = IO_lc(IO_stringValue(line,positions,2_pInt))
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case (tag)
case('orthotropic')
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
texture_symmetry(section) = 4_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('monoclinic')
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
texture_symmetry(section) = 2_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case default
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
texture_symmetry(section) = 1_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
general polishing
2013-03-28 19:20:20 +05:30
case ('(random)') textureType
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
gauss = gauss + 1_pInt
reworked restarting for compatibility with abaqus (not yet fully working) added new orientation feature for direct simulation: component type (random) asigns random orientation to an entire grain
2011-05-28 15:14:43 +05:30
texture_Gauss(1:3,gauss,section) = math_sampleRandomOri()
general polishing
2013-03-28 19:20:20 +05:30
do j = 2_pInt,4_pInt,2_pInt
tag = IO_lc(IO_stringValue(line,positions,j))
reworked restarting for compatibility with abaqus (not yet fully working) added new orientation feature for direct simulation: component type (random) asigns random orientation to an entire grain
2011-05-28 15:14:43 +05:30
select case (tag)
case('scatter')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(4,gauss,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
reworked restarting for compatibility with abaqus (not yet fully working) added new orientation feature for direct simulation: component type (random) asigns random orientation to an entire grain
2011-05-28 15:14:43 +05:30
case('fraction')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(5,gauss,section) = IO_floatValue(line,positions,j+1_pInt)
reworked restarting for compatibility with abaqus (not yet fully working) added new orientation feature for direct simulation: component type (random) asigns random orientation to an entire grain
2011-05-28 15:14:43 +05:30
end select
enddo
general polishing
2013-03-28 19:20:20 +05:30
case ('(gauss)') textureType
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
gauss = gauss + 1_pInt
general polishing
2013-03-28 19:20:20 +05:30
do j = 2_pInt,10_pInt,2_pInt
tag = IO_lc(IO_stringValue(line,positions,j))
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case (tag)
case('phi1')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(1,gauss,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('phi')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(2,gauss,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('phi2')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(3,gauss,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('scatter')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(4,gauss,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('fraction')
general polishing
2013-03-28 19:20:20 +05:30
texture_Gauss(5,gauss,section) = IO_floatValue(line,positions,j+1_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
enddo
reworked restarting for compatibility with abaqus (not yet fully working) added new orientation feature for direct simulation: component type (random) asigns random orientation to an entire grain
2011-05-28 15:14:43 +05:30
general polishing
2013-03-28 19:20:20 +05:30
case ('(fiber)') textureType
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
fiber = fiber + 1_pInt
general polishing
2013-03-28 19:20:20 +05:30
do j = 2_pInt,12_pInt,2_pInt
tag = IO_lc(IO_stringValue(line,positions,j))
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
select case (tag)
corrected spelling in subroutine "material_parseTexture"
2009-04-28 15:15:52 +05:30
case('alpha1')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(1,fiber,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('alpha2')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(2,fiber,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('beta1')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(3,fiber,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('beta2')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(4,fiber,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('scatter')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(5,fiber,section) = IO_floatValue(line,positions,j+1_pInt)*inRad
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
case('fraction')
general polishing
2013-03-28 19:20:20 +05:30
texture_Fiber(6,fiber,section) = IO_floatValue(line,positions,j+1_pInt)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
end select
enddo
general polishing
2013-03-28 19:20:20 +05:30
end select textureType
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
enddo
added possibility for multi-level inclusion of files in *.config and loadcase files. include subfiles by stating {path/to/include}
2013-06-27 00:49:00 +05:30
end subroutine material_parseTexture
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief populates the grains
!> @details populates the grains by identifying active microstructure/homogenization pairs,
!! calculates the volume of the grains and deals with texture components and hybridIA
!--------------------------------------------------------------------------------------------------
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 material_populateGrains
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
use math, only: &
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
math_RtoEuler, &
math_EulerToR, &
math_mul33x33, &
math_range, &
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
math_sampleRandomOri, &
math_sampleGaussOri, &
math_sampleFiberOri, &
math_symmetricEulers
use mesh, only: &
mesh_element, &
mesh_maxNips, &
mesh_NcpElems, &
mesh_ipVolume, &
FE_Nips, &
FE_geomtype
use IO, only: &
IO_error, &
IO_hybridIA
use FEsolving, only: &
FEsolving_execIP
use debug, only: &
debug_level, &
debug_material, &
debug_levelBasic
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
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
implicit none
integer(pInt), dimension (:,:), allocatable :: Ngrains
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
integer(pInt), dimension (microstructure_maxNconstituents) :: &
NgrainsOfConstituent, &
currentGrainOfConstituent, &
randomOrder
real(pReal), dimension (microstructure_maxNconstituents) :: &
rndArray
(1) change phaseOfGrain from real to integer and (2) introducing a condition to switch off the grain reshuffling when element homogeneous is requested.
2010-02-18 21:24:10 +05:30
real(pReal), dimension (:), allocatable :: volumeOfGrain
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
real(pReal), dimension (:,:), allocatable :: orientationOfGrain
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
real(pReal), dimension (3) :: orientation
real(pReal), dimension (3,3) :: symOrientation
added "texture " as possible crystallite output => reports ID for later discrimination...
2011-03-22 19:10:27 +05:30
integer(pInt), dimension (:), allocatable :: phaseOfGrain, textureOfGrain
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
integer(pInt) :: t,e,i,g,j,m,c,r,homog,micro,sgn,hme, myDebug, &
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
phaseID,textureID,dGrains,myNgrains,myNorientations,myNconstituents, &
moved field initialization from homogenization to material since it is needed in crystallite
2014-09-23 02:06:55 +05:30
grain,constituentGrain,ipGrain,symExtension, ip
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
real(pReal) :: extreme,rnd
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
integer(pInt), dimension (:,:), allocatable :: Nelems ! counts number of elements in homog, micro array
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
type(p_intvec), dimension (:,:), allocatable :: elemsOfHomogMicro ! lists element number in homog, micro array
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
changed variable name "debug_what" to "debug_level"
2012-07-05 15:24:50 +05:30
myDebug = debug_level(debug_material)
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
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(material_volume(homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems), source=0.0_pReal)
allocate(material_phase(homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems), source=0_pInt)
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
allocate(material_homog(mesh_maxNips,mesh_NcpElems), source=0_pInt)
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(material_texture(homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems), source=0_pInt)
allocate(material_EulerAngles(3,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems),source=0.0_pReal)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(Ngrains(material_Nhomogenization,material_Nmicrostructure), source=0_pInt)
allocate(Nelems(material_Nhomogenization,material_Nmicrostructure), source=0_pInt)
removed spectralPicturemode
2012-01-12 16:06:17 +05:30
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
! populating homogenization schemes in each
!--------------------------------------------------------------------------------------------------
do e = 1_pInt, mesh_NcpElems
prepared for using new state structure in homogenization
2014-09-18 20:57:52 +05:30
material_homog(1_pInt:FE_Nips(FE_geomtype(mesh_element(2,e))),e) = mesh_element(3,e)
started introducing new state structure in homogenisation
2014-08-21 23:18:20 +05:30
enddo
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
bugfix release
2012-01-12 22:01:23 +05:30
! precounting of elements for each homog/micro pair
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
do e = 1_pInt, mesh_NcpElems
bugfix release
2012-01-12 22:01:23 +05:30
homog = mesh_element(3,e)
micro = mesh_element(4,e)
Nelems(homog,micro) = Nelems(homog,micro) + 1_pInt
removed spectralPicturemode
2012-01-12 16:06:17 +05:30
enddo
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
allocate(elemsOfHomogMicro(material_Nhomogenization,material_Nmicrostructure))
do homog = 1,material_Nhomogenization
do micro = 1,material_Nmicrostructure
if (Nelems(homog,micro) > 0_pInt) then
allocate(elemsOfHomogMicro(homog,micro)%p(Nelems(homog,micro)))
elemsOfHomogMicro(homog,micro)%p = 0_pInt
endif
enddo
enddo
bugfix release
2012-01-12 22:01:23 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
bugfix release
2012-01-12 22:01:23 +05:30
! identify maximum grain count per IP (from element) and find grains per homog/micro pair
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
Nelems = 0_pInt ! reuse as counter
general polishing
2013-03-28 19:20:20 +05:30
elementLooping: do e = 1_pInt,mesh_NcpElems
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
t = FE_geomtype(mesh_element(2,e))
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
homog = mesh_element(3,e)
micro = mesh_element(4,e)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (homog < 1_pInt .or. homog > material_Nhomogenization) & ! out of bounds
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
call IO_error(154_pInt,e,0_pInt,0_pInt)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (micro < 1_pInt .or. micro > material_Nmicrostructure) & ! out of bounds
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore. all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 23:11:27 +05:30
call IO_error(155_pInt,e,0_pInt,0_pInt)
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
if (microstructure_elemhomo(micro)) then ! how many grains are needed at this element?
dGrains = homogenization_Ngrains(homog) ! only one set of Ngrains (other IPs are plain copies)
introduction of a possibility to have homogeneous element (all ips in one element are identical, sort of reduced integration) and bugs-fixing in crystallite.f90, homogenization_RGC.f90, numerics.f90.
2009-11-24 20:30:25 +05:30
else
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
dGrains = homogenization_Ngrains(homog) * FE_Nips(t) ! each IP has Ngrains
introduction of a possibility to have homogeneous element (all ips in one element are identical, sort of reduced integration) and bugs-fixing in crystallite.f90, homogenization_RGC.f90, numerics.f90.
2009-11-24 20:30:25 +05:30
endif
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
Ngrains(homog,micro) = Ngrains(homog,micro) + dGrains ! total grain count
Nelems(homog,micro) = Nelems(homog,micro) + 1_pInt ! total element count
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
elemsOfHomogMicro(homog,micro)%p(Nelems(homog,micro)) = e ! remember elements active in this homog/micro pair
general polishing
2013-03-28 19:20:20 +05:30
enddo elementLooping
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
using sourced allocation
2013-12-13 04:33:37 +05:30
allocate(volumeOfGrain(maxval(Ngrains)), source=0.0_pReal) ! reserve memory for maximum case
allocate(phaseOfGrain(maxval(Ngrains)), source=0_pInt) ! reserve memory for maximum case
allocate(textureOfGrain(maxval(Ngrains)), source=0_pInt) ! reserve memory for maximum case
allocate(orientationOfGrain(3,maxval(Ngrains)),source=0.0_pReal) ! reserve memory for maximum case
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
if (iand(myDebug,debug_levelBasic) /= 0_pInt) then
debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output) 0 : only version infos and all from "hypela2"/"umat" 1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info 2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines 3 : basic outputs from "homogenization.f90" 4 : extensive outputs from "homogenization.f90" 5 : basic outputs from "crystallite.f90" 6 : extensive outputs from "crystallite.f90" 7 : basic outputs from the constitutive files 8 : extensive outputs from the constitutive files If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 16:01:17 +05:30
!$OMP CRITICAL (write2out)
general polishing
2013-03-28 19:20:20 +05:30
write(6,'(/,a/)') ' MATERIAL grain population'
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 18:23:25 +05:30
write(6,'(a32,1x,a32,1x,a6)') 'homogenization_name','microstructure_name','grain#'
debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output) 0 : only version infos and all from "hypela2"/"umat" 1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info 2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines 3 : basic outputs from "homogenization.f90" 4 : extensive outputs from "homogenization.f90" 5 : basic outputs from "crystallite.f90" 6 : extensive outputs from "crystallite.f90" 7 : basic outputs from the constitutive files 8 : extensive outputs from the constitutive files If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 16:01:17 +05:30
!$OMP END CRITICAL (write2out)
endif
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
do homog = 1_pInt,material_Nhomogenization ! loop over homogenizations
dGrains = homogenization_Ngrains(homog) ! grain number per material point
do micro = 1_pInt,material_Nmicrostructure ! all pairs of homog and micro
if (Ngrains(homog,micro) > 0_pInt) then ! an active pair of homog and micro
myNgrains = Ngrains(homog,micro) ! assign short name for total number of grains to populate
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
myNconstituents = microstructure_Nconstituents(micro) ! assign short name for number of constituents
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
if (iand(myDebug,debug_levelBasic) /= 0_pInt) then
debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output) 0 : only version infos and all from "hypela2"/"umat" 1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info 2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines 3 : basic outputs from "homogenization.f90" 4 : extensive outputs from "homogenization.f90" 5 : basic outputs from "crystallite.f90" 6 : extensive outputs from "crystallite.f90" 7 : basic outputs from the constitutive files 8 : extensive outputs from the constitutive files If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 16:01:17 +05:30
!$OMP CRITICAL (write2out)
general polishing
2013-03-28 19:20:20 +05:30
write(6,'(/,a32,1x,a32,1x,i6)') homogenization_name(homog),microstructure_name(micro),myNgrains
debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output) 0 : only version infos and all from "hypela2"/"umat" 1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info 2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines 3 : basic outputs from "homogenization.f90" 4 : extensive outputs from "homogenization.f90" 5 : basic outputs from "crystallite.f90" 6 : extensive outputs from "crystallite.f90" 7 : basic outputs from the constitutive files 8 : extensive outputs from the constitutive files If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 16:01:17 +05:30
!$OMP END CRITICAL (write2out)
endif
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
! calculate volume of each grain
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
same as for rev327
2009-05-26 22:54:42 +05:30
volumeOfGrain = 0.0_pReal
removed spectralPicturemode
2012-01-12 16:06:17 +05:30
grain = 0_pInt
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
bugfix release
2012-01-12 22:01:23 +05:30
do hme = 1_pInt, Nelems(homog,micro)
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
e = elemsOfHomogMicro(homog,micro)%p(hme) ! my combination of homog and micro, only perform calculations for elements with homog, micro combinations which is indexed in cpElemsindex
switched element library to geomType based. saves to copy same geometry description for different elements that are essentially similar regarding the IP number but differ in total node count. introduced quadratic tetrahedron (Marc element 127 -- element 157 might also work, but did not perform well in fully elastic calc so far)
2012-11-16 04:15:20 +05:30
t = FE_geomtype(mesh_element(2,e))
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (microstructure_elemhomo(micro)) then ! homogeneous distribution of grains over each element's IPs
switched element library to geomType based. saves to copy same geometry description for different elements that are essentially similar regarding the IP number but differ in total node count. introduced quadratic tetrahedron (Marc element 127 -- element 157 might also work, but did not perform well in fully elastic calc so far)
2012-11-16 04:15:20 +05:30
volumeOfGrain(grain+1_pInt:grain+dGrains) = sum(mesh_ipVolume(1:FE_Nips(t),e))/&
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
real(dGrains,pReal) ! each grain combines size of all IPs in that element
grain = grain + dGrains ! wind forward by Ngrains@IP
removed spectralPicturemode
2012-01-12 16:06:17 +05:30
else
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
forall (i = 1_pInt:FE_Nips(t)) & ! loop over IPs
polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard. removed cut_off parameter for damask_spectral removed outpot of derived divergence measures and added RMS output in brackets added comments and options to the makefile
2012-02-16 00:28:38 +05:30
volumeOfGrain(grain+(i-1)*dGrains+1_pInt:grain+i*dGrains) = &
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
mesh_ipVolume(i,e)/dGrains ! assign IPvolume/Ngrains@IP to all grains of IP
grain = grain + FE_Nips(t) * dGrains ! wind forward by Nips*Ngrains@IP
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
enddo
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
if (grain /= myNgrains) &
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_error(0,el = homog,ip = micro,ext_msg = 'inconsistent grain count after volume calc')
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
! divide myNgrains as best over constituents
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
!
! example: three constituents with fractions of 0.25, 0.25, and 0.5 distributed over 20 (microstructure) grains
!
! ***** ***** **********
! NgrainsOfConstituent: 5, 5, 10
! counters:
! |-----> grain (if constituent == 2)
! |--> constituentGrain (of constituent 2)
!
NgrainsOfConstituent = 0_pInt ! reset counter of grains per constituent
forall (i = 1_pInt:myNconstituents) &
general polishing
2013-03-28 19:20:20 +05:30
NgrainsOfConstituent(i) = nint(microstructure_fraction(i,micro) * myNgrains, pInt) ! do rounding integer conversion
do while (sum(NgrainsOfConstituent) /= myNgrains) ! total grain count over constituents wrong?
sgn = sign(1_pInt, myNgrains - sum(NgrainsOfConstituent)) ! direction of required change
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
extreme = 0.0_pReal
t = 0_pInt
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
do i = 1_pInt,myNconstituents ! find largest deviator
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
if (real(sgn,pReal)*log(NgrainsOfConstituent(i)/myNgrains/microstructure_fraction(i,micro)) > extreme) then
extreme = real(sgn,pReal)*log(NgrainsOfConstituent(i)/myNgrains/microstructure_fraction(i,micro))
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
t = i
endif
enddo
general polishing
2013-03-28 19:20:20 +05:30
NgrainsOfConstituent(t) = NgrainsOfConstituent(t) + sgn ! change that by one
added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
2009-06-15 18:41:21 +05:30
enddo
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
!--------------------------------------------------------------------------------------------------
! assign phase and texture info
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
phaseOfGrain = 0_pInt
added "texture " as possible crystallite output => reports ID for later discrimination...
2011-03-22 19:10:27 +05:30
textureOfGrain = 0_pInt
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
orientationOfGrain = 0.0_pReal
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
texture: do i = 1_pInt,myNconstituents ! loop over constituents
grain = sum(NgrainsOfConstituent(1_pInt:i-1_pInt)) ! set microstructure grain index of current constituent
! "grain" points to start of this constituent's grain population
constituentGrain = 0_pInt ! constituent grain index
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
phaseID = microstructure_phase(i,micro)
textureID = microstructure_texture(i,micro)
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
phaseOfGrain (grain+1_pInt:grain+NgrainsOfConstituent(i)) = phaseID ! assign resp. phase
general polishing
2013-03-28 19:20:20 +05:30
textureOfGrain(grain+1_pInt:grain+NgrainsOfConstituent(i)) = textureID ! assign resp. texture
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
fixing inconsistencies in variable assignments
2012-02-10 17:26:05 +05:30
myNorientations = ceiling(real(NgrainsOfConstituent(i),pReal)/&
general polishing
2013-03-28 19:20:20 +05:30
real(texture_symmetry(textureID),pReal),pInt) ! max number of unique orientations (excl. symmetry)
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
! ...has texture components
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (texture_ODFfile(textureID) == '') then
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
gauss: do t = 1_pInt,texture_Ngauss(textureID) ! loop over Gauss components
general polishing
2013-03-28 19:20:20 +05:30
do g = 1_pInt,int(myNorientations*texture_Gauss(5,t,textureID),pInt) ! loop over required grain count
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
orientationOfGrain(:,grain+constituentGrain+g) = &
math_sampleGaussOri(texture_Gauss(1:3,t,textureID),&
texture_Gauss( 4,t,textureID))
enddo
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
constituentGrain = &
constituentGrain + int(myNorientations*texture_Gauss(5,t,textureID)) ! advance counter for grains of current constituent
enddo gauss
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
fiber: do t = 1_pInt,texture_Nfiber(textureID) ! loop over fiber components
general polishing
2013-03-28 19:20:20 +05:30
do g = 1_pInt,int(myNorientations*texture_Fiber(6,t,textureID),pInt) ! loop over required grain count
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
orientationOfGrain(:,grain+constituentGrain+g) = &
math_sampleFiberOri(texture_Fiber(1:2,t,textureID),&
texture_Fiber(3:4,t,textureID),&
texture_Fiber( 5,t,textureID))
enddo
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
constituentGrain = &
constituentGrain + int(myNorientations*texture_fiber(6,t,textureID),pInt) ! advance counter for grains of current constituent
enddo fiber
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
random: do constituentGrain = constituentGrain+1_pInt,myNorientations ! fill remainder with random
orientationOfGrain(:,grain+constituentGrain) = math_sampleRandomOri()
enddo random
!--------------------------------------------------------------------------------------------------
! ...has hybrid IA
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
else
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
orientationOfGrain(1:3,grain+1_pInt:grain+myNorientations) = &
IO_hybridIA(myNorientations,texture_ODFfile(textureID))
if (all(orientationOfGrain(1:3,grain+1_pInt) == -1.0_pReal)) call IO_error(156_pInt)
endif
general polishing
2013-03-28 19:20:20 +05:30
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
!--------------------------------------------------------------------------------------------------
renamed texture_rotation into texture_transformation changed keyword for material.config to axes in line with geom_fromAng.py former keyword rotation is still recognized for compatibility
2013-07-24 16:39:39 +05:30
! ...texture transformation
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
do j = 1_pInt,myNorientations ! loop over each "real" orientation
orientationOfGrain(1:3,grain+j) = math_RtoEuler( & ! translate back to Euler angles
math_mul33x33( & ! pre-multiply
math_EulertoR(orientationOfGrain(1:3,grain+j)), & ! face-value orientation
renamed texture_rotation into texture_transformation changed keyword for material.config to axes in line with geom_fromAng.py former keyword rotation is still recognized for compatibility
2013-07-24 16:39:39 +05:30
texture_transformation(1:3,1:3,textureID) & ! and transformation matrix
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
) &
)
enddo
!--------------------------------------------------------------------------------------------------
! ...sample symmetry
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
symExtension = texture_symmetry(textureID) - 1_pInt
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
if (symExtension > 0_pInt) then ! sample symmetry (number of additional equivalent orientations)
constituentGrain = myNorientations ! start right after "real" orientations
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
do j = 1_pInt,myNorientations ! loop over each "real" orientation
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
symOrientation = math_symmetricEulers(texture_symmetry(textureID), &
orientationOfGrain(1:3,grain+j)) ! get symmetric equivalents
e = min(symExtension,NgrainsOfConstituent(i)-constituentGrain) ! do not overshoot end of constituent grain array
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
if (e > 0_pInt) then
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
orientationOfGrain(1:3,grain+constituentGrain+1: &
grain+constituentGrain+e) = &
symOrientation(1:3,1:e)
constituentGrain = constituentGrain + e ! remainder shrinks by e
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
enddo
endif
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
general polishing
2013-03-28 19:20:20 +05:30
!--------------------------------------------------------------------------------------------------
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
! shuffle grains within current constituent
do j = 1_pInt,NgrainsOfConstituent(i)-1_pInt ! walk thru grains of current constituent
added "texture " as possible crystallite output => reports ID for later discrimination...
2011-03-22 19:10:27 +05:30
call random_number(rnd)
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
t = nint(rnd*(NgrainsOfConstituent(i)-j)+j+0.5_pReal,pInt) ! select a grain in remaining list
m = phaseOfGrain(grain+t) ! exchange current with random
phaseOfGrain(grain+t) = phaseOfGrain(grain+j)
phaseOfGrain(grain+j) = m
m = textureOfGrain(grain+t) ! exchange current with random
textureOfGrain(grain+t) = textureOfGrain(grain+j)
textureOfGrain(grain+j) = m
orientation = orientationOfGrain(1:3,grain+t) ! exchange current with random
orientationOfGrain(1:3,grain+t) = orientationOfGrain(1:3,grain+j)
orientationOfGrain(1:3,grain+j) = orientation
added "texture " as possible crystallite output => reports ID for later discrimination...
2011-03-22 19:10:27 +05:30
enddo
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
enddo texture
fixed bug in include_linux64 for Marc 2012, tests should run again after updating the copy in /msc/.... polished two other versions of this file unified output format for /echo/ in material.f90
2013-06-11 12:58:08 +05:30
!< @todo calc fraction after weighing with volumePerGrain, exchange in MC steps to improve result (humbug at the moment)
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
!--------------------------------------------------------------------------------------------------
! distribute grains of all constituents as accurately as possible to given constituent fractions
ip = 0_pInt
currentGrainOfConstituent = 0_pInt
bugfix release
2012-01-12 22:01:23 +05:30
do hme = 1_pInt, Nelems(homog,micro)
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
e = elemsOfHomogMicro(homog,micro)%p(hme) ! only perform calculations for elements with homog, micro combinations which is indexed in cpElemsindex
polishing
2012-11-21 22:27:57 +05:30
t = FE_geomtype(mesh_element(2,e))
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
if (microstructure_elemhomo(micro)) then ! homogeneous distribution of grains over each element's IPs
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
m = 1_pInt ! process only first IP
removed spectralPicturemode
2012-01-12 16:06:17 +05:30
else
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
m = FE_Nips(t) ! process all IPs
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
endif
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
do i = 1_pInt, m ! loop over necessary IPs
ip = ip + 1_pInt ! keep track of total ip count
ipGrain = 0_pInt ! count number of grains assigned at this IP
randomOrder = math_range(microstructure_maxNconstituents) ! start out with ordered sequence of constituents
call random_number(rndArray) ! as many rnd numbers as (max) constituents
do j = 1_pInt, myNconstituents - 1_pInt ! loop over constituents ...
r = nint(rndArray(j)*(myNconstituents-j)+j+0.5_pReal,pInt) ! ... select one in remaining list
c = randomOrder(r) ! ... call it "c"
randomOrder(r) = randomOrder(j) ! ... and exchange with present position in constituent list
grain = sum(NgrainsOfConstituent(1:c-1_pInt)) ! figure out actual starting index in overall/consecutive grain population
do g = 1_pInt, min(dGrains-ipGrain, & ! leftover number of grains at this IP
max(0_pInt, & ! no negative values
nint(real(ip * dGrains * NgrainsOfConstituent(c)) / & ! fraction of grains scaled to this constituent...
real(myNgrains),pInt) - & ! ...minus those already distributed
currentGrainOfConstituent(c)))
ipGrain = ipGrain + 1_pInt ! advance IP grain counter
currentGrainOfConstituent(c) = currentGrainOfConstituent(c) + 1_pInt ! advance index of grain population for constituent c
material_volume(ipGrain,i,e) = volumeOfGrain(grain+currentGrainOfConstituent(c)) ! assign properties
material_phase(ipGrain,i,e) = phaseOfGrain(grain+currentGrainOfConstituent(c))
material_texture(ipGrain,i,e) = textureOfGrain(grain+currentGrainOfConstituent(c))
material_EulerAngles(1:3,ipGrain,i,e) = orientationOfGrain(1:3,grain+currentGrainOfConstituent(c))
enddo; enddo
c = randomOrder(microstructure_Nconstituents(micro)) ! look up constituent remaining after random shuffling
grain = sum(NgrainsOfConstituent(1:c-1_pInt)) ! figure out actual starting index in overall/consecutive grain population
do ipGrain = ipGrain + 1_pInt, dGrains ! ensure last constituent fills up to dGrains
currentGrainOfConstituent(c) = currentGrainOfConstituent(c) + 1_pInt
material_volume(ipGrain,i,e) = volumeOfGrain(grain+currentGrainOfConstituent(c))
material_phase(ipGrain,i,e) = phaseOfGrain(grain+currentGrainOfConstituent(c))
material_texture(ipGrain,i,e) = textureOfGrain(grain+currentGrainOfConstituent(c))
material_EulerAngles(1:3,ipGrain,i,e) = orientationOfGrain(1:3,grain+currentGrainOfConstituent(c))
enddo
enddo
do i = i, FE_Nips(t) ! loop over IPs to (possibly) distribute copies from first IP
material_volume (1_pInt:dGrains,i,e) = material_volume (1_pInt:dGrains,1,e)
material_phase (1_pInt:dGrains,i,e) = material_phase (1_pInt:dGrains,1,e)
material_texture(1_pInt:dGrains,i,e) = material_texture(1_pInt:dGrains,1,e)
material_EulerAngles(1:3,1_pInt:dGrains,i,e) = material_EulerAngles(1:3,1_pInt:dGrains,1,e)
enddo
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 17:00:52 +05:30
endif ! active homog,micro pair
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
enddo
enddo
same as for rev327
2009-05-26 22:54:42 +05:30
deallocate(volumeOfGrain)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
deallocate(phaseOfGrain)
added "texture " as possible crystallite output => reports ID for later discrimination...
2011-03-22 19:10:27 +05:30
deallocate(textureOfGrain)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
deallocate(orientationOfGrain)
bugfix release
2012-01-12 22:01:23 +05:30
deallocate(Nelems)
made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables
2013-06-11 22:05:04 +05:30
!> @todo - causing segmentation fault: needs looking into
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
!do homog = 1,material_Nhomogenization
! do micro = 1,material_Nmicrostructure
reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents. removed error for volume fractions not equalling 1. implemented capability to rotate the texture given in material.config.
2013-05-29 22:53:49 +05:30
! if (Nelems(homog,micro) > 0_pInt) deallocate(elemsOfHomogMicro(homog,micro)%p)
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 18:43:52 +05:30
! enddo
!enddo
bugfix release
2012-01-12 22:01:23 +05:30
deallocate(elemsOfHomogMicro)
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
end subroutine material_populateGrains
overhaul to allow for multiphase/homogenization setup.
2009-03-04 17:18:54 +05:30
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
#ifdef HDF
small changes to HDF5, put link to homepage into README not to have double effort
2014-03-25 22:51:47 +05:30
integer(pInt) pure function material_NconstituentsPhase(matID)
implicit none
integer(pInt), intent(in) :: matID
material_NconstituentsPhase = count(microstructure_phase == matID)
end function
some changes related to new state fixed assigning of wrong output in J2 model
2014-05-08 20:25:19 +05:30
#endif
small changes to HDF5, put link to homepage into README not to have double effort
2014-03-25 22:51:47 +05:30
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-20 03:42:05 +05:30
end module material