齿轮接触动力学分析

lt511888

有没有人有齿轮接触动力学分析的资料呀，小弟急需。。。谢谢

warshipcj

/FILNAME,Spur Gear Based on ANSYS,0
/PREP7
!定义小齿轮齿轮参数
M=2E-3                  !齿轮模数
PI=3.1415               !定义π
Z1=36                   !小齿轮齿数
ALFA1=20/180*PI         !定义小齿轮的标准压力角
Ha=1                    !定义齿顶高系数
C=0.25                  !定义顶隙系数
D1=M*Z1                 !定义小齿轮分度圆直径
Dk1=M*Z1*COS(ALFA1)     !定义小齿轮基圆直径
Rk1=Dk1/2               !定义小齿轮基圆半径
Df1=D1-2*(Ha+C)*M       !定义小齿轮齿根圆直径
Rf1=Df1/2               !定义小齿轮齿根圆半径
Da1=D1+2*Ha*M           !定义小齿轮齿顶圆直径
Ra1=Da1/2               !定义小齿轮齿顶圆半径
ALFA_Ding1=ACOS(Rk1/Ra1)!定义小齿轮齿顶压力角   
SITA_Ding1=(TAN(ALFA_Ding1)-ALFA_Ding1)*180/PI !定义小齿轮齿顶展角
SITA1=(tan(ALFA1)-ALFA1)*180/PI                !定义小齿轮分度圆展角
!定义大齿轮几何参数
Z2=45                   !大齿轮齿数
ALFA2=20/180*PI         !定义大齿轮的标准压力角
D2=M*Z2                 !定义大齿轮分度圆直径
Dk2=M*Z2*COS(ALFA2)     !定义大齿轮基圆直径
Rk2=Dk2/2               !定义大齿轮基圆半径
Df2=D2-2*(Ha+C)*M       !定义大齿轮齿根圆直径
Rf2=Df2/2               !定义大齿轮齿根圆半径
Da2=D2+2*Ha*M           !定义大齿轮齿顶圆直径
Ra2=Da2/2               !定义大齿轮齿顶圆半径
ALFA_Ding2=ACOS(Rk2/Ra2)!定义大齿轮齿顶压力角   
SITA_Ding2=(TAN(ALFA_Ding2)-ALFA_Ding2)*180/PI !定义大齿轮齿顶展角
SITA2=(tan(ALFA2)-ALFA2)*180/PI                !定义大齿轮分度圆展角

ET,1,SOLID185
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06E11
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7850


!定义小齿轮的渐开线
ALFA_Dingn1=alfa_ding1/10 !定义压力角份数
*DIM,ALFA,ARRAY,10,1      !定义10×1数组ALFA
*DIM,SITA,ARRAY,9,1       !定义9×1数组SITA
*DIM,R,ARRAY,9,1          !定义9×1数组R
CSYS,1                    !当前坐标系为柱坐标
ALFA(1,1)=0               !数组数ALFA(1,1)=0
K,100,Rk1,0               !生成关键点100
*DO,i,1,9,1               !循环生成其他关键点
ALFA(i+1,1)=ALFA(i,1)+ALFA_Dingn1
SITA(i,1)=(tan(alfa(i+1,1))-alfa(i+1,1))*180/PI
R(i,1)=Rk1/cos(alfa(i+1,1))
/prep7
K,i+100,R(i,1),sita(i,1)
*enddo                          !循环结束生成9个关键点
k,110,ra1,sita_ding1            !生成齿顶圆上的关键点
bsplin,all                      !用样条曲线连接所有点
!定义小齿轮齿根过渡曲线
!定义小齿轮齿根圆
FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,Rf1,0,0  
CIRCLE,P51X, , , ,360, ,
!定义线
L,100,1
!定义齿根过渡曲线  
LCOMB,1,6

LFILLT,1,5,0.38*m, ,
LDELE,2, , ,1
LDELE,3, , ,1
LDELE,4, , ,1
LDELE,5, , ,1
LCOMB,1,6
!定义单个齿的平面模型
LOCAL,11,0,0,0,0,SITA1+90/Z1, , ,1,1,  
CSYS,11,
LSYMM,Y,1,  

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,RF1,0,0  
CIRCLE,P51X

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,RA1,0,0  
CIRCLE,P51X

CSYS,1

FLST,3,2,4,ORDE,2   
FITEM,3,1   
FITEM,3,-2  
LGEN,2,P51X, , , ,360/Z1, , ,0  

LOVLAP,ALL  !搭接所有的线

