47
48
49
50#include "implicit_f.inc"
51
52
53
54#include "mvsiz_p.inc"
55
56
57
58
59
60
61 INTEGER NEL
62
64 . px(mvsiz,10),py(mvsiz,10),pz(mvsiz,10),
65 . vol(*),sig(nel,6)
67 . k11(9,*) ,k12(9,*) ,k13(9,*) ,k14(9,*) ,k15(9,*) ,
68 . k16(9,*) ,k17(9,*) ,k18(9,*) ,k19(9,*) ,k10(9,*) ,
69 . k22(9,*) ,k23(9,*) ,k24(9,*) ,k25(9,*) ,k26(9,*) ,
70 . k27(9,*) ,k28(9,*) ,k29(9,*) ,k20(9,*) ,k33(9,*) ,
71 . k34(9,*) ,k35(9,*) ,k36(9,*) ,k37(9,*) ,k38(9,*) ,
72 . k39(9,*) ,k30(9,*) ,k44(9,*) ,k45(9,*) ,k46(9,*) ,
73 . k47(9,*) ,k48(9,*) ,k49(9,*) ,k40(9,*) ,k55(9,*) ,
74 . k56(9,*) ,k57(9,*) ,k58(9,*) ,k59(9,*) ,k50(9,*) ,
75 . k66(9,*) ,k67(9,*) ,k68(9,*) ,k69(9,*) ,k60(9,*) ,
76 . k77(9,*) ,k78(9,*) ,k79(9,*) ,k70(9,*) ,k88(9,*) ,
77 . k89(9,*) ,k80(9,*) ,k99(9,*) ,k90(9,*) ,k00(9,*)
78
79
80
81 INTEGER ,J
83 . fxx(mvsiz),fyy(mvsiz),fzz(mvsiz),fxy(mvsiz),fxz(mvsiz),
84 . fyz(mvsiz)
85
86 DO i=1,nel
87 fxx(i)=sig(i,1)*vol(i)
88 fyy(i)=sig(i,2)*vol(i)
89 fzz(i)=sig(i,3)*vol(i)
90 fxy(i)=sig(i,4)*vol(i)
91 fyz(i)=sig(i,5)*vol(i)
92 fxz(i)=sig(i,6)*vol(i)
93 ENDDO
94
96 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
97 . fxx, fyy ,fzz ,fxy, fyz, fxz,k11 )
99 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
100 . fxx, fyy ,fzz ,fxy, fyz, fxz,k22 )
102 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
103 . fxx, fyy ,fzz ,fxy, fyz, fxz,k33 )
105 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
106 . fxx, fyy ,fzz ,fxy, fyz, fxz,k44 )
108 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
109 . fxx, fyy ,fzz ,fxy, fyz, fxz,k55 )
111 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
112 . fxx, fyy ,fzz ,fxy, fyz, fxz,k66 )
114 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
115 . fxx, fyy ,fzz ,fxy, fyz, fxz,k77 )
117 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
118 . fxx, fyy ,fzz ,fxy, fyz, fxz,k88 )
120 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
121 . fxx, fyy ,fzz ,fxy, fyz, fxz,k99 )
123 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
124 . fxx, fyy ,fzz ,fxy, fyz, fxz,k00 )
125
127 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
128 . fxx, fyy ,fzz ,fxy, fyz, fxz,k12 )
130 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
131 . fxx, fyy ,fzz ,fxy, fyz, fxz,k13 )
133 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
134 . fxx, fyy ,fzz ,fxy, fyz, fxz,k14 )
136 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
137 . fxx, fyy ,fzz ,fxy, fyz, fxz,k15 )
139 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
140 . fxx, fyy ,fzz ,fxy, fyz, fxz,k16 )
142 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
143 . fxx, fyy ,fzz ,fxy, fyz, fxz,k17 )
145 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
146 . fxx, fyy ,fzz ,fxy, fyz, fxz,k18 )
148 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
149 . fxx, fyy ,fzz ,fxy, fyz, fxz,k19 )
151 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
152 . fxx, fyy ,fzz ,fxy, fyz, fxz,k10 )
154 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
155 . fxx, fyy ,fzz ,fxy, fyz, fxz,k23 )
157 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
158 . fxx, fyy ,fzz ,fxy, fyz, fxz,k24 )
160 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
161 . fxx, fyy ,fzz ,fxy, fyz, fxz,k25 )
163 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
164 . fxx, fyy ,fzz ,fxy, fyz, fxz,k26 )
166 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
167 . fxx, fyy ,fzz ,fxy, fyz, fxz,k27 )
169 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
170 . fxx, fyy ,fzz ,fxy, fyz, fxz,k28 )
172 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
173 . fxx, fyy ,fzz ,fxy, fyz, fxz,k29 )
175 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
176 . fxx, fyy ,fzz ,fxy, fyz, fxz,k20 )
178 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
179 . fxx, fyy ,fzz ,fxy, fyz, fxz,k34 )
181 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
182 . fxx, fyy ,fzz ,fxy, fyz, fxz,k35 )
184 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
185 . fxx, fyy ,fzz ,fxy, fyz, fxz,k36 )
187 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
188 . fxx, fyy ,fzz ,fxy, fyz, fxz,k37 )
190 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
191 . fxx, fyy ,fzz ,fxy, fyz, fxz,k38 )
193 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
194 . fxx, fyy ,fzz ,fxy, fyz, fxz,k39 )
196 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
197 . fxx, fyy ,fzz ,fxy, fyz, fxz,k30 )
199 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
200 . fxx, fyy ,fzz ,fxy, fyz, fxz,k45 )
202 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
203 . fxx, fyy ,fzz ,fxy, fyz, fxz,k46 )
205 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
206 . fxx, fyy ,fzz ,fxy, fyz, fxz,k47 )
208 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
209 . fxx, fyy ,fzz ,fxy, fyz, fxz,k48 )
211 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
212 . fxx, fyy ,fzz ,fxy, fyz, fxz,k49 )
214 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
215 . fxx, fyy ,fzz ,fxy, fyz, fxz,k40 )
217 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
218 . fxx, fyy ,fzz ,fxy, fyz, fxz,k56 )
220 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
221 . fxx, fyy ,fzz ,fxy, fyz, fxz,k57 )
223 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
224 . fxx, fyy ,fzz ,fxy, fyz, fxz,k58 )
226 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
227 . fxx, fyy ,fzz ,fxy, fyz, fxz,k59 )
229 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
230 . fxx, fyy ,fzz ,fxy, fyz, fxz,k50 )
232 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
233 . fxx, fyy ,fzz ,fxy, fyz, fxz,k67 )
235 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
236 . fxx, fyy ,fzz ,fxy, fyz, fxz,k68 )
238 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
239 . fxx, fyy ,fzz ,fxy, fyz, fxz,k69 )
241 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
242 . fxx, fyy ,fzz ,fxy, fyz, fxz,k60 )
244 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
245 . fxx, fyy ,fzz ,fxy, fyz, fxz,k78 )
247 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
248 . fxx, fyy ,fzz ,fxy, fyz, fxz,k79 )
250 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
251 . fxx, fyy ,fzz ,fxy, fyz, fxz,k70 )
253 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
254 . fxx, fyy ,fzz ,fxy, fyz, fxz,k89 )
256 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
257 . fxx, fyy ,fzz ,fxy, fyz, fxz,k80 )
259 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
260 . fxx, fyy ,fzz ,fxy, fyz, fxz,k90 )
261
262 RETURN
subroutine s4keh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, fxx, fyy, fzz, fxy, fyz, fxz, kij)