Traffic Sign Deflection — Welcome to LS-DYNA Examples

Home » Knowledge Base » CESE / DUAL-CESE » Examples DUAL-CESE » Traffic Sign Deflection

In this 3D FSI test case, the reaction of a traffic light plate due to a strong wind gust is simulated.

Description

In this 3D FSI test case, the reaction of a traffic light plate due to a strong wind gust is simulated.

Animated Result

traffic_light_plate_vibrating_v2.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_SOLID_WALL_SEGMENT_SET

*DUALCESE_INITIAL

*DUALCESE_PART

*DUALCESE_MAT_GAS

*DUALCESE_EOS_IDEAL_GAS

*DUALCESE_D3PLOT

*CONTROL_TERMINATION

*DATABASE_BINARY_D3PLOT

*DATABASE_BINARY_D3DUMP

*INCLUDE

*PART

*SECTION_SOLID

*SECTION_SHELL

*MAT_ELASTIC

*CONTACT_TIED_NODES_TO_SURFACE_OFFSET

*HOURGLASS

*CONTACT_AUTOMATIC_GENERAL_INTERIOR_MPP

*CONTROL_CONTACT

*SET_PART_LIST

*LOAD_BODY_Z

*LOAD_BODY_PARTS

*BOUNDARY_SPC_SET

*DEFINE_CURVE

Reduced Input

*KEYWORD

*TITLE

traffic_sign_vibration

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

$ Units:

$ rho (density):  Kg/m^3

$ V   (velocity): m/s

$ P   (pressure): Pa

$ T   (temperature): K

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

*CONTROL_TERMINATION

$   endtim    endcyc     dtmin     endeng   endmas

0.5     20000

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

$ 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

131

$                 a1                  o1

0.000               0.000

0.010               100.0

0.030               200.0

0.080               200.0

0.100                 0.0

0.160                 0.0

0.180               150.0

0.205                 0.0

0.265                 0.0

0.280               150.0

0.330               150.0

0.350                 0.0

0.420                 0.0

0.450               150.0

1000000.0               150.0

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

$ Handle output of state and restart data

*DATABASE_BINARY_D3PLOT

$  dt/cycl      lcdt      beam     npltc

5.0e-3                   0

$    5.0e-4                   0

*DATABASE_BINARY_D3DUMP

$  dt/cycl      lcdt      beam     npltc

5000

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

*INCLUDE

struct_setup.k

*DUALCESE_MODEL

dcese_setup.k

*END

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

 

*KEYWORD

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

*DUALCESE_CONTROL_SOLVER

$    ieqns     igeom iframeref

Euler        3D     fixed

*DUALCESE_CONTROL_TIMESTEP

$     iddt       cfl     dtint

2       0.9    0.0001

*DUALCESE_CONTROL_LIMITER

$    idlmt      alfa      beta      epsr

0       2.0       1.0       .01

*DUALCESE_INCLUDE_MODEL

fluid_mesh3.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        -1        -1        -1                  -1

131        -1        -1        -1        -1        -1

$     sf_u      sf_v      sf_w      sf_d      sf_p      sf_t

2.51988e+5

$ open BCs

*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET

$     ssid

2

3

5

10

*DUALCESE_BOUNDARY_SOLID_WALL_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

5         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

sign

111,17,19,0,0,0,0,0

*SECTION_SHELL

17,16,0.0,0.0,0.0

$ 0.004,0.004,0.004,0.004

0.03,0.03,0.03,0.03

*MAT_ELASTIC

19,5.0e+3,6.15e+11,0.1,0.0,0.0,0.0

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

*CONTACT_TIED_NODES_TO_SURFACE_OFFSET

48,89,4,0

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

 

1

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

*CONTACT_TIED_NODES_TO_SURFACE_OFFSET

58,89,4,0

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

 

1

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

*CONTACT_TIED_NODES_TO_SURFACE_OFFSET

49,111,4,3

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

 

1

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

*CONTACT_TIED_NODES_TO_SURFACE_OFFSET

59,111,4,3

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

 

1

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

*PART

pillar

$      PID     SECID       MID     EOSID      HGID      GRAV    ADPOPT      TMID

15         6         5         0         1

*SECTION_SOLID

6         1

*MAT_ELASTIC

5,6.4e+3,7.15e+11,0.1,0.0,0.0,0.0

*HOURGLASS

1         4

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

*PART

stand_1

$      PID     SECID       MID     EOSID      HGID      GRAV    ADPOPT      TMID

25         7         5         0         1

*SECTION_SOLID

7         1

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

*PART

stand_1

$      PID     SECID       MID     EOSID      HGID      GRAV    ADPOPT      TMID

35         8         5         0         1

*SECTION_SOLID

8         1

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

*CONTACT_AUTOMATIC_GENERAL_INTERIOR_MPP

,10

$*CONTACT_AUTOMATIC_SINGLE_SURFACE

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

15,25,35,111

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

*LOAD_BODY_Z

$ lcid     sf    lciddr  xc  yc  zc

101,9.81

*LOAD_BODY_PARTS

910

*BOUNDARY_SPC_SET

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

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

*INCLUDE

struct_mesh.k

*END

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

Download

DUALCESE_example_Trafficsign.zip