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

removed some leftovers of the shearbanding part

This commit is contained in:
Martin Diehl 2014-11-05 09:53:08 +00:00
parent c2799585d9
commit a21f4d80db
1 changed files with 18 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
code/constitutive_dislokmc.f90
View File

@ -1376,10 +1376,13 @@ subroutine constitutive_dislokmc_dotState(Tstar_v,Temperature,ipc,ip,el)
use lattice, only: & use lattice, only: &
lattice_Sslip_v, & lattice_Sslip_v, &
lattice_Stwin_v, & lattice_Stwin_v, &
lattice_Sslip, &
lattice_Stwin, &
lattice_maxNslipFamily, & lattice_maxNslipFamily, &
lattice_maxNtwinFamily, & lattice_maxNtwinFamily, &
lattice_NslipSystem, & lattice_NslipSystem, &
lattice_NtwinSystem, & lattice_NtwinSystem, &
lattice_NnonSchmid, &
lattice_sheartwin, & lattice_sheartwin, &
lattice_mu, & lattice_mu, &
lattice_structure, & lattice_structure, &
@ -1397,7 +1400,7 @@ subroutine constitutive_dislokmc_dotState(Tstar_v,Temperature,ipc,ip,el)
ip, & !< integration point ip, & !< integration point
el !< element el !< element
integer(pInt) :: instance,ns,nt,f,i,j,index_myFamily,s1,s2, & integer(pInt) :: instance,ns,nt,f,i,j,k,index_myFamily,s1,s2, &
ph, & ph, &
of of
real(pReal) :: & real(pReal) :: &
@ -1421,6 +1424,10 @@ subroutine constitutive_dislokmc_dotState(Tstar_v,Temperature,ipc,ip,el)
ClimbVelocity, & ClimbVelocity, &
DotRhoEdgeDipClimb, & DotRhoEdgeDipClimb, &
DotRhoDipFormation DotRhoDipFormation
real(pReal), dimension(3,3,2) :: &
nonSchmid_tensor
real(pReal), dimension(3,3,3,3) :: &
dLp_dTstar3333
real(pReal), dimension(constitutive_dislokmc_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: & real(pReal), dimension(constitutive_dislokmc_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_slip_pos, tau_slip_neg gdot_slip_pos, tau_slip_neg
real(pReal), dimension(constitutive_dislokmc_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: & real(pReal), dimension(constitutive_dislokmc_totalNtwin(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
@ -1555,9 +1562,9 @@ subroutine constitutive_dislokmc_dotState(Tstar_v,Temperature,ipc,ip,el)
!* Twin volume fraction evolution !* Twin volume fraction evolution
j = 0_pInt j = 0_pInt
slipFamilies: do f = 1_pInt,lattice_maxNtwinFamily twinFamilies: do f = 1_pInt,lattice_maxNtwinFamily
index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family index_myFamily = sum(lattice_NtwinSystem(1:f-1_pInt,ph)) ! at which index starts my family
slipSystems: do i = 1_pInt,constitutive_dislokmc_Ntwin(f,instance) twinSystems: do i = 1_pInt,constitutive_dislokmc_Ntwin(f,instance)
j = j+1_pInt j = j+1_pInt
!* Resolved shear stress on twin system !* Resolved shear stress on twin system
@ -1590,8 +1597,8 @@ subroutine constitutive_dislokmc_dotState(Tstar_v,Temperature,ipc,ip,el)
plasticState(ph)%dotState(3_pInt*ns+nt+j, of) = plasticState(ph)%dotState(3_pInt*ns+j, of) * & plasticState(ph)%dotState(3_pInt*ns+nt+j, of) = plasticState(ph)%dotState(3_pInt*ns+j, of) * &
lattice_sheartwin(index_myfamily+i,ph) lattice_sheartwin(index_myfamily+i,ph)
endif endif
enddo slipSystems enddo twinSystems
enddo slipFamilies enddo twinFamilies
end subroutine constitutive_dislokmc_dotState end subroutine constitutive_dislokmc_dotState
@ -1692,9 +1699,7 @@ function constitutive_dislokmc_postResults(Tstar_v,Temperature,ipc,ip,el)
use prec, only: & use prec, only: &
tol_math_check tol_math_check
use math, only: & use math, only: &
pi, & pi
math_Mandel6to33, &
math_spectralDecompositionSym33
use mesh, only: & use mesh, only: &
mesh_NcpElems, & mesh_NcpElems, &
mesh_maxNips mesh_maxNips
@ -1708,10 +1713,13 @@ function constitutive_dislokmc_postResults(Tstar_v,Temperature,ipc,ip,el)
use lattice, only: & use lattice, only: &
lattice_Sslip_v, & lattice_Sslip_v, &
lattice_Stwin_v, & lattice_Stwin_v, &
lattice_Sslip, &
lattice_Stwin, &
lattice_maxNslipFamily, & lattice_maxNslipFamily, &
lattice_maxNtwinFamily, & lattice_maxNtwinFamily, &
lattice_NslipSystem, & lattice_NslipSystem, &
lattice_NtwinSystem, & lattice_NtwinSystem, &
lattice_NnonSchmid, &
lattice_shearTwin, & lattice_shearTwin, &
lattice_mu, & lattice_mu, &
lattice_structure, & lattice_structure, &
@ -1734,7 +1742,7 @@ function constitutive_dislokmc_postResults(Tstar_v,Temperature,ipc,ip,el)
integer(pInt) :: & integer(pInt) :: &
instance,& instance,&
ns,nt,& ns,nt,&
f,o,i,c,j,index_myFamily,& f,o,i,c,j,k,index_myFamily,&
s1,s2, & s1,s2, &
ph, & ph, &
of of
@ -1743,8 +1751,6 @@ function constitutive_dislokmc_postResults(Tstar_v,Temperature,ipc,ip,el)
real(pReal) :: StressRatio_u,StressRatio_uminus1 real(pReal) :: StressRatio_u,StressRatio_uminus1
real(preal), dimension(constitutive_dislokmc_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: & real(preal), dimension(constitutive_dislokmc_totalNslip(phase_plasticityInstance(material_phase(ipc,ip,el)))) :: &
gdot_slip_pos, vel_slip gdot_slip_pos, vel_slip
real(pReal), dimension(3,3) :: eigVectors
real(pReal), dimension (3) :: eigValues
logical :: error logical :: error
!* Shortened notation !* Shortened notation
@ -1761,8 +1767,6 @@ function constitutive_dislokmc_postResults(Tstar_v,Temperature,ipc,ip,el)
c = 0_pInt c = 0_pInt
constitutive_dislokmc_postResults = 0.0_pReal constitutive_dislokmc_postResults = 0.0_pReal
!* Spectral decomposition of stress
call math_spectralDecompositionSym33(math_Mandel6to33(Tstar_v),eigValues,eigVectors, error)
do o = 1_pInt,constitutive_dislokmc_Noutput(instance) do o = 1_pInt,constitutive_dislokmc_Noutput(instance)
select case(constitutive_dislokmc_outputID(o,instance)) select case(constitutive_dislokmc_outputID(o,instance))