切换到宽版

 找回密码
 注册
Simdroid-非首页
Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)»论坛 › ..:: CAE专区 › A02：LS-DYNA › dyna流固耦合 出现穿透现象
返回列表 发新帖
查看: 597|回复: 17

[耦合分析] dyna流固耦合 出现穿透现象

[复制链接]
gujianjian
发表于 2011-7-12 09:37:41 | 显示全部楼层 |阅读模式 来自 湖南长沙
用dyna模型土的切削，以前也有做成功的k文件，但是这次做的时候，设置的都一样，但是每次都发生穿透，土没有发生变形和流动，一直在没有变化，但是应力还有，请问怎么解决！感谢各位大哥哥，大姐姐的帮忙，本人不胜感激！
不胜感激, 成功, 姐姐, 模型

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

×

相关帖子
回复

使用道具 举报

wyc412721
发表于 2011-7-12 09:42:02 | 显示全部楼层 来自 陕西咸阳
Simdroid开发平台
土没有变形和流动，是用流固耦合还是定义接触呢？
回复 赞 不支持

使用道具 举报

弓长三石
发表于 2011-7-12 09:45:45 | 显示全部楼层 来自 江苏南京
在ls-prepost后处理器中的选中Selpart中的fluid试试看

评分

1

查看全部评分

回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-12 14:42:54 | 显示全部楼层 来自 湖南长沙
用的事流固耦合  没有定义接触  请各位大侠帮帮忙！
回复 赞 不支持

使用道具 举报

wyc412721
发表于 2011-7-12 15:28:52 | 显示全部楼层 来自 陕西咸阳
是不是土壤和空气共节点，用ALE算法？
回复 赞 不支持

使用道具 举报

fanny2324
发表于 2011-7-12 16:05:05 | 显示全部楼层 来自 湖北武汉
本帖最后由 fanny2324 于 2011-7-12 16:06 编辑

介绍穿透
http://forum.simwe.com/viewthread.php?tid=844423&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D555
回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-12 17:15:55 | 显示全部楼层 来自 湖南长沙
问题补充   用的是欧拉算法  定义的是LAGRAGE-IN-SOLID，请问下   共节点和LAGRAGE-IN-SOLID选其中一个就行了
回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-12 17:20:09 | 显示全部楼层 来自 湖南长沙
共节点的话   是不是只要画的网格一样就行了    如果是这样的  话  那我后来就是设置的同样的网格  但是穿透依旧
回复 赞 不支持

使用道具 举报

wyc412721
发表于 2011-7-12 18:09:57 | 显示全部楼层 来自 陕西咸阳
K文件不大的话就上传一下，方法查看
回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-13 09:01:12 | 显示全部楼层 来自 湖南长沙
k文件如下，请大侠费心！
$# LS-DYNA Keyword file created by LS-PREPOST 3.0 - 05Aug010(16:13)
$# Created on Jul-11-2011 (21:31:05)
*KEYWORD  
*TITLE
IUnits                                                                          
*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.000     0.000  1.000000  0.001000         0         0     0.000         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
  0.300000         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.001500         1         0         1
*DATABASE_MATSUM
$#      dt    binary      lcur     ioopt
  0.001500         1         0         1
*DATABASE_RCFORC
$#      dt    binary      lcur     ioopt
  0.001500         1         0         1
*DATABASE_BINARY_D3DUMP
$#      dt      lcdt      beam     npltc    psetid
2.5000E+5         0         0         0         0
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid
  0.003000         0         0         0         0
$#   ioopt
         0
*DATABASE_BINARY_D3THDT
$#      dt      lcdt      beam     npltc    psetid
  0.015000         0         0         0         0
*DATABASE_FORMAT
$#   iform   ibinary
         0         0
*DATABASE_EXTENT_BINARY
$#   neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
         0         0         3         1         0         0         0         0
$#  cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat
         0         0         4         0         0         0         0         0
$# nintsld   pkp_sen      sclp    unused     msscl     therm    intout    nodout
         0         0  1.000000         0         0         0STRESS    STRESS   
*DATABASE_FSI
$#      dt
  0.001000
