OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
airbag2.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!|| airbagb ../engine/source/airbag/airbag2.F
25!||--- called by ------------------------------------------------------
26!|| monvol0 ../engine/source/airbag/monvol0.F
27!||--- calls -----------------------------------------------------
28!|| get_u_func ../engine/source/user_interface/ufunc.F
29!|| spmd_exch_fr6 ../engine/source/mpi/kinematic_conditions/spmd_exch_fr6.F
30!|| sum_6_float ../engine/source/system/parit.F
31!||--- uses -----------------------------------------------------
32!|| groupdef_mod ../common_source/modules/groupdef_mod.F
33!||====================================================================
34 SUBROUTINE airbagb(
35 1 IVOLU ,ICBAG ,NJET ,IBAGJET ,NVENT ,
36 2 IBAGHOL ,RVOLU ,RVOLUV ,RCBAG ,RBAGJET ,
37 3 RBAGHOL ,FSAV ,NORMAL ,NN ,
38 4 IGRSURF ,PORO ,IVOLUV ,RBAGVJET,
39 4 FR_MV )
40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE groupdef_mod
44C-----------------------------------------------
45C STRUCTURES AIRBAG, INPUT STARTER 4.4
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "param_c.inc"
54#include "units_c.inc"
55#include "com01_c.inc"
56#include "com04_c.inc"
57#include "com08_c.inc"
58#include "task_c.inc"
59C-----------------------------------------------
60C D u m m y A r g u m e n t s
61C-----------------------------------------------
62 INTEGER IVOLU(*),ICBAG(NICBAG,*),NJET,IBAGJET(NIBJET,*),
63 . NVENT,IBAGHOL(NIBHOL,*),
64 . NN,IVOLUV(NIMV,*),FR_MV(*)
65C REAL
66 my_real
67 . rvolu(*), rvoluv(nrvolu,*),rcbag(nrcbag,*),poro(*),
68 . rbagjet(nrbjet,*),rbaghol(nrbhol,*),fsav(*),normal(3,*),rbagvjet(*)
69 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
70C-----------------------------------------------
71C L o c a l V a r i a b l e s
72C-----------------------------------------------
73 INTEGER I,II, NAV, K, IDEF, KK, IPVENT, NNC, IAD,
74 . iport,iporp,ipora,iport1,iporp1,ipora1,ivdp,
75 . ij ,iv, radvois,pmain,
76 . in1,in2,in3,in4,ittf,idtpdef
77C REAL
78 my_real
79 . gama, cv, cp, pext, pdef, dtpdefi, dtpdefc, tvent, tstope,
80 . apvent, avent, bvent,
81 . amtot, p, ro, vol, hspec,
82 . gmtot, cpa, cpb, cpc, gmi, cpai, cpbi, cpci, tbag,
83 . u, deout, dmout, area, pcrit, pvois, tvois, aa, veps,
84 . aout, aout1, aoutot, flout, de, vvois,
85 . dgeout, dgmout, rnm, rmwi, rnmi, rmwg, rnmg,
86 . deri, temp(4),aisent,achemk,fchemk,vmax,
87 . fport,fporp,fpora,fport1,fporp1,fpora1,scalt,scalp,scals,
88 . fvdp, roex, uisent,tt1,
89 . f1(nn), f2(nn),ttf
90 my_real get_u_func
91 EXTERNAL get_u_func
92 DOUBLE PRECISION
93 . FRMV6(2,6), FRMV6B(6)
94C-----------------------------------------------
95 pmain = fr_mv(nspmd+2)
96 nav = ivolu(3)
97 ittf = ivolu(17)
98 ttf =rvolu(60)
99 gama =rvolu(1)
100 pext =rvolu(3)
101 vol =rvolu(16)
102 veps =rvolu(17)
103 tbag =rvolu(13)
104 amtot =rvolu(20)
105 p =rvolu(12)
106 area =rvolu(18)
107 scalt =rvolu(26)
108 scalp =rvolu(27)
109 scals =rvolu(28)
110C
111 ro = amtot/vol
112 pcrit = p*(two/(gama+one))**(gama/(gama-one))
113C
114C contribution du gaz a l'initial.
115 cpai=rvolu(7)
116 cpbi=rvolu(8)
117 cpci=rvolu(9)
118 gmi =rvolu(11)
119 hspec=gmi*tbag*(cpai+half*cpbi*tbag+third*cpci*tbag*tbag)
120C
121 DO ij=1,njet
122 cpa = rbagjet(2,ij)
123 cpb = rbagjet(3,ij)
124 cpc = rbagjet(4,ij)
125 gmtot= rbagjet(8,ij)
126 hspec=hspec+gmtot*tbag*(cpa+half*cpb*tbag+ third*cpc*tbag*tbag)
127 ENDDO
128 hspec=hspec/max(em20,amtot)
129C--------------------------------
130C FLUX SORTANT PAR LES TROUS
131C--------------------------------
132 aisent=zero
133 achemk=zero
134 fchemk=zero
135 DO iv=1,nvent
136 idef = ibaghol(1,iv)
137 ipvent = ibaghol(2,iv)
138 idtpdef =ibaghol(11,iv)
139C
140 pdef = rbaghol(1,iv)
141 dtpdefi= rbaghol(4,iv)
142 dtpdefc= rbaghol(5,iv)
143 avent = rbaghol(2,iv)
144 tvent = rbaghol(3,iv)
145 bvent = rbaghol(6,iv)
146 tstope = rbaghol(14,iv)
147 rbaghol(16,iv)=zero
148 rbaghol(17,iv)=zero
149 rbaghol(18,iv)=zero
150 rbaghol(21,iv)=zero
151 rbaghol(22,iv)=zero
152C
153 IF(ittf==11.OR.ittf==12.OR.ittf==13) THEN
154 IF(idef==0.AND.p>pdef+pext.
155 . and.dtpdefc>dtpdefi.
156 . and.vol>em3*area**three_half.
157 . and.tt<tstope+ttf
158 . .AND.idtpdef==0) THEN
159 idef=1
160 IF(ispmd+1==pmain) THEN
161 WRITE(iout,*)' *** AIRBAG VENT HOLES MEMBRANE IS DEFLATED ***'
162 WRITE(iout,*)' *** MONITORED VOLUME ',ivolu(1),
163 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
164 WRITE(istdo,*)' *** VENT HOLES MEMBRANE IS DEFLATED ***'
165 ENDIF
166 ENDIF
167 IF(idef==0.AND.dtpdefc>dtpdefi.
168 . and.tt<tstope+ttf
169 . .AND.idtpdef==1) THEN
170 idef=1
171 WRITE(iout,*)' *** AIRBAG VENT HOLES MEMBRANE IS DEFLATED ***'
172 WRITE(iout,*)' *** MONITORED VOLUME ',ivolu(1),
173 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
174 WRITE(istdo,*)' *** VENT HOLES MEMBRANE IS DEFLATED ***'
175 ENDIF
176 IF(idef==0 .AND. tt>tvent+ttf
177 . .AND. tt<tstope+ttf) THEN
178 idef=1
179 IF(ispmd+1==pmain) THEN
180 WRITE(iout,*) ' *** AIRBAG VENTING STARTS ***'
181 WRITE(iout,*) ' *** MONITORED VOLUME ',ivolu(1),
182 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
183 WRITE(istdo,*)' *** VENTING STARTS ***'
184 ENDIF
185 ENDIF
186 IF(idef==1 .AND. tt>=tstope+ttf) THEN
187 idef=0
188 IF(ispmd+1==pmain) THEN
189 WRITE(iout,*) ' *** AIRBAG VENTING STOPS ***'
190 WRITE(iout,*) ' *** MONITORED VOLUME ',ivolu(1),
191 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
192 WRITE(istdo,*)' *** VENTING STOPS ***'
193 END IF
194 END IF
195 ELSE IF(ittf==0) THEN
196 IF(idef==0.AND.p>pdef+pext.
197 . and.dtpdefc>dtpdefi.
198 . and.vol>em3*area**three_half.
199 . and.tt<tstope
200 . .AND.idtpdef==0) THEN
201 idef=1
202 IF(ispmd+1==pmain) THEN
203 WRITE(iout,*)' *** AIRBAG VENT HOLES MEMBRANE IS DEFLATED ***'
204 WRITE(iout,*)' *** MONITORED VOLUME ',ivolu(1),
205 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
206 WRITE(istdo,*)' *** VENT HOLES MEMBRANE IS DEFLATED ***'
207 ENDIF
208 ENDIF
209 IF(idef==0.AND.dtpdefc>dtpdefi.
210 . and.tt<tstope
211 . .AND.idtpdef==1) THEN
212 idef=1
213 WRITE(iout,*)' *** AIRBAG VENT HOLES MEMBRANE IS DEFLATED ***'
214 WRITE(iout,*)' *** MONITORED VOLUME ',ivolu(1),
215 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
216 WRITE(istdo,*)' *** VENT HOLES MEMBRANE IS DEFLATED ***'
217 ENDIF
218 IF(idef==0 .AND. tt>tvent
219 . .AND. tt<tstope) THEN
220 idef=1
221 IF(ispmd+1==pmain) THEN
222 WRITE(iout,*) ' *** AIRBAG VENTING STARTS ***'
223 WRITE(iout,*) ' *** MONITORED VOLUME ',ivolu(1),
224 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
225 WRITE(istdo,*)' *** VENTING STARTS ***'
226 ENDIF
227 ENDIF
228 IF(idef==1 .AND. tt>=tstope) THEN
229 idef=0
230 IF(ispmd+1==pmain) THEN
231 WRITE(iout,*) ' *** AIRBAG VENTING STOPS ***'
232 WRITE(iout,*) ' *** MONITORED VOLUME ',ivolu(1),
233 . ' VENT HOLES MEMBRANE NUMBER ',iv,' ***'
234 WRITE(istdo,*)' *** VENTING STOPS ***'
235 END IF
236 END IF
237 END IF
238C
239 IF(ipvent/=0)THEN
240 IF(intbag==0)THEN
241 nnc=igrsurf(ipvent)%NSEG
242 DO kk=1,nnc
243 IF(igrsurf(ipvent)%ELTYP(kk)==3)THEN
244 k=igrsurf(ipvent)%ELEM(kk)
245 ELSEIF(igrsurf(ipvent)%ELTYP(kk)==7)THEN
246 k=igrsurf(ipvent)%ELEM(kk) + numelc
247 ELSE
248 k=igrsurf(ipvent)%ELEM(kk) + numelc + numeltg
249 ENDIF
250 f1(kk) = sqrt( normal(1,k)**2+normal(2,k)**2+normal(3,k)**2 )
251 f2(kk) = zero
252 ENDDO
253 ELSE
254 nnc=igrsurf(ipvent)%NSEG
255 DO kk=1,nnc
256 IF(igrsurf(ipvent)%ELTYP(kk)==3)THEN
257 k=igrsurf(ipvent)%ELEM(kk)
258 ELSEIF(igrsurf(ipvent)%ELTYP(kk)==7)THEN
259 k=igrsurf(ipvent)%ELEM(kk) + numelc
260 ELSE
261 k=igrsurf(ipvent)%ELEM(kk) + numelc + numeltg
262 ENDIF
263 aa = sqrt( normal(1,k)**2+normal(2,k)**2+normal(3,k)**2 )
264 f1(kk) = (one - poro(k))*aa
265 f2(kk) = poro(k)*aa
266 ENDDO
267 ENDIF
268C
269C somme parith/on
270C
271 DO k = 1, 6
272 frmv6(1,k) = zero
273 frmv6(2,k) = zero
274 END DO
275 CALL sum_6_float(1, nnc, f1, frmv6(1,1),2)
276 CALL sum_6_float(1, nnc, f2, frmv6(2,1),2)
277C comm si necessaire
278 IF(nspmd > 1) THEN
279 CALL spmd_exch_fr6(fr_mv,frmv6,2*6)
280 ENDIF
281C
282 aout = frmv6(1,1)+frmv6(1,2)+frmv6(1,3)+
283 . frmv6(1,4)+frmv6(1,5)+frmv6(1,6)
284 aout1 = frmv6(2,1)+frmv6(2,2)+frmv6(2,3)+
285 . frmv6(2,4)+frmv6(2,5)+frmv6(2,6)
286 ELSE
287 aout =avent
288 avent=one
289 aout1=zero
290 ENDIF
291C
292 IF(idef==1 .AND. p>pext.
293 . and.vol>em3*area**three_half)THEN
294 iport =ibaghol(3,iv)
295 iporp =ibaghol(4,iv)
296 ipora =ibaghol(5,iv)
297 iport1=ibaghol(6,iv)
298 iporp1=ibaghol(7,iv)
299 ipora1=ibaghol(8,iv)
300 fport = rbaghol(7,iv)
301 fporp = rbaghol(8,iv)
302 fpora = rbaghol(9,iv)
303 fport1= rbaghol(10,iv)
304 fporp1= rbaghol(11,iv)
305 fpora1= rbaghol(12,iv)
306 IF(ipora/=0.AND.ipvent/=0)THEN
307 aout=fpora*avent*get_u_func(ipora,aout*scals,deri)
308 ELSE
309 aout=avent*aout
310 ENDIF
311 tt1=tt
312 IF (ittf==12) tt1=tt-ttf
313 IF (ittf==13) tt1=tt-ttf-tvent
314 IF(iport/=0)aout=fport*aout*
315 . get_u_func(iport,tt1*scalt,deri)
316 IF(iporp/=0)aout=fporp*aout*
317 . get_u_func(iporp,(p-pext)*scalp,deri)
318 IF(ipora1/=0.AND.ipvent/=0)THEN
319 aout1=fpora*bvent*get_u_func(ipora1,aout1*scals,deri)
320 ELSE
321 aout1=bvent*aout1
322 ENDIF
323 IF(iport1/=0)aout1=fport*aout1*
324 . get_u_func(iport1,tt1*scalt,deri)
325 IF(iporp1/=0)aout1=fporp*aout1*
326 . get_u_func(iporp1,(p-pext)*scalp,deri)
327 IF(ispmd+1==pmain) THEN
328 rbaghol(16,iv)=aout
329 rbaghol(17,iv)=aout1
330 END IF
331 ivdp=ibaghol(9,iv)
332 IF(ivdp==0)THEN
333 aisent=aisent+aout+aout1
334 ELSE
335 achemk=achemk+aout+aout1
336 fvdp=rbaghol(13,iv)
337 u=fvdp*get_u_func(ivdp,(p-pext)*scalp,deri)
338 fchemk= fchemk+(aout+aout1)*u
339 IF(ispmd+1==pmain) rbaghol(18,iv)=u
340 ENDIF
341 ENDIF
342 ibaghol(1,iv)=idef
343 ENDDO
344C-------
345 aoutot=aisent+achemk
346 uisent=zero
347 flout =zero
348 dmout =zero
349 IF(aoutot>zero)THEN
350 roex =ro*(pext/p)**(one/gama)
351 vmax =half*(p-pext)*vol/(gama-one)
352 . /max(em20,hspec*(roex*aisent+ro*achemk)*dt1)
353 vmax =min(vmax,half*vol/max(em20,aoutot*dt1))
354C
355 IF(aisent>zero)THEN
356 pext = max(pext,pcrit)
357 roex =ro*(pext/p)**(one/gama)
358 uisent=two*gama/(gama-one)*p/ro*(one-(pext/p)**((gama-one)/gama))
359 uisent=max(uisent,zero)
360 uisent=sqrt(uisent)
361 uisent=min(uisent,vmax)
362 flout=aisent*uisent
363 dmout=flout*roex
364 ENDIF
365C
366 IF(achemk>zero)THEN
367 fchemk=min(fchemk,vmax*achemk)
368 flout =flout +fchemk
369 dmout =dmout +ro*fchemk
370 ENDIF
371 IF(ispmd+1==pmain)THEN
372 DO iv=1,nvent
373 idef=ibaghol(1,iv)
374 ivdp=ibaghol(9,iv)
375 IF(idef==1)THEN
376 IF(ivdp==0)THEN
377 rbaghol(18,iv)=uisent
378 rbaghol(21,iv)= roex*uisent
379 . *(rbaghol(16,iv)+rbaghol(17,iv))
380 rbaghol(22,iv)=rbaghol(21,iv)*hspec
381 ELSE
382 rbaghol(18,iv)=min(rbaghol(18,iv),vmax)
383 rbaghol(21,iv)= ro*rbaghol(18,iv)
384 . *(rbaghol(16,iv)+rbaghol(17,iv))
385 rbaghol(22,iv)=rbaghol(21,iv)*hspec
386 END IF
387 END IF
388 END DO
389 END IF
390 ENDIF
391C---------------------------------------------
392 rnm =rvolu(14)
393 cv =rnm/amtot/(gama-one)
394 cp =gama*cv
395 IF(ispmd+1==pmain) THEN
396 fsav(1) =amtot
397 fsav(2) =vol
398 fsav(3) =p
399 fsav(4) =area
400 fsav(5) =tbag
401 fsav(6) =aoutot
402 fsav(7) =flout/max(em20,aoutot)
403 fsav(8)=zero
404 fsav(9)=zero
405 fsav(10)=cp
406 fsav(11)=cv
407 fsav(12)=gama
408 fsav(15)=zero
409 DO ij=1,njet
410 fsav(15)=fsav(15)+rbagjet(7,ij)
411 ENDDO
412 ENDIF
413C---------------------------------------------
414C MASSE et TRAVAIL par GAZ
415C---------------------------------------------
416 rmwi=rvolu(10)
417 rnmi=gmi*rmwi
418C VOLG/VOL=fraction molaire=RNMG/RNM
419 dgmout=rnmi/max(em20,rnm)*dmout
420 dgeout=dgmout*tbag*(cpai+half*cpbi*tbag+third*cpci*tbag*tbag)
421 rvolu(22)=rvolu(22) + dgeout
422 rvolu(24)=rvolu(24) + dgmout
423C------
424 DO ij=1,njet
425 cpa = rbagjet(2,ij)
426 cpb = rbagjet(3,ij)
427 cpc = rbagjet(4,ij)
428 gmtot =rbagjet(8,ij)
429 rmwg =rbagjet(1,ij)
430 rnmg =gmtot*rmwg
431 dgmout=rnmg/max(em20,rnm)*dmout
432 dgeout=dgmout*tbag*(cpa+half*cpb*tbag+third*cpc*tbag*tbag)
433 rbagjet( 9,ij)=rbagjet( 9,ij)+dgmout
434 rbagjet(10,ij)=rbagjet(10,ij)+dgeout
435 ENDDO
436C---------------------------------------------
437C AIRBAG COMMUNIQUANTS
438C---------------------------------------------
439 DO i=1,nav
440 ii = icbag(1,i)
441 ipvent = icbag(2,i)
442 idef = icbag(3,i)
443 iport = icbag(4,i)
444 iporp = icbag(5,i)
445 pdef = rcbag(1,i)
446 avent = rcbag(2,i)
447 tvent = rcbag(3,i)
448 dtpdefi= rcbag(4,i)
449 dtpdefc= rcbag(5,i)
450 fport = rcbag(6,i)
451 fporp = rcbag(7,i)
452 pvois=rvoluv(12,ii)
453 vvois=rvoluv(16,ii)
454 IF(ittf==0.OR.ittf==11.OR.ittf==12.OR.ittf==13)THEN
455 IF(idef==0.AND.p>pdef+pvois
456 . .AND.dtpdefc>dtpdefi
457 . .AND.vol>em3*area**three_half)THEN
458 idef=1
459 IF(ispmd+1==pmain) THEN
460 WRITE(iout,*)
461 . ' *** CHAMBER VENT HOLES MEMBRANE IS DEFLATED ***'
462 WRITE(iout,*)
463 . ' *** MONITORED VOLUME ',ivolu(1),' ***'
464 WRITE(istdo,*)
465 . ' *** CHAMBER VENT HOLES MEMBRANE IS DEFLATED ***'
466 ENDIF
467 ENDIF
468 IF(idef==0 .AND. tt>tvent+ttf) THEN
469 idef=1
470 IF(ispmd+1==pmain) THEN
471 WRITE(iout,*) ' *** CHAMBER COMMUNICATION STARTS ***'
472 WRITE(iout,*) ' *** MONITORED VOLUME ',ivolu(1),' ***'
473 WRITE(istdo,*)' *** COMMUNICATION STARTS ***'
474 ENDIF
475 ENDIF
476 ENDIF
477C
478 IF(ipvent/=0)THEN
479 nnc=igrsurf(ipvent)%NSEG
480 DO kk=1,nnc
481 IF(igrsurf(ipvent)%ELTYP(kk)==3)THEN
482 k=igrsurf(ipvent)%ELEM(kk)
483 ELSEIF(igrsurf(ipvent)%ELTYP(kk)==7)THEN
484 k=igrsurf(ipvent)%ELEM(kk) + numelc
485 ELSE
486 k=igrsurf(ipvent)%ELEM(kk) + numelc + numeltg
487 ENDIF
488 f1(kk) = sqrt( normal(1,k)**2+normal(2,k)**2+normal(3,k)**2 )
489 ENDDO
490C
491C Sommation p/on
492C
493 DO k = 1, 6
494 frmv6b(k) = zero
495 END DO
496 CALL sum_6_float(1, nnc, f1, frmv6b, 1)
497C comm si necessaire
498 IF(nspmd > 1) THEN
499 CALL spmd_exch_fr6(fr_mv,frmv6b,6)
500 ENDIF
501 apvent = frmv6b(1)+frmv6b(2)+frmv6b(3)+
502 . frmv6b(4)+frmv6b(5)+frmv6b(6)
503 ELSE
504 apvent = one
505 ENDIF
506C
507 aout=avent*apvent
508C
509 IF(idef==1 .AND. p>pvois.
510 . and.vol>em3*area**three_half)THEN
511 pvois = max(pvois,pcrit)
512 u=two*gama/(gama-one)*p/ro*(one-(pvois/p)**((gama-one)/gama))
513 u=sqrt(u)
514 u=min(u,half*vol/max(em20,aout*dt1))
515 de=ro*(pvois/p)**(one/gama)*hspec
516 u=min(u,(p-pvois)*half*min(vol,vvois)
517 . /(gama-one)/de/max(em20,aout*dt1))
518 flout=aout*u
519 dmout=flout*ro*(pvois/p)**(one/gama)
520 ELSE
521 dmout=zero
522 flout=zero
523 ENDIF
524 icbag(3,i) = idef
525C---------------------------------------------
526C MASSE et TRAVAIL par GAZ
527C---------------------------------------------
528C VOLG/VOL=fraction molaire=RNMG/RNM
529 dgmout=rnmi/max(em20,rnm)*dmout
530 dgeout=dgmout*tbag*(cpai+half*cpbi*tbag+third*cpci*tbag*tbag)
531C OUT
532 rvolu(22)=rvolu(22) + dgeout
533 rvolu(24)=rvolu(24) + dgmout
534C IN
535 rvoluv(22,ii)=rvoluv(22,ii) - dgeout
536 rvoluv(24,ii)=rvoluv(24,ii) - dgmout
537C------
538 radvois= ivoluv(10,ii)
539 DO ij=1,njet
540 gmtot =rbagjet(8,ij)
541 rmwg =rbagjet(1,ij)
542 rnmg =gmtot*rmwg
543 dgmout=rnmg/max(em20,rnm)*dmout
544 cpa = rbagjet(2,ij)
545 cpb = rbagjet(3,ij)
546 cpc = rbagjet(4,ij)
547 dgeout=dgmout*tbag*(cpa+half*cpb*tbag+third*cpc*tbag*tbag)
548C OUT
549 rbagjet(9,ij)=
550 . rbagjet(9,ij)+dgmout
551 rbagjet(10,ij)=
552 . rbagjet(10,ij)+dgeout
553C IN
554 rbagvjet(radvois+nrbjet*(ij-1)+9)=
555 . rbagvjet(radvois+nrbjet*(ij-1)+9)-dgmout
556 rbagvjet(radvois+nrbjet*(ij-1)+10)=
557 . rbagvjet(radvois+nrbjet*(ij-1)+10)-dgeout
558C
559 ENDDO
560 IF(ispmd+1==pmain) THEN
561 fsav(8)=fsav(8)+aout
562 fsav(9)=fsav(9)+flout
563 ENDIF
564 ENDDO
565 IF(ispmd+1==pmain) THEN
566 fsav(9)=fsav(9)/max(em20,fsav(8))
567 ENDIF
568C---------------------------------------------
569 RETURN
570 END
airbagb
subroutine airbagb(ivolu, icbag, njet, ibagjet, nvent, ibaghol, rvolu, rvoluv, rcbag, rbagjet, rbaghol, fsav, normal, nn, igrsurf, poro, ivoluv, rbagvjet, fr_mv)
Definition airbag2.F:40
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:662
sum_6_float
subroutine sum_6_float(jft, jlt, f, f6, n)
Definition parit.F:64
poro
subroutine poro(geo, nodpor, ms, x, v, w, af, am, skew, weight, nporgeo)
Definition poro.F:40
spmd_exch_fr6
subroutine spmd_exch_fr6(fr, fs6, len)
Definition spmd_exch_fr6.F:51