请教齿轮接触的问题
我如果想要得到齿轮啮合的最大接触应力,作静力分析的话,是不是要进行一系列的从啮入到退出的分析?还是有其它较为方便的办法,可以得到整个啮合过程的应力变化?还请不吝赐教!
这是个问题!
是啊,齿轮分析很麻烦!斑竹\各位高手帮帮小学生吧!我也有这方面的问题!
比如,能不能让齿轮转动起来分析接触啊?这样更好!
斑竹大人,能就这个问题仔细说说吗?需要解决这类问题的人很多?
拜托......... http://www.ugcn.cn/common/download.asp?module=upload&id=79785 一个例子:
/GST,ON
!
*DIM,ANG,,53 !定义53个计算点
ANG(1)=-5.625,-5.5,-5.25,-5,-4.75,-4.5,-4.25
ANG=-4,-3.75,-3.5,-3.25,-3,-2.75,-2.625,-2.5,-2.375
ANG(17)=-2.25,-2.125,-2,-1.75,-1.5
ANG(22)=-1.25,-1,-0.75,-0.5,-0.25
ANG(27)=0.00001,0.25,0.5,0.75,1
ANG(32)=1.25,1.5,1.75,2,2.125
ANG(37)=2.25,2.375,2.5,2.625,2.75
ANG(42)=3,3.25,3.5,3.75,4
ANG(47)=4.25,4.5,4.75,5,5.25
ANG(52)=5.5,5.625
FNAM=1
PARSAV,ALL,ANGLE_A,,
OLOOP
RESUME,Cornell_5,db,,0,0
!Cornell_5是个2D齿轮模型接触点在Pitch point
PARRES,,ANGLE_A,, ,,
/PREP7
!*
TOQ=150 !Nm
!*
!*
CSYS,4
FLST,3,52,5,ORDE,2
FITEM,3,53
FITEM,3,-104
AGEN, ,P51X, , , ,ANG(%FNAM%), , , ,1 !转被动轮
CSYS,1
FLST,3,52,5,ORDE,2
FITEM,3,1
FITEM,3,-52
AGEN, ,P51X, , , ,-ANG(%FNAM%), , , ,1 !转主动轮
SMRT,1 !定义接触区内smart mesh
ESIZE,0.2,0, !接触区单元大小0.2
MSHKEY,0
!*
FLST,5,6,5,ORDE,6
FITEM,5,2
FITEM,5,4
FITEM,5,6
FITEM,5,54
FITEM,5,56
FITEM,5,58
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
/UI,MESH,OFF
!*
!*
/COM, CONTACT PAIR CREATION - START !生成接触单元
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_LINECM,LINE
CM,_AREACM,AREA
/GSAV,cwz,gsav,,temp
MP,MU,1,0.1
MAT,1
R,3
REAL,3
ET,2,169
ET,3,172
KEYOPT,3,9,0
R,3,
RMORE,
RMORE,,0
! Generate the target surface
LSEL,S,,,188
LSEL,A,,,195
LSEL,A,,,202
CM,_TARGET,LINE
TYPE,2
NSLL,S,1
ESLN,S,0
ESURF,ALL
CMSEL,S,_ELEMCM
! Generate the contact surface
LSEL,S,,,110
LSEL,A,,,114
LSEL,A,,,128
CM,_CONTACT,LINE
TYPE,3
NSLL,S,1
ESLN,S,0
ESURF,ALL
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,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
!*
CSYS,1 !激活元拄坐标
FLST,2,74,1,ORDE,18
FITEM,2,2572
FITEM,2,2588
FITEM,2,-2591
FITEM,2,2636
FITEM,2,-2639
FITEM,2,2684
FITEM,2,2700
FITEM,2,-2702
FITEM,2,2732
FITEM,2,-2733
FITEM,2,2849
FITEM,2,2857
FITEM,2,3086
FITEM,2,-3087
FITEM,2,3105
FITEM,2,-3106
FITEM,2,3324
FITEM,2,-3376
NROTAT,P51X !转主动轮HUB的节点坐标为元拄坐标
!*
CSYS,4 !激活WORKING坐标(working plane is pola, it already been set at ! the centre of drven gear)
FLST,2,74,1,ORDE,17
FITEM,2,6350
FITEM,2,6366
FITEM,2,-6369
FITEM,2,6414
FITEM,2,-6417
FITEM,2,6462
FITEM,2,6478
FITEM,2,-6480
FITEM,2,6510
FITEM,2,6518
FITEM,2,6855
FITEM,2,-6856
FITEM,2,6960
FITEM,2,-6961
FITEM,2,6979
FITEM,2,-7032
FITEM,2,7384
NROTAT,P51X !转被动轮HUB的节点坐标为WORKING坐标
FLST,2,74,1,ORDE,17
FITEM,2,6350
FITEM,2,6366
FITEM,2,-6369
FITEM,2,6414
FITEM,2,-6417
FITEM,2,6462
FITEM,2,6478
FITEM,2,-6480
FITEM,2,6510
FITEM,2,6518
FITEM,2,6855
FITEM,2,-6856
FITEM,2,6960
FITEM,2,-6961
FITEM,2,6979
FITEM,2,-7032
FITEM,2,7384
!*
/GO
D,P51X, , , , , ,ALL, , , , , !加约束
FLST,2,74,1,ORDE,18
FITEM,2,2572
FITEM,2,2588
FITEM,2,-2591
FITEM,2,2636
FITEM,2,-2639
FITEM,2,2684
FITEM,2,2700
FITEM,2,-2702
FITEM,2,2732
FITEM,2,-2733
FITEM,2,2849
FITEM,2,2857
FITEM,2,3086
FITEM,2,-3087
FITEM,2,3105
FITEM,2,-3106
FITEM,2,3324
FITEM,2,-3376
!*
/GO
D,P51X, , , , , ,UX, , , , , !加约束
FLST,2,74,1,ORDE,18
FITEM,2,2572
FITEM,2,2588
FITEM,2,-2591
FITEM,2,2636
FITEM,2,-2639
FITEM,2,2684
FITEM,2,2700
FITEM,2,-2702
FITEM,2,2732
FITEM,2,-2733
FITEM,2,2849
FITEM,2,2857
FITEM,2,3086
FITEM,2,-3087
FITEM,2,3105
FITEM,2,-3106
FITEM,2,3324
FITEM,2,-3376
!*
/GO
F,P51X,FY,-TOQ*1.0640561821664183868908278357097 !加载
FLST,4,74,1,ORDE,18
FITEM,4,2572
FITEM,4,2588
FITEM,4,-2591
FITEM,4,2636
FITEM,4,-2639
FITEM,4,2684
FITEM,4,2700
FITEM,4,-2702
FITEM,4,2732
FITEM,4,-2733
FITEM,4,2849
FITEM,4,2857
FITEM,4,3086
FITEM,4,-3087
FITEM,4,3105
FITEM,4,-3106
FITEM,4,3324
FITEM,4,-3376
CP,1,UY,P51X !主动轮HUB加转向的绑定方程
FINISH
/SOLU
!*
NLGEOM,1
NROPT,AUTO, ,
LUMPM,0
EQSLV,ICCG,1.0e-5,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF
PSTRES
TOFFST,0,
!*
!*
TIME,1
AUTOTS,1
NSUBST,4,8,2,0
KBC,0
!*
TSRES,ERASE
AUTOTS,-1.0
NEQIT,200
/STATUS,SOLU
SOLVE
!*
!*
/POST1
RSYS,1
PRNSOL,DOF,U,Y
*CFOPEN,Cornell_5-%TOQ%Nm,lis,L:\Results2\,APPEND !打开*lis文件, L硬盘名
*CFWRITE,,%ANG(FNAM)%,UY%FNAM%,UY(2572) !写入计算结果UY
*CFCLOSE
SAVE,Cornell_5-%TOQ%Nm%ANG(FNAM)%,DB,L:\Results2\,0,0 !存盘以便进一步提取其他数据
FINISH
FNAM=FNAM+1
*IF,FNAM,GT,53,:ENDLOOP
PARSAV,ALL,ANGLE_A,,
*GO,OLOOP
:ENDLOOP
FINISH
!/CLEAR,START !如有必要在这里开始新一轮计算
!/INPUT,Cornell_5_250Nm,log,,, 0 楼上的朋友,能说明一下,你这个APDL程序,是关于什么的啊?
齿轮啮合接触?
希望你能和我联系.QQ110944020
邮箱:lwt5841314520@163.com 算不了啊
第五行"ANG=-4,-3.75,-3.5,-3.25,-3,-2.75,-2.625,-2.5,-2.375"ANG(?)括号里的数字是多少?
Cornell_5是个2D齿轮模型,模型在那?
[ 本帖最后由 syzl 于 2007-7-9 09:42 编辑 ]
没有这么复杂吧
我用ANSYS11做了一个简单的例子,只要计算机容量够大,够耐心,接触的过程和结果都能做得出来. 这是动画,是否是这样楼上,如何让齿轮转起来的啊?
这个太受用了~
多谢分享了 ~ 回复 8 # fengfzg 的帖子请教一下,这个是如何实现的?由于要做论文很需要这方面的知识,我做的是少齿差传动,希望能指导一下,多谢! 回复 8 # fengfzg 的帖子
怎么让它动起来,还能显示应力变化呢?请指教! 回复 8 # fengfzg 的帖子
能否详细介绍下怎么实现的? 谢谢 回复 8 # fengfzg 的帖子
我是求啮合刚度啊,你这个是怎么转起来的?
页:
[1]