Traffic Sign Deflection
In this 3D FSI test case, the reaction of a traffic light plate due to a strong wind gust is simulated.
*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
*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
$============================================================================