$#dbsfi_id       sid     stype      swid    convid
         1         1         1         0         0
*BOUNDARY_PRESCRIBED_MOTION_RIGID
         1         1         0         1  1.000000         01.0000E+28     0.000
*BOUNDARY_SPC_SET
         1         0         1         1         1         1         1         1
*SET_NODE_LIST
         1     0.000     0.000     0.000     0.000MECH      
     81991     81992     81993     81994     81995     81996     81997     81998
     81999     82000     82001     82002     82003     82004     82005     82006
     82007     82008     82009     82010     82011     82012     82013     82014
     82015     82016     82017     82018     82019     82020     82021     82022
     82023     82024     82025     82026     82027     82028     82029     82030
     82031     82032     82033     82034     82035     82036     82037     82038
.................................
     85303     85304     85305     85306     85307     85308     85309     85310
     85311         0         0         0         0         0         0         0
*BOUNDARY_SPC_SET
         2         0         0         0         1         1         1         0
*SET_NODE_LIST
         2     0.000     0.000     0.000     0.000MECH      
     89134     89135     89136     89137     89138     89139     89140     89141
     89142     89143     89144     89145     89146     89147     89148     89149
     89150     89151     89152     89153     91554     91555     91556     91557
     91558     91559     91560     91561     91562     91563     91564     91565
     91566     91567     91568     91569     91570     91571     91572     91573
     91574     91575     91576     91577     91578     91579     91580     91581
     91582     91583     91584     91585     91586     91587     91588     91589
     91590     91591     91592     91593     91594     91595     91596     91597
     91598     91599     91600     91601     91602     91603     91604     91605
...................................
    275986    275987    275988    275989    275990    275991    275992    275993
    275994         0         0         0         0         0         0         0
*BOUNDARY_SPC_SET
         3         0         0         0         1         1         1         0
*SET_NODE_LIST
         3     0.000     0.000     0.000     0.000MECH      
     89154     89155     89156     89157     89158     89159     89160     89161
     89162     89163     89164     89165     89166     89167     89168     89169
     89170     89171     89172     89173     89174     89175     89176     89177
......................
*PART
$# title
Part          1 for Mat         1 and Elem Type         1                       
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         1         0         0         0         0         0
*SECTION_SOLID
$#   secid    elform       aet
         1         1         0
*MAT_RIGID
$#     mid        ro         e        pr         n    couple         m     alias
         1 7800.00002.1000E+11  0.300000     0.000     0.000     0.000         
$#     cmo      con1      con2
  1.000000         5         7
$# lco or a1      a2        a3        v1        v2        v3
     0.000     0.000     0.000     0.000     0.000     0.000
*PART
$# title
Part          2 for Mat       147 and Elem Type         1                       
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2       147       147         0         0         0         0         0
*SECTION_SOLID_ALE
$#   secid    elform       aet
       147        12         1
$#    afac      bfac      cfac      dfac     start       end     aafac
     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_FHWA_SOIL
$#     mid        ro     nplot    spgrav    rhowat        vn    gammar    intrmx
       147 1850.0000         3  2.650000 1000.0000  1.100000     0.000         0
$#       k         g    phimax      ahyp       coh     eccen        an        et
1.7560E+7 5.8500E+6  0.314000 160.00000 6200.0000  0.700000     0.000     0.000
$#   mcont      pwd1     pwksk      pwd2    phires      dint      vdfm    damlev
  0.340000     0.000     0.000     0.000  0.001000 1.0000E-5 6.0000E-8     0.000
$#  epsmax
     0.000
*PART
$# title
Part          3 for Mat       148 and Elem Type         1                       
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3       147       147         0         0         0         0         0
*INITIAL_VOID_PART
$#     pid
         3
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         1         0  1.000000  1.000000     0.000     0.000         0
$#                a1                  o1
               0.000         -10.0000000
           3.0000000         -10.0000000
*SET_PART_LIST
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH      
$#    pid1      pid2      pid3      pid4      pid5      pid6      pid7      pid8
         2         3         0         0         0         0         0         0
