DAMASK_EICMD/src/phase_mechanical_eigen_ther...

116 lines
4.7 KiB
Fortran
Raw Normal View History

!--------------------------------------------------------------------------------------------------
!> @author Pratheek Shanthraj, Max-Planck-Institut für Eisenforschung GmbH
!> @brief material subroutine incorporating kinematics resulting from thermal expansion
!> @details to be done
!--------------------------------------------------------------------------------------------------
2021-02-13 23:27:41 +05:30
submodule(phase:eigen) thermalexpansion
2020-02-28 14:55:07 +05:30
2020-02-29 14:06:42 +05:30
integer, dimension(:), allocatable :: kinematics_thermal_expansion_instance
2019-05-17 02:44:47 +05:30
type :: tParameters
2020-02-29 15:40:23 +05:30
real(pReal) :: &
T_ref
2020-02-29 16:51:03 +05:30
real(pReal), dimension(3,3,3) :: &
2020-09-23 04:46:12 +05:30
A = 0.0_pReal
end type tParameters
2020-02-28 14:55:07 +05:30
type(tParameters), dimension(:), allocatable :: param
2020-02-28 14:55:07 +05:30
contains
!--------------------------------------------------------------------------------------------------
!> @brief module initialization
!> @details reads in material parameters, allocates arrays, and does sanity checks
!--------------------------------------------------------------------------------------------------
2021-02-13 23:11:30 +05:30
module function thermalexpansion_init(kinematics_length) result(myKinematics)
2020-02-28 14:55:07 +05:30
integer, intent(in) :: kinematics_length
2020-08-15 19:32:10 +05:30
logical, dimension(:,:), allocatable :: myKinematics
2020-02-28 14:55:07 +05:30
integer :: Ninstances,p,i,k
2020-08-15 19:32:10 +05:30
class(tNode), pointer :: &
phases, &
phase, &
mech, &
2020-08-15 19:32:10 +05:30
kinematics, &
kinematic_type
2021-02-13 23:27:41 +05:30
print'(/,a)', ' <<<+- phase:mechanical:eigen:thermalexpansion init -+>>>'
2020-08-15 19:32:10 +05:30
myKinematics = kinematics_active('thermalexpansion',kinematics_length)
Ninstances = count(myKinematics)
print'(a,i2)', ' # phases: ',Ninstances; flush(IO_STDOUT)
if(Ninstances == 0) return
2020-02-28 14:55:07 +05:30
phases => config_material%get('phase')
allocate(param(Ninstances))
2020-08-15 19:32:10 +05:30
allocate(kinematics_thermal_expansion_instance(phases%length), source=0)
do p = 1, phases%length
if(any(myKinematics(:,p))) kinematics_thermal_expansion_instance(p) = count(myKinematics(:,1:p))
phase => phases%get(p)
2020-08-15 19:32:10 +05:30
if(count(myKinematics(:,p)) == 0) cycle
mech => phase%get('mechanical')
kinematics => mech%get('eigen')
2020-08-15 19:32:10 +05:30
do k = 1, kinematics%length
if(myKinematics(k,p)) then
associate(prm => param(kinematics_thermal_expansion_instance(p)))
kinematic_type => kinematics%get(k)
prm%T_ref = kinematic_type%get_asFloat('T_ref', defaultVal=0.0_pReal)
prm%A(1,1,1) = kinematic_type%get_asFloat('A_11')
prm%A(1,1,2) = kinematic_type%get_asFloat('A_11,T',defaultVal=0.0_pReal)
prm%A(1,1,3) = kinematic_type%get_asFloat('A_11,T^2',defaultVal=0.0_pReal)
if (any(phase_lattice(p) == ['hP','tI'])) then
prm%A(3,3,1) = kinematic_type%get_asFloat('A_33')
prm%A(3,3,2) = kinematic_type%get_asFloat('A_33,T',defaultVal=0.0_pReal)
prm%A(3,3,3) = kinematic_type%get_asFloat('A_33,T^2',defaultVal=0.0_pReal)
endif
do i=1, size(prm%A,3)
prm%A(1:3,1:3,i) = lattice_applyLatticeSymmetry33(prm%A(1:3,1:3,i),phase_lattice(p))
enddo
2020-08-15 19:32:10 +05:30
end associate
endif
2020-02-29 17:27:19 +05:30
enddo
enddo
2021-02-13 23:11:30 +05:30
end function thermalexpansion_init
!--------------------------------------------------------------------------------------------------
2020-06-26 15:14:17 +05:30
!> @brief constitutive equation for calculating the velocity gradient
!--------------------------------------------------------------------------------------------------
module subroutine thermalexpansion_LiAndItsTangent(Li, dLi_dTstar, ph,me)
2020-02-28 14:55:07 +05:30
2021-01-24 15:49:10 +05:30
integer, intent(in) :: ph, me
real(pReal), intent(out), dimension(3,3) :: &
2019-06-16 00:02:53 +05:30
Li !< thermal velocity gradient
real(pReal), intent(out), dimension(3,3,3,3) :: &
2019-06-16 00:02:53 +05:30
dLi_dTstar !< derivative of Li with respect to Tstar (4th-order tensor defined to be zero)
2020-02-29 14:06:42 +05:30
2021-01-24 15:49:10 +05:30
real(pReal) :: T, dot_T
2020-02-28 14:55:07 +05:30
2021-05-24 00:32:35 +05:30
T = thermal_T(ph,me)
dot_T = thermal_dot_T(ph,me)
2020-02-28 14:55:07 +05:30
2021-01-24 15:49:10 +05:30
associate(prm => param(kinematics_thermal_expansion_instance(ph)))
Li = dot_T * ( &
prm%A(1:3,1:3,1)*(T - prm%T_ref)**0 & ! constant coefficient
+ prm%A(1:3,1:3,2)*(T - prm%T_ref)**1 & ! linear coefficient
+ prm%A(1:3,1:3,3)*(T - prm%T_ref)**2 & ! quadratic coefficient
) / &
(1.0_pReal &
+ prm%A(1:3,1:3,1)*(T - prm%T_ref)**1 / 1. &
+ prm%A(1:3,1:3,2)*(T - prm%T_ref)**2 / 2. &
+ prm%A(1:3,1:3,3)*(T - prm%T_ref)**3 / 3. &
)
end associate
2020-02-28 14:55:07 +05:30
dLi_dTstar = 0.0_pReal
end subroutine thermalexpansion_LiAndItsTangent
2021-01-26 05:50:45 +05:30
end submodule thermalexpansion