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

 作者:gjvvhh https://www.bilibili.com/read/cv27122162/ 出处:bilibili

(1条)
默认 最新
👍
评论 点赞
点赞 1 评论 1 收藏 8
关注