!删除指定的线
LDELE,21, , ,1
LDELE,8, , ,1
LDELE,9, , ,1
LDELE,10, , ,1
LDELE,18, , ,1
LDELE,4, , ,1
LDELE,5, , ,1
LDELE,19, , ,1
LDELE,12, , ,1
LDELE,22, , ,1
LDELE,23, , ,1
LDELE,11, , ,1
LDELE,13, , ,1
LDELE,20, , ,1
LDELE,15, , ,1
LDELE,17, , ,1
!通过复制产生完整的齿轮轮廓
LGEN,Z1,ALL, , , ,360/Z1, , ,0
LGLUE,ALL !粘结所有的线
!通过齿轮轮廓线建立齿轮平面模型
FLST,2,144,4
FITEM,2,160
FITEM,2,60  
FITEM,2,58  
FITEM,2,59  
FITEM,2,159
FITEM,2,56  
FITEM,2,54  
FITEM,2,55  
FITEM,2,158
FITEM,2,52  
FITEM,2,50  
FITEM,2,51  
FITEM,2,157
FITEM,2,48  
FITEM,2,46  
FITEM,2,47  
FITEM,2,156
FITEM,2,44  
FITEM,2,42  
FITEM,2,43  
FITEM,2,155
FITEM,2,40  
FITEM,2,38  
FITEM,2,39  
FITEM,2,154
FITEM,2,36  
FITEM,2,34  
FITEM,2,35  
FITEM,2,153
FITEM,2,32  
FITEM,2,30  
FITEM,2,31  
FITEM,2,152
FITEM,2,28  
FITEM,2,26  
FITEM,2,27  
FITEM,2,151
FITEM,2,24  
FITEM,2,22  
FITEM,2,23  
FITEM,2,150
FITEM,2,20  
FITEM,2,18  
FITEM,2,19  
FITEM,2,149
FITEM,2,15  
FITEM,2,12  
FITEM,2,13  
FITEM,2,148
FITEM,2,10  
FITEM,2,8   
FITEM,2,9   
FITEM,2,147
FITEM,2,6   
FITEM,2,4   
FITEM,2,5   
FITEM,2,3   
FITEM,2,145
FITEM,2,2   
FITEM,2,14  
FITEM,2,1   
FITEM,2,146
FITEM,2,142
FITEM,2,143
FITEM,2,180
FITEM,2,140
FITEM,2,138
FITEM,2,139
FITEM,2,179
FITEM,2,136
FITEM,2,134
FITEM,2,135
FITEM,2,178
FITEM,2,132
FITEM,2,130
FITEM,2,131
FITEM,2,177
FITEM,2,128
FITEM,2,126
FITEM,2,127
FITEM,2,176
FITEM,2,124
FITEM,2,122
FITEM,2,123
FITEM,2,175
FITEM,2,120
FITEM,2,118
FITEM,2,119
FITEM,2,174
FITEM,2,116
FITEM,2,114
FITEM,2,115
FITEM,2,173
FITEM,2,112
FITEM,2,110
FITEM,2,111
FITEM,2,172
FITEM,2,108
FITEM,2,106
FITEM,2,107
FITEM,2,171
FITEM,2,104
FITEM,2,102
FITEM,2,103
FITEM,2,170
FITEM,2,100
FITEM,2,98  
FITEM,2,99  
FITEM,2,169
FITEM,2,96  
FITEM,2,94  
FITEM,2,95  
FITEM,2,168
FITEM,2,92  
FITEM,2,90  
FITEM,2,91  
FITEM,2,167
FITEM,2,88  
FITEM,2,86  
FITEM,2,87  
FITEM,2,166
FITEM,2,84  
FITEM,2,82  
FITEM,2,83  
FITEM,2,165
FITEM,2,80  
FITEM,2,78  
FITEM,2,79  
FITEM,2,164
FITEM,2,76  
FITEM,2,74  
FITEM,2,75  
FITEM,2,163
FITEM,2,72  
FITEM,2,70  
FITEM,2,71  
FITEM,2,162
FITEM,2,68  
FITEM,2,66  
FITEM,2,67  
FITEM,2,161
FITEM,2,64  
FITEM,2,62  
FITEM,2,63  
AL,P51X

PCIRC,0.4*RA1, ,0,360,  
ASBA,       1,       2  
VOFFST,3,1E-2, ,
!建立大齿轮模型
!定义大齿轮的渐开线
ALFA_Dingn2=alfa_ding2/10    !定义压力角份数
csys,1                     !当前坐标系为柱坐标
ALFA(1,1)=0                !数组数alfa_12(1,1)=0
K,1000,Rk2,0                !生成关键点100
*DO,i,1,9,1                !循环生成其他关键点
ALFA(i+1,1)=ALFA(i,1)+ALFA_Dingn2
SITA(i,1)=(tan(ALFA(i+1,1))-ALFA(i+1,1))*180/PI
R(i,1)=Rk2/cos(alfa(i+1,1))
/prep7
K,i+1000,R(i,1),sita(i,1)
*enddo                          !循环结束生成9个关键点
K,1010,Ra2,sita_ding2            !生成齿顶圆上的关键点
KSEL,S, , ,1000,1010
bsplin,all                      !用样条曲线连接所有点
ALLSEL,ALL
!定义小齿轮齿根过渡曲线
FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,Rf2,0,0  
CIRCLE,P51X, , , ,360, ,


