https://ftp.lstc.com/anonymous/outgoing/hao/sale/models/waterlanding/
这是网络上的一个速度边界的SALE造波方法,包含k文。
速度边界造波(应该也叫无板造波吧)
BOUNDARY_PRESCRIBED_MOTION_SET设置输入端和输出端x方向速度,以及输入端y方向速度,去模拟波浪。
但直接下载的k文件只能使用高版本lsrun进行仿真还需要开启双精度,还有就是他进行了流固耦合,可能是计算量比较大,我直接使用仿真的话结果是很不理想的
根据以上例子删除相关结构,只保留水后进行了仿真还算是可以,边界流速和波高的关系还有待商榷,一起探讨。
main.k
$# LS-DYNA Keyword file created by LS-PrePost(R) V4.9.16-04Aug2023
$# Created on Oct-17-2023 (10:34:01)
$$ HM_OUTPUT_DECK created 10:36:08 07-07-2017 by HyperMesh Version 13.0.0.119
$$ Ls-dyna Input Deck Generated by HyperMesh Version : 13.0.0.119
$$ Generated using HyperMesh-Ls-dyna 971_R6.1 Template Version : 13.0
*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
*CONTROL_ALE
$# dct nadv meth afac bfac cfac dfac efac
0 0 1 -1.0 0.0 0.0 0.0 0.0
$# start end aafac vfact prit ebc pref nsidebc
0.01.00000E20 1.01.00000E-6 0 0 101300.0 0
$# ncpl nbkt imascl checkr beamin mmgpref pdifmx dtmufac
1 50 0 0.0
$# optimpp ialedr bndflx minmas
0 0 01.00000E-5
$-----------------------------------ALE-CARDS-----------------------------------
*CONTROL_MPP_DECOMPOSITION_DISTRIBUTE_ALE_ELEMENTS
*CONTROL_OUTPUT
$# npopt neecho nrefup iaccop opifs ipnint ikedit iflush
0 0 0 0 0.0 0 100 100
$# iprtf ierode tet10s8 msgmax ipcurv gmdt ip1dblt eocs
0 0 2 50 0 0.0 0 0
$# tolev newleg frfreq minfo solsig msgflg cdetol
2 0 1 0 0 0 10.0
$# phschng demden icrfile spc2bnd penout shlsig hisnout engout
0 0
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
0.4 0 0.0 0.0 0.0
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.7 0 0.0 0.0 0 0 0
$# dt2msf dt2mslc imscl unused unused rmscl emscl ihdo
0.0 0 0 0.0
$--------------------------------BOUNDARY-CONDITIONS----------------------------
*BOUNDARY_NON_REFLECTING
$# ssid ad as
1 0.0 0.0
*BOUNDARY_PRESCRIBED_MOTION_SET
$----wave making (3) and flow-in(2)/flow-out(12)
$# nsid dof vad lcid sf vid death birth
2 1 0 3 1.0 01.00000E28 0.0
$# nsid dof vad lcid sf vid death birth
12 1 0 3 1.0 01.00000E28 0.0
$# nsid dof vad lcid sf vid death birth
3 2 0 2 1.0 01.00000E28 0.0
*BOUNDARY_SPC_SET_ID
$# id heading
0cemian
$----no side flow-in/flow-out
$# nsid cid dofx dofy dofz dofrx dofry dofrz
11 0 0 0 1 1 1 1
$# id heading
0down
$----bottom flow along x direction only
$# nsid cid dofx dofy dofz dofrx dofry dofrz
4 0 0 1 1 1 1 1
*NODE
$# nid x y z tc rc
886400 -6.0 -3.0 0.0 0 0
2000000 -6.0 0.0 4.0 0 0
2000001 6.0 0.0 0.0 0 0
*ALE_AMBIENT_HYDROSTATIC
$# alesid stype vecid grav pbase ramptlc
2 2 1 9.8065 101300.0 0
$# nid mmgbl
886400 1
2000000 2
*ALE_MULTI-MATERIAL_GROUP
$# sid idtype
1 1
2 1
$-------------------------------------GRAVITY/LOADING---------------------------
*LOAD_BODY_Y
$# lcid sf lciddr xc yc zc cid
5 -9.80665 0 0.0 0.0 0.0 0
$-------------------------------------------------------------------------------
*PART
$# title
air
$# pid secid mid eosid hgid grav adpopt tmid
1 1 3 1 0 0 0 0
*SECTION_SOLID
$# secid elform aet unused unused unused cohoff gaskeit
1 11
*MAT_NULL_TITLE
Air
$# mid ro pc mu terod cerod ym pr
3 1.29 -10.0 0.0 0.0 0.0 0.0 0.0
*EOS_LINEAR_POLYNOMIAL_TITLE
Air
$# eosid c0 c1 c2 c3 c4 c5 c6
1 0.0 0.0 0.0 0.0 0.4 0.4 0.0
$# e0 v0
250000.0 0.0
*PART
$# title
water
$# pid secid mid eosid hgid grav adpopt tmid
2 1 2 2 0 0 0 0
*MAT_NULL_TITLE
Water
$# mid ro pc mu terod cerod ym pr
2 998.21 -10.0 0.0 0.0 0.0 0.0 0.0
*EOS_GRUNEISEN_TITLE
Water
$# eosid c s1 s2 s3 gamao a e0
2 1647.0 1.921 -0.096 0.0 0.35 0.0 289500.0
$# v0 unused lcid
0.0
$----------------------------------ASCII-DATABASE-------------------------------
*DATABASE_GLSTAT
$# dt binary lcur ioopt
0.004 1 0 1
*DATABASE_MATSUM
$# dt binary lcur ioopt
0.004 1 0 1
*DATABASE_NCFORC
$# dt binary lcur ioopt
0.004 0 0 1
*DATABASE_NODOUT
$# dt binary lcur ioopt option1 option2
0.004 1 0 1 0.0 0
*DATABASE_RBDOUT
$# dt binary lcur ioopt
0.004 1 0 1
*DATABASE_BINARY_D3DUMP
$# cycl lcdt beam npltc psetid
0.0 0 0 0 0
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.005 0 0 0 0
$# ioopt rate cutoff window type pset
0
*DATABASE_EXTENT_BINARY
$# neiph neips maxint strflg sigflg epsflg rltflg engflg
0 0 3 0 1 1 1 1
$# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat
0 0 0 1 1 3 1 0
$# nintsld pkp_sen sclp hydro msscl therm intout nodout
1 0 1.0 0 0 0STRESS STRESS
$# dtdt resplt neipb quadsld cubsld deleres
1 0 0
*DEFINE_CURVE_TITLE
sin
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 0.04 -0.15 0.0 0.0 0 0
$# a1 o1
0.0 0.0
0.1 6.18034
0.2 11.7557
0.3 16.18034
0.4 19.02113
0.5 20.0
0.6 19.02113
0.7 16.18034
0.8 11.7557
0.9 6.18034
1.0 2.4503000862e-15
1.1 -6.18034
1.2 -11.7557
1.3 -16.18034
1.4 -19.02113
1.5 -20.0
1.6 -19.02113
1.7 -16.18034
1.8 -11.7557
1.9 -6.18034
2.0 -4.9005900080e-15
2.1 6.18034
2.2 11.7557
2.3 16.18034
2.4 19.02113
2.5 20.0
2.6 19.02113
2.7 16.18034
2.8 11.7557
2.9 6.18034
3.0 7.3508898824e-15
3.1 -6.18034
3.2 -11.7557
3.3 -16.18034
3.4 -19.02113
3.5 -20.0
3.6 -19.02113
3.7 -16.18034
3.8 -11.7557
3.9 -6.18034
4.0 -9.8011901804e-15
4.1 6.18034
4.2 11.7557
4.3 16.18034
4.4 19.02113
4.5 20.0
4.6 19.02113
4.7 16.18034
4.8 11.7557
4.9 6.18034
5.0 1.2251499796e-14
5.1 -6.18034
5.2 -11.7557
5.3 -16.18034
5.4 -19.02113
5.5 -20.0
5.6 -19.02113
5.7 -16.18034
5.8 -11.7557
5.9 -6.18034
6.0 -1.4701800094e-14
6.1 6.18034
6.2 11.7557
6.3 16.18034
6.4 19.02113
6.5 20.0
6.6 19.02113
6.7 16.18034
6.8 11.7557
6.9 6.18034
7.0 1.7152100392e-14
7.1 -6.18034
7.2 -11.7557
7.3 -16.18034
7.4 -19.02113
7.5 -20.0
7.6 -19.02113
7.7 -16.18034
7.8 -11.7557
7.9 -6.18034
8.0 -1.9602400690e-14
8.1 6.18034
8.2 11.7557
8.3 16.18034
8.4 19.02113
8.5 20.0
8.6 19.02113
8.7 16.18034
8.8 11.7557
8.9 6.18034
9.0 2.2052699293e-14
9.1 -6.18034
9.2 -11.7557
9.3 -16.18034
9.4 -19.02113
9.5 -20.0
9.6 -19.02113
9.7 -16.18034
9.8 -11.7557
9.9 -6.18034
10.0 -2.4502999591e-14
*DEFINE_CURVE_TITLE
bo su
$# lcid sidr sfa sfo offa offo dattyp lcint
3 0 1.0 0.3 0.0 0.0 0 0
$# a1 o1
0.0 20.0
1.0000000000e+07 20.0
*DEFINE_CURVE_TITLE
load body y
$# lcid sidr sfa sfo offa offo dattyp lcint
5 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 1.0
200000.0 1.0
*DEFINE_CURVE_TITLE
x vel
$# lcid sidr sfa sfo offa offo dattyp lcint
6 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
0.2 0.0
0.2001 -9.144
0.8 -9.144
*DEFINE_CURVE_TITLE
y vel
$# lcid sidr sfa sfo offa offo dattyp lcint
7 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
0.2 0.0
0.2001 0.3
0.8 0.3
*DEFINE_VECTOR
$# vid xt yt zt xh yh zh cid
1 0.0 -1.0 0.0 0.0 0.0 0.0 0
*INCLUDE
mesh_sale.k
*INITIAL_VOLUME_FRACTION_GEOMETRY
$# fmsid fmidtyp bammg ntrace
6 1 1 0
$# conttyp fillopt fammg vx vy vz unused unused
3 0 2 0.0 0.0 0.0
$# x0 y0 z0 xcos ycos zcos unused unused
1.0 0.0 0.0 0.0 1.0 0.0
*INITIAL_HYDROSTATIC_ALE
$# alesid stype vecid grav pbase -
1 2 1 9.8065 101300.0
$# nid mmgblo
886400 1
2000000 2
*INITIAL_VELOCITY
$# nsid nsidex boxid irigid icid
1 0 0 0 0
$# vx vy vz vxr vyr vzr
6.0 0.0 0.0 0.0 0.0 0.0
$---------------------------------INITIAL-CONDITIONS----------------------------
$*INITIAL_VELOCITY
$ 1 0 0 0 0
$ 6.0 0.0 0.0 0.0 0.0 0.0
*MAT_RIGID_TITLE
plane
$# mid ro e pr n couple m alias
8 7800.02.10000E11 0.28 0.0 0.0 0.0
$# cmo con1 con2
0.0 0 0
$#lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*SECTION_SHELL_TITLE
plane
$# secid elform shrf nip propt qr/irid icomp setyp
8 2 1.0 2 1.0 0 0 1
$# t1 t2 t3 t4 nloc marea idof edgset
0.01 0.01 0.01 0.01 0.0 0.0 0.0 0
*SET_NODE_LIST
$# sid da1 da2 da3 da4 solver its -
3 0.0 0.0 0.0 0.0MECH 1
$# nid1 nid2 nid3 nid4 nid5 nid6 nid7 nid8
2000000 0 0 0 0 0 0 0
*SET_MULTI-MATERIAL_GROUP_LIST_TITLE
water
$# ammsid
1
$# ammgid1 ammgid2 ammgid3 ammgid4 ammgid5 ammgid6 ammgid7 ammgid8
2 0 0 0 0 0 0 0
*END
mesh_sale.k
$# LS-DYNA Keyword file created by LS-PrePost(R) V4.9.16-04Aug2023
$# Created on Oct-17-2023 (10:34:01)
*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
*ALE_STRUCTURED_MESH
$ mshid pid nbid ebid
$# mshid dpid nbid ebid unused unused unused tdeath
1 6 5000001 3000001
$ nptx npty nptz nid0 lcsid
$# cpidx cpidy cpidz nid0 lcsid
1001 1002 1003
*ALE_STRUCTURED_MESH_CONTROL_POINTS
$# cpid unused icase sfo unused offo
1001 0 1.0 0.0
$# n x ratio
1 -6.0 0.0
121 6.0 0.0
*ALE_STRUCTURED_MESH_CONTROL_POINTS
$# cpid unused icase sfo unused offo
1002 0 1.0 0.0
$# n x ratio
1 -3.0 0.0
61 3.0 0.0
*ALE_STRUCTURED_MESH_CONTROL_POINTS
$# cpid unused icase sfo unused offo
1003 0 0.0 0.0
$# n x ratio
1 -4.0 0.0
81 4.0 0.0
*SET_NODE_GENERAL
$# sid da1 da2 da3 da4 solver its -
1 0.0 0.0 0.0 0.0MECH 1
$# option mshid imin imax jmin jmax kmin kmax
SALECPT 1 1 120 1 61 1 81
*SET_NODE_GENERAL
$# sid da1 da2 da3 da4 solver its -
2 0.0 0.0 0.0 0.0MECH 1
$# option mshid imin imax jmin jmax kmin kmax
SALECPT 1 1 2 1 61 1 81
*SET_NODE_GENERAL
$# sid da1 da2 da3 da4 solver its -
4 0.0 0.0 0.0 0.0MECH 1
$# option mshid -x +x -y +y -z +z
SALEFAC 1 0 0 1 0 0 0
*SET_NODE_GENERAL
$ SID
$# sid da1 da2 da3 da4 solver its -
11 MECH 1
$ OPTION MSHID XMN XMX YMN YMX ZMN ZMX
$# option mshid -x +x -y +y -z +z
SALEFAC 1 1 1
*SET_NODE_GENERAL
$# sid da1 da2 da3 da4 solver its -
12 0.0 0.0 0.0 0.0MECH 1
$# option mshid -x +x -y +y -z +z
SALEFAC 1 0 1 0 0 0 0
*SET_SEGMENT_GENERAL
$ SID
$# sid da1 da2 da3 da4 solver its -
1 MECH
$ OPTION MSHID XMN XMX YMN YMX ZMN ZMX
$# option mshid -x +x -y +y -z +z
SALEFAC 1 1 1 1 1 1
*SET_SOLID_GENERAL
$ SID
$# sid solver
1MECH
$# option mshid imin imax jmin jmax kmin kmax
SALECPT 1 1 121 1 61 1 81
*SET_SOLID_GENERAL
$ SID
$# sid solver
2MECH
$# option mshid imin imax jmin jmax kmin kmax
SALECPT 1 1 2 1 61 1 81
*END
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删