【讨论】接触分析，简单的命令流如下：

bigbird

^^......我以為會很多人幫忙說........
先試試在SOLVE前ALLSEL一次

bigbird

^^...OK~~..基本上不難....只是ansys使用觀念問題.......
通常會沒有選擇到線的ERROR表示
有使用SELECT COMMEND..........
當使用xSEL時(x為K,N,A,V or CM).........
就已宣告接下來的動作......將只作用在被選擇的K,N,A,V or CM上.........
因此.....^^..接下來應該不用說了吧........
一個良好的習慣是
使用SELECT COMMEND後..........其後再加上ALLSEL.........
這樣懂了嗎^^

wuhanwangtao

我做了下
只有等效应力
无接触应力，磨察力
你的键模和约束有问题
你能不能把问题说具体？？？？？？？？？？？
我在帮你看看
图为：

wuhanwangtao

应力图为

wuhanwangtao

你的约束应该这样
六个颇面用对称约束，建易用菜单的方法加）
外套外圆用约束
csys,0  
nsel,s,loc,x  
d,all,ux  
asel,s,,,28  
nsla,,1  
d,all,all （！！！！！！！！！！！有问题，这不是对称约束！！！！！！！！！）
asel,s,,,29,35,6  
nsla,,1  
nrotat,all  
d,all,all  
asel,s,,,2  
nsla,,1  
cp,1,uy,all  
*get,nc,node,,num,min  
  
不知是不是我的ansys有问题
用plns 画接触压力，种说有问题，但用apdl有可提取数据
我做ANSYS5.7VM191，也一样，用plns 画接触压力，种说有问题，但用apdl有可提取数据
等效应力都可画出来，你试试。
你的问题有两个：
1）是约束
2）是没alls

wuhanwangtao

tu

wuhanwangtao

tu:

wuhanwangtao

关于约束有两种：
1）自身的约束，对称约束
你是提取实体的一部份，要用到对称约束对剖面约束
apply->displacement->symmetry B.C.->on area  选中六个剖面即可，它不让你选ux,uy
2）施加的约束
如限制中空圆柱棍子和高速轴横向移动和固定外套
要分开加。
  
ANSYS57中的VM191 不能显示摩察力，这可能与单元有关，48，49（点对面）单员不能显示摩察力
但169（面到面）可以。
如果你只是对系统分析，等效应力就可以了。
如果你想得出磨察力
可用编程提取
具体你可看王国强的接触分析！

wuhanwangtao

系统间的等效应力可理解为系统间的接触应力。
而且ansys只能给出每个节点的切向和法向力，但不能给个合力。
求磨察力，可通过和磨察力平横的力来求
也可通过编程来求
具体是：
先得到接触单元的状态
ETABLE,NSTAT,NMISC,1   
在求里
fn:法向力
fs:滑动力
具体看单元输出表。
  
EL Element Number Y Y  
  
FS1 Tangential force in element direction 1 (elastic or sliding) 2 2  
FS2 Tangential force in element direction 2 (elastic or sliding) 2 2  
FSLIM Friction limit force 2 2  
MU Coefficient of friction 2 2  
  
王国强的有限元分析（非常基础的那本），第130页的例子可看看。
要说的是求种的求磨察力和法向力，还没有一个好的方法。
你只想进行结构分析，等效应力够了，它是磨察力和法向力纵合的。
你要想具体得出种的磨察力和种的法向力，要自已编程。
我给你一个别人做过的东西

wuhanwangtao

