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

CPFEM_general call now uses CPFEM_en instead of cp_en -- the required mapping is not in place on first call!

This commit is contained in:
Philip Eisenlohr 2007-04-23 13:25:09 +00:00
parent 471d53e3ba
commit 83d2dea8c6
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
trunk/mpie_cpfem_marc.f90
View File

@ -169,13 +169,13 @@
! mpie_en element number ! mpie_en element number
! mpie_in intergration point number ! mpie_in intergration point number
!******************************************************************** !********************************************************************
cp_en = mesh_FEasCP('elem', n(1))
if ((lovl==6).or.(inc==0)) then if ((lovl==6).or.(inc==0)) then
call CPFEM_general(ffn, ffn1, inc, incsub, ncycle, timinc, cp_en, nn) call CPFEM_general(ffn, ffn1, inc, incsub, ncycle, timinc, n(1), nn)
endif endif
! return stress and jacobi ! return stress and jacobi
! Mandel: 11, 22, 33, 12, 23, 13 ! Mandel: 11, 22, 33, 12, 23, 13
! Marc: 11, 22, 33, 12, 23, 13 ! Marc: 11, 22, 33, 12, 23, 13
cp_en = mesh_FEasCP('elem', n(1))
s(1:ngens)=invnrmMandel(1:ngens)*CPFEM_stress_all(1:ngens, nn, cp_en) s(1:ngens)=invnrmMandel(1:ngens)*CPFEM_stress_all(1:ngens, nn, cp_en)
d(1:ngens,1:ngens)=CPFEM_jaco_old(1:ngens,1:ngens, nn, cp_en) d(1:ngens,1:ngens)=CPFEM_jaco_old(1:ngens,1:ngens, nn, cp_en)
forall(i=1:ngens) d(i,1:ngens)=d(i,1:ngens)*invnrmMandel(1:ngens) forall(i=1:ngens) d(i,1:ngens)=d(i,1:ngens)*invnrmMandel(1:ngens)