找回密码
 注册
Simdroid-非首页
查看: 108|回复: 1

[耦合分析] 梁单元与实体单元的流固耦合

[复制链接]
发表于 2013-8-2 23:10:47 | 显示全部楼层 |阅读模式 来自 天津
如题,请教前辈,梁单元能否参加流固耦合,我自己建了个小例子,算是能算,但非常慢,不知道是不是我的设置有问题还是梁单元不能参与流固耦合?以下是我的部分K文件内容,请大家指点!


$# LS-DYNA Keyword file created by LS-PrePost 3.2 - 12Mar2012(08:25)
$# Created on Jul-31-2013 (12:44:17)
*KEYWORD  
*TITLE
$# title
FSI
*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 1.0000E-6         0         0     0.000         0
$#    ncpl      nbkt    imascl    checkr
         1        50         0     0.000
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen
         2         2         2         2
*CONTROL_HOURGLASS
$#     ihq        qh
         1  0.100000
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
  2.000000         0     0.000     0.000     0.000
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
     0.000  0.900000         0     0.000     0.000         0         0         0
$#  dt2msf   dt2mslc     imscl
     0.000         0         0
*DATABASE_GLSTAT
$#      dt    binary      lcur     ioopt
  0.010000         0         0         1
*DATABASE_MATSUM
$#      dt    binary      lcur     ioopt
  0.010000         0         0         1
*DATABASE_SECFORC
$#      dt    binary      lcur     ioopt
  0.010000         0         0         1
*DATABASE_SLEOUT
$#      dt    binary      lcur     ioopt
  0.010000         0         0         1
*DATABASE_FORMAT
$#   iform   ibinary
         0         0
*DATABASE_EXTENT_BINARY
$#   neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
         0         0         3         1         1         1         1         1
$#  cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat
         0         0         4         1         1         1         2         0
$# nintsld   pkp_sen      sclp    unused     msscl     therm    intout    nodout
         0         0  1.000000         0         0         0STRESS    STRESS
*DATABASE_CROSS_SECTION_SET_ID
$#    csid                                                                 title
         1                                                                     
$#    nsid      hsid      bsid      ssid      tsid      dsid        id     itype
         1         1         0         0         0         0         0         0
*PART
$# title
L_r
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         1         0         0         0         0         0
*SECTION_BEAM
$#   secid    elform      shrf   qr/irid       cst     scoor       nsm
         1         1  1.000000         2         1     0.000     0.000
$#     ts1       ts2       tt1       tt2     nsloc     ntloc
  0.020000  0.020000     0.000     0.000     0.000     0.000
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         1 7800.00002.0600E+11  0.300000 3.3500E+8 1.2000E+9     0.000
$#     src       srp        fs        vp
40.000000  5.000000  0.120000     0.000
*PART
$# title
C_r
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         2         2         0         0         0         0         0
*SECTION_BEAM
$#   secid    elform      shrf   qr/irid       cst     scoor       nsm
         2         1  1.000000         2         1     0.000     0.000
$#     ts1       ts2       tt1       tt2     nsloc     ntloc
  0.010000  0.010000     0.000     0.000     0.000     0.000
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         2 7800.00002.0600E+11  0.300000 2.3500E+8 1.2000E+9     0.000
$#     src       srp        fs        vp
40.000000  5.000000  0.150000     0.000
*PART
$# title
Concrete
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         3        72         0         1         0         0         0
*SECTION_SOLID
$#   secid    elform       aet
         3        11         0
*SET_PART_LIST_TITLE
all_rebar
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#    pid1      pid2      pid3      pid4      pid5      pid6      pid7      pid8
         1         2         0         0         0         0         0         0

*CONSTRAINED_LAGRANGE_IN_SOLID
$#   slave    master     sstyp     mstyp     nquad     ctype     direc     mcoup
         1         3         0         1         3         2         1         0
$#   start       end      pfac      fric    frcmin      norm   normtyp      damp
     0.0001.0000E+10  0.100000     0.000  0.300000         0         0     0.000
$#      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
*END
发表于 2013-10-22 21:07:31 | 显示全部楼层 来自 福建福州
Simdroid开发平台
您好,我模拟框架爆破后倒塌,用流固耦合,结果可以算完,框架也可以倒塌,但是结果显示框架上的混凝土穿透刚体地面,是那种无视地面的穿透,就是刚体地面压根没起到作用。         框架的混凝土是拉格朗日单元,刚体地面也是拉格朗日单元,流固耦合是否只是定义了流体与固体的接触作用,但是没有定义固体与固体的接触,怎么定义固体与固体的作用??
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 14:14 , Processed in 0.027772 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表