OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
czlkec3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com04_c.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine czlkec3 (jft, jlt, vol, thk0, thk2, hm, hf, hz, a_i, z1, px1, px2, py1, py2, nplat, iplat, dhz, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, iorth, hmor, hfor, hmfor)
subroutine czlkecr3 (jft, jlt, vol, thk0, thk2, hm, hf, hz, a_i, z1, px1, px2, py1, py2, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, iorth, hmor, hfor, iplat, dhz, prx, pry, prxy, prz, hmfor, nplat)

Function/Subroutine Documentation

◆ czlkec3()

subroutine czlkec3 ( integer jft,
integer jlt,
vol,
thk0,
thk2,
hm,
hf,
hz,
a_i,
z1,
px1,
px2,
py1,
py2,
integer nplat,
integer, dimension(*) iplat,
dhz,
k11,
k12,
k13,
k14,
k22,
k23,
k24,
k33,
k34,
k44,
m11,
m12,
m13,
m14,
m22,
m23,
m24,
m33,
m34,
m44,
mf11,
mf12,
mf13,
mf14,
mf22,
mf23,
mf24,
mf33,
mf34,
mf44,
fm12,
fm13,
fm14,
fm23,
fm24,
fm34,
integer iorth,
hmor,
hfor,
hmfor )

Definition at line 28 of file czlkec3.F.

