Taylor撞击用Euler算法时，Taylor杆不边缘不流动？

heliling

一直在调整*Constraint_lagrangian_in_solid和Control_ALE里的参数，但是都没有任何改善。附k文件如下，请高手看看倒底哪个参数设置不当？

heliling

*CONTROL_ALE
$#     dct      nadv      meth      afac      bfac      cfac      dfac      efac
         2         1         2 -1.000000     0.000     0.000     0.000     0.000
$#   start       end     aafac     vfact      prit       ebc      pref   nsidebc
     0.0001.0000E+20  1.000000     0.000         0         0     0.000         0

*CONSTRAINED_LAGRANGE_IN_SOLID
$#   slave    master     sstyp     mstyp     nquad     ctype     direc     mcoup
         1         1         1         0         0         4         2         0
$#   start       end      pfac      fric    frcmin      norm   normtyp      damp
     0.0001.0000E+10  0.100000     0.000  0.500000         0         0         0
$#      cq      hmin      hmax     ileak     pleak   lcidpor     nvent  blockage
     0.000     0.000     0.000         0  0.010000         0         0         0
$#  iboxid   ipenchk   intforc   ialesof    lagmul    pfacmm      thkf
         0         0         0         0     0.000         0     0.000

1# heliling

heliling

2# heliling

*MAT_JOHNSON_COOK
$#     mid        ro         g         e        pr       dtf        vp
         2 7850.00008.0000E+10     0.000     0.000     0.000     0.000
$#       a         b         n         c         m        tm        tr      epso
1.0000E+9 5.1000E+8  0.260000  0.007000  1.030000 1811.0000 300.00000 1.0000E-6
$#      cp        pc     spall        it        d1        d2        d3        d4
477.00000     0.000  2.000000     0.000  3.000000  3.440000 -2.210000  0.002000
$#      d5
  0.610000
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         2 4680.0000  1.330000     0.000     0.000  1.670000  0.470000     0.000
$#      v0
     0.000


杆用JC模型模拟，钢材

wyc412721

泰勒杆撞击一般不用流固耦合算法呀，楼主想要研究什么呢

heliling

我在学习流固耦合，自己设计的算例，用纯拉格朗日、ALE和纯欧拉分别计算taylor撞击，看看有何不同





4# wyc412721