Simulation of the Effects of an Air Blast Wave - ELSA - Europa

More documents

Recommendations

Info

l1 = p2 d p3 d p4 d p1;l3 = p1 d p2 d p3 d p4 d p1;a1 = surf l3 plane;p10 = 0 0 (sizex);*a2 = l1 tran p10;a3 = a1 plus p10;elim (a2 et a3);l2 = p1 d p2;a4 = l2 tran p10;elim (a4 et a2);p11 = 0 0 (0-sizex);a1 = orie a1 p11;*elim (a1 et a2 et a3 et a4);geom1= (a1 et a2 et a3 et a4) coul roug;opti elem cub8;v1 = (geom1) volu;* Erstellen der Huelle fuer IMPEgeom2 = (geom1 et v1);*Explosiveep1 = p1;ep2 = (size) 0 0;ep3 = (size) (size) 0;ep4 = 0 (size) 0;ep5 = 0 0 (size);el1 = ep1 d ep2 d ep3 d ep4 d;ea1 = surf el1 plane;ve1 = ea1 volu (size/deel) tran ep5;elim 0.001 (v1 et ve1);air1=diff v1 ve1;exp1=ve1;vges= (air1 et ve1);nxpl = pxpdroi1 vges p1 (p1 plus (sizex 0 0)) 0.0005;nxpl = pxordpoi nxpl p1;nxp2 = pxpdroi1 vges p1 (p1 plus (sizex sizex sizex)) 0.0005;nxp2 = pxordpoi nxp2 p1;geom_new = (air1 et geom1 et exp1 et nxpl et nxp2);TASS geom_new;OPTI sauv form 'cubX.msh';sauv form geom_new;list (nbno geom1);list(mesu(exp1));con.epsCUBX - cubical model$ECHOCONV WINCAST geom_newTRID NONL ALEOPTI NF34OPTI TOLC 1e-1$DIMEPT6L X ! CUB1:100000, CUB2, CUB3, CUB4:300000,CUB4:FL38 840000 FL35 1 ZONE 4ECROU 2200000NBLO 100000 NALE 10000NBLE X ! CUB1:20000, CUB2, CUB3:150000,CUB4:230000TERM$GEOMFL38 air1FL38 exp1TERM*GRIL EULE LECT exp1 TERMALE TOUSAUTO AUTR$MATE$ airflut RO 1.3 EINT 2.1978E5 GAMM 1.35 PB 0ITER 1 ALF0 1 BET0 1 KINT 0 AHGF 0 CL 0.5CQ 2.56 PMIN 0 PREF 1.e5 NUM 11a 3.738e11 b 3.747e9 r1 4.15 r2 0.90ros 1630LECT air1 TERM$ explosiveflut ro 1630 eint 3.68e6 gamm 1.35 PB 0ITER 1 ALF0 1 BET0 1 KINT 0 AHGF 0 CL 0.5CQ 2.56 PMIN 0 PREF 1.e5 NUM 11a 3.738e11 b 3.747e9 r1 4.15 r2 0.90d 6930 TDET 0.0 pini 1e5xdet 0.0 ydet 0.0 zdet 0.0LECT exp1 TERM$LINK COUPFSR LECT geom1 TERMECRI DEPL VITE CONT ECRO TFRE 10.E-3FICH ALIC TFRE 2E-6$OPTI NOTE LOG 1$CALC TINI 0 TEND 4.75e-4finexplosive.dgibi*Conical Model for calculations only with the explosiveopti donn 'D:\Users\larchma\cast3m\pxpdroi1.procedur';opti donn 'D:\Users\larchma\cast3m\pxordpoi.procedur';OPTI echo 1;OPTI dime 3 elem cub8;dex=0.124; ! length of the cone**********************************************Parameter *******************************************************************************opex=2e-2; ! openingdexin= 0.001; ! element size in the centerdexfi= 0.01; ! element size at the enddpy=0.001; ! Length of the pyramid element*********************************************DENS 100;oppy=dpy*opex/dex;p1 = 0 0 0;*Points of the pyramidppy2 = (dpy) (oppy) (oppy);ppy3 = (dpy) (0-oppy) (oppy);ppy4 = (dpy) (0-oppy) (0-oppy);ppy5 = (dpy) (oppy) (0-oppy);lpyr1 = p1 d ppy2;80
lpyr2 = p1 d ppy3;lpyr3 = p1 d ppy4;lpyr4 = p1 d ppy5;pyra = manu pyr5 p1 ppy5 ppy4 ppy3 ppy2;*Points of the explosivepe2 = (dex) (opex) (opex);pe3 = (dex) (0-opex) (opex);pe4 = (dex) (0-opex) (0-opex);pe5 = (dex) (opex) (0-opex);*defining the lines between the pyramid and the explosivelpe1 = ppy2 d ppy3;lpe2 = ppy3 d ppy4;lpe3 = ppy4 d ppy5;lpe4 = ppy5 d ppy2;*defining the lines around explosiveDENS dexin;le1 = ppy2 d 'DINI' dexin 'DFIN' dexfi pe2;le2 = ppy3 d 'DINI' dexin 'DFIN' dexfi pe3;le3 = ppy4 d 'DINI' dexin 'DFIN' dexfi pe4;le4 = ppy5 d 'DINI' dexin 'DFIN' dexfi pe5;*defining the lines between the explosive and the airDENS 100.;lea1 = pe2 d pe3;lea2 = pe3 d pe4;lea3 = pe4 d pe5;lea4 = pe5 d pe2;*defining the surfaces around the explosiveae1=dall lpe1 le1 lea1 le2;ae2=dall lpe2 le2 lea2 le3;ae3=dall lpe3 le3 lea3 le4;ae4=dall lpe4 le4 lea4 le1;ae5=dall lpe1 lpe2 lpe3 lpe4;ae6=dall lea1 lea2 lea3 lea4;*defining the surfaces around the pyramideapyr1 = surf (lpyr1 et lpyr2 et lpe1) plane;apyr2 = surf (lpyr2 et lpyr3 et lpe2) plane;apyr3 = surf (lpyr3 et lpyr4 et lpe3) plane;apyr4 = surf (lpyr4 et lpyr1 et lpe4) plane;apyrsum = (apyr1 et apyr2 et apyr3 et apyr4 et ae6);elim 1e-10 apyrsum;*defining the volume of the explosivegeomex= (ae1 et ae2 et ae3 et ae4 et ae5 et ae6) coul roug;elim 1e-10 (geomex);vex = (geomex) volu;vex1=vex;*Points for the controllfp1 = vex1 poin proche (0.1 0 0);pp0 = fp1 et fp1;REPE I0 (NBEL vex1);TEST0 = pp0 INCL (vex1 ELEM CUB8 &I0) 'VOLU';SI ((NBEL TEST0)> 0);quit I0;FINS;FIN I0;MESS &I0;fe1 = vex1 elem CUB8 &I0;*lines for the output of the pressurevges = (vex1);nxpl = pxpdroi1 vges p1 (p1 plus pe2) 2e-5;nxpl = pxordpoi nxpl p1;*areas on the sidesasum1 = (apyr1 et ae1);asum2 = (apyr2 et ae2);asum3 = (apyr3 et ae3);asum4 = (apyr4 et ae4);asum = (asum1 et asum2 et asum3 et asum4);geom_new = (vex1 et pyra et nxpl et fe1 et asum);TASS geom_new;OPTI sauv form 'explosiveX.msh';sauv form geom_new;cv8.dgibi* Construction d'une sphere a partir d'un cube* control volume modelopti donn 'D:\Users\larchma\cast3m\pxpdroi1.procedur';opti donn 'D:\Users\larchma\cast3m\pxordpoi.procedur';opti dime 3 elem cub8;*Nombre de bissectionsnel0 = 30;sizeex = 1.0;sizeai = 4.0;*Cote du cube intermediairer0 = .25;dini = 3.141*sizeex/(4.*nel0);dfin = 3.141*sizeai/(4.*nel0);*Referenceo0 = 0. 0. 0.;x0 = (sizeex) 0. 0.;xa0 = 0 (sizeex) (sizeex);xb0 = (sizeex) (sizeex) (sizeex);xc0 = (sizeex) 0 (sizeex);xd0 = 0 0 (sizeex);y0 = 0. (sizeex) 0.;z0 = 0. 0. (sizeex);x1 = (sizeai) 0. 0.;y1 = 0. (sizeai) 0.;z1 = 0. 0. (sizeai);c0 = x0 plus y0 plus z0 / 2.;c1 = x1 plus y1 plus z1 / 2.;symp1 = (sizeex/2.) (sizeex/2.) (sizeex/2.);symp2 = (sizeai/2.) (sizeai/2.) (sizeai/2.);*Cube intermediaire (centre=o0 et arete=r0)cub0 = (o0 droi nel0 y0 tran nel0 x0)et (o0 droi nel0 z0 tran nel0 y0) et(o0 droi nel0 x0 tran nel0 z0)syme 'POINT' symp1 homo o0 r0;cub2 = (o0 droi nel0 y1 tran nel0 x1)et (o0 droi nel0 z1 tran nel0 y1) et(o0 droi nel0 x1 tran nel0 z1)syme 'POINT' symp2;*Pojection sur la sphere de rayon unitairespe1 = cub0 proj 'CONI' o0 'SPHE' o0 x0;spe2 = cub2 proj 'CONI' o0 'SPHE' o0 x1;a_press = spe1 coul 'BLEU';*Remplissageair1 = spe1 volu 'DINI' dini 'DFIN' dfin spe2 coul bleu;ages = enve air1;a_abso = ages diff a_press;*Points for the controllnxp1 = pxpdroi1 air1 x0 (x0 plus x1) 0.0005;nxp1 = pxordpoi nxp1 x0;81
Page 1 and 2:
Simulation of the Effectsof an Air
Page 3 and 4:
Distribution ListLechner S.Anthoine
Page 5 and 6:
1 IntroductionThis work is being co
Page 7 and 8:
• Far zone. The blast wave result
Page 9 and 10:
Figure 2: Model of Kingery [14] wit
Page 11 and 12:
2.2.3 ImpulseThe impulse of the air
Page 13 and 14:
Figure 5: Different parameters for
Page 15 and 16:
250200Impulse from KingeryImpulse w
Page 17 and 18:
the specific heat ratio and the spe
Page 19 and 20:
determined with a fluid pre-calcula
Page 21 and 22:
From this asymptotic form it can al
Page 23 and 24:
2.0E+101.6E+10without burn mass fra
Page 25 and 26:
Case opening d pyra d ex,in d ex,en
Page 27 and 28:
Pressure [Pa]2.4E+102.0E+101.6E+101
Page 29 and 30: 5.0E+094.0E+09EUROPLEXUS x=0.105EUR
Page 31 and 32: Cased pyr /d ex,in /d ex,end /d air
Page 33 and 34: Figure 21: Comparison of a coarse w
Page 35 and 36: CON9This model has approximately th
Page 37 and 38: 4E+7Max. pressure [Pa]3E+72E+7Kinge
Page 39 and 40: Case Size of the model Description
Page 41 and 42: CUB4This model uses a finer mesh th
Page 43 and 44: 6. Complete the volumes by adding t
Page 45 and 46: Figure 33: Maximum pressures versus
Page 47 and 48: experimental values is obtained wit
Page 49 and 50: 4.6.3 Arrival TimeThe arrival time
Page 51 and 52: value in [13]). These results show
Page 53 and 54: Modell airl bBubbleBubbleBubbleTNT
Page 55 and 56: Several calculations are done with
Page 57 and 58: 3. Calculation of the factor α bub
Page 59 and 60: quadratic. Therefore, the size of t
Page 61 and 62: 400Impulse [Pa sec]300200.Kingery 5
Page 63 and 64: 5 = hemispherical detonation (half
Page 65 and 66: 9 References[1] Alia, A.; Souli, M.
Page 67 and 68: 10 Apendix10.1 EUROPLEXUS Codedyms.
Page 69 and 70: CALLMECVAL(GET_FMT(6),NEWMAT%MATVAL
Page 71 and 72: 10.2 Miscellaneous codeairblastresu
Page 73 and 74: else{ //Kingerydouble t=log10(z);do
Page 75 and 76: maxar=0;REPE I0 (NBEL vges);ar = as
Page 77 and 78: *defining the surfaces around the a
Page 79: lpyr1 = p1 d ppy2;lpyr2 = p1 d ppy3
Page 83: *Cube intermediaire (centre=o0 et a
Page 86: The mission of the JRC is to provid
show all

Simulation of the Effects of an Air Blast Wave - ELSA - Europa

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?