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_dislotwin.f90

2357 lines
135 KiB
Fortran
Raw Normal View History

added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! $Id$
!--------------------------------------------------------------------------------------------------
!> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
!> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
!> @brief material subroutine incoprorating dislocation and twinning physics
!> @details to be done
!--------------------------------------------------------------------------------------------------
module constitutive_dislotwin
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
use prec, only: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
pReal, &
pInt
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
implicit none
private
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
integer(pInt), dimension(:), allocatable, public, protected :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_sizePostResults !< cumulative size of post results
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
integer(pInt), dimension(:,:), allocatable, target, public :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_sizePostResult !< size of each post result output
edited lines exceeding 132 chars (before any comment started)
2012-01-11 22:26:35 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
character(len=64), dimension(:,:), allocatable, target, public :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_output !< name of each post result output
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
character(len=12), dimension(3), parameter, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
CONSTITUTIVE_DISLOTWIN_listBasicSlipStates = &
['rhoEdge ', 'rhoEdgeDip ', 'accshearslip']
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
character(len=12), dimension(2), parameter, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
CONSTITUTIVE_DISLOTWIN_listBasicTwinStates = &
['twinFraction', 'accsheartwin']
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
character(len=19), dimension(2), parameter, private :: &
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
CONSTITUTIVE_DISLOTWIN_listBasicTransStates = &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
['stressTransFraction', 'strainTransFraction']
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
character(len=17), dimension(4), parameter, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
CONSTITUTIVE_DISLOTWIN_listDependentSlipStates = &
['invLambdaSlip ', 'invLambdaSlipTwin', 'meanFreePathSlip ', 'tauSlipThreshold ']
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
character(len=16), dimension(4), parameter, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
CONSTITUTIVE_DISLOTWIN_listDependentTwinStates = &
['invLambdaTwin ', 'meanFreePathTwin', 'tauTwinThreshold', 'twinVolume ']
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), parameter, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
kB = 1.38e-23_pReal !< Boltzmann constant in J/Kelvin
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 :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_Noutput !< number of outputs per instance of this plasticity
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
integer(pInt), dimension(:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_totalNslip, & !< total number of active slip systems for each instance
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
constitutive_dislotwin_totalNtwin, & !< total number of active twin systems for each instance
constitutive_dislotwin_totalNtrans !< number of active transformation systems
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
integer(pInt), dimension(:,:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_Nslip, & !< number of active slip systems for each family and instance
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
constitutive_dislotwin_Ntwin, & !< number of active twin systems for each family and instance
constitutive_dislotwin_Ntrans !< number of active transformation systems for each family and instance
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_CAtomicVolume, & !< atomic volume in Bugers vector unit
constitutive_dislotwin_D0, & !< prefactor for self-diffusion coefficient
constitutive_dislotwin_Qsd, & !< activation energy for dislocation climb
constitutive_dislotwin_GrainSize, & !< grain size
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_pShearBand, & !< p-exponent in shearband velocity
constitutive_dislotwin_qShearBand, & !< q-exponent in shearband velocity
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_MaxTwinFraction, & !< maximum allowed total twin volume fraction
constitutive_dislotwin_CEdgeDipMinDistance, & !<
constitutive_dislotwin_Cmfptwin, & !<
constitutive_dislotwin_Cthresholdtwin, & !<
constitutive_dislotwin_SolidSolutionStrength, & !< Strength due to elements in solid solution
constitutive_dislotwin_L0, & !< Length of twin nuclei in Burgers vectors
constitutive_dislotwin_xc, & !< critical distance for formation of twin nucleus
constitutive_dislotwin_VcrossSlip, & !< cross slip volume
constitutive_dislotwin_sbResistance, & !< value for shearband resistance (might become an internal state variable at some point)
constitutive_dislotwin_sbVelocity, & !< value for shearband velocity_0
constitutive_dislotwin_sbQedge, & !< value for shearband systems Qedge
constitutive_dislotwin_SFE_0K, & !< stacking fault energy at zero K
constitutive_dislotwin_dSFE_dT, & !< temperature dependance of stacking fault energy
added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given
2014-03-14 05:20:55 +05:30
constitutive_dislotwin_dipoleFormationFactor, & !< scaling factor for dipole formation: 0: off, 1: on. other values not useful
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_aTolRho, & !< absolute tolerance for integration of dislocation density
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
constitutive_dislotwin_aTolTwinFrac, & !< absolute tolerance for integration of twin volume fraction
Renaming transformation input variables to more logical names
2014-09-23 21:08:41 +05:30
constitutive_dislotwin_Cnuc, & !< Coefficient for strain-induced martensite nucleation
constitutive_dislotwin_Cdwp, & !< Coefficient for double well potential
constitutive_dislotwin_Cgro, & !< Coefficient for stress-assisted martensite growth
constitutive_dislotwin_deltaG !< Free energy difference between austensite and martensite [MPa]
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:,:,:,:), allocatable, private :: &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
constitutive_dislotwin_Ctwin66 !< twin elasticity matrix in Mandel notation for each instance
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:,:,:,:,:,:), allocatable, private :: &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
constitutive_dislotwin_Ctwin3333 !< twin elasticity matrix for each instance
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
real(pReal), dimension(:,:,:,:), allocatable, private :: &
constitutive_dislotwin_Ctrans66 !< trans elasticity matrix in Mandel notation for each instance
real(pReal), dimension(:,:,:,:,:,:), allocatable, private :: &
constitutive_dislotwin_Ctrans3333 !< trans elasticity matrix for each instance
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:,:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_rhoEdge0, & !< initial edge dislocation density per slip system for each family and instance
constitutive_dislotwin_rhoEdgeDip0, & !< initial edge dipole density per slip system for each family and instance
constitutive_dislotwin_burgersPerSlipFamily, & !< absolute length of burgers vector [m] for each slip family and instance
constitutive_dislotwin_burgersPerSlipSystem, & !< absolute length of burgers vector [m] for each slip system and instance
constitutive_dislotwin_burgersPerTwinFamily, & !< absolute length of burgers vector [m] for each twin family and instance
constitutive_dislotwin_burgersPerTwinSystem, & !< absolute length of burgers vector [m] for each twin system and instance
constitutive_dislotwin_QedgePerSlipFamily, & !< activation energy for glide [J] for each slip family and instance
constitutive_dislotwin_QedgePerSlipSystem, & !< activation energy for glide [J] for each slip system and instance
constitutive_dislotwin_v0PerSlipFamily, & !< dislocation velocity prefactor [m/s] for each family and instance
constitutive_dislotwin_v0PerSlipSystem, & !< dislocation velocity prefactor [m/s] for each slip system and instance
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_tau_peierlsPerSlipFamily, & !< Peierls stress [Pa] for each family and instance
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_Ndot0PerTwinFamily, & !< twin nucleation rate [1/m³s] for each twin family and instance
constitutive_dislotwin_Ndot0PerTwinSystem, & !< twin nucleation rate [1/m³s] for each twin system and instance
constitutive_dislotwin_tau_r, & !< stress to bring partial close together for each twin system and instance
constitutive_dislotwin_twinsizePerTwinFamily, & !< twin thickness [m] for each twin family and instance
constitutive_dislotwin_twinsizePerTwinSystem, & !< twin thickness [m] for each twin system and instance
constitutive_dislotwin_CLambdaSlipPerSlipFamily, & !< Adj. parameter for distance between 2 forest dislocations for each slip family and instance
constitutive_dislotwin_CLambdaSlipPerSlipSystem, & !< Adj. parameter for distance between 2 forest dislocations for each slip system and instance
constitutive_dislotwin_interaction_SlipSlip, & !< coefficients for slip-slip interaction for each interaction type and instance
constitutive_dislotwin_interaction_SlipTwin, & !< coefficients for slip-twin interaction for each interaction type and instance
constitutive_dislotwin_interaction_TwinSlip, & !< coefficients for twin-slip interaction for each interaction type and instance
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_interaction_TwinTwin, & !< coefficients for twin-twin interaction for each interaction type and instance
constitutive_dislotwin_pPerSlipFamily, & !< p-exponent in glide velocity
constitutive_dislotwin_qPerSlipFamily, & !< q-exponent in glide velocity
constitutive_dislotwin_rPerTwinFamily !< r-exponent in twin nucleation rate
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:,:,:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_interactionMatrix_SlipSlip, & !< interaction matrix of the different slip systems for each instance
constitutive_dislotwin_interactionMatrix_SlipTwin, & !< interaction matrix of slip systems with twin systems for each instance
constitutive_dislotwin_interactionMatrix_TwinSlip, & !< interaction matrix of twin systems with slip systems for each instance
constitutive_dislotwin_interactionMatrix_TwinTwin, & !< interaction matrix of the different twin systems for each instance
Projection matrices for transformation
2014-09-24 18:01:16 +05:30
constitutive_dislotwin_forestProjectionEdge, & !< matrix of forest projections of edge dislocations for each instance
constitutive_dislotwin_projectionMatrix_Trans !< matrix for projection of slip system shear on fault band (twin) systems for each instance
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
real(pReal), dimension(:,:,:,:,:), allocatable, private :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_sbSv
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
enum, bind(c)
enumerator :: undefined_ID, &
edge_density_ID, &
dipole_density_ID, &
shear_rate_slip_ID, &
accumulated_shear_slip_ID, &
mfp_slip_ID, &
resolved_stress_slip_ID, &
threshold_stress_slip_ID, &
edge_dipole_distance_ID, &
stress_exponent_ID, &
twin_fraction_ID, &
shear_rate_twin_ID, &
accumulated_shear_twin_ID, &
mfp_twin_ID, &
resolved_stress_twin_ID, &
threshold_stress_twin_ID, &
resolved_stress_shearband_ID, &
shear_rate_shearband_ID, &
sb_eigenvalues_ID, &
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
sb_eigenvectors_ID, &
stress_trans_fraction_ID, &
strain_trans_fraction_ID, &
trans_fraction_ID
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
end enum
integer(kind(undefined_ID)), 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_dislotwin_outputID !< ID of each post result output
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
public :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_init, &
constitutive_dislotwin_homogenizedC, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_microstructure, &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_LpAndItsTangent, &
constitutive_dislotwin_dotState, &
added helper functions to communicate accumulated slip between damage and plasticity modules
2014-10-10 22:04:51 +05:30
constitutive_dislotwin_getAccumulatedSlip, &
- added get slip rate helper function - added vacancy generation rate - cleaned up some bugs
2014-10-11 15:15:30 +05:30
constitutive_dislotwin_getSlipRate, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_postResults
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
private :: &
constitutive_dislotwin_stateInit, &
constitutive_dislotwin_aTolState
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
contains
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
!> @brief module initialization
!> @details reads in material parameters, allocates arrays, and does sanity checks
!--------------------------------------------------------------------------------------------------
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
subroutine constitutive_dislotwin_init(fileUnit)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
use, intrinsic :: iso_fortran_env ! to get compiler_version and compiler_options (at least for gfortran 4.6 at the moment)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use debug, only: &
debug_level,&
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
debug_constitutive,&
debug_levelBasic
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use math, only: &
math_Mandel3333to66, &
math_Voigt66to3333, &
math_mul3x3
use mesh, only: &
mesh_maxNips, &
mesh_NcpElems
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
use IO, only: &
moved reading in of lattice stru
2014-02-10 20:01:19 +05:30
IO_read, &
IO_lc, &
IO_getTag, &
IO_isBlank, &
IO_stringPos, &
IO_stringValue, &
IO_floatValue, &
IO_intValue, &
IO_warning, &
IO_error, &
IO_timeStamp, &
IO_EOF
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
use material, only: &
homogenization_maxNgrains, &
phase_plasticity, &
phase_plasticityInstance, &
phase_Noutput, &
PLASTICITY_DISLOTWIN_label, &
PLASTICITY_DISLOTWIN_ID, &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
material_phase, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState, &
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 14:19:47 +05:30
MATERIAL_partPhase
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use lattice
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
use numerics,only: &
only output from the root processor for parallel runs
2014-10-10 01:53:06 +05:30
worldrank, &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
numerics_integrator
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
implicit none
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
integer(pInt), intent(in) :: fileUnit
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
integer(pInt), dimension(1+2*MAXNCHUNKS) :: positions
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
integer(pInt) :: maxNinstance,mySize=0_pInt,phase,maxTotalNslip,maxTotalNtwin,maxTotalNtrans,&
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
f,instance,j,k,l,m,n,o,p,q,r,s,ns,nt,nr, &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
Nchunks_SlipSlip, Nchunks_SlipTwin, Nchunks_TwinSlip, Nchunks_TwinTwin, &
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
Nchunks_SlipFamilies, Nchunks_TwinFamilies, Nchunks_TransFamilies, &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
index_myFamily, index_otherFamily
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
integer(pInt) :: sizeState, sizeDotState
integer(pInt) :: NofMyPhase
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
character(len=65536) :: &
tag = '', &
line = ''
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
real(pReal), dimension(:), allocatable :: tempPerSlip, tempPerTwin, tempPerTrans
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
reinstating changes lost in previous commit
2014-10-10 22:15:14 +05:30
mainProcess: if (worldrank == 0) then
write(6,'(/,a)') ' <<<+- constitutive_'//PLASTICITY_DISLOTWIN_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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
maxNinstance = int(count(phase_plasticity == PLASTICITY_DISLOTWIN_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
if (iand(debug_level(debug_constitutive),debug_levelBasic) /= 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
write(6,'(a16,1x,i5,/)') '# instances:',maxNinstance
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_sizePostResults(maxNinstance), source=0_pInt)
allocate(constitutive_dislotwin_sizePostResult(maxval(phase_Noutput),maxNinstance),source=0_pInt)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
allocate(constitutive_dislotwin_output(maxval(phase_Noutput),maxNinstance))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_output = ''
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_outputID(maxval(phase_Noutput),maxNinstance), source=undefined_ID)
allocate(constitutive_dislotwin_Noutput(maxNinstance), source=0_pInt)
allocate(constitutive_dislotwin_Nslip(lattice_maxNslipFamily,maxNinstance), source=0_pInt)
allocate(constitutive_dislotwin_Ntwin(lattice_maxNtwinFamily,maxNinstance), source=0_pInt)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
allocate(constitutive_dislotwin_Ntrans(lattice_maxNtransFamily,maxNinstance), source=0_pInt)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_totalNslip(maxNinstance), source=0_pInt)
allocate(constitutive_dislotwin_totalNtwin(maxNinstance), source=0_pInt)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
allocate(constitutive_dislotwin_totalNtrans(maxNinstance), source=0_pInt)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_CAtomicVolume(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_D0(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Qsd(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_GrainSize(maxNinstance), source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(constitutive_dislotwin_pShearBand(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_qShearBand(maxNinstance), source=0.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_MaxTwinFraction(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_CEdgeDipMinDistance(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Cmfptwin(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Cthresholdtwin(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_SolidSolutionStrength(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_L0(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_xc(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_VcrossSlip(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_aTolRho(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_aTolTwinFrac(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_sbResistance(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_sbVelocity(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_sbQedge(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_SFE_0K(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_dSFE_dT(maxNinstance), source=0.0_pReal)
added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given
2014-03-14 05:20:55 +05:30
allocate(constitutive_dislotwin_dipoleFormationFactor(maxNinstance), source=1.0_pReal) !should be on by default
Renaming transformation input variables to more logical names
2014-09-23 21:08:41 +05:30
allocate(constitutive_dislotwin_Cnuc(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Cdwp(maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Cgro(maxNinstance), source=0.0_pReal)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
allocate(constitutive_dislotwin_deltaG(maxNinstance), source=0.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_rhoEdge0(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_rhoEdgeDip0(lattice_maxNslipFamily,maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_burgersPerSlipFamily(lattice_maxNslipFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_burgersPerTwinFamily(lattice_maxNtwinFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_QedgePerSlipFamily(lattice_maxNslipFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_v0PerSlipFamily(lattice_maxNslipFamily,maxNinstance), &
source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(constitutive_dislotwin_tau_peierlsPerSlipFamily(lattice_maxNslipFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_pPerSlipFamily(lattice_maxNslipFamily,maxNinstance),source=0.0_pReal)
allocate(constitutive_dislotwin_qPerSlipFamily(lattice_maxNslipFamily,maxNinstance),source=0.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_Ndot0PerTwinFamily(lattice_maxNtwinFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_twinsizePerTwinFamily(lattice_maxNtwinFamily,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_CLambdaSlipPerSlipFamily(lattice_maxNslipFamily,maxNinstance), &
source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(constitutive_dislotwin_rPerTwinFamily(lattice_maxNtwinFamily,maxNinstance),source=0.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_interaction_SlipSlip(lattice_maxNinteraction,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_interaction_SlipTwin(lattice_maxNinteraction,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_interaction_TwinSlip(lattice_maxNinteraction,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_interaction_TwinTwin(lattice_maxNinteraction,maxNinstance), &
source=0.0_pReal)
allocate(constitutive_dislotwin_sbSv(6,6,homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems), &
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
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
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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>
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
line = IO_read(fileUnit)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
line = IO_read(fileUnit)
moved reading in of lattice type and elastic constants to lattice module removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config constitutive modules will only be initialized if needed homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-09 02:20:31 +05:30
if (IO_isBlank(line)) cycle ! skip empty lines
if (IO_getTag(line,'<','>') /= '') then ! stop at next part
line = IO_read(fileUnit, .true.) ! reset IO_read
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
exit
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
if (IO_getTag(line,'[',']') /= '') then ! next phase section
phase = phase + 1_pInt ! advance phase section counter
if (phase_plasticity(phase) == PLASTICITY_DISLOTWIN_ID) then
Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt)
Nchunks_TwinFamilies = count(lattice_NtwinSystem(:,phase) > 0_pInt)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
Nchunks_TransFamilies =count(lattice_NtransSystem(:,phase)> 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
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 introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if(allocated(tempPerSlip)) deallocate(tempPerSlip)
if(allocated(tempPerTwin)) deallocate(tempPerTwin)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
if(allocated(tempPerTrans)) deallocate(tempPerTrans)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(tempPerSlip(Nchunks_SlipFamilies))
allocate(tempPerTwin(Nchunks_TwinFamilies))
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
allocate(tempPerTrans(Nchunks_TransFamilies))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
endif
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +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
if (phase > 0_pInt ) then; if (phase_plasticity(phase) == PLASTICITY_DISLOTWIN_ID) then ! do not short-circuit here (.and. with next if statemen). 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 ('edge_density')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = edge_density_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_dislotwin_output(constitutive_dislotwin_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 ('dipole_density')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = dipole_density_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_dislotwin_output(constitutive_dislotwin_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_rate_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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = shear_rate_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_dislotwin_output(constitutive_dislotwin_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 ('accumulated_shear_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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = accumulated_shear_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_dislotwin_output(constitutive_dislotwin_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 ('mfp_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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = mfp_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_dislotwin_output(constitutive_dislotwin_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 ('resolved_stress_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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = resolved_stress_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_dislotwin_output(constitutive_dislotwin_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 ('edge_dipole_distance')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = edge_dipole_distance_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_dislotwin_output(constitutive_dislotwin_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 ('stress_exponent')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = stress_exponent_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_dislotwin_output(constitutive_dislotwin_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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_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_dislotwin_output(constitutive_dislotwin_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_rate_twin')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = shear_rate_twin_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_dislotwin_output(constitutive_dislotwin_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 ('accumulated_shear_twin')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = accumulated_shear_twin_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_dislotwin_output(constitutive_dislotwin_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 ('mfp_twin')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = mfp_twin_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_dislotwin_output(constitutive_dislotwin_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 ('resolved_stress_twin')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = resolved_stress_twin_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_dislotwin_output(constitutive_dislotwin_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 ('threshold_stress_twin')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = threshold_stress_twin_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_dislotwin_output(constitutive_dislotwin_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 ('resolved_stress_shearband')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = resolved_stress_shearband_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_dislotwin_output(constitutive_dislotwin_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_rate_shearband')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = shear_rate_shearband_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_dislotwin_output(constitutive_dislotwin_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 ('sb_eigenvalues')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = sb_eigenvalues_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_dislotwin_output(constitutive_dislotwin_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 ('sb_eigenvectors')
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_dislotwin_Noutput(instance) = constitutive_dislotwin_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_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = sb_eigenvectors_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_dislotwin_output(constitutive_dislotwin_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
case ('stress_trans_fraction')
constitutive_dislotwin_Noutput(instance) = constitutive_dislotwin_Noutput(instance) + 1_pInt
constitutive_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = stress_trans_fraction_ID
constitutive_dislotwin_output(constitutive_dislotwin_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
case ('strain_trans_fraction')
constitutive_dislotwin_Noutput(instance) = constitutive_dislotwin_Noutput(instance) + 1_pInt
constitutive_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = strain_trans_fraction_ID
constitutive_dislotwin_output(constitutive_dislotwin_Noutput(instance),instance) = &
IO_lc(IO_stringValue(line,positions,2_pInt))
case ('trans_fraction','total_trans_fraction')
constitutive_dislotwin_Noutput(instance) = constitutive_dislotwin_Noutput(instance) + 1_pInt
constitutive_dislotwin_outputID(constitutive_dislotwin_Noutput(instance),instance) = trans_fraction_ID
constitutive_dislotwin_output(constitutive_dislotwin_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
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
!--------------------------------------------------------------------------------------------------
some minor polishing on the way
2014-04-04 20:10:30 +05:30
! parameters depending on number of slip system families
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 ('nslip')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (positions(1) < Nchunks_SlipFamilies + 1_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_warning(50_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (positions(1) > Nchunks_SlipFamilies + 1_pInt) &
call IO_error(150_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_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
Nchunks_SlipFamilies = positions(1) - 1_pInt
do j = 1_pInt, Nchunks_SlipFamilies
some minor polishing on the way
2014-04-04 20:10:30 +05:30
constitutive_dislotwin_Nslip(j,instance) = IO_intValue(line,positions,1_pInt+j)
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
enddo
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
case ('rhoedge0','rhoedgedip0','slipburgers','qedge','v0','clambdaslip','tau_peierls','p_slip','q_slip')
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, Nchunks_SlipFamilies
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
tempPerSlip(j) = IO_floatValue(line,positions,1_pInt+j)
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
enddo
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
select case(tag)
case ('rhoedge0')
constitutive_dislotwin_rhoEdge0(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('rhoedgedip0')
constitutive_dislotwin_rhoEdgeDip0(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('slipburgers')
constitutive_dislotwin_burgersPerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('qedge')
constitutive_dislotwin_QedgePerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('v0')
constitutive_dislotwin_v0PerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('clambdaslip')
constitutive_dislotwin_CLambdaSlipPerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('tau_peierls')
if (lattice_structure(phase) /= LATTICE_bcc_ID) &
call IO_warning(42_pInt,ext_msg=trim(tag)//' for non-bcc ('//PLASTICITY_DISLOTWIN_label//')')
constitutive_dislotwin_tau_peierlsPerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('p_slip')
constitutive_dislotwin_pPerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
case ('q_slip')
constitutive_dislotwin_qPerSlipFamily(1:Nchunks_SlipFamilies,instance) = tempPerSlip(1:Nchunks_SlipFamilies)
end select
!--------------------------------------------------------------------------------------------------
! parameters depending on slip number of twin families
case ('ntwin')
if (positions(1) < Nchunks_TwinFamilies + 1_pInt) &
fixed warning message in constitutive
2014-06-17 20:54:44 +05:30
call IO_warning(51_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (positions(1) > Nchunks_TwinFamilies + 1_pInt) &
call IO_error(150_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
Nchunks_TwinFamilies = positions(1) - 1_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
do j = 1_pInt, Nchunks_TwinFamilies
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_Ntwin(j,instance) = IO_intValue(line,positions,1_pInt+j)
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
enddo
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
case ('ndot0','twinsize','twinburgers','r_twin')
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, Nchunks_TwinFamilies
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
tempPerTwin(j) = IO_floatValue(line,positions,1_pInt+j)
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
enddo
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
select case(tag)
case ('ndot0')
if (lattice_structure(phase) == LATTICE_fcc_ID) &
call IO_warning(42_pInt,ext_msg=trim(tag)//' for fcc ('//PLASTICITY_DISLOTWIN_label//')')
constitutive_dislotwin_Ndot0PerTwinFamily(1:Nchunks_TwinFamilies,instance) = tempPerTwin(1:Nchunks_TwinFamilies)
case ('twinsize')
constitutive_dislotwin_twinsizePerTwinFamily(1:Nchunks_TwinFamilies,instance) = tempPerTwin(1:Nchunks_TwinFamilies)
case ('twinburgers')
constitutive_dislotwin_burgersPerTwinFamily(1:Nchunks_TwinFamilies,instance) = tempPerTwin(1:Nchunks_TwinFamilies)
case ('r_twin')
constitutive_dislotwin_rPerTwinFamily(1:Nchunks_TwinFamilies,instance) = tempPerTwin(1:Nchunks_TwinFamilies)
end select
some minor polishing on the way
2014-04-04 20:10:30 +05:30
!--------------------------------------------------------------------------------------------------
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
! parameters depending on number of transformation system families
case ('ntrans')
if (positions(1) < Nchunks_TransFamilies + 1_pInt) &
call IO_warning(53_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
if (positions(1) > Nchunks_TransFamilies + 1_pInt) &
call IO_error(150_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
Nchunks_TransFamilies = positions(1) - 1_pInt
do j = 1_pInt, Nchunks_TransFamilies
constitutive_dislotwin_Ntrans(j,instance) = IO_intValue(line,positions,1_pInt+j)
enddo
!--------------------------------------------------------------------------------------------------
some minor polishing on the way
2014-04-04 20:10:30 +05:30
! parameters depending on number of interactions
case ('interaction_slipslip','interactionslipslip')
if (positions(1) < 1_pInt + Nchunks_SlipSlip) &
call IO_warning(52_pInt,ext_msg=trim(tag)//' ('//PLASTICITY_DISLOTWIN_label//')')
do j = 1_pInt, Nchunks_SlipSlip
constitutive_dislotwin_interaction_SlipSlip(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_DISLOTWIN_label//')')
do j = 1_pInt, Nchunks_SlipTwin
constitutive_dislotwin_interaction_SlipTwin(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_DISLOTWIN_label//')')
do j = 1_pInt, Nchunks_TwinSlip
constitutive_dislotwin_interaction_TwinSlip(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_DISLOTWIN_label//')')
do j = 1_pInt, Nchunks_TwinTwin
constitutive_dislotwin_interaction_TwinTwin(j,instance) = IO_floatValue(line,positions,1_pInt+j)
enddo
!--------------------------------------------------------------------------------------------------
! parameters independent of number of slip/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
case ('grainsize')
constitutive_dislotwin_GrainSize(instance) = IO_floatValue(line,positions,2_pInt)
case ('maxtwinfraction')
constitutive_dislotwin_MaxTwinFraction(instance) = IO_floatValue(line,positions,2_pInt)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
case ('p_shearband')
constitutive_dislotwin_pShearBand(instance) = IO_floatValue(line,positions,2_pInt)
case ('q_shearband')
constitutive_dislotwin_qShearBand(instance) = IO_floatValue(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 ('d0')
constitutive_dislotwin_D0(instance) = IO_floatValue(line,positions,2_pInt)
case ('qsd')
constitutive_dislotwin_Qsd(instance) = IO_floatValue(line,positions,2_pInt)
case ('atol_rho')
constitutive_dislotwin_aTolRho(instance) = IO_floatValue(line,positions,2_pInt)
case ('atol_twinfrac')
constitutive_dislotwin_aTolTwinFrac(instance) = IO_floatValue(line,positions,2_pInt)
case ('cmfptwin')
constitutive_dislotwin_Cmfptwin(instance) = IO_floatValue(line,positions,2_pInt)
case ('cthresholdtwin')
constitutive_dislotwin_Cthresholdtwin(instance) = IO_floatValue(line,positions,2_pInt)
case ('solidsolutionstrength')
constitutive_dislotwin_SolidSolutionStrength(instance) = IO_floatValue(line,positions,2_pInt)
case ('l0')
constitutive_dislotwin_L0(instance) = IO_floatValue(line,positions,2_pInt)
case ('xc')
constitutive_dislotwin_xc(instance) = IO_floatValue(line,positions,2_pInt)
case ('vcrossslip')
constitutive_dislotwin_VcrossSlip(instance) = IO_floatValue(line,positions,2_pInt)
case ('cedgedipmindistance')
constitutive_dislotwin_CEdgeDipMinDistance(instance) = IO_floatValue(line,positions,2_pInt)
case ('catomicvolume')
constitutive_dislotwin_CAtomicVolume(instance) = IO_floatValue(line,positions,2_pInt)
case ('sfe_0k')
constitutive_dislotwin_SFE_0K(instance) = IO_floatValue(line,positions,2_pInt)
case ('dsfe_dt')
constitutive_dislotwin_dSFE_dT(instance) = IO_floatValue(line,positions,2_pInt)
added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given
2014-03-14 05:20:55 +05:30
case ('dipoleformationfactor')
constitutive_dislotwin_dipoleFormationFactor(instance) = IO_floatValue(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 ('shearbandresistance')
constitutive_dislotwin_sbResistance(instance) = IO_floatValue(line,positions,2_pInt)
case ('shearbandvelocity')
constitutive_dislotwin_sbVelocity(instance) = IO_floatValue(line,positions,2_pInt)
case ('qedgepersbsystem')
constitutive_dislotwin_sbQedge(instance) = IO_floatValue(line,positions,2_pInt)
Renaming transformation input variables to more logical names
2014-09-23 21:08:41 +05:30
case ('cnuc')
constitutive_dislotwin_Cnuc(instance) = IO_floatValue(line,positions,2_pInt)
case ('cdwp')
constitutive_dislotwin_Cdwp(instance) = IO_floatValue(line,positions,2_pInt)
case ('cgro')
constitutive_dislotwin_Cgro(instance) = IO_floatValue(line,positions,2_pInt)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
case ('deltag')
constitutive_dislotwin_deltaG(instance) = IO_floatValue(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
endif; endif
enddo parsingFile
sanityChecks: do phase = 1_pInt, size(phase_plasticity)
myPhase: if (phase_plasticity(phase) == PLASTICITY_dislotwin_ID) then
instance = phase_plasticityInstance(phase)
Projection matrices for transformation
2014-09-24 18:01:16 +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
if (sum(constitutive_dislotwin_Nslip(:,instance)) < 0_pInt) &
call IO_error(211_pInt,el=instance,ext_msg='Nslip ('//PLASTICITY_DISLOTWIN_label//')')
if (sum(constitutive_dislotwin_Ntwin(:,instance)) < 0_pInt) &
call IO_error(211_pInt,el=instance,ext_msg='Ntwin ('//PLASTICITY_DISLOTWIN_label//')')
do f = 1_pInt,lattice_maxNslipFamily
if (constitutive_dislotwin_Nslip(f,instance) > 0_pInt) then
if (constitutive_dislotwin_rhoEdge0(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='rhoEdge0 ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_rhoEdgeDip0(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='rhoEdgeDip0 ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_burgersPerSlipFamily(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='slipBurgers ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_v0PerSlipFamily(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='v0 ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='tau_peierls ('//PLASTICITY_DISLOTWIN_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
endif
enddo
do f = 1_pInt,lattice_maxNtwinFamily
if (constitutive_dislotwin_Ntwin(f,instance) > 0_pInt) then
if (constitutive_dislotwin_burgersPerTwinFamily(f,instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='twinburgers ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_Ndot0PerTwinFamily(f,instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='ndot0 ('//PLASTICITY_DISLOTWIN_label//')')
endif
enddo
if (constitutive_dislotwin_CAtomicVolume(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='cAtomicVolume ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_D0(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='D0 ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_Qsd(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='Qsd ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (sum(constitutive_dislotwin_Ntwin(:,instance)) > 0_pInt) then
if (constitutive_dislotwin_SFE_0K(instance) == 0.0_pReal .and. &
constitutive_dislotwin_dSFE_dT(instance) == 0.0_pReal .and. &
lattice_structure(phase) == LATTICE_fcc_ID) &
call IO_error(211_pInt,el=instance,ext_msg='SFE0K ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_aTolRho(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='aTolRho ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_aTolTwinFrac(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='aTolTwinFrac ('//PLASTICITY_DISLOTWIN_label//')')
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 (constitutive_dislotwin_sbResistance(instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='sbResistance ('//PLASTICITY_DISLOTWIN_label//')')
if (constitutive_dislotwin_sbVelocity(instance) < 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='sbVelocity ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (constitutive_dislotwin_sbVelocity(instance) > 0.0_pReal .and. &
constitutive_dislotwin_pShearBand(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='pShearBand ('//PLASTICITY_DISLOTWIN_label//')')
added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given
2014-03-14 05:20:55 +05:30
if (constitutive_dislotwin_dipoleFormationFactor(instance) /= 0.0_pReal .and. &
constitutive_dislotwin_dipoleFormationFactor(instance) /= 1.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='dipoleFormationFactor ('//PLASTICITY_DISLOTWIN_label//')')
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (constitutive_dislotwin_sbVelocity(instance) > 0.0_pReal .and. &
constitutive_dislotwin_qShearBand(instance) <= 0.0_pReal) &
call IO_error(211_pInt,el=instance,ext_msg='qShearBand ('//PLASTICITY_DISLOTWIN_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
!--------------------------------------------------------------------------------------------------
! Determine total number of active slip or twin systems
constitutive_dislotwin_Nslip(:,instance) = min(lattice_NslipSystem(:,phase),constitutive_dislotwin_Nslip(:,instance))
constitutive_dislotwin_Ntwin(:,instance) = min(lattice_NtwinSystem(:,phase),constitutive_dislotwin_Ntwin(:,instance))
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
constitutive_dislotwin_Ntrans(:,instance)= min(lattice_NtransSystem(:,phase),constitutive_dislotwin_Ntrans(:,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_dislotwin_totalNslip(instance) = sum(constitutive_dislotwin_Nslip(:,instance))
constitutive_dislotwin_totalNtwin(instance) = sum(constitutive_dislotwin_Ntwin(:,instance))
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
constitutive_dislotwin_totalNtrans(instance) = sum(constitutive_dislotwin_Ntrans(:,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 myPhase
enddo sanityChecks
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!--------------------------------------------------------------------------------------------------
! allocation of variables whose size depends on the total number of active slip systems
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
maxTotalNslip = maxval(constitutive_dislotwin_totalNslip)
maxTotalNtwin = maxval(constitutive_dislotwin_totalNtwin)
maxTotalNtrans = maxval(constitutive_dislotwin_totalNtrans)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_burgersPerSlipSystem(maxTotalNslip, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_burgersPerTwinSystem(maxTotalNtwin, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_QedgePerSlipSystem(maxTotalNslip, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_v0PerSlipSystem(maxTotalNslip, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Ndot0PerTwinSystem(maxTotalNtwin, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_tau_r(maxTotalNtwin, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_twinsizePerTwinSystem(maxTotalNtwin, maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_CLambdaSlipPerSlipSystem(maxTotalNslip, maxNinstance),source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(constitutive_dislotwin_interactionMatrix_SlipSlip(maxval(constitutive_dislotwin_totalNslip),& ! slip resistance from slip activity
maxval(constitutive_dislotwin_totalNslip),&
maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_interactionMatrix_SlipTwin(maxval(constitutive_dislotwin_totalNslip),& ! slip resistance from twin activity
maxval(constitutive_dislotwin_totalNtwin),&
maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_interactionMatrix_TwinSlip(maxval(constitutive_dislotwin_totalNtwin),& ! twin resistance from slip activity
maxval(constitutive_dislotwin_totalNslip),&
maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_interactionMatrix_TwinTwin(maxval(constitutive_dislotwin_totalNtwin),& ! twin resistance from twin activity
maxval(constitutive_dislotwin_totalNtwin),&
maxNinstance), source=0.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
allocate(constitutive_dislotwin_forestProjectionEdge(maxTotalNslip,maxTotalNslip,maxNinstance), &
source=0.0_pReal)
Projection matrices for transformation
2014-09-24 18:01:16 +05:30
allocate(constitutive_dislotwin_projectionMatrix_Trans(maxTotalNtrans,maxTotalNtrans,maxNinstance), &
source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
allocate(constitutive_dislotwin_Ctwin66(6,6,maxTotalNtwin,maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Ctwin3333(3,3,3,3,maxTotalNtwin,maxNinstance), source=0.0_pReal)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
allocate(constitutive_dislotwin_Ctrans66(6,6,maxTotalNtrans,maxNinstance), source=0.0_pReal)
allocate(constitutive_dislotwin_Ctrans3333(3,3,3,3,maxTotalNtrans,maxNinstance), source=0.0_pReal)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25: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
initializeInstances: do phase = 1_pInt, size(phase_plasticity)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
myPhase2: if (phase_plasticity(phase) == PLASTICITY_dislotwin_ID) then
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
NofMyPhase=count(material_phase==phase)
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 = phase_plasticityInstance(phase)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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
ns = constitutive_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
nr = constitutive_dislotwin_totalNtrans(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
!--------------------------------------------------------------------------------------------------
! Determine size of postResults array
outputsLoop: do o = 1_pInt,constitutive_dislotwin_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_dislotwin_outputID(o,instance))
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case(edge_density_ID, &
dipole_density_ID, &
shear_rate_slip_ID, &
accumulated_shear_slip_ID, &
mfp_slip_ID, &
resolved_stress_slip_ID, &
threshold_stress_slip_ID, &
edge_dipole_distance_ID, &
stress_exponent_ID &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
)
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 = ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case(twin_fraction_ID, &
shear_rate_twin_ID, &
accumulated_shear_twin_ID, &
mfp_twin_ID, &
resolved_stress_twin_ID, &
threshold_stress_twin_ID &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
)
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 = nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case(resolved_stress_shearband_ID, &
shear_rate_shearband_ID &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
)
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 = 6_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case(sb_eigenvalues_ID)
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 = 3_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case(sb_eigenvectors_ID)
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 = 9_pInt
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
case(stress_trans_fraction_ID, &
strain_trans_fraction_ID, &
trans_fraction_ID &
)
mySize = nr
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
end select
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
if (mySize > 0_pInt) then ! any meaningful output found
constitutive_dislotwin_sizePostResult(o,instance) = mySize
constitutive_dislotwin_sizePostResults(instance) = constitutive_dislotwin_sizePostResults(instance) + mySize
endif
enddo outputsLoop
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
!--------------------------------------------------------------------------------------------------
! allocate state arrays
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
sizeDotState = int(size(CONSTITUTIVE_DISLOTWIN_listBasicSlipStates),pInt) * ns &
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
+ int(size(CONSTITUTIVE_DISLOTWIN_listBasicTwinStates),pInt) * nt &
+ int(size(CONSTITUTIVE_DISLOTWIN_listBasicTransStates),pInt) * nr
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
sizeState = sizeDotState &
+ int(size(CONSTITUTIVE_DISLOTWIN_listDependentSlipStates),pInt) * ns &
+ int(size(CONSTITUTIVE_DISLOTWIN_listDependentTwinStates),pInt) * nt
plasticState(phase)%sizeState = sizeState
plasticState(phase)%sizeDotState = sizeDotState
Small fix for missing constitutive dislotwin output
2014-09-16 13:32:19 +05:30
plasticState(phase)%sizePostResults = constitutive_dislotwin_sizePostResults(instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
allocate(plasticState(phase)%aTolState (sizeState), source=0.0_pReal)
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)
updated to read in TRIP parameters
2014-07-22 13:13:03 +05:30
allocate(plasticState(phase)%deltaState (sizeDotState,NofMyPhase), source=0.0_pReal)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
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)
Projection matrices for transformation
2014-09-24 18:01:16 +05:30
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
!* Process slip related parameters ------------------------------------------------
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
slipFamiliesLoop: do f = 1_pInt,lattice_maxNslipFamily
index_myFamily = sum(constitutive_dislotwin_Nslip(1:f-1_pInt,instance)) ! index in truncated slip system list
slipSystemsLoop: do j = 1_pInt,constitutive_dislotwin_Nslip(f,instance)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Burgers vector,
! dislocation velocity prefactor,
! mean free path prefactor,
! and minimum dipole distance
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_dislotwin_burgersPerSlipSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_burgersPerSlipFamily(f,instance)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
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_dislotwin_QedgePerSlipSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_QedgePerSlipFamily(f,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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_dislotwin_v0PerSlipSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_v0PerSlipFamily(f,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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_dislotwin_CLambdaSlipPerSlipSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_CLambdaSlipPerSlipFamily(f,instance)
!* Calculation of forest projections for edge dislocations
!* Interaction matrices
do o = 1_pInt,lattice_maxNslipFamily
index_otherFamily = sum(constitutive_dislotwin_Nslip(1:o-1_pInt,instance))
do k = 1_pInt,constitutive_dislotwin_Nslip(o,instance) ! loop over (active) systems in other family (slip)
constitutive_dislotwin_forestProjectionEdge(index_myFamily+j,index_otherFamily+k,instance) = &
abs(math_mul3x3(lattice_sn(:,sum(lattice_NslipSystem(1:f-1,phase))+j,phase), &
lattice_st(:,sum(lattice_NslipSystem(1:o-1,phase))+k,phase)))
constitutive_dislotwin_interactionMatrix_SlipSlip(index_myFamily+j,index_otherFamily+k,instance) = &
constitutive_dislotwin_interaction_SlipSlip(lattice_interactionSlipSlip( &
sum(lattice_NslipSystem(1:f-1,phase))+j, &
sum(lattice_NslipSystem(1:o-1,phase))+k, &
phase), instance )
enddo; enddo
do o = 1_pInt,lattice_maxNtwinFamily
index_otherFamily = sum(constitutive_dislotwin_Ntwin(1:o-1_pInt,instance))
do k = 1_pInt,constitutive_dislotwin_Ntwin(o,instance) ! loop over (active) systems in other family (twin)
constitutive_dislotwin_interactionMatrix_SlipTwin(index_myFamily+j,index_otherFamily+k,instance) = &
constitutive_dislotwin_interaction_SlipTwin(lattice_interactionSlipTwin( &
sum(lattice_NslipSystem(1:f-1_pInt,phase))+j, &
sum(lattice_NtwinSystem(1:o-1_pInt,phase))+k, &
phase), instance )
enddo; enddo
enddo slipSystemsLoop
enddo slipFamiliesLoop
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Process twin related parameters ------------------------------------------------
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
twinFamiliesLoop: do f = 1_pInt,lattice_maxNtwinFamily
index_myFamily = sum(constitutive_dislotwin_Ntwin(1:f-1_pInt,instance)) ! index in truncated twin system list
twinSystemsLoop: do j = 1_pInt,constitutive_dislotwin_Ntwin(f,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Burgers vector,
! nucleation rate prefactor,
! and twin size
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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_dislotwin_burgersPerTwinSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_burgersPerTwinFamily(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_dislotwin_Ndot0PerTwinSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_Ndot0PerTwinFamily(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_dislotwin_twinsizePerTwinSystem(index_myFamily+j,instance) = &
constitutive_dislotwin_twinsizePerTwinFamily(f,instance)
!* Rotate twin elasticity matrices
index_otherFamily = sum(lattice_NtwinSystem(1:f-1_pInt,phase)) ! index in full lattice twin list
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
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_dislotwin_Ctwin3333(l,m,n,o,index_myFamily+j,instance) = &
constitutive_dislotwin_Ctwin3333(l,m,n,o,index_myFamily+j,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
lattice_C3333(p,q,r,s,instance) * &
lattice_Qtwin(l,p,index_otherFamily+j,phase) * &
lattice_Qtwin(m,q,index_otherFamily+j,phase) * &
lattice_Qtwin(n,r,index_otherFamily+j,phase) * &
lattice_Qtwin(o,s,index_otherFamily+j,phase)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
enddo; enddo; enddo; enddo
enddo; enddo; enddo; enddo
constitutive_dislotwin_Ctwin66(1:6,1:6,index_myFamily+j,instance) = &
math_Mandel3333to66(constitutive_dislotwin_Ctwin3333(1:3,1:3,1:3,1:3,index_myFamily+j,instance))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Interaction matrices
do o = 1_pInt,lattice_maxNslipFamily
index_otherFamily = sum(constitutive_dislotwin_Nslip(1:o-1_pInt,instance))
do k = 1_pInt,constitutive_dislotwin_Nslip(o,instance) ! loop over (active) systems in other family (slip)
constitutive_dislotwin_interactionMatrix_TwinSlip(index_myFamily+j,index_otherFamily+k,instance) = &
constitutive_dislotwin_interaction_TwinSlip(lattice_interactionTwinSlip( &
sum(lattice_NtwinSystem(1:f-1_pInt,phase))+j, &
sum(lattice_NslipSystem(1:o-1_pInt,phase))+k, &
phase), instance )
enddo; enddo
do o = 1_pInt,lattice_maxNtwinFamily
index_otherFamily = sum(constitutive_dislotwin_Ntwin(1:o-1_pInt,instance))
do k = 1_pInt,constitutive_dislotwin_Ntwin(o,instance) ! loop over (active) systems in other family (twin)
constitutive_dislotwin_interactionMatrix_TwinTwin(index_myFamily+j,index_otherFamily+k,instance) = &
constitutive_dislotwin_interaction_TwinTwin(lattice_interactionTwinTwin( &
sum(lattice_NtwinSystem(1:f-1_pInt,phase))+j, &
sum(lattice_NtwinSystem(1:o-1_pInt,phase))+k, &
phase), instance )
enddo; enddo
enddo twinSystemsLoop
enddo twinFamiliesLoop
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* Process transformation related parameters ------------------------------------------------
Projection matrices for transformation
2014-09-24 18:01:16 +05:30
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
transFamiliesLoop: do f = 1_pInt,lattice_maxNtransFamily
index_myFamily = sum(constitutive_dislotwin_Ntrans(1:f-1_pInt,instance)) ! index in truncated trans system list
transSystemsLoop: do j = 1_pInt,constitutive_dislotwin_Ntrans(f,instance)
!* Rotate trans elasticity matrices
index_otherFamily = sum(lattice_NtransSystem(1:f-1_pInt,phase)) ! index in full lattice trans list
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_dislotwin_Ctrans3333(l,m,n,o,index_myFamily+j,instance) = &
constitutive_dislotwin_Ctrans3333(l,m,n,o,index_myFamily+j,instance) + &
lattice_C3333(p,q,r,s,instance) * &
lattice_Qtrans(l,p,index_otherFamily+j,phase) * &
lattice_Qtrans(m,q,index_otherFamily+j,phase) * &
lattice_Qtrans(n,r,index_otherFamily+j,phase) * &
lattice_Qtrans(o,s,index_otherFamily+j,phase)
enddo; enddo; enddo; enddo
enddo; enddo; enddo; enddo
constitutive_dislotwin_Ctrans66(1:6,1:6,index_myFamily+j,instance) = &
math_Mandel3333to66(constitutive_dislotwin_Ctrans3333(1:3,1:3,1:3,1:3,index_myFamily+j,instance))
Projection matrices for transformation
2014-09-24 18:01:16 +05:30
!* Projection matrices for shear from slip systems to fault-band (twin) systems for strain-induced martensite nucleation
do o = 1_pInt,lattice_maxNtransFamily
index_otherFamily = sum(constitutive_dislotwin_Ntrans(1:o-1_pInt,instance))
do k = 1_pInt,constitutive_dislotwin_Ntrans(o,instance) ! loop over (active) systems in other family (trans)
constitutive_dislotwin_projectionMatrix_Trans(index_myFamily+j,index_otherFamily+k,instance) = &
lattice_projectionTrans( sum(lattice_NtransSystem(1:f-1,phase))+j, &
sum(lattice_NtransSystem(1:o-1,phase))+k, phase)
enddo; enddo
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
enddo transSystemsLoop
enddo transFamiliesLoop
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
call constitutive_dislotwin_stateInit(phase,instance)
call constitutive_dislotwin_aTolState(phase,instance)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
endif myPhase2
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
enddo initializeInstances
replaced missued family index by system index
2013-02-06 14:15:08 +05:30
end subroutine constitutive_dislotwin_init
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
!> @brief sets the relevant state values for a given instance of this plasticity
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_dislotwin_stateInit(ph,instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
use math, only: &
pi
use lattice, only: &
lattice_maxNslipFamily, &
lattice_structure, &
lattice_mu, &
lattice_bcc_ID
use material, only: &
plasticState
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
implicit none
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
integer(pInt), intent(in) :: &
instance, & !< number specifying the instance of the plasticity
ph
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
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
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
integer(pInt) :: i,j,f,ns,nt,nr, index_myFamily
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
real(pReal), dimension(constitutive_dislotwin_totalNslip(instance)) :: &
rhoEdge0, &
rhoEdgeDip0, &
invLambdaSlip0, &
MeanFreePathSlip0, &
tauSlipThreshold0
real(pReal), dimension(constitutive_dislotwin_totalNtwin(instance)) :: &
MeanFreePathTwin0,TwinVolume0
tempState = 0.0_pReal
ns = constitutive_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
!--------------------------------------------------------------------------------------------------
! initialize basic slip state variables
do f = 1_pInt,lattice_maxNslipFamily
index_myFamily = sum(constitutive_dislotwin_Nslip(1:f-1_pInt,instance)) ! index in truncated slip system list
rhoEdge0(index_myFamily+1_pInt: &
index_myFamily+constitutive_dislotwin_Nslip(f,instance)) = &
constitutive_dislotwin_rhoEdge0(f,instance)
rhoEdgeDip0(index_myFamily+1_pInt: &
index_myFamily+constitutive_dislotwin_Nslip(f,instance)) = &
constitutive_dislotwin_rhoEdgeDip0(f,instance)
enddo
tempState(1_pInt:ns) = rhoEdge0
tempState(ns+1_pInt:2_pInt*ns) = rhoEdgeDip0
!--------------------------------------------------------------------------------------------------
! initialize dependent slip microstructural variables
forall (i = 1_pInt:ns) &
invLambdaSlip0(i) = sqrt(dot_product((rhoEdge0+rhoEdgeDip0),constitutive_dislotwin_forestProjectionEdge(1:ns,i,instance)))/ &
constitutive_dislotwin_CLambdaSlipPerSlipSystem(i,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
tempState(3_pInt*ns+2_pInt*nt+2_pInt*nr+1:4_pInt*ns+2_pInt*nt+2_pInt*nr) = invLambdaSlip0
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
forall (i = 1_pInt:ns) &
MeanFreePathSlip0(i) = &
constitutive_dislotwin_GrainSize(instance)/(1.0_pReal+invLambdaSlip0(i)*constitutive_dislotwin_GrainSize(instance))
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
tempState(5_pInt*ns+3_pInt*nt+2_pInt*nr+1:6_pInt*ns+3_pInt*nt+2_pInt*nr) = MeanFreePathSlip0
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
forall (i = 1_pInt:ns) &
tauSlipThreshold0(i) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_mu(ph)*constitutive_dislotwin_burgersPerSlipSystem(i,instance) * &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
sqrt(dot_product((rhoEdge0+rhoEdgeDip0),constitutive_dislotwin_interactionMatrix_SlipSlip(i,1:ns,instance)))
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
tempState(6_pInt*ns+4_pInt*nt+2_pInt*nr+1:7_pInt*ns+4_pInt*nt+2_pInt*nr) = tauSlipThreshold0
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
!--------------------------------------------------------------------------------------------------
! initialize dependent twin microstructural variables
forall (j = 1_pInt:nt) &
MeanFreePathTwin0(j) = constitutive_dislotwin_GrainSize(instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
tempState(6_pInt*ns+3_pInt*nt+2_pInt*nr+1_pInt:6_pInt*ns+4_pInt*nt+2_pInt*nr) = MeanFreePathTwin0
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
forall (j = 1_pInt:nt) &
TwinVolume0(j) = &
(pi/4.0_pReal)*constitutive_dislotwin_twinsizePerTwinSystem(j,instance)*MeanFreePathTwin0(j)**(2.0_pReal)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
tempState(7_pInt*ns+5_pInt*nt+2_pInt*nr+1_pInt:7_pInt*ns+6_pInt*nt+2_pInt*nr) = TwinVolume0
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
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,:)))
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
end subroutine constitutive_dislotwin_stateInit
!--------------------------------------------------------------------------------------------------
!> @brief sets the relevant state values for a given instance of this plasticity
!--------------------------------------------------------------------------------------------------
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
subroutine constitutive_dislotwin_aTolState(ph,instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
use material, only: &
plasticState
implicit none
integer(pInt), intent(in) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
ph, &
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
instance ! number specifying the current instance of the plasticity
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
integer(pInt) :: ns, nt, nr
ns = constitutive_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
nr = constitutive_dislotwin_totalNtrans(instance)
Adding accumulated shear due to "slip" and "twin" as two outputs
2013-05-20 02:02:06 +05:30
! Tolerance state for dislocation densities
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
plasticState(ph)%aTolState(1_pInt: &
2_pInt*ns) = constitutive_dislotwin_aTolRho(instance)
implementation of twin nucleation criteria according to Davids PhD thesis applies for fcc crystal structure only
2013-07-01 18:36:01 +05:30
Adding accumulated shear due to "slip" and "twin" as two outputs
2013-05-20 02:02:06 +05:30
! Tolerance state for accumulated shear due to slip
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
plasticState(ph)%aTolState(2_pInt*ns+1_pInt: &
3_pInt*ns)=1.0e6_pReal
Adding accumulated shear due to "slip" and "twin" as two outputs
2013-05-20 02:02:06 +05:30
! Tolerance state for twin volume fraction
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
plasticState(ph)%aTolState(3_pInt*ns+1_pInt: &
3_pInt*ns+nt) = constitutive_dislotwin_aTolTwinFrac(instance)
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
! Tolerance state for accumulated shear due to twin
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
plasticState(ph)%aTolState(3_pInt*ns+nt+1_pInt: &
3_pInt*ns+2_pInt*nt) = 1.0e6_pReal
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
! Tolerance state for stress-assisted martensite volume fraction
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
plasticState(ph)%aTolState(3_pInt*ns+2_pInt*nt+1_pInt: &
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
3_pInt*ns+2_pInt*nt+nr) = 1.0e-6
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
! Tolerance state for strain-induced martensite volume fraction
plasticState(ph)%aTolState(3_pInt*ns+2_pInt*nt+nr+1_pInt: &
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
3_pInt*ns+2_pInt*nt+2_pInt*nr) = 1.0e-6
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
end subroutine constitutive_dislotwin_aTolState
updated to read in TRIP parameters
2014-07-22 13:13:03 +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_dislotwin_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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
implicit none
real(pReal), dimension(6,6) :: &
constitutive_dislotwin_homogenizedC
integer(pInt), intent(in) :: &
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
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
integer(pInt) :: instance,ns,nt,nr,i, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
ph, &
of
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
real(pReal) :: sumf, sumftr
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_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
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 to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sumf = sum(plasticState(ph)%state((3_pInt*ns+1_pInt):(3_pInt*ns+nt),of)) ! safe for nt == 0
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* Total transformed volume fraction
sumftr = sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of)) + &
sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of))
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
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
constitutive_dislotwin_homogenizedC = (1.0_pReal-sumf-sumftr)*lattice_C66(1:6,1:6,ph)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
do i=1_pInt,nt
bufgix for last commit, Franz please rather check this
2014-08-08 17:46:51 +05:30
constitutive_dislotwin_homogenizedC = constitutive_dislotwin_homogenizedC &
+ plasticState(ph)%state(3_pInt*ns+i, of)*constitutive_dislotwin_Ctwin66(1:6,1:6,i,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
do i=1_pInt,nr
constitutive_dislotwin_homogenizedC = constitutive_dislotwin_homogenizedC &
+ (plasticState(ph)%state(3_pInt*ns+2_pInt*nt+i, of) + plasticState(ph)%state(3_pInt*ns+2_pInt*nt+nr+i, of))*&
constitutive_dislotwin_Ctrans66(1:6,1:6,i,instance)
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
end function constitutive_dislotwin_homogenizedC
!--------------------------------------------------------------------------------------------------
!> @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_dislotwin_microstructure(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 math, only: &
pi
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: &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
lattice_structure, &
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_mu, &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
lattice_nu, &
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
lattice_bcc_ID, &
lattice_maxNslipFamily
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
implicit none
integer(pInt), intent(in) :: &
ipc, & !< component-ID of integration point
ip, & !< integration point
el !< element
real(pReal), intent(in) :: &
temperature !< temperature at IP
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
integer(pInt) :: &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
instance, &
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
ns,nt,nr,s,t, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
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) :: &
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
sumf,sfe,x0,sumftr
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(constitutive_dislotwin_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: fOverStacksize
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_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
!BASIC STATES
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* State: 1 : ns rho_edge
!* State: ns+1 : 2*ns rho_dipole
!* State: 2*ns+1 : 3*ns accumulated shear due to slip
!* State: 3*ns+1 : 3*ns+nt f
!* State: 3*ns+nt+1 : 3*ns+2*nt accumulated shear due to twin
!* State: 3*ns+2*nt+1 : 3*ns+2*nt+nr stress-assisted martensite volume fraction
!* State: 3*ns+2*nt+nr+1 : 3*ns+2*nt+2*nr strain-induced martensite volume fraction
started to introduce state for phase transformation (adding one set of additional basic states)
2014-08-08 17:50:43 +05:30
!DEPENDENT STATES
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* State: 3*ns+2*nt+2*nr+1 : 4*ns+2*nt+2*nr 1/lambda_slip
!* State: 4*ns+2*nt+2*nr+1 : 5*ns+2*nt+2*nr 1/lambda_sliptwin
!* State: 5*ns+2*nt+2*nr+1 : 5*ns+3*nt+2*nr 1/lambda_twin
!* State: 5*ns+3*nt+2*nr+1 : 6*ns+3*nt+2*nr mfp_slip
!* State: 6*ns+3*nt+2*nr+1 : 6*ns+4*nt+2*nr mfp_twin
!* State: 6*ns+4*nt+2*nr+1 : 7*ns+4*nt+2*nr threshold_stress_slip
!* State: 7*ns+4*nt+2*nr+1 : 7*ns+5*nt+2*nr threshold_stress_twin
!* State: 7*ns+5*nt+2*nr+1 : 7*ns+6*nt+2*nr twin volume
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 to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sumf = sum(plasticState(ph)%state((3*ns+1):(3*ns+nt), of)) ! 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
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
!* Total transformed volume fraction
sumftr = sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of)) + &
sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* Stacking fault energy
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
sfe = constitutive_dislotwin_SFE_0K(instance) + &
constitutive_dislotwin_dSFE_dT(instance) * Temperature
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* rescaled twin volume fraction for topology
forall (t = 1_pInt:nt) &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
fOverStacksize(t) = &
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+t, of)/constitutive_dislotwin_twinsizePerTwinSystem(t,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* 1/mean free distance between 2 forest dislocations seen by a moving dislocation
forall (s = 1_pInt:ns) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(3_pInt*ns+2_pInt*nt+2_pInt*nr+s, of) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sqrt(dot_product((plasticState(ph)%state(1:ns,of)+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_dislotwin_forestProjectionEdge(1:ns,s,instance)))/ &
constitutive_dislotwin_CLambdaSlipPerSlipSystem(s,instance)
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* 1/mean free distance between 2 twin stacks from different systems seen by a moving dislocation
!$OMP CRITICAL (evilmatmul)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state((4_pInt*ns+2_pInt*nt+2_pInt*nr+1_pInt):(5_pInt*ns+2_pInt*nt+2_pInt*nr), of) = 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
if (nt > 0_pInt .and. ns > 0_pInt) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state((4_pInt*ns+2_pInt*nt+2_pInt*nr+1):(5_pInt*ns+2_pInt*nt+2_pInt*nr), of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
matmul(constitutive_dislotwin_interactionMatrix_SlipTwin(1:ns,1:nt,instance),fOverStacksize(1:nt))/(1.0_pReal-sumf)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!$OMP END CRITICAL (evilmatmul)
!* 1/mean free distance between 2 twin stacks from different systems seen by a growing twin
!$OMP CRITICAL (evilmatmul)
if (nt > 0_pInt) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state((5_pInt*ns+2_pInt*nt+2_pInt*nr+1_pInt):(5_pInt*ns+3_pInt*nt+2_pInt*nr), of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
matmul(constitutive_dislotwin_interactionMatrix_TwinTwin(1:nt,1:nt,instance),fOverStacksize(1:nt))/(1.0_pReal-sumf)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!$OMP END CRITICAL (evilmatmul)
!* mean free path between 2 obstacles seen by a moving dislocation
do s = 1_pInt,ns
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
if (nt > 0_pInt) then
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(5_pInt*ns+3_pInt*nt+2_pInt*nr+s, of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_dislotwin_GrainSize(instance)/(1.0_pReal+constitutive_dislotwin_GrainSize(instance)*&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
(plasticState(ph)%state(3_pInt*ns+2_pInt*nt+2_pInt*nr+s, of) + &
plasticState(ph)%state(4_pInt*ns+2_pInt*nt+2_pInt*nr+s, of)))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
else
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(5_pInt*ns+2_pInt*nr+s, of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_dislotwin_GrainSize(instance)/&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
(1.0_pReal+constitutive_dislotwin_GrainSize(instance)*(plasticState(ph)%state(3_pInt*ns+s, of)))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
endif
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25: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
!* mean free path between 2 obstacles seen by a growing twin
forall (t = 1_pInt:nt) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(6_pInt*ns+3_pInt*nt+2_pInt*nr+t, of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_dislotwin_Cmfptwin(instance)*constitutive_dislotwin_GrainSize(instance))/&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
(1.0_pReal+constitutive_dislotwin_GrainSize(instance)*plasticState(ph)%state(5_pInt*ns+2_pInt*nt+2_pInt*nr+t, of))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* threshold stress for dislocation motion
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
forall (s = 1_pInt:ns) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(6_pInt*ns+4_pInt*nt+2_pInt*nr+s, of) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_mu(ph)*constitutive_dislotwin_burgersPerSlipSystem(s,instance)*&
sqrt(dot_product((plasticState(ph)%state(1:ns, of)+plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of)),&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_interactionMatrix_SlipSlip(s,1:ns,instance)))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* threshold stress for growing twin
forall (t = 1_pInt:nt) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(7_pInt*ns+4_pInt*nt+2_pInt*nr+t, of) = &
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_dislotwin_Cthresholdtwin(instance)*&
(sfe/(3.0_pReal*constitutive_dislotwin_burgersPerTwinSystem(t,instance))+&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
3.0_pReal*constitutive_dislotwin_burgersPerTwinSystem(t,instance)*lattice_mu(ph)/&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_dislotwin_L0(instance)*constitutive_dislotwin_burgersPerSlipSystem(t,instance)))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* final twin volume after growth
forall (t = 1_pInt:nt) &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(7_pInt*ns+5_pInt*nt+2_pInt*nr+t, of) = &
(pi/4.0_pReal)*constitutive_dislotwin_twinsizePerTwinSystem(t,instance)*&
plasticState(ph)%state(6_pInt*ns+3_pInt*nt+2_pInt*nr+t, of)**(2.0_pReal)
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* equilibrium seperation of partial dislocations
do t = 1_pInt,nt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
x0 = lattice_mu(ph)*constitutive_dislotwin_burgersPerTwinSystem(t,instance)**(2.0_pReal)/&
(sfe*8.0_pReal*pi)*(2.0_pReal+lattice_nu(ph))/(1.0_pReal-lattice_nu(ph))
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
constitutive_dislotwin_tau_r(t,instance)= &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
lattice_mu(ph)*constitutive_dislotwin_burgersPerTwinSystem(t,instance)/(2.0_pReal*pi)*&
(1/(x0+constitutive_dislotwin_xc(instance))+cos(pi/3.0_pReal)/x0) !!! used where??
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
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
end subroutine constitutive_dislotwin_microstructure
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +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 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_dislotwin_LpAndItsTangent(Lp,dLp_dTstar,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 prec, only: &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
tol_math_check
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use math, only: &
math_Plain3333to99, &
math_Mandel6to33, &
math_Mandel33to6, &
math_spectralDecompositionSym33, &
math_tensorproduct, &
math_symmetric33, &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
math_mul33x3, &
math_norm33
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
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_Sslip, &
lattice_Sslip_v, &
lattice_Stwin, &
lattice_Stwin_v, &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
lattice_NItrans, &
lattice_NItrans_v, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
lattice_maxNslipFamily,&
lattice_maxNtwinFamily, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_maxNtransFamily, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
lattice_NslipSystem, &
lattice_NtwinSystem, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_NtransSystem, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
lattice_shearTwin, &
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_structure, &
lattice_fcc_twinNucleationSlipPair, &
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
LATTICE_fcc_ID
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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) :: ipc,ip,el
real(pReal), intent(in) :: Temperature
real(pReal), dimension(6), intent(in) :: Tstar_v
real(pReal), dimension(3,3), intent(out) :: Lp
real(pReal), dimension(9,9), intent(out) :: dLp_dTstar
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
integer(pInt) :: instance,ph,of,ns,nt,nr,f,i,j,k,l,m,n,index_myFamily,s1,s2
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
real(pReal) :: sumf,sumftr,StressRatio_p,StressRatio_pminus1,StressRatio_r,BoltzmannRatio,DotGamma0,Ndot0
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(3,3,3,3) :: dLp_dTstar3333
real(pReal), dimension(constitutive_dislotwin_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
gdot_slip,dgdot_dtauslip,tau_slip
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(constitutive_dislotwin_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
gdot_twin,dgdot_dtautwin,tau_twin
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
real(pReal), dimension(constitutive_dislotwin_totalNtrans(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
V_trans,Vdot_trans,dVdot_dtautrans,tau_trans,fstress_trans
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) :: gdot_sb,dgdot_dtausb,tau_sb
real(pReal), dimension(3,3) :: eigVectors, sb_Smatrix
real(pReal), dimension(3) :: eigValues, sb_s, sb_m
real(pReal), dimension(3,6), parameter :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
sb_sComposition = &
reshape(real([&
1, 0, 1, &
1, 0,-1, &
1, 1, 0, &
1,-1, 0, &
0, 1, 1, &
0, 1,-1 &
],pReal),[ 3,6]), &
sb_mComposition = &
reshape(real([&
1, 0,-1, &
1, 0,+1, &
1,-1, 0, &
1, 1, 0, &
0, 1,-1, &
0, 1, 1 &
],pReal),[ 3,6])
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
logical error
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_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +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 to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sumf = sum(plasticState(ph)%state((3_pInt*ns+1_pInt):(3_pInt*ns+nt), of)) ! safe for nt == 0
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* Total transformed volume fraction
sumftr = sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of)) + &
sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
Lp = 0.0_pReal
dLp_dTstar3333 = 0.0_pReal
dLp_dTstar = 0.0_pReal
!* Dislocation glide part
gdot_slip = 0.0_pReal
dgdot_dtauslip = 0.0_pReal
j = 0_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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
slipSystemsLoop: do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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))
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if((abs(tau_slip(j))-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of)) > tol_math_check) then
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_p = ((abs(tau_slip(j))- plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
**constitutive_dislotwin_pPerSlipFamily(f,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_pminus1 = ((abs(tau_slip(j))-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
**(constitutive_dislotwin_pPerSlipFamily(f,instance)-1.0_pReal)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Boltzmann ratio
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
BoltzmannRatio = constitutive_dislotwin_QedgePerSlipSystem(j,instance)/(kB*Temperature)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Initial shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
DotGamma0 = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)*&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_v0PerSlipSystem(j,instance)
!* Shear rates due to slip
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
gdot_slip(j) = (1.0_pReal - sumf - sumftr) * DotGamma0 &
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
* exp(-BoltzmannRatio*(1-StressRatio_p) ** constitutive_dislotwin_qPerSlipFamily(f,instance)) &
* sign(1.0_pReal,tau_slip(j))
!* Derivatives of shear rates
dgdot_dtauslip(j) = &
abs(gdot_slip(j))*BoltzmannRatio*constitutive_dislotwin_pPerSlipFamily(f,instance)&
*constitutive_dislotwin_qPerSlipFamily(f,instance)/&
(constitutive_dislotwin_SolidSolutionStrength(instance)+constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance))*&
StressRatio_pminus1*(1-StressRatio_p)**(constitutive_dislotwin_qPerSlipFamily(f,instance)-1.0_pReal)
endif
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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 + gdot_slip(j)*lattice_Sslip(:,:,1,index_myFamily+i,ph)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Calculation of the tangent of Lp
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
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)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
enddo slipSystemsLoop
enddo slipFamiliesLoop
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* Shear banding (shearband) part
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
if(constitutive_dislotwin_sbVelocity(instance) /= 0.0_pReal .and. &
constitutive_dislotwin_sbResistance(instance) /= 0.0_pReal) then
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
gdot_sb = 0.0_pReal
dgdot_dtausb = 0.0_pReal
call math_spectralDecompositionSym33(math_Mandel6to33(Tstar_v),eigValues,eigVectors, error)
do j = 1_pInt,6_pInt
sb_s = 0.5_pReal*sqrt(2.0_pReal)*math_mul33x3(eigVectors,sb_sComposition(1:3,j))
sb_m = 0.5_pReal*sqrt(2.0_pReal)*math_mul33x3(eigVectors,sb_mComposition(1:3,j))
sb_Smatrix = math_tensorproduct(sb_s,sb_m)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
constitutive_dislotwin_sbSv(1:6,j,ipc,ip,el) = math_Mandel33to6(math_symmetric33(sb_Smatrix))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Calculation of Lp
!* Resolved shear stress on shear banding system
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
tau_sb(j) = dot_product(Tstar_v,constitutive_dislotwin_sbSv(1:6,j,ipc,ip,el))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Stress ratios
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (abs(tau_sb(j)) < tol_math_check) then
StressRatio_p = 0.0_pReal
StressRatio_pminus1 = 0.0_pReal
else
StressRatio_p = (abs(tau_sb(j))/constitutive_dislotwin_sbResistance(instance))&
**constitutive_dislotwin_pShearBand(instance)
StressRatio_pminus1 = (abs(tau_sb(j))/constitutive_dislotwin_sbResistance(instance))&
**(constitutive_dislotwin_pShearBand(instance)-1.0_pReal)
endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Boltzmann ratio
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatio = constitutive_dislotwin_sbQedge(instance)/(kB*Temperature)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Initial shear rates
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
DotGamma0 = constitutive_dislotwin_sbVelocity(instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Shear rates due to shearband
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
gdot_sb(j) = DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)**&
constitutive_dislotwin_qShearBand(instance))*sign(1.0_pReal,tau_sb(j))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Derivatives of shear rates
dgdot_dtausb(j) = &
((abs(gdot_sb(j))*BoltzmannRatio*&
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_pShearBand(instance)*constitutive_dislotwin_qShearBand(instance))/&
constitutive_dislotwin_sbResistance(instance))*&
StressRatio_pminus1*(1_pInt-StressRatio_p)**(constitutive_dislotwin_qShearBand(instance)-1.0_pReal)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Plastic velocity gradient for shear banding
Lp = Lp + gdot_sb(j)*sb_Smatrix
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Calculation of the tangent of Lp
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
dLp_dTstar3333(k,l,m,n) = &
dLp_dTstar3333(k,l,m,n) + dgdot_dtausb(j)*&
sb_Smatrix(k,l)*&
sb_Smatrix(m,n)
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
end if
all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState).
2012-05-16 20:13:26 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* Mechanical twinning part
gdot_twin = 0.0_pReal
dgdot_dtautwin = 0.0_pReal
j = 0_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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
twinSystemsLoop: do i = 1_pInt,constitutive_dislotwin_Ntwin(f,instance)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
j = j+1_pInt
!* 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))
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Stress ratios
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (tau_twin(j) > tol_math_check) then
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_r = (plasticState(ph)%state(7*ns+4*nt+2*nr+j, of)/&
tau_twin(j))**constitutive_dislotwin_rPerTwinFamily(f,instance)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Shear rates and their derivatives due to twin
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
select case(lattice_structure(ph))
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (LATTICE_fcc_ID)
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
s1=lattice_fcc_twinNucleationSlipPair(1,index_myFamily+i)
s2=lattice_fcc_twinNucleationSlipPair(2,index_myFamily+i)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
if (tau_twin(j) < constitutive_dislotwin_tau_r(j,instance)) then
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Ndot0=(abs(gdot_slip(s1))*(plasticState(ph)%state(s2,of)+plasticState(ph)%state(ns+s2, of))+&
abs(gdot_slip(s2))*(plasticState(ph)%state(s1, of)+plasticState(ph)%state(ns+s1, of)))/&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_dislotwin_L0(instance)*constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
(1.0_pReal-exp(-constitutive_dislotwin_VcrossSlip(instance)/(kB*Temperature)*&
(constitutive_dislotwin_tau_r(j,instance)-tau_twin(j))))
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
else
Ndot0=0.0_pReal
end if
case default
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
Ndot0=constitutive_dislotwin_Ndot0PerTwinSystem(j,instance)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
end select
gdot_twin(j) = &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
(1.0_pReal-sumf-sumftr)*lattice_shearTwin(index_myFamily+i,ph)*&
plasticState(ph)%state(7*ns+5*nt+2*nr+j, of)*Ndot0*exp(-StressRatio_r)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
dgdot_dtautwin(j) = ((gdot_twin(j)*constitutive_dislotwin_rPerTwinFamily(f,instance))/tau_twin(j))*StressRatio_r
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
endif
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +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 + gdot_twin(j)*lattice_Stwin(:,:,index_myFamily+i,ph)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
!* Calculation of the tangent of Lp
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
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)
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
enddo twinSystemsLoop
enddo twinFamiliesLoop
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
!* Phase transformation part
Vdot_trans = 0.0_pReal
dVdot_dtautrans = 0.0_pReal
j = 0_pInt
transFamiliesLoop: do f = 1_pInt,lattice_maxNtransFamily
index_myFamily = sum(lattice_NtransSystem(1:f-1_pInt,ph)) ! at which index starts my family
transSystemsLoop: do i = 1_pInt,constitutive_dislotwin_Ntrans(f,instance)
j = j+1_pInt
!* Resolved shear stress on transformation system
tau_trans(j) = dot_product(Tstar_v,lattice_NItrans_v(:,index_myFamily+i,ph))
!* Total martensite volume fraction for transformation system
V_trans(j) = plasticState(ph)%state(3*ns+2*nt+j, of) + plasticState(ph)%state(3*ns+2*nt+nr+j, of)
!* Driving force for stress-assisted martensite growth on transformation system
fstress_trans(j) = tau_trans(j) - &
constitutive_dislotwin_Cdwp(instance)*(2*V_trans(j) - 6*V_trans(j)**2 + 4*V_trans(j)**3) - &
12*constitutive_dislotwin_deltaG(instance)*(V_trans(j)**2 - V_trans(j)**3)
!* Transformation rate of stress-assisted martensite and its derivative
if (fstress_trans(j) > tol_math_check) then
Vdot_trans(j) = constitutive_dislotwin_Cgro(instance)*(1.0_pReal - sumf - sumftr)*fstress_trans(j)
dVdot_dtautrans(j) = constitutive_dislotwin_Cgro(instance)*(1.0_pReal - sumf - sumftr)
endif
!* Plastic velocity gradient for phase transformation
Lp = Lp + Vdot_trans(j)*lattice_NItrans(:,:,index_myFamily+i,ph)
!* Calculation of the tangent of Lp
forall (k=1_pInt:3_pInt,l=1_pInt:3_pInt,m=1_pInt:3_pInt,n=1_pInt:3_pInt) &
dLp_dTstar3333(k,l,m,n) = &
dLp_dTstar3333(k,l,m,n) + dVdot_dtautrans(j)*&
lattice_NItrans(k,l,index_myFamily+i,ph)*&
lattice_NItrans(m,n,index_myFamily+i,ph)
enddo transSystemsLoop
enddo transFamiliesLoop
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
dLp_dTstar = math_Plain3333to99(dLp_dTstar3333)
end subroutine constitutive_dislotwin_LpAndItsTangent
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +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 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_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
use prec, only: &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
tol_math_check
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
use math, only: &
pi
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
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
use lattice, only: &
lattice_Sslip_v, &
lattice_Stwin_v, &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
lattice_NItrans_v, &
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
lattice_maxNslipFamily, &
lattice_maxNtwinFamily, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_maxNtransFamily, &
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
lattice_NslipSystem, &
lattice_NtwinSystem, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_NtransSystem, &
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
lattice_sheartwin, &
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_mu, &
lattice_structure, &
lattice_fcc_twinNucleationSlipPair, &
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
lattice_fcc_transNucleationTwinPair, &
lattice_fcc_shearCritTrans, &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
LATTICE_fcc_ID, &
LATTICE_bcc_ID
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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):: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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) :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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) :: &
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
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
integer(pInt) :: instance,ns,nt,nr,f,i,j,index_myFamily,s1,s2,a,b,sa,sb,ssa,ssb, &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
ph, &
of
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
real(pReal) :: sumf,sumftr,StressRatio_p,StressRatio_pminus1,BoltzmannRatio,DotGamma0,&
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
EdgeDipMinDistance,AtomicVolume,VacancyDiffusion,StressRatio_r,Ndot0
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(constitutive_dislotwin_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_slip,tau_slip,DotRhoMultiplication,EdgeDipDistance,DotRhoEdgeEdgeAnnihilation,DotRhoEdgeDipAnnihilation,&
ClimbVelocity,DotRhoEdgeDipClimb,DotRhoDipFormation
real(pReal), dimension(constitutive_dislotwin_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
tau_twin
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
real(pReal), dimension(constitutive_dislotwin_totalNtrans(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
V_trans,Vdot_trans,dVdot_dtautrans,tau_trans,fstress_trans,shear_trans,shearrate_trans,probrate_trans
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +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_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
finished updating to new state of dislokmc
2014-08-07 21:20:52 +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 to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sumf = sum(plasticState(ph)%state((3_pInt*ns+1_pInt):(3_pInt*ns+nt), of)) ! safe for nt == 0
plasticState(ph)%dotState(:,of) = 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
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
!* Total transformed volume fraction
sumftr = sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of)) + &
sum(plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
!* Dislocation density evolution
gdot_slip = 0.0_pReal
j = 0_pInt
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
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
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 i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in family
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
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
!* 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))
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if((abs(tau_slip(j))-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of)) > tol_math_check) then
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_p = ((abs(tau_slip(j))-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**constitutive_dislotwin_pPerSlipFamily(f,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_pminus1 = ((abs(tau_slip(j))-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**(constitutive_dislotwin_pPerSlipFamily(f,instance)-1.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
!* Boltzmann ratio
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
BoltzmannRatio = constitutive_dislotwin_QedgePerSlipSystem(j,instance)/(kB*Temperature)
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
!* Initial shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
DotGamma0 = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)*&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_v0PerSlipSystem(j,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Shear rates due to slip
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
gdot_slip(j) = DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)** &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_qPerSlipFamily(f,instance))*sign(1.0_pReal,tau_slip(j))
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +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
!* Multiplication
DotRhoMultiplication(j) = abs(gdot_slip(j))/&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(constitutive_dislotwin_burgersPerSlipSystem(j,instance)* &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(5*ns+3*nt+2*nr+j, of))
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
!* Dipole formation
EdgeDipMinDistance = &
constitutive_dislotwin_CEdgeDipMinDistance(instance)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)
if (tau_slip(j) == 0.0_pReal) then
DotRhoDipFormation(j) = 0.0_pReal
else
EdgeDipDistance(j) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(3.0_pReal*lattice_mu(ph)*constitutive_dislotwin_burgersPerSlipSystem(j,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
(16.0_pReal*pi*abs(tau_slip(j)))
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if (EdgeDipDistance(j)>plasticState(ph)%state(5*ns+3*nt+2*nr+j, of)) &
EdgeDipDistance(j)=plasticState(ph)%state(5*ns+3*nt+2*nr+j, of)
reintroduced dislocation hardening for bcc
2014-03-17 21:12:46 +05:30
if (EdgeDipDistance(j)<EdgeDipMinDistance) EdgeDipDistance(j)=EdgeDipMinDistance
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
DotRhoDipFormation(j) = &
((2.0_pReal*EdgeDipDistance(j))/constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*abs(gdot_slip(j))*constitutive_dislotwin_dipoleFormationFactor(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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Spontaneous annihilation of 2 single edge dislocations
DotRhoEdgeEdgeAnnihilation(j) = &
((2.0_pReal*EdgeDipMinDistance)/constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*abs(gdot_slip(j))
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
!* Spontaneous annihilation of a single edge dislocation with a dipole constituent
DotRhoEdgeDipAnnihilation(j) = &
((2.0_pReal*EdgeDipMinDistance)/constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
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)*abs(gdot_slip(j))
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
!* Dislocation dipole climb
AtomicVolume = &
constitutive_dislotwin_CAtomicVolume(instance)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)**(3.0_pReal)
VacancyDiffusion = &
constitutive_dislotwin_D0(instance)*exp(-constitutive_dislotwin_Qsd(instance)/(kB*Temperature))
if (tau_slip(j) == 0.0_pReal) then
DotRhoEdgeDipClimb(j) = 0.0_pReal
else
ClimbVelocity(j) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
((3.0_pReal*lattice_mu(ph)*VacancyDiffusion*AtomicVolume)/(2.0_pReal*pi*kB*Temperature))*&
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
(1/(EdgeDipDistance(j)+EdgeDipMinDistance))
DotRhoEdgeDipClimb(j) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(4.0_pReal*ClimbVelocity(j)*plasticState(ph)%state(ns+j, of))/(EdgeDipDistance(j)-EdgeDipMinDistance)
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
!* Edge dislocation density rate of change
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(j, of) = &
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
DotRhoMultiplication(j)-DotRhoDipFormation(j)-DotRhoEdgeEdgeAnnihilation(j)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Edge dislocation dipole density rate of change
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) = &
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
DotRhoDipFormation(j)-DotRhoEdgeDipAnnihilation(j)-DotRhoEdgeDipClimb(j)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
!* Dotstate for accumulated shear due to slip
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(2_pInt*ns+j, of) = gdot_slip(j)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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
enddo
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
!* Twin volume fraction evolution
j = 0_pInt
do f = 1_pInt,lattice_maxNtwinFamily ! loop over all twin families
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
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 i = 1_pInt,constitutive_dislotwin_Ntwin(f,instance) ! process each (active) twin system in family
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
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
!* 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))
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
!* Stress ratios
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (tau_twin(j) > tol_math_check) then
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_r = (plasticState(ph)%state(7*ns+4*nt+2*nr+j, of)/&
tau_twin(j))**constitutive_dislotwin_rPerTwinFamily(f,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
!* Shear rates and their derivatives due to twin
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
select case(lattice_structure(ph))
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 (LATTICE_fcc_ID)
s1=lattice_fcc_twinNucleationSlipPair(1,index_myFamily+i)
s2=lattice_fcc_twinNucleationSlipPair(2,index_myFamily+i)
if (tau_twin(j) < constitutive_dislotwin_tau_r(j,instance)) then
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Ndot0=(abs(gdot_slip(s1))*(plasticState(ph)%state(s2, of)+plasticState(ph)%state(ns+s2, of))+&
abs(gdot_slip(s2))*(plasticState(ph)%state(s1, of)+plasticState(ph)%state(ns+s1, of)))/&
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_dislotwin_L0(instance)*constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
(1.0_pReal-exp(-constitutive_dislotwin_VcrossSlip(instance)/(kB*Temperature)*&
(constitutive_dislotwin_tau_r(j,instance)-tau_twin(j))))
else
Ndot0=0.0_pReal
end if
case default
Ndot0=constitutive_dislotwin_Ndot0PerTwinSystem(j,instance)
end select
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(3_pInt*ns+j, of) = &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
(1.0_pReal-sumf-sumftr)*&
plasticState(ph)%state(7_pInt*ns+5_pInt*nt+2*nr+j, of)*Ndot0*exp(-StressRatio_r)
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
!* Dotstate for accumulated shear due to twin
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%dotState(3_pInt*ns+nt+j, of) = plasticState(ph)%dotState(3_pInt*ns+j, of) * &
lattice_sheartwin(index_myfamily+i,ph)
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
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
prepared new State for dislotwin
2014-06-11 17:41:14 +05:30
TRIP contribution to Lp and dotState
2014-10-01 13:41:39 +05:30
!* Transformation volume fraction evolution
shear_trans = 0.0_pReal
shearrate_trans = 0.0_pReal
probrate_trans = 0.0_pReal
j = 0_pInt
do f = 1_pInt,lattice_maxNtransFamily ! loop over all trans families
index_myFamily = sum(lattice_NtransSystem(1:f-1_pInt,ph)) ! at which index starts my family
!* Projection of shear and shear rate on fault band (twin) systems
shear_trans = matmul(constitutive_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%state(2_pInt*ns+1_pInt:3_pInt*ns, of))
shearrate_trans = matmul(constitutive_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%dotState(2_pInt*ns+1_pInt:3_pInt*ns, of))
do i = 1_pInt,constitutive_dislotwin_Ntrans(f,instance) ! process each (active) trans system in family
j = j+1_pInt
!* Resolved shear stress on transformation system
tau_trans(j) = dot_product(Tstar_v,lattice_NItrans_v(:,index_myFamily+i,ph))
!* Total martensite volume fraction for transformation system
V_trans(j) = plasticState(ph)%state(3*ns+2*nt+j, of) + plasticState(ph)%state(3*ns+2*nt+nr+j, of)
!* Driving force for stress-assisted martensite growth
fstress_trans(j) = tau_trans(j) - &
constitutive_dislotwin_Cdwp(instance)*(2*V_trans(j) - 6*V_trans(j)**2 + 4*V_trans(j)**3) - &
12*constitutive_dislotwin_deltaG(instance)*(V_trans(j)**2 - V_trans(j)**3)
!* Dotstate for stress-assisted martensite volume fraction
if (fstress_trans(j) > tol_math_check) then
plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+j, of) = constitutive_dislotwin_Cgro(instance)*&
(1.0_pReal - sumf - sumftr)*fstress_trans(j)
else
plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+j, of) = 0.0_pReal
endif
!* Probability rate of fault band intersection for strain-induced martensite nucleation
select case(lattice_structure(ph))
case (LATTICE_fcc_ID)
a = lattice_fcc_transNucleationTwinPair(1,j)
b = lattice_fcc_transNucleationTwinPair(2,j)
sa = sign(1, a)
sb = sign(1, b)
ssa = sign(1.0, shearrate_trans(a))
ssb = sign(1.0, shearrate_trans(b))
if (sa == ssa .and. sb == ssb) then
probrate_trans(j) = (abs(shear_trans(a)*shearrate_trans(b)) + abs(shear_trans(b)*shearrate_trans(a)))
else
probrate_trans(j) = 0.0_pReal
endif
probrate_trans(j) = probrate_trans(j)/lattice_fcc_shearCritTrans**2
!* Dotstate for strain-induced martensite volume fraction
plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+nr+j, of) = constitutive_dislotwin_Cnuc(instance)*&
(1.0_pReal - sumf - sumftr)*probrate_trans(j)
case default
plasticState(ph)%dotState(3_pInt*ns+2_pInt*nt+nr+j, of) = 0.0_pReal
end select
enddo
enddo
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
end subroutine constitutive_dislotwin_dotState
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
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_dislotwin_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_dislotwin_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_dislotwin_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_dislotwin_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 rate
!--------------------------------------------------------------------------------------------------
subroutine constitutive_dislotwin_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_dislotwin_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_dislotwin_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_dislotwin_getSlipRate
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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_dislotwin_postResults(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 prec, only: &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
tol_math_check
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
use math, only: &
pi, &
math_Mandel6to33, &
math_spectralDecompositionSym33
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
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_Sslip_v, &
lattice_Stwin_v, &
lattice_maxNslipFamily, &
lattice_maxNtwinFamily, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_maxNtransFamily, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
lattice_NslipSystem, &
lattice_NtwinSystem, &
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
lattice_NtransSystem, &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
lattice_shearTwin, &
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_mu, &
lattice_structure, &
lattice_fcc_twinNucleationSlipPair, &
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
LATTICE_fcc_ID
sorry, last version did not compile. now with correct encoding.
2011-04-13 17:21:46 +05:30
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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) :: &
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +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) :: &
removed temperature integration and corresponding data structures and debugging options temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver. introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver. went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 18:34:59 +05:30
temperature !< 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) :: &
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
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
real(pReal), dimension(constitutive_dislotwin_sizePostResults(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
constitutive_dislotwin_postResults
integer(pInt) :: &
removed unused variables
2014-09-23 02:04:42 +05:30
instance,&
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
ns,nt,nr,&
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
f,o,i,c,j,index_myFamily,&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
s1,s2, &
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,tau,StressRatio_p,StressRatio_pminus1,BoltzmannRatio,DotGamma0,StressRatio_r,Ndot0,dgdot_dtauslip
real(preal), dimension(constitutive_dislotwin_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_slip
real(pReal), dimension(3,3) :: eigVectors
real(pReal), dimension (3) :: eigValues
logical :: error
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_dislotwin_totalNslip(instance)
nt = constitutive_dislotwin_totalNtwin(instance)
Modification of basic states for phase transformation
2014-08-14 18:48:33 +05:30
nr = constitutive_dislotwin_totalNtrans(instance)
prepared new State for dislotwin
2014-06-11 17:41:14 +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 to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
sumf = sum(plasticState(ph)%state((3_pInt*ns+1_pInt):(3_pInt*ns+nt), of)) ! 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
c = 0_pInt
constitutive_dislotwin_postResults = 0.0_pReal
!* Spectral decomposition of stress
call math_spectralDecompositionSym33(math_Mandel6to33(Tstar_v),eigValues,eigVectors, error)
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_dislotwin_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_dislotwin_outputID(o,instance))
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (edge_density_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(1_pInt:ns, of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (dipole_density_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+ns) = plasticState(ph)%state(ns+1_pInt:2_pInt*ns, of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (shear_rate_slip_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
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
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* 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 = dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of)) > tol_math_check) then
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_p = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**constitutive_dislotwin_pPerSlipFamily(f,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_pminus1 = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**(constitutive_dislotwin_pPerSlipFamily(f,instance)-1.0_pReal)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Boltzmann ratio
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
BoltzmannRatio = constitutive_dislotwin_QedgePerSlipSystem(j,instance)/(kB*Temperature)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Initial shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
DotGamma0 = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)* &
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_v0PerSlipSystem(j,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Shear rates due to slip
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_postResults(c+j) = &
DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)**&
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_qPerSlipFamily(f,instance))*sign(1.0_pReal,tau)
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
else
constitutive_dislotwin_postResults(c+j) = 0.0_pReal
endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo ; enddo
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (accumulated_shear_slip_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+ns) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state((2_pInt*ns+1_pInt):(3_pInt*ns), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (mfp_slip_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+ns) =&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state((5_pInt*ns+3_pInt*nt+2*nr+1_pInt):(6_pInt*ns+3_pInt*nt+2*nr), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (resolved_stress_slip_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = j + 1_pInt
constitutive_dislotwin_postResults(c+j) =&
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo; enddo
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (threshold_stress_slip_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+ns) = &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state((6_pInt*ns+4_pInt*nt+2*nr+1_pInt):(7_pInt*ns+4_pInt*nt+2*nr), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (edge_dipole_distance_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = j + 1_pInt
constitutive_dislotwin_postResults(c+j) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(3.0_pReal*lattice_mu(ph)*constitutive_dislotwin_burgersPerSlipSystem(j,instance))/&
(16.0_pReal*pi*abs(dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))))
constitutive_dislotwin_postResults(c+j)=min(constitutive_dislotwin_postResults(c+j),&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(5*ns+3*nt+2*nr+j, of))
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
! constitutive_dislotwin_postResults(c+j)=max(constitutive_dislotwin_postResults(c+j),&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
! plasticState(ph)%state(4*ns+2*nt+2*nr+j, of))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo; enddo
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (resolved_stress_shearband_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do j = 1_pInt,6_pInt ! loop over all shearband families
constitutive_dislotwin_postResults(c+j) = dot_product(Tstar_v, &
constitutive_dislotwin_sbSv(1:6,j,ipc,ip,el))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo
c = c + 6_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (shear_rate_shearband_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do j = 1_pInt,6_pInt ! loop over all shearbands
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Resolved shear stress on shearband system
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
tau = dot_product(Tstar_v,constitutive_dislotwin_sbSv(1:6,j,ipc,ip,el))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Stress ratios
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
if (abs(tau) < tol_math_check) then
StressRatio_p = 0.0_pReal
StressRatio_pminus1 = 0.0_pReal
else
StressRatio_p = (abs(tau)/constitutive_dislotwin_sbResistance(instance))**&
constitutive_dislotwin_pShearBand(instance)
StressRatio_pminus1 = (abs(tau)/constitutive_dislotwin_sbResistance(instance))**&
(constitutive_dislotwin_pShearBand(instance)-1.0_pReal)
endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Boltzmann ratio
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
BoltzmannRatio = constitutive_dislotwin_sbQedge(instance)/(kB*Temperature)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Initial shear rates
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
DotGamma0 = constitutive_dislotwin_sbVelocity(instance)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
! Shear rate due to shear band
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
constitutive_dislotwin_postResults(c+j) = &
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)**constitutive_dislotwin_qShearBand(instance))*&
sign(1.0_pReal,tau)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo
c = c + 6_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (twin_fraction_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+nt) = plasticState(ph)%state((3_pInt*ns+1_pInt):(3_pInt*ns+nt), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (shear_rate_twin_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
if (nt > 0_pInt) then
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
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
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* 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 = dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of)) > tol_math_check) then
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_p = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**constitutive_dislotwin_pPerSlipFamily(f,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_pminus1 = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**(constitutive_dislotwin_pPerSlipFamily(f,instance)-1.0_pReal)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Boltzmann ratio
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
BoltzmannRatio = constitutive_dislotwin_QedgePerSlipSystem(j,instance)/(kB*Temperature)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Initial shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
DotGamma0 = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)* &
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_v0PerSlipSystem(j,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Shear rates due to slip
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
gdot_slip(j) = DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)**&
constitutive_dislotwin_qPerSlipFamily(f,instance))*sign(1.0_pReal,tau)
else
gdot_slip(j) = 0.0_pReal
endif
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo;enddo
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNtwinFamily ! loop over all twin families
index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1,constitutive_dislotwin_Ntwin(f,instance) ! process each (active) twin system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
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
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* 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 = dot_product(Tstar_v,lattice_Stwin_v(:,index_myFamily+i,ph))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_r = (plasticState(ph)%state(7_pInt*ns+4_pInt*nt+2*nr+j, of)/ &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
tau)**constitutive_dislotwin_rPerTwinFamily(f,instance)
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
!* Shear rates due to twin
if ( tau > 0.0_pReal ) then
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
select case(lattice_structure(ph))
change latticeName to latticeID
2013-11-27 21:50:27 +05:30
case (LATTICE_fcc_ID)
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
s1=lattice_fcc_twinNucleationSlipPair(1,index_myFamily+i)
s2=lattice_fcc_twinNucleationSlipPair(2,index_myFamily+i)
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
if (tau < constitutive_dislotwin_tau_r(j,instance)) then
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
Ndot0=(abs(gdot_slip(s1))*(plasticState(ph)%state(s2, of)+plasticState(ph)%state(ns+s2, of))+&
abs(gdot_slip(s2))*(plasticState(ph)%state(s1, of)+plasticState(ph)%state(ns+s1, of)))/&
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
(constitutive_dislotwin_L0(instance)*&
constitutive_dislotwin_burgersPerSlipSystem(j,instance))*&
(1.0_pReal-exp(-constitutive_dislotwin_VcrossSlip(instance)/(kB*Temperature)*&
(constitutive_dislotwin_tau_r(j,instance)-tau)))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
else
Ndot0=0.0_pReal
end if
case default
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 15:48:40 +05:30
Ndot0=constitutive_dislotwin_Ndot0PerTwinSystem(j,instance)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
end select
constitutive_dislotwin_postResults(c+j) = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
(constitutive_dislotwin_MaxTwinFraction(instance)-sumf)*lattice_shearTwin(index_myFamily+i,ph)*&
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
plasticState(ph)%state(7_pInt*ns+5_pInt*nt+2*nr+j, of)*Ndot0*exp(-StressRatio_r)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
endif
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo ; enddo
endif
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (accumulated_shear_twin_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+nt) = plasticState(ph)% &
state((3_pInt*ns+nt+1_pInt) :(3_pInt*ns+2_pInt*nt), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (mfp_twin_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+nt) = plasticState(ph)% &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
state((6_pInt*ns+3_pInt*nt+2_pInt*nr+1_pInt):(6_pInt*ns+4_pInt*nt+2_pInt*nr), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (resolved_stress_twin_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
if (nt > 0_pInt) then
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNtwinFamily ! loop over all slip families
index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Ntwin(f,instance) ! process each (active) slip system in family
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = j + 1_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+j) = dot_product(Tstar_v,lattice_Stwin_v(:,index_myFamily+i,ph))
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
enddo; enddo
endif
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (threshold_stress_twin_ID)
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+nt) = plasticState(ph)% &
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
state((7_pInt*ns+4_pInt*nt+2_pInt*nr+1_pInt):(7_pInt*ns+5_pInt*nt+2_pInt*nr), of)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + nt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (stress_exponent_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
j = 0_pInt
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
do f = 1_pInt,lattice_maxNslipFamily ! loop over all slip families
index_myFamily = sum(lattice_NslipSystem(1:f-1_pInt,ph)) ! at which index starts my family
do i = 1_pInt,constitutive_dislotwin_Nslip(f,instance) ! process each (active) slip system in 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
j = j + 1_pInt
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25: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
!* 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 = dot_product(Tstar_v,lattice_Sslip_v(:,1,index_myFamily+i,ph))
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
if((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of)) > tol_math_check) then
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
!* Stress ratios
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_p = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**constitutive_dislotwin_pPerSlipFamily(f,instance)
Added Ctrans to homogenizedC and strain-induced martensite part to basic states
2014-09-10 17:42:17 +05:30
StressRatio_pminus1 = ((abs(tau)-plasticState(ph)%state(6*ns+4*nt+2*nr+j, of))/&
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance)))&
**(constitutive_dislotwin_pPerSlipFamily(f,instance)-1.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
!* Boltzmann ratio
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
BoltzmannRatio = constitutive_dislotwin_QedgePerSlipSystem(j,instance)/(kB*Temperature)
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
!* Initial shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
DotGamma0 = &
changed to new state, please report bugs to Luv or Martin
2014-07-02 17:57:39 +05:30
plasticState(ph)%state(j, of)*constitutive_dislotwin_burgersPerSlipSystem(j,instance)* &
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
constitutive_dislotwin_v0PerSlipSystem(j,instance)
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25: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
!* Shear rates due to slip
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
gdot_slip(j) = DotGamma0*exp(-BoltzmannRatio*(1_pInt-StressRatio_p)**&
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25:40 +05:30
constitutive_dislotwin_qPerSlipFamily(f,instance))*sign(1.0_pReal,tau)
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
!* Derivatives of shear rates
modified constitutive description in line with other dislocation density based models
2014-03-30 20:34:06 +05:30
dgdot_dtauslip = &
abs(gdot_slip(j))*BoltzmannRatio*constitutive_dislotwin_pPerSlipFamily(f,instance)&
*constitutive_dislotwin_qPerSlipFamily(f,instance)/&
(constitutive_dislotwin_SolidSolutionStrength(instance)+&
constitutive_dislotwin_tau_peierlsPerSlipFamily(f,instance))*&
StressRatio_pminus1*(1-StressRatio_p)**(constitutive_dislotwin_qPerSlipFamily(f,instance)-1.0_pReal)
else
gdot_slip(j) = 0.0_pReal
dgdot_dtauslip = 0.0_pReal
endif
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code forgot to commit dislotwin last time, now seems to work
2014-03-12 05:25: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
!* Stress exponent
if (gdot_slip(j)==0.0_pReal) then
constitutive_dislotwin_postResults(c+j) = 0.0_pReal
else
constitutive_dislotwin_postResults(c+j) = (tau/gdot_slip(j))*dgdot_dtauslip
endif
enddo ; enddo
c = c + ns
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (sb_eigenvalues_ID)
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
forall (j = 1_pInt:3_pInt) &
constitutive_dislotwin_postResults(c+j) = eigValues(j)
c = c + 3_pInt
added enum for dislotwin output, fixed bug when using recursive file input function
2013-12-12 05:12:33 +05:30
case (sb_eigenvectors_ID)
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 13:34:05 +05:30
constitutive_dislotwin_postResults(c+1_pInt:c+9_pInt) = reshape(eigVectors,[9])
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
c = c + 9_pInt
Outputs for TRIP volume fraction
2014-09-24 12:56:27 +05:30
case (stress_trans_fraction_ID)
constitutive_dislotwin_postResults(c+1_pInt:c+nr) = &
plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of)
c = c + nr
case (strain_trans_fraction_ID)
constitutive_dislotwin_postResults(c+1_pInt:c+nr) = &
plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of)
c = c + nr
case (trans_fraction_ID)
constitutive_dislotwin_postResults(c+1_pInt:c+nr) = &
plasticState(ph)%state((3_pInt*ns+2_pInt*nt+1_pInt):(3_pInt*ns+2_pInt*nt+nr), of) + &
plasticState(ph)%state((3_pInt*ns+2_pInt*nt+nr+1_pInt):(3_pInt*ns+2_pInt*nt+2_pInt*nr), of)
c = c + nr
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 16:24:45 +05:30
end select
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 21:57:26 +05:30
enddo
end function constitutive_dislotwin_postResults
end module constitutive_dislotwin