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

Go to the source code of this file.

Functions/Subroutines

subroutine cbadef (jft, jlt, ng, vcore, area, cdet, vqn, vq, vjfi, vxyz, rxyz, vdef, vnrm, vastn, hx, hy, veta, vksi, bm, bmf, bf, bc, tc, nplat, iplat, isrot, brz)
subroutine cbadefsh (jft, jlt, x13, x24, y13, y24, bm, vdef, vxyz, nplat, iplat)
subroutine cbaderirz (jft, jlt, area, x13, x24, y13, y24, bm0rz, bmkrz, bmerz, vcore, nplat, iplat, ismstr)
subroutine cbadefrz (jft, jlt, area, rlz, vdef, vxyz, bm0rz, bmkrz, bmerz, vrlz, bmrz, brz, bm, nplat, iplat, ng)
subroutine cbaderirz0 (jft, jlt, area, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, bm0rz, bmkrz, bmerz, ismstr)
subroutine cbpatch (jft, jlt, bm0rz, bmkrz, bmerz, area, mx13, mx23, mx34, my13, my23, my34)
subroutine cbadef1 (jft, jlt, ng, vcore, vxyz, vdef, hx, hy, bm, nplat, iplat, isrot)
subroutine cbadeft (jft, jlt, vxyz, rlz, vdef, bm, nplat, iplat, isrot, bmrz)
subroutine cbadeft1 (jft, jlt, ng, vcore, vxyz, vdef, hx, hy, bm, nplat, iplat, isrot, bmrz, rxyz, wxy)
subroutine cbaderit1 (jft, jlt, ng, vcore, vq, vjfi, hx, hy, veta, vksi, bm, nplat, iplat, isrot)
subroutine cbaderirzt (jft, jlt, ng, bm0rz, bmkrz, bmerz, bmrz)
subroutine cbadeftw (jft, jlt, vxyz, rxyz, bm, bmf, bf, nplat, iplat, wxy)

Function/Subroutine Documentation

◆ cbadef()

subroutine cbadef ( integer jft,
integer jlt,
integer ng,
vcore,
area,
cdet,
vqn,
vq,
vjfi,
vxyz,
rxyz,
vdef,
vnrm,
vastn,
hx,
hy,
veta,
vksi,
bm,
bmf,
bf,
bc,
tc,
integer nplat,
integer, dimension(*) iplat,
integer isrot,
brz )

Definition at line 28 of file cbadef.F.

