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

2019 lines
117 KiB
Fortran
Raw Normal View History

comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
! $Id$
!--------------------------------------------------------------------------------------------------
!> @author Alankar Alankar, Max-Planck-Institut für Eisenforschung GmbH
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief material subroutine for titanium
!--------------------------------------------------------------------------------------------------
added "public" (and where possible "protected") statements to functions/subroutines/variables/parameters where needed added -fmodule-private to compile options for gfortran, i.e. gfortran won't compile anymore if an entity is not explicitly declared as public but used in another module to ensure stronger encapsulation
2013-03-28 13:10:30 +05:30
module constitutive_titanmod
use prec, only: &
pReal, &
pInt
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
implicit none
private
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
character(len=18), dimension(3), parameter, private :: &
CONSTITUTIVE_TITANMOD_listBasicSlipStates = &
['rho_edge ', 'rho_screw ', 'shear_system']
character(len=18), dimension(1), parameter, private :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
CONSTITUTIVE_TITANMOD_listBasicTwinStates = ['gdot_twin']
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
character(len=19), dimension(11), parameter, private :: &
CONSTITUTIVE_TITANMOD_listDependentSlipStates = &
['segment_edge ', 'segment_screw ', &
'resistance_edge ', 'resistance_screw ', &
'tau_slip ', &
'velocity_edge ', 'velocity_screw ', &
'gdot_slip_edge ', 'gdot_slip_screw ', &
'stressratio_edge_p ', 'stressratio_screw_p' ]
character(len=18), dimension(2), parameter, private :: &
constitutive_titanmod_listDependentTwinStates = &
['twin_fraction', 'tau_twin ']
real(pReal), parameter, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
kB = 1.38e-23_pReal !< Boltzmann constant in J/Kelvin
comments for doxygen
2013-07-12 12:27:15 +05:30
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
integer(pInt), dimension(:), allocatable, public, protected :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_sizePostResults !< cumulative size of post results
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
comments for doxygen
2013-07-12 12:27:15 +05:30
integer(pInt), dimension(:,:), allocatable, target, public :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_sizePostResult !< size of each post result output
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
comments for doxygen
2013-07-12 12:27:15 +05:30
character(len=64), dimension(:,:), allocatable, target, public :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_output !< name of each post result output
added new, flexible debugging scheme. now all modules have their own debug specification. compiles and runs, I hope nothing is broken did a lot of polishing
2012-03-09 01:55:28 +05:30
use constitutive_xxx_Noutput instead of phase_Noutput as number of outputs for more consistency with constitutive multi physics output
2014-09-26 15:55:26 +05:30
integer(pInt), dimension(:), allocatable, target, public :: &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_Noutput !< number of outputs per instance of this plasticity !< ID of the lattice structure
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 and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt), dimension(:), allocatable, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_totalNslip, & !< total number of active slip systems for each instance
constitutive_titanmod_totalNtwin !< total number of active twin systems for each instance
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 and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt), dimension(:,:), allocatable, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_Nslip, & !< number of active slip systems for each family and instance
constitutive_titanmod_Ntwin, & !< number of active twin systems for each family and instance
constitutive_titanmod_slipFamily, & !< lookup table relating active slip system to slip family for each instance
constitutive_titanmod_twinFamily, & !< lookup table relating active twin system to twin family for each instance
constitutive_titanmod_slipSystemLattice, & !< lookup table relating active slip system index to lattice slip system index for each instance
constitutive_titanmod_twinSystemLattice !< lookup table relating active twin system index to lattice twin system index for each instance
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 and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(:), allocatable, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_debyefrequency, & !< Debye frequency
constitutive_titanmod_kinkf0, & !<
constitutive_titanmod_CAtomicVolume, & !< atomic volume in Bugers vector unit
constitutive_titanmod_dc, & !< prefactor for self-diffusion coefficient
constitutive_titanmod_twinhpconstant, & !< activation energy for dislocation climb
constitutive_titanmod_GrainSize, & !< grain size - Not being used
constitutive_titanmod_MaxTwinFraction, & !< maximum allowed total twin volume fraction
constitutive_titanmod_r, & !< r-exponent in twin nucleation rate
constitutive_titanmod_CEdgeDipMinDistance, & !< Not being used
constitutive_titanmod_Cmfptwin, & !< Not being used
constitutive_titanmod_Cthresholdtwin, & !< Not being used
constitutive_titanmod_aTolRho !< absolute tolerance for integration of dislocation density
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 and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(:,:), allocatable, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_rho_edge0, & !< initial edge dislocation density per slip system for each family and instance
constitutive_titanmod_rho_screw0, & !< initial screw dislocation density per slip system for each family and instance
constitutive_titanmod_shear_system0, & !< accumulated shear on each system
constitutive_titanmod_burgersPerSlipFam, & !< absolute length of burgers vector [m] for each slip family and instance
constitutive_titanmod_burgersPerSlipSys, & !< absolute length of burgers vector [m] for each slip system and instance
constitutive_titanmod_burgersPerTwinFam, & !< absolute length of burgers vector [m] for each twin family and instance
constitutive_titanmod_burgersPerTwinSys, & !< absolute length of burgers vector [m] for each twin system and instance
constitutive_titanmod_f0_PerSlipFam, & !< activation energy for glide [J] for each slip family and instance
constitutive_titanmod_f0_PerSlipSys, & !< activation energy for glide [J] for each slip system and instance
constitutive_titanmod_twinf0_PerTwinFam, & !< activation energy for glide [J] for each twin family and instance
constitutive_titanmod_twinf0_PerTwinSys, & !< activation energy for glide [J] for each twin system and instance
constitutive_titanmod_twinshearconstant_PerTwinFam, & !< activation energy for glide [J] for each twin family and instance
constitutive_titanmod_twinshearconstant_PerTwinSys, & !< activation energy for glide [J] for each twin system and instance
constitutive_titanmod_tau0e_PerSlipFam, & !< Initial yield stress for edge dislocations per slip family
constitutive_titanmod_tau0e_PerSlipSys, & !< Initial yield stress for edge dislocations per slip system
constitutive_titanmod_tau0s_PerSlipFam, & !< Initial yield stress for screw dislocations per slip family
constitutive_titanmod_tau0s_PerSlipSys, & !< Initial yield stress for screw dislocations per slip system
constitutive_titanmod_twintau0_PerTwinFam, & !< Initial yield stress for edge dislocations per twin family
constitutive_titanmod_twintau0_PerTwinSys, & !< Initial yield stress for edge dislocations per twin system
constitutive_titanmod_capre_PerSlipFam, & !< Capture radii for edge dislocations per slip family
constitutive_titanmod_capre_PerSlipSys, & !< Capture radii for edge dislocations per slip system
constitutive_titanmod_caprs_PerSlipFam, & !< Capture radii for screw dislocations per slip family
constitutive_titanmod_caprs_PerSlipSys, & !< Capture radii for screw dislocations per slip system
constitutive_titanmod_pe_PerSlipFam, & !< p-exponent in glide velocity
constitutive_titanmod_ps_PerSlipFam, & !< p-exponent in glide velocity
constitutive_titanmod_qe_PerSlipFam, & !< q-exponent in glide velocity
constitutive_titanmod_qs_PerSlipFam, & !< q-exponent in glide velocity
constitutive_titanmod_pe_PerSlipSys, & !< p-exponent in glide velocity
constitutive_titanmod_ps_PerSlipSys, & !< p-exponent in glide velocity
constitutive_titanmod_qe_PerSlipSys, & !< q-exponent in glide velocity
constitutive_titanmod_qs_PerSlipSys, & !< q-exponent in glide velocity
constitutive_titanmod_twinp_PerTwinFam, & !< p-exponent in glide velocity
constitutive_titanmod_twinq_PerTwinFam, & !< q-exponent in glide velocity
constitutive_titanmod_twinp_PerTwinSys, & !< p-exponent in glide velocity
constitutive_titanmod_twinq_PerTwinSys, & !< p-exponent in glide velocity
constitutive_titanmod_v0e_PerSlipFam, & !< edge dislocation velocity prefactor [m/s] for each family and instance
constitutive_titanmod_v0e_PerSlipSys, & !< screw dislocation velocity prefactor [m/s] for each slip system and instance
constitutive_titanmod_v0s_PerSlipFam, & !< edge dislocation velocity prefactor [m/s] for each family and instance
constitutive_titanmod_v0s_PerSlipSys, & !< screw dislocation velocity prefactor [m/s] for each slip system and instance
constitutive_titanmod_twingamma0_PerTwinFam, & !< edge dislocation velocity prefactor [m/s] for each family and instance
constitutive_titanmod_twingamma0_PerTwinSys, & !< screw dislocation velocity prefactor [m/s] for each slip system and instance
constitutive_titanmod_kinkcriticallength_PerSlipFam, & !< screw dislocation mobility prefactor for kink-pairs per slip family
constitutive_titanmod_kinkcriticallength_PerSlipSys, & !< screw dislocation mobility prefactor for kink-pairs per slip system
constitutive_titanmod_twinsizePerTwinFam, & !< twin thickness [m] for each twin family and instance
constitutive_titanmod_twinsizePerTwinSys, & !< twin thickness [m] for each twin system and instance
constitutive_titanmod_CeLambdaSlipPerSlipFam, & !< Adj. parameter for distance between 2 forest dislocations for each slip family and instance
constitutive_titanmod_CeLambdaSlipPerSlipSys, & !< Adj. parameter for distance between 2 forest dislocations for each slip system and instance
constitutive_titanmod_CsLambdaSlipPerSlipFam, & !< Adj. parameter for distance between 2 forest dislocations for each slip family and instance
constitutive_titanmod_CsLambdaSlipPerSlipSys, & !< Adj. parameter for distance between 2 forest dislocations for each slip system and instance
constitutive_titanmod_twinLambdaSlipPerTwinFam, & !< Adj. parameter for distance between 2 forest dislocations for each slip family and instance
constitutive_titanmod_twinLambdaSlipPerTwinSys, & !< Adj. parameter for distance between 2 forest dislocations for each slip system and instance
constitutive_titanmod_interactionSlipSlip, & !< coefficients for slip-slip interaction for each interaction type and instance
constitutive_titanmod_interaction_ee, & !< coefficients for e-e interaction for each interaction type and instance
constitutive_titanmod_interaction_ss, & !< coefficients for s-s interaction for each interaction type and instance
constitutive_titanmod_interaction_es, & !< coefficients for e-s-twin interaction for each interaction type and instance
constitutive_titanmod_interactionSlipTwin, & !< coefficients for twin-slip interaction for each interaction type and instance
constitutive_titanmod_interactionTwinSlip, & !< coefficients for twin-slip interaction for each interaction type and instance
constitutive_titanmod_interactionTwinTwin !< coefficients for twin-twin interaction for each interaction type and instance
comments for doxygen
2013-07-12 12:27:15 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(:,:,:), allocatable, private :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
constitutive_titanmod_interactionMatrixSlipSlip, & !< interaction matrix of the different slip systems for each instance
constitutive_titanmod_interactionMatrix_ee, & !< interaction matrix of e-e for each instance
constitutive_titanmod_interactionMatrix_ss, & !< interaction matrix of s-s for each instance
constitutive_titanmod_interactionMatrix_es, & !< interaction matrix of e-s for each instance
constitutive_titanmod_interactionMatrixSlipTwin, & !< interaction matrix of slip systems with twin systems for each instance
constitutive_titanmod_interactionMatrixTwinSlip, & !< interaction matrix of twin systems with slip systems for each instance
constitutive_titanmod_interactionMatrixTwinTwin, & !< interaction matrix of the different twin systems for each instance
constitutive_titanmod_forestProjectionEdge, & !< matrix of forest projections of edge dislocations for each instance
constitutive_titanmod_forestProjectionScrew, & !< matrix of forest projections of screw dislocations for each instance
constitutive_titanmod_TwinforestProjectionEdge, & !< matrix of forest projections of edge dislocations in twin system for each instance
constitutive_titanmod_TwinforestProjectionScrew !< matrix of forest projections of screw dislocations in twin system for each instance
comments for doxygen
2013-07-12 12:27:15 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(:,:,:,:), allocatable, private :: &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_Ctwin66 !< twin elasticity matrix in Mandel notation for each instance
comments for doxygen
2013-07-12 12:27:15 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(:,:,:,:,:,:), allocatable, private :: &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_Ctwin3333 !< twin elasticity matrix for each instance
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
enum, bind(c)
enumerator :: undefined_ID, &
rhoedge_ID, rhoscrew_ID, &
segment_edge_ID, segment_screw_ID, &
resistance_edge_ID, resistance_screw_ID, &
velocity_edge_ID, velocity_screw_ID, &
tau_slip_ID, &
gdot_slip_edge_ID, gdot_slip_screw_ID, &
gdot_slip_ID, &
stressratio_edge_p_ID, stressratio_screw_p_ID, &
shear_system_ID, &
twin_fraction_ID, &
shear_basal_ID, shear_prism_ID, shear_pyra_ID, shear_pyrca_ID, &
rhoedge_basal_ID, rhoedge_prism_ID, rhoedge_pyra_ID, rhoedge_pyrca_ID, &
rhoscrew_basal_ID, rhoscrew_prism_ID, rhoscrew_pyra_ID, rhoscrew_pyrca_ID, &
shear_total_ID
end enum
integer(kind(undefined_ID)), dimension(:,:), allocatable, private :: &
constitutive_titanmod_outputID !< ID of each post result output
comments for doxygen
2013-07-12 12:27:15 +05:30
public :: &
constitutive_titanmod_microstructure, &
constitutive_titanmod_stateInit, &
constitutive_titanmod_init, &
constitutive_titanmod_LpAndItsTangent, &
constitutive_titanmod_dotState, &
added helper functions to communicate accumulated slip between damage and plasticity modules
2014-10-10 22:04:51 +05:30
constitutive_titanmod_getAccumulatedSlip, &
- added get slip rate helper function - added vacancy generation rate - cleaned up some bugs
2014-10-11 15:15:30 +05:30
constitutive_titanmod_getSlipRate, &
comments for doxygen
2013-07-12 12:27:15 +05:30
constitutive_titanmod_postResults, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_homogenizedC
comments for doxygen
2013-07-12 12:27:15 +05:30
contains
!--------------------------------------------------------------------------------------------------
!> @brief module initialization
!> @details reads in material parameters, allocates arrays, and does sanity checks
!--------------------------------------------------------------------------------------------------
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 constitutive_titanmod_init(fileUnit)
comments for doxygen
2013-07-12 12:27:15 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
use debug, only: &
debug_level,&
debug_constitutive,&
debug_levelBasic
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
use math, only: &
math_Mandel3333to66,&
math_Voigt66to3333,&
math_mul3x3
use IO, only: &
IO_read, &
IO_lc, &
IO_getTag, &
IO_isBlank, &
IO_stringPos, &
IO_stringValue, &
IO_floatValue, &
IO_intValue, &
IO_warning, &
IO_error, &
IO_timeStamp, &
IO_EOF
use material, only: &
homogenization_maxNgrains, &
phase_plasticity, &
phase_plasticityInstance, &
phase_Noutput, &
PLASTICITY_TITANMOD_label, &
PLASTICITY_TITANMOD_ID, &
new state related changes
2014-06-11 17:57:41 +05:30
plasticState, &
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
MATERIAL_partPhase
comments for doxygen
2013-07-12 12:27:15 +05:30
use lattice
new state related changes
2014-06-11 17:57:41 +05:30
use numerics,only: &
only output from the root processor for parallel runs
2014-10-10 01:53:06 +05:30
worldrank, &
new state related changes
2014-06-11 17:57:41 +05:30
numerics_integrator
remove omp statements during initialization, "line" for reading in from material.config now has a default (empty) value to prevent conditional jump depend on this value
2013-01-09 03:41:59 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
implicit none
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 22:39:59 +05:30
integer(pInt), intent(in) :: fileUnit
comments for doxygen
2013-07-12 12:27:15 +05:30
remove dotTemperature (returned always 0 anyway) to make it easier to replace it with a heating rate funtction remove potentially dangerous short circuiting in homogenization files added doxygen comment for nonlocal (giving the functions a name) fixed wrong definition of maxnchungs introcuded into dislotwin with last commit, causes trouble with hex
2013-10-09 11:42:16 +05:30
integer(pInt), parameter :: MAXNCHUNKS = LATTICE_maxNinteraction + 1_pInt
comments for doxygen
2013-07-12 12:27:15 +05:30
integer(pInt), dimension(1_pInt+2_pInt*MAXNCHUNKS) :: positions
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt) :: &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
phase, &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
instance, j, k, l, m, n, p, q, r, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
f, o, &
s, s1, s2, &
t, t1, t2, &
ns, nt, &
remove dotTemperature (returned always 0 anyway) to make it easier to replace it with a heating rate funtction remove potentially dangerous short circuiting in homogenization files added doxygen comment for nonlocal (giving the functions a name) fixed wrong definition of maxnchungs introcuded into dislotwin with last commit, causes trouble with hex
2013-10-09 11:42:16 +05:30
Nchunks_SlipSlip, Nchunks_SlipTwin, Nchunks_TwinSlip, Nchunks_TwinTwin, &
Nchunks_SlipFamilies, Nchunks_TwinFamilies, &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
mySize, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
maxTotalNslip,maxTotalNtwin, maxNinstance
new state related changes
2014-06-11 17:57:41 +05:30
integer(pInt) :: sizeState, sizeDotState
integer(pInt) :: NofMyPhase
comments for doxygen
2013-07-12 12:27:15 +05:30
character(len=65536) :: &
tag = '', &
fixed small bug in detecting the lattice structure. only first 3 letters count (e.g. 'ort' and 'orthography' will give the orthorombic lattice type)
2013-11-28 14:26:02 +05:30
line = ''
comments for doxygen
2013-07-12 12:27:15 +05:30
reinstating changes lost in previous commit
2014-10-10 22:15:14 +05:30
mainProcess: if (worldrank == 0) then
write(6,'(/,a)') ' <<<+- constitutive_'//PLASTICITY_TITANMOD_label//' 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"
reinstating changes lost in previous commit
2014-10-10 22:15:14 +05:30
endif mainProcess
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
maxNinstance = int(count(phase_plasticity == PLASTICITY_TITANMOD_ID),pInt)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
if (maxNinstance == 0_pInt) return
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 0_pInt) &
write(6,'(a16,1x,i5,/)') '# instances:',maxNinstance
parsing of material.config file now expects adequate number of (slip/twin) family entries according to given lattice type. used to read up to maxN, but that caused unnecessary IO_warnings...
2013-02-15 03:54:55 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_sizePostResults(maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_sizePostResult(maxval(phase_Noutput),maxNinstance), source=0_pInt)
comments for doxygen
2013-07-12 12:27:15 +05:30
allocate(constitutive_titanmod_output(maxval(phase_Noutput),maxNinstance))
constitutive_titanmod_output = ''
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_outputID(maxval(phase_Noutput),maxNinstance), source=undefined_ID)
allocate(constitutive_titanmod_Noutput(maxNinstance), source=0_pInt)
comments for doxygen
2013-07-12 12:27:15 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_Nslip(lattice_maxNslipFamily,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_Ntwin(lattice_maxNtwinFamily,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_slipFamily(lattice_maxNslip,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_twinFamily(lattice_maxNtwin,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_slipSystemLattice(lattice_maxNslip,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_twinSystemLattice(lattice_maxNtwin,maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_totalNslip(maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_totalNtwin(maxNinstance), source=0_pInt)
allocate(constitutive_titanmod_debyefrequency(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_kinkf0(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CAtomicVolume(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_dc(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinhpconstant(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_GrainSize(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_MaxTwinFraction(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_r(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CEdgeDipMinDistance(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_Cmfptwin(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_Cthresholdtwin(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_aTolRho(maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_rho_edge0(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_rho_screw0(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_shear_system0(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_burgersPerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_burgersPerTwinFam(lattice_maxNtwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_f0_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_tau0e_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_tau0s_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_capre_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_caprs_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_pe_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_ps_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_qe_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_qs_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_v0e_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_v0s_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_kinkcriticallength_PerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinsizePerTwinFam(lattice_maxNtwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CeLambdaSlipPerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CsLambdaSlipPerSlipFam(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
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
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_twinf0_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinshearconstant_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twintau0_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinp_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinq_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twingamma0_PerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinLambdaSlipPerTwinFam(lattice_maxNTwinFamily,maxNinstance), source=0.0_pReal)
comments for doxygen
2013-07-12 12:27:15 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_interactionSlipSlip(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interaction_ee(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interaction_ss(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interaction_es(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionSlipTwin(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionTwinSlip(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionTwinTwin(lattice_maxNinteraction,maxNinstance), source=0.0_pReal)
comments for doxygen
2013-07-12 12:27:15 +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)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
phase = 0_pInt
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>')) /= material_partPhase) ! 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)
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
parsingFile: do while (trim(line) /= IO_EOF) ! read through sections of phase part
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)
comments for doxygen
2013-07-12 12:27:15 +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
comments for doxygen
2013-07-12 12:27:15 +05:30
if (IO_getTag(line,'[',']') /= '') then ! next section
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
phase = phase + 1_pInt ! advance section counter
if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then
Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt)
Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt)
Nchunks_SlipSlip = maxval(lattice_interactionSlipSlip(:,:,phase))
Nchunks_SlipTwin = maxval(lattice_interactionSlipTwin(:,:,phase))
Nchunks_TwinSlip = maxval(lattice_interactionTwinSlip(:,:,phase))
Nchunks_TwinTwin = maxval(lattice_interactionTwinTwin(:,:,phase))
fixed bug in initialization when reading in material.config. Fortran does not short-circuit, e.g. (if a>0 .and. b(a) ==c) might cause an out-of-bounds error when a=-1 it'll try to access b(-1)!
2013-06-12 01:46:40 +05:30
endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
cycle ! skip to next line
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
endif
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
if (phase > 0_pInt ) then; if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then ! one of my sections. Do not short-circuit here (.and. between if-statements), it's not safe in Fortran
instance = phase_plasticityInstance(phase) ! which instance of my plasticity is present phase
positions = IO_stringPos(line,MAXNCHUNKS)
tag = IO_lc(IO_stringValue(line,positions,1_pInt)) ! extract key
select case(tag)
case ('(output)')
select case(IO_lc(IO_stringValue(line,positions,2_pInt)))
case ('rhoedge')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoedge_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoscrew')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoscrew_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('segment_edge')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = segment_edge_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('segment_screw')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = segment_screw_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('resistance_edge')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = resistance_edge_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('resistance_screw')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = resistance_screw_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('velocity_edge')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = velocity_edge_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('velocity_screw')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = velocity_screw_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('tau_slip')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = tau_slip_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('gdot_slip_edge')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = gdot_slip_edge_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('gdot_slip_screw')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = gdot_slip_screw_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('gdot_slip')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = gdot_slip_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('stressratio_edge_p')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = stressratio_edge_p_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('stressratio_screw_p')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = stressratio_screw_p_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_system')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_system_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('twin_fraction')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = twin_fraction_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_basal')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_basal_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_prism')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_prism_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_pyra')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_pyra_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_pyrca')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_pyrca_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoedge_basal')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoedge_basal_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoedge_prism')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoedge_prism_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoedge_pyra')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoedge_pyra_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoedge_pyrca')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoedge_pyrca_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoscrew_basal')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoscrew_basal_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoscrew_prism')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoscrew_prism_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoscrew_pyra')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoscrew_pyra_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('rhoscrew_pyrca')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = rhoscrew_pyrca_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
case ('shear_total')
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_Noutput(instance) = constitutive_titanmod_Noutput(instance) + 1_pInt
fixed error with Noutput counter and few other newstate related changes (nonlocal)
2014-06-26 19:23:12 +05:30
constitutive_titanmod_outputID(constitutive_titanmod_Noutput(instance),instance) = shear_total_ID
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
constitutive_titanmod_output(constitutive_titanmod_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
end select
case ('debyefrequency')
constitutive_titanmod_debyefrequency(instance) = IO_floatValue(line,positions,2_pInt)
case ('kinkf0')
constitutive_titanmod_kinkf0(instance) = IO_floatValue(line,positions,2_pInt)
case ('nslip')
if (positions(1) < 1_pInt + Nchunks_SlipFamilies) &
call IO_warning(50_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_Nslip(j,instance) = IO_intValue(line,positions,1_pInt+j)
enddo
case ('ntwin')
if (positions(1) < 1_pInt + Nchunks_TwinFamilies) &
call IO_warning(51_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_Ntwin(j,instance) = IO_intValue(line,positions,1_pInt+j)
enddo
case ('rho_edge0')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_rho_edge0(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('rho_screw0')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_rho_screw0(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('slipburgers')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_burgersPerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinburgers')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_burgersPerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('f0')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_f0_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinf0')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinf0_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('tau0e')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_tau0e_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twintau0')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twintau0_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('tau0s')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_tau0s_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('capre')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_capre_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('caprs')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_caprs_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('v0e')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_v0e_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twingamma0')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twingamma0_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('v0s')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_v0s_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('kinkcriticallength')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_kinkcriticallength_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinsize')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinsizePerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('celambdaslip')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_CeLambdaSlipPerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinlambdaslip')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinlambdaslipPerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('cslambdaslip')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_CsLambdaSlipPerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('grainsize')
constitutive_titanmod_GrainSize(instance) = IO_floatValue(line,positions,2_pInt)
case ('maxtwinfraction')
constitutive_titanmod_MaxTwinFraction(instance) = IO_floatValue(line,positions,2_pInt)
case ('pe')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_pe_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinp')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinp_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('ps')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_ps_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('qe')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_qe_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinq')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinq_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('qs')
do j = 1_pInt, Nchunks_SlipFamilies
constitutive_titanmod_qs_PerSlipFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('twinshearconstant')
do j = 1_pInt, Nchunks_TwinFamilies
constitutive_titanmod_twinshearconstant_PerTwinFam(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('dc')
constitutive_titanmod_dc(instance) = IO_floatValue(line,positions,2_pInt)
case ('twinhpconstant')
constitutive_titanmod_twinhpconstant(instance) = IO_floatValue(line,positions,2_pInt)
case ('atol_rho')
constitutive_titanmod_aTolRho(instance) = IO_floatValue(line,positions,2_pInt)
case ('interactionee')
do j = 1_pInt, lattice_maxNinteraction
constitutive_titanmod_interaction_ee(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interactionss')
do j = 1_pInt, lattice_maxNinteraction
constitutive_titanmod_interaction_ss(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interactiones')
do j = 1_pInt, lattice_maxNinteraction
constitutive_titanmod_interaction_es(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interaction_slipslip','interactionslipslip')
if (positions(1) < 1_pInt + Nchunks_SlipSlip) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_SlipSlip
constitutive_titanmod_interactionSlipSlip(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interaction_sliptwin','interactionsliptwin')
if (positions(1) < 1_pInt + Nchunks_SlipTwin) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_SlipTwin
constitutive_titanmod_interactionSlipTwin(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interaction_twinslip','interactiontwinslip')
if (positions(1) < 1_pInt + Nchunks_TwinSlip) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_TwinSlip
constitutive_titanmod_interactionTwinSlip(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
case ('interaction_twintwin','interactiontwintwin')
if (positions(1) < 1_pInt + Nchunks_TwinTwin) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_TITANMOD_label//')')
do j = 1_pInt, Nchunks_TwinTwin
constitutive_titanmod_interactionTwinTwin(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
end select
endif; endif
enddo parsingFile
sanityChecks: do phase = 1_pInt, size(phase_plasticity)
myPhase: if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then
instance = phase_plasticityInstance(phase)
added kmc example material.config added elastic test for mkc, titanmod, and J2
2014-08-08 19:21:11 +05:30
if (sum(constitutive_titanmod_Nslip(:,instance)) < 0_pInt) &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
call IO_error(211_pInt,el=instance,ext_msg='nslip ('//PLASTICITY_TITANMOD_label//')')
if (sum(constitutive_titanmod_Ntwin(:,instance)) < 0_pInt) &
call IO_error(211_pInt,el=instance,ext_msg='ntwin ('//PLASTICITY_TITANMOD_label//')')
do f = 1_pInt,lattice_maxNslipFamily
if (constitutive_titanmod_Nslip(f,instance) > 0_pInt) then
if (constitutive_titanmod_rho_edge0(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='rho_edge0 ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_rho_screw0(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='rho_screw0 ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_burgersPerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='slipburgers ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_f0_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='f0 ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_tau0e_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='tau0e ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_tau0s_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='tau0s ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_capre_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='capre ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_caprs_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='caprs ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_v0e_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='v0e ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_v0s_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='v0s ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_kinkcriticallength_PerSlipFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='kinkCriticalLength ('//PLASTICITY_TITANMOD_label//')')
endif
enddo
do f = 1_pInt,lattice_maxNtwinFamily
if (constitutive_titanmod_Ntwin(f,instance) > 0_pInt) then
if (constitutive_titanmod_burgersPerTwinFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinburgers ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_twinf0_PerTwinFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinf0 ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_twinshearconstant_PerTwinFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinshearconstant ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_twintau0_PerTwinFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twintau0 ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_twingamma0_PerTwinFam(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twingamma0 ('//PLASTICITY_TITANMOD_label//')')
endif
enddo
if (constitutive_titanmod_dc(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='dc ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_twinhpconstant(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinhpconstant ('//PLASTICITY_TITANMOD_label//')')
if (constitutive_titanmod_aTolRho(instance) <= 0.0_pReal) &
added kmc example material.config added elastic test for mkc, titanmod, and J2
2014-08-08 19:21:11 +05:30
call IO_error(211_pInt,el=instance,ext_msg='aTol_rho ('//PLASTICITY_TITANMOD_label//')')
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! determine total number of active slip or twin systems
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_Nslip(:,instance) = min(lattice_NslipSystem(:,phase),constitutive_titanmod_Nslip(:,instance))
constitutive_titanmod_Ntwin(:,instance) = min(lattice_NtwinSystem(:,phase),constitutive_titanmod_Ntwin(:,instance))
constitutive_titanmod_totalNslip(instance) = sum(constitutive_titanmod_Nslip(:,instance))
constitutive_titanmod_totalNtwin(instance) = sum(constitutive_titanmod_Ntwin(:,instance))
endif myPhase
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo sanityChecks
!--------------------------------------------------------------------------------------------------
! allocation of variables whose size depends on the total number of active slip systems
maxTotalNslip = maxval(constitutive_titanmod_totalNslip)
maxTotalNtwin = maxval(constitutive_titanmod_totalNtwin)
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_burgersPerSlipSys(maxTotalNslip, maxNinstance), source=0.0_pReal)
comments for doxygen
2013-07-12 12:27:15 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_f0_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_tau0e_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_tau0s_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_capre_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_caprs_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_pe_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_ps_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_qe_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_qs_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_v0e_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_v0s_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_kinkcriticallength_PerSlipSys(maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CeLambdaSlipPerSlipSys(maxTotalNslip, maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_CsLambdaSlipPerSlipSys(maxTotalNslip, maxNinstance), source=0.0_pReal)
comments for doxygen
2013-07-12 12:27:15 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_burgersPerTwinSys(maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinf0_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinshearconstant_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twintau0_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinp_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinq_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twingamma0_PerTwinSys(maxTotalNTwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinsizePerTwinSys(maxTotalNtwin, maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_twinLambdaSlipPerTwinSys(maxTotalNtwin, maxNinstance), source=0.0_pReal)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
allocate(constitutive_titanmod_Ctwin66 (6,6,maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_Ctwin3333 (3,3,3,3,maxTotalNtwin,maxNinstance), source=0.0_pReal)
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
allocate(constitutive_titanmod_interactionMatrixSlipSlip(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrix_ee(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrix_ss(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrix_es(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrixSlipTwin(maxTotalNslip,maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrixTwinSlip(maxTotalNtwin,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_interactionMatrixTwinTwin(maxTotalNtwin,maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_forestProjectionEdge(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_forestProjectionScrew(maxTotalNslip,maxTotalNslip,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_TwinforestProjectionEdge(maxTotalNtwin,maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_titanmod_TwinforestProjectionScrew(maxTotalNtwin,maxTotalNtwin,maxNinstance), source=0.0_pReal)
comments for doxygen
2013-07-12 12:27:15 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
initializeInstances: do phase = 1_pInt, size(phase_plasticity)
if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then
instance = phase_plasticityInstance(phase)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! inverse lookup of slip system family
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
l = 0_pInt
do f = 1_pInt,lattice_maxNslipFamily
do s = 1_pInt,constitutive_titanmod_Nslip(f,instance)
l = l + 1_pInt
constitutive_titanmod_slipFamily(l,instance) = f
constitutive_titanmod_slipSystemLattice(l,instance) = sum(lattice_NslipSystem(1:f-1_pInt,phase)) + s
enddo; enddo
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! inverse lookup of twin system family
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
l = 0_pInt
do f = 1_pInt,lattice_maxNtwinFamily
do t = 1_pInt,constitutive_titanmod_Ntwin(f,instance)
l = l + 1_pInt
constitutive_titanmod_twinFamily(l,instance) = f
constitutive_titanmod_twinSystemLattice(l,instance) = sum(lattice_NtwinSystem(1:f-1_pInt,phase)) + t
enddo; enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! determine size of state array
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
new state related changes
2014-06-11 17:57:41 +05:30
sizeDotState = &
size(constitutive_titanmod_listBasicSlipStates)*ns + &
size(constitutive_titanmod_listBasicTwinStates)*nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sizeState = sizeDotState+ &
new state related changes
2014-06-11 17:57:41 +05:30
size(constitutive_titanmod_listDependentSlipStates)*ns + &
size(constitutive_titanmod_listDependentTwinStates)*nt
parsing of material.config file now expects adequate number of (slip/twin) family entries according to given lattice type. used to read up to maxN, but that caused unnecessary IO_warnings...
2013-02-15 03:54:55 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! determine size of postResults array
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
outputsLoop: do o = 1_pInt,constitutive_titanmod_Noutput(instance)
mySize = 0_pInt
select case(constitutive_titanmod_outputID(o,instance))
case(rhoedge_ID, rhoscrew_ID, &
segment_edge_ID, segment_screw_ID, &
resistance_edge_ID, resistance_screw_ID, &
velocity_edge_ID, velocity_screw_ID, &
tau_slip_ID, &
gdot_slip_edge_ID, gdot_slip_screw_ID, &
gdot_slip_ID, &
stressratio_edge_p_ID, stressratio_screw_p_ID, &
shear_system_ID)
mySize = constitutive_titanmod_totalNslip(instance)
case(twin_fraction_ID)
mySize = constitutive_titanmod_totalNtwin(instance)
case(shear_basal_ID, shear_prism_ID, shear_pyra_ID, shear_pyrca_ID, & ! use only if all 4 slip families in hex are considered
rhoedge_basal_ID, rhoedge_prism_ID, rhoedge_pyra_ID, rhoedge_pyrca_ID, &
rhoscrew_basal_ID, rhoscrew_prism_ID, rhoscrew_pyra_ID, rhoscrew_pyrca_ID, &
shear_total_ID)
mySize = 1_pInt
case default
call IO_error(105_pInt,ext_msg=constitutive_titanmod_output(o,instance)// &
' ('//PLASTICITY_TITANMOD_label//')')
end select
outputFound: if (mySize > 0_pInt) then
constitutive_titanmod_sizePostResult(o,instance) = mySize
constitutive_titanmod_sizePostResults(instance) = constitutive_titanmod_sizePostResults(instance) + mySize
endif outputFound
enddo outputsLoop
new state related changes
2014-06-11 17:57:41 +05:30
! Determine size of state array
plasticState(phase)%sizeState = sizeState
plasticState(phase)%sizeDotState = sizeDotState
no need to allocate a large array for constitutive(thermal,damage)_sizePostResults (changes apply only to new state)
2014-06-30 20:17:30 +05:30
plasticState(phase)%sizePostResults = constitutive_titanmod_sizePostResults(instance)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
allocate(plasticState(phase)%aTolState (sizeState), source=constitutive_titanmod_aTolRho(instance))
new state related changes
2014-06-11 17:57:41 +05:30
allocate(plasticState(phase)%state0 (sizeState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%partionedState0 (sizeState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%subState0 (sizeState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%state (sizeState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%state_backup (sizeState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%dotState (sizeDotState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%deltaState (sizeDotState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%dotState_backup (sizeDotState,NofMyPhase), source=0.0_pReal)
if (any(numerics_integrator == 1_pInt)) then
allocate(plasticState(phase)%previousDotState (sizeDotState,NofMyPhase), source=0.0_pReal)
allocate(plasticState(phase)%previousDotState2 (sizeDotState,NofMyPhase), source=0.0_pReal)
endif
if (any(numerics_integrator == 4_pInt)) &
allocate(plasticState(phase)%RK4dotState (sizeDotState,NofMyPhase), source=0.0_pReal)
if (any(numerics_integrator == 5_pInt)) &
allocate(plasticState(phase)%RKCK45dotState (6,sizeDotState,NofMyPhase),source=0.0_pReal)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! construction of the twin elasticity matrices
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
do j=1_pInt,lattice_maxNtwinFamily
do k=1_pInt,constitutive_titanmod_Ntwin(j,instance)
do l=1_pInt,3_pInt ; do m=1_pInt,3_pInt ; do n=1_pInt,3_pInt ; do o=1_pInt,3_pInt
do p=1_pInt,3_pInt ; do q=1_pInt,3_pInt ; do r=1_pInt,3_pInt ; do s=1_pInt,3_pInt
constitutive_titanmod_Ctwin3333(l,m,n,o,sum(constitutive_titanmod_Nslip(1:j-1_pInt,instance))+k,instance) = &
constitutive_titanmod_Ctwin3333(l,m,n,o,sum(constitutive_titanmod_Nslip(1:j-1_pInt,instance))+k,instance) + &
lattice_C3333(p,q,r,s,phase)*&
lattice_Qtwin(l,p,sum(lattice_NslipSystem(1:j-1_pInt,phase))+k,phase)* &
lattice_Qtwin(m,q,sum(lattice_NslipSystem(1:j-1_pInt,phase))+k,phase)* &
lattice_Qtwin(n,r,sum(lattice_NslipSystem(1:j-1_pInt,phase))+k,phase)* &
lattice_Qtwin(o,s,sum(lattice_NslipSystem(1:j-1_pInt,phase))+k,phase)
enddo; enddo; enddo; enddo
enddo; enddo; enddo ; enddo
constitutive_titanmod_Ctwin66(1:6,1:6,k,instance) = &
math_Mandel3333to66(constitutive_titanmod_Ctwin3333(1:3,1:3,1:3,1:3,k,instance))
enddo; enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! Burgers vector, dislocation velocity prefactor for each slip system
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
do s = 1_pInt,constitutive_titanmod_totalNslip(instance)
f = constitutive_titanmod_slipFamily(s,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_burgersPerSlipSys(s,instance) = &
constitutive_titanmod_burgersPerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_f0_PerSlipSys(s,instance) = &
constitutive_titanmod_f0_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_tau0e_PerSlipSys(s,instance) = &
constitutive_titanmod_tau0e_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_tau0s_PerSlipSys(s,instance) = &
constitutive_titanmod_tau0s_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_capre_PerSlipSys(s,instance) = &
constitutive_titanmod_capre_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_caprs_PerSlipSys(s,instance) = &
constitutive_titanmod_caprs_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_v0e_PerSlipSys(s,instance) = &
constitutive_titanmod_v0e_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_v0s_PerSlipSys(s,instance) = &
constitutive_titanmod_v0s_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_kinkcriticallength_PerSlipSys(s,instance) = &
constitutive_titanmod_kinkcriticallength_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_pe_PerSlipSys(s,instance) = &
constitutive_titanmod_pe_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_ps_PerSlipSys(s,instance) = &
constitutive_titanmod_ps_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_qe_PerSlipSys(s,instance) = &
constitutive_titanmod_qe_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_qs_PerSlipSys(s,instance) = &
constitutive_titanmod_qs_PerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_CeLambdaSlipPerSlipSys(s,instance) = &
constitutive_titanmod_CeLambdaSlipPerSlipFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_CsLambdaSlipPerSlipSys(s,instance) = &
constitutive_titanmod_CsLambdaSlipPerSlipFam(f,instance)
enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! Burgers vector, nucleation rate prefactor and twin size for each twin system
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
do t = 1_pInt,constitutive_titanmod_totalNtwin(instance)
f = constitutive_titanmod_twinFamily(t,instance)
constitutive_titanmod_burgersPerTwinSys(t,instance) = &
constitutive_titanmod_burgersPerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinsizePerTwinSys(t,instance) = &
constitutive_titanmod_twinsizePerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinf0_PerTwinSys(t,instance) = &
constitutive_titanmod_twinf0_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(t,instance) = &
constitutive_titanmod_twinshearconstant_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twintau0_PerTwinSys(t,instance) = &
constitutive_titanmod_twintau0_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twingamma0_PerTwinSys(t,instance) = &
constitutive_titanmod_twingamma0_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinp_PerTwinSys(t,instance) = &
constitutive_titanmod_twinp_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinq_PerTwinSys(t,instance) = &
constitutive_titanmod_twinq_PerTwinFam(f,instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinLambdaSlipPerTwinSys(t,instance) = &
constitutive_titanmod_twinLambdaSlipPerTwinFam(f,instance)
enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! Construction of interaction matrices
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
do s1 = 1_pInt,constitutive_titanmod_totalNslip(instance)
do s2 = 1_pInt,constitutive_titanmod_totalNslip(instance)
constitutive_titanmod_interactionMatrixSlipSlip(s1,s2,instance) = &
constitutive_titanmod_interactionSlipSlip(lattice_interactionSlipSlip( &
constitutive_titanmod_slipSystemLattice(s1,instance),&
constitutive_titanmod_slipSystemLattice(s2,instance),phase),instance)
constitutive_titanmod_interactionMatrix_ee(s1,s2,instance) = &
constitutive_titanmod_interaction_ee(lattice_interactionSlipSlip ( &
constitutive_titanmod_slipSystemLattice(s1,instance), &
constitutive_titanmod_slipSystemLattice(s2,instance), phase),instance)
constitutive_titanmod_interactionMatrix_ss(s1,s2,instance) = &
constitutive_titanmod_interaction_ss(lattice_interactionSlipSlip( &
constitutive_titanmod_slipSystemLattice(s1,instance), &
constitutive_titanmod_slipSystemLattice(s2,instance), phase),instance)
constitutive_titanmod_interactionMatrix_es(s1,s2,instance) = &
constitutive_titanmod_interaction_es(lattice_interactionSlipSlip( &
constitutive_titanmod_slipSystemLattice(s1,instance), &
constitutive_titanmod_slipSystemLattice(s2,instance), phase),instance)
enddo; enddo
do s1 = 1_pInt,constitutive_titanmod_totalNslip(instance)
do t2 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
constitutive_titanmod_interactionMatrixSlipTwin(s1,t2,instance) = &
constitutive_titanmod_interactionSlipTwin(lattice_interactionSlipTwin( &
constitutive_titanmod_slipSystemLattice(s1,instance), &
constitutive_titanmod_twinSystemLattice(t2,instance), phase),instance)
enddo; enddo
do t1 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
do s2 = 1_pInt,constitutive_titanmod_totalNslip(instance)
constitutive_titanmod_interactionMatrixTwinSlip(t1,s2,instance) = &
constitutive_titanmod_interactionTwinSlip(lattice_interactionTwinSlip( &
constitutive_titanmod_twinSystemLattice(t1,instance), &
constitutive_titanmod_slipSystemLattice(s2,instance), phase),instance)
enddo; enddo
do t1 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
do t2 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
constitutive_titanmod_interactionMatrixTwinTwin(t1,t2,instance) = &
constitutive_titanmod_interactionTwinTwin(lattice_interactionTwinTwin( &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinSystemLattice(t1,instance), &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_twinSystemLattice(t2,instance), phase),instance)
enddo; enddo
do s1 = 1_pInt,constitutive_titanmod_totalNslip(instance)
do s2 = 1_pInt,constitutive_titanmod_totalNslip(instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! calculation of forest projections for edge dislocations
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_forestProjectionEdge(s1,s2,instance) = &
abs(math_mul3x3(lattice_sn(:,constitutive_titanmod_slipSystemLattice(s1,instance),phase), &
lattice_st(:,constitutive_titanmod_slipSystemLattice(s2,instance),phase)))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! calculation of forest projections for screw dislocations
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_forestProjectionScrew(s1,s2,instance) = &
abs(math_mul3x3(lattice_sn(:,constitutive_titanmod_slipSystemLattice(s1,instance),phase), &
lattice_sd(:,constitutive_titanmod_slipSystemLattice(s2,instance),phase)))
enddo; enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! calculation of forest projections for edge dislocations in twin system
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
do t1 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
do t2 = 1_pInt,constitutive_titanmod_totalNtwin(instance)
constitutive_titanmod_TwinforestProjectionEdge(t1,t2,instance) = &
abs(math_mul3x3(lattice_tn(:,constitutive_titanmod_twinSystemLattice(t1,instance),phase), &
lattice_tt(:,constitutive_titanmod_twinSystemLattice(t2,instance),phase)))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! calculation of forest projections for screw dislocations in twin system
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_TwinforestProjectionScrew(t1,t2,instance) = &
abs(math_mul3x3(lattice_tn(:,constitutive_titanmod_twinSystemLattice(t1,instance),phase), &
lattice_td(:,constitutive_titanmod_twinSystemLattice(t2,instance),phase)))
enddo; enddo
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
call constitutive_titanmod_stateInit(phase,instance)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
endif
enddo initializeInstances
comments for doxygen
2013-07-12 12:27:15 +05:30
end subroutine constitutive_titanmod_init
!--------------------------------------------------------------------------------------------------
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!> @brief sets the initial microstructural state for a given instance of this plasticity
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_titanmod_stateInit(ph,instance)
new state related changes
2014-06-11 17:57:41 +05:30
use lattice, only: &
lattice_maxNslipFamily, &
lattice_maxNtwinFamily, &
lattice_mu
use material, only: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState, &
mappingConstitutive
new state related changes
2014-06-11 17:57:41 +05:30
implicit none
integer(pInt), intent(in) :: instance !< number specifying the instance of the plasticity
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
integer(pInt), intent(in) :: ph !< number specifying the phase of the plasticity
new state related changes
2014-06-11 17:57:41 +05:30
integer(pInt) :: &
s,s0,s1, &
t,t0,t1, &
ns,nt,f
real(pReal), dimension(constitutive_titanmod_totalNslip(instance)) :: &
rho_edge0, &
rho_screw0, &
shear_system0, &
segment_edge0, &
segment_screw0, &
resistance_edge0, &
resistance_screw0
real(pReal), dimension(constitutive_titanmod_totalNtwin(instance)) :: &
twingamma_dot0, &
resistance_twin0
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
real(pReal), dimension(plasticState(ph)%sizeState) :: tempState !!!!!!!!!????????? check
new state related changes
2014-06-11 17:57:41 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
tempState = 0.0_pReal
!--------------------------------------------------------------------------------------------------
! initialize basic slip state variables for slip
s1 = 0_pInt
do f = 1_pInt,lattice_maxNslipFamily
s0 = s1 + 1_pInt
s1 = s0 + constitutive_titanmod_Nslip(f,instance) - 1_pInt
do s = s0,s1
rho_edge0(s) = constitutive_titanmod_rho_edge0(f,instance)
rho_screw0(s) = constitutive_titanmod_rho_screw0(f,instance)
shear_system0(s) = 0.0_pReal
enddo
enddo
!--------------------------------------------------------------------------------------------------
! initialize basic slip state variables for twin
t1 = 0_pInt
do f = 1_pInt,lattice_maxNtwinFamily
t0 = t1 + 1_pInt
t1 = t0 + constitutive_titanmod_Ntwin(f,instance) - 1_pInt
do t = t0,t1
twingamma_dot0(t)=0.0_pReal
enddo
enddo
!--------------------------------------------------------------------------------------------------
! initialize dependent slip microstructural variables
forall (s = 1_pInt:ns)
segment_edge0(s) = constitutive_titanmod_CeLambdaSlipPerSlipSys(s,instance)/ &
sqrt(dot_product((rho_edge0),constitutive_titanmod_forestProjectionEdge(1:ns,s,instance))+ &
dot_product((rho_screw0),constitutive_titanmod_forestProjectionScrew(1:ns,s,instance)))
segment_screw0(s) = constitutive_titanmod_CsLambdaSlipPerSlipSys(s,instance)/ &
sqrt(dot_product((rho_edge0),constitutive_titanmod_forestProjectionEdge(1:ns,s,instance))+ &
dot_product((rho_screw0),constitutive_titanmod_forestProjectionScrew(1:ns,s,instance)))
resistance_edge0(s) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_mu(ph)*constitutive_titanmod_burgersPerSlipSys(s,instance)* &
new state related changes
2014-06-11 17:57:41 +05:30
sqrt(dot_product((rho_edge0),constitutive_titanmod_interactionMatrix_ee(1:ns,s,instance))+ &
dot_product((rho_screw0),constitutive_titanmod_interactionMatrix_es(1:ns,s,instance)))
resistance_screw0(s) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_mu(ph)*constitutive_titanmod_burgersPerSlipSys(s,instance)* &
new state related changes
2014-06-11 17:57:41 +05:30
sqrt(dot_product((rho_edge0),constitutive_titanmod_interactionMatrix_es(1:ns,s,instance))+ &
dot_product((rho_screw0), constitutive_titanmod_interactionMatrix_ss(1:ns,s,instance)))
end forall
forall (t = 1_pInt:nt) &
resistance_twin0(t) = 0.0_pReal
tempState = 0.0_pReal
tempState (1:ns) = rho_edge0
tempState (1_pInt*ns+1_pInt:2_pInt*ns) = rho_screw0
tempState (2_pInt*ns+1_pInt:3_pInt*ns) = shear_system0
tempState (3_pInt*ns+1_pInt:3_pInt*ns+nt) = twingamma_dot0
tempState (3_pInt*ns+nt+1_pInt:4_pInt*ns+nt) = segment_edge0
tempState (4_pInt*ns+nt+1_pInt:5_pInt*ns+nt) = segment_screw0
tempState (5_pInt*ns+nt+1_pInt:6_pInt*ns+nt) = resistance_edge0
tempState (6_pInt*ns+nt+1_pInt:7_pInt*ns+nt) = resistance_screw0
tempState (7_pInt*ns+nt+1_pInt:7_pInt*ns+2_pInt*nt)=resistance_twin0
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state0 = spread(tempState,2,size(plasticState(ph)%state(1,:)))
end subroutine constitutive_titanmod_stateInit
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns the homogenized elasticity matrix
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
function constitutive_titanmod_homogenizedC(ipc,ip,el)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_plasticityInstance, &
plasticState, &
mappingConstitutive
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
use lattice, only: &
lattice_C66
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal), dimension(6,6) :: &
constitutive_titanmod_homogenizedC
integer(pInt), intent(in) :: &
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
comments for doxygen
2013-07-12 12:27:15 +05:30
real(pReal), dimension(constitutive_titanmod_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
volumefraction_PerTwinSys
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
ph, &
of, &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
instance, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ns, nt, &
i
real(pReal) :: &
sumf
new state related changes
2014-06-11 17:57:41 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! shortened notation
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! ph = material_phase(ipc,ip,el)
of = mappingConstitutive(1,ipc,ip,el)
ph = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(ph)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! total twin volume fraction
comments for doxygen
2013-07-12 12:27:15 +05:30
do i=1_pInt,nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
volumefraction_PerTwinSys(i)=plasticState(ph)%state(3_pInt*ns+i, of)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(i,instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
sumf = sum(abs(volumefraction_PerTwinSys(1:nt))) ! safe for nt == 0
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! homogenized elasticity matrix
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_homogenizedC = (1.0_pReal-sumf)*lattice_C66(1:6,1:6,ph)
comments for doxygen
2013-07-12 12:27:15 +05:30
do i=1_pInt,nt
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
constitutive_titanmod_homogenizedC = constitutive_titanmod_homogenizedC &
+ volumefraction_PerTwinSys(i)*&
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
constitutive_titanmod_Ctwin66(1:6,1:6,i,instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
end function constitutive_titanmod_homogenizedC
Brushed up accountability of twinning to Lp
2010-09-07 20:14:37 +05:30
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief calculates derived quantities from state
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_titanmod_microstructure(temperature,ipc,ip,el)
comments for doxygen
2013-07-12 12:27:15 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase,&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_plasticityInstance, &
plasticState, &
mappingConstitutive
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
use lattice, only: &
lattice_mu
comments for doxygen
2013-07-12 12:27:15 +05:30
implicit none
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt), intent(in) :: &
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
real(pReal), intent(in) :: &
temperature !< temperature at IP
integer(pInt) :: &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
instance, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ns, nt, s, t, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
i, &
ph, &
of
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal) :: &
sumf, &
sfe ! stacking fault energy
comments for doxygen
2013-07-12 12:27:15 +05:30
real(pReal), dimension(constitutive_titanmod_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
volumefraction_PerTwinSys
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
new state related changes
2014-06-11 17:57:41 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!Shortened notation
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
of = mappingConstitutive(1,ipc,ip,el)
ph = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(ph)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! total twin volume fraction
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
forall (i = 1_pInt:nt) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
volumefraction_PerTwinSys(i)=plasticState(ph)%state(3_pInt*ns+i, of)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(i,instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
sumf = sum(abs(volumefraction_PerTwinSys(1:nt))) ! safe for nt == 0
sfe = 0.0002_pReal*Temperature-0.0396_pReal
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
! average segment length for edge dislocations in matrix
comments for doxygen
2013-07-12 12:27:15 +05:30
forall (s = 1_pInt:ns) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(3_pInt*ns+nt+s, of) = constitutive_titanmod_CeLambdaSlipPerSlipSys(s,instance)/ &
sqrt(dot_product(plasticState(ph)%state(1:ns, of), &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_forestProjectionEdge(1:ns,s,instance))+ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
dot_product(plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of), &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_forestProjectionScrew(1:ns,s,instance)))
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
! average segment length for screw dislocations in matrix
comments for doxygen
2013-07-12 12:27:15 +05:30
forall (s = 1_pInt:ns) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(4_pInt*ns+nt+s, of) = constitutive_titanmod_CsLambdaSlipPerSlipSys(s,instance)/ &
sqrt(dot_product(plasticState(ph)%state(1:ns, of), &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_forestProjectionEdge(1:ns,s,instance))+ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
dot_product(plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of), &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_forestProjectionScrew(1:ns,s,instance)))
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
! threshold stress or slip resistance for edge dislocation motion
forall (s = 1_pInt:ns) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(5_pInt*ns+nt+s, of) = &
lattice_mu(ph)*constitutive_titanmod_burgersPerSlipSys(s,instance)*&
sqrt(dot_product((plasticState(ph)%state(1:ns, of)),&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_interactionMatrix_ee(1:ns,s,instance))+ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
dot_product((plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of)),&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_interactionMatrix_es(1:ns,s,instance)))
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
! threshold stress or slip resistance for screw dislocation motion
forall (s = 1_pInt:ns) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(6_pInt*ns+nt+s, of) = &
lattice_mu(ph)*constitutive_titanmod_burgersPerSlipSys(s,instance)*&
sqrt(dot_product((plasticState(ph)%state(1:ns, of)),&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_interactionMatrix_es(1:ns,s,instance))+ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
dot_product((plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of)),&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_interactionMatrix_ss(1:ns,s,instance)))
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
! threshold stress or slip resistance for dislocation motion in twin
forall (t = 1_pInt:nt) &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(7_pInt*ns+nt+t, of) = &
lattice_mu(ph)*constitutive_titanmod_burgersPerTwinSys(t,instance)*&
(dot_product((abs(plasticState(ph)%state(2_pInt*ns+1_pInt:2_pInt*ns+nt, of))),&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_interactionMatrixTwinTwin(1:nt,t,instance)))
new state related changes
2014-06-11 17:57:41 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! state=tempState
comments for doxygen
2013-07-12 12:27:15 +05:30
end subroutine constitutive_titanmod_microstructure
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief calculates plastic velocity gradient and its tangent
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_titanmod_LpAndItsTangent(Lp,dLp_dTstar99,Tstar_v,temperature,ipc,ip,el)
comments for doxygen
2013-07-12 12:27:15 +05:30
use math, only: &
math_Plain3333to99, &
math_Mandel6to33
use lattice, only: &
lattice_Sslip, &
lattice_Sslip_v, &
lattice_Stwin, &
lattice_Stwin_v, &
lattice_maxNslipFamily, &
lattice_maxNtwinFamily, &
lattice_NslipSystem, &
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
lattice_NtwinSystem, &
lattice_structure, &
LATTICE_hex_ID
comments for doxygen
2013-07-12 12:27:15 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_plasticityInstance, &
plasticState, &
mappingConstitutive
comments for doxygen
2013-07-12 12:27:15 +05:30
implicit none
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), dimension(3,3), intent(out) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
Lp !< plastic velocity gradient
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), dimension(9,9), intent(out) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
dLp_dTstar99 !< derivative of Lp with respect to 2nd Piola Kirchhoff stress
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), dimension(6), intent(in) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
Tstar_v !< 2nd Piola Kirchhoff stress tensor in Mandel notation
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), intent(in) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
temperature !< temperature at IP
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
integer(pInt), intent(in) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily, instance, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ns,nt, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
f,i,j,k,l,m,n, &
ph, &
of
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal) :: sumf, &
StressRatio_edge_p, minusStressRatio_edge_p, StressRatio_edge_pminus1, BoltzmannRatioedge, &
StressRatio_screw_p, minusStressRatio_screw_p, StressRatio_screw_pminus1, BoltzmannRatioscrew, &
twinStressRatio_p, twinminusStressRatio_p, twinStressRatio_pminus1, BoltzmannRatiotwin, &
twinDotGamma0, bottomstress_edge, bottomstress_screw, screwvelocity_prefactor
comments for doxygen
2013-07-12 12:27:15 +05:30
real(pReal), dimension(3,3,3,3) :: dLp_dTstar3333
real(pReal), dimension(constitutive_titanmod_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
gdot_slip,dgdot_dtauslip,tau_slip, &
edge_velocity, screw_velocity, &
gdot_slip_edge, gdot_slip_screw
comments for doxygen
2013-07-12 12:27:15 +05:30
real(pReal), dimension(constitutive_titanmod_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_twin,dgdot_dtautwin,tau_twin, volumefraction_PerTwinSys
new state related changes
2014-06-11 17:57:41 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! tempState=state
new state related changes
2014-06-11 17:57:41 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! shortened notation
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
of = mappingConstitutive(1,ipc,ip,el)
ph = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(ph)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
do i=1_pInt,nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
volumefraction_PerTwinSys(i)=plasticState(ph)%state(3_pInt*ns+i, of)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(i,instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
sumf = sum(abs(volumefraction_PerTwinSys(1:nt))) ! safe for nt == 0
Lp = 0.0_pReal
dLp_dTstar3333 = 0.0_pReal
dLp_dTstar99 = 0.0_pReal
!* Dislocation glide part
gdot_slip = 0.0_pReal
gdot_slip_edge = 0.0_pReal
gdot_slip_screw = 0.0_pReal
dgdot_dtauslip = 0.0_pReal
j = 0_pInt
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
slipFamiliesLoop: do f = 1_pInt,lattice_maxNslipFamily
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
do i = 1_pInt,constitutive_titanmod_Nslip(f,instance) ! process each (active) slip system in family
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
j = j+1_pInt
!* Calculation of Lp
!* Resolved shear stress on slip system
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
tau_slip(j) = dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))
if(lattice_structure(ph)==LATTICE_hex_ID) then ! only for prismatic and pyr <a> systems in hex
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
screwvelocity_prefactor=constitutive_titanmod_debyefrequency(instance)* &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(4_pInt*ns+nt+j, of)*(constitutive_titanmod_burgersPerSlipSys(j,instance)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_kinkcriticallength_PerSlipSys(j,instance))**2
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
!* Stress ratio for screw ! No slip resistance for screw dislocations, only Peierls stress
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
bottomstress_screw=constitutive_titanmod_tau0s_PerSlipSys(j,instance)
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
StressRatio_screw_p = ((abs(tau_slip(j)))/ &
Final version of titanmod
2011-03-31 14:51:43 +05:30
( bottomstress_screw) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**constitutive_titanmod_ps_PerSlipSys(j,instance)
Final version of titanmod
2011-03-31 14:51:43 +05:30
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
if((1.0_pReal-StressRatio_screw_p)>0.001_pReal) then
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
minusStressRatio_screw_p=1.0_pReal-StressRatio_screw_p
else
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
minusStressRatio_screw_p=0.001_pReal
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
bottomstress_screw=constitutive_titanmod_tau0s_PerSlipSys(j,instance)
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
StressRatio_screw_pminus1 = ((abs(tau_slip(j)))/ &
Final version of titanmod
2011-03-31 14:51:43 +05:30
( bottomstress_screw) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**(constitutive_titanmod_ps_PerSlipSys(j,instance)-1.0_pReal)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
!* Boltzmann ratio for screw
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatioscrew = constitutive_titanmod_kinkf0(instance)/(kB*Temperature)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
Final version of titanmod
2011-03-31 14:51:43 +05:30
else ! if the structure is not hex or the slip family is basal
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
screwvelocity_prefactor=constitutive_titanmod_v0s_PerSlipSys(j,instance)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
bottomstress_screw=constitutive_titanmod_tau0s_PerSlipSys(j,instance)+ &
plasticState(ph)%state(6*ns+nt+j, of)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
StressRatio_screw_p = ((abs(tau_slip(j)))/( bottomstress_screw ))**constitutive_titanmod_ps_PerSlipSys(j,instance)
Final version of titanmod
2011-03-31 14:51:43 +05:30
if((1.0_pReal-StressRatio_screw_p)>0.001_pReal) then
minusStressRatio_screw_p=1.0_pReal-StressRatio_screw_p
else
minusStressRatio_screw_p=0.001_pReal
endif
StressRatio_screw_pminus1 = ((abs(tau_slip(j)))/( bottomstress_screw))** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_titanmod_ps_PerSlipSys(j,instance)-1.0_pReal)
Final version of titanmod
2011-03-31 14:51:43 +05:30
!* Boltzmann ratio for screw
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatioscrew = constitutive_titanmod_f0_PerSlipSys(j,instance)/(kB*Temperature)
Final version of titanmod
2011-03-31 14:51:43 +05:30
endif
!* Stress ratio for edge
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
bottomstress_edge=constitutive_titanmod_tau0e_PerSlipSys(j,instance)+ &
plasticState(ph)%state(5*ns+nt+j, of)
Final version of titanmod
2011-03-31 14:51:43 +05:30
StressRatio_edge_p = ((abs(tau_slip(j)))/ &
( bottomstress_edge) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**constitutive_titanmod_pe_PerSlipSys(j,instance)
Final version of titanmod
2011-03-31 14:51:43 +05:30
if((1.0_pReal-StressRatio_edge_p)>0.001_pReal) then
minusStressRatio_edge_p=1.0_pReal-StressRatio_edge_p
else
minusStressRatio_edge_p=0.001_pReal
endif
StressRatio_edge_pminus1 = ((abs(tau_slip(j)))/( bottomstress_edge))** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_titanmod_pe_PerSlipSys(j,instance)-1.0_pReal)
Final version of titanmod
2011-03-31 14:51:43 +05:30
!* Boltzmann ratio for edge. For screws it is defined above
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatioedge = constitutive_titanmod_f0_PerSlipSys(j,instance)/(kB*Temperature)
Final version of titanmod
2011-03-31 14:51:43 +05:30
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
screw_velocity(j) =screwvelocity_prefactor * & ! there is no v0 for screw now because it is included in the prefactor
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
exp(-BoltzmannRatioscrew*(minusStressRatio_screw_p)** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_qs_PerSlipSys(j,instance))
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
edge_velocity(j) =constitutive_titanmod_v0e_PerSlipSys(j,instance)*exp(-BoltzmannRatioedge* &
Final version of titanmod
2011-03-31 14:51:43 +05:30
(minusStressRatio_edge_p)** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_qe_PerSlipSys(j,instance))
Final version of titanmod
2011-03-31 14:51:43 +05:30
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
!* Shear rates due to edge slip
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
gdot_slip_edge(j) = constitutive_titanmod_burgersPerSlipSys(j,instance)*(plasticState(ph)%state(j, of)* &
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
edge_velocity(j))* sign(1.0_pReal,tau_slip(j))
!* Shear rates due to screw slip
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
gdot_slip_screw(j) = constitutive_titanmod_burgersPerSlipSys(j,instance)*(plasticState(ph)%state(ns+j, of) * &
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
screw_velocity(j))* sign(1.0_pReal,tau_slip(j))
!Total shear rate
gdot_slip(j) = gdot_slip_edge(j) + gdot_slip_screw(j)
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state( 7*ns+2*nt+j, of)= edge_velocity(j)
plasticState(ph)%state( 8*ns+2*nt+j, of)= screw_velocity(j)
plasticState(ph)%state( 9*ns+2*nt+j, of)= tau_slip(j)
plasticState(ph)%state(10*ns+2*nt+j, of)= gdot_slip_edge(j)
plasticState(ph)%state(11*ns+2*nt+j, of)= gdot_slip_screw(j)
plasticState(ph)%state(12*ns+2*nt+j, of)= StressRatio_edge_p
plasticState(ph)%state(13*ns+2*nt+j, of)= StressRatio_screw_p
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Derivatives of shear rates
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
dgdot_dtauslip(j) = constitutive_titanmod_burgersPerSlipSys(j,instance)*(( &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
( &
( &
( &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(edge_velocity(j)*plasticState(ph)%state(j, of))) * &
Final version of titanmod
2011-03-31 14:51:43 +05:30
BoltzmannRatioedge*&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_pe_PerSlipSys(j,instance)* &
constitutive_titanmod_qe_PerSlipSys(j,instance) &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
)/ &
Final version of titanmod
2011-03-31 14:51:43 +05:30
bottomstress_edge &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
)*&
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
StressRatio_edge_pminus1*(minusStressRatio_edge_p)** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_titanmod_qe_PerSlipSys(j,instance)-1.0_pReal) &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
) + &
( &
( &
( &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(plasticState(ph)%state(ns+j, of) * screw_velocity(j)) * &
Final version of titanmod
2011-03-31 14:51:43 +05:30
BoltzmannRatioscrew* &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_ps_PerSlipSys(j,instance)* &
constitutive_titanmod_qs_PerSlipSys(j,instance) &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
)/ &
Final version of titanmod
2011-03-31 14:51:43 +05:30
bottomstress_screw &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
)*&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
StressRatio_screw_pminus1*(minusStressRatio_screw_p)**(constitutive_titanmod_qs_PerSlipSys(j,instance)-1.0_pReal) &
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
) &
) !* sign(1.0_pReal,tau_slip(j))
Final version of titanmod
2011-03-31 14:51:43 +05:30
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
!*************************************************
Brushed up accountability of twinning to Lp
2010-09-07 20:14:37 +05:30
!sumf=0.0_pReal
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Plastic velocity gradient for dislocation glide
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Lp = Lp + (1.0_pReal - sumf)*gdot_slip(j)*lattice_Sslip(1:3,1:3,1,index_myFamily+i,ph)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Calculation of the tangent of Lp
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
dLp_dTstar3333(k,l,m,n) = &
dLp_dTstar3333(k,l,m,n) + dgdot_dtauslip(j)*&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_Sslip(k,l,1,index_myFamily+i,ph)*&
lattice_Sslip(m,n,1,index_myFamily+i,ph)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo slipFamiliesLoop
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Mechanical twinning part
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
gdot_twin = 0.0_pReal
dgdot_dtautwin = 0.0_pReal
j = 0_pInt
twinFamiliesLoop: do f = 1_pInt,lattice_maxNtwinFamily
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
do i = 1_pInt,constitutive_titanmod_Ntwin(f,instance) ! process each (active) slip system in family
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
j = j+1_pInt
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Calculation of Lp
!* Resolved shear stress on twin system
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
tau_twin(j) = dot_product(Tstar_v,lattice_Stwin_v(:,index_myFamily+i,ph))
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!**************************************************************************************
!* Stress ratios
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! StressRatio_r = (plasticState(ph)%state6*ns+3*nt+j, of)/tau_twin(j))**constitutive_titanmod_r(instance)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
!* Shear rates and their derivatives due to twin
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
! if ( tau_twin(j) > 0.0_pReal ) !then
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
! gdot_twin(j) = 0.0_pReal!&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! (constitutive_titanmod_MaxTwinFraction(instance)-sumf)*lattice_shearTwin(index_myFamily+i,ph)*&
! plasticState(ph)%state(6*ns+4*nt+j, of)*constitutive_titanmod_Ndot0PerTwinSys(f,instance)*exp(-StressRatio_r)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
! dgdot_dtautwin(j) = ((gdot_twin(j)*constitutive_titanmod_r(instance))/tau_twin(j))*StressRatio_r
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
! endif
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!**************************************************************************************
!* Stress ratio for edge
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinStressRatio_p = ((abs(tau_twin(j)))/ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
( constitutive_titanmod_twintau0_PerTwinSys(j,instance)+plasticState(ph)%state(7*ns+nt+j, of)) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**constitutive_titanmod_twinp_PerTwinSys(j,instance)
Final version of titanmod
2011-03-31 14:51:43 +05:30
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
if((1.0_pReal-twinStressRatio_p)>0.001_pReal) then
twinminusStressRatio_p=1.0_pReal-twinStressRatio_p
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
else
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinminusStressRatio_p=0.001_pReal
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
endif
Final version of titanmod
2011-03-31 14:51:43 +05:30
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinStressRatio_pminus1 = ((abs(tau_twin(j)))/ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
( constitutive_titanmod_twintau0_PerTwinSys(j,instance)+plasticState(ph)%state(7*ns+nt+j, of)) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**(constitutive_titanmod_twinp_PerTwinSys(j,instance)-1.0_pReal)
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!* Boltzmann ratio
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatiotwin = constitutive_titanmod_twinf0_PerTwinSys(j,instance)/(kB*Temperature)
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
!* Initial twin shear rates
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
TwinDotGamma0 = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twingamma0_PerTwinSys(j,instance)
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
!* Shear rates due to twin
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
gdot_twin(j) =sign(1.0_pReal,tau_twin(j))*constitutive_titanmod_twingamma0_PerTwinSys(j,instance)* &
exp(-BoltzmannRatiotwin*(twinminusStressRatio_p)**constitutive_titanmod_twinq_PerTwinSys(j,instance))
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
Final version of titanmod
2011-03-31 14:51:43 +05:30
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!* Derivatives of shear rates in twin
dgdot_dtautwin(j) = ( &
( &
( &
(abs(gdot_twin(j))) * &
Final version of titanmod
2011-03-31 14:51:43 +05:30
BoltzmannRatiotwin*&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinp_PerTwinSys(j,instance)* &
constitutive_titanmod_twinq_PerTwinSys(j,instance) &
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twintau0_PerTwinSys(j,instance) &
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
)*&
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinStressRatio_pminus1*(twinminusStressRatio_p)** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_titanmod_twinq_PerTwinSys(j,instance)-1.0_pReal) &
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
) !* sign(1.0_pReal,tau_slip(j))
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 17:06:51 +05:30
!* Plastic velocity gradient for mechanical twinning
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! Lp = Lp + sumf*gdot_twin(j)*lattice_Stwin(:,:,index_myFamily+i,ph)
Lp = Lp + gdot_twin(j)*lattice_Stwin(:,:,index_myFamily+i,ph)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
!* Calculation of the tangent of Lp
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
dLp_dTstar3333(k,l,m,n) = &
dLp_dTstar3333(k,l,m,n) + dgdot_dtautwin(j)*&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_Stwin(k,l,index_myFamily+i,ph)*&
lattice_Stwin(m,n,index_myFamily+i,ph)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo twinFamiliesLoop
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
dLp_dTstar99 = math_Plain3333to99(dLp_dTstar3333)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! tempState=state
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
end subroutine constitutive_titanmod_LpAndItsTangent
!--------------------------------------------------------------------------------------------------
!> @brief calculates the rate of change of microstructure
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_titanmod_dotState(Tstar_v,temperature,ipc,ip,el)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use lattice, only: &
lattice_Stwin_v, &
lattice_maxNslipFamily, &
lattice_maxNtwinFamily, &
lattice_NslipSystem, &
lattice_NtwinSystem
comments for doxygen
2013-07-12 12:27:15 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_plasticityInstance, &
plasticState, &
mappingConstitutive
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +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
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), dimension(6), intent(in):: &
comments for doxygen
2013-07-12 12:27:15 +05:30
Tstar_v !< 2nd Piola Kirchhoff stress tensor in Mandel notation
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
real(pReal), intent(in) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
temperature !< temperature at integration point
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
integer(pInt), intent(in) :: &
comments for doxygen
2013-07-12 12:27:15 +05:30
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily, instance, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ns,nt,&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
f,i,j, &
ph, &
of
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal) :: &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
sumf,BoltzmannRatio, &
twinStressRatio_p,twinminusStressRatio_p
real(pReal), dimension(constitutive_titanmod_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
DotRhoEdgeGeneration, &
DotRhoEdgeAnnihilation, &
DotRhoScrewGeneration, &
DotRhoScrewAnnihilation
real(pReal), dimension(constitutive_titanmod_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_twin, &
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
tau_twin, &
volumefraction_PerTwinSys
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! shortened notation
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
of = mappingConstitutive(1,ipc,ip,el)
ph = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(ph)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
do i=1_pInt,nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
volumefraction_PerTwinSys(i)=plasticState(ph)%state(3_pInt*ns+i, of)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(i,instance)
Brushed up accountability of twinning to Lp
2010-09-07 20:14:37 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
enddo
Brushed up accountability of twinning to Lp
2010-09-07 20:14:37 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
sumf = sum(abs(volumefraction_PerTwinSys(1_pInt:nt))) ! safe for nt == 0
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(:,of) = 0.0_pReal
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
j = 0_pInt
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
slipFamiliesLoop: do f = 1_pInt,lattice_maxNslipFamily
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
do i = 1_pInt,constitutive_titanmod_Nslip(f,instance) ! process each (active) slip system in family
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
j = j+1_pInt
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
DotRhoEdgeGeneration(j) = & ! multiplication of edge dislocations
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(ns+j, of)*plasticState(ph)%state(8*ns+2*nt+j, of)/plasticState(ph)%state(4*ns+nt+j, of)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
DotRhoScrewGeneration(j) = & ! multiplication of screw dislocations
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*plasticState(ph)%state(7*ns+2*nt+j, of)/plasticState(ph)%state(3*ns+nt+j, of)
DotRhoEdgeAnnihilation(j) = -((plasticState(ph)%state(j, of))**2)* & ! annihilation of edge dislocations
constitutive_titanmod_capre_PerSlipSys(j,instance)*plasticState(ph)%state(7*ns+2*nt+j, of)*0.5_pReal
DotRhoScrewAnnihilation(j) = -((plasticState(ph)%state(ns+j, of))**2)* & ! annihilation of screw dislocations
constitutive_titanmod_caprs_PerSlipSys(j,instance)*plasticState(ph)%state(8*ns+2*nt+j, of)*0.5_pReal
plasticState(ph)%dotState(j, of) = & ! edge dislocation density rate of change
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
DotRhoEdgeGeneration(j)+DotRhoEdgeAnnihilation(j)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(ns+j, of) = & ! screw dislocation density rate of change
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
DotRhoScrewGeneration(j)+DotRhoScrewAnnihilation(j)
Calibration (Prismatic slip) of the new dislocation density based CPFE model for single crystal alpha-Ti (Titanmod). Formulated edge and screw dislocations mobilities in presence of kink pairs. Updated comments. The model shows the plateau in stress evolution in prismatic orientation which is observed in experiments.
2010-07-15 12:46:15 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(2*ns+j, of) = & ! sum of shear due to edge and screw
plasticState(ph)%state(10*ns+2*nt+j, of)+plasticState(ph)%state(11*ns+2*nt+j, of)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo slipFamiliesLoop
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!* Twin fraction evolution
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
j = 0_pInt
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
twinFamiliesLoop: do f = 1_pInt,lattice_maxNtwinFamily
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
do i = 1_pInt,constitutive_titanmod_Ntwin(f,instance) ! process each (active) twin system in family
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
j = j+1_pInt
!* Resolved shear stress on twin system
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
tau_twin(j) = dot_product(Tstar_v,lattice_Stwin_v(:,index_myFamily+i,ph))
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
!* Stress ratio for edge
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinStressRatio_p = ((abs(tau_twin(j)))/ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
( constitutive_titanmod_twintau0_PerTwinSys(j,instance)+plasticState(ph)%state(7*ns+nt+j, of)) &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
)**(constitutive_titanmod_twinp_PerTwinSys(j,instance))
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
if((1.0_pReal-twinStressRatio_p)>0.001_pReal) then
twinminusStressRatio_p=1.0_pReal-twinStressRatio_p
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
else
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
twinminusStressRatio_p=0.001_pReal
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 15:37:52 +05:30
endif
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatio = constitutive_titanmod_twinf0_PerTwinSys(j,instance)/(kB*Temperature)
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
gdot_twin(j) =constitutive_titanmod_twingamma0_PerTwinSys(j,instance)*exp(-BoltzmannRatio* &
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
(twinminusStressRatio_p)** &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinq_PerTwinSys(j,instance))*sign(1.0_pReal,tau_twin(j))
Cleaned-up, rearranged state variables
2010-09-22 14:30:40 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(3*ns+j, of)=gdot_twin(j)
Brushed up accountability of twinning to Lp
2010-09-07 20:14:37 +05:30
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo twinFamiliesLoop
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
end subroutine constitutive_titanmod_dotState
added helper functions to communicate accumulated slip between damage and plasticity modules
2014-10-10 22:04:51 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns accumulated slip
!--------------------------------------------------------------------------------------------------
subroutine constitutive_titanmod_getAccumulatedSlip(nSlip,accumulatedSlip,ipc, ip, el)
use lattice, only: &
lattice_maxNslipFamily
use material, only: &
mappingConstitutive, &
plasticState, &
phase_plasticityInstance
implicit none
real(pReal), dimension(:), allocatable :: &
accumulatedSlip
integer(pInt) :: &
nSlip
integer(pInt), intent(in) :: &
ipc, & !< grain number
ip, & !< integration point number
el !< element number
integer(pInt) :: &
offset, &
phase, &
instance, &
offset_accshear_slip, &
f, j, i
offset = mappingConstitutive(1,ipc,ip,el)
phase = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(phase)
nSlip = constitutive_titanmod_totalNslip(instance)
allocate(accumulatedSlip(nSlip))
offset_accshear_slip = 2_pInt*nSlip
j = 0_pInt
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
do i = 1_pInt,constitutive_titanmod_Nslip(f,instance) ! process each (active) slip system in family
j = j+1_pInt
accumulatedSlip(j) = plasticState(phase)%state(offset_accshear_slip+j,offset)
enddo
enddo
end subroutine constitutive_titanmod_getAccumulatedSlip
- added get slip rate helper function - added vacancy generation rate - cleaned up some bugs
2014-10-11 15:15:30 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief returns accumulated slip
!--------------------------------------------------------------------------------------------------
subroutine constitutive_titanmod_getSlipRate(nSlip,slipRate,ipc, ip, el)
use lattice, only: &
lattice_maxNslipFamily
use material, only: &
mappingConstitutive, &
plasticState, &
phase_plasticityInstance
implicit none
real(pReal), dimension(:), allocatable :: &
slipRate
integer(pInt) :: &
nSlip
integer(pInt), intent(in) :: &
ipc, & !< grain number
ip, & !< integration point number
el !< element number
integer(pInt) :: &
offset, &
phase, &
instance, &
offset_accshear_slip, &
f, j, i
offset = mappingConstitutive(1,ipc,ip,el)
phase = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(phase)
nSlip = constitutive_titanmod_totalNslip(instance)
allocate(slipRate(nSlip))
offset_accshear_slip = 2_pInt*nSlip
j = 0_pInt
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
do i = 1_pInt,constitutive_titanmod_Nslip(f,instance) ! process each (active) slip system in family
j = j+1_pInt
slipRate(j) = plasticState(phase)%dotState(offset_accshear_slip+j,offset)
enddo
enddo
end subroutine constitutive_titanmod_getSlipRate
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief return array of constitutive results
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
function constitutive_titanmod_postResults(ipc,ip,el)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use mesh, only: &
mesh_NcpElems, &
mesh_maxNips
use material, only: &
homogenization_maxNgrains, &
material_phase, &
phase_plasticityInstance, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
phase_Noutput, &
plasticState, &
mappingConstitutive
comments for doxygen
2013-07-12 12:27:15 +05:30
implicit none
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 12:13:49 +05:30
integer(pInt), intent(in) :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
real(pReal), dimension(constitutive_titanmod_sizePostResults(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
constitutive_titanmod_postResults
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
integer(pInt) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
instance, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
ns,nt,&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
o,i,c, &
ph, &
of
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
real(pReal) :: sumf
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
comments for doxygen
2013-07-12 12:27:15 +05:30
real(pReal), dimension(constitutive_titanmod_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
volumefraction_PerTwinSys
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! shortened notation
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
of = mappingConstitutive(1,ipc,ip,el)
ph = mappingConstitutive(2,ipc,ip,el)
instance = phase_plasticityInstance(ph)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
ns = constitutive_titanmod_totalNslip(instance)
nt = constitutive_titanmod_totalNtwin(instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
do i=1_pInt,nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
volumefraction_PerTwinSys(i)=plasticState(ph)%state(3_pInt*ns+i, of)/ &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_titanmod_twinshearconstant_PerTwinSys(i,instance)
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
sumf = sum(abs(volumefraction_PerTwinSys(1:nt))) ! safe for nt == 0
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! required output
comments for doxygen
2013-07-12 12:27:15 +05:30
c = 0_pInt
constitutive_titanmod_postResults = 0.0_pReal
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....)
2014-06-25 04:23:25 +05:30
do o = 1_pInt,constitutive_titanmod_Noutput(instance)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
select case(constitutive_titanmod_outputID(o,instance))
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoedge_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(1_pInt:ns, of)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoscrew_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of)
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (segment_edge_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(3_pInt*ns+nt+1_pInt:4_pInt*ns+nt, of)
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (segment_screw_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(4_pInt*ns+nt+1_pInt:5_pInt*ns+nt, of)
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (resistance_edge_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(5_pInt*ns+nt+1_pInt:6_pInt*ns+nt, of)
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 21:23:45 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (resistance_screw_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(6_pInt*ns+nt+1_pInt:7_pInt*ns+nt, of)
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (velocity_edge_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(7*ns+2*nt+1:8*ns+2*nt, of)
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (velocity_screw_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(8*ns+2*nt+1:9*ns+2*nt, of)
Final version of titanmod
2011-03-31 14:51:43 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (tau_slip_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(9*ns+2*nt+1:10*ns+2*nt, of))
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (gdot_slip_edge_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(10*ns+2*nt+1:11*ns+2*nt, of))
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (gdot_slip_screw_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(11*ns+2*nt+1:12*ns+2*nt, of))
Final version of titanmod
2011-03-31 14:51:43 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (gdot_slip_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(10*ns+2*nt+1:11*ns+2*nt, of)) + &
abs(plasticState(ph)%state(11*ns+2*nt+1:12*ns+2*nt, of))
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (stressratio_edge_p_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(12*ns+2*nt+1:13*ns+2*nt, of))
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (stressratio_screw_p_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(13*ns+2*nt+1:14*ns+2*nt, of))
Added more state dependent variables in post results.
2010-09-29 12:05:08 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_system_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+ns) = abs(plasticState(ph)%state(2*ns+1:3*ns, of))
Final version of titanmod
2011-03-31 14:51:43 +05:30
c = c + ns
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_basal_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(abs(plasticState(ph)%state(2*ns+1:2*ns+3, of)))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_prism_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(abs(plasticState(ph)%state(2*ns+4:2*ns+6, of)))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_pyra_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(abs(plasticState(ph)%state(2*ns+7:2*ns+12, of)))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_pyrca_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(abs(plasticState(ph)%state(2*ns+13:2*ns+24, of)))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
Final version of titanmod
2011-03-31 14:51:43 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoedge_basal_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(1:3, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoedge_prism_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(4:6, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoedge_pyra_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(7:12,of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoedge_pyrca_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(13:24, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
Final version of titanmod
2011-03-31 14:51:43 +05:30
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoscrew_basal_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(ns+1:ns+3, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoscrew_prism_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(ns+4:ns+6, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoscrew_pyra_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(ns+7:ns+12, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (rhoscrew_pyrca_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(plasticState(ph)%state(ns+13:ns+24, of))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (shear_total_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+1_pInt) = sum(abs(plasticState(ph)%state(2*ns+1:3*ns, of)))
polishing: adding pInt, removing unused use-statements etc
2012-02-21 21:30:00 +05:30
c = c + 1_pInt
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 22:38:13 +05:30
case (twin_fraction_ID)
comments for doxygen
2013-07-12 12:27:15 +05:30
constitutive_titanmod_postResults(c+1_pInt:c+nt) = abs(volumefraction_PerTwinSys(1:nt))
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
c = c + nt
end select
comments for doxygen
2013-07-12 12:27:15 +05:30
enddo
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
end function constitutive_titanmod_postResults
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 13:49:25 +05:30
comments for doxygen
2013-07-12 12:27:15 +05:30
end module constitutive_titanmod