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

快乐神仙

我现在做了一个钻动的钻头钻岩石的受力分析，岩石固定，钻头受7t的压力和200r/min的钻速。我定义钻头为目标面，岩石为接触面，采用3D面－面接触。由于模型比较复杂，网格划分上我一直没有好的办法，四个刮刀长度不一，还要想办法只允许钻头在Y方向有位移，在X和Z方向上不会偏移。结果老是不收敛。下边是我的模型的命令流。
/BATCH  
/COM,ANSYS RELEASE  9.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 RELEASE  9.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 编辑 ]

快乐神仙

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

快乐神仙

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

快乐神仙

版主帮忙啊！