OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbabe3.F File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine cbabe3 (jft, jlt, ng, vcore, area, cdet, vqn, vq, vjfi, vnrm, vastn, hx, hy, veta, vksi, bm, bmf, bf, bc, tc, bzz, nplat, iplat, isrot, brz)
subroutine cbabec3 (jft, jlt, x13, x24, y13, y24, bm, nplat, iplat)
subroutine cbaber3 (jft, jlt, bm0rz, bmkrz, bmerz, bmrz, brz, bm, nplat, iplat, ng)

Function/Subroutine Documentation

◆ cbabe3()

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

Definition at line 34 of file cbabe3.F.

38C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
39C CALCUL [B] AU POINT D'INTEGRATION Index changed*****EP-->I*****
40C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
41C-----------------------------------------------
42C I M P L I C I T T Y P E S
43C-----------------------------------------------
44#include "implicit_f.inc"
45#include "mvsiz_p.inc"
46C-----------------------------------------------
47C D U M M Y A R G U M E N T S
48C-----------------------------------------------
49 INTEGER NPG,NG,JFT,JLT,NPLAT,IPLAT(*),ISROT
50 parameter(npg = 4)
51 my_real
52 . vcore(mvsiz,12),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),area(*)
57 my_real
58 . cdet(*),tc(mvsiz,2,2),bzz(mvsiz,8),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,NG1,K
63 my_real
64 . vpg(2,npg),thk(mvsiz),detj,det,
65 . tfn(3,2),bcx,bcy,v1(2),v2(2),c1,c2,bc1,bc2,vb1,
66 . vjf1(2,3),vjf(3,3),tbi(2,2),tbc(mvsiz,2,2),v11(4)
67 my_real
68 . a_1,c11,c12,c21,c22,cc,beta1,ksi1,ksiy1,beta2,ksi2,ksiy2,
69 . bmt(8),bct(12),bxy(3),byx(3)
70 my_real
71 . pg1,pg
72 parameter(pg=.577350269189626)
73 parameter(pg1=-.577350269189626)
74C--------------------------
75C INITIALISATION
76C--------------------------
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-------BM(1-4)->BxI,BM(5-8)->ByI
84C------- (I=1:BMT(1)+BMT(3);I=2:BMT(2)+BMT(4);I=3:-BMT(1)+BMT(3);I=4:-BMT(2)+BMT(4))
85 bmt(1)=vcore(ep,1)+hx(ep,ng)*vcore(ep,5)
86 bmt(2)=vcore(ep,2)+hx(ep,ng)*vcore(ep,6)
87 bmt(3)=hx(ep,ng)*fourth
88 bmt(4)=-bmt(3)
89 bm(i,1)=bmt(1)+bmt(3)
90 bm(i,2)=bmt(2)+bmt(4)
91 bm(i,3)=-bmt(1)+bmt(3)
92 bm(i,4)=-bmt(2)+bmt(4)
93 bmt(5)=vcore(ep,3)+hy(ep,ng)*vcore(ep,5)
94 bmt(6)=vcore(ep,4)+hy(ep,ng)*vcore(ep,6)
95 bmt(7)=hy(ep,ng)*fourth
96 bmt(8)=-bmt(7)
97 bm(i,5)=bmt(5)+bmt(7)
98 bm(i,6)=bmt(6)+bmt(8)
99 bm(i,7)=-bmt(5)+bmt(7)
100 bm(i,8)=-bmt(6)+bmt(8)
101 ENDDO
102 DO i=jft,nplat
103 ep=iplat(i)
104 a_1 = one/max(cdet(ep),em20)
105 c11=(vcore(ep,10)+vcore(ep,12)*vpg(1,ng))*a_1
106 c12=-(vcore(ep,8)+vcore(ep,12)*vpg(2,ng))*a_1
107 c21=-(vcore(ep,9)+vcore(ep,11)*vpg(1,ng))*a_1
108 c22=( vcore(ep,7)+vcore(ep,11)*vpg(2,ng))*a_1
109 bzz(i,1)=c11*vksi(1,ng)+c12*veta(1,ng)
110 bzz(i,2)=c21*vksi(1,ng)+c22*veta(1,ng)
111 bzz(i,3)=c11*vksi(2,ng)+c12*veta(2,ng)
112 bzz(i,4)=c21*vksi(2,ng)+c22*veta(2,ng)
113 bzz(i,5)=c11*vksi(3,ng)+c12*veta(3,ng)
114 bzz(i,6)=c21*vksi(3,ng)+c22*veta(3,ng)
115 bzz(i,7)=c11*vksi(4,ng)+c12*veta(4,ng)
116 bzz(i,8)=c21*vksi(4,ng)+c22*veta(4,ng)
117 c11=fourth*c11
118 c12=fourth*c12
119 c21=fourth*c21
120 c22=fourth*c22
121C--------C.T-----
122C-----PARTIE SYM---------
123 beta1=vcore(ep,8)+vcore(ep,12)*vpg(2,ng)
124 ksi1=vcore(ep,10)+vcore(ep,12)*vpg(1,ng)
125 beta2=vcore(ep,7)+vcore(ep,11)*vpg(2,ng)
126 ksi2=vcore(ep,9)+vcore(ep,11)*vpg(1,ng)
127C-------NOEUD--- 3
128 bct(1)=c11*vpg(2,ng)+c12*vpg(1,ng)
129 bct(2)=c21*vpg(2,ng)+c22*vpg(1,ng)
130 bct(3)=-beta1*c11-ksi1*c12
131 bct(4)=-beta1*c21-ksi1*c22
132 bct(5)=c11*beta2+c12*ksi2
133 bct(6)=c21*beta2+c22*ksi2
134C-------NOEUD--- 4
135 bct(7)=-bct(1)
136 bct(8)=-bct(2)
137 bct(9)=bct(3)
138 bct(10)=bct(4)
139 bct(11)=bct(5)
140 bct(12)=bct(6)
141C
142 beta1=vcore(ep,12)+vcore(ep,8)*vpg(2,ng)
143 ksi1=vcore(ep,12)+vcore(ep,10)*vpg(1,ng)
144 beta2=vcore(ep,11)+vcore(ep,7)*vpg(2,ng)
145 ksi2=vcore(ep,11)+vcore(ep,9)*vpg(1,ng)
146C-------NOEUD--- 1
147 bc(i,1)=-c11-c12+bct(1)
148 bc(i,2)=-c21-c22+bct(2)
149 bc(i,3)= beta1*c11+ksi1*c12+bct(3)
150 bc(i,4)= beta1*c21+ksi1*c22+bct(4)
151 bc(i,5)= -beta2*c11-ksi2*c12+bct(5)
152 bc(i,6)= -beta2*c21-ksi2*c22+bct(6)
153C-------NOEUD--- 2
154 bc(i,7)= c11-c12+bct(7)
155 bc(i,8)= c21-c22+bct(8)
156 bc(i,9)= beta1*c11-ksi1*c12+bct(9)
157 bc(i,10)=beta1*c21-ksi1*c22+bct(10)
158 bc(i,11)=-beta2*c11+ksi2*c12+bct(11)
159 bc(i,12)=-beta2*c21+ksi2*c22+bct(12)
160C-------NOEUD--- 3
161 bc(i,13)=-bc(i,1)+bct(1)+bct(1)
162 bc(i,14)=-bc(i,2)+bct(2)+bct(2)
163 bc(i,15)=-bc(i,3)+bct(3)+bct(3)
164 bc(i,16)=-bc(i,4)+bct(4)+bct(4)
165 bc(i,17)=-bc(i,5)+bct(5)+bct(5)
166 bc(i,18)=-bc(i,6)+bct(6)+bct(6)
167C-------NOEUD--- 4
168 bc(i,19)=-bc(i,7)+bct(7)+bct(7)
169 bc(i,20)=-bc(i,8)+bct(8)+bct(8)
170 bc(i,21)=-bc(i,9)+bct(9)+bct(9)
171 bc(i,22)=-bc(i,10)+bct(10)+bct(10)
172 bc(i,23)=-bc(i,11)+bct(11)+bct(11)
173 bc(i,24)=-bc(i,12)+bct(12)+bct(12)
174 ENDDO
175C
176#include "vectorize.inc"
177 DO i=nplat+1,jlt
178 ep=iplat(i)
179C---------------------------------------------------
180C CALCUL DE [FN]
181C---------------------------------------------------
182 tfn(1,1)=vksi(1,ng)*vqn(ep,7,1)+vksi(2,ng)*vqn(ep,7,2)
183 1 +vksi(3,ng)*vqn(ep,7,3)+vksi(4,ng)*vqn(ep,7,4)
184 tfn(2,1)=vksi(1,ng)*vqn(ep,8,1)+vksi(2,ng)*vqn(ep,8,2)
185 1 +vksi(3,ng)*vqn(ep,8,3)+vksi(4,ng)*vqn(ep,8,4)
186 tfn(3,1)=vksi(1,ng)*vqn(ep,9,1)+vksi(2,ng)*vqn(ep,9,2)
187 1 +vksi(3,ng)*vqn(ep,9,3)+vksi(4,ng)*vqn(ep,9,4)
188 tfn(1,2)=veta(1,ng)*vqn(ep,7,1)+veta(2,ng)*vqn(ep,7,2)
189 1 +veta(3,ng)*vqn(ep,7,3)+veta(4,ng)*vqn(ep,7,4)
190 tfn(2,2)=veta(1,ng)*vqn(ep,8,1)+veta(2,ng)*vqn(ep,8,2)
191 1 +veta(3,ng)*vqn(ep,8,3)+veta(4,ng)*vqn(ep,8,4)
192 tfn(3,2)=veta(1,ng)*vqn(ep,9,1)+veta(2,ng)*vqn(ep,9,2)
193 1 +veta(3,ng)*vqn(ep,9,3)+veta(4,ng)*vqn(ep,9,4)
194C--------------------------------------------------------
195C CALCUL [BN]=[F]^-1 [FN]
196C--------------------------------------------------------
197 tbi(2,2)=vjfi(ep,1,1,ng)*tfn(1,1)+vjfi(ep,2,1,ng)*tfn(2,1)
198 1 + vjfi(ep,3,1,ng)*tfn(3,1)
199 tbi(2,1)=vjfi(ep,1,2,ng)*tfn(1,1)+vjfi(ep,2,2,ng)*tfn(2,1)
200 1 + vjfi(ep,3,2,ng)*tfn(3,1)
201 tbi(1,2)=vjfi(ep,1,1,ng)*tfn(1,2)+vjfi(ep,2,1,ng)*tfn(2,2)
202 1 + vjfi(ep,3,1,ng)*tfn(3,2)
203 tbi(1,1)=vjfi(ep,1,2,ng)*tfn(1,2)+vjfi(ep,2,2,ng)*tfn(2,2)
204 1 + vjfi(ep,3,2,ng)*tfn(3,2)
205C
206 thk(ep) =-(tbi(1,1)+tbi(2,2))
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(ep,1,1)=tbi(1,1)*tc(ep,1,1)+tbi(1,2)*tc(ep,2,1)
224 tbc(ep,2,1)=tbi(2,1)*tc(ep,1,1)+tbi(2,2)*tc(ep,2,1)
225 tbc(ep,1,2)=tbi(1,1)*tc(ep,1,2)+tbi(1,2)*tc(ep,2,2)
226 tbc(ep,2,2)=tbi(2,1)*tc(ep,1,2)+tbi(2,2)*tc(ep,2,2)
227 ENDDO
228C
229C--------------------------
230C--------------------------
231C DEFORMATION IN-PLANE
232C--------------------------
233C [BM](3,3*NPG),[BMF](3,3*NPG),[BF](3,2*NPG)
234C--------------------------
235C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
236C---------- NOEUDS J=1 II =0 KK = 0
237#include "vectorize.inc"
238 DO i=nplat+1,jlt
239 ep=iplat(i)
240C--------TERMES DE [B0]=<T>CI-------
241 c1=vksi(1,ng)*tc(ep,1,1)+veta(1,ng)*tc(ep,2,1)
242 c2=vksi(1,ng)*tc(ep,1,2)+veta(1,ng)*tc(ep,2,2)
243C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
244 bc1=vksi(1,ng)*tbc(ep,1,1)+veta(1,ng)*tbc(ep,2,1)
245 bc2=vksi(1,ng)*tbc(ep,1,2)+veta(1,ng)*tbc(ep,2,2)
246C I = 1; JJ = 3*KK =0
247C---------CONTRIBUTION DU <T1> :EPS-X
248 IF (isrot>0) THEN
249 bxy(1)=vq(ep,1,1,ng)*c2
250 bxy(2)=vq(ep,2,1,ng)*c2
251 bxy(3)=vq(ep,3,1,ng)*c2
252 byx(1)=vq(ep,1,2,ng)*c1
253 byx(2)=vq(ep,2,2,ng)*c1
254 byx(3)=vq(ep,3,2,ng)*c1
255 brz(i,1,1)=-bxy(1)+byx(1)
256 brz(i,2,1)=-bxy(2)+byx(2)
257 brz(i,3,1)=-bxy(3)+byx(3)
258 bm(i,3)=bxy(1)+byx(1)
259 bm(i,6)=bxy(2)+byx(2)
260 bm(i,9)=bxy(3)+byx(3)
261 ELSE
262 bm(i,9)=zero
263 END IF
264 bm(i,1)=vq(ep,1,1,ng)*c1
265 bmf(i,1)=thk(ep)*bm(i,1)+vq(ep,1,1,ng)*bc1
266C---------CONTRIBUTION DU <T2> :EPS-Y
267 bm(i,2)=vq(ep,1,2,ng)*c2
268 bmf(i,2)=thk(ep)*bm(i,2)+vq(ep,1,2,ng)*bc2
269C---------CONTRIBUTION DU <T> :EPS-XY
270 bmf(i,3)=thk(ep)*bm(i,3)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
271C I = 2; JJ = 3*(KK+1)=3
272C---------:EPS-X :EPS-Y :EPS-XY
273 bm(i,4)=vq(ep,2,1,ng)*c1
274 bmf(i,4)=thk(ep)*bm(i,4)+vq(ep,2,1,ng)*bc1
275 bm(i,5)=vq(ep,2,2,ng)*c2
276 bmf(i,5)=thk(ep)*bm(i,5)+vq(ep,2,2,ng)*bc2
277 bmf(i,6)=thk(ep)*bm(i,6)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
278C I = 3; JJ = 3*(KK+2)=6
279C---------:EPS-X:EPS-Y:EPS-XY
280 bm(i,7)=vq(ep,3,1,ng)*c1
281 bmf(i,7)=thk(ep)*bm(i,7)+vq(ep,3,1,ng)*bc1
282 bm(i,8)=vq(ep,3,2,ng)*c2
283 bmf(i,8)=thk(ep)*bm(i,8)+vq(ep,3,2,ng)*bc2
284 bmf(i,9)=thk(ep)*bm(i,9)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
285C---------- TERMES ASSOCIEES AU BETA ----------
286C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
287 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,1)+vq(ep,2,1,ng)*vqn(ep,2,1)
288 1 +vq(ep,3,1,ng)*vqn(ep,3,1)
289 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,1)+vq(ep,2,1,ng)*vqn(ep,5,1)
290 1 +vq(ep,3,1,ng)*vqn(ep,6,1))
291 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,1)+vq(ep,2,2,ng)*vqn(ep,2,1)
292 1 +vq(ep,3,2,ng)*vqn(ep,3,1)
293 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,1)+vq(ep,2,2,ng)*vqn(ep,5,1)
294 1 +vq(ep,3,2,ng)*vqn(ep,6,1))
295C----------TERMES ASSOCIEES AUX BETA DE [B01]--
296C JJ = 3*II =0
297 bf(i,1)=v1(1)*c1
298 bf(i,2)=v2(1)*c2
299 bf(i,3)=v1(1)*c2+v2(1)*c1
300 bf(i,4)=v1(2)*c1
301 bf(i,5)=v2(2)*c2
302 bf(i,6)=v1(2)*c2+v2(2)*c1
303 ENDDO
304C
305C--------J=2---II=(J-1)*2 =2 KK = 3*(J-1) =3
306#include "vectorize.inc"
307 DO i=nplat+1,jlt
308 ep=iplat(i)
309C--------- [B0], [B01]---------
310C--------TERMES DE [B0]=<T>CI-------
311 c1=vksi(2,ng)*tc(ep,1,1)+veta(2,ng)*tc(ep,2,1)
312 c2=vksi(2,ng)*tc(ep,1,2)+veta(2,ng)*tc(ep,2,2)
313C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
314 bc1=vksi(2,ng)*tbc(ep,1,1)+veta(2,ng)*tbc(ep,2,1)
315 bc2=vksi(2,ng)*tbc(ep,1,2)+veta(2,ng)*tbc(ep,2,2)
316C
317C I = 1; JJ = 3*KK =9
318C---------:EPS-X:EPS-Y:EPS-XY
319 IF (isrot>0) THEN
320 bxy(1)=vq(ep,1,1,ng)*c2
321 bxy(2)=vq(ep,2,1,ng)*c2
322 bxy(3)=vq(ep,3,1,ng)*c2
323 byx(1)=vq(ep,1,2,ng)*c1
324 byx(2)=vq(ep,2,2,ng)*c1
325 byx(3)=vq(ep,3,2,ng)*c1
326 brz(i,1,2)=-bxy(1)+byx(1)
327 brz(i,2,2)=-bxy(2)+byx(2)
328 brz(i,3,2)=-bxy(3)+byx(3)
329 bm(i,12)=bxy(1)+byx(1)
330 bm(i,15)=bxy(2)+byx(2)
331 bm(i,18)=bxy(3)+byx(3)
332 ELSE
333 bm(i,18)=zero
334 END IF
335 bm(i,10)=vq(ep,1,1,ng)*c1
336 bmf(i,10)=thk(ep)*bm(i,10)+vq(ep,1,1,ng)*bc1
337 bm(i,11)=vq(ep,1,2,ng)*c2
338 bmf(i,11)=thk(ep)*bm(i,11)+vq(ep,1,2,ng)*bc2
339 bmf(i,12)=thk(ep)*bm(i,12)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
340C
341C I = 2; JJ = 3*(KK+1)=12
342C---------:EPS-X:EPS-Y:EPS-XY
343 bm(i,13)=vq(ep,2,1,ng)*c1
344 bmf(i,13)=thk(ep)*bm(i,13)+vq(ep,2,1,ng)*bc1
345 bm(i,14)=vq(ep,2,2,ng)*c2
346 bmf(i,14)=thk(ep)*bm(i,14)+vq(ep,2,2,ng)*bc2
347 bmf(i,15)=thk(ep)*bm(i,15)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
348C
349C I = 3; JJ = 3*(KK+2)=15
350C---------:EPS-X:EPS-Y:EPS-XY
351 bm(i,16)=vq(ep,3,1,ng)*c1
352 bmf(i,16)=thk(ep)*bm(i,16)+vq(ep,3,1,ng)*bc1
353 bm(i,17)=vq(ep,3,2,ng)*c2
354 bmf(i,17)=thk(ep)*bm(i,17)+vq(ep,3,2,ng)*bc2
355 bmf(i,18)=thk(ep)*bm(i,18)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
356C
357C---------- TERMES ASSOCIEES AU BETA ----------
358C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
359 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,2)+vq(ep,2,1,ng)*vqn(ep,2,2)
360 1 +vq(ep,3,1,ng)*vqn(ep,3,2)
361 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,2)+vq(ep,2,1,ng)*vqn(ep,5,2)
362 1 +vq(ep,3,1,ng)*vqn(ep,6,2))
363 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,2)+vq(ep,2,2,ng)*vqn(ep,2,2)
364 1 +vq(ep,3,2,ng)*vqn(ep,3,2)
365 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,2)+vq(ep,2,2,ng)*vqn(ep,5,2)
366 1 +vq(ep,3,2,ng)*vqn(ep,6,2))
367C----------TERMES ASSOCIEES AUX BETA DE [B01]--
368C JJ = 3*II=6
369 bf(i,7)=v1(1)*c1
370 bf(i,8)=v2(1)*c2
371 bf(i,9)=v1(1)*c2+v2(1)*c1
372 bf(i,10)=v1(2)*c1
373 bf(i,11)=v2(2)*c2
374 bf(i,12)=v1(2)*c2+v2(2)*c1
375 ENDDO
376C
377C--------J=3---II=(J-1)*2 =4 KK = 3*(J-1)=6
378#include "vectorize.inc"
379 DO i=nplat+1,jlt
380 ep=iplat(i)
381C--------- [B0], [B01]---------
382C--------TERMES DE [B0]=<T>CI-------
383 c1=vksi(3,ng)*tc(ep,1,1)+veta(3,ng)*tc(ep,2,1)
384 c2=vksi(3,ng)*tc(ep,1,2)+veta(3,ng)*tc(ep,2,2)
385C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
386 bc1=vksi(3,ng)*tbc(ep,1,1)+veta(3,ng)*tbc(ep,2,1)
387 bc2=vksi(3,ng)*tbc(ep,1,2)+veta(3,ng)*tbc(ep,2,2)
388C I = 1;JJ = 3*KK=18
389C---------:EPS-X:EPS-Y:EPS-XY
390 IF (isrot>0) THEN
391 bxy(1)=vq(ep,1,1,ng)*c2
392 bxy(2)=vq(ep,2,1,ng)*c2
393 bxy(3)=vq(ep,3,1,ng)*c2
394 byx(1)=vq(ep,1,2,ng)*c1
395 byx(2)=vq(ep,2,2,ng)*c1
396 byx(3)=vq(ep,3,2,ng)*c1
397 brz(i,1,3)=-bxy(1)+byx(1)
398 brz(i,2,3)=-bxy(2)+byx(2)
399 brz(i,3,3)=-bxy(3)+byx(3)
400 bm(i,21)=bxy(1)+byx(1)
401 bm(i,24)=bxy(2)+byx(2)
402 bm(i,27)=bxy(3)+byx(3)
403 ELSE
404 bm(i,27)=zero
405 END IF
406 bm(i,19)=vq(ep,1,1,ng)*c1
407 bmf(i,19)=thk(ep)*bm(i,19)+vq(ep,1,1,ng)*bc1
408 bm(i,20)=vq(ep,1,2,ng)*c2
409 bmf(i,20)=thk(ep)*bm(i,20)+vq(ep,1,2,ng)*bc2
410 bmf(i,21)=thk(ep)*bm(i,21)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
411C I = 2; JJ = 3*(KK+1)=21
412 bm(i,22)=vq(ep,2,1,ng)*c1
413 bmf(i,22)=thk(ep)*bm(i,22)+vq(ep,2,1,ng)*bc1
414 bm(i,23)=vq(ep,2,2,ng)*c2
415 bmf(i,23)=thk(ep)*bm(i,23)+vq(ep,2,2,ng)*bc2
416 bmf(i,24)=thk(ep)*bm(i,24)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
417C I = 3; JJ = 3*(KK+2)=24
418 bm(i,25)=vq(ep,3,1,ng)*c1
419 bmf(i,25)=thk(ep)*bm(i,25)+vq(ep,3,1,ng)*bc1
420 bm(i,26)=vq(ep,3,2,ng)*c2
421 bmf(i,26)=thk(ep)*bm(i,26)+vq(ep,3,2,ng)*bc2
422 bmf(i,27)=thk(ep)*bm(i,27)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
423C
424C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
425C----------[B01]---------
426C---------- TERMES ASSOCIEES AU BETA ----------
427C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
428 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,3)+vq(ep,2,1,ng)*vqn(ep,2,3)
429 1 +vq(ep,3,1,ng)*vqn(ep,3,3)
430 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,3)+vq(ep,2,1,ng)*vqn(ep,5,3)
431 1 +vq(ep,3,1,ng)*vqn(ep,6,3))
432 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,3)+vq(ep,2,2,ng)*vqn(ep,2,3)
433 1 +vq(ep,3,2,ng)*vqn(ep,3,3)
434 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,3)+vq(ep,2,2,ng)*vqn(ep,5,3)
435 1 +vq(ep,3,2,ng)*vqn(ep,6,3))
436C----------TERMES ASSOCIEES AUX BETA DE [B01]--
437C JJ = 3*II =12
438 bf(i,13)=v1(1)*c1
439 bf(i,14)=v2(1)*c2
440 bf(i,15)=v1(1)*c2+v2(1)*c1
441 bf(i,16)=v1(2)*c1
442 bf(i,17)=v2(2)*c2
443 bf(i,18)=v1(2)*c2+v2(2)*c1
444 ENDDO
445C
446C--------J=4---II=(J-1)*2 =6 KK = 3*(J-1)=9
447#include "vectorize.inc"
448 DO i=nplat+1,jlt
449 ep=iplat(i)
450C--------- [B0], [B01]---------
451C--------TERMES DE [B0]=<T>CI-------
452 c1=vksi(4,ng)*tc(ep,1,1)+veta(4,ng)*tc(ep,2,1)
453 c2=vksi(4,ng)*tc(ep,1,2)+veta(4,ng)*tc(ep,2,2)
454C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
455 bc1=vksi(4,ng)*tbc(ep,1,1)+veta(4,ng)*tbc(ep,2,1)
456 bc2=vksi(4,ng)*tbc(ep,1,2)+veta(4,ng)*tbc(ep,2,2)
457C I = 1; JJ = 3*KK =27
458 bm(i,28)=vq(ep,1,1,ng)*c1
459 bmf(i,28)=thk(ep)*bm(i,28)+vq(ep,1,1,ng)*bc1
460 bm(i,29)=vq(ep,1,2,ng)*c2
461 bmf(i,29)=thk(ep)*bm(i,29)+vq(ep,1,2,ng)*bc2
462C---------CONTRIBUTION DU <T> :EPS-XY
463 IF (isrot>0) THEN
464 bxy(1)=vq(ep,1,1,ng)*c2
465 bxy(2)=vq(ep,2,1,ng)*c2
466 bxy(3)=vq(ep,3,1,ng)*c2
467 byx(1)=vq(ep,1,2,ng)*c1
468 byx(2)=vq(ep,2,2,ng)*c1
469 byx(3)=vq(ep,3,2,ng)*c1
470 brz(i,1,4)=-bxy(1)+byx(1)
471 brz(i,2,4)=-bxy(2)+byx(2)
472 brz(i,3,4)=-bxy(3)+byx(3)
473 bm(i,30)=bxy(1)+byx(1)
474 bm(i,33)=bxy(2)+byx(2)
475 bm(i,36)=bxy(3)+byx(3)
476 ELSE
477 bm(i,36)=zero
478 END IF
479 bmf(i,30)=thk(ep)*bm(i,30)+vq(ep,1,1,ng)*bc2+vq(ep,1,2,ng)*bc1
480C I = 2; JJ = 3*(KK+1)=30
481 bm(i,31)=vq(ep,2,1,ng)*c1
482 bmf(i,31)=thk(ep)*bm(i,31)+vq(ep,2,1,ng)*bc1
483 bm(i,32)=vq(ep,2,2,ng)*c2
484 bmf(i,32)=thk(ep)*bm(i,32)+vq(ep,2,2,ng)*bc2
485 bmf(i,33)=thk(ep)*bm(i,33)+vq(ep,2,1,ng)*bc2+vq(ep,2,2,ng)*bc1
486C I = 3; JJ = 3*(KK+2)=33
487 bm(i,34)=vq(ep,3,1,ng)*c1
488 bmf(i,34)=thk(ep)*bm(i,34)+vq(ep,3,1,ng)*bc1
489 bm(i,35)=vq(ep,3,2,ng)*c2
490 bmf(i,35)=thk(ep)*bm(i,35)+vq(ep,3,2,ng)*bc2
491 bmf(i,36)=thk(ep)*bm(i,36)+vq(ep,3,1,ng)*bc2+vq(ep,3,2,ng)*bc1
492C
493C---------- TERMES ASSOCIEES AU BETA ----------
494C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
495 v1(2)=vq(ep,1,1,ng)*vqn(ep,1,4)+vq(ep,2,1,ng)*vqn(ep,2,4)
496 1 +vq(ep,3,1,ng)*vqn(ep,3,4)
497 v1(1)=-(vq(ep,1,1,ng)*vqn(ep,4,4)+vq(ep,2,1,ng)*vqn(ep,5,4)
498 1 +vq(ep,3,1,ng)*vqn(ep,6,4))
499 v2(2)=vq(ep,1,2,ng)*vqn(ep,1,4)+vq(ep,2,2,ng)*vqn(ep,2,4)
500 1 +vq(ep,3,2,ng)*vqn(ep,3,4)
501 v2(1)=-(vq(ep,1,2,ng)*vqn(ep,4,4)+vq(ep,2,2,ng)*vqn(ep,5,4)
502 1 +vq(ep,3,2,ng)*vqn(ep,6,4))
503C----------TERMES ASSOCIEES AUX BETA DE [B01]--
504C JJ = 3*II=18
505 bf(i,19)=v1(1)*c1
506 bf(i,20)=v2(1)*c2
507 bf(i,21)=v1(1)*c2+v2(1)*c1
508 bf(i,22)=v1(2)*c1
509 bf(i,23)=v2(2)*c2
510 bf(i,24)=v1(2)*c2+v2(2)*c1
511C
512 ENDDO
513C--------------------------
514C DEFORMATION OUT-PLANE (C.T)
515C--------------------------
516C [BC](2,5*NPG),VDEF(4,5)
517C--------------------------
518#include "vectorize.inc"
519 DO i=nplat+1,jlt
520 ep=iplat(i)
521 v11(1)=vksi(2,ng)
522 v11(2)=vksi(3,ng)
523 v11(3)=veta(4,ng)
524 v11(4)=veta(3,ng)
525C--- J=1---------------------
526C-----------POINTEURS POUR BC(II,J1~J2)
527C J1=(KA1(1,J)-1)*5=0
528C J2=(KA1(2,J)-1)*5=5
529C II=KIND(J) =1
530C JJ = 3*(J-1) =0
531C--------TERMES <NAI>-------
532 c1=v11(1)*vnrm(ep,1)
533 bc(i,1)=-c1
534 bc(i,11)=c1
535C
536 c1=v11(1)*vnrm(ep,2)
537 bc(i,3)=-c1
538 bc(i,13)=c1
539C
540 c1=v11(1)*vnrm(ep,3)
541 bc(i,5)=-c1
542 bc(i,15)=c1
543C--------TERMES <ASI>-------
544 bc(i,7)=v11(1)*vastn(ep,1)
545 bc(i,9)=v11(1)*vastn(ep,2)
546 bc(i,17)=v11(1)*vastn(ep,3)
547 bc(i,19)=v11(1)*vastn(ep,4)
548C--- J=2---------------------
549C-----------POINTEURS POUR BC(II,J1~J2)
550C J1=(KA1(1,J)-1)*5 =15
551C J2=(KA1(2,J)-1)*5 =10
552C II=KIND(J) =1
553C JJ = 3*(J-1) =3
554C--------TERMES <NAI>-------
555 c1=v11(2)*vnrm(ep,4)
556 bc(i,31)=-c1
557 bc(i,21)=c1
558C
559 c1=v11(2)*vnrm(ep,5)
560 bc(i,33)=-c1
561 bc(i,23)=c1
562C
563 c1=v11(2)*vnrm(ep,6)
564 bc(i,35)=-c1
565 bc(i,25)=c1
566C--------TERMES <ASI>-------
567 bc(i,37)=v11(2)*vastn(ep,5)
568 bc(i,39)=v11(2)*vastn(ep,6)
569 bc(i,27)=v11(2)*vastn(ep,7)
570 bc(i,29)=v11(2)*vastn(ep,8)
571C--- J=3---------------------
572C-----------POINTEURS POUR BC(II,J1~J2)
573C J1=(KA1(1,J)-1)*5 =0
574C J2=(KA1(2,J)-1)*5 =15
575C II=KIND(J) =2
576C JJ = 3*(J-1) =6
577C--------TERMES <NAI>-------
578 c1=v11(3)*vnrm(ep,7)
579 bc(i,2)=-c1
580 bc(i,32)=c1
581C
582 c1=v11(3)*vnrm(ep,8)
583 bc(i,4)=-c1
584 bc(i,34)=c1
585C
586 c1=v11(3)*vnrm(ep,9)
587 bc(i,6)=-c1
588 bc(i,36)=c1
589C--------TERMES <ASI>-------
590 bc(i,8)=v11(3)*vastn(ep,9)
591 bc(i,10)=v11(3)*vastn(ep,10)
592 bc(i,38)=v11(3)*vastn(ep,11)
593 bc(i,40)=v11(3)*vastn(ep,12)
594C--- J=4---------------------
595C-----------POINTEURS POUR BC(II,J1~J2)
596C J1=(KA1(1,J)-1)*5 =5
597C J2=(KA1(2,J)-1)*5 =10
598C II=KIND(J) =2
599C JJ = 3*(J-1) =9
600C--------TERMES <NAI>-------
601 c1=v11(4)*vnrm(ep,10)
602 bc(i,12)=-c1
603 bc(i,22)=c1
604C
605 c1=v11(4)*vnrm(ep,11)
606 bc(i,14)=-c1
607 bc(i,24)=c1
608C
609 c1=v11(4)*vnrm(ep,12)
610 bc(i,16)=-c1
611 bc(i,26)=c1
612C--------TERMES <ASI>-------
613 bc(i,18)=v11(4)*vastn(ep,13)
614 bc(i,20)=v11(4)*vastn(ep,14)
615 bc(i,28)=v11(4)*vastn(ep,15)
616 bc(i,30)=v11(4)*vastn(ep,16)
617 ENDDO
618C--------pour rigidite en rzz-------
619 IF (isrot==0) THEN
620#include "vectorize.inc"
621 DO i=nplat+1,jlt
622 ep=iplat(i)
623 bzz(i,1)=tc(ep,1,1)*vksi(1,ng)+tc(ep,2,1)*veta(1,ng)
624 bzz(i,2)=tc(ep,1,2)*vksi(1,ng)+tc(ep,2,2)*veta(1,ng)
625 bzz(i,3)=tc(ep,1,1)*vksi(2,ng)+tc(ep,2,1)*veta(2,ng)
626 bzz(i,4)=tc(ep,1,2)*vksi(2,ng)+tc(ep,2,2)*veta(2,ng)
627 bzz(i,5)=tc(ep,1,1)*vksi(3,ng)+tc(ep,2,1)*veta(3,ng)
628 bzz(i,6)=tc(ep,1,2)*vksi(3,ng)+tc(ep,2,2)*veta(3,ng)
629 bzz(i,7)=tc(ep,1,1)*vksi(4,ng)+tc(ep,2,1)*veta(4,ng)
630 bzz(i,8)=tc(ep,1,2)*vksi(4,ng)+tc(ep,2,2)*veta(4,ng)
631 ENDDO
632 END IF !(ISROT==0) THEN
633 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

