39
40
41
42#include "implicit_f.inc"
43
44
45
46#include "mvsiz_p.inc"
47
48
49
50
51
52
53 INTEGER NEL
54
56 . pxc1(*), pxc2(*), pxc3(*), pxc4(*),
57 . pyc1(*), pyc2(*), pyc3(*), pyc4(*),
58 . k11(9,*),k12(9,*),k13(9,*),k14(9,*) ,k15(9,*) ,
59 . k16(9,*),k22(9,*) ,k23(9,*) ,k24(9,*),k25(9,*),
60 . k26(9,*),k33(9,*),k34(9,*),k35(9,*),k36(9,*) ,
61 . k44(9,*),k45(9,*),k46(9,*) ,k55(9,*),k56(9,*),
62 . k66(9,*),vol(*),sig(nel,6)
63
64
65
66 INTEGER I,J
67
69 . fxx(mvsiz),fyy(mvsiz),fxy(mvsiz),
70 . px1(mvsiz),px2(mvsiz),px3(mvsiz),
71 . px4(mvsiz),px5(mvsiz),px6(mvsiz),
72 . py1(mvsiz),py2(mvsiz),py3(mvsiz),
73 . py4(mvsiz),py5(mvsiz),py6(mvsiz)
74
75 DO i=1,nel
76 fxx(i)=sig(i,1)*vol(i)
77 fyy(i)=sig(i,2)*vol(i)
78 fxy(i)=sig(i,4)*vol(i)
79 ENDDO
80
81 DO i=1,nel
82 px1(i) = pxc1(i)-pxc4(i)
83 px2(i) = pxc2(i)-pxc4(i)
84 px3(i) = pxc3(i)-pxc4(i)
85 px4(i) = pxc1(i)+pxc4(i)
86 px5(i) = pxc2(i)+pxc4(i)
87 px6(i) = pxc3(i)+pxc4(i)
88 ENDDO
89 DO i=1,nel
90 py1(i) = pyc1(i)-pyc4(i)
91 py2(i) = pyc2(i)-pyc4(i)
92 py3(i) = pyc3(i)-pyc4(i)
93 py4(i) = pyc1(i)+pyc4(i)
94 py5(i) = pyc2(i)+pyc4(i)
95 py6(i) = pyc3(i)+pyc4(i)
96 ENDDO
98 . px1, py1 ,px1, py1, fxx, fyy ,fxy,k11)
100 . px1, py1 ,px2, py2, fxx, fyy ,fxy,k12)
102 . px1, py1 ,px3, py3, fxx, fyy ,fxy,k13)
104 . px1, py1 ,px4, py4, fxx, fyy ,fxy,k14)
106 . px1, py1 ,px5, py5, fxx, fyy ,fxy,k15)
108 . px1, py1 ,px6, py6, fxx, fyy ,fxy,k16)
110 . px2, py2 ,px2, py2, fxx, fyy ,fxy,k22)
112 . px2, py2 ,px3, py3, fxx, fyy ,fxy,k23)
114 . px2, py2 ,px4, py4, fxx, fyy ,fxy,k24)
116 . px2, py2 ,px5, py5, fxx, fyy ,fxy,k25)
118 . px2, py2 ,px6, py6, fxx, fyy ,fxy,k26)
120 . px3, py3 ,px3, py3, fxx, fyy ,fxy,k33)
122 . px3, py3 ,px4, py4, fxx, fyy ,fxy,k34)
124 . px3, py3 ,px5, py5, fxx, fyy ,fxy,k35)
126 . px3, py3 ,px6, py6, fxx, fyy ,fxy,k36)
128 . px4, py4 ,px4, py4, fxx, fyy ,fxy,k44)
130 . px4, py4 ,px5, py5, fxx, fyy ,fxy,k45)
132 . px4, py4 ,px6, py6, fxx, fyy ,fxy,k46)
134 . px5, py5 ,px5, py5, fxx, fyy ,fxy,k55)
136 . px5, py5 ,px6, py6, fxx, fyy ,fxy,k56)
138 . px6, py6 ,px6, py6, fxx, fyy ,fxy,k66)
139
140 RETURN
subroutine s6ckeh3(lft, llt, pxi, pyi, pxj, pyj, fxx, fyy, fxy, kij)