The proposed test case consists of a flexible gate, clamped at one end and free at the other which interacts with a mass of water initially confined in a free surface tank behind the gate. Due to the uncertainty of the material used, an elastic material as well as hyperelastic models (Ogden, Mooney-Rivlin) will be tested.
Fluid pressure fringes
*KEYWORD *TITLE *BOUNDARY_SPC_SET_BIRTH_DEATH *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *DATABASE_BINARY_D3PLOT *DATABASE_HISTORY_NODE *DATABASE_NODOUT *DEFINE_CURVE_TITLE *ELEMENT_SHELL *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_FSI *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_CONTROL_ADAPT_SIZE *ICFD_CONTROL_FSI *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_TIMESTEP *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *LOAD_BODY_Y *MAT_ELASTIC *MAT_MOONEY-RIVLIN_RUBBER *MAT_OGDEN_RUBBER *MESH_INTERF *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *NODE *PARAMETER *PART *SECTION_SHELL *SET_NODE_LIST *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R9 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Dam break with elastic gate opening
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 0.4
R dt_plot 5e-3
$
$--- Fluid
$
Rrho_fluid 1000
R mu_fluid 0.001
R dt_fluid 5e-4
R grav 9.81
R fsi_bt 1e-3
*INCLUDE
mesh.k
*INCLUDE
struc.k
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$# ttm dt
&T_end &dt_fluid
*ICFD_CONTROL_FSI
$# owc bt
0 &fsi_bt
*ICFD_CONTROL_OUTPUT
$# msglv
3
*ICFD_CONTROL_ADAPT_SIZE
$# asize
1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SECTION
$# sid
1
*ICFD_MAT
$# mid flg ro vis
1 1&rho_fluid &mu_fluid
*ICFD_MAT
$# mid flg ro vis
2 0
*ICFD_PART
$# pid secid mid
1 1 1
*ICFD_PART
$# pid secid mid
2 1 1
*ICFD_PART
$# pid secid mid
3 1 1
*ICFD_PART
$# pid secid mid
4 1 1
*ICFD_PART
$# pid secid mid
5 1 1
*ICFD_PART
$# pid secid mid
6 1 2
*ICFD_PART
$# pid secid mid
7 1 2
*ICFD_PART
$# pid secid mid
8 1 2
*ICFD_PART
$# pid secid mid
9 1 2
*ICFD_PART
$# pid secid mid
10 1 2
*ICFD_PART
$# pid secid mid
11 1 2
*ICFD_PART
$# pid secid mid
12 1 2
*ICFD_PART
$# pid secid mid
13 1 2
*ICFD_PART
$# pid secid mid
14 1 2
*ICFD_PART_VOL
$# pid secid mid
15 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 4 5 6
*ICFD_PART_VOL
$# pid secid mid
16 1 2
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
2 3 6 7 8 9 10 11
$# spid9 spid10 spid11 spid12
12 13 14
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FREESLIP
$# pid
1
*ICFD_BOUNDARY_NONSLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
4
*ICFD_BOUNDARY_FREESLIP
$# pid
5
*ICFD_BOUNDARY_FREESLIP
$# pid
7
*ICFD_BOUNDARY_FREESLIP
$# pid
8
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
9 2
*ICFD_BOUNDARY_FREESLIP
$# pid
10
*ICFD_BOUNDARY_NONSLIP
$# pid
11
*ICFD_BOUNDARY_FREESLIP
$# pid
12
*ICFD_BOUNDARY_FREESLIP
$# pid
13
*ICFD_BOUNDARY_FREESLIP
$# pid
14
*ICFD_BOUNDARY_FSI
$# pid
3
*ICFD_BOUNDARY_FSI
$# pid
4
*ICFD_BOUNDARY_FSI
$# pid
11
*LOAD_BODY_Y
$# lcid sf
1 1
*DEFINE_CURVE_TITLE
Gravity force
$# lcid sidr sfa sfo offa offo dattyp
1 &grav
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_CURVE_TITLE
Pressure outlet
$# lcid sidr sfa sfo offa offo dattyp
2
$# a1 o1
0.0 0.0
10000.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$# volid
30
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 3 4 5 7 8 9 10
$# pid9 pid10 pid11 pid12
11 12 13 14
*MESH_INTERF
$# pid
30
$# spid1
2
*MESH_INTERF
$# pid
30
$# spid1
6
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_TIMESTEP
$# outlv
1
*DATABASE_BINARY_D3PLOT
&dt_plot
*END
The proposed test case consists of a flexible gate, clamped at one end and free at the other which interacts with a mass of water initially confined in a free surface tank behind the gate. Due to the uncertainty of the material used, an elastic material as well as hyperelastic models (Ogden, Mooney-Rivlin) will be tested.
References :
[1] C. Antoci, M. Gallati, and S. Sibilla, “Numerical simulation of fluid-structure interaction by sph,” Computers and Structures, vol. 85, no. 11–14, pp. 879 – 890, 2007. Fourth MIT Conference on Computational Fluid and Solid Mechanics.
[2] ASME 27th International Conference on Offshore Mechanics and Arctic Engineering, Fluid structure interaction simulation using an incompressible SPH method, (Estoril,Portugal), OMAE2008, June 2008.
[3 Q. Yang, SPH Simulation of Fluid-Structure Interaction Problems with Application to Hovercraft. PhD thesis, Virginia Polytechnic Institute and State University, 2011.