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

$ lcid_u lcid_v lcid_w lcid_d lcid_p lcid_t

-1

0.0 0.0 1.22 1.51988e5

$ open BCs

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

$ ssid

$ bottom solid BC_7

$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

$ ssid

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

$ 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

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

*CONTROL_CONTACT

0.10000

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

*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

$ a1 o1

0.0 0.98325e5

1000000.0 0.98325e5

*DEFINE_CURVE

$ lcid sidr sfa sfo offa offo dattyp

$ a1 o1

0.0 1.51988e5

1000000.0 1.51988e5

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

*INCLUDE

struct_m1h.k

*END

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