30
31
32
33#include "implicit_f.inc"
34
35
36
37#include "scr14_c.inc"
38
39
40
41 INTEGER :: JFT, JLT,NFT,ISTRAIN,NEL
42 my_real :: gstr(nel,8), exx(*), eyy(*), exy(*),kxx(*), kyy(*), kxy(*),epsdot(6,*),dt1,gstrpg(nel,8)
43
44
45
46 INTEGER I,J
47
48 IF(iepsdot /= 0)THEN
49 DO i=jft,jlt
50 j = i + nft
51 epsdot(4,j)=kxx(i)
52 epsdot(5,j)=kyy(i)
53 epsdot(6,j)=kxy(i)
54 ENDDO
55 ENDIF
56
57 DO i=jft,jlt
58 kxx(i)=kxx(i)*dt1
59 kyy(i)=kyy(i)*dt1
60 kxy(i)=kxy(i)*dt1
61 ENDDO
62
63 IF(istrain /= 0)THEN
64 DO i=jft,jlt
65 gstr(i,6)=gstr(i,6)+kxx(i)*third
66 gstr(i,7)=gstr(i,7)+kyy(i)*third
67 gstr(i,8)=gstr(i,8)+kxy(i)*third
68 gstrpg(i,1)=gstrpg(i,1)+exx(i)
69 gstrpg(i,2)=gstrpg(i,2)+eyy(i)
70 gstrpg(i,3)=gstrpg(i,3)+exy(i)
71 gstrpg(i,6)=gstrpg(i,6)+kxx(i)
72 gstrpg(i,7)=gstrpg(i,7)+kyy(i)
73 gstrpg(i,8)=gstrpg(i,8)+kxy(i)
74 ENDDO
75 ENDIF
76
77 RETURN