Wave-Structure Interaction
A container is partially filled with water subjected to gravity. As the wave propagates, it hits a rigid column placed roughly in the middle of the container. The total contact forces are recorded for the entire simulation, and compared to experimental data: Gómez-Gesteira M, Crespo AJC, Rogers BD, Dalrymple RA, Domínguez JM and Barreiro A. 2012. SPHysics - development of a free-surface fluid solver- Part 2: Efficiency and test cases. Computers & Geosciences, 48: 300-307. doi:10.1016/j.cageo.2012.02.028.
https://www.dynaexamples.com/sph/intermediate-examples/wavestructure
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Wave-Structure Interaction
A container is partially filled with water subjected to gravity. As the wave propagates, it hits a rigid column placed roughly in the middle of the container. The total contact forces are recorded for the entire simulation, and compared to experimental data: Gómez-Gesteira M, Crespo AJC, Rogers BD, Dalrymple RA, Domínguez JM and Barreiro A. 2012. SPHysics - development of a free-surface fluid solver- Part 2: Efficiency and test cases. Computers & Geosciences, 48: 300-307. doi:10.1016/j.cageo.2012.02.028.
Watch animated result on YouTube
$# LS-DYNA Keyword file created by LS-PrePost(R) V4.5.14 - 14Feb2018
$# Created on Apr-6-2018 (17:51:30)
$-----------------------------------------------------------------------------
$
$ Example provided by Edouard (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2018 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 R10 (or higher)
$X
$X------------------------------------------------------------------------------
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
$# title
Wave-Structure Interaction, 3D
*CONTROL_BULK_VISCOSITY
$# q1 q2 type btype
0.0011.0000E-12 1 0
*CONTROL_SPH
$# ncbs boxid dt idim memory form start maxv
1 11.00000E20 3 500 15 0.0 20.0
$# cont deriv ini ishow ierod icont iavis isymp
0 0 0 1 0 0 1 100
$# ithk
0
*DEFINE_BOX
$# boxid xmn xmx ymn ymx zmn zmx
1 -0.05 1.65 -0.31 0.31 -0.05 2.0
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
4.0 0 0.0 0.0 0.0 0
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.45 0 0.0 0.0 0 0 0
$# dt2msf dt2mslc imscl unused unused rmscl
0.0 0 0 0.0
*DATABASE_RCFORC
$# dt binary lcur ioopt
0.001 0 0 1
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.005 0 0 0 0
$# ioopt
0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ GRAVITY LOADING $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*LOAD_BODY_Z
$# lcid sf lciddr xc yc zc cid
1 1.0 0 0.0 0.0 0.0 0
*DEFINE_CURVE_TITLE
Gravity
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
0.05 9.81
1000.0 9.81
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ FEM STRUCTURE $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
Column
$# pid secid mid eosid hgid grav adpopt tmid
1 2 3 0 0 0 0 0
*SECTION_SOLID
$# secid elform aet
2 1 0
*MAT_RIGID
$# mid ro e pr n couple m alias
3 0.001 2.0 0.25 0.0 0.0 0.0
$# cmo con1 con2
1.0 7 7
$#lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*PART
$# title
Container
$# pid secid mid eosid hgid grav adpopt tmid
3 2 3 0 0 0 0 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ SPH FLUID (WATER) $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
SPH Water
$# pid secid mid eosid hgid grav adpopt tmid
2 1 1 1 0 0 0 0
*SECTION_SPH
$# secid cslh hmin hmax sphini death start
1 1.2 1.0 1.0 0.01.00000E20 0.0
*MAT_NULL
$# mid ro pc mu terod cerod ym pr
1 1000.0 -50000.0 0.001 0.0 0.0 0.0 0.0
*EOS_MURNAGHAN
$# eosid gamma k0 v0
1 7.0 150000.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ FLUID-STRUCTURE INTERACTION $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_MPP_ID
$# cid title
1SPH to Column
$# ignore bucket lcbucket ns2track inititer parmax unused cparm8
0 50 0 3 2 1.0005 0
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
2 1 3 3 0 0 0 0
$# fs fd dc vc vdc penchk bt dt
0.0 0.0 0.0 0.0 0.0 0 0.01.00000E20
$# sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.5 0.0 0.0 1.0 1.0 1.0 1.0
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
1 0.1 0 1.025 2.0 2 50 1
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_MPP_ID
$# cid title
2SPH to Container
$# ignore bucket lcbucket ns2track inititer parmax unused cparm8
0 50 0 3 2 1.0005 0
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
2 3 3 3 0 0 0 0
$# fs fd dc vc vdc penchk bt dt
0.0 0.0 0.0 0.0 0.0 0 0.01.00000E20
$# sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.5 0.0 0.0 1.0 1.0 1.0 1.0
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
1 0.1 0 1.025 2.0 2 50 1
*RIGIDWALL_PLANAR_ID
$# id title
1
$# nsid nsidex boxid offset birth death rwksf
0 0 0 0.0 0.01.00000E20 0.01
$# xt yt zt xh yh zh fric wvel
0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0
*END
*CONTROL_BULK_VISCOSITY *CONTROL_SPH *DEFINE_BOX *DATABASE_RCFORC *LOAD_BODY_Z *DEFINE_CURVE_TITLE *SECTION_SPH *MAT_NULL *EOS_MURNAGHAN *CONTACT_AUTOMATIC_NODES_TO_SURFACE_MPP_ID *RIGIDWALL_PLANAR_ID
- sph_wave_structure.zip
- sph_wave_structure.tar.gz