32C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
33C CALCUL DES DEFORMATIONS GENERALISEES AUX POINTS DE GAUSS ET [B]
34C ENTREES : NEL,NG,VCORE,VQN,VXYZ,RXYZ,VNRM,VASM,VASTN
35C SORTIES : VDEF,BM (MEMBRANE),BMF(COUPLAGE MEM-FLEXION),BF (FLEXION),BC(C.T)
36C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
37C-----------------------------------------------
38C I M P L I C I T T Y P E S
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C G l o b a l P a r a m e t e r s
43C-----------------------------------------------
44#include "mvsiz_p.inc"
45C-----------------------------------------------
46C D U M M Y A R G U M E N T S
47C-----------------------------------------------
48 INTEGER NPG,NG,JFT,JLT,NPLAT,IPLAT(*),ISROT
49 parameter(npg = 4)
50 my_real
51 . rxyz(mvsiz,8),vcore(mvsiz,12),vxyz(mvsiz,12),
52 . vqn(mvsiz,9,4),vksi(4,4),veta(4,4),
53 . bm(mvsiz,36),bmf(mvsiz,36),bf(mvsiz,24),bc(mvsiz,40),hx(mvsiz,4),hy(mvsiz,4)
54 my_real
55 . vnrm(mvsiz,12),vastn(mvsiz,16),vjfi(mvsiz,3,2,4),
56 . vq(mvsiz,3,3,4),vdef(mvsiz,8),dt1,area(*)
57 my_real
58 . cdet(*),tc(mvsiz,2,2),brz(mvsiz,4,4)
59C-----------------------------------------------
60C L O C A L V A R I A B L E S
61C-----------------------------------------------
62 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K
63 my_real
64 . vpg(2,npg),pg1,pg,thk,detj,det,
65 . tfn(3,2),bcx,bcy,bxy(3),byx(3),
66 . v1(2),v2(2),rv1,rv2,
67 . c1,c2,vt1,vt2,bc1,bc2,vb1,
68 . vjf1(2,3),vjf(3,3),tbi(2,2),tbc(2,2),v11(4)
69 my_real
70 . a_1,c11,c12,c21,c22,cc,beta1,ksi1,ksiy1,beta2,ksi2,ksiy2
71 parameter(pg=.577350269189626)
72 parameter(pg1=-.577350269189626)
73C--------------------------
74C INITIALISATION
75C--------------------------
76 DATA iun/1/
77 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
78C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
79#include "vectorize.inc"
80 DO i=jft,nplat
81 ep=iplat(i)
82C-------VCORE(12):BX0(2),BY0(2),GAMA(2),MX23,MY23,MX34,MY34,MX13,MY13
83C-------V13(I)->VXYZ(j,1,I),V24(I)->VXYZ(j,2,I),VHI(I)->VXYZ(j,3,I)------------
84 bm(ep,1)=vcore(ep,1)+hx(ep,ng)*vcore(ep,5)
85 bm(ep,2)=vcore(ep,2)+hx(ep,ng)*vcore(ep,6)
86 bm(ep,3)=hx(ep,ng)*fourth
87 bm(ep,4)=-bm(ep,3)
88 bm(ep,5)=vcore(ep,3)+hy(ep,ng)*vcore(ep,5)
89 bm(ep,6)=vcore(ep,4)+hy(ep,ng)*vcore(ep,6)
90 bm(ep,7)=hy(ep,ng)*fourth
91 bm(ep,8)=-bm(ep,7)
92C
93 vdef(ep,1)=bm(ep,1)*vxyz(ep,1)+bm(ep,2)*vxyz(ep,4)
94 1 +bm(ep,3)*vxyz(ep,7)
95 vdef(ep,2)=bm(ep,5)*vxyz(ep,2)+bm(ep,6)*vxyz(ep,5)
96 1 +bm(ep,7)*vxyz(ep,8)
97C--------FLEXION-----
98 vdef(ep,6)=bm(ep,1)*rxyz(ep,2)+bm(ep,2)*rxyz(ep,4)
99 1 +bm(ep,3)*rxyz(ep,6)
100 vdef(ep,7)=-(bm(ep,5)*rxyz(ep,1)+bm(ep,6)*rxyz(ep,3)
101 1 +bm(ep,7)*rxyz(ep,5))
102 vdef(ep,8)=-(bm(ep,1)*rxyz(ep,1)+bm(ep,2)*rxyz(ep,3)
103 1 +bm(ep,3)*rxyz(ep,5))
104 2 +bm(ep,5)*rxyz(ep,2)+bm(ep,6)*rxyz(ep,4)
105 3 +bm(ep,7)*rxyz(ep,6)
106C--------C.T-----
107C-------VCORE(12):BX0(2),BY0(2),GAMA(2),MX23,MY23,MX34,MY34,MX13,MY13
108 a_1 = 0.25/max(cdet(ep),em20)
109 c11=(vcore(ep,10)+vcore(ep,12)*vpg(1,ng))*a_1
110 c12=-(vcore(ep,8)+vcore(ep,12)*vpg(2,ng))*a_1
111 c21=-(vcore(ep,9)+vcore(ep,11)*vpg(1,ng))*a_1
112 c22=( vcore(ep,7)+vcore(ep,11)*vpg(2,ng))*a_1
113 beta1=vcore(ep,12)+vcore(ep,8)*vpg(2,ng)
114 ksi1=vcore(ep,12)+vcore(ep,10)*vpg(1,ng)
115 beta2=vcore(ep,11)+vcore(ep,7)*vpg(2,ng)
116 ksi2=vcore(ep,11)+vcore(ep,9)*vpg(1,ng)
117C-------NOEUD--- 1
118 bc(ep,1)=-c11-c12
119 bc(ep,2)=-c21-c22
120 bc(ep,3)= beta1*c11+ksi1*c12
121 bc(ep,4)= beta1*c21+ksi1*c22
122 bc(ep,5)= -beta2*c11-ksi2*c12
123 bc(ep,6)= -beta2*c21-ksi2*c22
124C-------NOEUD--- 2
125 bc(ep,7)= c11-c12
126 bc(ep,8)= c21-c22
127 bc(ep,9)= beta1*c11-ksi1*c12
128 bc(ep,10)=beta1*c21-ksi1*c22
129 bc(ep,11)=-beta2*c11+ksi2*c12
130 bc(ep,12)=-beta2*c21+ksi2*c22
131C
132 beta1=vcore(ep,8)+vcore(ep,12)*vpg(2,ng)
133 ksi1=vcore(ep,10)+vcore(ep,12)*vpg(1,ng)
134 beta2=vcore(ep,7)+vcore(ep,11)*vpg(2,ng)
135 ksi2=vcore(ep,9)+vcore(ep,11)*vpg(1,ng)
136C-----PARTIE SYM---------
137C-------NOEUD--- 1
138 bc(ep,13)=c11*vpg(2,ng)+c12*vpg(1,ng)
139 bc(ep,14)=c21*vpg(2,ng)+c22*vpg(1,ng)
140 bc(ep,15)=-beta1*c11-ksi1*c12
141 bc(ep,16)=-beta1*c21-ksi1*c22
142 bc(ep,17)=c11*beta2+c12*ksi2
143 bc(ep,18)=c21*beta2+c22*ksi2
144C-------NOEUD--- 2
145 bc(ep,19)=-bc(ep,13)
146 bc(ep,20)=-bc(ep,14)
147 bc(ep,21)=bc(ep,15)
148 bc(ep,22)=bc(ep,16)
149 bc(ep,23)=bc(ep,17)
150 bc(ep,24)=bc(ep,18)
151C
152 vdef(ep,4)=bc(ep,1)*vxyz(ep,3)+bc(ep,7)*vxyz(ep,6)
153 1 +bc(ep,13)*vxyz(ep,9)+bc(ep,3)*rxyz(ep,1)
154 2 +bc(ep,9)*rxyz(ep,3)+bc(ep,15)*rxyz(ep,7)
155 3 +bc(ep,5)*rxyz(ep,2)+bc(ep,11)*rxyz(ep,4)
156 4 +bc(ep,17)*rxyz(ep,8)
157 vdef(ep,5)=bc(ep,2)*vxyz(ep,3)+bc(ep,8)*vxyz(ep,6)
158 1 +bc(ep,14)*vxyz(ep,9)+bc(ep,4)*rxyz(ep,1)
159 2 +bc(ep,10)*rxyz(ep,3)+bc(ep,16)*rxyz(ep,7)
160 3 +bc(ep,6)*rxyz(ep,2)+bc(ep,12)*rxyz(ep,4)
161 4 +bc(ep,18)*rxyz(ep,8)
162 ENDDO
163 IF (isrot>0) THEN
164#include "vectorize.inc"
165 DO i=jft,nplat
166 ep=iplat(i)
167C----------no more constant shear-----------
168 vdef(ep,3)=bm(ep,1)*vxyz(ep,2)+bm(ep,2)*vxyz(ep,5)
169 1 +bm(ep,3)*vxyz(ep,8)
170 1 +bm(ep,5)*vxyz(ep,1)+bm(ep,6)*vxyz(ep,4)
171 1 +bm(ep,7)*vxyz(ep,7)
172 END DO
173 END IF !(ISROT>0) THEN
174C------------ELEMENT GAUCH-------------------
175#include "vectorize.inc"
176 DO 150 i=nplat+1,jlt
177 ep=iplat(i)
178C---------------------------------------------------
179C CALCUL DE [FN]
180C---------------------------------------------------
181 tfn(1,1)=vksi(1,ng)*vqn(ep,7,1)+vksi(2,ng)*vqn(ep,7,2)
182 1 +vksi(3,ng)*vqn(ep,7,3)+vksi(4,ng)*vqn(ep,7,4)
183 tfn(2,1)=vksi(1,ng)*vqn(ep,8,1)+vksi(2,ng)*vqn(ep,8,2)
184 1 +vksi(3,ng)*vqn(ep,8,3)+vksi(4,ng)*vqn(ep,8,4)
185 tfn(3,1)=vksi(1,ng)*vqn(ep,9,1)+vksi(2,ng)*vqn(ep,9,2)
186 1 +vksi(3,ng)*vqn(ep,9,3)+vksi(4,ng)*vqn(ep,9,4)
187 tfn(1,2)=veta(1,ng)*vqn(ep,7,1)+veta(2,ng)*vqn(ep,7,2)
188 1 +veta(3,ng)*vqn(ep,7,3)+veta(4,ng)*vqn(ep,7,4)
189 tfn(2,2)=veta(1,ng)*vqn(ep,8,1)+veta(2,ng)*vqn(ep,8,2)
190 1 +veta(3,ng)*vqn(ep,8,3)+veta(4,ng)*vqn(ep,8,4)
191 tfn(3,2)=veta(1,ng)*vqn(ep,9,1)+veta(2,ng)*vqn(ep,9,2)
192 1 +veta(3,ng)*vqn(ep,9,3)+veta(4,ng)*vqn(ep,9,4)
193C--------------------------------------------------------
194C CALCUL [BN]=[F]^-1 [FN]
195C--------------------------------------------------------
196 tbi(2,2)=vjfi(ep,1,1,ng)*tfn(1,1)+vjfi(ep,2,1,ng)*tfn(2,1)
197 1 + vjfi(ep,3,1,ng)*tfn(3,1)
198 tbi(2,1)=vjfi(ep,1,2,ng)*tfn(1,1)+vjfi(ep,2,2,ng)*tfn(2,1)
199 1 + vjfi(ep,3,2,ng)*tfn(3,1)
200 tbi(1,2)=vjfi(ep,1,1,ng)*tfn(1,2)+vjfi(ep,2,1,ng)*tfn(2,2)
201 1 + vjfi(ep,3,1,ng)*tfn(3,2)
202 tbi(1,1)=vjfi(ep,1,2,ng)*tfn(1,2)+vjfi(ep,2,2,ng)*tfn(2,2)
203 1 + vjfi(ep,3,2,ng)*tfn(3,2)
204C
205 thk =-(tbi(1,1)+tbi(2,2))
206C THK(2)=TBI(1,1)*TBI(2,2)-TBI(1,2)*TBI(2,1)
207 tbi(1,2)=-tbi(1,2)
208 tbi(2,1)=-tbi(2,1)
209C---------------------------
210C CALCUL DE TC=VJFI*VQ
211C---------------------------
212 tc(ep,1,1)=vjfi(ep,1,1,ng)*vq(ep,1,1,ng)+vjfi(ep,2,1,ng)
213 + *vq(ep,2,1,ng)+ vjfi(ep,3,1,ng)*vq(ep,3,1,ng)
214 tc(ep,2,1)=vjfi(ep,1,2,ng)*vq(ep,1,1,ng)+vjfi(ep,2,2,ng)
215 + *vq(ep,2,1,ng)+ vjfi(ep,3,2,ng)*vq(ep,3,1,ng)
216 tc(ep,1,2)=vjfi(ep,1,1,ng)*vq(ep,1,2,ng)+vjfi(ep,2,1,ng)
217 + *vq(ep,2,2,ng)+ vjfi(ep,3,1,ng)*vq(ep,3,2,ng)
218 tc(ep,2,2)=vjfi(ep,1,2,ng)*vq(ep,1,2,ng)+vjfi(ep,2,2,ng)
219 + *vq(ep,2,2,ng)+ vjfi(ep,3,2,ng)*vq(ep,3,2,ng)
220C------------------------------
221C CALCUL DE TBC=TBI*TC
222C------------------------------
223 tbc(1,1)=tbi(1,1)*tc(ep,1,1)+tbi(1,2)*tc(ep,2,1)
224 tbc(2,1)=tbi(2,1)*tc(ep,1,1)+tbi(2,2)*tc(ep,2,1)
225 tbc(1,2)=tbi(1,1)*tc(ep,1,2)+tbi(1,2)*tc(ep,2,2)
226 tbc(2,2)=tbi(2,1)*tc(ep,1,2)+tbi(2,2)*tc(ep,2,2)
227C--------------------------
228C CONSTRUIRE LA MATRICE [B] ET DEFORMATION
229C--------------------------
230C--------------------------
231C DEFORMATION IN-PLANE
232C--------------------------
233C [BM](3,3*NPG),[BMF](3,3*NPG),
234C [BF](3,2*NPG),VDEF(8),
235C--------------------------
236C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
237C---------- BOUCLE SUR LES NOEUDS---
238C--------J=1 II =0 KK = 0
239C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
240 vt1=vq(ep,1,1,ng)*vxyz(ep,1)+vq(ep,2,1,ng)*vxyz(ep,2)
241 1 +vq(ep,3,1,ng)*vxyz(ep,3)
242C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
243 vt2=vq(ep,1,2,ng)*vxyz(ep,1)+vq(ep,2,2,ng)*vxyz(ep,2)
244 1 +vq(ep,3,2,ng)*vxyz(ep,3)
245C--------TERMES DE [B0]=<T>CI-------
246 c1=vksi(1,ng)*tc(ep,1,1)+veta(1,ng)*tc(ep,2,1)
247 c2=vksi(1,ng)*tc(ep,1,2)+veta(1,ng)*tc(ep,2,2)
248C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
249 bc1=vksi(1,ng)*tbc(1,1)+veta(1,ng)*tbc(2,1)
250 bc2=vksi(1,ng)*tbc(1,2)+veta(1,ng)*tbc(2,2)
251C I = 1; JJ = 3*KK =0
252 IF (isrot>0) THEN
253 bxy(1)=vq(ep,1,1,ng)*c2
254 bxy(2)=vq(ep,2,1,ng)*c2
255 bxy(3)=vq(ep,3,1,ng)*c2
256 byx(1)=vq(ep,1,2,ng)*c1
257 byx(2)=vq(ep,2,2,ng)*c1
258 byx(3)=vq(ep,3,2,ng)*c1
259 brz(i,1,1)=-bxy(1)+byx(1)
260 brz(i,2,1)=-bxy(2)+byx(2)
261 brz(i,3,1)=-bxy(3)+byx(3)
262 bm(ep,3)=bxy(1)+byx(1)
263 bm(ep,6)=bxy(2)+byx(2)
264 bm(ep,9)=bxy(3)+byx(3)
265 vdef(ep,3)= c1*vt2+c2*vt1
266 ELSE
267 bm(ep,9)=zero
268 END IF
269C---------CONTRIBUTION DU <T1> :EPS-X
270 bm(ep,1)=vq(ep,1,1,ng)*c1
271 bmf(ep,1)=thk*bm(ep,1)+vq(ep,1,1,ng)*bc1
272C---------CONTRIBUTION DU <T2> :EPS-Y
273 bm(ep,2)=vq(ep,1,2,ng)*c2
274 bmf(ep,2)=thk*bm(ep,2)+vq(ep,1,2,ng)*bc2
275C---------CONTRIBUTION DU <T> :EPS-XY
276 bmf(ep,3)=thk*bm(ep,3)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
277C I = 2; JJ = 3*(KK+1)=3
278C---------:EPS-X :EPS-Y :EPS-XY
279 bm(ep,4)=vq(ep,2,1,ng)*c1
280 bmf(ep,4)=thk*bm(ep,4)+vq(ep,2,1,ng)*bc1
281 bm(ep,5)=vq(ep,2,2,ng)*c2
282 bmf(ep,5)=thk*bm(ep,5)+vq(ep,2,2,ng)*bc2
283 bmf(ep,6)=thk*bm(ep,6)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
284C I = 3; JJ = 3*(KK+2)=6
285C---------:EPS-X:EPS-Y:EPS-XY
286 bm(ep,7)=vq(ep,3,1,ng)*c1
287 bmf(ep,7)=thk*bm(ep,7)+vq(ep,3,1,ng)*bc1
288 bm(ep,8)=vq(ep,3,2,ng)*c2
289 bmf(ep,8)=thk*bm(ep,8)+vq(ep,3,2,ng)*bc2
290 bmf(ep,9)=thk*bm(ep,9)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
291C
292C---------DEF MEMBRANE------
293 vdef(ep,1)= c1*vt1
294 vdef(ep,2)= c2*vt2
295C
296C---------FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
297C----------[B01]---------
298 vdef(ep,6)= bc1*vt1
299 vdef(ep,7)= bc2*vt2
300 vdef(ep,8)= bc1*vt2+bc2*vt1
301C---------- TERMES ASSOCIEES AU BETA ----------
302C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
303 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,1)+vq(ep,2,1,ng)*vqn(ep,2,1)
304 1 +vq(ep,3,1,ng)*vqn(ep,3,1)
305 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,1)+vq(ep,2,1,ng)*vqn(ep,5,1)
306 1 +vq(ep,3,1,ng)*vqn(ep,6,1))
307 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,1)+vq(ep,2,2,ng)*vqn(ep,2,1)
308 1 +vq(ep,3,2,ng)*vqn(ep,3,1)
309 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,1)+vq(ep,2,2,ng)*vqn(ep,5,1)
310 1 +vq(ep,3,2,ng)*vqn(ep,6,1))
311C----------RV1=<V1>*RXYZ AU CHAQUE NOEUD---------
312 rv1=v1(1)*rxyz(ep,1)+v1(2)*rxyz(ep,2)
313C----------RV2=<V2>*RXYZ AU CHAQUE NOEUD---------
314 rv2=v2(1)*rxyz(ep,1)+v2(2)*rxyz(ep,2)
315C----------TERMES ASSOCIEES AUX BETA DE [B01]--
316C JJ = 3*II =0
317 bf(ep,1)=v1(1)*c1
318 bf(ep,2)=v2(1)*c2
319 bf(ep,3)=v1(1)*c2+v2(1)*c1
320 bf(ep,4)=v1(2)*c1
321 bf(ep,5)=v2(2)*c2
322 bf(ep,6)=v1(2)*c2+v2(2)*c1
323C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
324 vdef(ep,6)= vdef(ep,6)+c1*rv1
325 vdef(ep,7)= vdef(ep,7)+c2*rv2
326 vdef(ep,8)= vdef(ep,8)+c1*rv2+c2*rv1
327C
328C--------J=2---II=(J-1)*2 =2 KK = 3*(J-1) =3
329C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
330 vt1=vq(ep,1,1,ng)*vxyz(ep,4)+vq(ep,2,1,ng)*vxyz(ep,5)
331 1 +vq(ep,3,1,ng)*vxyz(ep,6)
332C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
333 vt2=vq(ep,1,2,ng)*vxyz(ep,4)+vq(ep,2,2,ng)*vxyz(ep,5)
334 1 +vq(ep,3,2,ng)*vxyz(ep,6)
335C--------- [B0], [B01]---------
336C--------TERMES DE [B0]=<T>CI-------
337 c1=vksi(2,ng)*tc(ep,1,1)+veta(2,ng)*tc(ep,2,1)
338 c2=vksi(2,ng)*tc(ep,1,2)+veta(2,ng)*tc(ep,2,2)
339C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
340 bc1=vksi(2,ng)*tbc(1,1)+veta(2,ng)*tbc(2,1)
341 bc2=vksi(2,ng)*tbc(1,2)+veta(2,ng)*tbc(2,2)
342C
343C I = 1; JJ = 3*KK =9
344 IF (isrot>0) THEN
345 bxy(1)=vq(ep,1,1,ng)*c2
346 bxy(2)=vq(ep,2,1,ng)*c2
347 bxy(3)=vq(ep,3,1,ng)*c2
348 byx(1)=vq(ep,1,2,ng)*c1
349 byx(2)=vq(ep,2,2,ng)*c1
350 byx(3)=vq(ep,3,2,ng)*c1
351 brz(i,1,2)=-bxy(1)+byx(1)
352 brz(i,2,2)=-bxy(2)+byx(2)
353 brz(i,3,2)=-bxy(3)+byx(3)
354 bm(ep,12)=bxy(1)+byx(1)
355 bm(ep,15)=bxy(2)+byx(2)
356 bm(ep,18)=bxy(3)+byx(3)
357 vdef(ep,3)= vdef(ep,3)+ c1*vt2+c2*vt1
358 ELSE
359 bm(ep,18)=zero
360 END IF
361C---------:EPS-X:EPS-Y:EPS-XY
362 bm(ep,10)=vq(ep,1,1,ng)*c1
363 bmf(ep,10)=thk*bm(ep,10)+vq(ep,1,1,ng)*bc1
364 bm(ep,11)=vq(ep,1,2,ng)*c2
365 bmf(ep,11)=thk*bm(ep,11)+vq(ep,1,2,ng)*bc2
366 bmf(ep,12)=thk*bm(ep,12)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
367C
368C I = 2; JJ = 3*(KK+1)=12
369C---------:EPS-X:EPS-Y:EPS-XY
370 bm(ep,13)=vq(ep,2,1,ng)*c1
371 bmf(ep,13)=thk*bm(ep,13)+vq(ep,2,1,ng)*bc1
372 bm(ep,14)=vq(ep,2,2,ng)*c2
373 bmf(ep,14)=thk*bm(ep,14)+vq(ep,2,2,ng)*bc2
374 bmf(ep,15)=thk*bm(ep,15)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
375C
376C I = 3; JJ = 3*(KK+2)=15
377C---------:EPS-X:EPS-Y:EPS-XY
378 bm(ep,16)=vq(ep,3,1,ng)*c1
379 bmf(ep,16)=thk*bm(ep,16)+vq(ep,3,1,ng)*bc1
380 bm(ep,17)=vq(ep,3,2,ng)*c2
381 bmf(ep,17)=thk*bm(ep,17)+vq(ep,3,2,ng)*bc2
382 bmf(ep,18)=thk*bm(ep,18)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
383C
384C---------DEF MEMBRANE------
385 vdef(ep,1)= vdef(ep,1)+ c1*vt1
386 vdef(ep,2)= vdef(ep,2)+ c2*vt2
387C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
388C----------[B01]---------
389 vdef(ep,6)= vdef(ep,6)+ bc1*vt1
390 vdef(ep,7)= vdef(ep,7)+ bc2*vt2
391 vdef(ep,8)= vdef(ep,8)+ bc1*vt2+bc2*vt1
392C---------- TERMES ASSOCIEES AU BETA ----------
393C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
394 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,2)+vq(ep,2,1,ng)*vqn(ep,2,2)
395 1 +vq(ep,3,1,ng)*vqn(ep,3,2)
396 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,2)+vq(ep,2,1,ng)*vqn(ep,5,2)
397 1 +vq(ep,3,1,ng)*vqn(ep,6,2))
398 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,2)+vq(ep,2,2,ng)*vqn(ep,2,2)
399 1 +vq(ep,3,2,ng)*vqn(ep,3,2)
400 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,2)+vq(ep,2,2,ng)*vqn(ep,5,2)
401 1 +vq(ep,3,2,ng)*vqn(ep,6,2))
402C----------RV1=<V1>*RXYZ AU CHAQUE NOEUD---------
403 rv1=v1(1)*rxyz(ep,2+1)+v1(2)*rxyz(ep,2+2)
404C----------RV2=<V2>*RXYZ AU CHAQUE NOEUD---------
405 rv2=v2(1)*rxyz(ep,2+1)+v2(2)*rxyz(ep,2+2)
406C----------TERMES ASSOCIEES AUX BETA DE [B01]--
407C JJ = 3*II=6
408 bf(ep,7)=v1(1)*c1
409 bf(ep,8)=v2(1)*c2
410 bf(ep,9)=v1(1)*c2+v2(1)*c1
411 bf(ep,10)=v1(2)*c1
412 bf(ep,11)=v2(2)*c2
413 bf(ep,12)=v1(2)*c2+v2(2)*c1
414C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
415 vdef(ep,6)= vdef(ep,6)+c1*rv1
416 vdef(ep,7)= vdef(ep,7)+c2*rv2
417 vdef(ep,8)= vdef(ep,8)+c1*rv2+c2*rv1
418C
419C--------J=3---II=(J-1)*2 =4 KK = 3*(J-1)=6
420C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
421 vt1=vq(ep,1,1,ng)*vxyz(ep,7)+vq(ep,2,1,ng)*vxyz(ep,8)
422 1 +vq(ep,3,1,ng)*vxyz(ep,9)
423C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
424 vt2=vq(ep,1,2,ng)*vxyz(ep,7)+vq(ep,2,2,ng)*vxyz(ep,8)
425 1 +vq(ep,3,2,ng)*vxyz(ep,9)
426C--------- [B0], [B01]---------
427C--------TERMES DE [B0]=<T>CI-------
428 c1=vksi(3,ng)*tc(ep,1,1)+veta(3,ng)*tc(ep,2,1)
429 c2=vksi(3,ng)*tc(ep,1,2)+veta(3,ng)*tc(ep,2,2)
430C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
431 bc1=vksi(3,ng)*tbc(1,1)+veta(3,ng)*tbc(2,1)
432 bc2=vksi(3,ng)*tbc(1,2)+veta(3,ng)*tbc(2,2)
433C I = 1;JJ = 3*KK=18
434 IF (isrot>0) THEN
435 bxy(1)=vq(ep,1,1,ng)*c2
436 bxy(2)=vq(ep,2,1,ng)*c2
437 bxy(3)=vq(ep,3,1,ng)*c2
438 byx(1)=vq(ep,1,2,ng)*c1
439 byx(2)=vq(ep,2,2,ng)*c1
440 byx(3)=vq(ep,3,2,ng)*c1
441 brz(i,1,3)=-bxy(1)+byx(1)
442 brz(i,2,3)=-bxy(2)+byx(2)
443 brz(i,3,3)=-bxy(3)+byx(3)
444 bm(ep,21)=bxy(1)+byx(1)
445 bm(ep,24)=bxy(2)+byx(2)
446 bm(ep,27)=bxy(3)+byx(3)
447 vdef(ep,3)= vdef(ep,3)+ c1*vt2+c2*vt1
448 ELSE
449 bm(ep,27)=zero
450 END IF
451C---------:EPS-X:EPS-Y:EPS-XY
452 bm(ep,19)=vq(ep,1,1,ng)*c1
453 bmf(ep,19)=thk*bm(ep,19)+vq(ep,1,1,ng)*bc1
454 bm(ep,20)=vq(ep,1,2,ng)*c2
455 bmf(ep,20)=thk*bm(ep,20)+vq(ep,1,2,ng)*bc2
456 bmf(ep,21)=thk*bm(ep,21)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
457C
458C I = 2; JJ = 3*(KK+1)=21
459 bm(ep,22)=vq(ep,2,1,ng)*c1
460 bmf(ep,22)=thk*bm(ep,22)+vq(ep,2,1,ng)*bc1
461 bm(ep,23)=vq(ep,2,2,ng)*c2
462 bmf(ep,23)=thk*bm(ep,23)+vq(ep,2,2,ng)*bc2
463 bmf(ep,24)=thk*bm(ep,24)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
464C
465C I = 3; JJ = 3*(KK+2)=24
466 bm(ep,25)=vq(ep,3,1,ng)*c1
467 bmf(ep,25)=thk*bm(ep,25)+vq(ep,3,1,ng)*bc1
468 bm(ep,26)=vq(ep,3,2,ng)*c2
469 bmf(ep,26)=thk*bm(ep,26)+vq(ep,3,2,ng)*bc2
470 bmf(ep,27)=thk*bm(ep,27)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
471C
472C---------DEF MEMBRANE------
473 vdef(ep,1)= vdef(ep,1)+ c1*vt1
474 vdef(ep,2)= vdef(ep,2)+ c2*vt2
475C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
476C----------[B01]---------
477 vdef(ep,6)= vdef(ep,6)+ bc1*vt1
478 vdef(ep,7)= vdef(ep,7)+ bc2*vt2
479 vdef(ep,8)= vdef(ep,8)+ bc1*vt2+bc2*vt1
480C---------- TERMES ASSOCIEES AU BETA ----------
481C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
482 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,3)+vq(ep,2,1,ng)*vqn(ep,2,3)
483 1 +vq(ep,3,1,ng)*vqn(ep,3,3)
484 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,3)+vq(ep,2,1,ng)*vqn(ep,5,3)
485 1 +vq(ep,3,1,ng)*vqn(ep,6,3))
486 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,3)+vq(ep,2,2,ng)*vqn(ep,2,3)
487 1 +vq(ep,3,2,ng)*vqn(ep,3,3)
488 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,3)+vq(ep,2,2,ng)*vqn(ep,5,3)
489 1 +vq(ep,3,2,ng)*vqn(ep,6,3))
490C----------RV1=<V1>*RXYZ AU CHAQUE NOEUD---------
491 rv1=v1(1)*rxyz(ep,4+1)+v1(2)*rxyz(ep,4+2)
492C----------RV2=<V2>*RXYZ AU CHAQUE NOEUD---------
493 rv2=v2(1)*rxyz(ep,4+1)+v2(2)*rxyz(ep,4+2)
494C----------TERMES ASSOCIEES AUX BETA DE [B01]--
495C JJ = 3*II =12
496 bf(ep,13)=v1(1)*c1
497 bf(ep,14)=v2(1)*c2
498 bf(ep,15)=v1(1)*c2+v2(1)*c1
499 bf(ep,16)=v1(2)*c1
500 bf(ep,17)=v2(2)*c2
501 bf(ep,18)=v1(2)*c2+v2(2)*c1
502C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
503 vdef(ep,6)= vdef(ep,6)+c1*rv1
504 vdef(ep,7)= vdef(ep,7)+c2*rv2
505 vdef(ep,8)= vdef(ep,8)+c1*rv2+c2*rv1
506C--------J=4---II=(J-1)*2 =6 KK = 3*(J-1)=9
507C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
508 vt1=vq(ep,1,1,ng)*vxyz(ep,10)+vq(ep,2,1,ng)*vxyz(ep,11)
509 1 +vq(ep,3,1,ng)*vxyz(ep,12)
510C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
511 vt2=vq(ep,1,2,ng)*vxyz(ep,10)+vq(ep,2,2,ng)*vxyz(ep,11)
512 1 +vq(ep,3,2,ng)*vxyz(ep,12)
513C--------- [B0], [B01]---------
514C--------TERMES DE [B0]=<T>CI-------
515 c1=vksi(4,ng)*tc(ep,1,1)+veta(4,ng)*tc(ep,2,1)
516 c2=vksi(4,ng)*tc(ep,1,2)+veta(4,ng)*tc(ep,2,2)
517C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
518 bc1=vksi(4,ng)*tbc(1,1)+veta(4,ng)*tbc(2,1)
519 bc2=vksi(4,ng)*tbc(1,2)+veta(4,ng)*tbc(2,2)
520C I = 1; JJ = 3*KK =27
521 IF (isrot>0) THEN
522 bxy(1)=vq(ep,1,1,ng)*c2
523 bxy(2)=vq(ep,2,1,ng)*c2
524 bxy(3)=vq(ep,3,1,ng)*c2
525 byx(1)=vq(ep,1,2,ng)*c1
526 byx(2)=vq(ep,2,2,ng)*c1
527 byx(3)=vq(ep,3,2,ng)*c1
528 brz(i,1,4)=-bxy(1)+byx(1)
529 brz(i,2,4)=-bxy(2)+byx(2)
530 brz(i,3,4)=-bxy(3)+byx(3)
531 bm(ep,30)=bxy(1)+byx(1)
532 bm(ep,33)=bxy(2)+byx(2)
533 bm(ep,36)=bxy(3)+byx(3)
534 vdef(ep,3)= vdef(ep,3)+ c1*vt2+c2*vt1
535 ELSE
536 bm(ep,36)=zero
537 END IF
538 bm(ep,28)=vq(ep,1,1,ng)*c1
539 bmf(ep,28)=thk*bm(ep,28)+vq(ep,1,1,ng)*bc1
540 bm(ep,29)=vq(ep,1,2,ng)*c2
541 bmf(ep,29)=thk*bm(ep,29)+vq(ep,1,2,ng)*bc2
542C---------CONTRIBUTION DU <T> :EPS-XY
543 bmf(ep,30)=thk*bm(ep,30)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
544C
545C I = 2; JJ = 3*(KK+1)=30
546 bm(ep,31)=vq(ep,2,1,ng)*c1
547 bmf(ep,31)=thk*bm(ep,31)+vq(ep,2,1,ng)*bc1
548 bm(ep,32)=vq(ep,2,2,ng)*c2
549 bmf(ep,32)=thk*bm(ep,32)+vq(ep,2,2,ng)*bc2
550 bmf(ep,33)=thk*bm(ep,33)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
551C
552C I = 3; JJ = 3*(KK+2)=33
553 bm(ep,34)=vq(ep,3,1,ng)*c1
554 bmf(ep,34)=thk*bm(ep,34)+vq(ep,3,1,ng)*bc1
555 bm(ep,35)=vq(ep,3,2,ng)*c2
556 bmf(ep,35)=thk*bm(ep,35)+vq(ep,3,2,ng)*bc2
557 bmf(ep,36)=thk*bm(ep,36)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
558C
559C---------DEF MEMBRANE------
560 vdef(ep,1)= vdef(ep,1)+ c1*vt1
561 vdef(ep,2)= vdef(ep,2)+ c2*vt2
562C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
563C----------SEULEMENT LA PARTIE DE [B01]---------
564 vdef(ep,6)= vdef(ep,6)+ bc1*vt1
565 vdef(ep,7)= vdef(ep,7)+ bc2*vt2
566 vdef(ep,8)= vdef(ep,8)+ bc1*vt2+bc2*vt1
567C---------- TERMES ASSOCIEES AU BETA ----------
568C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
569 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,4)+vq(ep,2,1,ng)*vqn(ep,2,4)
570 1 +vq(ep,3,1,ng)*vqn(ep,3,4)
571 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,4)+vq(ep,2,1,ng)*vqn(ep,5,4)
572 1 +vq(ep,3,1,ng)*vqn(ep,6,4))
573 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,4)+vq(ep,2,2,ng)*vqn(ep,2,4)
574 1 +vq(ep,3,2,ng)*vqn(ep,3,4)
575 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,4)+vq(ep,2,2,ng)*vqn(ep,5,4)
576 1 +vq(ep,3,2,ng)*vqn(ep,6,4))
577C----------RV1=<V1>*RXYZ AU CHAQUE NOEUD---------
578 rv1=v1(1)*rxyz(ep,6+1)+v1(2)*rxyz(ep,6+2)
579C----------RV2=<V2>*RXYZ AU CHAQUE NOEUD---------
580 rv2=v2(1)*rxyz(ep,6+1)+v2(2)*rxyz(ep,6+2)
581C----------TERMES ASSOCIEES AUX BETA DE [B01]--
582C JJ = 3*II=18
583 bf(ep,19)=v1(1)*c1
584 bf(ep,20)=v2(1)*c2
585 bf(ep,21)=v1(1)*c2+v2(1)*c1
586 bf(ep,22)=v1(2)*c1
587 bf(ep,23)=v2(2)*c2
588 bf(ep,24)=v1(2)*c2+v2(2)*c1
589C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
590 vdef(ep,6)= vdef(ep,6)+thk*vdef(ep,1)+c1*rv1
591 vdef(ep,7)= vdef(ep,7)+thk*vdef(ep,2)+c2*rv2
592 vdef(ep,8)= vdef(ep,8)+thk*vdef(ep,3)+c1*rv2+c2*rv1
593C
594C--------------------------
595C DEFORMATION OUT-PLANE (C.T)
596C--------------------------
597C [BC](2,5*NPG),VDEF(4,5)
598C--------------------------
599C
600 v11(1)=vksi(2,ng)
601 v11(2)=vksi(3,ng)
602 v11(3)=veta(4,ng)
603 v11(4)=veta(3,ng)
604C--- J=1---------------------
605C-----------POINTEURS POUR BC(II,J1~J2)
606C J1=(KA1(1,J)-1)*5=0
607C J2=(KA1(2,J)-1)*5=5
608C II=KIND(J) =1
609C JJ = 3*(J-1) =0
610C--------TERMES <NAI>-------
611 c1=v11(1)*vnrm(ep,1)
612 bc(ep,1)=-c1
613 bc(ep,11)=c1
614C
615 c1=v11(1)*vnrm(ep,2)
616 bc(ep,3)=-c1
617 bc(ep,13)=c1
618C
619 c1=v11(1)*vnrm(ep,3)
620 bc(ep,5)=-c1
621 bc(ep,15)=c1
622C--------TERMES <ASI>-------
623 bc(ep,7)=v11(1)*vastn(ep,1)
624 bc(ep,9)=v11(1)*vastn(ep,2)
625 bc(ep,17)=v11(1)*vastn(ep,3)
626 bc(ep,19)=v11(1)*vastn(ep,4)
627C--- J=2---------------------
628C-----------POINTEURS POUR BC(II,J1~J2)
629C J1=(KA1(1,J)-1)*5 =15
630C J2=(KA1(2,J)-1)*5 =10
631C II=KIND(J) =1
632C JJ = 3*(J-1) =3
633C--------TERMES <NAI>-------
634 c1=v11(2)*vnrm(ep,4)
635 bc(ep,31)=-c1
636 bc(ep,21)=c1
637C
638 c1=v11(2)*vnrm(ep,5)
639 bc(ep,33)=-c1
640 bc(ep,23)=c1
641C
642 c1=v11(2)*vnrm(ep,6)
643 bc(ep,35)=-c1
644 bc(ep,25)=c1
645C--------TERMES <ASI>-------
646 bc(ep,37)=v11(2)*vastn(ep,5)
647 bc(ep,39)=v11(2)*vastn(ep,6)
648 bc(ep,27)=v11(2)*vastn(ep,7)
649 bc(ep,29)=v11(2)*vastn(ep,8)
650C--- J=3---------------------
651C-----------POINTEURS POUR BC(II,J1~J2)
652C J1=(KA1(1,J)-1)*5 =0
653C J2=(KA1(2,J)-1)*5 =15
654C II=KIND(J) =2
655C JJ = 3*(J-1) =6
656C--------TERMES <NAI>-------
657 c1=v11(3)*vnrm(ep,7)
658 bc(ep,2)=-c1
659 bc(ep,32)=c1
660C
661 c1=v11(3)*vnrm(ep,8)
662 bc(ep,4)=-c1
663 bc(ep,34)=c1
664C
665 c1=v11(3)*vnrm(ep,9)
666 bc(ep,6)=-c1
667 bc(ep,36)=c1
668C--------TERMES <ASI>-------
669 bc(ep,8)=v11(3)*vastn(ep,9)
670 bc(ep,10)=v11(3)*vastn(ep,10)
671 bc(ep,38)=v11(3)*vastn(ep,11)
672 bc(ep,40)=v11(3)*vastn(ep,12)
673C--- J=4---------------------
674C-----------POINTEURS POUR BC(II,J1~J2)
675C J1=(KA1(1,J)-1)*5 =5
676C J2=(KA1(2,J)-1)*5 =10
677C II=KIND(J) =2
678C JJ = 3*(J-1) =9
679C--------TERMES <NAI>-------
680 c1=v11(4)*vnrm(ep,10)
681 bc(ep,12)=-c1
682 bc(ep,22)=c1
683C
684 c1=v11(4)*vnrm(ep,11)
685 bc(ep,14)=-c1
686 bc(ep,24)=c1
687C
688 c1=v11(4)*vnrm(ep,12)
689 bc(ep,16)=-c1
690 bc(ep,26)=c1
691C--------TERMES <ASI>-------
692 bc(ep,18)=v11(4)*vastn(ep,13)
693 bc(ep,20)=v11(4)*vastn(ep,14)
694 bc(ep,28)=v11(4)*vastn(ep,15)
695 bc(ep,30)=v11(4)*vastn(ep,16)
696C
697C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
698C----------J=1--II=(J-1)*3 =0 JJ=(J-1)*2 =0 KK=10*(J-1)=0
699 bcx=
700 1 bc(ep,1)*vxyz(ep,1)+bc(ep,3)*vxyz(ep,2)
701 2 +bc(ep,5)*vxyz(ep,3)+bc(ep,7)*rxyz(ep,1)
702 3 +bc(ep,9)*rxyz(ep,2)
703 bcy=
704 1 bc(ep,2)*vxyz(ep,1)+bc(ep,4)*vxyz(ep,2)
705 2 +bc(ep,6)*vxyz(ep,3)+bc(ep,8)*rxyz(ep,1)
706 3 +bc(ep,10)*rxyz(ep,2)
707C----------J=2 II=(J-1)*3 =3 JJ=(J-1)*2 =2 KK=10*(J-1) =10
708 bcx=bcx
709 1 +bc(ep,11)*vxyz(ep,4)+bc(ep,13)*vxyz(ep,5)
710 2 +bc(ep,15)*vxyz(ep,6)+bc(ep,17)*rxyz(ep,3)
711 3 +bc(ep,19)*rxyz(ep,4)
712 bcy=bcy
713 1 +bc(ep,12)*vxyz(ep,4)+bc(ep,14)*vxyz(ep,5)
714 2 +bc(ep,16)*vxyz(ep,6)+bc(ep,18)*rxyz(ep,3)
715 3 +bc(ep,20)*rxyz(ep,4)
716C----------J=3-- II=(J-1)*3=6 JJ=(J-1)*2=4 KK=10*(J-1)=20
717 bcx=bcx
718 1 +bc(ep,21)*vxyz(ep,7)+bc(ep,23)*vxyz(ep,8)
719 2 +bc(ep,25)*vxyz(ep,9)+bc(ep,27)*rxyz(ep,5)
720 3 +bc(ep,29)*rxyz(ep,6)
721 bcy=bcy
722 1 +bc(ep,22)*vxyz(ep,7)+bc(ep,24)*vxyz(ep,8)
723 2 +bc(ep,26)*vxyz(ep,9)+bc(ep,28)*rxyz(ep,5)
724 3 +bc(ep,30)*rxyz(ep,6)
725C----------J=4-- II=(J-1)*3 =9 JJ=(J-1)*2 =6 KK=10*(J-1)=30
726 bcx=bcx
727 1 +bc(ep,31)*vxyz(ep,10)+bc(ep,33)*vxyz(ep,11)
728 2 +bc(ep,35)*vxyz(ep,12)+bc(ep,37)*rxyz(ep,7)
729 3 +bc(ep,39)*rxyz(ep,8)
730 bcy=bcy
731 1 +bc(ep,32)*vxyz(ep,10)+bc(ep,34)*vxyz(ep,11)
732 2 +bc(ep,36)*vxyz(ep,12)+bc(ep,38)*rxyz(ep,7)
733 3 +bc(ep,40)*rxyz(ep,8)
734 vdef(ep,4)=tc(ep,1,1)*bcx+tc(ep,2,1)*bcy
735 vdef(ep,5)=tc(ep,1,2)*bcx+tc(ep,2,2)*bcy
736C
737 150 CONTINUE
738 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21

