浅谈齿轮强度设计几个问题的探讨论文
0 引言
齿轮传动是机械传动中最重要的传动之一。公元前300 多年，古希腊哲学家亚里士多德在《机械问题》中，就阐述了用青铜或铸铁齿轮传递旋转运动的问题。17 世纪末到18 世纪初，人们开始对齿轮的强度问题进行研究。欧洲工业革命以后，齿轮技术得到高速发展，齿轮传动在机械传动及整个机械领域中的应用极其广泛。齿轮设计成为机械设计中重要的设计内容之一。目前国际上比较常见的有关齿轮强度设计公式，除了我国的国家标准( GB) 有关齿轮强度的计算方法以外主要有: 国际标准化组织( ISO) 计算方法; 美国齿轮制造商协会( AGMA) 标准计算方法;德国工业标准( DIN) 计算方法; 日本齿轮工业会( JGMA)计算方法; 英国BS 计算方法等。作者在从事机械设计特别对齿轮设计的教学中，发现不少地方的知识点描述比较简单，不容易理解，为此，在文中对齿轮设计的几个问题如齿轮的失效方式、齿轮强度设计的历史、现状进行了深入分析，探讨我国齿轮强度设计的历史来源以及在齿轮设计中的一些困惑。通过深入的分析，有助于大家更好地理解齿轮设计公式的意义和来龙去脉。
1 齿轮失效方式的探讨
齿轮在传动过程中会出现各种形式的失效，甚至丧失传动能力。齿轮传动的失效方式与齿轮的材料、热处理方式、润滑条件、载荷大小、载荷变化规律以及转动速度等有关。人们对齿轮失效的认识是一个发展的过程。18 世纪中叶人们就开始对齿轮的失效进行研究。对齿轮摩擦磨损、点蚀形成和齿面胶合有了初步的认识。1928 年，白金汉发表了有关齿轮磨损的论文，并将齿面失效分为点蚀、磨粒磨损、胶合、剥落、擦伤和咬死等6 种失效形式。1939 年，Rideout 将齿轮损伤分为正常磨损、点蚀、剥落、胶合、擦伤、切伤、滚轧和锤击等8 种形式。1953 年Borsoff 和Sorem 将齿轮损伤分为6 类。1967 年尼曼根据大量试验，对渐开线齿轮的4 种失效形式画出了承载能力的限制关系图，并指出当齿轮转速较低时，影响软齿面齿轮承载能力的主要因素是点蚀，影响硬齿面齿轮承载能力的是断齿; 而对于高速重载传动齿轮，影响因素往往是胶合。自上世纪50 年代以来，一些国家以标准的形式对齿轮损伤形式进行分类，对名词术语、表现特征、引发原因等都有规定。如1951 年美国将齿轮损伤分为两大类，一类是齿面损坏，包括磨损、塑性变形、胶合、表面疲劳等，另一类是轮齿的折断。前一大类齿面损坏是齿轮作为高副由于摩擦学原因而引起的表面损伤; 后一大类轮齿的折断是轮齿作为受力构件由于体积强度不够而发生的破坏。1968 年奥地利国家标准规定了齿轮损伤的名词术语。
1983 年，我国颁布了齿轮轮齿损伤的术语、特征和原因国家标准( GB /T3481 - 83) ，将齿轮损伤形式分为5 大类，即磨损、齿面疲劳( 包括点蚀和剥落) 、塑性变形、轮齿折断和其他损伤，共26 种失效形式。1997 年，我国颁布了对GB/T3481 - 1983 修订的GB/T3481 -1997 国家标准。目前我国在大多数的机械设计教材和机械设计手册中齿轮失效方式都进行了简化，一般分为5 大类，即轮齿折断、齿面疲劳点蚀、齿面胶合、齿面磨损和塑性变形。
2 齿轮强度设计的探讨
2. 1 轮齿弯曲强度计算
1785 年，英国瓦特提出了齿根弯曲强度的计算方法，把轮齿看成为矩形截面的板状悬臂梁，随后出现多种弯曲强度计算公式。1893年，路易斯发表了轮齿弯曲强度计算式，而且用内切抛物线法找齿轮的危险截面，这一方法称为“抛物线法”[12]，如图1 所示。路易斯以载荷作用于齿顶推导出齿根弯曲应力公式，但是对于重合度大于1 小于2 的齿轮传动，理论上只有当单对齿啮合时，载荷才全部由一个齿承受。对于重合度大于2 小于3 的足够精密的齿轮，因为同时有2 对以上的齿轮在啮合，其最大弯曲应力的作用点要低。
在此之后，又出现30°切线法、尼曼法、白金汉法等。1980 年， ISO 提出“渐开线圆柱齿轮承载能力的基本原理”( ISO 6336 - 1980) ，公布了轮齿弯曲强度、齿面接触强度的计算方法。
过去，我国的齿轮强度计算方法一直比较混乱，没有统一的标准，对生产、科研以及教学带来诸多问题。于是， 1981 年我国成立了“渐开线圆柱齿轮承载能力计算方法”国家标准课题组，以ISO6336—1980为根据，开展全面的研究工作。1983 年颁布了渐开线圆柱齿轮承载能力计算方法的国家标准( GB /T3480—1983) 。
目前，我国有关齿轮弯曲强度的设计公式基本上采用30° 切线法，即作与轮齿对称中心线成30°夹角并与齿根圆角相切的斜线，两切点的连线是齿根危险截面位置。而且以单对齿啮合区的最高点作为最不利载荷作用点，这时产生的弯曲应力最大，如图2 所示。另外，弯曲疲劳强度计算公式中，齿形系数在许多机械设计中只是说明与齿数有关，与模数无关，并未做详细说明，不容易理解。下面对相关问题进行详细分析。如图2 所示，齿根弯曲应力为σF =MW= FnhFcosαFbS2F /6 = 6KFthFcosαFbS2Fcosα= KFtbm6( hFm) cosαF( SFm)2cosα( 1)式中，αF为齿顶圆压力角。令式( 1) 中的YF =6( hFm) cos αF( SFm)2cos α式中，YF称为齿形系数，由路易斯在其轮齿弯曲强度计算式中首次引用。可以看出，YF是与齿轮形状的几何参数有关的一个系数。因为，根据齿轮形成原理，齿数的变化将引起轮齿上hF、SF、aF等参数的变化，由于hF、SF、aF均与齿轮模数成正比，致使齿形系数中的模数可以约去。因此，齿形系数不受模数的影响，而只与齿数有关，齿数越多YF越小，反之YF越大。这就是在机械设计的教材中经常会看到“标准齿轮的齿形系数只与齿数有关而与模数无关”的原因。
2. 2 齿轮压应力对弯曲应力的影响
根据30°切线法及齿轮受力分析。将法向力Fn移至轮齿中线并分解成相互垂直的两个分力，即圆周力Ft和径向力Fr。根据力学理论，Ft使齿根产生弯曲应力为σF，Fr则产生压应力σy。因此齿根危险截面上受到的应力为弯曲和压缩组成的组合应力，并导致齿根两边的应力大小不相等。然而，在相关的机械设计资料中都没有将由于径向力产生的压应力计算在齿轮的弯曲强度计算公式中，而且在大多数的相关教材中都认为: 压应力相对于齿根最大弯曲应力比较小，可以忽略不计。但是压应力到底多少，为什么可以忽略不计，很少有人进行计算，下面对压应力与弯曲应力进行探讨。如图2 中，Ft产生其弯曲应力σF如式( 1) 所示。由Fr产生压应力σy为σy = Fnsin αFbSF( 2)由式( 1) 及式( 2) 可得σyσF= SF6hFtan αF设OD = h'，则SF = 2h' tan30°，因此σyσF= tan 30tan αF3h'hF假设标准齿轮模数为m，齿数z。则齿顶圆压力角为cos αF = rbra= zz + 2cos α，由于h'hF< 1，因此，当不考虑h'hF的影响时，σyσF的大小取决于齿轮的齿数。为了便于讨论，取ξ = σyσF称为压应力对弯曲应力的影响系数。则根据计算可以得到ξ 与齿数的对应关系，如图3 所示。可见，压应力对弯曲应力的影响与齿数有关，而模数无关，而且随着齿数的变化而变化，齿数越少其影响越大，反之影响就越小，最终趋于一水平线。最小约为最大弯曲应力的8%，特别当h'hF< 1 时，压应力更小，可以忽略不计。这就是为了简化计算，在计算轮齿弯曲强度时一般只考虑弯曲应力的原因。从图2 可知，弯曲应力分为拉伸侧的拉应力和压缩侧的压应力。实际证明，拉伸侧是危险侧，因拉伸侧的`裂纹扩展速度较大。压缩侧有时虽裂纹出现较早，但发展速度较慢。所以大多数的公式以拉伸侧的应力作为设计时的计算应力。而且根据齿轮弯曲疲劳实验分析证明，考虑弯曲应力、压应力与只考虑弯曲应力的结果，实际上没有多大差别。因此，在齿轮弯曲疲劳强度计算中只考虑弯曲应力。
2. 3 齿面接触疲劳强度计算
图4 赫兹接触应力模型齿面接触疲劳强度计算是针对齿轮齿面疲劳点蚀失效进行计算的强度计算。1881 年，赫兹提出两个圆柱体接触时接触面上载荷分布公式，该式作为齿面强度计算的理论基础，如图4 所示。根据赫兹接触应力理论，在载荷作用下接触区产生的最大接触应力为σH = Fnπb·1ρ1± 1ρ21 - μ21E1+ 1 - μ22槡 E2( 3)式中，Fn为作用在圆柱体上的载荷; b 为接触长度;μ1、μ2分别为两圆柱体材料的泊松比; E1、E2为两圆柱体材料的弹性模量。ρ1、ρ2为两圆柱体接触处的半径，式中“+”号用于外接触，“-”号用于内接触。1898 年，拉塞根据法向力应用“压强”原理研究齿面的接触疲劳强度问题。1908 年，奥地利的维德基将赫兹的两个圆柱体的接触应力理论应用于计算轮齿齿面应力，并绘出了沿啮合线最大接触应力变化图。1932 年，英国BS 根据实验数据提出基础表面应力作为齿面强度计算方法。1940 年，美国AGMA 采用齿面强度最重负荷点的接触应力最大值计算方法。
1949 年，白金汉提出节圆上齿面接触应力不超过许用值的计算方法，后来该方法被许多计算方法所采用。1954 年，尼曼采用最大负荷点上滚动压力。至今，我国皆以赫兹公式作为计算齿面接触疲劳强度的理论基础，即以赫兹应力作为点蚀的判断指标。通常令1ρΣ= 1ρ1± 1ρ2，ρΣ称为综合曲率，对于标准齿轮，1ρΣ= 2d1 sin αi ± 1i 。并令式( 3 ) 中的ZE =1π 1 - μ21E1+ 1 - μ22E 槡为弹性影响系数。从而，获得渐开线直齿圆柱齿轮接触疲劳强度的基本公式为σH = ZEZH2KT1bd21i ± 1槡 i #[ σ ] H( 4) 式中，ZH = 2槡sin αcos α，称为区域系数，对于压力角α= 20°的标准齿轮，ZH≈2. 5。在机械设计手册或机械设计教材中，有关齿轮接触疲劳强度公式有很多版本，其中最常见的是将一对钢制标准齿轮齿面接触强度校核公式进行简化，取钢制齿轮的E1 = E2 =2. 06 ×105MPa，μ1 =μ2 =0. 3，便获得机械设计中常用的校核公式。σH = 671 KT1bd21i ± 1槡 i ≤[ σ ] H( 5)
2. 4 齿面胶合强度计算
齿轮另外一个常见的失效是齿面胶合。有关齿轮胶合比较统一的说法是: 相互啮合的两金属齿面，在一定的压力下直接接触发生黏着，同时又随着齿面运动而使金属从齿面上撕落而引起的黏着磨损现象。胶合分为冷胶合和热胶合。对于高速重载的齿轮传动，齿面瞬时温度较高，相对滑动速度较大，则容易发生热胶合。对于低速重载的重型齿轮传动，由于齿面间压力过大，导致齿面油膜被破坏，尽管齿面温度不高，但也容易产生胶合，称为冷胶合。
对于齿轮齿面胶合强度计算的研究，目前主要基于两种理论，一是基于Pv 值( 压力与速度的乘积) 或PTv ( T 为啮合点到节点的距离) 值作为计算胶合的指标。另一种是以齿面温度作为判定胶合的准则的布洛克算法。1975 年，温特提出积分温度法。现在ISO 的标准中主要以这两种方法为主。2003年，我国颁布“圆柱齿轮、锥齿轮和准双曲面齿轮胶合承载能力计算方法”国家标准( GB - Z 6413. 1 - 2003和GB - Z 6413. 2 - 2003)。该标准等同采用了ISO/TR 13989 - 2000“圆柱齿轮、锥齿轮和准双曲面齿轮胶合承载能力计算方法”。曾经有人试图以按弹性流体动力润滑理论计算齿面间的油膜厚度作为胶合的评判依据。
我国多数的机械设计教材中齿轮强度设计一般只提供齿面接触疲劳强度和齿根弯曲疲劳强度两种计算方法，并未提供有关齿面胶合的强度计算公式。
3 结束语
文中分别对机械设计教学中有关齿轮的强度设计问题进行了分析和探讨，详细解读我国齿轮强度设计的历史沿革及现状，以及齿轮强度设计计算过程中让人困惑的问题及解决方法。研究指出，在齿轮弯曲疲劳强度的计算中，压应力对弯曲应力的影响是有限的，一般可忽略不计，只有当需要精确计算时，应当考虑其影响。论文的研究可以帮助齿轮设计人员和学生更好地理解齿轮设计中的相关内容，为将来从事机械设计工作打下良好的基础。
小模数齿轮齿形误差图像测量法权转菊 (西安东风仪表厂计量处 710065) 摘要:本文提出了一种小模数齿轮齿形误差测量的新方法,该方法在极坐标系下采集齿廓边缘 摘要 点的坐标值,通过测量模型计算获得齿形误差,符合齿形误差定义,具有较高的精度. 关键字:齿形误差 关键字 光学测量 极坐标 数学模型 引言小模数齿轮尤其是模数在 0.05mm～0.5mm 的小模数齿轮广泛应用于航天航空,国防,IT,钟 表等领域的精密仪器仪表制造中.作为关键的运动传动件,其质量直接影响到仪器仪表的运动精度, 噪音,寿命等.因此,实现对小模数齿轮的高精度测量是保证仪器仪表质量的一个关键技术问题. 小模数齿轮由于其模数小而齿数通常较多,齿槽空间很小,很难采用传统的齿轮测量技术和仪 器.目前,普遍采用的测量方法是轮廓投影比较法和分度盘展成法.轮廓投影比较法即在轮廓投影 仪上,将齿轮与标准放大图进行比较,从而判定加工齿轮的齿廓精度,这种方法显然不能实现精确 检测.分度盘展成法测量效率低,受找正误差,分度误差的影响精度并不是很高. 近年来,随着光 学坐标测量机的应用和发展,基于 CCD 技术的齿轮测量方法的研究不断增多.本文作者研究了在光 学坐标测量机上对小模数齿轮齿形误差进行精密测量的一种新方法. 1 数学模型的建立 1.1 展开角增量与展开弧长增量的关系 按渐开线形成原理,渐开线上某一点的曲率半径 ρ 等于基圆上形成渐开线的起点 A 到曲率半 径 ρ 与基圆切点 B 间的弧长,ρ 也即展开弧长, 其展开角 w 与 ρ 之间的关系为: w =ρ/ r0 (式 1) ? 式中: r0 为齿轮的基圆半径 当展开角 w1 增加w 转角时,展开弧长的增 量为ρ. 与ρ 之间有一定的比例关系, w 如当 齿轮转动一度, ρ1 应增加齿轮基圆圆周长度为 则 1/360,所以得如下关系式: ρ=2πr0w/360=π/180wr0 式中:ρ 为展开弧长增量 ( 2 B1 A1 A2 ?A 式 2) 图 (1) 渐 开 线 形 成 原 理 1 ρ=ρ2-ρ1 1.2 极坐标系下展开角与极径的关系 按照几何关系,从图中可以看出: ρ= R2 r 0 2 2 R12 r 0 2 (式 3) wx=w2-w1=B2+cos-1 r0 r -B1- cos-1 0 R2 R1 r0 r - cos-1 0 R2 R1 (式 4) 也即:wx=B2-B1+cos-1 将式 3 和式 4 代入式 2 得: 2 R2 r 0 R12 r 0 2 2 =π/180r0 2-B1+cos-1 (B r0 r - cos-1 0 ) (式 5) R2 R1 从式 5 中可以看出, 如果我们以齿轮中心为极坐标中心, 靠近渐开线起始点测量一点作为极坐 标起点,建立极坐标系,在此坐标系下齿形上各点极径与极角应满足式 5 中的关系. 2 齿形误差的测量 由于齿形误差的影响,实际齿廓上各点的坐标值与理论值有差异,即相对一展开角实际齿廓上 展开弧长与理论值有差异.变换式 5,我们可以求出实际齿轮左右齿廓上这种差异. f= Ri r 0 2 2 R12 r 0 -π/180r0 i-B1|+cos-1 (|B 2 r0 r - cos-1 0 ) Ri R1 (式 6) 这种差异也是齿廓上各点曲率半径与渐开线上相应的理论曲率半径的差异. 最小值与最大值之 差即为包容实际齿形的两条最近的理论渐开线间的法向距离,符合 GB10095-85 规定齿形误差ff 定 义(见图 2) ,则齿形误差为: ff=fmax-fmin (式 7) 齿 顶倒 棱高 度 ? 工 作部 分 设 计 齿形 齿 根 起始 工作 圆 图 1 齿形 误差 示意 图 2 3 测量实例选用有背光照明系统和 CCD 视像头的坐标测量机对一模数为 0.5mm,齿数为 50 的小模数齿轮 进行实际测量.首先以齿轮中心为坐标原点,以齿廓上大于基圆半径一点为起点建立极坐标,然后 对齿廓进行测量,求得各点的极径及极角,通过数学模型计算齿形误差.测量数据及结果见表 1. 表1 右齿 R(mm) ω(°) 9.518 9.655 9.711 9.806 9.898 9.958 360 359.813 359.743 359.570 359.404 359.284 f (mm) 0 -0.0037 -0.0016 -0.0046 -0.0043 -0.0043 R(mm) 9.539 9.602 9.683 9.750 9.924 9.938 齿形误差测量结果 左齿 ω(°) 0 0.072 0.192 0.299 0.613 0.627 f (mm) 0 0.0002 -0.0013 -0.0019 -0.0020 0.0002 齿形误差:0-(-0.0046)=0.0046 齿形误差:0.0002-(-0.0020)=0.0022 4 结论本文结合渐开线展成原理对极坐标系下渐开线上各点坐标关系进行了分析, 并给出了数学模型, 由此得出齿形误差测量方法.通过测量实例对测量方法进行了说明. 这种方法与传统的使用分度盘 测量齿形的方法相比,同样是图像测量法,但由于少了分度盘的找正误差,分度误差等误差影响因 素,测量精度大为提高,并且可利用坐标测量机柔性定位功能,形成测量程序进行批量测量,实现 对小模数齿轮齿形误差的高精度,高效率测量. 3
郭隐彪,黄元庆,粱锡昌. 免修整柔性金刚石砂轮的磨削性能研究. 中国机械工程, Vol.12, No.7, 2001: 775~777郭隐彪,黄元庆,田波,庄司克雄. 非轴对称非球面平行磨削误差补偿技术研究. 机械工程学报, Vol.38,No.5,2002:118~121 （EI收录02447178250）郭隐彪,黄元庆,粱锡昌. 超精密研磨砂轮的优化设计. 重庆大学学报, Vol.25, No.2, 2002 : 8~10郭隐彪,张翊,庄司克雄. 研磨轮结构参数对工件加工精度的影响.中国机械工程, Vol.13,No.20,2002: 1724~1726 （EI收录02507271598）张翊,郭隐彪,庄司克雄. 超精密平面研磨加工参数对精度的影响. 制造技术与机床, 2002(10):37~39厨川常元，郭隐彪，庄司克雄.Study of Ultra-precision Plane Honing. Japan Society of Abrasive Technology, Sizioka, Japan, 2002: 177~178郭隐彪,徐雅珍. 机电学科新课程建设的几点思考. 厦门大学学报（哲学社会科学版） 增刊, 教学论文专辑: 168~170洪永强,郭隐彪,蒋文芳. 基于CCD的高精度非球面测量与控制技术. 仪器仪表学报 2002, No.3（增刊）:184~185仲进安,郭隐彪,张翊. 超声波微加工间隙在线测控技术研究. 仪器仪表学报 2002, No.3（增刊）: 367~368张翊,郭隐彪,庄司克雄. 微小振动影响超精密非球面加工精度的研究. 金刚石与磨料磨具工程, No3,2003:17~20（EI收录03327585136）郭隐彪,张翊,徐雅珍,庄司克雄. 圆弧金刚石砂轮的修形及修锐技术研究. 制造技术与机床, No4, 2003: 39~40王振忠,郑琳,郭隐彪. 超精密磨床系统振动监控及平衡技术研究. 机电技术 增刊, 省第三届学术年会装备制造业分会论文集, 2003: 128~131纪跃波,郭隐彪. 虚拟式多功能时频分析系统的研究与开发. 机电技术 增刊, 省第三届学术年会装备制造业分会论文集, 2003: 132~137仲进安,郭隐彪. 光学非球面自动对焦测量系统研究. 机电技术 增刊, 省第三届学术年会装备制造业分会论文集, 2003: 138~141谢仁宁,郭隐彪,黄浩. 基于高精度非球面测量的数据采集与处理系统研究. 机电技术 增刊, 省第三届学术年会装备制造业分会论文集, 2003: 205~208黄浩,郭隐彪,谢仁宁. 高精度非球面制造系统控制软件的设计. 机电技术 增刊, 省第三届学术年会装备制造业分会论文集, 2003: 209~212黄浩,郭隐彪,王振忠,谢仁宁. 轴对称非球面加工进给速度控制技术研究. 金刚石与磨料磨具工程, No.2, 2004: 20~24（EI收录04308283557）黄浩,郭隐彪,谢仁宁. 高精度非球面制造系统控制软件的设计. 厦门大学学报(自然科学版), Vol.43, No.2, 2004: 179~182范剑红,周学才,郭隐彪,邹琴. 激光共焦扫描显微镜运动控制系统设计. 重庆大学学报, Vol.27, No.3, 2004: 9~12谢仁宁,郭隐彪,黄浩. 基于高精度非球面测量的数据采集与处理系统研究. 厦门大学学报（自然科学版）, Vol.43, No.3, 2004: 339~343朱立秒,龚子伟,郭隐彪. 基于远程制造的多功能数控接口设计及实现. 厦门大学学报（自然科学版）, Vol.43, No.5, 2004: 628~631杨小潘,郑琳,郭隐彪. 超精密磨削加工微小振动模拟系统研究. 组合机床与自动化加工技术,No.6,2004:1~张翊,郭隐彪,庄司克雄. 超精密平面研磨加工速度对精度的影响. 制造技术与机床, No.7, 2004: 56~59魏丽珍,郭隐彪,苏秋丽,龚子伟. 非球面镜超精密加工系统控制软件的研究. 制造技术与机床,No.8, 2004: 111~114谢仁宁,郭隐彪,黄浩. 基于非球面加工的砂轮表面误差测量技术研究. 中国机械工程,Vol.15,2004:4~6黄浩,郭隐彪,王振忠,谢仁宁. 高精度非球面加工砂轮均匀磨损技术研究. 中国机械工程,Vol.15,2004:25~29郭隐彪,兰劲,龚子伟. 远程制造的多功能数控接口技术. 先进制造技术与齿轮学术会议中国齿轮制造研究会第七届学术年会论文集2004年11月6－7日重庆,pp:10~14郑杰,郭隐彪. 嵌入式系统中的低功耗设计. 机电技术, 2004增刊, 省第四届学术年会论文集, 2004: 69~72林荣川,郭隐彪. 液压缸临界载荷计算. 机电技术, 2004增刊, 省第四届学术年会论文集, 2004: 89~92王振忠,郭隐彪. 虚拟制造技术浅析. 机电技术, 2004增刊, 省第四届学术年会论文集, 2004: 140~145董亮,郭隐彪. 空调冷水机组参数的分析. 机电技术, 2004增刊, 省第四届学术年会论文集, 2004: 245~249邹琴,郭隐彪,范剑红. 一种开放式智能化数控系统研究. 厦门大学学报(自然科学版), Vol.43, No.6, 2004: 797~801邹琴,郭隐彪. 多功能数控砂轮修整系统. 2004年中国机械工程学会年会论文集.2004：142~146.郭隐彪,兰劲,龚子伟. 远程制造的多功能数控接口技术. 现代制造工程 No.1, 2005:42~45仲进安，郭隐彪. 高精度非球面表面测量技术研究. 机械工程学报，Vol.41，No.1，2005：176~180（EI收录：05149025347）邹琴, 郭隐彪. 多功能数控砂轮修整系统. 机电21世纪，No.3，2005：34~36
