Wave-Structure Interaction — Welcome to LS-DYNA Examples

Home » Knowledge Base » ALE » Intermediate Example » 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

Keywords

*KEYWORD
*TITLE
*CONTROL_ALE
*CONTROL_ENERGY
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_FSI
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DATABASE_BINARY_D3PLOT
*PART
*SECTION_SOLID_TITLE
*MAT_VACUUM_TITLE
*HOURGLASS_TITLE
*MAT_NULL_TITLE
*EOS_MURNAGHAN_TITLE
*MAT_RIGID_TITLE
*LOAD_BODY_Z
*DEFINE_CURVE_TITLE
*ALE_MULTI-MATERIAL_GROUP
*SET_MULTI-MATERIAL_GROUP_LIST
*CONSTRAINED_LAGRANGE_IN_SOLID_TITLE
*CONSTRAINED_GLOBAL
*ELEMENT_SOLID
*NODE
*END

Reduced Input

$# LS-DYNA Keyword file created by James M. Kennedy - KBS2 Inc. - March 2020
$#
$# units are kg, m, s, N, Pa, N-m
$# (kilogram, meter, second, Newton, Pascal, Newton-meter)
$#
*KEYWORD
*TITLE
$# title
ALE - Wave-Structure Interaction
$
*CONTROL_ALE
$# dct nadv meth afac bfac cfac dfac efac
-1 1 1 -1.0 0.0 0.0 0.0 0.0
$# start end aafac vfact vlimit ebc pref nsidebc
0.0 1.0000E20 1.0 1.0000E-6 0 0 0.0 0
$# ncpl nbkt imascl checkr beamin mmgpref pdfimx dtmufac
1 50 0 0.0 0.0 0 0.0 0.0
$# optimpp ialedr bndflx minmas
0 0 0 1.0000E-5
*CONTROL_ENERGY
$# hgen rwen slnten rylen irgen
2 2 2 2 2
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
4.0 0 0.0 0.0 1.00000E8 0
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.70 0 0.0 0.0 0 0 0
$# dt2msf dt2mslc imscl unused unused rmscl unused ihdo
0.0 0 0 0.0 0
$
*DATABASE_FSI
$# dt binary
0.001 1
$#dbfsi_id sid sidtype
101 4 1
*DATABASE_GLSTAT
$# dt binary lcur ioopt
0.001 0 0 1
*DATABASE_MATSUM
$# dt binary lcur ioopt
0.001 0 0 1
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.040 0 0 0 0
$# ioopt rate cutoff window type pset
0 0.0 0.0 0.0 0 0
$
*PART
Air
$# pid secid mid eosid hgid grav adpopt tmid
1 1 1 0 1 0 0 0
*SECTION_SOLID_TITLE
Air
$# secid elform aet
1 11 0
*MAT_VACUUM_TITLE
Air
$ mid ro
1
*HOURGLASS_TITLE
Air
$# hgid ihq qm ibq q1 q2 qb qw
1 1 1.00E-06
$
*PART
Water
$# pid secid mid eosid hgid grav adpopt tmid
2 2 2 2 2 0 0 0
*SECTION_SOLID_TITLE
Water
$# secid elform aet
2 11 0
*MAT_NULL_TITLE
Water
$ mid ro pc mu terod cerod ym pr
2 998.21 -50000.0 0.0 0.0 0.0 0.0 0.0
*EOS_MURNAGHAN_TITLE
Water
$# eosid gamma k0 v0
2 7.0 150000.0 0.0
*HOURGLASS_TITLE
Water
$# hgid ihq qm ibq q1 q2 qb qw
2 1 1.0000E-6
$
*PART
Column, rigid, fixed
$# pid secid mid eosid hgid grav adpopt tmid
4 4 4 0 0 0 0 0
*SECTION_SOLID_TITLE
Rigid column
$# secid elform aet
4 1 0
*MAT_RIGID_TITLE
Rigid column
$# mid ro e pr n couple m alias
4 100.0 2.0E4 0.25 0.0 0.0 0.0
$# cmo con1 con2
1.0 7 7
$#lco / a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
$
*LOAD_BODY_Z
$# lcid sf lciddr xc yc zc cid
1 9.81 0 0.0 0.0 0.0 0
*DEFINE_CURVE_TITLE
Gravity
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 0.0 0.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
0.050 1.0
100.0 1.0
$
*ALE_MULTI-MATERIAL_GROUP
$# sid idtype
1 1 Air
2 1 Water
*SET_MULTI-MATERIAL_GROUP_LIST
$# ammsid
101
$# ammgid1
1
*SET_MULTI-MATERIAL_GROUP_LIST
$# ammsid
102
$# ammgid1
2
*CONSTRAINED_LAGRANGE_IN_SOLID_TITLE
1 Column-to-Air
$# slave master sstyp mstyp nquad ctype direc mcoup
4 1 1 1 4 4 2 -102
$# start end pfac fric frcmin norm normtyp damp
0.0 1.0000E20 0.1 0.0 0.5 0 0 0.0
$# k hmin hmax ileak pleak lcidpor nvent iblock
0.0 0.0 0.0 0 0.1 0 0 0
$# boxid ipenchk intforc ialesof lagmul pfacmm thkf
0 0 0 0 0.0 0 0.0
*CONSTRAINED_LAGRANGE_IN_SOLID_TITLE
2 Column-to-Water
$# slave master sstyp mstyp nquad ctype direc mcoup
4 2 1 1 4 4 2 -102
$# start end pfac fric frcmin norm normtyp damp
0.0 1.0000E20 0.1 0.0 0.5 0 0 0.0
$# k hmin hmax ileak pleak lcidpor nvent iblock
0.0 0.0 0.0 0 0.1 0 0 0
$# boxid ipenchk intforc ialesof lagmul pfacmm thkf
0 0 0 0 0.0 0 0.0
$
*CONSTRAINED_GLOBAL
$# ALE - global boundary constraint planes
$# tc rc dir x y z tol
1 0 1 0.0 0.0 0.0 0.0
1 0 1 1.60 0.0 0.0 0.0
2 0 2 0.0 0.0 0.0 0.0
2 0 2 0.0 0.61 0.0 0.0
3 0 3 0.0 0.0 0.0 0.0
$
*END

Figures

ale_vs_sale_0.010_cells_fx.jpgsale_0.010_vs_0.005_cells_fx.jpg

Animated Result

wavestructure_ani2.gif

Download