45
46
47
48#include "implicit_f.inc"
49
50
51
52#include "mvsiz_p.inc"
53
54
55
56
57
58
59 INTEGER NEL
60
62 . px1(*), px2(*), px3(*), px4(*),
63 . py1(*), py2(*), py3(*), py4(*),
64 . pz1(*), pz2(*), pz3(*), pz4(*),
65 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*) ,k15(3,3,*) ,
66 . k16(3,3,*),k17(3,3,*),k18(3,3,*),k22(3,3,*) ,k23(3,3,*) ,
67 . k24(3,3,*),k25(3,3,*),k26(3,3,*),k27(3,3,*) ,k28(3,3,*) ,
68 . k33(3,3,*),k34(3,3,*),k35(3,3,*),k36(3,3,*) ,k37(3,3,*) ,
69 . k38(3,3,*),k44(3,3,*),k45(3,3,*),k46(3,3,*) ,k47(3,3,*) ,
70 . k48(3,3,*),k55(3,3,*),k56(3,3,*),k57(3,3,*) ,k58(3,3,*) ,
71 . k66(3,3,*),k67(3,3,*),k68(3,3,*),k77(3,3,*) ,k78(3,3,*) ,
72 . k88(3,3,*),vol(*),sig(nel,6)
73
74
75
76 INTEGER I,J
77
79 . fxx(mvsiz),fyy(mvsiz),fzz(mvsiz),fxy(mvsiz),fxz(mvsiz),
80 . fyz(mvsiz),hij(mvsiz)
81
82 DO i=1,nel
83 fxx(i)=sig(i,1)*vol(i)
84 fyy(i)=sig(i,2)*vol(i)
85 fzz(i)=sig(i,3)*vol(i)
86 fxy(i)=sig(i,4)*vol(i)
87 fyz(i)=sig(i,5)*vol(i)
88 fxz(i)=sig(i,6)*vol(i)
89 ENDDO
90
92 . px1, py1 ,pz1 ,px1, py1, pz1, hij ,
93 . fxx, fyy ,fzz ,fxy, fyz, fxz)
94 DO j = 1,3
95 DO i=1,nel
96 k11(j,j,i)=k11(j,j,i)+hij(i)
97 k77(j,j,i)=k77(j,j,i)+hij(i)
98 k17(j,j,i)=k17(j,j,i)-hij(i)
99 ENDDO
100 ENDDO
101
103 . px1, py1 ,pz1 ,px2, py2, pz2, hij ,
104 . fxx, fyy ,fzz ,fxy, fyz, fxz)
105 DO j = 1,3
106 DO i=1,nel
107 k12(j,j,i)=k12(j,j,i)+hij(i)
108 k18(j,j,i)=k18(j,j,i)-hij(i)
109 k27(j,j,i)=k27(j,j,i)-hij(i)
110 k78(j,j,i)=k78(j,j,i)+hij(i)
111 ENDDO
112 ENDDO
113
115 . px1, py1 ,pz1 ,px3, py3, pz3, hij ,
116 . fxx, fyy ,fzz ,fxy, fyz, fxz)
117 DO j = 1,3
118 DO i=1,nel
119 k13(j,j,i)=k13(j,j,i)+hij(i)
120 k15(j,j,i)=k15(j,j,i)-hij(i)
121 k37(j,j,i)=k37(j,j,i)-hij(i)
122 k57(j,j,i)=k57(j,j,i)+hij(i)
123 ENDDO
124 ENDDO
125
127 . px1, py1 ,pz1 ,px4, py4, pz4, hij ,
128 . fxx, fyy ,fzz ,fxy, fyz, fxz)
129 DO j = 1,3
130 DO i=1,nel
131 k14(j,j,i)=k14(j,j,i)+hij(i)
132 k16(j,j,i)=k16(j,j,i)-hij(i)
133 k47(j,j,i)=k47(j,j,i)-hij(i)
134 k67(j,j,i)=k67(j,j,i)+hij(i)
135 ENDDO
136 ENDDO
137
139 . px2, py2 ,pz2 ,px2, py2, pz2, hij ,
140 . fxx, fyy ,fzz ,fxy, fyz, fxz)
141 DO j = 1,3
142 DO i=1,nel
143 k22(j,j,i)=k22(j,j,i)+hij(i)
144 k88(j,j,i)=k88(j,j,i)+hij(i)
145 k28(j,j,i)=k28(j,j,i)-hij(i)
146 ENDDO
147 ENDDO
148
150 . px2, py2 ,pz2 ,px3, py3, pz3, hij ,
151 . fxx, fyy ,fzz ,fxy, fyz, fxz)
152 DO j = 1,3
153 DO i=1,nel
154 k23(j,j,i)=k23(j,j,i)+hij(i)
155 k25(j,j,i)=k25(j,j,i)-hij(i)
156 k38(j,j,i)=k38(j,j,i)-hij(i)
157 k58(j,j,i)=k58(j,j,i)+hij(i)
158 ENDDO
159 ENDDO
160
162 . px2, py2 ,pz2 ,px4, py4, pz4, hij ,
163 . fxx, fyy ,fzz ,fxy, fyz, fxz)
164 DO j = 1,3
165 DO i=1,nel
166 k24(j,j,i)=k24(j,j,i)+hij(i)
167 k26(j,j,i)=k26(j,j,i)-hij(i)
168 k48(j,j,i)=k48(j,j,i)-hij(i)
169 k68(j,j,i)=k68(j,j,i)+hij(i)
170 ENDDO
171 ENDDO
172
174 . px3, py3 ,pz3 ,px3, py3, pz3, hij ,
175 . fxx, fyy ,fzz ,fxy, fyz, fxz)
176 DO j = 1,3
177 DO i=1,nel
178 k33(j,j,i)=k33(j,j,i)+hij(i)
179 k55(j,j,i)=k55(j,j,i)+hij(i)
180 k35(j,j,i)=k35(j,j,i)-hij(i)
181 ENDDO
182 ENDDO
183
185 . px3, py3 ,pz3 ,px4, py4, pz4, hij ,
186 . fxx, fyy ,fzz ,fxy, fyz, fxz)
187 DO j = 1,3
188 DO i=1,nel
189 k34(j,j,i)=k34(j,j,i)+hij(i)
190 k36(j,j,i)=k36(j,j,i)-hij(i)
191 k56(j,j,i)=k56(j,j,i)+hij(i)
192 k45(j,j,i)=k45(j,j,i)-hij(i)
193 ENDDO
194 ENDDO
195
197 . px4, py4 ,pz4 ,px4, py4, pz4, hij ,
198 . fxx, fyy ,fzz ,fxy, fyz, fxz)
199 DO j = 1,3
200 DO i=1,nel
201 k44(j,j,i)=k44(j,j,i)+hij(i)
202 k66(j,j,i)=k66(j,j,i)+hij(i)
203 k46(j,j,i)=k46(j,j,i)-hij(i)
204 ENDDO
205 ENDDO
206
207 RETURN
subroutine s8zkeh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, hij, fxx, fyy, fzz, fxy, fyz, fxz)