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

generic interfaces makes sense

This commit is contained in:
Sharan Roongta 2020-07-12 13:27:28 +02:00
parent 3a5e3b36c1
commit 5602abe690
6 changed files with 25 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
src/constitutive.f90
View File

@ -140,7 +140,7 @@ module constitutive
el !< element
end function plastic_dislotwin_homogenizedC
module subroutine constitutive_dependentState(F, Fp, ipc, ip, el)
module subroutine constitutive_plastic_dependentState(F, Fp, ipc, ip, el)
integer, intent(in) :: &
ipc, & !< component-ID of integration point
@ -149,9 +149,9 @@ module constitutive
real(pReal), intent(in), dimension(3,3) :: &
F, & !< elastic deformation gradient
Fp !< plastic deformation gradient
end subroutine constitutive_dependentState
end subroutine constitutive_plastic_dependentState
module subroutine constitutive_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, &
module subroutine constitutive_plastic_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, &
S, Fi, ipc, ip, el)
integer, intent(in) :: &
ipc, & !< component-ID of integration point
@ -166,7 +166,7 @@ module constitutive
dLp_dS, &
dLp_dFi !< derivative of Lp with respect to Fi
end subroutine constitutive_LpAndItsTangents
end subroutine constitutive_plastic_LpAndItsTangents
pure module function kinematics_thermal_expansion_initialStrain(homog,phase,offset) result(initialStrain)
@ -351,6 +351,18 @@ module constitutive
end interface
interface constitutive_LpAndItsTangents
module procedure :: constitutive_plastic_LpAndItsTangents
end interface constitutive_LpAndItsTangents
interface constitutive_dependentState
module procedure :: constitutive_plastic_dependentState
end interface constitutive_dependentState
interface constitutive_getRateAndItsTangents
module procedure :: constitutive_damage_getRateAndItsTangents , &
constitutive_thermal_getRateAndItsTangents
end interface constitutive_getRateAndItsTangents
type :: tDebugOptions
logical :: &
@ -376,8 +388,7 @@ module constitutive
constitutive_collectDotState, &
constitutive_deltaState, &
plastic_nonlocal_updateCompatibility, &
constitutive_damage_getRateAndItsTangents, &
constitutive_thermal_getRateAndItsTangents, &
constitutive_getRateAndItsTangents, &
constitutive_results
contains

8
src/constitutive_plastic.f90
View File

@ -171,7 +171,7 @@ end subroutine plastic_init
!--------------------------------------------------------------------------------------------------
!> @brief calls microstructure function of the different constitutive models
!--------------------------------------------------------------------------------------------------
module procedure constitutive_dependentState
module procedure constitutive_plastic_dependentState
integer :: &
ho, & !< homogenization
@ -192,14 +192,14 @@ module procedure constitutive_dependentState
call plastic_nonlocal_dependentState (F,Fp,instance,of,ip,el)
end select plasticityType
end procedure constitutive_dependentState
end procedure constitutive_plastic_dependentState
!--------------------------------------------------------------------------------------------------
!> @brief contains the constitutive equation for calculating the velocity gradient
! ToDo: Discuss whether it makes sense if crystallite handles the configuration conversion, i.e.
! Mp in, dLp_dMp out
!--------------------------------------------------------------------------------------------------
module procedure constitutive_LpAndItsTangents
module procedure constitutive_plastic_LpAndItsTangents
real(pReal), dimension(3,3,3,3) :: &
dLp_dMp !< derivative of Lp with respect to Mandel stress
@ -250,7 +250,7 @@ module procedure constitutive_LpAndItsTangents
dLp_dS(i,j,1:3,1:3) = matmul(matmul(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3)) ! ToDo: @PS: why not: dLp_dMp:(FiT Fi)
enddo; enddo
end procedure constitutive_LpAndItsTangents
end procedure constitutive_plastic_LpAndItsTangents
end submodule constitutive_plastic

2
src/damage_local.f90
View File

@ -131,7 +131,7 @@ subroutine damage_local_getSourceAndItsTangent(phiDot, dPhiDot_dPhi, phi, ip, el
phiDot = 0.0_pReal
dPhiDot_dPhi = 0.0_pReal
call constitutive_damage_getRateAndItsTangents(phiDot, dPhiDot_dPhi, phi, ip, el)
call constitutive_getRateAndItsTangents(phiDot, dPhiDot_dPhi, phi, ip, el)
phiDot = phiDot/real(homogenization_Ngrains(material_homogenizationAt(el)),pReal)
dPhiDot_dPhi = dPhiDot_dPhi/real(homogenization_Ngrains(material_homogenizationAt(el)),pReal)

2
src/damage_nonlocal.f90
View File

@ -100,7 +100,7 @@ subroutine damage_nonlocal_getSourceAndItsTangent(phiDot, dPhiDot_dPhi, phi, ip,
phiDot = 0.0_pReal
dPhiDot_dPhi = 0.0_pReal
call constitutive_damage_getRateAndItsTangents(phiDot, dPhiDot_dPhi, phi, ip, el)
call constitutive_getRateAndItsTangents(phiDot, dPhiDot_dPhi, phi, ip, el)
phiDot = phiDot/real(homogenization_Ngrains(material_homogenizationAt(el)),pReal)
dPhiDot_dPhi = dPhiDot_dPhi/real(homogenization_Ngrains(material_homogenizationAt(el)),pReal)

2
src/thermal_adiabatic.f90
View File

@ -131,7 +131,7 @@ subroutine thermal_adiabatic_getSourceAndItsTangent(Tdot, dTdot_dT, T, ip, el)
dTdot_dT = 0.0_pReal
homog = material_homogenizationAt(el)
call constitutive_thermal_getRateAndItsTangents(TDot, dTDot_dT, T, crystallite_S, crystallite_Lp, ip, el)
call constitutive_getRateAndItsTangents(TDot, dTDot_dT, T, crystallite_S, crystallite_Lp, ip, el)
Tdot = Tdot/real(homogenization_Ngrains(homog),pReal)
dTdot_dT = dTdot_dT/real(homogenization_Ngrains(homog),pReal)

2
src/thermal_conduction.f90
View File

@ -90,7 +90,7 @@ subroutine thermal_conduction_getSourceAndItsTangent(Tdot, dTdot_dT, T, ip, el)
dTdot_dT = 0.0_pReal
homog = material_homogenizationAt(el)
call constitutive_thermal_getRateAndItsTangents(TDot, dTDot_dT, T, crystallite_S,crystallite_Lp ,ip, el)
call constitutive_getRateAndItsTangents(TDot, dTDot_dT, T, crystallite_S,crystallite_Lp ,ip, el)
Tdot = Tdot/real(homogenization_Ngrains(homog),pReal)
dTdot_dT = dTdot_dT/real(homogenization_Ngrains(homog),pReal)