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

GRIL EULE LECT vex1 pyra 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 vai1 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 vex1 pyra TERM$LINK COUPFSR LECT asum1 asum2 asum3 asum4 p1 a6 TERMECRI DEPL VITE CONT ECRO TFRE 10.E-3FICH ALIC TFRE 2E-6$OPTI NOTELOG 1$CALC TINI 0 TEND 1.5e-3*============================================SUITPost-treatment (time curves from alice temps file)ECHO*RESU ALIC GARD PSCR*SORT GRAP*AXTE 1.0 'Time [s]'*COUR 1 'press_1' ECRO COMP 1 GAUSS 1 ELEM LECT fe1TERMCOUR 2 'press_2' ECRO COMP 1 GAUSS 1 ELEM LECT fe2TERMCOUR 3 'press_3' ECRO COMP 1 GAUSS 1 ELEM LECT fe3TERMCOUR 4 'press_4' ECRO COMP 1 GAUSS 1 ELEM LECT fe4TERMCOUR 5 'press_5' ECRO COMP 1 GAUSS 1 ELEM LECT fe5TERMCOUR 6 'press_6' ECRO COMP 1 GAUSS 1 ELEM LECT fe6TERMCOUR 7 'int1' subc 1 1e5COUR 8 'int2' subc 2 1e5COUR 11 'int3' INT 7COUR 12 'int4' INT 8COUR 21 'press_1' ECRO COMP 1 GAUSS 1 ELEM LECT 1TERMCOUR 22 'press_2' ECRO COMP 1 GAUSS 1 ELEM LECT 2TERMCOUR 23 'press_3' ECRO COMP 1 GAUSS 1 ELEM LECT 3TERMCOUR 24 'vi_1_160' VITE LECT 160 TERMCOUR 25 'vi_2_160' VITE LECT 160 TERMCOUR 26 'vi_1_156' VITE LECT 156 TERMCOUR 27 'vi_1_155' VITE LECT 155 TERMSCOURBE 33 '1E-4' NSTO 51 SAXE 1 'x' LECT nxpl TERMECRO COMP 1SCOURBE 34 '5E-4' NSTO 251 SAXE 1 'x' LECT nxpl TERMECRO COMP 1trac 1 2 3 4 5 6 TEXT axes 1.0 'PRESS [MPa]'trac 11 12 TEXT axes 1.0 'PRESS*time [MPa]'trac 21 22 23 TEXT axes 1.0 'PRESS [MPa]' xmin 0 xmax 2e-4nx 10trac 24 25 26 27 TEXT axes 1.0 'Velocity [m/sec]'trac 33 34 TEXT axes 1.0 'distance [m]'SCOURBE 101 ' 2.00E-06 ' NSTO 2 SAXE 1 'x' LECT nxplTERM ECRO COMP 1SCOURBE 102 ' 4.00E-06 ' NSTO 3 SAXE 1 'x' LECT nxplTERM ECRO COMP 1** and so on****************************************SCOURBE 849 ' 1.500E-03 ' NSTO 751 SAXE 1 'x' LECT nxplTERM ECRO COMP 1*COUR 1 'press_1' ECRO COMP 1 GAUSS 1 ELEM LECT fe1TERMlist 101 102 103 104 105 106 107 108 109 110 axes 1 'PRESS[MPa]'** and so on****************************************list 841 842 843 844 845 846 847 848 849 axes 1 'PRESS[MPa]'*=============================================FINcon22.dgibi*CON22 Conical model with 12.8 kg TNT*quadratic elements*opti donn 'D:\Users\larchma\cast3m\pxpdroi1.procedur';opti donn 'D:\Users\larchma\cast3m\pxordpoi.procedur';OPTI echo 1;OPTI dime 3 elem cub8;***************************************************** Parameters *******************************************************************************************opening=0.1;length=10.0;dex=0.124;dexin= 0.0001;*Length of the pyramiddpy=0.001;****************************************************DENS 100;opex=dex*opening/length;dexfi= opex*2.;daiin= opex*2.;daifi= opening*2.;oppy=dpy*opening/length;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);78