◆ cbabec3()

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

Definition at line 640 of file cbabe3.F.

642C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
643C CALCUL 'membrane shear traitement' use only PARTIE CONSTANTE
644C-----------------------------------------------
645C I M P L I C I T T Y P E S
646C-----------------------------------------------
647#include "implicit_f.inc"
648#include "mvsiz_p.inc"
649C-----------------------------------------------
650C D U M M Y A R G U M E N T S
651C-----------------------------------------------
652 INTEGER JFT ,JLT,NPLAT,IPLAT(*)
653 my_real
654 . x13(*),x24(*),y13(*),y24(*),bm(mvsiz,36)
655C-----------------------------------------------
656C L O C A L V A R I A B L E S
657C-----------------------------------------------
658 INTEGER I ,J,EP
659#include "vectorize.inc"
660 DO ep=nplat+1,jlt
661 i=iplat(ep)
662 bm(ep,3)=-x24(i)
663 bm(ep,6)=y24(i)
664 bm(ep,12)=x13(i)
665 bm(ep,15)=-y13(i)
666 bm(ep,21)=-bm(ep,3)
667 bm(ep,24)=-bm(ep,6)
668 bm(ep,30)=-bm(ep,12)
669 bm(ep,33)=-bm(ep,15)
670 bm(ep,9)=zero
671 bm(ep,18)=zero
672 bm(ep,27)=zero
673 bm(ep,36)=zero
674 ENDDO
675 RETURN

