sph接触问题

huang_zhi_long · 发表于 2010-10-2 14:20:22

本帖最后由 huang_zhi_long 于 2010-10-2 14:21 编辑

侵蝕接觸
*CONTACT_ERODING_NODES_TO_SURFACE. E& {, d* R" J) |6 e
2,1,4,3% M" c1 I$ H3 }# U: C
0,0,0,0,0,1* D0 u# g8 h V- p a
1,1,0,0,1,1,1,14 P6 K6 Z* O+ H3 h- R; E5 ?2 A; m7 Q
1,1,1
10 a! Z+ Q5 z" ?; B

第一个接触视频输出如下d3plot

出视频如下现穿透现象
* v7 G+ R7 h( u3 j$ \! V
2.做侵蝕接觸! c4 r/ A4 I0 d2 X* Z4 Z
*CONTACT_ERODING_SURFACE_TO_SURFACE
2,1,4,3
0,0,0,0,0,1' p9 L) H( d7 a2 X4 k2 k9 o3 m
1,1,0,0,1,1,1,1
1,1,1 K+ _# u# M Z4 Q
1
出现下图错误，大家帮忙看看啊

nwpuzyq · 发表于 2010-10-2 14:20:23

我记得我用网格生成SPH粒子时也出现过这样的问题，后来用prepost直接生成粒子就好了

huang_zhi_long · 发表于 2010-10-2 14:20:59

我的是用有网格直接生成的sph，
$# LS-DYNA Keyword file created by LS-PREPOST 3.1 (Beta) - 19Aug2010(08:33)
$# Created on Oct-01-2010 (20:00:07)
*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
*CONTROL_ENERGY
$# hgen    rwen slnten    rylen
      2       2       1       1
*CONTROL_OUTPUT
$# npopt neecho nrefup iaccop    opifs ipnint ikedit iflush
      0       0       0       0    0.000       0    100    3000
$# iprtf ierode    tet10 msgmax ipcurv
      0       0       2       50       0
*CONTROL_SHELL
$#  wrpang    esort    irnxx istupd theory    bwc    miter    proj
20.000000       1       -1       1       2       2       1       0
$# rotascl intgrd lamsht cstyp6 tshell nfail1 nfail4 psnfail
  1.000000       0       0       1       0       0       0       0
$# psstupd irquad    cntco
      0       0       0
*CONTROL_SPH
$# ncbs    boxid       dt    idim memory    form    start    maxv
      1       11.0000E+20       0    150       0    0.0001.0000E+15
$# cont    deriv             iact
      0       0       0       0
*CONTROL_TERMINATION
$#  endtim endcyc    dtmin endeng endmas
  0.400000       0    0.000    0.000    0.000
*CONTROL_TIMESTEP
$#  dtinit tssfac    isdo tslimt    dt2ms    lctm    erode    ms1st
   0.000  0.700000       0    0.000    0.000       0       0       0
$#  dt2msf dt2mslc    imscl
   0.000       0       0
*DATABASE_GLSTAT
$#    dt binary    lcur    ioopt
  0.001000       0       0       1
*DATABASE_MATSUM
$#    dt binary    lcur    ioopt
  0.010000       0       0       1
*DATABASE_RCFORC
$#    dt binary    lcur    ioopt
  0.001000       0       0       1
*DATABASE_BINARY_D3PLOT
$#    dt    lcdt    beam    npltc psetid
  0.001000       0       0       0       0
$# ioopt
      0
*DATABASE_BINARY_D3THDT
$#    dt    lcdt    beam    npltc psetid
  0.001000       0       0       0       0
*DATABASE_BINARY_INTFOR
$#    dt    lcdt    beam    npltc psetid
  0.010000       0       0       0       0
*DATABASE_FORMAT
$# iform ibinary
      0       1
*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       1       2       1
$# nintsld pkp_sen    sclp unused    msscl    therm intout nodout
      0       0  1.000000       0       0       0STRESS STRESS
*BOUNDARY_PRESCRIBED_MOTION_RIGID
$#    pid    dof    vad    lcid       sf    vid    death    birth
      1       3       0       1  1.000000       01.0000E+28    0.000
*BOUNDARY_SPC_SET
$# nsid    cid    dofx    dofy    dofz    dofrx    dofry    dofrz
      3       0       1       1       1       1       1       1
*SET_NODE_LIST_TITLE
NODESET(SPC) 3
$#    sid    da1    da2    da3    da4 solver
      3    0.000    0.000    0.000    0.000MECH
$# nid1    nid2    nid3    nid4    nid5    nid6    nid7    nid8
   48441    48443    48444    48445    48446    48447    48448    48449
   48450    48451    48452    48453    48454    48455    48456    50142
   50143    50144    50145    50146    50147    50148    50149    50150
、、、、、、
   81513    81514    81515    81516    81517    81518    81519    81520
   81521    81522    81523    81524    81525    81526    81527    81528
   81529    81530    81531    81532    81533    81534    81535    81536
   81537    81538    81539    81540    81541    81542    81543       0
*SET_SEGMENT
$#    sid    da1    da2    da3    da4 solver
      1    0.000    0.000    0.000    0.000MECH
$#    n1       n2       n3       n4       a1       a2       a3       a4
   4590    14462    14463    4591    0.000    0.000    0.000    0.000
   4590    4591    14759    14760    0.000    0.000    0.000    0.000
   14537    14529    14544    14542    0.000    0.000    0.000    0.000
、、、、、
  15238    1 4756 15239 4726    332    0    0    0    0
15239    1 4727 15239 4756    127    0    0    0    0
*ELEMENT_SPH
$# nid    pid          mass
48441    3    0.0030000
48443    3    0.0030000
48444    3    0.0030000
48445    3    0.0030000
48446    3    0.0030000
、、、
  81637    3    0.0030000
81638    3    0.0030000
81639    3    0.0030000
81640    3    0.0030000
81641    3    0.0030000
*NODE
$# nid             x             y             z    tc    rc
   1 -277.1166992    -4.0000000    687.0463867    0    0
   2 -277.1166992    -10.5000000    687.0463867    0    0
   3 -282.1018677    1.0000000    686.6613770    0    0
   4 -294.0662231    1.0000000    685.7371826    0    0
   5 -298.9763489    -3.1357410    685.3579102    0    0
、、、、、、
81639 -140.0000000    -90.0000000    470.0000000    0    0
81640 -140.0000000    -90.0000000    480.0000000    0    0
81641 -140.0000000    -90.0000000    490.0000000    0    0
*END
*COMPONENT
$# clid color1 color2 color3 color4
      1  0.769000  0.004000  0.110000    0.000       0       0       0
$# name
Part 1
*COMPONENT_PART
$#    pid    clid
      1       1
$#    pid    clid
      3       1
*COMPONENT_END

huang_zhi_long · 发表于 2010-10-5 12:19:59

b] 3# nwpuzyq
我的这个问题我自己已经解决了，三楼的兄的，粒子的体积在哪里能设置？

nwpuzyq · 发表于 2010-10-5 15:56:33

我不知道你用的什么软件。我用的是prepost，只要输入密度，软件直接生成的粒子的质量。不过你要想算粒子体积的话可以自己计算，也很容易，只要知道模型的总体积除以粒子总数，然后用Excel或者UltraEdit编辑一下就行了。另外，你的问题是怎么解决的？希望交流一下

huang_zhi_long · 发表于 2010-10-6 10:48:33

5# nwpuzyq
我的QQ263048847

sccxhrq · 发表于 2013-4-9 21:05:18

请问楼主你这个问题是怎么解决的，我也出现这个问题了，希望楼主解决之后跟大家分享下啊，呵呵

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[接触碰撞] sph接触问题

最佳答案

相关帖子