37C--------------------------------------------------------------------------------------------------
38C-----------------------------------------------
39C I M P L I C I T T Y P E S
40C-----------------------------------------------
41#include "implicit_f.inc"
42#include "mvsiz_p.inc"
43#include "com04_c.inc"
44C-----------------------------------------------
45C D U M M Y A R G U M E N T S
46C-----------------------------------------------
47C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
48 INTEGER JFT,JLT,NPLAT,IPLAT(*),IKGEO,IORTH
49 my_real
50 . vol(*),thk0(*),thk2(*),a_i(*),z1(*),
51 . hm(mvsiz,4),hf(mvsiz,4),hz(*),
52 . px1(*) ,px2(*) ,py1(*) ,py2(*) ,
53 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*),
54 . k22(3,3,*),k23(3,3,*),k24(3,3,*),k33(3,3,*),
55 . m11(3,3,*),m12(3,3,*),m13(3,3,*),m14(3,3,*),
56 . m22(3,3,*),m23(3,3,*),m24(3,3,*),m33(3,3,*),
57 . mf11(3,3,*),mf12(3,3,*),mf13(3,3,*),mf14(3,3,*),
58 . mf22(3,3,*),mf23(3,3,*),mf24(3,3,*),mf33(3,3,*),
59 . fm12(3,3,*),fm13(3,3,*),fm14(3,3,*),
60 . fm23(3,3,*),fm24(3,3,*),fm34(3,3,*),
61 . k34(3,3,*),k44(3,3,*),m34(3,3,*),m44(3,3,*),
62 . mf34(3,3,*),mf44(3,3,*),dhz(*),hmor(mvsiz,2),hfor(mvsiz,2),
63 . hmfor(mvsiz,6)
64C---------------|[KIJ][MFIJ]|----
65C-----KE(6x6)= | |
66C---------------|[FMIJ]{MIJ]|----
67C-----------------------------------------------
68C L O C A L V A R I A B L E S
69C-----------------------------------------------
70 INTEGER EP,I,J,NF,M
71 my_real
72 . dm(2,2,mvsiz),df(2,2,mvsiz),
73 . c1,c2,gm(mvsiz),gf(mvsiz),dg(mvsiz),
74 . px1px1(mvsiz),px1px2(mvsiz) ,px2px2(mvsiz),
75 . px1py1(mvsiz),px1py2(mvsiz) ,px2py1(mvsiz),
76 . px2py2(mvsiz),py1py1(mvsiz),py1py2(mvsiz),py2py2(mvsiz),
77 . facz(mvsiz),facz2(mvsiz),dfz,dm5(mvsiz),dm6(mvsiz),df5(mvsiz),
78 . df6(mvsiz),dm5_2(mvsiz),dm6_2(mvsiz),df5_2(mvsiz),
79 . df6_2(mvsiz),dmf(3,3,mvsiz),
80 . bxibxj(2,2,mvsiz),bxibyj(2,2,mvsiz),byibxj(2,2,mvsiz),
81 . byibyj(2,2,mvsiz),fac2z
82C-----------Attention Matrice sym Kii ne calcul que la moitie---------72
83 nf=nplat+1
84#include "vectorize.inc"
85C--------Constante parts---------------
86 DO m=jft,jlt
87 ep=iplat(m)
88 c2=vol(ep)
89 c1=thk2(ep)*c2
90 dm(1,1,m)=hm(ep,1)*c2
91 dm(2,2,m)=hm(ep,2)*c2
92 dm(1,2,m)=hm(ep,3)*c2
93 dm(2,1,m)=dm(1,2,m)
94 gm(m) =hm(ep,4)*c2
95 df(1,1,m)=hf(ep,1)*c1
96 df(2,2,m)=hf(ep,2)*c1
97 df(1,2,m)=hf(ep,3)*c1
98 df(2,1,m)=df(1,2,m)
99 gf(m) =hf(ep,4)*c1
100 dhz(m)=hz(ep)*c1
101 ENDDO
102C---------px3=-px1,px4=-px2-------------
103 DO ep=jft,jlt
104 px1px1(ep) = px1(ep)*px1(ep)
105 px1px2(ep) = px1(ep)*px2(ep)
106 px2px2(ep) = px2(ep)*px2(ep)
107 px1py1(ep) = px1(ep)*py1(ep)
108 px1py2(ep) = px1(ep)*py2(ep)
109 px2py1(ep) = px2(ep)*py1(ep)
110 px2py2(ep) = px2(ep)*py2(ep)
111 py1py1(ep) = py1(ep)*py1(ep)
112 py1py2(ep) = py1(ep)*py2(ep)
113 py2py2(ep) = py2(ep)*py2(ep)
114 ENDDO
115C
116 DO ep=jft,jlt
117 k11(1,1,ep)=dm(1,1,ep)*px1px1(ep)+gm(ep)*py1py1(ep)
118 k12(1,1,ep)=dm(1,1,ep)*px1px2(ep)+gm(ep)*py1py2(ep)
119 k22(1,1,ep)=dm(1,1,ep)*px2px2(ep)+gm(ep)*py2py2(ep)
120C
121 m11(2,2,ep)=df(1,1,ep)*px1px1(ep)+gf(ep)*py1py1(ep)
122 m12(2,2,ep)=df(1,1,ep)*px1px2(ep)+gf(ep)*py1py2(ep)
123 m22(2,2,ep)=df(1,1,ep)*px2px2(ep)+gf(ep)*py2py2(ep)
124 ENDDO
125C
126 DO ep=jft,jlt
127 k11(1,2,ep)=(dm(1,2,ep)+gm(ep))*px1py1(ep)
128 k12(1,2,ep)=dm(1,2,ep)*px1py2(ep)+gm(ep)*px2py1(ep)
129 k22(1,2,ep)=(dm(1,2,ep)+gm(ep))*px2py2(ep)
130 k12(2,1,ep)=dm(1,2,ep)*px2py1(ep)+gm(ep)*px1py2(ep)
131C
132 m11(1,2,ep)=-(df(1,2,ep)+gf(ep))*px1py1(ep)
133 m12(1,2,ep)=-df(1,2,ep)*px2py1(ep)-gf(ep)*px1py2(ep)
134 m22(1,2,ep)=-df(1,2,ep)*px2py2(ep)-gf(ep)*px2py2(ep)
135 m12(2,1,ep)=-df(1,2,ep)*px1py2(ep)-gf(ep)*px2py1(ep)
136 ENDDO
137C
138 DO ep=jft,jlt
139 k11(2,2,ep)=dm(2,2,ep)*py1py1(ep)+gm(ep)*px1px1(ep)
140 k12(2,2,ep)=dm(2,2,ep)*py1py2(ep)+gm(ep)*px1px2(ep)
141 k22(2,2,ep)=dm(2,2,ep)*py2py2(ep)+gm(ep)*px2px2(ep)
142C
143 m11(1,1,ep)=df(2,2,ep)*py1py1(ep)+gf(ep)*px1px1(ep)
144 m12(1,1,ep)=df(2,2,ep)*py1py2(ep)+gf(ep)*px1px2(ep)
145 m22(1,1,ep)=df(2,2,ep)*py2py2(ep)+gf(ep)*px2px2(ep)
146 ENDDO
147 IF (iorth >0 ) THEN
148#include "vectorize.inc"
149 DO m=jft,jlt
150 ep=iplat(m)
151 c2=vol(ep)
152 c1=thk2(ep)*c2
153 dm5(m)=hmor(ep,1)*c2
154 dm6(m)=hmor(ep,2)*c2
155 df5(m)=hfor(ep,1)*c1
156 df6(m)=hfor(ep,2)*c1
157 ENDDO
158 DO m=jft,jlt
159 dm5_2(m)=two*dm5(m)
160 dm6_2(m)=two*dm6(m)
161 df5_2(m)=two*df5(m)
162 df6_2(m)=two*df6(m)
163 ENDDO
164 DO ep=jft,jlt
165 k11(1,1,ep)=k11(1,1,ep)+px1py1(ep)*dm5_2(ep)
166 k11(1,2,ep)=k11(1,2,ep)+
167 1 px1px1(ep)*dm5(ep)+py1py1(ep)*dm6(ep)
168 k11(2,2,ep)=k11(2,2,ep)+px1py1(ep)*dm6_2(ep)
169 k12(1,1,ep)=k12(1,1,ep)+(px1py2(ep)+px2py1(ep))*dm5(ep)
170 c1 = px1px2(ep)*dm5(ep)+py1py2(ep)*dm6(ep)
171 k12(1,2,ep)=k12(1,2,ep)+c1
172 k12(2,1,ep)=k12(2,1,ep)+c1
173 k12(2,2,ep)=k12(2,2,ep)+(px1py2(ep)+px2py1(ep))*dm6(ep)
174 k22(1,1,ep)=k22(1,1,ep)+px2py2(ep)*dm5_2(ep)
175 k22(1,2,ep)=k22(1,2,ep)+
176 1 px2px2(ep)*dm5(ep)+py2py2(ep)*dm6(ep)
177 k22(2,2,ep)=k22(2,2,ep)+px2py2(ep)*dm6_2(ep)
178C
179 m11(1,1,ep)=m11(1,1,ep)+px1py1(ep)*df6_2(ep)
180 m11(1,2,ep)=m11(1,2,ep)-
181 1 px1px1(ep)*df5(ep)-py1py1(ep)*df6(ep)
182 m11(2,2,ep)=m11(2,2,ep)+px1py1(ep)*df5_2(ep)
183 m12(1,1,ep)=m12(1,1,ep)+(px1py2(ep)+px2py1(ep))*df6(ep)
184 c2 = px1px2(ep)*df5(ep)+py1py2(ep)*df6(ep)
185 m12(1,2,ep)=m12(1,2,ep)-c2
186 m12(2,1,ep)=m12(2,1,ep)-c2
187 m12(2,2,ep)=m12(2,2,ep)+(px1py2(ep)+px2py1(ep))*df5(ep)
188 m22(1,1,ep)=m22(1,1,ep)+px2py2(ep)*df6_2(ep)
189 m22(1,2,ep)=m22(1,2,ep)-
190 1 px2px2(ep)*df5(ep)-py2py2(ep)*df6(ep)
191 m22(2,2,ep)=m22(2,2,ep)+px2py2(ep)*df5_2(ep)
192 ENDDO
193 DO m=jft,jlt
194 ep=iplat(m)
195 c2=vol(ep)*thk0(ep)
196 dmf(1,1,m)=hmfor(ep,1)*c2
197 dmf(2,2,m)=hmfor(ep,2)*c2
198 dmf(1,2,m)=hmfor(ep,3)*c2
199 dmf(1,3,m)=hmfor(ep,5)*c2
200 dmf(2,3,m)=hmfor(ep,6)*c2
201 dmf(2,1,m)=dmf(1,2,m)
202 dmf(3,1,m)=dmf(1,3,m)
203 dmf(3,2,m)=dmf(2,3,m)
204 dmf(3,3,m)=hmfor(ep,4)*c2
205 ENDDO
206C ----add mem-bending coupling terms of orthotrope--coplanar first---
207 DO ep=jft,nplat
208 mf11(1,1,ep)=
209 1 -dmf(1,2,ep)*px1py1(ep)-dmf(1,3,ep)*px1px1(ep)
210 2 -dmf(2,3,ep)*py1py1(ep)-dmf(3,3,ep)*px1py1(ep)
211 mf11(1,2,ep)=
212 1 dmf(1,1,ep)*px1px1(ep)+dmf(1,3,ep)*px1py1(ep)
213 2 +dmf(1,3,ep)*px1py1(ep)+dmf(3,3,ep)*py1py1(ep)
214 mf11(2,1,ep)=
215 1 -dmf(2,2,ep)*py1py1(ep)-dmf(2,3,ep)*px1py1(ep)
216 2 -dmf(2,3,ep)*px1py1(ep)-dmf(3,3,ep)*px1px1(ep)
217 mf11(2,2,ep)=
218 1 dmf(1,2,ep)*px1py1(ep)+dmf(2,3,ep)*py1py1(ep)
219 2 +dmf(1,3,ep)*px1px1(ep)+dmf(3,3,ep)*px1py1(ep)
220 ENDDO
221 DO ep=jft,nplat
222 mf12(1,1,ep)=
223 1 -dmf(1,2,ep)*px1py2(ep)-dmf(1,3,ep)*px1px2(ep)
224C--------------------------------------------------------YIXJ-----
225 2 -dmf(2,3,ep)*py1py2(ep)-dmf(3,3,ep)*px2py1(ep)
226 mf12(1,2,ep)=
227 1 dmf(1,1,ep)*px1px2(ep)+dmf(1,3,ep)*px1py2(ep)
228C-------------------------YIXJ-----
229 2 +dmf(1,3,ep)*px2py1(ep)+dmf(3,3,ep)*py1py2(ep)
230 mf12(2,1,ep)=
231C---------------------------------------------------------YIXJ-----
232 1 -dmf(2,2,ep)*py1py2(ep)-dmf(2,3,ep)*px2py1(ep)
233 2 -dmf(2,3,ep)*px1py2(ep)-dmf(3,3,ep)*px1px2(ep)
234 mf12(2,2,ep)=
235C---------------------------------YIXJ-----
236 1 dmf(1,2,ep)*px2py1(ep)+dmf(2,3,ep)*py1py2(ep)
237 2 +dmf(1,3,ep)*px1px2(ep)+dmf(3,3,ep)*px1py2(ep)
238 ENDDO
239 DO ep=jft,nplat
240 mf22(1,1,ep)=
241 1 -dmf(1,2,ep)*px2py2(ep)-dmf(1,3,ep)*px2px2(ep)
242C--------------------------------------------------------YIXJ-----
243 2 -dmf(2,3,ep)*py2py2(ep)-dmf(3,3,ep)*px2py2(ep)
244 mf22(1,2,ep)=
245 1 dmf(1,1,ep)*px2px2(ep)+dmf(1,3,ep)*px2py2(ep)
246C-------------------------YIXJ-----
247 2 +dmf(1,3,ep)*px2py2(ep)+dmf(3,3,ep)*py2py2(ep)
248 mf22(2,1,ep)=
249C---------------------------------------------------------YIXJ-----
250 1 -dmf(2,2,ep)*py2py2(ep)-dmf(2,3,ep)*px2py2(ep)
251 2 -dmf(2,3,ep)*px2py2(ep)-dmf(3,3,ep)*px2px2(ep)
252 mf22(2,2,ep)=
253C---------------------------------YIXJ-----
254 1 dmf(1,2,ep)*px2py2(ep)+dmf(2,3,ep)*py2py2(ep)
255 2 +dmf(1,3,ep)*px2px2(ep)+dmf(3,3,ep)*px2py2(ep)
256 ENDDO
257 DO ep=jft,nplat
258 mf23(1,1,ep)=
259 1 dmf(1,2,ep)*px2py1(ep)+dmf(1,3,ep)*px1px2(ep)
260C--------------------------------------------------------YIXJ-----
261 2 +dmf(2,3,ep)*py1py2(ep)+dmf(3,3,ep)*px1py2(ep)
262 mf23(1,2,ep)=
263 1 -dmf(1,1,ep)*px1px2(ep)-dmf(1,3,ep)*px2py1(ep)
264C-------------------------YIXJ-----
265 2 -dmf(1,3,ep)*px1py2(ep)-dmf(3,3,ep)*py1py2(ep)
266 mf23(2,1,ep)=
267C---------------------------------------------------------YIXJ-----
268 1 dmf(2,2,ep)*py1py2(ep)+dmf(2,3,ep)*px1py2(ep)
269 2 +dmf(2,3,ep)*px2py1(ep)+dmf(3,3,ep)*px1px2(ep)
270 mf23(2,2,ep)=
271C---------------------------------YIXJ-----
272 1 -dmf(1,2,ep)*px1py2(ep)-dmf(2,3,ep)*py1py2(ep)
273 2 -dmf(1,3,ep)*px1px2(ep)-dmf(3,3,ep)*px2py1(ep)
274 ENDDO
275 DO i=1,2
276 DO j=1,2
277 DO ep=jft,nplat
278 mf13(i,j,ep)=-mf11(i,j,ep)
279 mf14(i,j,ep)=-mf12(i,j,ep)
280 mf24(i,j,ep)=-mf22(i,j,ep)
281 mf33(i,j,ep)=mf11(i,j,ep)
282 mf34(i,j,ep)=mf12(i,j,ep)
283 mf44(i,j,ep)=mf22(i,j,ep)
284 ENDDO
285 ENDDO
286 ENDDO
287C
288 DO i=1,2
289 DO j=1,2
290 DO ep=jft,nplat
291 fm12(i,j,ep)=-mf23(j,i,ep)
292 fm13(i,j,ep)=-mf11(j,i,ep)
293 fm14(i,j,ep)= mf23(j,i,ep)
294 fm23(i,j,ep)=-mf12(j,i,ep)
295 fm24(i,j,ep)=-mf22(j,i,ep)
296 fm34(i,j,ep)=-mf23(j,i,ep)
297 ENDDO
298 ENDDO
299 ENDDO
300 ENDIF !(IORTH==1)
301C---+----1----+----2----+warped elements----+----5----+----6----+----7----+----8
302#include "vectorize.inc"
303 DO m=nf,jlt
304 ep=iplat(m)
305 facz(m) = four*z1(ep)*a_i(ep)
306 facz2(m) = facz(m)*facz(m)
307 dg(m)=gf(m)*facz2(m)
308 ENDDO
309C---+----1----+----2----+warped elements----+----5----+----6----+----7----+----8
310C--------px1'=py2;px2'=py1;py1'=-px2;py2'=-px1;-------------
311 DO ep=nf,jlt
312 dfz =df(1,1,ep)*facz2(ep)
313 k11(1,1,ep)=k11(1,1,ep)+dfz*py2py2(ep)+dg(ep)*px2px2(ep)
314 k12(1,1,ep)=k12(1,1,ep)+dfz*py1py2(ep)+dg(ep)*px1px2(ep)
315 k22(1,1,ep)=k22(1,1,ep)+dfz*py1py1(ep)+dg(ep)*px1px1(ep)
316 ENDDO
317C
318 DO ep=nf,jlt
319 dfz =df(1,2,ep)*facz2(ep)
320 k11(1,2,ep)=k11(1,2,ep)-(dfz+dg(ep))*px2py2(ep)
321 k12(1,2,ep)=k12(1,2,ep)-dfz*px1py2(ep)-dg(ep)*px2py1(ep)
322 k22(1,2,ep)=k22(1,2,ep)-(dfz+dg(ep))*px1py1(ep)
323 k12(2,1,ep)=k12(2,1,ep)-dfz*px2py1(ep)-dg(ep)*px1py2(ep)
324 ENDDO
325C
326 DO ep=nf,jlt
327 dfz =df(2,2,ep)*facz2(ep)
328 k11(2,2,ep)=k11(2,2,ep)+dfz*px2px2(ep)+dg(ep)*py2py2(ep)
329 k12(2,2,ep)=k12(2,2,ep)+dfz*px1px2(ep)+dg(ep)*py1py2(ep)
330 k22(2,2,ep)=k22(2,2,ep)+dfz*px1px1(ep)+dg(ep)*py1py1(ep)
331 ENDDO
332C--------(px1px1)'=py2py2;(px1px2)'=py1py2;(px2px2)'=py1py1;-------------
333C--------(py1py1)'=px2px2;(py1py2)'=px1px2;(py2py2)'=px1px1;-------------
334C--------(px1py1)'=-px2py2;(px2py1)'=-px2py1;(px2py2)'=-px1py1;(px1py2)'=-px1py2;-------------
335 IF (iorth > 0) THEN
336 DO ep=nf,jlt
337 k11(1,1,ep)=k11(1,1,ep)-px2py2(ep)*df5_2(ep)*facz2(ep)
338 k11(1,2,ep)=k11(1,2,ep)+facz2(ep)*(
339 1 py2py2(ep)*df5(ep)+px2px2(ep)*df6(ep))
340 k11(2,2,ep)=k11(2,2,ep)-px2py2(ep)*df6_2(ep)*facz2(ep)
341 k12(1,1,ep)=k12(1,1,ep)-facz2(ep)*(
342 1 px1py2(ep)+px2py1(ep))*df5(ep)
343 c1 = (py1py2(ep)*df5(ep)+px1px2(ep)*df6(ep))*facz2(ep)
344 k12(1,2,ep)=k12(1,2,ep)+c1
345 k12(2,1,ep)=k12(2,1,ep)+c1
346 k12(2,2,ep)=k12(2,2,ep)-facz2(ep)*(
347 1 px1py2(ep)+px2py1(ep))*df6(ep)
348 k22(1,1,ep)=k22(1,1,ep)-px1py1(ep)*df5_2(ep)*facz2(ep)
349 k22(1,2,ep)=k22(1,2,ep)+facz2(ep)*(
350 1 py1py1(ep)*df5(ep)+px1px1(ep)*df6(ep))
351 k22(2,2,ep)=k22(2,2,ep)-px1py1(ep)*df6_2(ep)*facz2(ep)
352 ENDDO
353C ----add mem-bending coupling terms of orthotrope--warped elements-suite---
354C--------px1px1'=px1'px1=px1py2;px1px2'=px1py1;px1'px2=px2py2;px2px2'=px2'px2=px2py1;-------------
355C--------py1py1'=py1'py1=-px2py1;py1py2'=-px1py1;py1'py2=-px2py2;py2py2'=py2'py2=-px1py2;-------------
356C--------px1py1'=-px2px1;px1'py1=py1py2;px2py1'=-px2px2;px2'py1=py1py1;px2py2'=-px1px2;px2'py2=py1py2;px1py2'=-px1px1;px1'py2=py2py2;-------------
357 DO ep=nf,jlt
358 fac2z=two*facz(ep)
359 bxibxj(1,1,ep)=fac2z*px1py2(ep)
360 bxibxj(1,2,ep)=facz(ep)*(px1py1(ep)+px2py2(ep))
361 bxibxj(2,1,ep)=bxibxj(1,2,ep)
362 bxibxj(2,2,ep)=fac2z*px2py1(ep)
363C
364 bxibyj(1,1,ep)=facz(ep)*(-px1px2(ep)+py1py2(ep))
365 bxibyj(1,2,ep)=facz(ep)*(-px1px1(ep)+py2py2(ep))
366 bxibyj(2,1,ep)=facz(ep)*(-px2px2(ep)+py1py1(ep))
367 bxibyj(2,2,ep)=bxibyj(1,1,ep)
368C
369 byibyj(1,1,ep)=-fac2z*px2py1(ep)
370 byibyj(1,2,ep)=facz(ep)*(-px1py1(ep)-px2py2(ep))
371 byibyj(2,1,ep)=byibyj(1,2,ep)
372 byibyj(2,2,ep)=-fac2z*px1py2(ep)
373C
374 byibxj(1,1,ep)=facz(ep)*(py1py2(ep)-px1px2(ep))
375 byibxj(1,2,ep)=facz(ep)*(py1py1(ep)-px2px2(ep))
376 byibxj(2,1,ep)=facz(ep)*(py2py2(ep)-px1px1(ep))
377 byibxj(2,2,ep)=byibxj(1,1,ep)
378C
379 ENDDO
380 DO ep=nf,jlt
381 k11(1,1,ep)=k11(1,1,ep)+dmf(1,1,ep)*bxibxj(1,1,ep)+
382 1 dmf(1,3,ep)*(bxibyj(1,1,ep)+byibxj(1,1,ep))+
383 2 dmf(3,3,ep)*byibyj(1,1,ep)
384 k11(1,2,ep)=k11(1,2,ep)+dmf(1,2,ep)*bxibyj(1,1,ep)+
385 1 dmf(1,3,ep)*bxibxj(1,1,ep)+
386 2 dmf(2,3,ep)*byibyj(1,1,ep)+
387 3 dmf(3,3,ep)*byibxj(1,1,ep)
388 k11(2,2,ep)=k11(2,2,ep)+dmf(2,2,ep)*byibyj(1,1,ep)+
389 1 dmf(2,3,ep)*(bxibyj(1,1,ep)+byibxj(1,1,ep))+
390 2 dmf(3,3,ep)*bxibxj(1,1,ep)
391C
392 k12(1,1,ep)=k12(1,1,ep)+dmf(1,1,ep)*bxibxj(1,2,ep)+
393 1 dmf(1,3,ep)*(bxibyj(1,2,ep)+byibxj(1,2,ep))+
394 2 dmf(3,3,ep)*byibyj(1,2,ep)
395 k12(1,2,ep)=k12(1,2,ep)+dmf(1,2,ep)*bxibyj(1,2,ep)+
396 1 dmf(1,3,ep)*bxibxj(1,2,ep)+
397 2 dmf(2,3,ep)*byibyj(1,2,ep)+
398 3 dmf(3,3,ep)*byibxj(1,2,ep)
399 k12(2,1,ep)=k12(2,1,ep)+dmf(1,2,ep)*byibxj(1,2,ep)+
400 1 dmf(1,3,ep)*bxibxj(1,2,ep)+
401 2 dmf(2,3,ep)*byibyj(1,2,ep)+
402 3 dmf(3,3,ep)*bxibyj(1,2,ep)
403 k12(2,2,ep)=k12(2,2,ep)+dmf(2,2,ep)*byibyj(1,2,ep)+
404 1 dmf(2,3,ep)*(bxibyj(1,2,ep)+byibxj(1,2,ep))+
405 2 dmf(3,3,ep)*bxibxj(1,2,ep)
406C
407 k22(1,1,ep)=k22(1,1,ep)+dmf(1,1,ep)*bxibxj(2,2,ep)+
408 1 dmf(1,3,ep)*(bxibyj(2,2,ep)+byibxj(2,2,ep))+
409 2 dmf(3,3,ep)*byibyj(2,2,ep)
410 k22(1,2,ep)=k22(1,2,ep)+dmf(1,2,ep)*bxibyj(2,2,ep)+
411 1 dmf(1,3,ep)*bxibxj(2,2,ep)+
412 2 dmf(2,3,ep)*byibyj(2,2,ep)+
413 3 dmf(3,3,ep)*byibxj(2,2,ep)
414 k22(2,2,ep)=k22(2,2,ep)+dmf(2,2,ep)*byibyj(2,2,ep)+
415 1 dmf(2,3,ep)*(bxibyj(2,2,ep)+byibxj(2,2,ep))+
416 2 dmf(3,3,ep)*bxibxj(2,2,ep)
417 ENDDO
418C
419 END IF !(IORTH==1)
420
421C
422 DO i=1,2
423 DO j=i,2
424 DO ep=jft,jlt
425 k13(i,j,ep)=-k11(i,j,ep)
426 k14(i,j,ep)=-k12(i,j,ep)
427 k23(i,j,ep)=-k12(j,i,ep)
428 k24(i,j,ep)=-k22(i,j,ep)
429 k33(i,j,ep)=k11(i,j,ep)
430 k34(i,j,ep)=k12(i,j,ep)
431 k44(i,j,ep)=k22(i,j,ep)
432 m13(i,j,ep)=-m11(i,j,ep)
433 m14(i,j,ep)=-m12(i,j,ep)
434 m23(i,j,ep)=-m12(j,i,ep)
435 m24(i,j,ep)=-m22(i,j,ep)
436 m33(i,j,ep)=m11(i,j,ep)
437 m34(i,j,ep)=m12(i,j,ep)
438 m44(i,j,ep)=m22(i,j,ep)
439 ENDDO
440 ENDDO
441 ENDDO
442C
443 DO ep=jft,jlt
444 k13(2,1,ep)=-k11(1,2,ep)
445 k14(2,1,ep)=-k12(2,1,ep)
446 k23(2,1,ep)=-k12(1,2,ep)
447 k24(2,1,ep)=-k22(1,2,ep)
448 k34(2,1,ep)=k12(2,1,ep)
449 m13(2,1,ep)=-m11(1,2,ep)
450 m14(2,1,ep)=-m12(2,1,ep)
451 m23(2,1,ep)=-m12(1,2,ep)
452 m24(2,1,ep)=-m22(1,2,ep)
453 m34(2,1,ep)=m12(2,1,ep)
454 ENDDO
455C---+---------+---------+warped elements-suite------------
456C--------px1'=py2;px2'=py1;py1'=-px2;py2'=-px1;-------------
457C--------px1'px1=px1py2; px1'px2=px2py2; px2'px2=px2py1;-------------
458C--------py1'py1=-px2py1; py1'py2=-px2py2; py2'py2=-px1py2;-------------
459C--------px1'py1=py1py2;px2'py1=py1py1;px2'py2=py1py2;px1'py2=py2py2;-------------
460 DO ep=nf,jlt
461 dg(ep)=gf(ep)*facz(ep)
462 dfz =df(1,1,ep)*facz(ep)
463 mf11(1,2,ep)=dfz*px1py2(ep)-dg(ep)*px2py1(ep)
464 mf12(1,2,ep)=(dfz-dg(ep))*px2py2(ep)
465 mf22(1,2,ep)=dfz*px2py1(ep)-dg(ep)*px1py2(ep)
466 mf23(1,2,ep)=-(dfz-dg(ep))*px1py1(ep)
467 ENDDO
468 DO ep=nf,jlt
469 dfz =df(1,2,ep)*facz(ep)
470 mf11(1,1,ep)=-dfz*py1py2(ep)+dg(ep)*px1px2(ep)
471 mf12(1,1,ep)=-dfz*py2py2(ep)+dg(ep)*px2px2(ep)
472 mf22(1,1,ep)=mf11(1,1,ep)
473 mf23(1,1,ep)=dfz*py1py1(ep)-dg(ep)*px1px1(ep)
474 mf11(2,2,ep)=-dfz*px1px2(ep)+dg(ep)*py1py2(ep)
475 mf12(2,2,ep)=-dfz*px2px2(ep)+dg(ep)*py2py2(ep)
476 mf22(2,2,ep)= mf11(2,2,ep)
477 mf23(2,2,ep)=dfz*px1px1(ep)-dg(ep)*py1py1(ep)
478 ENDDO
479 DO ep=nf,jlt
480 dfz =df(2,2,ep)*facz(ep)
481 mf11(2,1,ep)=dfz*px2py1(ep)-dg(ep)*px1py2(ep)
482 mf12(2,1,ep)=dfz*px2py2(ep)-dg(ep)*px2py2(ep)
483 mf22(2,1,ep)=dfz*px1py2(ep)-dg(ep)*px2py1(ep)
484 mf23(2,1,ep)=-(dfz-dg(ep))*px1py1(ep)
485 ENDDO
486 IF (iorth>0) THEN
487 DO m=nf,jlt
488 df5(m)=facz(m)*df5(m)
489 df6(m)=facz(m)*df6(m)
490 ENDDO
491C
492 DO ep=nf,jlt
493 dfz=df5(ep)*px1py2(ep)-df6(ep)*px2py1(ep)
494 mf11(1,1,ep)=mf11(1,1,ep)-dfz
495 mf11(1,2,ep)=mf11(1,2,ep)+df5(ep)*(py1py2(ep)-px1px2(ep))
496 mf11(2,1,ep)=mf11(2,1,ep)+df6(ep)*(px1px2(ep)-py1py2(ep))
497 mf11(2,2,ep)=mf11(2,2,ep)+dfz
498 ENDDO
499 DO ep=nf,jlt
500 dfz=(df5(ep)-df6(ep))*px2py2(ep)
501 mf12(1,1,ep)=mf12(1,1,ep)-dfz
502 mf12(1,2,ep)=mf12(1,2,ep)+df5(ep)*(py2py2(ep)-px2px2(ep))
503 mf12(2,1,ep)=mf12(2,1,ep)+df6(ep)*(px2px2(ep)-py2py2(ep))
504 mf12(2,2,ep)=mf12(2,2,ep)+dfz
505 ENDDO
506 DO ep=nf,jlt
507 dfz=df5(ep)*px2py1(ep)-df6(ep)*px1py2(ep)
508 mf22(1,1,ep)=mf22(1,1,ep)-dfz
509 mf22(1,2,ep)=mf22(1,2,ep)+df5(ep)*(py1py2(ep)-px1px2(ep))
510 mf22(2,1,ep)=mf22(2,1,ep)+df6(ep)*(px1px2(ep)-py1py2(ep))
511 mf22(2,2,ep)=mf22(2,2,ep)+dfz
512 ENDDO
513 DO ep=nf,jlt
514 dfz=(df6(ep)-df5(ep))*px1py1(ep)
515 mf23(1,1,ep)=mf23(1,1,ep)-dfz
516 mf23(1,2,ep)=mf23(1,2,ep)-df5(ep)*(py1py1(ep)-px1px1(ep))
517 mf23(2,1,ep)=mf23(2,1,ep)-df6(ep)*(px1px1(ep)-py1py1(ep))
518 mf23(2,2,ep)=mf23(2,2,ep)+dfz
519 ENDDO
520C ----add mem-bending coupling terms of orthotrope--warped(same than coplanar---
521 DO ep=nf,jlt
522 mf11(1,1,ep)=mf11(1,1,ep)
523 1 -dmf(1,2,ep)*px1py1(ep)-dmf(1,3,ep)*px1px1(ep)
524 2 -dmf(2,3,ep)*py1py1(ep)-dmf(3,3,ep)*px1py1(ep)
525 mf11(1,2,ep)=mf11(1,2,ep)
526 1 +dmf(1,1,ep)*px1px1(ep)+dmf(1,3,ep)*px1py1(ep)
527 2 +dmf(1,3,ep)*px1py1(ep)+dmf(3,3,ep)*py1py1(ep)
528 mf11(2,1,ep)=mf11(2,1,ep)
529 1 -dmf(2,2,ep)*py1py1(ep)-dmf(2,3,ep)*px1py1(ep)
530 2 -dmf(2,3,ep)*px1py1(ep)-dmf(3,3,ep)*px1px1(ep)
531 mf11(2,2,ep)=mf11(2,2,ep)
532 1 +dmf(1,2,ep)*px1py1(ep)+dmf(2,3,ep)*py1py1(ep)
533 2 +dmf(1,3,ep)*px1px1(ep)+dmf(3,3,ep)*px1py1(ep)
534 ENDDO
535 DO ep=nf,jlt
536 mf12(1,1,ep)=mf12(1,1,ep)
537 1 -dmf(1,2,ep)*px1py2(ep)-dmf(1,3,ep)*px1px2(ep)
538 2 -dmf(2,3,ep)*py1py2(ep)-dmf(3,3,ep)*px2py1(ep)
539 mf12(1,2,ep)=mf12(1,2,ep)
540 1 +dmf(1,1,ep)*px1px2(ep)+dmf(1,3,ep)*px1py2(ep)
541 2 +dmf(1,3,ep)*px2py1(ep)+dmf(3,3,ep)*py1py2(ep)
542 mf12(2,1,ep)=mf12(2,1,ep)
543 1 -dmf(2,2,ep)*py1py2(ep)-dmf(2,3,ep)*px2py1(ep)
544 2 -dmf(2,3,ep)*px1py2(ep)-dmf(3,3,ep)*px1px2(ep)
545 mf12(2,2,ep)=mf12(2,2,ep)
546 1 +dmf(1,2,ep)*px2py1(ep)+dmf(2,3,ep)*py1py2(ep)
547 2 +dmf(1,3,ep)*px1px2(ep)+dmf(3,3,ep)*px1py2(ep)
548 ENDDO
549 DO ep=nf,jlt
550 mf22(1,1,ep)=mf22(1,1,ep)
551 1 -dmf(1,2,ep)*px2py2(ep)-dmf(1,3,ep)*px2px2(ep)
552 2 -dmf(2,3,ep)*py2py2(ep)-dmf(3,3,ep)*px2py2(ep)
553 mf22(1,2,ep)=mf22(1,2,ep)
554 1 +dmf(1,1,ep)*px2px2(ep)+dmf(1,3,ep)*px2py2(ep)
555 2 +dmf(1,3,ep)*px2py2(ep)+dmf(3,3,ep)*py2py2(ep)
556 mf22(2,1,ep)=mf22(2,1,ep)
557 1 -dmf(2,2,ep)*py2py2(ep)-dmf(2,3,ep)*px2py2(ep)
558 2 -dmf(2,3,ep)*px2py2(ep)-dmf(3,3,ep)*px2px2(ep)
559 mf22(2,2,ep)=mf22(2,2,ep)
560 1 +dmf(1,2,ep)*px2py2(ep)+dmf(2,3,ep)*py2py2(ep)
561 2 +dmf(1,3,ep)*px2px2(ep)+dmf(3,3,ep)*px2py2(ep)
562 ENDDO
563 DO ep=nf,jlt
564 mf23(1,1,ep)=mf23(1,1,ep)
565 1 +dmf(1,2,ep)*px2py1(ep)+dmf(1,3,ep)*px1px2(ep)
566 2 +dmf(2,3,ep)*py1py2(ep)+dmf(3,3,ep)*px1py2(ep)
567 mf23(1,2,ep)=mf23(1,2,ep)
568 1 -dmf(1,1,ep)*px1px2(ep)-dmf(1,3,ep)*px2py1(ep)
569 2 -dmf(1,3,ep)*px1py2(ep)-dmf(3,3,ep)*py1py2(ep)
570 mf23(2,1,ep)=mf23(2,1,ep)
571 1 +dmf(2,2,ep)*py1py2(ep)+dmf(2,3,ep)*px1py2(ep)
572 2 +dmf(2,3,ep)*px2py1(ep)+dmf(3,3,ep)*px1px2(ep)
573 mf23(2,2,ep)=mf23(2,2,ep)
574 1 -dmf(1,2,ep)*px1py2(ep)-dmf(2,3,ep)*py1py2(ep)
575 2 -dmf(1,3,ep)*px1px2(ep)-dmf(3,3,ep)*px2py1(ep)
576 ENDDO
577 END IF !(IORTH==1)
578C
579 DO i=1,2
580 DO j=1,2
581 DO ep=nf,jlt
582 mf13(i,j,ep)=-mf11(i,j,ep)
583 mf14(i,j,ep)=-mf12(i,j,ep)
584 mf24(i,j,ep)=-mf22(i,j,ep)
585 mf33(i,j,ep)=mf11(i,j,ep)
586 mf34(i,j,ep)=mf12(i,j,ep)
587 mf44(i,j,ep)=mf22(i,j,ep)
588 ENDDO
589 ENDDO
590 ENDDO
591C
592 DO i=1,2
593 DO j=1,2
594 DO ep=nf,jlt
595 fm12(i,j,ep)=-mf23(j,i,ep)
596 fm13(i,j,ep)=-mf11(j,i,ep)
597 fm14(i,j,ep)= mf23(j,i,ep)
598 fm23(i,j,ep)=-mf12(j,i,ep)
599 fm24(i,j,ep)=-mf22(j,i,ep)
600 fm34(i,j,ep)=-mf23(j,i,ep)
601 ENDDO
602 ENDDO
603 ENDDO
604C------ Mzz ----------
605 DO ep=jft,jlt
606 m11(3,3,ep)=dhz(ep)*(px1px1(ep)+py1py1(ep))
607 m12(3,3,ep)=dhz(ep)*(px1px2(ep)+py1py2(ep))
608 m13(3,3,ep)=-m11(3,3,ep)
609 m14(3,3,ep)=-m12(3,3,ep)
610 m22(3,3,ep)=dhz(ep)*(px2px2(ep)+py2py2(ep))
611 m23(3,3,ep)=-m12(3,3,ep)
612 m24(3,3,ep)=-m22(3,3,ep)
613 m33(3,3,ep)=m11(3,3,ep)
614 m34(3,3,ep)=m12(3,3,ep)
615 m44(3,3,ep)=m22(3,3,ep)
616 ENDDO
617C------ renforce positive ----------
618 IF (neig==0) THEN
619 DO ep=jft,jlt
620 c1 = em8*m11(3,3,ep)
621 c2 = em8*m22(3,3,ep)
622 m11(3,3,ep)=m11(3,3,ep)+c1
623 m22(3,3,ep)=m22(3,3,ep)+c2
624 m33(3,3,ep)=m11(3,3,ep)
625 m44(3,3,ep)=m22(3,3,ep)
626 ENDDO
627 ENDIF
628C
629 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ czlkecr3()

subroutine czlkecr3 ( integer jft,
integer jlt,
vol,
thk0,
thk2,
hm,
hf,
hz,
a_i,
z1,
px1,
px2,
py1,
py2,
k11,
k12,
k13,
k14,
k22,
k23,
k24,
k33,
k34,
k44,
m11,
m12,
m13,
m14,
m22,
m23,
m24,
m33,
m34,
m44,
mf11,
mf12,
mf13,
mf14,
mf22,
mf23,
mf24,
mf33,
mf34,
mf44,
fm12,
fm13,
fm14,
fm23,
fm24,
fm34,
integer iorth,
hmor,
hfor,
integer, dimension(*) iplat,
dhz,
prx,
pry,
prxy,
prz,
hmfor,
integer nplat )

Definition at line 636 of file czlkec3.F.

645C--------------------------------------------------------------------------------------------------
646C-----------------------------------------------
647C I M P L I C I T T Y P E S
648C-----------------------------------------------
649#include "implicit_f.inc"
650#include "mvsiz_p.inc"
651C-----------------------------------------------
652C D U M M Y A R G U M E N T S
653C-----------------------------------------------
654C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
655 INTEGER JFT,JLT,IPLAT(*),IORTH,NPLAT
656 my_real
657 . vol(*),thk0(*),thk2(*),a_i(*),z1(*),
658 . hm(mvsiz,4),hf(mvsiz,4),hz(*),
659 . px1(*) ,px2(*) ,py1(*) ,py2(*) ,
660 . k11(3,3,*),k12(3,3,*),k13(3,3,*),k14(3,3,*),
661 . k22(3,3,*),k23(3,3,*),k24(3,3,*),k33(3,3,*),
662 . m11(3,3,*),m12(3,3,*),m13(3,3,*),m14(3,3,*),
663 . m22(3,3,*),m23(3,3,*),m24(3,3,*),m33(3,3,*),
664 . mf11(3,3,*),mf12(3,3,*),mf13(3,3,*),mf14(3,3,*),
665 . mf22(3,3,*),mf23(3,3,*),mf24(3,3,*),mf33(3,3,*),
666 . fm12(3,3,*),fm13(3,3,*),fm14(3,3,*),
667 . fm23(3,3,*),fm24(3,3,*),fm34(3,3,*),
668 . k34(3,3,*),k44(3,3,*),m34(3,3,*),m44(3,3,*),
669 . mf34(3,3,*),mf44(3,3,*),dhz(*),hmor(mvsiz,2),hfor(mvsiz,2)
670 my_real
671 . prx(4,*),pry(4,*),prxy(4,*),prz(4,*) ,hmfor(mvsiz,6)
672C---------------|[KIJ][MFIJ]|----
673C-----KE(6x6)= | |
674C---------------|[FMIJ]{MIJ]|----
675C-----------------------------------------------
676C L O C A L V A R I A B L E S
677C-----------------------------------------------
678 INTEGER EP,I,J,NF,M
679 my_real
680 . dm(3,3,mvsiz),dprz(4,mvsiz),
681 . c1,c2,cbrx(4,mvsiz),cbry(4,mvsiz),cbrz(4,mvsiz),
682 . kz11(2,2,mvsiz),kz12(2,2,mvsiz),kz22(2,2,mvsiz),
683 . dmf(3,3,mvsiz),pcx1(mvsiz),pcx2(mvsiz),pcy1(mvsiz),
684 . pcy2(mvsiz),facz(mvsiz)
685C-----------Attention Matrice sym Kii ne calcul que la moitie---------72
686#include "vectorize.inc"
687C--------Constante parts---------------
688 DO m=jft,jlt
689 ep=iplat(m)
690 c2=vol(ep)
691 dm(1,1,m)=hm(ep,1)*c2
692 dm(2,2,m)=hm(ep,2)*c2
693 dm(1,2,m)=hm(ep,3)*c2
694 dm(3,3,m)=hm(ep,4)*c2
695 dhz(m)=hz(ep)*fourth*c2
696 dm(1,3,m)=zero
697 dm(2,3,m)=zero
698 ENDDO
699 IF (iorth==1) THEN
700 DO m=jft,jlt
701 ep=iplat(m)
702 c2=vol(ep)
703 dm(1,3,m)=hmor(ep,1)*c2
704 dm(2,3,m)=hmor(ep,2)*c2
705 ENDDO
706 ENDIF
707C
708 DO m=jft,jlt
709 dm(2,1,m)=dm(1,2,m)
710 dm(3,1,m)=dm(1,3,m)
711 dm(3,2,m)=dm(2,3,m)
712 ENDDO
713C-------0.5*[-By Bx 0 0 0 BRZ]^tKG[-By Bx 0 0 0 BRZ]*0.5
714C------ BRZ = (-BRXY+BRYX)-0.5
715 DO ep=jft,jlt
716 kz11(1,1,ep)=py1(ep)*dhz(ep)*py1(ep)
717 kz11(1,2,ep)=-py1(ep)*dhz(ep)*px1(ep)
718 kz11(2,2,ep)=px1(ep)*dhz(ep)*px1(ep)
719 kz12(1,1,ep)=py1(ep)*dhz(ep)*py2(ep)
720 kz12(1,2,ep)=-py1(ep)*dhz(ep)*px2(ep)
721 kz12(2,1,ep)=-px1(ep)*dhz(ep)*py2(ep)
722 kz12(2,2,ep)=px1(ep)*dhz(ep)*px2(ep)
723 kz22(1,1,ep)=py2(ep)*dhz(ep)*py2(ep)
724 kz22(1,2,ep)=-py2(ep)*dhz(ep)*px2(ep)
725 kz22(2,2,ep)=px2(ep)*dhz(ep)*px2(ep)
726 ENDDO
727C
728 DO i=1,2
729 DO j=i,2
730 DO ep=jft,jlt
731 k11(i,j,ep)=k11(i,j,ep)+kz11(i,j,ep)
732 k12(i,j,ep)=k12(i,j,ep)+kz12(i,j,ep)
733 k22(i,j,ep)=k22(i,j,ep)+kz22(i,j,ep)
734 k13(i,j,ep)=k13(i,j,ep)-kz11(i,j,ep)
735 k14(i,j,ep)=k14(i,j,ep)-kz12(i,j,ep)
736 k23(i,j,ep)=k23(i,j,ep)-kz12(j,i,ep)
737 k24(i,j,ep)=k24(i,j,ep)-kz22(i,j,ep)
738 k33(i,j,ep)=k33(i,j,ep)+kz11(i,j,ep)
739 k34(i,j,ep)=k34(i,j,ep)+kz12(i,j,ep)
740 k44(i,j,ep)=k44(i,j,ep)+kz22(i,j,ep)
741 ENDDO
742 ENDDO
743 ENDDO
744C
745 DO ep=jft,jlt
746 k12(2,1,ep)=k12(2,1,ep)+kz12(2,1,ep)
747 k13(2,1,ep)=k13(2,1,ep)-kz11(1,2,ep)
748 k14(2,1,ep)=k14(2,1,ep)-kz12(2,1,ep)
749 k23(2,1,ep)=k23(2,1,ep)-kz12(1,2,ep)
750 k24(2,1,ep)=k24(2,1,ep)-kz22(1,2,ep)
751 k34(2,1,ep)=k34(2,1,ep)+kz12(2,1,ep)
752 ENDDO
753C
754 DO j=1,4
755 DO m=jft,jlt
756 dprz(j,m)=prz(j,m)*dhz(m)
757 ENDDO
758 ENDDO
759C
760 DO m=jft,jlt
761 c1 = -py2(m)*dprz(1,m)
762 c2 = px2(m)*dprz(1,m)
763 fm12(3,1,m)=fm12(3,1,m)+ c1
764 fm14(3,1,m)=fm14(3,1,m)- c1
765 fm12(3,2,m)=fm12(3,2,m)+ c2
766 fm14(3,2,m)=fm14(3,2,m)- c2
767 c1 = -py1(m)*dprz(1,m)
768 c2 = px1(m)*dprz(1,m)
769 mf11(1,3,m)=mf11(1,3,m)+ c1
770 fm13(3,1,m)=fm13(3,1,m)- c1
771 mf11(2,3,m)=mf11(2,3,m)+ c2
772 fm13(3,2,m)=fm13(3,2,m)- c2
773 c1 = -py1(m)*dprz(2,m)
774 c2 = px1(m)*dprz(2,m)
775 mf12(1,3,m)=mf12(1,3,m)+ c1
776 fm23(3,1,m)=fm23(3,1,m)- c1
777 mf12(2,3,m)=mf12(2,3,m)+ c2
778 fm23(3,2,m)=fm23(3,2,m)- c2
779 c1 = -py1(m)*dprz(3,m)
780 c2 = px1(m)*dprz(3,m)
781 mf13(1,3,m)=mf13(1,3,m)+ c1
782 mf33(1,3,m)=mf33(1,3,m)- c1
783 mf13(2,3,m)=mf13(2,3,m)+ c2
784 mf33(2,3,m)=mf33(2,3,m)- c2
785 c1 = -py1(m)*dprz(4,m)
786 c2 = px1(m)*dprz(4,m)
787 mf14(1,3,m)=mf14(1,3,m)+ c1
788 mf34(1,3,m)=mf34(1,3,m)- c1
789 mf14(2,3,m)=mf14(2,3,m)+ c2
790 mf34(2,3,m)=mf34(2,3,m)- c2
791 c1 = -py2(m)*dprz(2,m)
792 c2 = px2(m)*dprz(2,m)
793 mf22(1,3,m)=mf22(1,3,m)+ c1
794 fm24(3,1,m)=fm24(3,1,m)- c1
795 mf22(2,3,m)=mf22(2,3,m)+ c2
796 fm24(3,2,m)=fm24(3,2,m)- c2
797 c1 = -py2(m)*dprz(3,m)
798 c2 = px2(m)*dprz(3,m)
799 mf23(1,3,m)=mf23(1,3,m)+ c1
800 fm34(3,1,m)=fm34(3,1,m)- c1
801 mf23(2,3,m)=mf23(2,3,m)+ c2
802 fm34(3,2,m)=fm34(3,2,m)- c2
803 c1 = -py2(m)*dprz(4,m)
804 c2 = px2(m)*dprz(4,m)
805 mf24(1,3,m)=mf24(1,3,m)+ c1
806 mf44(1,3,m)=mf44(1,3,m)- c1
807 mf24(2,3,m)=mf24(2,3,m)+ c2
808 mf44(2,3,m)=mf44(2,3,m)- c2
809 ENDDO
810C------ Mzz ----------
811 DO ep=jft,jlt
812 m11(3,3,ep)=prz(1,ep)*dprz(1,ep)
813 m12(3,3,ep)=prz(1,ep)*dprz(2,ep)
814 m13(3,3,ep)=prz(1,ep)*dprz(3,ep)
815 m14(3,3,ep)=prz(1,ep)*dprz(4,ep)
816 m22(3,3,ep)=prz(2,ep)*dprz(2,ep)
817 m23(3,3,ep)=prz(2,ep)*dprz(3,ep)
818 m24(3,3,ep)=prz(2,ep)*dprz(4,ep)
819 m33(3,3,ep)=prz(3,ep)*dprz(3,ep)
820 m34(3,3,ep)=prz(3,ep)*dprz(4,ep)
821 m44(3,3,ep)=prz(4,ep)*dprz(4,ep)
822 ENDDO
823C-------[MFIJ]=[Bm]^t[C][BRm]; [MIJ]=[BRm]^t[C][BRm];----
824C---------------|0 0 BRX |----
825C-----BR(3x3)= |0 0 BRY |
826C---------------|0 0 BRXY+BRYX |----
827 DO j=1,4
828 DO m=jft,jlt
829 cbrx(j,m) =dm(1,1,m)*prx(j,m)+dm(1,2,m)*pry(j,m)+
830 . dm(1,3,m)*prxy(j,m)
831 cbry(j,m) =dm(2,1,m)*prx(j,m)+dm(2,2,m)*pry(j,m)+
832 . dm(2,3,m)*prxy(j,m)
833 cbrz(j,m) =dm(3,1,m)*prx(j,m)+dm(3,2,m)*pry(j,m)+
834 . dm(3,3,m)*prxy(j,m)
835 ENDDO
836 ENDDO
837C
838 DO m=jft,jlt
839 c1 = px2(m)*cbrx(1,m)+py2(m)*cbrz(1,m)
840 c2 = py2(m)*cbry(1,m)+px2(m)*cbrz(1,m)
841 fm12(3,1,m)=fm12(3,1,m)+ c1
842 fm14(3,1,m)=fm14(3,1,m)- c1
843 fm12(3,2,m)=fm12(3,2,m)+ c2
844 fm14(3,2,m)=fm14(3,2,m)- c2
845C
846 c1 = px1(m)*cbrx(1,m)+py1(m)*cbrz(1,m)
847 c2 = py1(m)*cbry(1,m)+px1(m)*cbrz(1,m)
848 mf11(1,3,m)=mf11(1,3,m)+ c1
849 fm13(3,1,m)=fm13(3,1,m)- c1
850 mf11(2,3,m)=mf11(2,3,m)+ c2
851 fm13(3,2,m)=fm13(3,2,m)- c2
852 c1 = px1(m)*cbrx(2,m)+py1(m)*cbrz(2,m)
853 c2 = py1(m)*cbry(2,m)+px1(m)*cbrz(2,m)
854 mf12(1,3,m)=mf12(1,3,m)+ c1
855 fm23(3,1,m)=fm23(3,1,m)- c1
856 mf12(2,3,m)=mf12(2,3,m)+ c2
857 fm23(3,2,m)=fm23(3,2,m)- c2
858 c1 = px1(m)*cbrx(3,m)+py1(m)*cbrz(3,m)
859 c2 = py1(m)*cbry(3,m)+px1(m)*cbrz(3,m)
860 mf13(1,3,m)=mf13(1,3,m)+ c1
861 mf33(1,3,m)=mf33(1,3,m)- c1
862 mf13(2,3,m)=mf13(2,3,m)+ c2
863 mf33(2,3,m)=mf33(2,3,m)- c2
864 c1 = px1(m)*cbrx(4,m)+py1(m)*cbrz(4,m)
865 c2 = py1(m)*cbry(4,m)+px1(m)*cbrz(4,m)
866 mf14(1,3,m)=mf14(1,3,m)+ c1
867 mf34(1,3,m)=mf34(1,3,m)- c1
868 mf14(2,3,m)=mf14(2,3,m)+ c2
869 mf34(2,3,m)=mf34(2,3,m)- c2
870 c1 = px2(m)*cbrx(2,m)+py2(m)*cbrz(2,m)
871 c2 = py2(m)*cbry(2,m)+px2(m)*cbrz(2,m)
872 mf22(1,3,m)=mf22(1,3,m)+ c1
873 fm24(3,1,m)=fm24(3,1,m)- c1
874 mf22(2,3,m)=mf22(2,3,m)+ c2
875 fm24(3,2,m)=fm24(3,2,m)- c2
876 c1 = px2(m)*cbrx(3,m)+py2(m)*cbrz(3,m)
877 c2 = py2(m)*cbry(3,m)+px2(m)*cbrz(3,m)
878 mf23(1,3,m)=mf23(1,3,m)+ c1
879 fm34(3,1,m)=fm34(3,1,m)- c1
880 mf23(2,3,m)=mf23(2,3,m)+ c2
881 fm34(3,2,m)=fm34(3,2,m)- c2
882 c1 = px2(m)*cbrx(4,m)+py2(m)*cbrz(4,m)
883 c2 = py2(m)*cbry(4,m)+px2(m)*cbrz(4,m)
884 mf24(1,3,m)=mf24(1,3,m)+ c1
885 mf44(1,3,m)=mf44(1,3,m)- c1
886 mf24(2,3,m)=mf24(2,3,m)+ c2
887 mf44(2,3,m)=mf44(2,3,m)- c2
888 ENDDO
889C------ Mzz ----------
890 DO m=jft,jlt
891 m11(3,3,m)=m11(3,3,m)+prx(1,m)*cbrx(1,m)+
892 . pry(1,m)*cbry(1,m)+prxy(1,m)*cbrz(1,m)
893 m12(3,3,m)=m12(3,3,m)+prx(1,m)*cbrx(2,m)+
894 . pry(1,m)*cbry(2,m)+prxy(1,m)*cbrz(2,m)
895 m13(3,3,m)=m13(3,3,m)+prx(1,m)*cbrx(3,m)+
896 . pry(1,m)*cbry(3,m)+prxy(1,m)*cbrz(3,m)
897 m14(3,3,m)=m14(3,3,m)+prx(1,m)*cbrx(4,m)+
898 . pry(1,m)*cbry(4,m)+prxy(1,m)*cbrz(4,m)
899 m22(3,3,m)=m22(3,3,m)+prx(2,m)*cbrx(2,m)+
900 . pry(2,m)*cbry(2,m)+prxy(2,m)*cbrz(2,m)
901 m23(3,3,m)=m23(3,3,m)+prx(2,m)*cbrx(3,m)+
902 . pry(2,m)*cbry(3,m)+prxy(2,m)*cbrz(3,m)
903 m24(3,3,m)=m24(3,3,m)+prx(2,m)*cbrx(4,m)+
904 . pry(2,m)*cbry(4,m)+prxy(2,m)*cbrz(4,m)
905 m33(3,3,m)=m33(3,3,m)+prx(3,m)*cbrx(3,m)+
906 . pry(3,m)*cbry(3,m)+prxy(3,m)*cbrz(3,m)
907 m34(3,3,m)=m34(3,3,m)+prx(3,m)*cbrx(4,m)+
908 . pry(3,m)*cbry(4,m)+prxy(3,m)*cbrz(4,m)
909 m44(3,3,m)=m44(3,3,m)+prx(4,m)*cbrx(4,m)+
910 . pry(4,m)*cbry(4,m)+prxy(4,m)*cbrz(4,m)
911 ENDDO
912 IF (iorth>0) THEN
913C ----add mem-bending coupling terms of orthotrope--coplanar first---
914 nf=nplat+1
915#include "vectorize.inc"
916 DO m=jft,jlt
917 ep=iplat(m)
918 c2=vol(ep)*thk0(ep)
919 dmf(1,1,m)=hmfor(ep,1)*c2
920 dmf(2,2,m)=hmfor(ep,2)*c2
921 dmf(1,2,m)=hmfor(ep,3)*c2
922 dmf(1,3,m)=hmfor(ep,5)*c2
923 dmf(2,3,m)=hmfor(ep,6)*c2
924 dmf(2,1,m)=dmf(1,2,m)
925 dmf(3,1,m)=dmf(1,3,m)
926 dmf(3,2,m)=dmf(2,3,m)
927 dmf(3,3,m)=hmfor(ep,4)*c2
928 ENDDO
929C-------[C][BRm];----
930 DO j=1,4
931 DO ep=jft,jlt
932 cbrx(j,ep) =dmf(1,1,ep)*prx(j,ep)+dmf(1,2,ep)*pry(j,ep)+
933 . dmf(1,3,ep)*prxy(j,ep)
934 cbry(j,ep) =dmf(2,1,ep)*prx(j,ep)+dmf(2,2,ep)*pry(j,ep)+
935 . dmf(2,3,ep)*prxy(j,ep)
936 cbrz(j,ep) =dmf(3,1,ep)*prx(j,ep)+dmf(3,2,ep)*pry(j,ep)+
937 . dmf(3,3,ep)*prxy(j,ep)
938 ENDDO
939 ENDDO
940C
941 DO m=nf,jlt
942 ep=iplat(m)
943 facz(m) = four*z1(ep)*a_i(ep)
944 ENDDO
945C--------BMF : px1'=py2;px2'=py1;py1'=-px2;py2'=-px1;-------------
946 DO ep=nf,jlt
947 pcx1(ep)=facz(ep)*py2(ep)
948 pcx2(ep)=facz(ep)*py1(ep)
949 pcy1(ep)=-facz(ep)*px2(ep)
950 pcy2(ep)=-facz(ep)*px1(ep)
951 ENDDO
952C ----add Rz-bending coupling terms of orthotrope-----[Bf]^t[C][BRm]+[Bmf]^t[C][BRm]
953C MIJ(1,3)=-BY(I)*CBR2(J)-BX(I)*CBR3(J)
954C MIJ(2,3)=BX(I)*CBR1(J)+BY(I)*CBR3(J)
955C MIJ(3,1)=-BY(J)*CBR2(I)-BX(J)*CBR3(I)
956C MIJ(3,2)=BX(J)*CBR1(I)+BY(J)*CBR3(I)
957C
958 DO ep=jft,jlt
959 c1=-py1(ep)*cbry(1,ep)-px1(ep)*cbrz(1,ep)
960 c2= px1(ep)*cbrx(1,ep)+py1(ep)*cbrz(1,ep)
961 m11(1,3,ep)=m11(1,3,ep)+c1
962 m11(2,3,ep)=m11(2,3,ep)+c2
963 m13(3,1,ep)=m13(3,1,ep)-c1
964 m13(3,2,ep)=m13(3,2,ep)-c2
965 c1=-py1(ep)*cbry(2,ep)-px1(ep)*cbrz(2,ep)
966 c2= px1(ep)*cbrx(2,ep)+py1(ep)*cbrz(2,ep)
967 m12(1,3,ep)=m12(1,3,ep)+c1
968 m12(2,3,ep)=m12(2,3,ep)+c2
969 m23(3,1,ep)=m23(3,1,ep)-c1
970 m23(3,2,ep)=m23(3,2,ep)-c2
971 c1=-py1(ep)*cbry(4,ep)-px1(ep)*cbrz(4,ep)
972 c2= px1(ep)*cbrx(4,ep)+py1(ep)*cbrz(4,ep)
973 m14(1,3,ep)=m14(1,3,ep)+c1
974 m14(2,3,ep)=m14(2,3,ep)+c2
975 m34(1,3,ep)=m34(1,3,ep)-c1
976 m34(2,3,ep)=m34(2,3,ep)-c2
977 c1=-py2(ep)*cbry(2,ep)-px2(ep)*cbrz(2,ep)
978 c2= px2(ep)*cbrx(2,ep)+py2(ep)*cbrz(2,ep)
979 m22(1,3,ep)=m22(1,3,ep)+c1
980 m22(2,3,ep)=m22(2,3,ep)+c2
981 m24(3,1,ep)=m24(3,1,ep)-c1
982 m24(3,2,ep)=m24(3,2,ep)-c2
983 c1=-py2(ep)*cbry(3,ep)-px2(ep)*cbrz(3,ep)
984 c2= px2(ep)*cbrx(3,ep)+py2(ep)*cbrz(3,ep)
985 m23(1,3,ep)=m23(1,3,ep)+c1
986 m23(2,3,ep)=m23(2,3,ep)+c2
987 m34(3,1,ep)=m34(3,1,ep)-c1
988 m34(3,2,ep)=m34(3,2,ep)-c2
989 c1= py1(ep)*cbry(3,ep)+px1(ep)*cbrz(3,ep)
990 c2=-px1(ep)*cbrx(3,ep)-py1(ep)*cbrz(3,ep)
991 m33(1,3,ep)=m33(1,3,ep)+c1
992 m33(2,3,ep)=m33(2,3,ep)+c2
993 m13(1,3,ep)=m13(1,3,ep)-c1
994 m13(2,3,ep)=m13(2,3,ep)-c2
995 c1= py2(ep)*cbry(4,ep)+px2(ep)*cbrz(4,ep)
996 c2=-px2(ep)*cbrx(4,ep)-py2(ep)*cbrz(4,ep)
997 m44(1,3,ep)=m44(1,3,ep)+c1
998 m44(2,3,ep)=m44(2,3,ep)+c2
999 m24(1,3,ep)=m24(1,3,ep)-c1
1000 m24(2,3,ep)=m24(2,3,ep)-c2
1001 c1=-py2(ep)*cbry(1,ep)-px2(ep)*cbrz(1,ep)
1002 c2= px2(ep)*cbrx(1,ep)+py2(ep)*cbrz(1,ep)
1003 m12(3,1,ep)=m12(3,1,ep)+c1
1004 m12(3,2,ep)=m12(3,2,ep)+c2
1005 m14(3,1,ep)=m14(3,1,ep)-c1
1006 m14(3,2,ep)=m14(3,2,ep)-c2
1007 ENDDO
1008C---+----1----+----2----+warped elements----+----5----+----6----+----7----+----8
1009C MFIJ(1,3)=BCX(I)*CBR1(J)+BCY(I)*CBR3(J)
1010C MFIJ(2,3)=BCY(I)*CBR2(J)+BCX(I)*CBR3(J)
1011C FMIJ(3,1)=BCX(J)*CBR1(I)+BCY(J)*CBR3(I)
1012C FMIJ(3,2)=BCY(J)*CBR2(I)+BCX(J)*CBR3(I)
1013C
1014 DO ep=nf,jlt
1015 c1= pcx1(ep)*cbrx(1,ep)+pcy1(ep)*cbrz(1,ep)
1016 c2= pcy1(ep)*cbry(1,ep)+pcx1(ep)*cbrz(1,ep)
1017 mf11(1,3,ep)=mf11(1,3,ep)+c1
1018 mf11(2,3,ep)=mf11(2,3,ep)+c2
1019 fm13(3,1,ep)=fm13(3,1,ep)-c1
1020 fm13(3,2,ep)=fm13(3,2,ep)-c2
1021 c1= pcx1(ep)*cbrx(2,ep)+pcy1(ep)*cbrz(2,ep)
1022 c2= pcy1(ep)*cbry(2,ep)+pcx1(ep)*cbrz(2,ep)
1023 mf12(1,3,ep)=mf12(1,3,ep)+c1
1024 mf12(2,3,ep)=mf12(2,3,ep)+c2
1025 fm23(3,1,ep)=fm23(3,1,ep)-c1
1026 fm23(3,2,ep)=fm23(3,2,ep)-c2
1027 c1= pcx1(ep)*cbrx(4,ep)+pcy1(ep)*cbrz(4,ep)
1028 c2= pcy1(ep)*cbry(4,ep)+pcx1(ep)*cbrz(4,ep)
1029 mf14(1,3,ep)=mf14(1,3,ep)+c1
1030 mf14(2,3,ep)=mf14(2,3,ep)+c2
1031 mf34(1,3,ep)=mf34(1,3,ep)-c1
1032 mf34(2,3,ep)=mf34(2,3,ep)-c2
1033 c1= pcx2(ep)*cbrx(2,ep)+pcy2(ep)*cbrz(2,ep)
1034 c2= pcy2(ep)*cbry(2,ep)+pcx2(ep)*cbrz(2,ep)
1035 mf22(1,3,ep)=mf22(1,3,ep)+c1
1036 mf22(2,3,ep)=mf22(2,3,ep)+c2
1037 fm24(3,1,ep)=fm24(3,1,ep)-c1
1038 fm24(3,2,ep)=fm24(3,2,ep)-c2
1039 c1= pcx2(ep)*cbrx(3,ep)+pcy2(ep)*cbrz(3,ep)
1040 c2= pcy2(ep)*cbry(3,ep)+pcx2(ep)*cbrz(3,ep)
1041 mf23(1,3,ep)=mf23(1,3,ep)+c1
1042 mf23(2,3,ep)=mf23(2,3,ep)+c2
1043 fm34(3,1,ep)=fm34(3,1,ep)-c1
1044 fm34(3,2,ep)=fm34(3,2,ep)-c2
1045 c1=-pcx1(ep)*cbrx(3,ep)-pcy1(ep)*cbrz(3,ep)
1046 c2=-pcy1(ep)*cbry(3,ep)-pcx1(ep)*cbrz(3,ep)
1047 mf33(1,3,ep)=mf33(1,3,ep)+c1
1048 mf33(2,3,ep)=mf33(2,3,ep)+c2
1049 mf13(1,3,ep)=mf13(1,3,ep)-c1
1050 mf13(2,3,ep)=mf13(2,3,ep)-c2
1051 c1=-pcx2(ep)*cbrx(4,ep)-pcy2(ep)*cbrz(4,ep)
1052 c2=-pcy2(ep)*cbry(4,ep)-pcx2(ep)*cbrz(4,ep)
1053 mf44(1,3,ep)=mf44(1,3,ep)+c1
1054 mf44(2,3,ep)=mf44(2,3,ep)+c2
1055 mf24(1,3,ep)=mf24(1,3,ep)-c1
1056 mf24(2,3,ep)=mf24(2,3,ep)-c2
1057 c1= pcx2(ep)*cbrx(1,ep)+pcy2(ep)*cbrz(1,ep)
1058 c2= pcy2(ep)*cbry(1,ep)+pcx2(ep)*cbrz(1,ep)
1059 fm12(3,1,ep)=fm12(3,1,ep)+c1
1060 fm12(3,2,ep)=fm12(3,2,ep)+c2
1061 fm14(3,1,ep)=fm14(3,1,ep)-c1
1062 fm14(3,2,ep)=fm14(3,2,ep)-c2
1063 ENDDO
1064 ENDIF !IF (IORTH==1)
1065C
1066 RETURN