32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "mvsiz_p.inc"
40
41
42
43 INTEGER JFT, JLT, NFT, IEPSDOT, ISTRAIN,NEL
45 . vdef(mvsiz,8), gstr(nel,8),
46 . exx(mvsiz),eyy(mvsiz),exy(mvsiz),exz(mvsiz),eyz(mvsiz),
47 . kxx(mvsiz),kyy(mvsiz),kxy(mvsiz),dt1c(mvsiz),
48 . epsdot(6,*)
49
50
51
52 INTEGER I, J
53
54 IF(iepsdot/=0)THEN
55 DO i=jft,jlt
56 j = i + nft
57 epsdot(1,j) = vdef(i,1)
58 epsdot(2,j) = vdef(i,2)
59 epsdot(3,j) = vdef(i,3)
60 epsdot(4,j) = vdef(i,6)
61 epsdot(5,j) = vdef(i,7)
62 epsdot(6,j) = vdef(i,8)
63 ENDDO
64 ENDIF
65
66 DO i=jft,jlt
67 exx(i) = vdef(i,1) * dt1c(i)
68 eyy(i) = vdef(i,2) * dt1c(i)
69 exy(i) = vdef(i,3) * dt1c(i)
70 eyz(i) = 0.0
71 exz(i) = 0.0
72 kxx(i) = vdef(i,6) * dt1c(i)
73 kyy(i) = vdef(i,7) * dt1c(i)
74 kxy(i) = vdef(i,8) * dt1c(i)
75 ENDDO
76
77 IF(istrain/=0)THEN
78 DO i=jft,jlt
79 gstr(i,1)=gstr(i,1)+exx(i)
80 gstr(i,2)=gstr(i,2)+eyy(i)
81 gstr(i,3)=gstr(i,3)+exy(i)
82
83
84 gstr(i,6)=gstr(i,6)+kxx(i)
85 gstr(i,7)=gstr(i,7)+kyy(i)
86 gstr(i,8)=gstr(i,8)+kxy(i)
87 ENDDO
88 ENDIF
89
90 RETURN