快乐神仙 发表于 2006-11-24 21:20

救急：一个难办的接触问题，大家帮帮忙吧！

我现在做了一个钻动的钻头钻岩石的受力分析，岩石固定，钻头受7t的压力和200r/min的钻速。我定义钻头为目标面，岩石为接触面，采用3D面－面接触。由于模型比较复杂，网格划分上我一直没有好的办法，四个刮刀长度不一，还要想办法只允许钻头在Y方向有位移，在X和Z方向上不会偏移。结果老是不收敛。下边是我的模型的命令流。
/BATCH
/COM,ANSYS RELEASE9.0    UP20041104       13:50:55    11/24/2006
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE
/COM,ANSYS RELEASE9.0    UP20041104       13:32:57    11/24/2006
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE
/PREP7
!*
ET,1,SOLID92
!*
ET,2,TARGE170   
!*
ET,3,CONTA174   
!*
KEYOPT,1,5,2
KEYOPT,1,6,0
!*
KEYOPT,2,1,1
KEYOPT,2,2,0
KEYOPT,2,3,0
!*
KEYOPT,2,4,111111   
KEYOPT,2,5,0
!*
KEYOPT,3,1,0
KEYOPT,3,2,0
KEYOPT,3,4,0
KEYOPT,3,5,3
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,1   
KEYOPT,3,11,0   
KEYOPT,3,12,0   
!*
R,1,,,30000, 0.1,0.1,   
RMORE,0.1,0.01,, ,,
RMORE,,,, ,,
RMORE,,,, ,,
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.1e6
MPDATA,PRXY,1,,0.0691   
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,1,,0.2
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2.6e4
MPDATA,PRXY,2,,0.24
TB,DP,2,,,
TBMODIF,1,1,3   
TBMODIF,1,2,45
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
K,1,0,0,,   
K,2,44.56,-11.94,,
K,3,107.5,36.36,,   
K,4,0,36.36,,   
K,5,0,20,,
K,6,57.18,25,,
K,7,80,49.41,,
K,8,80,180,,
K,9,0,180,,
K,10,107.5,60,,
LSTR,       1,       2
LSTR,       3,      10
LSTR,       3,       4
LSTR,       4,       1
LSTR,       5,       6
LSTR,       7,       8
LSTR,       8,       9
LSTR,       9,       5
!*
L2TAN,5,6   
!*
L2TAN,1,2   
LDELE,       2
KDELE,      10
LPLOT   
FLST,2,5,4
FITEM,2,8   
FITEM,2,5   
FITEM,2,9   
FITEM,2,6   
FITEM,2,7   
AL,P51X
FLST,2,4,4
FITEM,2,4   
FITEM,2,1   
FITEM,2,10
FITEM,2,3   
AL,P51X
FLST,2,1,5,ORDE,1   
FITEM,2,1   
FLST,8,2,3
FITEM,8,9   
FITEM,8,5   
VROTAT,P51X, , , , , ,P51X, ,360, ,
LPLOT   
FLST,2,1,5,ORDE,1   
FITEM,2,2   
FLST,8,2,3
FITEM,8,4   
FITEM,8,1   
VROTAT,P51X, , , , , ,P51X, ,360, ,
LPLOT   
/USER,1   
/VIEW,1,0.318376480369    ,0.341225726606E-01,0.947350023372
/ANG,   1,   1.48108757190
/REPLO
/VIEW,1,0.433742317120    ,0.459994397550    ,0.774772712840
/ANG,   1,   3.38417201420
/REPLO
KWPAVE,       4
wprot,0,-90
/VIEW,1, -0.864640048957E-02,0.208684905800    ,0.977944706949
/ANG,   1,-7.04404814856
/REPLO
/FOC,   1,-9.70925176500    ,   42.4445315627    ,   8.67878127022
/REPLO
CSYS,4
BLC5, , ,500,500,-200   
FLST,5,4,6,ORDE,2   
FITEM,5,1   
FITEM,5,-4
VSEL,U, , ,P51X
VPLOT   
FLST,2,4,6,ORDE,2   
FITEM,2,5   
FITEM,2,-8
VADD,P51X   
FLST,3,1,6,ORDE,1   
FITEM,3,10
VGEN,2,P51X, , , , , , ,0      
VSBV,       9,      10   
ALLSEL,ALL   
VSEL,U, , ,       6   
ALLSEL,ALL   
VSEL,S, , ,       6   
ALLSEL,ALL   
VSEL,U, , ,       6      
KWPAVE,      24
wpoff,-10   
wprot,0,-90
CSYS,4
BLC5, , ,20,160,-107.5
FLST,2,2,6,ORDE,2   
FITEM,2,5   
FITEM,2,7   
VINV,P51X   
VPLOT   
VOFFST,56,80, ,
FLST,2,2,6,ORDE,2   
FITEM,2,5   
FITEM,2,8   
VADD,P51X   
VPLOT   
FLST,3,2,6,ORDE,2   
FITEM,3,1   
FITEM,3,3   
VGEN,2,P51X, , , , , , ,0   
VSBV,       7,       8
VPLOT   
CSYS,5
FLST,3,1,6,ORDE,1   
FITEM,3,9   
VGEN,4,P51X, , , ,90, , ,0
VPLOT   
KWPAVE,       1
wprot,0,-90,0   
CSYS,4
wpoff,0,0,-50   
BLC5, , ,20,60,100
VSBV,       7,      11
VPLOT   
BLC5, , ,60,40,100
VSBV,       8,       7
VPLOT   
BLC5, , ,70,70,100
VSBV,       9,       7
VPLOT   
/VIEW,1,0.838302860588E-01, -0.993155305187    , -0.813327911660E-01
/ANG,   1,   176.871680595
/REPLO
NUMCMP,ALL
ALLSEL,ALL
请哪位大侠给我划分一下合理的网格，施加力后给我调试一下，这是我的毕业论文，真是急死我了，大家帮帮忙吧！
我的邮箱weixiuyan1981@126.com,谢谢啦！

[ 本帖最后由 快乐神仙 于 2006-11-28 09:29 编辑 ]

快乐神仙 发表于 2006-11-24 21:26

加力时只要下边岩石固定，钻头顶端施加一个7吨的均匀压力，在对钻头施加角位移或是角速度，同时要考虑钻头不会在XZ面上偏移，最终要考察的是钻头的刮刀底面的受力情况。
大家帮帮忙吧，谢谢啦！

快乐神仙 发表于 2006-11-26 15:43

自己顶一下，请高手赐教啊！焦急等待中。。。。。。

快乐神仙 发表于 2006-11-28 09:31

版主帮忙啊！
页: [1]
查看完整版本: 救急：一个难办的接触问题，大家帮帮忙吧！