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

put recalculation of elastic stress back in (was accidentally deleted) 

changed check of residuum improvement for small residuum
This commit is contained in:
Franz Roters 2007-08-13 07:12:24 +00:00
parent 85df3b95ee
commit 28fc9cc221
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
trunk/CPFEM.f90
View File

@ -397,7 +397,8 @@ stress: do ! inner iteration: stress
forall(i=1:3) T_elastic(i)=T_elastic(i)-p_hydro forall(i=1:3) T_elastic(i)=T_elastic(i)-p_hydro
Rstress = Tstar_v - T_elastic Rstress = Tstar_v - T_elastic
! step size control: if residuum does not improve redo iteration with reduced step size ! step size control: if residuum does not improve redo iteration with reduced step size
if(maxval(abs(Rstress)) > maxval(abs(Rstress_old)) .and. iStress > 1) then if(maxval(abs(Rstress)) > maxval(abs(Rstress_old)) .and. &
maxval(abs(Rstress)) > 1.0e-6 .and. iStress > 1) then
! write(6,*) 'Hallo', iStress ! write(6,*) 'Hallo', iStress
Tstar_v=Tstar_v+0.5*dTstar_v Tstar_v=Tstar_v+0.5*dTstar_v
dTstar_v=0.5*dTstar_v dTstar_v=0.5*dTstar_v
@ -437,6 +438,7 @@ stress: do ! inner iteration: stress
enddo stress enddo stress
! write(6,*) 'istress', istress ! write(6,*) 'istress', istress
Tstar_v = 0.5_pReal*matmul(C_66,math_Mandel33to6(matmul(transpose(B),AB)-math_I3))
dstate = dt*constitutive_dotState(Tstar_v,state_new,grain,CPFEM_in,cp_en) ! evolution of microstructure dstate = dt*constitutive_dotState(Tstar_v,state_new,grain,CPFEM_in,cp_en) ! evolution of microstructure
Rstate = state_new - (state_old+dstate) Rstate = state_new - (state_old+dstate)
RstateS = 0.0_pReal RstateS = 0.0_pReal
@ -458,7 +460,6 @@ stress: do ! inner iteration: stress
endif endif
Fp_new = Fp_new*det**(1.0_pReal/3.0_pReal) ! det = det(InvFp_new) !! Fp_new = Fp_new*det**(1.0_pReal/3.0_pReal) ! det = det(InvFp_new) !!
Fe_new = matmul(Fg_new,invFp_new) Fe_new = matmul(Fg_new,invFp_new)
return return
END SUBROUTINE END SUBROUTINE