2D Pole Deformation — Welcome to LS-DYNA Examples

Home » Knowledge Base » CESE / DUAL-CESE » Examples DUAL-CESE » 2D Pole Deformation

This is a 2D fluid-structure interaction (FSI) problem. In this test, a structure pole (shell element) stands in the middle of a flow field and its bottom end is fixed. When a strong wind flows from left to right, it pushes the pole and forces it to deform. It is simulated by the dual CESE FSI-mmm solver (this FSI solver is based on the moving mesh method (mmm), and can be more accurate than the immersed boundary method (ibm) in capturing small deformation problems).

Description

This is a 2D fluid-structure interaction (FSI) problem. In this test, a structure pole (shell element) stands in the middle of a flow field and its bottom end is fixed. When a strong wind flows from left to right, it pushes the pole and forces it to deform. It is simulated by the dual CESE FSI-mmm solver (this FSI solver is based on the moving mesh method (mmm), and can be more accurate than the immersed boundary method (ibm) in capturing small deformation problems).

Animated Result

pole_deform_pre1.gif

Keywords

*DUALCESE_MODEL

*DUALCESE_CONTROL_SOLVER

*DUALCESE_CONTROL_TIMESTEP

*DUALCESE_CONTROL_LIMITER

*DUALCESE_INCLUDE_MODEL

*DUALCESE_BOUNDARY_PRESCRIBED_SEGMENT_SET

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

*DUALCESE_INITIAL

*DUALCESE_PART

*DUALCESE_CONTROL_MESH_MOV

*DUALCESE_MAT_GAS

*DUALCESE_EOS_IDEAL_GAS

*DUALCESE_D3PLOT

*CONTROL_TERMINATION

*DATABASE_BINARY_D3PLOT

*DATABASE_BINARY_D3DUMP

*INCLUDE

*PART

*SECTION_SHELL

*MAT_ELASTIC

*CONTROL_CONTACT

*BOUNDARY_SPC_SET

*CONTROL_STRUCTURED

Reduced Input

*KEYWORD

*TITLE

shell_swing_in_flow

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

$ Units:

$ rho (density):  Kg/m^3

$ V   (velocity): m/s

$ P   (pressure): Pa

$ T   (temperature): K

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

*CONTROL_TERMINATION

$   endtim    endcyc     dtmin     endeng   endmas

0.015     20000

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

$ Handle output of state and restart data

*DATABASE_BINARY_D3PLOT

$  dt/cycl      lcdt      beam     npltc

1.0e-4                   0

*DATABASE_BINARY_D3DUMP

$  dt/cycl      lcdt      beam     npltc

10000

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

*INCLUDE

struct_setup.k

$

*DUALCESE_MODEL

fluid_setup.k

*END

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

 

*KEYWORD

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

*DUALCESE_CONTROL_SOLVER

$    ieqns     igeom iframeref

Euler        2D     fixed

$        NS        3D     fixed

*DUALCESE_CONTROL_TIMESTEP

$     iddt       cfl     dtint

2       0.9    0.1e-5

*DUALCESE_CONTROL_LIMITER

$    idlmt      alfa      beta      epsr

0       2.0       1.0      0.2

*DUALCESE_INCLUDE_MODEL

fluid_mesh1.k

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

$ 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

$     sf_u      sf_v      sf_w      sf_d      sf_p      sf_t

1.0e+2       0.0       0.0       1.5     2.0e5

$

$ open BCs

$

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

$     ssid

2

$         5

$

$ bottom solid BC_7

$

$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET

*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET

$     ssid

3

5

*DUALCESE_BOUNDARY_FSI_SEGMENT_SET

$     ssid

11

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

$ Setup the initial conditions for fluid

*DUALCESE_INITIAL

$      uic       vic       wic     rhoic       pic       tic       hic

0.0       0.0       0.0     1.225 1.01325e5

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

$ Setup fluid properties

*DUALCESE_PART

$      pid       mid     eosid  FSI-type MOVMSHALG

1         4         3       MMM      21

*DUALCESE_CONTROL_MESH_MOV

$       ID      IALG     NITER RELERRTOL  MXDISP_R

21         9                        1.5e-4

*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

vorticity

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

*END

Download

DUALCESE_example_2DPolemmm.zip