*CONSTRAINED_LAGRANGE_IN_SOLID
         1         1         1         0         3         4         2         1
     0.0001.0000E+10  0.100000  0.200000  0.300000         0         0     0.000
     0.000     0.000     0.000         0     0.000         0         0         0
         0         0         0         0     0.000         0     0.000
*ELEMENT_SOLID
$#   eid     pid      n1      n2      n3      n4      n5      n6      n7      n8
       1       1     657     424     452     439   10580    5221    5865    5566
       2       1   10580    5221    5865    5566   10581    5222    5866    5567
.....................
  297423       3  276032  276796  386908  386713  276013  276795  386909  386714
  297424       3  276013  276795  386909  386714  276736  276834  270593  270598
*NODE
$#   nid               x               y               z      tc      rc
       1      -0.0325454           0.000       0.0120624       0       0
       2      -0.0283121           0.000      -0.0016514       0       0
.......................

  386907      -0.0769101       0.1450000      -0.0360000       0       0
  386908      -0.0769101       0.1500000      -0.0360000       0       0
  386909      -0.0769101       0.1550000      -0.0360000       0       0
*END
*COMPONENT
         1  0.769000  0.004000  0.110000     0.000         0         0         0
Part 1                          
*COMPONENT_PART
         1         1
         2         1
         3         1
*COMPONENT_END
回复 赞 不支持

使用道具 举报

ywh362568
发表于 2011-7-13 20:33:46 | 显示全部楼层 来自 浙江嘉兴
如果不涉及机密的话，最好直接传k文件，这样可以运行起来，更能发现问题所在
回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-14 09:26:14 | 显示全部楼层 来自 湖南长沙
K文件太大了  上传不了啊 几十兆
回复 赞 不支持

使用道具 举报

wyc412721
发表于 2011-7-17 20:33:39 | 显示全部楼层 来自 陕西西安
帮忙上传一下附件

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

×
回复 赞 不支持

使用道具 举报

gujianjian
 楼主| 发表于 2011-7-18 08:25:18 | 显示全部楼层 来自 湖南长沙
请各位帮忙解决下问题  k文件已经上传   感谢版主的帮助
回复 赞 不支持

使用道具 举报

conlor
发表于 2011-9-30 17:23:27 | 显示全部楼层 来自 新加坡
本帖最后由 conlor 于 2011-9-30 17:24 编辑

我遇见和你一模一样的问题，土壤用欧拉，和固体耦合，也是有力，但是穿透，土壤看不见变形。

我的模型时对称的，使用整个模型建模的时候，无穿透，

当简化成四分之一模型的时候（在对称面设置对称边界条件），就发生这样的穿透现象。苦闷不已。
回复 赞 不支持

使用道具 举报

gaojing52402130
发表于 2011-10-15 12:29:54 | 显示全部楼层 来自 内蒙古呼和浩特
gujianjian 发表于 2011-7-13 09:01
k文件如下，请大侠费心！
$# LS-DYNA Keyword file created by LS-PREPOST 3.0 - 05Aug010(16:13)
$# Creat ...

你好，我有个问题跟你的问题很相似，还没有你的问题解决的好，我看到你在k文件中使用了*BOUNDARY_PRESCIBED_MOTIOM_RIGID这个关键字了，还有那个VOID_PART那个，我现在想实现一个边界流入一股流体，并且与流体域中的固体耦合，这个关键字必须配合solid12使用么？
回复 赞 不支持

使用道具 举报

护韶
发表于 2013-10-19 20:06:07 | 显示全部楼层 来自 福建福州
gujianjian 发表于 2011-7-12 14:42
用的事流固耦合  没有定义接触  请各位大侠帮帮忙！

我也是同样的问题。。。
回复 赞 不支持

使用道具 举报

flairy1988
发表于 2013-11-23 22:57:03 | 显示全部楼层 来自 天津
弓长三石 发表于 2011-7-12 09:45
在ls-prepost后处理器中的选中Selpart中的fluid试试看

猫猫，呵呵，没想到在论坛里看见你了
回复 赞 不支持

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

GMT+8, 2024-5-7 04:30 , Processed in 0.054205 second(s), 18 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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