LOVLAP,445,446
LCOMB,449,450
LFILLT,453,449,0.38*m, ,
LCOMB,445,453
LDELE,452, , ,1
LDELE,451, , ,1
LDELE,447, , ,1
LDELE,448, , ,1
LDELE,449, , ,1

!定义单个齿的平面模型
LOCAL,11,0,0,0,0,SITA2+90/Z2, , ,1,1,  
CSYS,11,
LSYMM,Y,445,

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,RF2,0,0  
CIRCLE,P51X

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,RA2,0,0  
CIRCLE,P51X

CSYS,1

FLST,3,2,4,ORDE,2   
FITEM,3,445
FITEM,3,-446
LGEN,2,P51X, , , ,360/Z2, , ,0  

LOVLAP,447,450,451,454,456,455,446,445
LDELE,456, , ,1
LDELE,452, , ,1
LDELE,453, , ,1
LDELE,452, , ,1
LDELE,465, , ,1
LDELE,462, , ,1
LDELE,448, , ,1
LDELE,449, , ,1
LDELE,463, , ,1
LDELE,467, , ,1
LDELE,466, , ,1
LDELE,457, , ,1
LDELE,461, , ,1
LDELE,459, , ,1
LDELE,455, , ,1
LDELE,464, , ,1
LDELE,454, , ,1

FLST,2,4,4,ORDE,4   
FITEM,2,445
FITEM,2,-446
FITEM,2,458
FITEM,2,460
LCOMB,P51X, ,0

LGEN,Z2,445, , , ,360/Z2, , ,0

LSEL,S, , ,445,489
LGLUE,ALL
ALLSEL,ALL
LSEL,S,LOC,X,Ra1*1.01,Ra2*1.01
NUMMRG,ALL, , , ,LOW

FLST,2,135,4
FITEM,2,519
FITEM,2,518
FITEM,2,611
FITEM,2,515
FITEM,2,514
FITEM,2,613
FITEM,2,511
FITEM,2,510
FITEM,2,615
FITEM,2,509
FITEM,2,508
FITEM,2,617
FITEM,2,505
FITEM,2,504
FITEM,2,619
FITEM,2,553
FITEM,2,552
FITEM,2,621
FITEM,2,521
FITEM,2,520
FITEM,2,623
FITEM,2,507
FITEM,2,506
FITEM,2,624
FITEM,2,503
FITEM,2,502
FITEM,2,622
FITEM,2,517
FITEM,2,516
FITEM,2,620
FITEM,2,513
FITEM,2,512
FITEM,2,618
FITEM,2,531
FITEM,2,530
FITEM,2,616
FITEM,2,549
FITEM,2,548
FITEM,2,614
FITEM,2,497
FITEM,2,496
FITEM,2,612
FITEM,2,545
FITEM,2,544
FITEM,2,610
FITEM,2,501
FITEM,2,500
FITEM,2,608
FITEM,2,493
FITEM,2,492
FITEM,2,606
FITEM,2,491
FITEM,2,490
FITEM,2,604
FITEM,2,499
FITEM,2,498
FITEM,2,602
FITEM,2,494
FITEM,2,495
FITEM,2,600
FITEM,2,579
FITEM,2,578
FITEM,2,598
FITEM,2,577
FITEM,2,576
FITEM,2,596
FITEM,2,575
FITEM,2,574
FITEM,2,594
FITEM,2,573
FITEM,2,572
FITEM,2,592
FITEM,2,571
FITEM,2,570
FITEM,2,590
FITEM,2,569
FITEM,2,568
FITEM,2,588
FITEM,2,567
FITEM,2,566
FITEM,2,586
FITEM,2,565
FITEM,2,564
FITEM,2,584
FITEM,2,563
FITEM,2,562
FITEM,2,582
FITEM,2,561
FITEM,2,560
FITEM,2,580
FITEM,2,559
FITEM,2,558
FITEM,2,581
FITEM,2,557
FITEM,2,556
FITEM,2,583
FITEM,2,555
FITEM,2,554
FITEM,2,585
FITEM,2,551
FITEM,2,550
FITEM,2,587
FITEM,2,547
FITEM,2,546
FITEM,2,589
FITEM,2,543
FITEM,2,542
FITEM,2,591
FITEM,2,541
FITEM,2,540
FITEM,2,593
FITEM,2,537
FITEM,2,535
FITEM,2,595
FITEM,2,536
FITEM,2,534
FITEM,2,597
FITEM,2,539
FITEM,2,538
FITEM,2,599
FITEM,2,533
FITEM,2,532
FITEM,2,601
FITEM,2,529
FITEM,2,528
FITEM,2,603
FITEM,2,527
FITEM,2,526
FITEM,2,605
FITEM,2,525
FITEM,2,524
FITEM,2,607
FITEM,2,523
FITEM,2,522
FITEM,2,609
AL,P51X
ALLSEL,ALL
!定义一个圆形
PCIRC,0.4*RA2, ,0,360,  
ASBA,       151,       152  !进行面的布尔减法
VOFFST,153,1E-2, ,  !把齿轮面沿着法向拉伸成体

