ansys齿轮接触分析命令流.docx
《ansys齿轮接触分析命令流.docx》由会员分享,可在线阅读,更多相关《ansys齿轮接触分析命令流.docx(8页珍藏版)》请在冰豆网上搜索。
ansys齿轮接触分析命令流
ansys齿轮接触分析命令流
/prep7
!
单位是n,和mm
!
先建立两个齿轮模型(没又齿根过渡线,只适合基圆半径小于齿根圆的情况)
m=4!
齿轮模数
z1=45!
齿轮齿数
pi=3.1415926
ang=20!
分度圆上的压力角
ha_c=1!
齿顶高系数
c_c=0.25!
顶隙系数
ha=ha_c*m!
齿顶高
hf=(ha_c+c_c)*m!
齿根高
d=m*z1!
分度圆直径
db=d*cos(ang*pi/180)!
基圆直径
da=d+2*ha!
齿顶圆直径
df=d-2*hf!
齿根圆直径
x=0
s=pi*m/2+2*x*m*tan(ang*pi/180)!
分度圆齿厚theta_s=tan(ang*pi/180)-ang*pi/180
fai_s=s/d
alfa_a=acos(db/da)*180/pi!
齿顶圆压力角(角度表示)alfa_f=acos(db/df)*180/pi!
齿根圆压力角(角度表示)dfr=0.38*m!
齿根圆角半径
!
b=0.012!
齿宽
deata_ang=360/z1!
齿轮两齿间的夹角
*dim,alfa,array,10
*dim,rk,array,10
*dim,theta,array,10
*dim,fai,array,10
*do,i,1,10
alfa(i)=(alfa_f+((alfa_a-alfa_f)/9)*(i-1))*pi/180!
每条渐开线上生成10各关键点所
对应的压力角大小(用弧表示)
*enddo
*do,i,1,10
rk(i)=(db/2)/cos(alfa(i))theta(i)=tan(alfa(i))-alfa(i)!
弧度表示
*enddo
*do,i,1,10
fai(i)=(theta_s+fai_s-theta(i))*180/pi!
角度
*enddo
CSYS,1
*do,i,1,10
k,i,rk(i),fai(i)*enddo
ksel,all
BSPLIN,ALL!
绘制齿廓线
csys,0
LSYMM,Y,1
K,30000
LARC,10,12,30000,da/2csys,1
LSEL,ALL
LGEN,Z1,ALL,,,0,-deata_ang
kpoint=10+1
*do,i,1,z1-1,1
larc,kpoint,kpoint+2,30000,df/2
kpoint=kpoint+4
*enddo
larc,187,1,30000,df/2LSEL,ALL
AL,ALL
CYL4,,,50
ASBA,1,2csys,0
AGEN,2,3,,,1.00*d,,,,0
!
*********************************************************
ET,1,PLANE82
MP,EX,1,500000
MP,NUXY,1,0.3
MP,EX,2,400000
MP,NUXY,2,0.3
TYPE,1
MAT,1REAL,
ESYS,0
SECNUM,
!
*
SMRT,3
MSHAPE,0,2DMSHKEY,0
AMESH,3
!
*
MAT,2
AMESH,1
FLST,5,90,2,ORDE,73
FITEM,5,5FITEM,5,-6FITEM,5,61FITEM,5,-62FITEM,5,83FITEM,5,-84FITEM,5,89FITEM,5,-90FITEM,5,145FITEM,5,-147FITEM,5,180FITEM,5,-181FITEM,5,187FITEM,5,211FITEM,5,-212FITEM,5,216FITEM,5,-217FITEM,5,242FITEM,5,-243FITEM,5,248FITEM,5,-249FITEM,5,393FITEM,5,-394FITEM,5,444FITEM,5,-446FITEM,5,461FITEM,5,-464FITEM,5,483FITEM,5,-484
FITEM,5,525FITEM,5,584FITEM,5,587FITEM,5,624FITEM,5,635FITEM,5,641FITEM,5,643FITEM,5,681FITEM,5,703FITEM,5,-704FITEM,5,732FITEM,5,737FITEM,5,774FITEM,5,-775FITEM,5,864FITEM,5,-865FITEM,5,882FITEM,5,-883FITEM,5,895FITEM,5,903FITEM,5,-904FITEM,5,922FITEM,5,941FITEM,5,961FITEM,5,967FITEM,5,-968FITEM,5,982FITEM,5,1000FITEM,5,-1003
FITEM,5,1114FITEM,5,-1117
FITEM,5,1127FITEM,5,1163FITEM,5,-1166
FITEM,5,1171FITEM,5,-1176
FITEM,5,1266FITEM,5,1270FITEM,5,1310FITEM,5,-1312
FITEM,5,1343FITEM,5,-1346
FITEM,5,1350FITEM,5,-1351
CM,_Y,ELEM
ESEL,,,,P51XCM,_Y1,ELEM
CMSEL,S,_Y
CMDELE,_Y
!
*
!
*
EREF,_Y1,,,1,0,1,1
CMDELE,_Y1
CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_KPCM,KP
CM,_LINECM,LINECM,_AREACM,AREACM,_VOLUCM,VOLU/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2R,3,
RMORE,
RMORE,,0
RMORE,0
!
Generatethetargetsurface
LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,3
LSEL,A,,,4
LSEL,A,,,5
LSEL,A,,,6
LSEL,A,,,7
LSEL,A,,,8
LSEL,A,,,9
LSEL,A,,,10
LSEL,A,,,128
LSEL,A,,,130
LSEL,A,,,131
LSEL,A,,,132LSEL,A,,,133LSEL,A,,,134LSEL,A,,,135LSEL,A,,,136LSEL,A,,,137LSEL,A,,,138LSEL,A,,,178LSEL,A,,,179LSEL,A,,,180CM,_TARGET,LINE
TYPE,2
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
!
Generatethecontactsurface
LSEL,S,,,264LSEL,A,,,265LSEL,A,,,266LSEL,A,,,267LSEL,A,,,268LSEL,A,,,269LSEL,A,,,270LSEL,A,,,271LSEL,A,,,272LSEL,A,,,273LSEL,A,,,274LSEL,A,,,275LSEL,A,,,276LSEL,A,,,277LSEL,A,,,278LSEL,A,,,279LSEL,A,,,280LSEL,A,,,281LSEL,A,,,282LSEL,A,,,283LSEL,A,,,284LSEL,A,,,285LSEL,A,,,286LSEL,A,,,287LSEL,A,,,288CM,_CONTACT,LINE
TYPE,3
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECMCMDEL,_NODECMCMSEL,A,_ELEMCMCMDEL,_ELEMCMCMSEL,S,_KPCMCMDEL,_KPCM
CMSEL,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACMCMSEL,S,_VOLUCMCMDEL,_VOLUCM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACTNROPT,UNSYM
LOCAL,11,1,180,0,0,,,,1,1,
EPLOT
CSYS,1
LSEL,ALL
FLST,5,4,4,ORDE,2
FITEM,5,181
FITEM,5,-184
LSEL,S,,,P51XNSEL,ALL
NSLL,S,1
NPLOT
FLST,2,88,1,ORDE,4FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
NROTAT,P51X
FINISH
/SOL
FLST,2,88,1,ORDE,4FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
!
*
/GO
D,P51X,,,,,,UX,,,,,
FLST,2,88,1,ORDE,4FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
!
*
/GO
D,P51X,,0.6,,,,UY,,,,,
ALLSEL,ALL
CSYS,11,
EPLOT
LSEL,ALL
FLST,5,4,4,ORDE,2FITEM,5,365
FITEM,5,-368
LSEL,S,,,P51XNSEL,ALL
NSLL,S,1
NPLOT
FINISH
/PREP7
FLST,2,88,1,ORDE,4FITEM,2,2815
FITEM,2,-2818FITEM,2,3269
FITEM,2,-3352NROTAT,P51X
FINISH
/SOL
FLST,2,88,1,ORDE,4FITEM,2,2815
FITEM,2,-2818FITEM,2,3269
FITEM,2,-3352!
*
/GO
D,P51X,,0,,,,UX,,,,,
ALLSEL,ALL
EPLOT
!
*
!
*
OUTRES,ALL,ALL,!
*
TIME,