OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbafori.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| cbafori ../engine/source/elements/shell/coqueba/cbafori.f
25!||--- called by ------------------------------------------------------
26!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
27!||====================================================================
28 SUBROUTINE cbafori(JFT ,JLT ,NG ,CDET ,THK0,
29 2 TH12 ,FF0 ,MM0 ,NEL ,BM ,
30 3 BMF ,BF ,BC ,TC ,VF ,
31 4 VM ,NPLAT ,IPLAT ,VOL )
32C--------------------------------------------------------------------------------------------------
33C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
34C CALCUL DES FORCES INTERNES AUX NOEUDS
35C ENTREES : NEL,MNT,NPT,NP,NG,NIXC,IXC,HZ,OFF,CDET,EPAIS,BUFCON,VQN,BM,BMF,BF,BC
36C SORTIES : VF,VM
37C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
38C--------------------------------------------------------------------------------------------------
39C-----------------------------------------------
40C I M P L I C I T T Y P E S
41C-----------------------------------------------
42#include "implicit_f.inc"
43#include "mvsiz_p.inc"
44C-----------------------------------------------
45C D U M M Y A R G U M E N T S
46C-----------------------------------------------
47 INTEGER JFT,JLT,NG,NPLAT,IPLAT(*),NEL
48 MY_REAL
49 . FF0(NEL,5),MM0(NEL,3),CDET(*),VOL(*),
50 . bm(mvsiz,36),bmf(mvsiz,36),bf(mvsiz,24),bc(mvsiz,40),
51 . thk0(*),vf(mvsiz,3,4),vm(mvsiz,2,4),th12(*),
52 . c1,c2,cm1,cm2,cc1,cc2,bc1,bc3,bc5,
53 . tc(mvsiz,2,2),bcx,bcy
54C-----------------------------------------------
55C L O C A L V A R I A B L E S
56C-----------------------------------------------
57 INTEGER EP,I,J
58 MY_REAL
59 . FF(MVSIZ,5),MM(MVSIZ,3)
60C--------------------------------------------------------------------------------------------------
61C---due to bug of convention-----
62 DO ep=jft,jlt
63 ff(ep,1) = ff0(ep,1)
64 ff(ep,2) = ff0(ep,2)
65 ff(ep,3) = ff0(ep,3)
66 ff(ep,5) = ff0(ep,4)
67 ff(ep,4) = ff0(ep,5)
68
69 mm(ep,1) = mm0(ep,1)
70 mm(ep,2) = mm0(ep,2)
71 mm(ep,3) = mm0(ep,3)
72 END DO !
73#include "vectorize.inc"
74 DO i=jft,nplat
75 ep=iplat(i)
76 c2=vol(ep)
77 c1=th12(ep)*cdet(ep)
78 cm1=c1*(bm(ep,7)*mm(ep,2)+bm(ep,3)*mm(ep,3))
79 cm2=c1*(bm(ep,3)*mm(ep,1)+bm(ep,7)*mm(ep,3))
80 cc1=c2*(bc(ep,15)*ff(ep,4)+bc(ep,16)*ff(ep,5))
81 cc2=c2*(bc(ep,17)*ff(ep,4)+bc(ep,18)*ff(ep,5))
82 bc1=bc(ep,1)*ff(ep,4)
83 bc3=bc(ep,3)*ff(ep,4)
84 bc5=bc(ep,5)*ff(ep,4)
85C----------NOEUD 1 ------------
86 vf(ep,1,1)=vf(ep,1,1)+c2*(bm(ep,1)*ff(ep,1))
87 vf(ep,2,1)=vf(ep,2,1)+c2*(bm(ep,5)*ff(ep,2))
88 vf(ep,3,1)=vf(ep,3,1)+c2*(bc1+bc(ep,2)*ff(ep,5))
89C
90 vm(ep,1,1)=vm(ep,1,1)+c2*(bc3+bc(ep,4)*ff(ep,5))
91 1 -c1*(bm(ep,5)*mm(ep,2)+bm(ep,1)*mm(ep,3))
92 vm(ep,2,1)=vm(ep,2,1)+c2*(bc5+bc(ep,6)*ff(ep,5))+
93 1 c1*(bm(ep,1)*mm(ep,1)+bm(ep,5)*mm(ep,3))
94C
95 vf(ep,1,3)=vf(ep,1,3)+c2*(bm(ep,3)*ff(ep,1))
96 vf(ep,2,3)=vf(ep,2,3)+c2*(bm(ep,7)*ff(ep,2))
97C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
98 vf(ep,3,3)=vf(ep,3,3)+c2*(bc(ep,13)*ff(ep,4)+bc(ep,14)*ff(ep,5))
99C
100 vm(ep,1,3)=vm(ep,1,3)+cc1-cm1
101 vm(ep,2,3)=vm(ep,2,3)+cc2+cm2
102C----------NOEUD 2 ------------
103 vf(ep,1,2)=vf(ep,1,2)+c2*(bm(ep,2)*ff(ep,1))
104 vf(ep,2,2)=vf(ep,2,2)+c2*(bm(ep,6)*ff(ep,2))
105 vf(ep,3,2)=vf(ep,3,2)+c2*(bc(ep,7)*ff(ep,4)+bc(ep,8)*ff(ep,5))
106 vm(ep,1,2)=vm(ep,1,2)+c2*(bc(ep,9)*ff(ep,4)+bc(ep,10)*ff(ep,5))
107 1 -c1*(bm(ep,6)*mm(ep,2)+bm(ep,2)*mm(ep,3))
108 vm(ep,2,2)=vm(ep,2,2)
109 1 +c2*(bc(ep,11)*ff(ep,4)+bc(ep,12)*ff(ep,5))+
110 1 c1*(bm(ep,2)*mm(ep,1)+bm(ep,6)*mm(ep,3))
111C
112 vf(ep,1,4)=-vf(ep,1,3)
113 vf(ep,2,4)=-vf(ep,2,3)
114 vf(ep,3,4)=-vf(ep,3,3)
115C
116 vm(ep,1,4)=vm(ep,1,4)+cc1+cm1
117 vm(ep,2,4)=vm(ep,2,4)+cc2-cm2
118C
119 ENDDO
120#include "vectorize.inc"
121 DO i=nplat+1,jlt
122 ep=iplat(i)
123 c2=vol(ep)
124 c1=th12(ep)*cdet(ep)
125 bcx=tc(ep,1,1)*ff(ep,4)+tc(ep,1,2)*ff(ep,5)
126 bcy=tc(ep,2,1)*ff(ep,4)+tc(ep,2,2)*ff(ep,5)
127C J=1 K0=0 K1=0 K2=0
128 vf(ep,1,1)=vf(ep,1,1) +c2*(bm(ep,1)*ff(ep,1)+
129 1 bm(ep,2)*ff(ep,2)+
130 3 bc(ep,1)*bcx+
131 4 bc(ep,2)*bcy )
132 5 +c1*(bmf(ep,1)*mm(ep,1)+
133 6 bmf(ep,2)*mm(ep,2)+
134 7 bmf(ep,3)*mm(ep,3))
135C I1=1 I2=4
136 vm(ep,1,1)=vm(ep,1,1)+c2*(bc(ep,7)*bcx+
137 1 bc(ep,8)*bcy)
138 2 +c1*(bf(ep,1)*mm(ep,1)+
139 3 bf(ep,2)*mm(ep,2)+
140 4 bf(ep,3)*mm(ep,3) )
141C I1=2 I2=2
142 vf(ep,2,1)=vf(ep,2,1)+c2*(bm(ep,4)*ff(ep,1)+
143 2 bm(ep,5)*ff(ep,2)+
144 4 bc(ep,3)*bcx+
145 5 bc(ep,4)*bcy )
146 6 +c1*(bmf(ep,4)*mm(ep,1)+
147 7 bmf(ep,5)*mm(ep,2)+
148 8 bmf(ep,6)*mm(ep,3))
149C I1=2 I2=5
150 vm(ep,2,1)=vm(ep,2,1)+c2*(bc(ep,9)*bcx+
151 1 bc(ep,10)*bcy)
152 2 +c1*(bf(ep,4)*mm(ep,1)+
153 3 bf(ep,5)*mm(ep,2)+
154 4 bf(ep,6)*mm(ep,3) )
155C I1=3 I2=3
156 vf(ep,3,1)=vf(ep,3,1)+c2*(bm(ep,7)*ff(ep,1)+
157 2 bm(ep,8)*ff(ep,2)+
158 4 bc(ep,5)*bcx+
159 5 bc(ep,6)*bcy )
160 6 +c1*(bmf(ep,7)*mm(ep,1)+
161 7 bmf(ep,8)*mm(ep,2)+
162 8 bmf(ep,9)*mm(ep,3))
163C J=2 K0=2 K1=3 K2=5 I1=4 I2=6
164 vf(ep,1,2)=vf(ep,1,2)+c2*(bm(ep,10)*ff(ep,1)+
165 2 bm(ep,11)*ff(ep,2)+
166 4 bc(ep,11)*bcx+
167 5 bc(ep,12)*bcy )
168 6 +c1*(bmf(ep,10)*mm(ep,1)+
169 7 bmf(ep,11)*mm(ep,2)+
170 8 bmf(ep,12)*mm(ep,3))
171C I1=3 I2=9
172 vm(ep,1,2)=vm(ep,1,2)+c2*(bc(ep,17)*bcx+
173 1 bc(ep,18)*bcy)
174 2 +c1*(bf(ep,7)*mm(ep,1)+
175 3 bf(ep,8)*mm(ep,2)+
176 4 bf(ep,9)*mm(ep,3) )
177C I1=5 I2=7
178 vf(ep,2,2)=vf(ep,2,2)+c2*(bm(ep,13)*ff(ep,1)+
179 2 bm(ep,14)*ff(ep,2)+
180 4 bc(ep,13)*bcx+
181 5 bc(ep,14)*bcy )
182 6 +c1*(bmf(ep,13)*mm(ep,1)+
183 7 bmf(ep,14)*mm(ep,2)+
184 8 bmf(ep,15)*mm(ep,3))
185C I1=4 I2=10
186 vm(ep,2,2)=vm(ep,2,2)+c2*(bc(ep,19)*bcx+
187 1 bc(ep,20)*bcy)
188 2 +c1*(bf(ep,10)*mm(ep,1)+
189 3 bf(ep,11)*mm(ep,2)+
190 4 bf(ep,12)*mm(ep,3) )
191C I1=6 I2=8
192 vf(ep,3,2)=vf(ep,3,2)+c2*(bm(ep,16)*ff(ep,1)+
193 2 bm(ep,17)*ff(ep,2)+
194 4 bc(ep,15)*bcx+
195 5 bc(ep,16)*bcy )
196 6 +c1*(bmf(ep,16)*mm(ep,1)+
197 7 bmf(ep,17)*mm(ep,2)+
198 8 bmf(ep,18)*mm(ep,3))
199C J=3 K0=4 K1=6 K2=10 I1=7 I2=11
200 vf(ep,1,3)=vf(ep,1,3)+c2*(bm(ep,19)*ff(ep,1)+
201 2 bm(ep,20)*ff(ep,2)+
202 4 bc(ep,21)*bcx+
203 5 bc(ep,22)*bcy )
204 6 +c1*(bmf(ep,19)*mm(ep,1)+
205 7 bmf(ep,20)*mm(ep,2)+
206 8 bmf(ep,21)*mm(ep,3))
207C I1=5 I2=14
208 vm(ep,1,3)=vm(ep,1,3)+c2*(bc(ep,27)*bcx+
209 1 bc(ep,28)*bcy)
210 2 +c1*(bf(ep,13)*mm(ep,1)+
211 3 bf(ep,14)*mm(ep,2)+
212 4 bf(ep,15)*mm(ep,3) )
213C I1=8 I2=12
214 vf(ep,2,3)=vf(ep,2,3)+c2*(bm(ep,22)*ff(ep,1)+
215 2 bm(ep,23)*ff(ep,2)+
216 4 bc(ep,23)*bcx+
217 5 bc(ep,24)*bcy )
218 6 +c1*(bmf(ep,22)*mm(ep,1)+
219 7 bmf(ep,23)*mm(ep,2)+
220 8 bmf(ep,24)*mm(ep,3))
221C I1=6 I2=15
222 vm(ep,2,3)=vm(ep,2,3)+c2*(bc(ep,29)*bcx+
223 1 bc(ep,30)*bcy)
224 2 +c1*(bf(ep,16)*mm(ep,1)+
225 3 bf(ep,17)*mm(ep,2)+
226 4 bf(ep,18)*mm(ep,3) )
227C I1=9 I2=13
228 vf(ep,3,3)=vf(ep,3,3)+c2*(bm(ep,25)*ff(ep,1)+
229 2 bm(ep,26)*ff(ep,2)+
230 4 bc(ep,25)*bcx+
231 5 bc(ep,26)*bcy )
232 6 +c1*(bmf(ep,25)*mm(ep,1)+
233 7 bmf(ep,26)*mm(ep,2)+
234 8 bmf(ep,27)*mm(ep,3))
235C J=4 K0=6 K1=9 K2=15 I=1 I1=10 I2=16
236 vf(ep,1,4)=vf(ep,1,4)+c2*(bm(ep,28)*ff(ep,1)+
237 2 bm(ep,29)*ff(ep,2)+
238 4 bc(ep,31)*bcx+
239 5 bc(ep,32)*bcy )
240 6 +c1*(bmf(ep,28)*mm(ep,1)+
241 7 bmf(ep,29)*mm(ep,2)+
242 8 bmf(ep,30)*mm(ep,3))
243C I1=7 I2=19
244 vm(ep,1,4)=vm(ep,1,4)+c2*(bc(ep,37)*bcx+
245 1 bc(ep,38)*bcy)
246 2 +c1*(bf(ep,19)*mm(ep,1)+
247 3 bf(ep,20)*mm(ep,2)+
248 4 bf(ep,21)*mm(ep,3) )
249C I1=11 I2=17
250 vf(ep,2,4)=vf(ep,2,4)+c2*(bm(ep,31)*ff(ep,1)+
251 2 bm(ep,32)*ff(ep,2)+
252 4 bc(ep,33)*bcx+
253 5 bc(ep,34)*bcy )
254 6 +c1*(bmf(ep,31)*mm(ep,1)+
255 7 bmf(ep,32)*mm(ep,2)+
256 8 bmf(ep,33)*mm(ep,3))
257C I1=8 I2=20
258 vm(ep,2,4)=vm(ep,2,4)+c2*(bc(ep,39)*bcx +
259 1 bc(ep,40)*bcy )
260 2 +c1*(bf(ep,22)*mm(ep,1) +
261 3 bf(ep,23)*mm(ep,2) +
262 4 bf(ep,24)*mm(ep,3) )
263C I1=12 I2=18
264 vf(ep,3,4)=vf(ep,3,4)+c2*(bm(ep,34)*ff(ep,1)+
265 2 bm(ep,35)*ff(ep,2)+
266 4 bc(ep,35)*bcx+
267 5 bc(ep,36)*bcy )
268 6 +c1*(bmf(ep,34)*mm(ep,1)+
269 7 bmf(ep,35)*mm(ep,2)+
270 8 bmf(ep,36)*mm(ep,3))
271C
272 ENDDO
273C
274 RETURN
275 END
276!||====================================================================
277!|| cbaforct ../engine/source/elements/shell/coqueba/cbafori.F
278!||--- called by ------------------------------------------------------
279!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
280!||====================================================================
281 SUBROUTINE cbaforct(JFT ,JLT ,VOL ,X13 ,X24 ,
282 2 Y13 ,Y24 ,VSTRE ,VF ,NPLAT,
283 3 IPLAT ,OFF ,NEL )
284C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
285C CALCUL 'membrane shear traitement' use only PARTIE CONSTANTE
286C VF(3,NODE),VM(2,NODE) : FORCES INTERNES GENERALISEES
287C ENTREES : THK,X13,X24,Y13,Y24,
288C SORTIES : VF,VM
289C-----------------------------------------------
290C I M P L I C I T T Y P E S
291C-----------------------------------------------
292#include "implicit_f.inc"
293C-----------------------------------------------
294C G l o b a l P a r a m e t e r s
295C-----------------------------------------------
296#include "mvsiz_p.inc"
297C-----------------------------------------------
298C D U M M Y A R G U M E N T S
299C-----------------------------------------------
300 INTEGER JFT ,JLT,NPLAT,IPLAT(*),NEL
301 my_real
302 . VOL(*), OFF(*), X13(*) ,X24(*) ,Y13(*) ,
303 . VSTRE(NEL,5),VF(MVSIZ,3,4) ,Y24(*)
304C-----------------------------------------------
305C L O C A L V A R I A B L E S
306C-----------------------------------------------
307 INTEGER I ,J,EP
308 my_real SX1, SX2, SY1, SY2, THOFF
309C
310#include "vectorize.inc"
311 DO EP=jft,nplat
312 i=iplat(ep)
313 thoff =vol(i)*vstre(i,3)*off(i)
314 sx1=-thoff*x24(i)
315 sy1= thoff*y24(i)
316C
317 sx2= thoff*x13(i)
318 sy2= -thoff*y13(i)
319C-------------------------------------------------
320C PARTIE ANTISYM STOCKEE DANS NOEUDS 1,2
321C PARTIE SYM NOEUDS STOCKEE DANS NOEUDS 3,4
322C-------------------------------------------------
323 vf(i,1,1)=vf(i,1,1)+sx1
324 vf(i,2,1)=vf(i,2,1)+sy1
325C
326 vf(i,1,2)=vf(i,1,2)+sx2
327 vf(i,2,2)=vf(i,2,2)+sy2
328 ENDDO
329#include "vectorize.inc"
330 DO ep=nplat+1, jlt
331 i=iplat(ep)
332 thoff =vol(i)*vstre(i,3)*off(i)
333 sx1=-thoff*x24(i)
334 sy1= thoff*y24(i)
335 sx2= thoff*x13(i)
336 sy2= -thoff*y13(i)
337 vf(i,1,1)=vf(i,1,1)+sx1
338 vf(i,2,1)=vf(i,2,1)+sy1
339C
340 vf(i,1,2)=vf(i,1,2)+sx2
341 vf(i,2,2)=vf(i,2,2)+sy2
342C
343 vf(i,1,3)=vf(i,1,3)-sx1
344 vf(i,2,3)=vf(i,2,3)-sy1
345C
346 vf(i,1,4)=vf(i,1,4)-sx2
347 vf(i,2,4)=vf(i,2,4)-sy2
348 ENDDO
349C
350 RETURN
351 END
352!||====================================================================
353!|| cbaforrz ../engine/source/elements/shell/coqueba/cbafori.F
354!||--- called by ------------------------------------------------------
355!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
356!||====================================================================
357 SUBROUTINE cbaforrz(JFT ,JLT ,VOL ,FF ,VSRZ ,
358 2 VF ,VMZ ,BM ,BMRZ ,BRZ ,
359 3 KRZ ,VRLZ ,EINT ,OFF ,DT1C ,
360 4 NPLAT,IPLAT,NG ,NEL )
361C-----------------------------------------------
362C I M P L I C I T T Y P E S
363C-----------------------------------------------
364#include "implicit_f.inc"
365#include "mvsiz_p.inc"
366C-----------------------------------------------
367C D U M M Y A R G U M E N T S
368C-----------------------------------------------
369 INTEGER JFT ,JLT,NPLAT,IPLAT(*),NG,NEL
370 my_real
371 . VOL(*) ,FF(NEL,5),VSRZ(NEL,5) ,VF(MVSIZ,3,4) ,VMZ(MVSIZ,4),
372 . BM(MVSIZ,36),BRZ(MVSIZ,4,4),BMRZ(MVSIZ,3,4),KRZ(*),VRLZ(*),
373 . EINT(NEL,2),OFF(*),DT1C
374C-----------------------------------------------
375C L O C A L V A R I A B L E S
376C-----------------------------------------------
377 INTEGER I ,J,EP
378 my_real
379 . s1, c2, a05,erz,dev,c3,fx3,fy3,csrz
380C -------- sig_rz added here----a verifier EINT(QEPH sans A)-----
381#include "vectorize.inc"
382 DO i=jft,jlt
383 ep =iplat(i)
384 c2=vol(ep)
385 erz= vrlz(i)*dt1c
386 dev = erz*c2*half*off(ep)
387 eint(ep,1) = eint(ep,1)+ vsrz(ep,ng)*dev
388 vsrz(ep,ng)= vsrz(ep,ng)+krz(ep)*erz*off(ep)
389 eint(ep,1) = eint(ep,1)+ vsrz(ep,ng)*dev
390 DO j=1,4
391 vmz(i,j)=vmz(i,j)+(bmrz(i,1,j)*ff(ep,1)+
392 . bmrz(i,2,j)*ff(ep,2)+bmrz(i,3,j)*ff(ep,3))*c2
393 ENDDO
394 ENDDO
395C
396#include "vectorize.inc"
397 DO i=jft,jlt
398 ep =iplat(i)
399 c2=vol(ep)*vsrz(ep,ng)*half
400 DO j=1,4
401 vmz(i,j)= vmz(i,j)+brz(i,4,j)*c2
402 ENDDO
403 ENDDO
404#include "vectorize.inc"
405 DO i=jft,nplat
406 ep=iplat(i)
407 c2=vol(ep)
408 c3=c2*ff(ep,3)
409 csrz=c2*vsrz(ep,ng)*half
410C-------------------------------------------------
411C PARTIE ANTISYM STOCKEE DANS NOEUDS 1,2
412C PARTIE SYM NOEUDS STOCKEE DANS NOEUDS 3,4
413C-------------------------------------------------
414C----------NOEUD 1 ------------
415 vf(ep,1,1)=vf(ep,1,1)+(c3-csrz)*bm(ep,5)
416 vf(ep,2,1)=vf(ep,2,1)+(c3+csrz)*bm(ep,1)
417C
418 fx3 = (c3-csrz)*bm(ep,7)
419 fy3 = (c3+csrz)*bm(ep,3)
420 vf(ep,1,3)=vf(ep,1,3)+fx3
421 vf(ep,2,3)=vf(ep,2,3)+fy3
422C----------NOEUD 2 ------------
423 vf(ep,1,2)=vf(ep,1,2)+(c3-csrz)*bm(ep,6)
424 vf(ep,2,2)=vf(ep,2,2)+(c3+csrz)*bm(ep,2)
425C
426 vf(ep,1,4)=vf(ep,1,4)-fx3
427 vf(ep,2,4)=vf(ep,2,4)-fy3
428 ENDDO
429#include "vectorize.inc"
430 DO i=nplat+1,jlt
431 ep=iplat(i)
432 c2=vol(ep)
433 c3=c2*ff(ep,3)
434 csrz=c2*vsrz(ep,ng)*half
435 vf(ep,1,1)=vf(ep,1,1)+c3*bm(ep,3)+csrz*brz(i,1,1)
436 vf(ep,2,1)=vf(ep,2,1)+c3*bm(ep,6)+csrz*brz(i,2,1)
437 vf(ep,3,1)=vf(ep,3,1)+c3*bm(ep,9)+csrz*brz(i,3,1)
438 vf(ep,1,2)=vf(ep,1,2)+c3*bm(ep,12)+csrz*brz(i,1,2)
439 vf(ep,2,2)=vf(ep,2,2)+c3*bm(ep,15)+csrz*brz(i,2,2)
440 vf(ep,3,2)=vf(ep,3,2)+c3*bm(ep,18)+csrz*brz(i,3,2)
441 vf(ep,1,3)=vf(ep,1,3)+c3*bm(ep,21)+csrz*brz(i,1,3)
442 vf(ep,2,3)=vf(ep,2,3)+c3*bm(ep,24)+csrz*brz(i,2,3)
443 vf(ep,3,3)=vf(ep,3,3)+c3*bm(ep,27)+csrz*brz(i,3,3)
444 vf(ep,1,4)=vf(ep,1,4)+c3*bm(ep,30)+csrz*brz(i,1,4)
445 vf(ep,2,4)=vf(ep,2,4)+c3*bm(ep,33)+csrz*brz(i,2,4)
446 vf(ep,3,4)=vf(ep,3,4)+c3*bm(ep,36)+csrz*brz(i,3,4)
447 ENDDO
448C
449 RETURN
450 END
451!||====================================================================
452!|| cbafori1 ../engine/source/elements/shell/coqueba/cbafori.F
453!||--- called by ------------------------------------------------------
454!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
455!||====================================================================
456 SUBROUTINE cbafori1(JFT ,JLT ,FF ,BM ,VF,
457 . NPLAT,IPLAT,VOL ,NEL )
458C--------------------------------------------------------------------------------------------------
459C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
460C CALCUL DES FORCES INTERNES AUX NOEUDS (MEMBRANE ONLY)
461C ENTREES : NEL,MNT,NPT,NP,NG,NIXC,IXC,HZ,OFF,CDET,EPAIS,BUFCON,VQN,BM,BMF,BF,BC
462C SORTIES : VF,VM
463C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
464C--------------------------------------------------------------------------------------------------
465C-----------------------------------------------
466C I M P L I C I T T Y P E S
467C-----------------------------------------------
468#include "implicit_f.inc"
469#include "mvsiz_p.inc"
470C-----------------------------------------------
471C D U M M Y A R G U M E N T S
472C-----------------------------------------------
473 INTEGER JFT,JLT,NPLAT,IPLAT(*),NEL
474 MY_REAL
475 . FF(NEL,5),VOL(*),BM(MVSIZ,36),VF(MVSIZ,3,4)
476C-----------------------------------------------
477C L O C A L V A R I A B L E S
478C-----------------------------------------------
479 INTEGER EP,I,J
480 MY_REAL
481 . C1,C2
482C--------------------------------------------------------------------------------------------------
483#include "vectorize.inc"
484 DO i=jft,nplat
485 ep=iplat(i)
486 c2=vol(ep)
487C----------NOEUD 1 ------------
488 vf(ep,1,1)=vf(ep,1,1)+c2*(bm(ep,1)*ff(ep,1))
489 vf(ep,2,1)=vf(ep,2,1)+c2*(bm(ep,5)*ff(ep,2))
490C
491 vf(ep,1,3)=vf(ep,1,3)+c2*(bm(ep,3)*ff(ep,1))
492 vf(ep,2,3)=vf(ep,2,3)+c2*(bm(ep,7)*ff(ep,2))
493C----------NOEUD 2 ------------
494 vf(ep,1,2)=vf(ep,1,2)+c2*(bm(ep,2)*ff(ep,1))
495 vf(ep,2,2)=vf(ep,2,2)+c2*(bm(ep,6)*ff(ep,2))
496C
497 vf(ep,1,4)=-vf(ep,1,3)
498 vf(ep,2,4)=-vf(ep,2,3)
499C
500 ENDDO
501 RETURN
502 END
cbafori1
subroutine cbafori1(jft, jlt, ff, bm, vf, nplat, iplat, vol, nel)
Definition cbafori.F:458
cbaforct
subroutine cbaforct(jft, jlt, vol, x13, x24, y13, y24, vstre, vf, nplat, iplat, off, nel)
Definition cbafori.F:284
cbaforrz
subroutine cbaforrz(jft, jlt, vol, ff, vsrz, vf, vmz, bm, bmrz, brz, krz, vrlz, eint, off, dt1c, nplat, iplat, ng, nel)
Definition cbafori.F:361
cbafori
subroutine cbafori(jft, jlt, ng, cdet, thk0, th12, ff0, mm0, nel, bm, bmf, bf, bc, tc, vf, vm, nplat, iplat, vol)
Definition cbafori.F:32