x
Our website uses cookies. By using the website you agree ot its use. More information can be found in our privacy policy.

Reduced Input

*KEYWORD 

*TITLE 

fluid_2d_bullet_interaction 

*CONTROL_TERMINATION 

$   endtim    endcyc     dtmin     endeng   endmas 

      .004    200000 

$=========================================================================== 

$ Units: 

$ rho (density):  Kg/m^3 

$ V   (velocity): m/s 

$ P   (pressure): Pa 

$ T   (temperature): K 

$============================================================================ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ Handle output of state and restart data 

*DATABASE_BINARY_D3PLOT 

$  dt/cycl      lcdt      beam     npltc 

    2.e-5                   0 

*DATABASE_BINARY_D3DUMP 

$  dt/cycl      lcdt      beam     npltc 

     10000 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

*INCLUDE 

struct_setup.k 

$ 

*DUALCESE_MODEL 

            dcese_setup.k 

*END 

$--------------------------------------------------------------------------------------------------------------------------------- 

*KEYWORD 

*INCLUDE 

struct_mesh.k 

$ Time step control for structure side 

*CONTROL_TIMESTEP 

$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st 

 1.0000E-7     0.900         0     0.000     0.000                   0         0 

$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl      

     0.000         0         0                         0.000 

*DEFINE_CURVE_TITLE 

minimum_time 

$#    lcid      sidr       sfa       sfo      offa      offo    dattyp       

         2         0     1.000     1.000     0.000     0.000         0 

$#                a1                  o1   

               0.000        1.000000e-07 

          100.000000        1.000000e-07 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

*PART 

$#                                                                         title 

Bullet_part 

$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid 

         2         2         1         0         0         0         0         0 

*SECTION_SHELL_TITLE 

bullet_motion 

$#   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp 

         2        12  1.000000         2         1         0         0         1 

$#      t1        t2        t3        t4      nloc     marea      idof    edgset 

  1.000000  1.000000  1.000000  1.000000     0.000     0.000     0.000         0 

*MAT_RIGID_TITLE 

bullet_material 

$#     mid        ro         e        pr         n    couple         m     alias 

         1    7850.0  2.10E+11     0.200     0.000     0.000     0.000           

$#     cmo      con1      con2     

     1.000        5.        7. 

$#lco or a1        a2        a3        v1        v2        v3   

     0.000     0.000     0.000     0.000     0.000     0.000 

$#                                                                         title 

*BOUNDARY_PRESCRIBED_MOTION_RIGID_ID 

$#      id                                                               heading 

         1bullet_part 

$#     pid       dof       vad      lcid        sf       vid     death     birth 

         2         1         0         3  1.000000         01.0000E+28     0.000 

*DEFINE_CURVE 

$#    lcid      sidr       sfa       sfo      offa      offo    dattyp       

         3         0   1.00E-6       1.0     0.000     0.000         0 

$#                a1                  o1   

                 0.0                 0.0 

                 1.0               340.0 

                 2.0               500.0 

                 5.0              1000.0 

                 7.5              1250.0 

                10.0              1500.0 

             1.0e+05              1500.0 

*end 

$--------------------------------------------------------------------------------------------------------------------------- 

*KEYWORD 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

*DUALCESE_CONTROL_SOLVER 

$    ieqns     igeom iframeref 

     Euler        2D     fixed 

*DUALCESE_CONTROL_TIMESTEP 

$     iddt       cfl     dtint 

         2       0.5    0.1e-6                

*DUALCESE_CONTROL_LIMITER 

$    idlmt      alfa      beta      epsr 

         0       2.0       1.0      0.05 

*DUALCESE_INCLUDE_MODEL 

dcese_mesh.k 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ Setup the boundary conditions for fluid 

$ open BCs 

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET 

$     ssid 

         2 

$ 

$ bottom solid BCs 

$ 

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET 

$     ssid 

         1 

         3 

         4 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ Setup the initial conditions for fluid 

$ 

*DUALCESE_INITIAL 

$      uic       vic       wic     rhoic       pic       tic       hic 

       0.0       0.0       0.0      1.12   101325.     300.0 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ Setup fluid properties 

$ 

*DUALCESE_PART 

$      pid       mid     eosid  FSI-type MOVMSHALG  

         1         7         5       IBM 

*DUALCESE_MAT_001 

$      mid        c1        c2      prnd 

         7  1.458e-6     110.4      0.72 

*DUALCESE_EOS_IDEAL_GAS 

$    eosid        cv        cp        t0        v0 

         5     717.5    1004.5 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ Handle output of state and restart data 

*DUALCESE_D3PLOT 

density 

pressure 

temperature 

velocity 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

*END