OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbafori.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 cbafori (jft, jlt, ng, cdet, thk0, th12, ff0, mm0, nel, bm, bmf, bf, bc, tc, vf, vm, nplat, iplat, vol)
subroutine cbaforct (jft, jlt, vol, x13, x24, y13, y24, vstre, vf, nplat, iplat, off, nel)
subroutine cbaforrz (jft, jlt, vol, ff, vsrz, vf, vmz, bm, bmrz, brz, krz, vrlz, eint, off, dt1c, nplat, iplat, ng, nel)
subroutine cbafori1 (jft, jlt, ff, bm, vf, nplat, iplat, vol, nel)

Function/Subroutine Documentation

◆ cbaforct()

subroutine cbaforct ( integer jft,
integer jlt,
vol,
x13,
x24,
y13,
y24,
vstre,
vf,
integer nplat,
integer, dimension(*) iplat,
off,
integer nel )

Definition at line 281 of file cbafori.F.

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
my_real
#define my_real
Definition cppsort.cpp:32

◆ cbafori()

subroutine cbafori ( integer jft,
integer jlt,
integer ng,
cdet,
thk0,
th12,
ff0,
mm0,
integer nel,
bm,
bmf,
bf,
bc,
tc,
vf,
vm,
integer nplat,
integer, dimension(*) iplat,
vol )

Definition at line 28 of file cbafori.F.

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

◆ cbafori1()

subroutine cbafori1 ( integer jft,
integer jlt,
ff,
bm,
vf,
integer nplat,
integer, dimension(*) iplat,
vol,
integer nel )

Definition at line 456 of file cbafori.F.

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

◆ cbaforrz()

subroutine cbaforrz ( integer jft,
integer jlt,
vol,
ff,
vsrz,
vf,
vmz,
bm,
bmrz,
brz,
krz,
vrlz,
eint,
off,
dt1c,
integer nplat,
integer, dimension(*) iplat,
integer ng,
integer nel )

Definition at line 357 of file cbafori.F.

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