◆ cbadef1()

subroutine cbadef1 ( integer jft,
integer jlt,
integer ng,
vcore,
vxyz,
vdef,
hx,
hy,
bm,
integer nplat,
integer, dimension(*) iplat,
integer isrot )

Definition at line 1220 of file cbadef.F.

1222C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1223C CALCUL DES DEFORMATIONS GENERALISEES AUX POINTS DE GAUSS ET [B] (MEMBRANE ONLY)
1224C ENTREES : NEL,NG,VCORE,VXYZ
1225C SORTIES : VDEF,BM (MEMBRANE)
1226C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1227C-----------------------------------------------
1228C I M P L I C I T T Y P E S
1229C-----------------------------------------------
1230#include "implicit_f.inc"
1231C-----------------------------------------------
1232C G l o b a l P a r a m e t e r s
1233C-----------------------------------------------
1234#include "mvsiz_p.inc"
1235C-----------------------------------------------
1236C D U M M Y A R G U M E N T S
1237C-----------------------------------------------
1238 INTEGER NG,JFT,JLT,NPLAT,IPLAT(*),ISROT
1239 my_real
1240 . vcore(mvsiz,12),vxyz(mvsiz,12),bm(mvsiz,36),hx(mvsiz,4),hy(mvsiz,4)
1241 my_real
1242 . vdef(mvsiz,8)
1243C-----------------------------------------------
1244C L O C A L V A R I A B L E S
1245C-----------------------------------------------
1246 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K
1247C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1248#include "vectorize.inc"
1249 DO i=jft,nplat
1250 ep=iplat(i)
1251C-------VCORE(12):BX0(2),BY0(2),GAMA(2),MX23,MY23,MX34,MY34,MX13,MY13
1252C-------V13(I)->VXYZ(j,1,I),V24(I)->VXYZ(j,2,I),VHI(I)->VXYZ(j,3,I)------------
1253 bm(ep,1)=vcore(ep,1)+hx(ep,ng)*vcore(ep,5)
1254 bm(ep,2)=vcore(ep,2)+hx(ep,ng)*vcore(ep,6)
1255 bm(ep,3)=hx(ep,ng)*fourth
1256 bm(ep,4)=-bm(ep,3)
1257 bm(ep,5)=vcore(ep,3)+hy(ep,ng)*vcore(ep,5)
1258 bm(ep,6)=vcore(ep,4)+hy(ep,ng)*vcore(ep,6)
1259 bm(ep,7)=hy(ep,ng)*fourth
1260 bm(ep,8)=-bm(ep,7)
1261C
1262 vdef(ep,1)=bm(ep,1)*vxyz(ep,1)+bm(ep,2)*vxyz(ep,4)
1263 1 +bm(ep,3)*vxyz(ep,7)
1264 vdef(ep,2)=bm(ep,5)*vxyz(ep,2)+bm(ep,6)*vxyz(ep,5)
1265 1 +bm(ep,7)*vxyz(ep,8)
1266 vdef(ep,4)=zero
1267 vdef(ep,5)=zero
1268 vdef(ep,6)=zero
1269 vdef(ep,7)=zero
1270 vdef(ep,8)=zero
1271 ENDDO
1272 IF (isrot>0) THEN
1273#include "vectorize.inc"
1274 DO i=jft,nplat
1275 ep=iplat(i)
1276C----------no more constant shear-----------
1277 vdef(ep,3)=bm(ep,1)*vxyz(ep,2)+bm(ep,2)*vxyz(ep,5)
1278 1 +bm(ep,3)*vxyz(ep,8)
1279 1 +bm(ep,5)*vxyz(ep,1)+bm(ep,6)*vxyz(ep,4)
1280 1 +bm(ep,7)*vxyz(ep,7)
1281 END DO
1282 END IF !(ISROT>0) THEN
1283C
1284 RETURN

