跪求接触分析命令流错误解决方案

那个秋天那地方

觉得应该是一个比较简单的问题啊！但算不出来！头皮发麻啊兄弟们！！给点建议！！！以下是命令流：
  fini
  /clear
  plate_r1=50
  plate_r2=60
  pin_r=45
  t=10
  l=100
  /prep7
  et,1,plane183
  et,3,targe169
  et,4,conta172,,,,,1
  r,1,,,1
  mp,ex,1,2e5
  mp,prxy,1,0.3
  mp,mu,0.15
  cyl4,,,plate_r2,-90,,90
  rectng,0,plate_r2,0,l
  aadd,all
  cyl4,,,plate_r1
  asba,3,1
  wpoffs,,pin_r-plate_r1
  csys,4
  cyl4,,,pin_r,-90
  
   asel,s,,,1
   lesize,2,,,20
   lesize,3,,,20
   lesize,4,,,20
   mshape,0
   mshkey,1,2d
   amesh,1
   allsel
   
   wpcsys,,0
   wprota,,90
   asbw,2
   lesize,10,,,20
   lesize,8,,,20
   lesize,1,,,20
   lesize,11,,,4
   lesize,7,,,4
   lesize,5,,,34
   lesize,12,,,4
   lesize,9,,,10
   lesize,6,,,10
   mshape,0
   mshkey,1,2d
   amesh,3
   mshkey,2,2d
   amesh,4
   allsel
   
   
   csys,1
   lsel,s,,,10
   nsll,s,1
   real,1
   type,3
   esurf
   allsel
   
    csys,4
    wpcsys,,0
   wpoffs,,pin_r-plate_r1
   wpstyl,,,,,,1
   lsel,s,,,2
   nsll,s,1
   real,1
   type,4
   esurf
   allsel
   /solu
    dl,6,,all
    dl,9,,ux
    dl,12,,ux
    dl,11,,ux
    dl,4,,ux
    lsel,s,,,3
  nsll,s,1
  sf,all,pres,1e4
  cp,1,uy,all
  alls
   lnsrch,on
  autots,on
  nsubst,20,100,5
  outres,all,all
  nropt,full,,off
  solve
   

sunxueshuai

你的命令流写得有点乱，很多没用的句子。例如：定义目标面和接触面前的csys语句就没有任何意义，划分网格的方法也不够简洁。不过这都不是主要问题，我在你的命令流基础上修改了接触参数，如下：
fini
  /clear
  /prep7
  plate_r1=50
  plate_r2=60
  pin_r=45
  t=10
  l=100
  et,1,plane183,,,3
  et,3,targe169
  et,4,conta172,,,,,1
  r,1,,,t           ！上面的t我猜测是板厚，在此加上了
  mp,ex,1,2e5
  mp,prxy,1,0.3
  mp,mu,0.15
  cyl4,,,plate_r2,-90,,90
  rectng,0,plate_r2,0,l
  aadd,all
  cyl4,,,plate_r1
  asba,3,1
  wpoffs,,pin_r-plate_r1
  csys,4
  cyl4,,,pin_r,-90
  
   asel,s,,,1
   lesize,2,,,20
   lesize,3,,,20
   lesize,4,,,20
   mshape,0
   mshkey,1,2d
   amesh,1
   allsel
   
   wpcsys,,0
   wprota,,90
   asbw,2
   lesize,10,,,20
   lesize,8,,,20
   lesize,1,,,20
   lesize,11,,,4
   lesize,7,,,4
   lesize,5,,,34
   lesize,12,,,4
   lesize,9,,,10
   lesize,6,,,10
   mshape,0
   mshkey,1,2d
   amesh,3
   mshkey,2,2d
   amesh,4
   allsel
   
   
   csys,1
   R,3,,,0.1,0.1,,     ！这里是关键，你定的参数有问题
   lsel,s,,,10
   nsll,s,1
   esln,s,0
   real,3
   type,3
   esurf
   allsel
   
    csys,4
    wpcsys,,0
   wpoffs,,pin_r-plate_r1
   wpstyl,,,,,,1
   lsel,s,,,2
   nsll,s,1
   esln,s,0
   real,3
   type,4
   esurf
   allsel
   /solu
    dl,6,,all
    dl,9,,ux
    dl,12,,ux
    dl,11,,ux
    dl,4,,ux
    lsel,s,,,3
  nsll,s,1
  sf,all,pres,1e4
  cp,1,uy,all
  alls
  pred,on
  autots,on
  nsubst,100,1000,10  ！我增大了迭代数，不过这个不修改应该也可以
  outres,all,all
  nropt,full,,off
  solve
最后能算出结果，你可以试试

那个秋天那地方

sunxueshuai 发表于 2013-3-26 17:23
你的命令流写得有点乱，很多没用的句子。例如：定义目标面和接触面前的csys语句就没有任何意义，划分网格的 ...

嗯 是有点乱 嘿嘿  嗯 其实那个t是原来我建模为实体的厚度，用面来建模的时候忘去掉了，现在看来是我设置参数有问题

sunxueshuai

网格划分一般优先采用映射网格划分，这样生成的网格比较规则。不过面映射网格划分需要满足一定的条件：必须是3条边或4条边组成的面，所以面4不能采用，不过可以用lcomb进行线合并使之满足条件。你的划分方法是对的，我觉得不够简洁是因为映射划分只需要用lesize定义任意两条相邻边的单元尺寸就行了，对边默认划分为相同数目的单元。你可以查看一些映射网格划分的资料。