◆ cbaber3()

subroutine cbaber3 ( integer jft,
integer jlt,
bm0rz,
bmkrz,
bmerz,
bmrz,
brz,
bm,
integer nplat,
integer, dimension(*) iplat,
integer ng )

Definition at line 682 of file cbabe3.F.

685C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
686#include "implicit_f.inc"
687#include "mvsiz_p.inc"
688C-----------------------------------------------
689C D U M M Y A R G U M E N T S
690C-----------------------------------------------
691 INTEGER JFT,JLT,IPLAT(*),NPLAT,NG
692 my_real
693 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
694 . bmrz(mvsiz,3,4),bm(mvsiz,36),brz(mvsiz,4,4)
695C-----------------------------------------------
696C L O C A L V A R I A B L E S
697C-----------------------------------------------
698 INTEGER I,J,EP,NPG
699 parameter(npg = 4)
700 my_real
701 . vpg(2,npg),pg,pg1
702 parameter(pg=.577350269189626)
703 parameter(pg1=-.577350269189626)
704C--------------------------
705C INITIALISATION
706C--------------------------
707 DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
708C---------------
709C--------pay attention QEPH: Bi*A--- RLZ: VRZ/A-- QBAT:Bi--RLZ:VRZ
710 DO j=1,4
711 DO i=jft,jlt
712 bmrz(i,1,j) = bm0rz(i,1,j)+bmkrz(i,1,j)*vpg(1,ng)+
713 . bmerz(i,1,j)*vpg(2,ng)
714 bmrz(i,2,j) = bm0rz(i,2,j)+bmkrz(i,2,j)*vpg(1,ng)+
715 . bmerz(i,2,j)*vpg(2,ng)
716 bmrz(i,3,j) = bm0rz(i,3,j)+bmkrz(i,3,j)*vpg(1,ng)+
717 . bmerz(i,3,j)*vpg(2,ng)
718C
719 brz(i,4,j) = bm0rz(i,4,j)+bmkrz(i,4,j)*vpg(1,ng)+
720 . bmerz(i,4,j)*vpg(2,ng)
721 ENDDO
722 ENDDO
723C
724 RETURN