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

corrected CPFEM_stressIP to use FFN1 at the end of the time increment 

calculate for ncycle==0 and not only inc==0 in hypela2
This commit is contained in:
Franz Roters 2007-05-14 12:13:36 +00:00
parent 6aece52d1c
commit 49886a072c
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
trunk/CPFEM.f90
View File

@ -193,7 +193,7 @@
Fg(:,:,i_then) = Fg(:,:,i_then)+deltaFg ! corresponding Fg Fg(:,:,i_then) = Fg(:,:,i_then)+deltaFg ! corresponding Fg
else ! full step solution else ! full step solution
t = CPFEM_dt ! final time t = CPFEM_dt ! final time
Fg(:,:,i_then) = CPFEM_ffn_all(:,:,CPFEM_in,cp_en) ! final Fg Fg(:,:,i_then) = CPFEM_ffn1_all(:,:,CPFEM_in,cp_en) ! final Fg
endif endif
call CPFEM_stressCrystallite(msg,cs,cd,Tstar_v,Fp(:,:,i_then),Fe,state(:,i_then),& call CPFEM_stressCrystallite(msg,cs,cd,Tstar_v,Fp(:,:,i_then),Fe,state(:,i_then),&

4
trunk/mpie_cpfem_marc.f90
View File

@ -156,7 +156,7 @@
dimension e(*),de(*),t(*),dt(*),g(*),d(ngens,*),s(*), n(2),coord(ncrd,*),disp(ndeg,*),matus(2),dispt(ndeg,*),ffn(itel,*),& dimension e(*),de(*),t(*),dt(*),g(*),d(ngens,*),s(*), n(2),coord(ncrd,*),disp(ndeg,*),matus(2),dispt(ndeg,*),ffn(itel,*),&
frotn(itel,*),strechn(itel),eigvn(itel,*),ffn1(itel,*),frotn1(itel,*),strechn1(itel),eigvn1(itel,*),kcus(2) frotn(itel,*),strechn(itel),eigvn(itel,*),ffn1(itel,*),frotn1(itel,*),strechn1(itel),eigvn1(itel,*),kcus(2)
! !
! call general material routine only in increment 0 and for lovl==6 (stress recovery) ! call general material routine only in cycle 0 and for lovl==6 (stress recovery)
! subroutine cpfem_general(mpie_ffn, mpie_ffn1, mpie_cn, mpie_tinc, mpie_enp, mpie_in) ! subroutine cpfem_general(mpie_ffn, mpie_ffn1, mpie_cn, mpie_tinc, mpie_enp, mpie_in)
!******************************************************************** !********************************************************************
@ -169,7 +169,7 @@
! mpie_en element number ! mpie_en element number
! mpie_in intergration point number ! mpie_in intergration point number
!******************************************************************** !********************************************************************
if ((lovl==6).or.(inc==0)) then if ((lovl==6).or.(ncycle==0)) then
call CPFEM_general(ffn, ffn1, inc, incsub, ncycle, timinc, n(1), nn) call CPFEM_general(ffn, ffn1, inc, incsub, ncycle, timinc, n(1), nn)
endif endif
! return stress and jacobi ! return stress and jacobi