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_solid_shell_interaction 

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

$ Units: 

$ rho (density):  Kg/m^3 

$ V   (velocity): m/s 

$ P   (pressure): Pa 

$ T   (temperature): K 

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

*CONTROL_TERMINATION 

$   endtim    endcyc     dtmin     endeng   endmas 

      0.03     20000 

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

$ Handle output of state and restart data 

*DATABASE_BINARY_D3PLOT 

$  dt/cycl      lcdt      beam     npltc 

    5.0e-4                   0 

$    1.0e-7                   0 

*DATABASE_BINARY_D3DUMP 

$  dt/cycl      lcdt      beam     npltc 

      1000 

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

*INCLUDE 

struct_setup.k 

*DUALCESE_MODEL 

           dcese_setup.k 

*END 

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

 

*KEYWORD 

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

*DUALCESE_CONTROL_SOLVER 

$    ieqns     igeom iframeref 

     Euler        3D     fixed 

$        NS        3D     fixed 

*DUALCESE_CONTROL_TIMESTEP 

$     iddt       cfl     dtint 

         2      .99    0.1e-5                

*DUALCESE_CONTROL_LIMITER 

$    idlmt      alfa      beta      epsr 

         0       1.0       1.0       .1 

*DUALCESE_INCLUDE_MODEL 

dcese_m1h.k 

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

*DUALCESE_FSI_EXCLUDE 

    10005 

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

$ Setup the boundary conditions for fluid 

$ Prescribed inlet BC (left) 

*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET 

$     ssid 

         1 

$   lcid_u    lcid_v    lcid_w    lcid_d    lcid_p    lcid_t 

                                                          -1 

  1. 0.0       0.0      1.22 1.51988e5 

$ open BCs 

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET 

$     ssid 

         2 

         3 

         5 

        10 

$ bottom solid BC_7 

$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET 

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET 

$     ssid 

         7 

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

$ Setup the initial conditions for fluid 

*DUALCESE_INITIAL 

$      uic       vic       wic     rhoic       pic       tic       hic 

       0.0       0.0       0.0      1.22 1.01325e5 

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

$ Setup fluid properties 

*DUALCESE_PART 

$      pid       mid     eosid FSI-type MOVMSHALG  

         1         4         3      IBM 

*DUALCESE_MAT_GAS 

$      mid        c1        c2      prnd 

         4  1.458e-6     110.4      0.72 

*DUALCESE_EOS_IDEAL_GAS 

$    eosid        cv        cp        t0        v0 

         3     713.5    1001.5 

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

$ Handle output of state and restart data 

*DUALCESE_D3PLOT 

density 

pressure 

velocity 

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

*END 

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

 

*KEYWORD 

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

*PART 

  solid 

$      PID     SECID       MID     EOSID      HGID      GRAV    ADPOPT      TMID 

         3         6         5         0         1 

*SECTION_SOLID 

         6         1 

*MAT_ELASTIC 

 5,1.4e+3,1.15e+10,0.1,0.0,0.0,0.0 

*HOURGLASS 

         1         4     

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

*PART  

   flimsy 

 15,17,19,0,0,0,0,0 

*SECTION_SHELL 

 17,2,0.0,0.0,0.0 

 0.003,0.003,0.003,0.003 

*MAT_ELASTIC 

 19,6.4e+2,5.15e+8,0.1,0.0,0.0,0.0 

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

*PART  

  bottom_table 

 10005,7,9,0,0,0,0,0 

*SECTION_SHELL 

 7,16,0.0,0.0,0.0 

 0.002,0.002,0.002,0.002 

*MAT_ELASTIC 

 9,5.0e+3,1.0e+11,0.1,0.0,0.0,0.0 

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

*CONTACT_AUTOMATIC_GENERAL_INTERIOR_MPP 

$*CONTACT_AUTOMATIC_GENERAL_INTERIOR 

$*CONTACT_AUTOMATIC_SINGLE_SURFACE 

,10 

 910,0,2 

 0.0,0.0,0.0,0.0,50.0,2,0.0,0.0 

 

 1 

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

*CONTROL_CONTACT 

   0.10000                    

 

 

         2 

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

*SET_PART_LIST 

 910,0.0,0.0,0.0,0.0 

 3,15,10005 

*SET_PART_LIST 

 922,0.0,0.0,0.0,0.0 

 3,15 

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

*LOAD_BODY_Z 

$ lcid     sf    lciddr  xc  yc  zc 

  101,9.81e-0 

*LOAD_BODY_PARTS 

 922 

*BOUNDARY_SPC_SET 

79,0,1,1,1,1,1,1 

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

$ Define the load curves 

*DEFINE_CURVE 

$     lcid      sidr       sfa       sfo      offa      offo    dattyp 

       101 

$                 a1                  o1 

                 0.0                 1.0 

           1000000.0                 1.0 

*DEFINE_CURVE 

$     lcid      sidr       sfa       sfo      offa      offo    dattyp 

       11 

$                 a1                  o1 

                 0.0           0.98325e5 

           1000000.0           0.98325e5 

*DEFINE_CURVE 

$     lcid      sidr       sfa       sfo      offa      offo    dattyp 

       21 

$                 a1                  o1 

                 0.0           1.51988e5 

           1000000.0           1.51988e5 

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

*INCLUDE 

struct_m1h.k 

*END 

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