CSYS,0
VGEN, ,2, , ,81e-3, , , , ,1
CSYS,1
VGEN, ,1, , , ,3.08, , , ,1

SMRT,8        !定义智能网格划分，等级为8
ESIZE,0.004,0, !定义单元的总体尺寸
VSWEEP,ALL    !通过体扫描进行网格划分



!定义接触节点组
ASEL,R, , ,50,100
NSLA,S,1
CM,CE,NODE
  
ALLSEL,ALL
!定义目标节点组  
ASEL,R, , ,154,200
NSLA,S,1
CM,TE,NODE
ALLSEL,ALL


!使用接触向导建立接触对
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.28
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3  
ET,2,170
ET,3,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
NROPT,UNSYM
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2   
KEYOPT,3,11,0   
KEYOPT,3,12,0   
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface   
NSEL,S,,,TE
CM,_TARGET,NODE
TYPE,2  
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,S,,,CE
CM,_CONTACT,NODE
TYPE,3  
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
ALLSEL,ALL



CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
! Define surface-based constraint type of pair  
MAT,1   
R,4
REAL,4  
ET,4,170
ET,5,174
KEYOPT,5,12,5   
KEYOPT,5,4,2
KEYOPT,5,2,2
KEYOPT,4,2,0
KEYOPT,4,4,11111
TYPE,4  
! Create a pilot node   
! At center of mass of contact geometric entities   
CSYS,0  
N,31196, 4.32192713462e-019,-5.87422651478e-018,0.005   
CSYS,1  
TSHAP,PILO  
E,31196
! Generate the contact surface  
ASEL,S,,,147
ASEL,A,,,148
ASEL,A,,,149
ASEL,A,,,150
CM,_CONTACT,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
!*  
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
! Define surface-based constraint type of pair  
MAT,1   
R,5
REAL,5  
ET,6,170
ET,7,174
KEYOPT,7,12,5   
KEYOPT,7,4,1
KEYOPT,7,2,2
KEYOPT,6,2,0
KEYOPT,6,4,11111
TYPE,6  
! Create a pilot node   
! At center of mass of contact geometric entities   
CSYS,0  
N,31197, 0.081,-5.78874346437e-018,0.005
CSYS,1  
TSHAP,PILO  
E,31197
! Generate the contact surface  
ASEL,S,,,288
ASEL,A,,,289
ASEL,A,,,290
ASEL,A,,,291
CM,_CONTACT,AREA
TYPE,7  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,5  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,5  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  



/SOL
ANTYPE,4   !设置分析类型为瞬态动力学分析
TRNOPT,FULL !使用完全法
LUMPM,0
!定义函数
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'NJ'  
*SET,_FNCCSYS,0
! /INPUT,NJ.func,,,1
*DIM,%_FNCNAME%,TABLE,6,3,1,,,,%_FNCCSYS%   
!   
! Begin of equation: 1200   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1200, 0, 0, 0   
*SET,%_FNCNAME%(0,2,1),   0
*SET,%_FNCNAME%(0,3,1),   0
! End of equation: 1200
!-->
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'RZ'  
*SET,_FNCCSYS,0
! /INPUT,RZ.func,,,1
*DIM,%_FNCNAME%,TABLE,6,3,1,,,,%_FNCCSYS%   
!   
! Begin of equation: 0.2*{TIME}
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.2, 0, 0, 1
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,3,1),   0, 99, 0, 1, -2, 0, 0
! End of equation: 0.2*{TIME}   

F,31197,MZ, %NJ% !施加阻力矩
D,31196, , %RZ% , , , ,ROTZ, , , , ,  !施加驱动转动位移
OUTRES,ALL,ALL,
TIME,1
NSUBST,10
SOLVE

circle_out

上面贴的命令流运行后结果，好像是0.4E10 Pa吧，这样的话就是4000MPa，这明显大于屈服极限啦啊，这应该不正常的吧，是什么原因呢？