S-ALE造波技术:无板造波与边界流速法详解

https://ftp.lstc.com/anonymous/outgoing/hao/sale/models/waterlanding/

这是网络上的一个速度边界的SALE造波方法,包含k文。

waterlanding.gif

速度边界造波(应该也叫无板造波吧)

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


免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空