ansys分析之接触面上的压力总和.doc

45.接触面上的压力总和单面接触用于当一个物体的外表面与自身接触或和另一个物体的外表面接触时使用.单面接触是ANSYS/LS-DYNA 中最通用的接触类型，因为程序将搜索模型中的所有外表面，检查其间是否相互发生穿透.由于所有的外表面都在搜索范围内, 不需要定义接触面与目标面.在预先不知接触情况时，单面接触非常有用.相对于ANSYS 隐式分析 , ANSYS/LS-DYNA 的单面接触不会非常耗时。大多数冲击与碰撞问题需要定义单面接触当接触面之间的穿透超过接触单元厚度的40%时，单面接触自动释放接触，从而对下列问题造成潜在的问题 :1. 超薄部分2. 具有低刚度的软体3. 高速运动物体之间的接触单面接触在 ASCII rcforc 文件中不记录所有的接触反作用力. 如果需要接触反力, 可以使用点到面或面到面接触。有效的接触类型有 Single Surface, Automatic Single Surface, Automatic General, Eroding Single Surface和 Single Edge第三版，欢迎指正！!calculate mean contact pressure!file name : pressure.txt!command: /input,pressure,txt!pave is the mean contact pressureesel,s,type,3 !contact element 172*get,nelm,elem,0,count*get,el,elem,0,num,min*dim,pel,nelm,1*dim,sfel,nelm,1pfsum=0233pfsumx=0pfsumy=0sfsum=0sfsumx=0sfsumy=0asum=0asumx=0asumy=0pave=0 !the average contact pressurenumcont=0 !the number of contact elementscnfxsum=0cnfysum=0pi=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 its open!*GET,ael,ELEM,el,AREA!For axisymmetric model*GET,lel,ELEM,el,LENGnd1=nelem(el,1)nd2=nelem(el,2)sintheta=(ny(nd2)+uy(nd2)-ny(nd1)-uy(nd1)/lelcostheta=(nx(nd2)+ux(nd2)-nx(nd1)-ux(nd1)/lelr1=nx(nd1)+ux(nd1)r2=nx(nd2)+ux(nd2)!ael=pi*(r1+r2)*lelael=(r1+r2)*lel!*!calculate force in x and y direction*GET,lelx,ELEM,el,LPROJ,X*GET,lely,ELEM,el,LPROJ,Yaely=(r1+r2)*lelyaelx=(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 fxsum=cnfxsum+cnfxcnfysum=cnfysum+cnfy!*asum=asum+aelpfsum=pfsum+ael*pel( i )sfsum=sfsum+ael*sfel( i )234!sum of the forces in x and y directionsasumx=asumx+aelxasumy=asumy+aelypfsumx=pfsumx+ael*pel( i )*sinthetapfsumy=pfsumy+ael*pel( i )*costhetasfsumx=sfsumx+aelx*sfel( i )*costhetasfsumy=sfsumy+aely*sfel( i )*sintheta!*numcont=numcont+1*endifel=elnext(el)*enddo!average and total forcespave=pfsum/asumpfsum=pi*pfsumsfave=sfsum/asumsfsum=pi*sfsum!average and total forces in x and y directionspavex=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 pressureesel,s,type,3 !contact element 172*get,nelm,elem,0,count*get,el,elem,0,num,min*dim,pel,nelm,1*dim,sfel,nelm,1pfsum=0pfsumx=0pfsumy=0sfsum=0235sfsumx=0sfsumy=0asum=0asumx=0asumy=0pave=0 !the average contact pressurenumcont=0 !the number of contact elementscnfxsum=0cnfysum=0pi=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 its open!*GET,ael,ELEM,el,AREA!For axisymmetric model*GET,lel,ELEM,el,LENGnd1=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/lelcostheta=dx1/lelr1=nx(nd1)+ux(nd1)r2=nx(nd2)+ux(nd2)!ael=pi*(r1+r2)*lelael=(r1+r2)*lel!*!calculate force in x and y directionlelx=abs(dx1)lely=abs(dy1)aely=(r1+r2)*lelyaelx=(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 fxsum=cnfxsum+cnfxcnfysum=cnfysum+cnfy!*asum=asum+aelpfsum=pfsum+ael*pel( i )sfsum=sfsum+ael*sfel( i )236!sum of the forces in x and y directionsasumx=asumx+aelxasumy=asumy+aelypfsumx=pfsumx+ael*pel( i )*sinthetapfsumy=pfsumy+ael*pel( i )*costhetasfsumx=sfsumx+ael*sfel( i )*costhetasfsumy=sfsumy+ael*sfel( i )*sintheta!*numcont=numcont+1*endifel=elnext(el)*enddo!average and total forcespave=pfsum/asumpfsum=pi*pfsumsfave=sfsum/asumsfsum=pi*