◆ cbadefrz()

subroutine cbadefrz ( integer jft,
integer jlt,
area,
rlz,
vdef,
vxyz,
bm0rz,
bmkrz,
bmerz,
vrlz,
bmrz,
brz,
bm,
integer nplat,
integer, dimension(*) iplat,
integer ng )

Definition at line 866 of file cbadef.F.

870C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
871#include "implicit_f.inc"
872#include "mvsiz_p.inc"
873C-----------------------------------------------
874C D U M M Y A R G U M E N T S
875C-----------------------------------------------
876 INTEGER JFT,JLT,NPLAT,IPLAT(*),NG
877 my_real
878 . rlz(mvsiz,4),vxyz(mvsiz,3,4),
879 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
880 . bmrz(mvsiz,3,4),bm(mvsiz,36)
881 my_real
882 . area(*),vdef(mvsiz,8),vrlz(*),brz(mvsiz,4,4)
883C-----------------------------------------------
884C L O C A L V A R I A B L E S
885C-----------------------------------------------
886 INTEGER I,J,EP,NPG
887 parameter(npg = 4)
888 my_real
889 . aa4,bxv2,byv1,erz,vpg(2,npg),pg,pg1
890 parameter(pg=.577350269189626)
891 parameter(pg1=-.577350269189626)
892C--------------------------
893C INITIALISATION
894C--------------------------
895 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
896C--------pay attention QEPH: Bi*A--- RLZ: VRZ/A-- QBAT:Bi--RLZ:VRZ
897 DO i=jft,jlt
898 DO j=1,4
899 bmrz(i,1,j) = bm0rz(i,1,j)+bmkrz(i,1,j)*vpg(1,ng)+
900 . bmerz(i,1,j)*vpg(2,ng)
901 bmrz(i,2,j) = bm0rz(i,2,j)+bmkrz(i,2,j)*vpg(1,ng)+
902 . bmerz(i,2,j)*vpg(2,ng)
903 bmrz(i,3,j) = bm0rz(i,3,j)+bmkrz(i,3,j)*vpg(1,ng)+
904 . bmerz(i,3,j)*vpg(2,ng)
905C
906 brz(i,4,j) = bm0rz(i,4,j)+bmkrz(i,4,j)*vpg(1,ng)+
907 . bmerz(i,4,j)*vpg(2,ng)
908 ENDDO
909 ENDDO
910C
911#include "vectorize.inc"
912 DO i=jft,jlt
913C---------------
914C MEMBRANE CONSTANT
915C---------------
916 ep=iplat(i)
917 vdef(ep,1)=vdef(ep,1)+bmrz(i,1,1)*rlz(i,1)+bmrz(i,1,2)*rlz(i,2)
918 1 +bmrz(i,1,3)*rlz(i,3)+bmrz(i,1,4)*rlz(i,4)
919 vdef(ep,2)=vdef(ep,2)+bmrz(i,2,1)*rlz(i,1)+bmrz(i,2,2)*rlz(i,2)
920 1 +bmrz(i,2,3)*rlz(i,3)+bmrz(i,2,4)*rlz(i,4)
921 vdef(ep,3)=vdef(ep,3)+bmrz(i,3,1)*rlz(i,1)+bmrz(i,3,2)*rlz(i,2)
922 1 +bmrz(i,3,3)*rlz(i,3)+bmrz(i,3,4)*rlz(i,4)
923 ENDDO
924#include "vectorize.inc"
925 DO i=jft,nplat
926 ep=iplat(i)
927 bxv2= bm(ep,1)*vxyz(ep,2,1)+bm(ep,2)*vxyz(ep,2,2)
928 1 +bm(ep,3)*vxyz(ep,2,3)
929 byv1=bm(ep,5)*vxyz(ep,1,1)+bm(ep,6)*vxyz(ep,1,2)
930 1 +bm(ep,7)*vxyz(ep,1,3)
931 vrlz(i)=(bxv2-byv1+
932 1 brz(i,4,1)*rlz(i,1)+brz(i,4,2)*rlz(i,2)
933 1 +brz(i,4,3)*rlz(i,3)+brz(i,4,4)*rlz(i,4))*half
934 ENDDO
935C
936#include "vectorize.inc"
937 DO i=nplat+1,jlt
938 ep=iplat(i)
939 erz=zero
940 DO j=1,4
941 erz=erz+brz(i,1,j)*vxyz(ep,1,j)+
942 1 brz(i,2,j)*vxyz(ep,2,j)+brz(i,3,j)*vxyz(ep,3,j)+
943 1 brz(i,4,j)*rlz(i,j)
944 ENDDO
945 vrlz(i)=erz*half
946 ENDDO
947C
948 RETURN

◆ cbadefsh()

subroutine cbadefsh ( integer jft,
integer jlt,
x13,
x24,
y13,
y24,
bm,
vdef,
vxyz,
integer nplat,
integer, dimension(*) iplat )

Definition at line 745 of file cbadef.F.