接触面上的压力总和怎样得到
第三版，欢迎指正！
!calculate mean contact pressure
!file name : pressure.txt
!command: /input,pressure,txt
!pave is the mean contact pressure
esel,s,type,,3 !contact element 172
*get,nelm,elem,0,count
*get,el,elem,0,num,min
*dim,pel,,nelm,1
*dim,sfel,,nelm,1
pfsum=0
pfsumx=0
pfsumy=0
sfsum=0
sfsumx=0
sfsumy=0
asum=0
asumx=0
asumy=0
pave=0 !the average contact pressure
numcont=0 !the number of contact elements
cnfxsum=0
cnfysum=0
pi=4*atan(1)
*do,i,1,nelm
*GET,stat1,ELEM,el,NMISC,19,1,2 !stat1, contact element status
*if,stat1,ge,2,then !if stat1>=2, it is in close contact,else it's open
!*GET,ael,ELEM,el,AREA  
!For axisymmetric model
*GET,lel,ELEM,el,LENG  
nd1=nelem(el,1)
nd2=nelem(el,2)
sintheta=(ny(nd2)+uy(nd2)-ny(nd1)-uy(nd1))/lel
costheta=(nx(nd2)+ux(nd2)-nx(nd1)-ux(nd1))/lel
r1=nx(nd1)+ux(nd1)
r2=nx(nd2)+ux(nd2)
!ael=pi*(r1+r2)*lel
ael=(r1+r2)*lel
!*
!calculate force in x and y direction
*GET,lelx,ELEM,el,LPROJ,X  
*GET,lely,ELEM,el,LPROJ,Y  
aely=(r1+r2)*lely
aelx=(r1+r2)*lelx
!*
*GET,pel( i ),ELEM,el,SMISC, 5,1,2 !get normal contact pressure of contact element el.
*GET,sfel( i ),ELEM,el,SMISC, 3,4 !get tangential contact stress of contact element el.
*GET,cnfx,ELEM,el,NMISC, 21 !get Contact element el force-x component.
*GET,cnfy,ELEM,el,NMISC, 22 !get Contact element el force-y component.
cnfxsum=cnfxsum+cnfx
cnfysum=cnfysum+cnfy
!*
asum=asum+ael
pfsum=pfsum+ael*pel( i )
sfsum=sfsum+ael*sfel( i )
!sum of the forces in x and y directions
asumx=asumx+aelx
asumy=asumy+aely
pfsumx=pfsumx+ael*pel( i )*sintheta
pfsumy=pfsumy+ael*pel( i )*costheta
sfsumx=sfsumx+aelx*sfel( i )*costheta
sfsumy=sfsumy+aely*sfel( i )*sintheta
!*
numcont=numcont+1  
*endif
el=elnext(el)
*enddo
!average and total forces
pave=pfsum/asum
pfsum=pi*pfsum
sfave=sfsum/asum
sfsum=pi*sfsum
!average and total forces in x and y directions
pavex=pfsumx/asumx ！平均法向接触力x向分量
pfsumx=pi*pfsumx ！法向接触力x向分量
pavey=pfsumy/asumy ！平均法向接触力y向分量
pfsumy=pi*pfsumy ！法向接触力y向分量
sfavex=sfsumx/asumx ！平均摩擦力x向分量
sfsumx=pi*sfsumx ！摩擦力x向分量
sfavey=sfsumy/asumy ！平均摩擦力y向分量
sfsumy=pi*sfsumy ！摩擦力y向分量
!*
cnfxsum=cnfxsum ！x向接触力(摩擦力和法向接触力的x向分量)
cnfysum=cnfysum ！y向接触力（摩擦力和法向接触力的y向分量）
  
第四版，欢迎指正！
!calculate mean contact pressure
!file nameressure.txt
!command: /input,pressure,txt
!pave is the mean contact pressure
esel,s,type,,3 !contact element 172
*get,nelm,elem,0,count
*get,el,elem,0,num,min
*dim,pel,,nelm,1
*dim,sfel,,nelm,1
pfsum=0
pfsumx=0
pfsumy=0
sfsum=0
sfsumx=0
sfsumy=0
asum=0
asumx=0
asumy=0
pave=0 !the average contact pressure
numcont=0 !the number of contact elements
cnfxsum=0
cnfysum=0
pi=4*atan(1)
*do,i,1,nelm
*GET,stat1,ELEM,el,NMISC,19,1,2 !stat1, contact element status
*if,stat1,ge,2,then !if stat1>=2, it is in close contact,else it's open
!*GET,ael,ELEM,el,AREA  
!For axisymmetric model
*GET,lel,ELEM,el,LENG  
nd1=nelem(el,1)
nd2=nelem(el,2)
dx1=nx(nd2)+ux(nd2)-nx(nd1)-ux(nd1)
dy1=ny(nd2)+uy(nd2)-ny(nd1)-uy(nd1)
lel=sqrt(dx1*dx1+dy1*dy1)
sintheta=dy1/lel
costheta=dx1/lel
r1=nx(nd1)+ux(nd1)
r2=nx(nd2)+ux(nd2)
!ael=pi*(r1+r2)*lel
ael=(r1+r2)*lel
!*
!calculate force in x and y direction
lelx=abs(dx1)
lely=abs(dy1)
aely=(r1+r2)*lely
aelx=(r1+r2)*lelx
!*
*GET,pel( i ),ELEM,el,SMISC, 5,1,2 !get normal contact pressure of contact element el.
*GET,sfel( i ),ELEM,el,SMISC, 3,4 !get tangential contact stress of contact element el.
*GET,cnfx,ELEM,el,NMISC, 21 !get Contact element el force-x component.
*GET,cnfy,ELEM,el,NMISC, 22 !get Contact element el force-y component.
cnfxsum=cnfxsum+cnfx
cnfysum=cnfysum+cnfy
!*
asum=asum+ael
pfsum=pfsum+ael*pel( i )
sfsum=sfsum+ael*sfel( i )
!sum of the forces in x and y directions
asumx=asumx+aelx
asumy=asumy+aely
pfsumx=pfsumx+ael*pel( i )*sintheta
pfsumy=pfsumy+ael*pel( i )*costheta
sfsumx=sfsumx+ael*sfel( i )*costheta
sfsumy=sfsumy+ael*sfel( i )*sintheta
!*
numcont=numcont+1  
*endif
el=elnext(el)
*enddo
!average and total forces
pave=pfsum/asum
pfsum=pi*pfsum
sfave=sfsum/asum
sfsum=pi*sfsum
!average and total forces in x and y directions
pavex=pfsumx/asumy
pfsumx=pi*pfsumx
pavey=pfsumy/asumx
pfsumy=pi*pfsumy
sfavex=sfsumx/asumy
sfsumx=pi*sfsumx
sfavey=sfsumy/asumx
sfsumy=pi*sfsumy
!*
cnfxsum=cnfxsum
cnfysum=cnfysum