【分享】接触面上的压力总和怎样得到

wuhanwangtao

第四版，欢迎指正！
!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

jacobi

wuhanwangtao 老师太强了！
再一次表示敬意！
谢谢！！

mondgao

不好意思，能不能大致讲一下您的思路？我刚学apdl，查了命令手册，但有些还是不明白，比如
*GET,stat1,ELEM,el,NMISC,19,1,2 !stat1, contact element status  
这里19，1，2是什么参数？
比较愚钝，多谢指点

hsmy

楼主能不能贴出命令流的解释阿，谢谢了

下午茶

是啊
只看命令流有点困难

crazywolf411

厉害

crazywolf411

厉害

lucychen

我最近做了一个电阻点焊预压力接触分析，在接触压力云图中有接触压力，可提取不出来是为什么

lucychen

我最近做了一个电阻点焊预压力接触分析，在接触压力云图中有接触压力，可提取不出来是为什么

galen_03

郁悶，怎麼看不懂

zero966

感谢!

hai00011

*GET,lel,ELEM,el,LENG
lel=sqrt(dx1*dx1+dy1*dy1)出现了两个lel,不知道怎么回事？
另外ael=(r1+r2)*lel ，这是表示接触面积吗？怎么得来的啊 ？
有没有人能回答这个问题啊？wuhanwangtao是不是已经不来这个论坛了？

zhouboshi

看不懂呀

suyanping333

我也想看接触压力，这个命令流好像看不太懂，能有解释就好了

battlestar

为什么楼主不用172单元的carea（nmisc,27,28）去计算接触面积ael呢？而是通过判断单元接触状态，去计算接触面积，有点不懂哦~
还有你的接触单元上节点的合力（cnfx，cnfy） 和 单元上的接触压力和切向力（pel，sel），好像前者没有在后面数据处理中使用，而后者用到了，如果是为了求某个方向合力的话，为何不用fsum呢？

strongsailor

楼主强大！！！！！！！！！！！！！！！！！