747C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
748C CALCUL 'membrane shear traitement' use only PARTIE CONSTANTE
749C-----------------------------------------------
750C I M P L I C I T T Y P E S
751C-----------------------------------------------
752#include "implicit_f.inc"
753C-----------------------------------------------
754C G l o b a l P a r a m e t e r s
755C-----------------------------------------------
756#include "mvsiz_p.inc"
757C-----------------------------------------------
758C D U M M Y A R G U M E N T S
759C-----------------------------------------------
760 INTEGER JFT ,JLT,NPLAT,IPLAT(*)
761 my_real
762 . x13(*) ,x24(*) ,y13(*) ,y24(*) ,
763 . bm(mvsiz,36),vdef(mvsiz,8) ,vxyz(mvsiz,3,4)
764C-----------------------------------------------
765C L O C A L V A R I A B L E S
766C-----------------------------------------------
767 INTEGER I ,J,EP
768C
769 DO ep=jft,nplat
770 i=iplat(ep)
771 vdef(i,3)=y24(i)*vxyz(i,2,1)-y13(i)*vxyz(i,2,2)
772 . -x24(i)*vxyz(i,1,1)+x13(i)*vxyz(i,1,2)
773 ENDDO
774 DO ep=nplat+1,jlt
775 i=iplat(ep)
776 bm(i,3)=-x24(i)
777 bm(i,6)=y24(i)
778 bm(i,12)=x13(i)
779 bm(i,15)=-y13(i)
780 bm(i,21)=-bm(i,3)
781 bm(i,24)=-bm(i,6)
782 bm(i,30)=-bm(i,12)
783 bm(i,33)=-bm(i,15)
784 vdef(i,3)=y24(i)*(vxyz(i,2,1)-vxyz(i,2,3))+
785 1 y13(i)*(-vxyz(i,2,2)+vxyz(i,2,4))
786 2 -x24(i)*(vxyz(i,1,1)-vxyz(i,1,3))
787 3 +x13(i)*(vxyz(i,1,2)-vxyz(i,1,4))
788 ENDDO
789 RETURN

◆ cbadeft()

subroutine cbadeft ( integer jft,
integer jlt,
vxyz,
rlz,
vdef,
bm,
integer nplat,
integer, dimension(*) iplat,
integer isrot,
bmrz )

Definition at line 1291 of file cbadef.F.

