48
49
50
51#include "implicit_f.inc"
52
53
54
55#include "mvsiz_p.inc"
56
57
58
59
60
61
62 INTEGER NEL
63
65 . px1(*), px2(*), px3(*), px4(*),
66 . py1(*), py2(*), py3(*), py4(*),
67 . pz1(*), pz2(*), pz3(*), pz4(*),
68 . px5(*), px6(*), px7(*), px8(*),
69 . py5(*), py6(*), py7(*), py8(*),
70 . pz5(*), pz6(*), pz7(*), pz8(*),
71 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*) ,k15(3,3,*) ,
72 . k16(3,3,*),k17(3,3,*),k18(3,3,*),k22(3,3,*) ,k23(3,3,*) ,
73 . k24(3,3,*),k25(3,3,*),k26(3,3,*),k27(3,3,*) ,k28(3,3,*) ,
74 . k33(3,3,*),k34(3,3,*),k35(3,3,*),k36(3,3,*) ,k37(3,3,*) ,
75 . k38(3,3,*),k44(3,3,*),k45(3,3,*),k46(3,3,*) ,k47(3,3,*) ,
76 . k48(3,3,*),k55(3,3,*),k56(3,3,*),k57(3,3,*) ,k58(3,3,*) ,
77 . k66(3,3,*),k67(3,3,*),k68(3,3,*),k77(3,3,*) ,k78(3,3,*) ,
78 . k88(3,3,*),vol(*),sig(nel,6)
79
80
81
82 INTEGER I,J
83
85 . fxx(mvsiz),fyy(mvsiz),fzz(mvsiz),fxy(mvsiz),fxz(mvsiz),
86 . fyz(mvsiz)
87
88 DO i=1,nel
89 fxx(i)=sig(i,1)*vol(i)
90 fyy(i)=sig(i,2)*vol(i)
91 fzz(i)=sig(i,3)*vol(i)
92 fxy(i)=sig(i,4)*vol(i)
93 fyz(i)=sig(i,5)*vol(i)
94 fxz(i)=sig(i,6)*vol(i)
95 ENDDO
96
98 . px1, py1 ,pz1 ,px1, py1, pz1, k11 ,
99 . fxx, fyy ,fzz ,fxy, fyz, fxz)
101 . px1, py1 ,pz1 ,px7, py7, pz7, k17 ,
102 . fxx, fyy ,fzz ,fxy, fyz, fxz)
104 . px7, py7, pz7 ,px7, py7, pz7, k77 ,
105 . fxx, fyy ,fzz ,fxy, fyz, fxz)
106
108 . px1, py1 ,pz1 ,px2, py2, pz2, k12 ,
109 . fxx, fyy ,fzz ,fxy, fyz, fxz)
111 . px1, py1 ,pz1 ,px8, py8, pz8, k18 ,
112 . fxx, fyy ,fzz ,fxy, fyz, fxz)
114 . px2, py2 ,pz2 ,px7, py7, pz7, k27 ,
115 . fxx, fyy ,fzz ,fxy, fyz, fxz)
117 . px7, py7 ,pz7 ,px8, py8, pz8, k78 ,
118 . fxx, fyy ,fzz ,fxy, fyz, fxz)
119
121 . px1, py1 ,pz1 ,px3, py3, pz3, k13 ,
122 . fxx, fyy ,fzz ,fxy, fyz, fxz)
124 . px1, py1 ,pz1 ,px5, py5, pz5, k15 ,
125 . fxx, fyy ,fzz ,fxy, fyz, fxz)
127 . px3, py3 ,pz3 ,px7, py7, pz7, k37 ,
128 . fxx, fyy ,fzz ,fxy, fyz, fxz)
130 . px5, py5 ,pz5 ,px7, py7, pz7, k57 ,
131 . fxx, fyy ,fzz ,fxy, fyz, fxz)
132
134 . px1, py1 ,pz1 ,px4, py4, pz4, k14 ,
135 . fxx, fyy ,fzz ,fxy, fyz, fxz)
137 . px1, py1 ,pz1 ,px6, py6, pz6, k16 ,
138 . fxx, fyy ,fzz ,fxy, fyz, fxz)
140 . px4, py4 ,pz4 ,px7, py7, pz7, k47 ,
141 . fxx, fyy ,fzz ,fxy, fyz, fxz)
143 . px6, py6 ,pz6 ,px7, py7, pz7, k67 ,
144 . fxx, fyy ,fzz ,fxy, fyz, fxz)
145
147 . px2, py2 ,pz2 ,px2, py2, pz2, k22 ,
148 . fxx, fyy ,fzz ,fxy, fyz, fxz)
150 . px2, py2 ,pz2 ,px8, py8, pz8, k28 ,
151 . fxx, fyy ,fzz ,fxy, fyz, fxz)
153 . px8, py8 ,pz8 ,px8, py8, pz8, k88 ,
154 . fxx, fyy ,fzz ,fxy, fyz, fxz)
155
157 . px2, py2 ,pz2 ,px3, py3, pz3, k23 ,
158 . fxx, fyy ,fzz ,fxy, fyz, fxz)
160 . px2, py2 ,pz2 ,px5, py5, pz5, k25 ,
161 . fxx, fyy ,fzz ,fxy, fyz, fxz)
163 . px3, py3 ,pz3 ,px8, py8, pz8, k38 ,
164 . fxx, fyy ,fzz ,fxy, fyz, fxz)
166 . px5, py5 ,pz5 ,px8, py8, pz8, k58 ,
167 . fxx, fyy ,fzz ,fxy, fyz, fxz)
168
170 . px2, py2 ,pz2 ,px4, py4, pz4, k24 ,
171 . fxx, fyy ,fzz ,fxy, fyz, fxz)
173 . px2, py2 ,pz2 ,px6, py6, pz6, k26 ,
174 . fxx, fyy ,fzz ,fxy, fyz, fxz)
176 . px4, py4 ,pz4 ,px8, py8, pz8, k48 ,
177 . fxx, fyy ,fzz ,fxy, fyz, fxz)
179 . px6, py6 ,pz6 ,px8, py8, pz8, k68 ,
180 . fxx, fyy ,fzz ,fxy, fyz, fxz)
181
183 . px3, py3 ,pz3 ,px3, py3, pz3, k33 ,
184 . fxx, fyy ,fzz ,fxy, fyz, fxz)
186 . px3, py3 ,pz3 ,px5, py5, pz5, k35 ,
187 . fxx, fyy ,fzz ,fxy, fyz, fxz)
189 . px5, py5 ,pz5 ,px5, py5, pz5, k55 ,
190 . fxx, fyy ,fzz ,fxy, fyz, fxz)
191
193 . px3, py3 ,pz3 ,px4, py4, pz4, k34 ,
194 . fxx, fyy ,fzz ,fxy, fyz, fxz)
196 . px3, py3 ,pz3 ,px6, py6, pz6, k36 ,
197 . fxx, fyy ,fzz ,fxy, fyz, fxz)
199 . px5, py5 ,pz5 ,px6, py6, pz6, k56 ,
200 . fxx, fyy ,fzz ,fxy, fyz, fxz)
202 . px4, py4 ,pz4 ,px5, py5, pz5, k45 ,
203 . fxx, fyy ,fzz ,fxy, fyz, fxz)
204
206 . px4, py4 ,pz4 ,px4, py4, pz4, k44 ,
207 . fxx, fyy ,fzz ,fxy, fyz, fxz)
209 . px4, py4 ,pz4 ,px6, py6, pz6, k46 ,
210 . fxx, fyy ,fzz ,fxy, fyz, fxz)
212 . px6, py6 ,pz6 ,px6, py6, pz6, k66 ,
213 . fxx, fyy ,fzz ,fxy, fyz, fxz)
214
215 RETURN
subroutine s8zkehij(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, kij, fxx, fyy, fzz, fxy, fyz, fxz)