1293C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1294C-----------------------------------------------
1295C I M P L I C I T T Y P E S
1296C-----------------------------------------------
1297#include "implicit_f.inc"
1298C-----------------------------------------------
1299C G l o b a l P a r a m e t e r s
1300C-----------------------------------------------
1301#include "mvsiz_p.inc"
1302C-----------------------------------------------
1303C D U M M Y A R G U M E N T S
1304C-----------------------------------------------
1305 INTEGER JFT,JLT,NPLAT,IPLAT(*),ISROT
1306 my_real
1307 . rlz(mvsiz,4),vxyz(mvsiz,3,4),
1308 . bm(mvsiz,9,4),vdef(mvsiz,8)
1309C . BM(9,4,*),BMF(9,4,*),BF(6,4,*),VDEF(8,*)
1310 my_real
1311 . bmrz(mvsiz,4,4)
1312C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1313C---VDEF: 1:8 exx eyy,exy,eyx,byx,bxx,byy,bxy
1314C-----------------------------------------------
1315C L O C A L V A R I A B L E S
1316C-----------------------------------------------
1317 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K
1318C
1319 my_real
1320 . a_1,c11,c12,c21,c22,cdemi
1321C--------------------------
1322C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1323#include "vectorize.inc"
1324 DO i=jft,nplat
1325 ep=iplat(i)
1326C-------BM(8):BXI(2),HX(1),not used,BYI(2),HY(1),not used
1327C-------V13(I)->VXYZ(j,1,I),V24(I)->VXYZ(j,2,I),VHI(I)->VXYZ(j,3,I)------------
1328C
1329 vdef(ep,1)=bm(ep,1,1)*vxyz(ep,1,1)+bm(ep,2,1)*vxyz(ep,1,2)
1330 1 +bm(ep,3,1)*vxyz(ep,1,3)
1331 vdef(ep,2)=bm(ep,5,1)*vxyz(ep,2,1)+bm(ep,6,1)*vxyz(ep,2,2)
1332 1 +bm(ep,7,1)*vxyz(ep,2,3)
1333C--------FLEXION-----
1334C VDEF(EP,6)=BM(EP,1,1)*RXYZ(EP,2,1)+BM(EP,2,1)*RXYZ(EP,2,2)
1335C 1 +BM(EP,3,1)*RXYZ(EP,2,3)
1336C VDEF(EP,7)=-(BM(EP,5,1)*RXYZ(EP,1,1)+BM(EP,6,1)*RXYZ(EP,1,2)
1337C 1 +BM(EP,7,1)*RXYZ(EP,1,3))
1338C VDEF(EP,5)=-(BM(EP,1,1)*RXYZ(EP,1,1)+BM(EP,2,1)*RXYZ(EP,1,2)
1339C 1 +BM(EP,3,1)*RXYZ(EP,1,3))
1340C VDEF(EP,8)=BM(EP,5,1)*RXYZ(EP,2,1)+BM(EP,6,1)*RXYZ(EP,2,2)
1341C 3 +BM(EP,7,1)*RXYZ(EP,2,3)
1342 ENDDO
1343c constant shear is removed
1344#include "vectorize.inc"
1345 DO i=jft,nplat
1346 ep=iplat(i)
1347C---------eyx-----------
1348 vdef(ep,4)=bm(ep,1,1)*vxyz(ep,2,1)+bm(ep,2,1)*vxyz(ep,2,2)
1349 1 +bm(ep,3,1)*vxyz(ep,2,3)
1350C-------- exy---
1351 vdef(ep,3)=bm(ep,5,1)*vxyz(ep,1,1)+bm(ep,6,1)*vxyz(ep,1,2)
1352 1 +bm(ep,7,1)*vxyz(ep,1,3)
1353 END DO
1354c
1355C------------ELEMENT GAUCH-------------------
1356#include "vectorize.inc"
1357 DO ep=nplat+1,jlt
1358 i=iplat(ep)
1359C I = 1; JJ = 3*KK =0 terms VXYZ(3,j,I) shared
1360 cdemi = half*(bm(i,9,1)*vxyz(i,3,1)+bm(i,9,2)*vxyz(i,3,2)+
1361 + bm(i,9,3)*vxyz(i,3,3)+bm(i,9,4)*vxyz(i,3,4))
1362 vdef(i,4)=bm(i,6,1)*vxyz(i,2,1)+bm(i,6,2)*vxyz(i,2,2)+
1363 1 bm(i,6,3)*vxyz(i,2,3)+bm(i,6,4)*vxyz(i,2,4)+cdemi
1364 vdef(i,3)=bm(i,3,1)*vxyz(i,1,1)+bm(i,3,2)*vxyz(i,1,2)+
1365 1 bm(i,3,3)*vxyz(i,1,3)+bm(i,3,4)*vxyz(i,1,4)+cdemi
1366C---------DEF MEMBRANE------
1367 vdef(i,1)=
1368 +bm(i,1,1)*vxyz(i,1,1)+bm(i,4,1)*vxyz(i,2,1)+bm(i,7,1)*vxyz(i,3,1)+
1369 +bm(i,1,2)*vxyz(i,1,2)+bm(i,4,2)*vxyz(i,2,2)+bm(i,7,2)*vxyz(i,3,2)+
1370 +bm(i,1,3)*vxyz(i,1,3)+bm(i,4,3)*vxyz(i,2,3)+bm(i,7,3)*vxyz(i,3,3)+
1371 +bm(i,1,4)*vxyz(i,1,4)+bm(i,4,4)*vxyz(i,2,4)+bm(i,7,4)*vxyz(i,3,4)
1372 vdef(i,2)=
1373 +bm(i,2,1)*vxyz(i,1,1)+bm(i,5,1)*vxyz(i,2,1)+bm(i,8,1)*vxyz(i,3,1)+
1374 +bm(i,2,2)*vxyz(i,1,2)+bm(i,5,2)*vxyz(i,2,2)+bm(i,8,2)*vxyz(i,3,2)+
1375 +bm(i,2,3)*vxyz(i,1,3)+bm(i,5,3)*vxyz(i,2,3)+bm(i,8,3)*vxyz(i,3,3)+
1376 +bm(i,2,4)*vxyz(i,1,4)+bm(i,5,4)*vxyz(i,2,4)+bm(i,8,4)*vxyz(i,3,4)
1377C
1378C---------FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
1379C----------[B01]---------
1380C VDEF(I,6)= BMF(1,1,I)*VXYZ(1,1,I)+
1381C + BMF(4,1,I)*VXYZ(2,1,I)+BMF(7,1,I)*VXYZ(3,1,I)+
1382C + BMF(1,2,I)*VXYZ(1,2,I)+
1383C + BMF(4,2,I)*VXYZ(2,2,I)+BMF(7,2,I)*VXYZ(3,2,I)+
1384C + BMF(1,3,I)*VXYZ(1,3,I)+
1385C + BMF(4,3,I)*VXYZ(2,3,I)+BMF(7,3,I)*VXYZ(3,3,I)+
1386C + BMF(1,4,I)*VXYZ(1,4,I)+
1387C + BMF(4,4,I)*VXYZ(2,4,I)+BMF(7,4,I)*VXYZ(3,4,I)
1388C VDEF(I,7)= BMF(2,1,I)*VXYZ(1,1,I)+
1389C + BMF(5,1,I)*VXYZ(2,1,I)+BMF(8,1,I)*VXYZ(3,1,I)+
1390C + BMF(2,2,I)*VXYZ(1,2,I)+
1391C + BMF(5,2,I)*VXYZ(2,2,I)+BMF(8,2,I)*VXYZ(3,2,I)+
1392C + BMF(2,3,I)*VXYZ(1,3,I)+
1393C + BMF(5,3,I)*VXYZ(2,3,I)+BMF(8,3,I)*VXYZ(3,3,I)+
1394C + BMF(2,4,I)*VXYZ(1,4,I)+
1395C + BMF(5,4,I)*VXYZ(2,4,I)+BMF(8,4,I)*VXYZ(3,4,I)
1396C VXYZ(3,j,I) shared
1397C CDEMI = HALF*(BMF(9,1,I)*VXYZ(3,1,I)+BMF(9,2,I)*VXYZ(3,2,I)+
1398C + BMF(9,3,I)*VXYZ(3,3,I)+BMF(9,4,I)*VXYZ(3,4,I))
1399C VDEF(I,8)= BMF(3,1,I)*VXYZ(1,1,I)+BMF(3,2,I)*VXYZ(1,2,I)+
1400C + BMF(3,3,I)*VXYZ(1,3,I)+BMF(3,4,I)*VXYZ(1,4,I)+CDEMI
1401C VDEF(I,5)= BMF(6,1,I)*VXYZ(2,1,I)+BMF(6,2,I)*VXYZ(2,2,I)+
1402C + BMF(6,3,I)*VXYZ(2,3,I)+BMF(6,4,I)*VXYZ(2,4,I)+CDEMI
1403C---------- TERMES ASSOCIEES AU BETA ----------
1404C----------TERMES ASSOCIEES AUX BETA DE [B01]--
1405C JJ = 3*II =0
1406C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
1407C VDEF(I,6)= VDEF(I,6)+
1408C + BF(1,1,I)*RXYZ(1,1,I)+BF(4,1,I)*RXYZ(2,1,I)+
1409C + BF(1,2,I)*RXYZ(1,2,I)+BF(4,2,I)*RXYZ(2,2,I)+
1410C + BF(1,3,I)*RXYZ(1,3,I)+BF(4,3,I)*RXYZ(2,3,I)+
1411C + BF(1,4,I)*RXYZ(1,4,I)+BF(4,4,I)*RXYZ(2,4,I)
1412C VDEF(I,7)= VDEF(I,7)+
1413C + BF(2,1,I)*RXYZ(1,1,I)+BF(5,1,I)*RXYZ(2,1,I)+
1414C + BF(2,2,I)*RXYZ(1,2,I)+BF(5,2,I)*RXYZ(2,2,I)+
1415C + BF(2,3,I)*RXYZ(1,3,I)+BF(5,3,I)*RXYZ(2,3,I)+
1416C + BF(2,4,I)*RXYZ(1,4,I)+BF(5,4,I)*RXYZ(2,4,I)
1417C VDEF(I,8)= VDEF(I,5+
1418C + BF(6,1,I)*RXYZ(2,1,I)+BF(6,2,I)*RXYZ(2,2,I)+
1419C + BF(6,3,I)*RXYZ(2,3,I)+BF(6,4,I)*RXYZ(2,4,I)
1420C VDEF(I,5)= VDEF(I,8+
1421C + BF(3,1,I)*RXYZ(1,1,I)+BF(3,2,I)*RXYZ(1,2,I)+
1422C + BF(3,3,I)*RXYZ(1,3,I)+BF(3,4,I)*RXYZ(1,4,I)
1423C
1424 END DO
1425 IF (isrot>0) THEN
1426#include "vectorize.inc"
1427 DO i=jft,jlt
1428C---------------
1429C MEMBRANE CONSTANT
1430C---------------
1431 ep=iplat(i)
1432 vdef(ep,1)=vdef(ep,1)+bmrz(i,1,1)*rlz(ep,1)+bmrz(i,1,2)*
1433 1 rlz(ep,2)+bmrz(i,1,3)*rlz(ep,3)+bmrz(i,1,4)*rlz(ep,4)
1434 vdef(ep,2)=vdef(ep,2)+bmrz(i,2,1)*rlz(ep,1)+bmrz(i,2,2)*
1435 1 rlz(ep,2)+bmrz(i,2,3)*rlz(ep,3)+bmrz(i,2,4)*rlz(ep,4)
1436 vdef(ep,3)=vdef(ep,3)+
1437 1 bmrz(i,3,1)*rlz(ep,1)+bmrz(i,3,2)*rlz(ep,2)
1438 2 +bmrz(i,3,3)*rlz(ep,3)+bmrz(i,3,4)*rlz(ep,4)
1439 vdef(ep,4)=vdef(ep,4)+
1440 1 bmrz(i,4,1)*rlz(ep,1)+bmrz(i,4,2)*rlz(ep,2)
1441 2 +bmrz(i,4,3)*rlz(ep,3)+bmrz(i,4,4)*rlz(ep,4)
1442 ENDDO
1443 END IF !(ISROT>0) THEN
1444C
1445 RETURN

◆ cbadeft1()

subroutine cbadeft1 ( integer jft,
integer jlt,
integer ng,
vcore,
vxyz,
vdef,
hx,
hy,
bm,
integer nplat,
integer, dimension(*) iplat,
integer isrot,
bmrz,
rxyz,
wxy )

Definition at line 1452 of file cbadef.F.

1455C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1456C-----------------------------------------------
1457C I M P L I C I T T Y P E S
1458C-----------------------------------------------
1459#include "implicit_f.inc"
1460C-----------------------------------------------
1461C G l o b a l P a r a m e t e r s
1462C-----------------------------------------------
1463#include "mvsiz_p.inc"
1464C-----------------------------------------------
1465C D U M M Y A R G U M E N T S
1466C-----------------------------------------------
1467 INTEGER NG,JFT,JLT,NPLAT,IPLAT(*),ISROT
1468 my_real
1469 . vcore(mvsiz,12),vxyz(mvsiz,3,4),bm(mvsiz,9,4),hx(mvsiz,4),hy(mvsiz,4)
1470 my_real
1471 . vdef(mvsiz,8),rxyz(mvsiz,4),bmrz(mvsiz,4,4),wxy(*)
1472C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1473C---VDEF: 1:8 exx eyy,exy,eyx,byx,bxx,byy,bxy
1474C-----------------------------------------------
1475C L O C A L V A R I A B L E S
1476C-----------------------------------------------
1477 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K
1478C
1479 my_real
1480 . a_1,c11,c12,c21,c22,cdemi
1481C--------------------------
1482C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1483 DO ep=jft,jlt
1484C--------FLEXION-----
1485 vdef(ep,6)=zero
1486 vdef(ep,7)=zero
1487 vdef(ep,8)=zero
1488 vdef(ep,5)=zero
1489 wxy(ep)=zero
1490 ENDDO
1491#include "vectorize.inc"
1492 DO i=jft,nplat
1493 ep=iplat(i)
1494C-------BM(8):BXI(2),HX(1),not used,BYI(2),HY(1),not used
1495 bm(ep,1,1)=vcore(ep,1)+hx(ep,ng)*vcore(ep,5)
1496 bm(ep,2,1)=vcore(ep,2)+hx(ep,ng)*vcore(ep,6)
1497 bm(ep,3,1)=hx(ep,ng)*fourth
1498 bm(ep,4,1)=-bm(ep,3,1)
1499 bm(ep,5,1)=vcore(ep,3)+hy(ep,ng)*vcore(ep,5)
1500 bm(ep,6,1)=vcore(ep,4)+hy(ep,ng)*vcore(ep,6)
1501 bm(ep,7,1)=hy(ep,ng)*fourth
1502 bm(ep,8,1)=-bm(ep,7,1)
1503C-------V13(I)->VXYZ(j,1,I),V24(I)->VXYZ(j,2,I),VHI(I)->VXYZ(j,3,I)------------
1504C
1505 vdef(ep,1)=bm(ep,1,1)*vxyz(ep,1,1)+bm(ep,2,1)*vxyz(ep,1,2)
1506 1 +bm(ep,3,1)*vxyz(ep,1,3)
1507 vdef(ep,2)=bm(ep,5,1)*vxyz(ep,2,1)+bm(ep,6,1)*vxyz(ep,2,2)
1508 1 +bm(ep,7,1)*vxyz(ep,2,3)
1509C---------eyx-----------
1510 vdef(ep,4)=bm(ep,1,1)*vxyz(ep,2,1)+bm(ep,2,1)*vxyz(ep,2,2)
1511 1 +bm(ep,3,1)*vxyz(ep,2,3)
1512C-------- exy---
1513 vdef(ep,3)=bm(ep,5,1)*vxyz(ep,1,1)+bm(ep,6,1)*vxyz(ep,1,2)
1514 1 +bm(ep,7,1)*vxyz(ep,1,3)
1515 ENDDO
1516 IF (isrot>0) THEN
1517#include "vectorize.inc"
1518 DO i=jft,nplat
1519C---------------
1520C MEMBRANE CONSTANT
1521C---------------
1522 ep=iplat(i)
1523 vdef(ep,1)=vdef(ep,1)+bmrz(i,1,1)*rxyz(ep,1)+bmrz(i,1,2)*
1524 1 rxyz(ep,2)+bmrz(i,1,3)*rxyz(ep,3)+bmrz(i,1,4)*rxyz(ep,4)
1525 vdef(ep,2)=vdef(ep,2)+bmrz(i,2,1)*rxyz(ep,1)+bmrz(i,2,2)*
1526 1 rxyz(ep,2)+bmrz(i,2,3)*rxyz(ep,3)+bmrz(i,2,4)*rxyz(ep,4)
1527 vdef(ep,3)=vdef(ep,3)+
1528 1 bmrz(i,3,1)*rxyz(ep,1)+bmrz(i,3,2)*rxyz(ep,2)
1529 2 +bmrz(i,3,3)*rxyz(ep,3)+bmrz(i,3,4)*rxyz(ep,4)
1530 vdef(ep,4)=vdef(ep,4)+
1531 1 bmrz(i,4,1)*rxyz(ep,1)+bmrz(i,4,2)*rxyz(ep,2)
1532 2 +bmrz(i,4,3)*rxyz(ep,3)+bmrz(i,4,4)*rxyz(ep,4)
1533 ENDDO
1534 END IF !(ISROT>0) THEN
1535C
1536 RETURN

◆ cbadeftw()

subroutine cbadeftw ( integer jft,
integer jlt,
vxyz,
rxyz,
bm,
bmf,
bf,
integer nplat,
integer, dimension(*) iplat,
wxy )

Definition at line 1784 of file cbadef.F.

1787C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1788C-----------------------------------------------
1789C I M P L I C I T T Y P E S
1790C-----------------------------------------------
1791#include "implicit_f.inc"
1792C-----------------------------------------------
1793C G l o b a l P a r a m e t e r s
1794C-----------------------------------------------
1795#include "mvsiz_p.inc"
1796C-----------------------------------------------
1797C D U M M Y A R G U M E N T S
1798C-----------------------------------------------
1799 INTEGER JFT,JLT,NPLAT,IPLAT(*)
1800 my_real
1801 . rxyz(mvsiz,2,4),vxyz(mvsiz,3,4),
1802 . bm(mvsiz,9,4),bmf(mvsiz,9,4),bf(mvsiz,6,4),wxy(*)
1803C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1804C---VDEF: 1:8 exx eyy,exy,eyx,byx,bxx,byy,bxy
1805C-----------------------------------------------
1806C L O C A L V A R I A B L E S
1807C-----------------------------------------------
1808 INTEGER I,J,EP,K
1809C
1810 my_real
1811 . a_1,c11,c12,c21,c22,cdemi,kxy,kyx
1812C--------------------------
1813C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1814#include "vectorize.inc"
1815 DO i=jft,nplat
1816 ep=iplat(i)
1817 kyx =-(bm(ep,1,1)*rxyz(ep,1,1)+bm(ep,2,1)*rxyz(ep,1,2)
1818 1 +bm(ep,3,1)*rxyz(ep,1,3))
1819 kxy =bm(ep,5,1)*rxyz(ep,2,1)+bm(ep,6,1)*rxyz(ep,2,2)
1820 3 +bm(ep,7,1)*rxyz(ep,2,3)
1821 wxy(ep) = kxy - kyx
1822 ENDDO
1823C------------ELEMENT GAUCH-------------------
1824#include "vectorize.inc"
1825 DO i=nplat+1,jlt
1826 ep=iplat(i)
1827C
1828 cdemi =half*(bmf(ep,9,1)*vxyz(ep,3,1)+bmf(i,9,2)*vxyz(ep,3,2)+
1829 + bmf(ep,9,3)*vxyz(ep,3,3)+bmf(i,9,4)*vxyz(ep,3,4))
1830 kxy = bmf(ep,3,1)*vxyz(ep,1,1)+bmf(ep,3,2)*vxyz(ep,1,2)+
1831 + bmf(ep,3,3)*vxyz(ep,1,3)+bmf(ep,3,4)*vxyz(ep,1,4)+cdemi
1832 kyx = bmf(ep,6,1)*vxyz(ep,2,1)+bmf(ep,6,2)*vxyz(ep,2,2)+
1833 + bmf(ep,6,3)*vxyz(ep,2,3)+bmf(ep,6,4)*vxyz(ep,2,4)+cdemi
1834 kxy = kxy +
1835 + bf(ep,6,1)*rxyz(ep,2,1)+bf(ep,6,2)*rxyz(ep,2,2)+
1836 + bf(ep,6,3)*rxyz(ep,2,3)+bf(ep,6,4)*rxyz(ep,2,4)
1837 kyx = kyx +
1838 + bf(ep,3,1)*rxyz(ep,1,1)+bf(ep,3,2)*rxyz(ep,1,2)+
1839 + bf(ep,3,3)*rxyz(ep,1,3)+bf(ep,3,4)*rxyz(ep,1,4)
1840C
1841 wxy(ep) = kxy - kyx
1842 END DO
1843C
1844 RETURN

◆ cbaderirz()

subroutine cbaderirz ( integer jft,
integer jlt,
area,
x13,
x24,
y13,
y24,
bm0rz,
bmkrz,
bmerz,
vcore,
integer nplat,
integer, dimension(*) iplat,
integer ismstr )

Definition at line 799 of file cbadef.F.

802C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
803#include "implicit_f.inc"
804#include "mvsiz_p.inc"
805C-----------------------------------------------
806C D U M M Y A R G U M E N T S
807C-----------------------------------------------
808 INTEGER JFT,JLT,NPLAT,IPLAT(*),ISMSTR
809 my_real
810 . vcore(mvsiz,3,4),x13(*),x24(*),y13(*),y24(*),
811 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
812 . area(*)
813C-----------------------------------------------
814C L O C A L V A R I A B L E S
815C-----------------------------------------------
816 INTEGER I,J,EP
817 my_real
818 . mx13(mvsiz),mx23(mvsiz),mx34(mvsiz),my13(mvsiz),
819 . my23(mvsiz),my34(mvsiz),a1(mvsiz),
820 . x13l(mvsiz),x24l(mvsiz),y13l(mvsiz),y24l(mvsiz)
821C-----------------------------------------------
822#include "vectorize.inc"
823C-------VCORE(12):BX0(2),BY0(2),GAMA(2),MX23,MY23,MX34,MY34,MX13,MY13
824 DO i=jft,nplat
825 ep =iplat(i)
826 mx23(i) = vcore(ep,1,3)
827 my23(i) = vcore(ep,2,3)
828 mx34(i) = vcore(ep,3,3)
829 my34(i) = vcore(ep,1,4)
830 mx13(i) = vcore(ep,2,4)
831 my13(i) = vcore(ep,3,4)
832 x13l(i) = x13(ep)
833 x24l(i) = x24(ep)
834 y13l(i) = y13(ep)
835 y24l(i) = y24(ep)
836 a1(i) = area(ep)
837 END DO !I=JFT,NPLAT
838C
839 DO i=nplat+1,jlt
840 ep =iplat(i)
841 mx13(i) = (vcore(ep,1,1)+vcore(ep,1,3))*half
842 my13(i) = (vcore(ep,2,1)+vcore(ep,2,3))*half
843 mx23(i) = (vcore(ep,1,2)+vcore(ep,1,3))*half
844 my23(i) = (vcore(ep,2,2)+vcore(ep,2,3))*half
845 mx34(i) = (vcore(ep,1,3)+vcore(ep,1,4))*half
846 my34(i) = (vcore(ep,2,3)+vcore(ep,2,4))*half
847 a1(i) = area(ep)
848 x13l(i) = x13(ep)
849 x24l(i) = x24(ep)
850 y13l(i) = y13(ep)
851 y24l(i) = y24(ep)
852 END DO !I=NPLAT+1,JLT
853C
854 CALL cbaderirz0(jft ,jlt ,a1 ,x13l ,x24l ,
855 2 y13l ,y24l ,mx13 ,mx23 ,mx34 ,
856 3 my13 ,my23 ,my34 ,bm0rz,bmkrz ,
857 4 bmerz,ismstr)
858C
859 RETURN
cbaderirz0
subroutine cbaderirz0(jft, jlt, area, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, bm0rz, bmkrz, bmerz, ismstr)
Definition cbadef.F:961

◆ cbaderirz0()

subroutine cbaderirz0 ( integer jft,
integer jlt,
area,
x13,
x24,
y13,
y24,
mx13,
mx23,
mx34,
my13,
my23,
my34,
bm0rz,
bmkrz,
bmerz,
integer ismstr )

Definition at line 957 of file cbadef.F.

961C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
962#include "implicit_f.inc"
963#include "mvsiz_p.inc"
964C-----------------------------------------------
965C D U M M Y A R G U M E N T S
966C-----------------------------------------------
967 INTEGER JFT,JLT,ISMSTR
968 my_real
969 . x13(*),x24(*),y13(*),y24(*),
970 . mx13(*),mx23(*),mx34(*),my13(*),my23(*),my34(*)
971 my_real
972 . area(*),
973 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4)
974C-----------------------------------------------
975C L O C A L V A R I A B L E S
976C-----------------------------------------------
977 INTEGER I,J
978 my_real
979 . aa4,bxv2,byv1,a05,nxy,nyx,a05k(4),a05e(4),
980 . xij(mvsiz,4),yij(mvsiz,4),xil(mvsiz,4),yil(mvsiz,4)
981C--------(14,23,32,41)/4A----------X13:X13/A-----------------------------
982 DO i=jft,jlt
983 a05=half/area(i)
984 xij(i,1) = a05*(mx13(i)-mx34(i))
985 xij(i,2) = half*x13(i)+a05*(mx23(i)-mx13(i))
986 xij(i,3) = -xij(i,2)
987 xij(i,4) = -xij(i,1)
988 yij(i,1) = a05*(my13(i)-my34(i))
989 yij(i,2) = half*y13(i)+a05*(my23(i)-my13(i))
990 yij(i,3) = -yij(i,2)
991 yij(i,4) = -yij(i,1)
992C--------(12,21,34,43)/4---------------------------------------
993 xil(i,1) = a05*(mx13(i)-mx23(i))
994 xil(i,2) = -xil(i,1)
995 xil(i,3) = a05*(mx13(i)-mx34(i))-half*x13(i)
996 xil(i,4) = -xil(i,3)
997 yil(i,1) = a05*(my13(i)-my23(i))
998 yil(i,2) = -yil(i,1)
999 yil(i,3) = a05*(my13(i)-my34(i))-half*y13(i)
1000 yil(i,4) = -yil(i,3)
1001 ENDDO
1002C
1003 DO i=jft,jlt
1004C--------NxI,x ------
1005C-------- MY34YIJ KSI - MY23YIL ETA
1006 bm0rz(i,1,1) = (-my34(i)*yij(i,1)+my23(i)*yil(i,1))
1007 bm0rz(i,1,2) = ( my34(i)*yij(i,2)+my23(i)*yil(i,2))
1008 bm0rz(i,1,3) = ( my34(i)*yij(i,3)-my23(i)*yil(i,3))
1009 bm0rz(i,1,4) = (-my34(i)*yij(i,4)-my23(i)*yil(i,4))
1010C-------- MY13YIJ KSI - 2MY34YIL
1011 bmkrz(i,1,1) = (-my13(i)*yij(i,1)-two*my34(i)*yil(i,1))
1012 bmkrz(i,1,2) = ( my13(i)*yij(i,2)-two*my34(i)*yil(i,2))
1013 bmkrz(i,1,3) = ( my13(i)*yij(i,3)-two*my34(i)*yil(i,3))
1014 bmkrz(i,1,4) = (-my13(i)*yij(i,4)-two*my34(i)*yil(i,4))
1015C-------- -MY13YIL ETA + 2MY23YIJ
1016 bmerz(i,1,1) = ( my13(i)*yil(i,1)+two*my23(i)*yij(i,1))
1017 bmerz(i,1,2) = ( my13(i)*yil(i,2)+two*my23(i)*yij(i,2))
1018 bmerz(i,1,3) = (-my13(i)*yil(i,3)+two*my23(i)*yij(i,3))
1019 bmerz(i,1,4) = (-my13(i)*yil(i,4)+two*my23(i)*yij(i,4))
1020C -------NyI,y -------------
1021C-------- MX34XIJ KSI - MX23XIL ETA
1022 bm0rz(i,2,1) = (-mx34(i)*xij(i,1)+mx23(i)*xil(i,1))
1023 bm0rz(i,2,2) = ( mx34(i)*xij(i,2)+mx23(i)*xil(i,2))
1024 bm0rz(i,2,3) = ( mx34(i)*xij(i,3)-mx23(i)*xil(i,3))
1025 bm0rz(i,2,4) = (-mx34(i)*xij(i,4)-mx23(i)*xil(i,4))
1026C-------- MX13XIJ KSI - 2MX34XIL
1027 bmkrz(i,2,1) = (-mx13(i)*xij(i,1)-two*mx34(i)*xil(i,1))
1028 bmkrz(i,2,2) = ( mx13(i)*xij(i,2)-two*mx34(i)*xil(i,2))
1029 bmkrz(i,2,3) = ( mx13(i)*xij(i,3)-two*mx34(i)*xil(i,3))
1030 bmkrz(i,2,4) = (-mx13(i)*xij(i,4)-two*mx34(i)*xil(i,4))
1031C-------- -MX13XIL ETA + 2MX23XIJ
1032 bmerz(i,2,1) = ( mx13(i)*xil(i,1)+two*mx23(i)*xij(i,1))
1033 bmerz(i,2,2) = ( mx13(i)*xil(i,2)+two*mx23(i)*xij(i,2))
1034 bmerz(i,2,3) = (-mx13(i)*xil(i,3)+two*mx23(i)*xij(i,3))
1035 bmerz(i,2,4) = (-mx13(i)*xil(i,4)+two*mx23(i)*xij(i,4))
1036C--------NxI,y
1037C-------- -MX34YIJ KSI + MX23YIL ETA
1038 bm0rz(i,3,1) = ( mx34(i)*yij(i,1)-mx23(i)*yil(i,1))
1039 bm0rz(i,3,2) = (-mx34(i)*yij(i,2)-mx23(i)*yil(i,2))
1040 bm0rz(i,3,3) = (-mx34(i)*yij(i,3)+mx23(i)*yil(i,3))
1041 bm0rz(i,3,4) = ( mx34(i)*yij(i,4)+mx23(i)*yil(i,4))
1042C-------- -MX13YIJ KSI + 2MX34YIL
1043 bmkrz(i,3,1) = ( mx13(i)*yij(i,1)+two*mx34(i)*yil(i,1))
1044 bmkrz(i,3,2) = (-mx13(i)*yij(i,2)+two*mx34(i)*yil(i,2))
1045 bmkrz(i,3,3) = (-mx13(i)*yij(i,3)+two*mx34(i)*yil(i,3))
1046 bmkrz(i,3,4) = ( mx13(i)*yij(i,4)+two*mx34(i)*yil(i,4))
1047C-------- MX13YIL ETA - 2MX23YIJ
1048 bmerz(i,3,1) = (-mx13(i)*yil(i,1)-two*mx23(i)*yij(i,1))
1049 bmerz(i,3,2) = (-mx13(i)*yil(i,2)-two*mx23(i)*yij(i,2))
1050 bmerz(i,3,3) = ( mx13(i)*yil(i,3)-two*mx23(i)*yij(i,3))
1051 bmerz(i,3,4) = ( mx13(i)*yil(i,4)-two*mx23(i)*yij(i,4))
1052C--------NyI,x
1053C-------- -MY34XIJ KSI + MY23XIJ ETA
1054 bm0rz(i,4,1) = ( my34(i)*xij(i,1)-my23(i)*xil(i,1))
1055 bm0rz(i,4,2) = (-my34(i)*xij(i,2)-my23(i)*xil(i,2))
1056 bm0rz(i,4,3) = (-my34(i)*xij(i,3)+my23(i)*xil(i,3))
1057 bm0rz(i,4,4) = ( my34(i)*xij(i,4)+my23(i)*xil(i,4))
1058C-------- -MY13XIJ KSI + 2MY34XIL
1059 bmkrz(i,4,1) = ( my13(i)*xij(i,1)+two*my34(i)*xil(i,1))
1060 bmkrz(i,4,2) = (-my13(i)*xij(i,2)+two*my34(i)*xil(i,2))
1061 bmkrz(i,4,3) = (-my13(i)*xij(i,3)+two*my34(i)*xil(i,3))
1062 bmkrz(i,4,4) = ( my13(i)*xij(i,4)+two*my34(i)*xil(i,4))
1063C-------- MY13XIL ETA - 2MY23XIJ
1064 bmerz(i,4,1) = (-my13(i)*xil(i,1)-two*my23(i)*xij(i,1))
1065 bmerz(i,4,2) = (-my13(i)*xil(i,2)-two*my23(i)*xij(i,2))
1066 bmerz(i,4,3) = ( my13(i)*xil(i,3)-two*my23(i)*xij(i,3))
1067 bmerz(i,4,4) = ( my13(i)*xil(i,4)-two*my23(i)*xij(i,4))
1068 ENDDO
1069C--------correction to pass the patch test
1070 CALL cbpatch(jft ,jlt ,bm0rz,bmkrz,bmerz ,area ,
1071 + mx13 ,mx23,mx34 ,my13 ,my23 ,my34)
1072C
1073 a05= half
1074 IF (ismstr==10) a05=zero
1075 DO i=jft,jlt
1076C--------BM0RZ(I,3,1)->(NxI,y+NyI,x) -;BM0RZ(I,4,1)->(-NxI,y+NyI,x -2NI)-----
1077 a05k(1)= -a05
1078 a05k(2)= a05
1079 a05k(3)= a05
1080 a05k(4)= -a05
1081 a05e(1)= -a05
1082 a05e(2)= -a05
1083 a05e(3)= a05
1084 a05e(4)= a05
1085 DO j=1,4
1086 nxy=bm0rz(i,3,j)
1087 nyx=bm0rz(i,4,j)
1088 bm0rz(i,3,j)=nxy+nyx
1089 bm0rz(i,4,j)=-nxy+nyx-a05
1090 nxy=bmkrz(i,3,j)
1091 nyx=bmkrz(i,4,j)
1092 bmkrz(i,3,j)=nxy+nyx
1093 bmkrz(i,4,j)=-nxy+nyx-a05k(j)
1094 nxy=bmerz(i,3,j)
1095 nyx=bmerz(i,4,j)
1096 bmerz(i,3,j)=nxy+nyx
1097 bmerz(i,4,j)=-nxy+nyx-a05e(j)
1098 ENDDO
1099 ENDDO
1100C
1101 RETURN
cbpatch
subroutine cbpatch(jft, jlt, bm0rz, bmkrz, bmerz, area, mx13, mx23, mx34, my13, my23, my34)
Definition cbadef.F:1112

◆ cbaderirzt()

subroutine cbaderirzt ( integer jft,
integer jlt,
integer ng,
bm0rz,
bmkrz,
bmerz,
bmrz )

Definition at line 1730 of file cbadef.F.

1731C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1732C CALCUL DES DEFORMATIONS GENERALISEES AUX POINTS DE GAUSS ET [B]
1733C ENTREES : JFT,JLT,NG,BM0RZ,BMKRZ,BMERZ
1734C SORTIES : BMRZ
1735C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1736C-----------------------------------------------
1737C I M P L I C I T T Y P E S
1738C-----------------------------------------------
1739#include "implicit_f.inc"
1740#include "mvsiz_p.inc"
1741C-----------------------------------------------
1742C D U M M Y A R G U M E N T S
1743C-----------------------------------------------
1744 INTEGER NG,JFT,JLT
1745 my_real
1746 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
1747 . bmrz(mvsiz,4,4)
1748C-----------------------------------------------
1749C L O C A L V A R I A B L E S
1750C-----------------------------------------------
1751 INTEGER NPG,I,J
1752 parameter(npg = 4)
1753 my_real
1754 . vpg(2,npg),pg1,pg,nxyp,nyxn
1755 parameter(pg=.577350269189626)
1756 parameter(pg1=-.577350269189626)
1757C--------------------------
1758C INITIALISATION
1759C--------------------------
1760 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
1761C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1762C---+BMRZ(3,J,I): NXY ,BMRZ(4,J,I):NYX
1763 DO i=jft,jlt
1764 DO j=1,4
1765 bmrz(i,1,j) = bm0rz(i,1,j)+bmkrz(i,1,j)*vpg(1,ng)+
1766 . bmerz(i,1,j)*vpg(2,ng)
1767 bmrz(i,2,j) = bm0rz(i,2,j)+bmkrz(i,2,j)*vpg(1,ng)+
1768 . bmerz(i,2,j)*vpg(2,ng)
1769 nxyp = bm0rz(i,3,j)+bmkrz(i,3,j)*vpg(1,ng)+
1770 . bmerz(i,3,j)*vpg(2,ng)
1771 nyxn = bm0rz(i,4,j)+bmkrz(i,4,j)*vpg(1,ng)+
1772 . bmerz(i,4,j)*vpg(2,ng)
1773 bmrz(i,4,j) = half*(nxyp + nyxn)
1774 bmrz(i,3,j) = half*(nxyp - nyxn)
1775 ENDDO
1776 ENDDO
1777 RETURN

◆ cbaderit1()

subroutine cbaderit1 ( integer jft,
integer jlt,
integer ng,
vcore,
vq,
vjfi,
hx,
hy,
veta,
vksi,
bm,
integer nplat,
integer, dimension(*) iplat,
integer isrot )

Definition at line 1543 of file cbadef.F.

1546C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1547C CALCUL DES DEFORMATIONS GENERALISEES AUX POINTS DE GAUSS ET [B]
1548C ENTREES : NEL,NG,VCORE,VQN,VXYZ,RXYZ,VNRM,VASM,VASTN
1549C SORTIES : VDEF,BM (MEMBRANE),BMF(COUPLAGE MEM-FLEXION),BF (FLEXION),BC(C.T)
1550C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1551C-----------------------------------------------
1552C I M P L I C I T T Y P E S
1553C-----------------------------------------------
1554#include "implicit_f.inc"
1555#include "mvsiz_p.inc"
1556C-----------------------------------------------
1557C D U M M Y A R G U M E N T S
1558C-----------------------------------------------
1559 INTEGER NPG,NG,JFT,JLT,NPLAT,IPLAT(*),ISROT
1560 parameter(npg = 4)
1561 my_real
1562 . vcore(mvsiz,12),vksi(4,4),veta(4,4),
1563 . bm(mvsiz,36),hx(mvsiz,4),hy(mvsiz,4)
1564 my_real
1565 . vjfi(mvsiz,3,2,4),vq(mvsiz,3,3,4)
1566C-----------------------------------------------
1567C L O C A L V A R I A B L E S
1568C-----------------------------------------------
1569 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K
1570 my_real
1571 . vpg(2,npg),pg1,pg,thk,detj,det,
1572 . tfn(3,2),bcx,bcy,bxy(3),byx(3),
1573 . v1(2),v2(2),rv1,rv2,
1574 . c1,c2,vt1,vt2,bc1,bc2,vb1,tc(mvsiz,2,2),
1575 . vjf1(2,3),vjf(3,3),tbi(2,2),tbc(2,2),v11(4)
1576 my_real
1577 . a_1,c11,c12,c21,c22,cc,beta1,ksi1,ksiy1,beta2,ksi2,ksiy2
1578 parameter(pg=.577350269189626)
1579 parameter(pg1=-.577350269189626)
1580C--------------------------
1581C INITIALISATION
1582C--------------------------
1583 DATA iun/1/
1584 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
1585C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1586#include "vectorize.inc"
1587 DO i=jft,nplat
1588 ep=iplat(i)
1589C-------VCORE(12):BX0(2),BY0(2),GAMA(2),MX23,MY23,MX34,MY34,MX13,MY13
1590 bm(ep,1)=vcore(ep,1)+hx(ep,ng)*vcore(ep,5)
1591 bm(ep,2)=vcore(ep,2)+hx(ep,ng)*vcore(ep,6)
1592 bm(ep,3)=hx(ep,ng)*fourth
1593 bm(ep,4)=-bm(ep,3)
1594 bm(ep,5)=vcore(ep,3)+hy(ep,ng)*vcore(ep,5)
1595 bm(ep,6)=vcore(ep,4)+hy(ep,ng)*vcore(ep,6)
1596 bm(ep,7)=hy(ep,ng)*fourth
1597 bm(ep,8)=-bm(ep,7)
1598 ENDDO
1599C------------ELEMENT GAUCH-------------------
1600#include "vectorize.inc"
1601 DO 150 i=nplat+1,jlt
1602 ep=iplat(i)
1603C---------------------------
1604C CALCUL DE TC=VJFI*VQ
1605C---------------------------
1606 tc(ep,1,1)=vjfi(ep,1,1,ng)*vq(ep,1,1,ng)+vjfi(ep,2,1,ng)
1607 + *vq(ep,2,1,ng)+ vjfi(ep,3,1,ng)*vq(ep,3,1,ng)
1608 tc(ep,2,1)=vjfi(ep,1,2,ng)*vq(ep,1,1,ng)+vjfi(ep,2,2,ng)
1609 + *vq(ep,2,1,ng)+ vjfi(ep,3,2,ng)*vq(ep,3,1,ng)
1610 tc(ep,1,2)=vjfi(ep,1,1,ng)*vq(ep,1,2,ng)+vjfi(ep,2,1,ng)
1611 + *vq(ep,2,2,ng)+ vjfi(ep,3,1,ng)*vq(ep,3,2,ng)
1612 tc(ep,2,2)=vjfi(ep,1,2,ng)*vq(ep,1,2,ng)+vjfi(ep,2,2,ng)
1613 + *vq(ep,2,2,ng)+ vjfi(ep,3,2,ng)*vq(ep,3,2,ng)
1614C--------------------------
1615C CONSTRUIRE LA MATRICE [B] ET DEFORMATION
1616C--------------------------
1617C--------------------------
1618C DEFORMATION IN-PLANE
1619C--------------------------
1620C [BM](3,3*NPG),[BMF](3,3*NPG),
1621C--------------------------
1622C--------TERMES DE [B0]=<T>CI-------
1623 c1=vksi(1,ng)*tc(ep,1,1)+veta(1,ng)*tc(ep,2,1)
1624 c2=vksi(1,ng)*tc(ep,1,2)+veta(1,ng)*tc(ep,2,2)
1625C I = 1; JJ = 3*KK =0
1626c IF (ISROT>0) THEN
1627 bxy(1)=vq(ep,1,1,ng)*c2
1628 bxy(2)=vq(ep,2,1,ng)*c2
1629 bxy(3)=vq(ep,3,1,ng)*c2
1630 byx(1)=vq(ep,1,2,ng)*c1
1631 byx(2)=vq(ep,2,2,ng)*c1
1632 byx(3)=vq(ep,3,2,ng)*c1
1633 bm(ep,3)=bxy(1)+byx(1)
1634 bm(ep,6)=bxy(2)+byx(2)
1635 bm(ep,9)=bxy(3)+byx(3)
1636C---------CONTRIBUTION DU <T1> :EPS-X
1637 bm(ep,1)=vq(ep,1,1,ng)*c1
1638C---------CONTRIBUTION DU <T2> :EPS-Y
1639 bm(ep,2)=vq(ep,1,2,ng)*c2
1640C---------:EPS-X :EPS-Y :EPS-XY
1641 bm(ep,4)=vq(ep,2,1,ng)*c1
1642 bm(ep,5)=vq(ep,2,2,ng)*c2
1643C---------:EPS-X:EPS-Y:EPS-XY
1644 bm(ep,7)=vq(ep,3,1,ng)*c1
1645 bm(ep,8)=vq(ep,3,2,ng)*c2
1646C
1647C--------J=2---II=(J-1)*2 =2 KK = 3*(J-1) =3
1648C--------- [B0], [B01]---------
1649C--------TERMES DE [B0]=<T>CI-------
1650 c1=vksi(2,ng)*tc(ep,1,1)+veta(2,ng)*tc(ep,2,1)
1651 c2=vksi(2,ng)*tc(ep,1,2)+veta(2,ng)*tc(ep,2,2)
1652C I = 1; JJ = 3*KK =9
1653c IF (ISROT>0) THEN
1654 bxy(1)=vq(ep,1,1,ng)*c2
1655 bxy(2)=vq(ep,2,1,ng)*c2
1656 bxy(3)=vq(ep,3,1,ng)*c2
1657 byx(1)=vq(ep,1,2,ng)*c1
1658 byx(2)=vq(ep,2,2,ng)*c1
1659 byx(3)=vq(ep,3,2,ng)*c1
1660 bm(ep,12)=bxy(1)+byx(1)
1661 bm(ep,15)=bxy(2)+byx(2)
1662 bm(ep,18)=bxy(3)+byx(3)
1663C---------:EPS-X:EPS-Y:EPS-XY
1664 bm(ep,10)=vq(ep,1,1,ng)*c1
1665 bm(ep,11)=vq(ep,1,2,ng)*c2
1666C I = 2; JJ = 3*(KK+1)=12
1667C---------:EPS-X:EPS-Y:EPS-XY
1668 bm(ep,13)=vq(ep,2,1,ng)*c1
1669 bm(ep,14)=vq(ep,2,2,ng)*c2
1670C I = 3; JJ = 3*(KK+2)=15
1671C---------:EPS-X:EPS-Y:EPS-XY
1672 bm(ep,16)=vq(ep,3,1,ng)*c1
1673 bm(ep,17)=vq(ep,3,2,ng)*c2
1674C---------DEF DE FLEXION ASSOCIEE AUX THETA1,THETA2------
1675C--------J=3---II=(J-1)*2 =4 KK = 3*(J-1)=6
1676C--------- [B0], [B01]---------
1677C--------TERMES DE [B0]=<T>CI-------
1678 c1=vksi(3,ng)*tc(ep,1,1)+veta(3,ng)*tc(ep,2,1)
1679 c2=vksi(3,ng)*tc(ep,1,2)+veta(3,ng)*tc(ep,2,2)
1680c IF (ISROT>0) THEN
1681 bxy(1)=vq(ep,1,1,ng)*c2
1682 bxy(2)=vq(ep,2,1,ng)*c2
1683 bxy(3)=vq(ep,3,1,ng)*c2
1684 byx(1)=vq(ep,1,2,ng)*c1
1685 byx(2)=vq(ep,2,2,ng)*c1
1686 byx(3)=vq(ep,3,2,ng)*c1
1687 bm(ep,21)=bxy(1)+byx(1)
1688 bm(ep,24)=bxy(2)+byx(2)
1689 bm(ep,27)=bxy(3)+byx(3)
1690C---------:EPS-X:EPS-Y:EPS-XY
1691 bm(ep,19)=vq(ep,1,1,ng)*c1
1692 bm(ep,20)=vq(ep,1,2,ng)*c2
1693C I = 2; JJ = 3*(KK+1)=21
1694 bm(ep,22)=vq(ep,2,1,ng)*c1
1695 bm(ep,23)=vq(ep,2,2,ng)*c2
1696C I = 3; JJ = 3*(KK+2)=24
1697 bm(ep,25)=vq(ep,3,1,ng)*c1
1698 bm(ep,26)=vq(ep,3,2,ng)*c2
1699C---------- TERMES ASSOCIEES AU BETA ----------
1700C--------TERMES DE [B0]=<T>CI-------
1701 c1=vksi(4,ng)*tc(ep,1,1)+veta(4,ng)*tc(ep,2,1)
1702 c2=vksi(4,ng)*tc(ep,1,2)+veta(4,ng)*tc(ep,2,2)
1703c IF (ISROT>0) THEN
1704 bxy(1)=vq(ep,1,1,ng)*c2
1705 bxy(2)=vq(ep,2,1,ng)*c2
1706 bxy(3)=vq(ep,3,1,ng)*c2
1707 byx(1)=vq(ep,1,2,ng)*c1
1708 byx(2)=vq(ep,2,2,ng)*c1
1709 byx(3)=vq(ep,3,2,ng)*c1
1710 bm(ep,30)=bxy(1)+byx(1)
1711 bm(ep,33)=bxy(2)+byx(2)
1712 bm(ep,36)=bxy(3)+byx(3)
1713 bm(ep,28)=vq(ep,1,1,ng)*c1
1714 bm(ep,29)=vq(ep,1,2,ng)*c2
1715C I = 2; JJ = 3*(KK+1)=30
1716 bm(ep,31)=vq(ep,2,1,ng)*c1
1717 bm(ep,32)=vq(ep,2,2,ng)*c2
1718C I = 3; JJ = 3*(KK+2)=33
1719 bm(ep,34)=vq(ep,3,1,ng)*c1
1720 bm(ep,35)=vq(ep,3,2,ng)*c2
1721C
1722 150 CONTINUE
1723 RETURN

◆ cbpatch()

subroutine cbpatch ( integer jft,
integer jlt,
bm0rz,
bmkrz,
bmerz,
area,
mx13,
mx23,
mx34,
my13,
my23,
my34 )

Definition at line 1110 of file cbadef.F.

1112C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
1113#include "implicit_f.inc"
1114#include "mvsiz_p.inc"
1115C-----------------------------------------------
1116C D U M M Y A R G U M E N T S
1117C-----------------------------------------------
1118 INTEGER JFT,JLT
1119 my_real
1120 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
1121 . mx13(*),mx23(*),mx34(*),my13(*),my23(*),my34(*),area(*)
1122C-----------------------------------------------
1123C L O C A L V A R I A B L E S
1124C-----------------------------------------------
1125 INTEGER I,J,EP,NPG,NG,IAC
1126 parameter(npg = 4)
1127 my_real
1128 . bm0rzm(mvsiz,4,4),j1,j2,j0,jac(mvsiz,4),det,jk_a,je_a
1129 my_real
1130 . vpg(2,npg),pg,pg1
1131 parameter(pg=.577350269189626)
1132 parameter(pg1=-.577350269189626)
1133C--------------------------
1134C INITIALISATION
1135C--------------------------
1136 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
1137C-------------Correction of BM0RZ to pass Patch test- IAC: analytical compute
1138 iac = 1
1139 IF (iac>0) THEN
1140 DO i=jft,jlt
1141 jk_a=four_over_3*(mx23(i)*my13(i)-mx13(i)*my23(i))/area(i)
1142 je_a=four_over_3*(mx13(i)*my34(i)-mx34(i)*my13(i))/area(i)
1143C--------NxI,x ------
1144 bm0rz(i,1,1) = -bmkrz(i,1,1)*jk_a-bmerz(i,1,1)*je_a
1145 bm0rz(i,1,2) = -bmkrz(i,1,2)*jk_a-bmerz(i,1,2)*je_a
1146 bm0rz(i,1,3) = -bmkrz(i,1,3)*jk_a-bmerz(i,1,3)*je_a
1147 bm0rz(i,1,4) = -bmkrz(i,1,4)*jk_a-bmerz(i,1,4)*je_a
1148C -------NyI,y -------------
1149 bm0rz(i,2,1) = -bmkrz(i,2,1)*jk_a-bmerz(i,2,1)*je_a
1150 bm0rz(i,2,2) = -bmkrz(i,2,2)*jk_a-bmerz(i,2,2)*je_a
1151 bm0rz(i,2,3) = -bmkrz(i,2,3)*jk_a-bmerz(i,2,3)*je_a
1152 bm0rz(i,2,4) = -bmkrz(i,2,4)*jk_a-bmerz(i,2,4)*je_a
1153C--------NxI,y ------
1154 bm0rz(i,3,1) = -bmkrz(i,3,1)*jk_a-bmerz(i,3,1)*je_a
1155 bm0rz(i,3,2) = -bmkrz(i,3,2)*jk_a-bmerz(i,3,2)*je_a
1156 bm0rz(i,3,3) = -bmkrz(i,3,3)*jk_a-bmerz(i,3,3)*je_a
1157 bm0rz(i,3,4) = -bmkrz(i,3,4)*jk_a-bmerz(i,3,4)*je_a
1158C -------NyI,y -------------
1159 bm0rz(i,4,1) = -bmkrz(i,4,1)*jk_a-bmerz(i,4,1)*je_a
1160 bm0rz(i,4,2) = -bmkrz(i,4,2)*jk_a-bmerz(i,4,2)*je_a
1161 bm0rz(i,4,3) = -bmkrz(i,4,3)*jk_a-bmerz(i,4,3)*je_a
1162 bm0rz(i,4,4) = -bmkrz(i,4,4)*jk_a-bmerz(i,4,4)*je_a
1163 ENDDO
1164 RETURN
1165 END IF !(ICN>0) THEN
1166C
1167 DO i=jft,jlt
1168 j1=(mx23(i)*my13(i)-mx13(i)*my23(i))*pg
1169 j2=-(mx13(i)*my34(i)-mx34(i)*my13(i))*pg
1170 j0=area(i)*fourth
1171 jac(i,1)=j0+j2-j1
1172 jac(i,2)=j0+j2+j1
1173 jac(i,3)=j0-j2+j1
1174 jac(i,4)=j0-j2-j1
1175 ENDDO
1176C
1177 DO i=jft,jlt
1178 DO j=1,4
1179 bm0rzm(i,1,j)=zero
1180 bm0rzm(i,2,j)=zero
1181 bm0rzm(i,3,j)=zero
1182 bm0rzm(i,4,j)=zero
1183 ENDDO
1184 ENDDO
1185C
1186 DO ng= 1,4
1187 DO i=jft,jlt
1188 det=jac(i,ng)/area(i)
1189 DO j=1,4
1190 bm0rzm(i,1,j) = bm0rzm(i,1,j)+ (bmkrz(i,1,j)*vpg(1,ng)+
1191 . bmerz(i,1,j)*vpg(2,ng))*det
1192 bm0rzm(i,2,j) = bm0rzm(i,2,j)+(bmkrz(i,2,j)*vpg(1,ng)+
1193 . bmerz(i,2,j)*vpg(2,ng))*det
1194 bm0rzm(i,3,j) = bm0rzm(i,3,j)+(bmkrz(i,3,j)*vpg(1,ng)+
1195 . bmerz(i,3,j)*vpg(2,ng))*det
1196 bm0rzm(i,4,j) = bm0rzm(i,4,j)+(bmkrz(i,4,j)*vpg(1,ng)+
1197 . bmerz(i,4,j)*vpg(2,ng))*det
1198C
1199 ENDDO
1200 ENDDO
1201 ENDDO
1202C
1203 DO i=jft,jlt
1204 DO j=1,4
1205 bm0rz(i,1,j) = -bm0rzm(i,1,j)
1206 bm0rz(i,2,j) = -bm0rzm(i,2,j)
1207 bm0rz(i,3,j) = -bm0rzm(i,3,j)
1208 bm0rz(i,4,j) = -bm0rzm(i,4,j)
1209C
1210 ENDDO
1211 ENDDO
1212C
1213 RETURN