机器人用RV传动必须具有高的运动精度和小的回差外,还必须具有很高的刚度。这也是它与谐波传动相比,在机器人传动中最突出的优点之一。RV传动的低速级是采用多齿啮合的摆线针轮传动,其接触刚度之高是众所周知的。高速级则采用有3个行星轮的渐开线行星齿轮传动,不仅提高了承载能力,同时也提高了刚度。除此之外,输出机构为刚架结构的行星架,也具有很高的刚度。机器人用RV传动对扭转刚度有严格的要求,通常在额定转矩下,由扭转弹性变形引起的弹性回差不超过4′~5′。
1)采用常规的力学方法对RV传动的每个组成部分进行刚度分析计算
①摆线针轮传动部分的弹性变形。摆线针轮传动部分的变形主要由针齿与摆线轮齿的接触变形及针齿与针齿壳上针齿孔的接触变形两部分组成。
由(2)的受力分析理论,可以得出在额定转矩的作用下,针齿与摆线轮及针齿壳这两部分的总接触变形的最大数值wmax,由此,假定摆线轮固定时,由于此接触变形而引起针轮产生的转角为
式中 r′p——针轮节圆半径。
将针轮弹性转角折算到输出轴上,可得当固定渐开线中心轮时,由于针齿与摆线轮齿的接触变形及针齿与针齿壳的接触变形所引起的输出轴的转角Δθc为
Δθc=Δθpi165 (11-139)
式中 i165——当渐开线中心轮固定时,输出轴相对于针轮的传动比,即
②渐开线传动部分的弹性变形 按照ISO刚度计算方法,首先计算单齿刚度,然后根据单齿刚度计算其啮合刚度。对于刚性啮合齿轮,在中等载荷作用下,其单对齿刚度按下式近似计算:
式中 q——单位齿宽柔度(mm·μm/N)。
当x1≥x2,-0.5≤x1+x2≤2.0时,q按下式计算:
式中 x1——小齿轮的变位系数;
x2——大齿轮的变位系数;
zv1——小齿轮的当量齿数,对直齿轮zv1=z1;
zv2——大齿轮的当量齿数,对直齿轮zv2=z2。
根据ISO啮合刚度计算B法,考虑齿轮啮合时的重合度的影响,齿轮的啮合刚度为
cr=(0.75εα+0.25)c′ (11-142)
式中 εα——端面重合度;
c′——单对齿刚度[N/(mm·μm)]。
根据刚度的定义,啮合点的位移为
式中 Fgt——中心轮受的啮合力的切向分力(N);
bg——渐开线齿轮宽度(mm);
α′——渐开线齿轮的啮合角(°)。
由啮合点位移使中心轮产生的转角为
式中 r′1——渐开线中心轮的节圆半径(mm)。
当固定针齿壳及中心轮时,由于渐开线齿轮传动部分的弹性变形而引起输出轴的转角为
式中 i516——当针轮固定时,输入轴相对于输出轴的传动比。
③行星架的弹性变形。行星架是由两块侧板1、2和三个截面近似梯形的连接柱销3组成。行星架的弹性变形是指在转矩T的作用下两块侧板的相对转角,如图11-61所示,它也可以表示为在半径为rn圆周方向上由切向力Fnt引起的位移Δ的大小,而切向力
式中 FG1t——曲柄作用在侧板上的切向力(N);
rn——连接柱销形心到行星架中心的距离(mm);
a0——曲柄支撑轴承中心到输出轴中心的距
离(mm),其数值等于渐开线行星齿轮传动中心距。
图11-61 行星架示意图
参照图11-62行星架展开图,切向位移量Δ可表示为
式中 αm1——切向力Fnt对侧板1的弯曲变形和剪切变形的影响系数;
αm2——切向力Fnt对侧板2的弯曲变形和剪切变形的影响系数;
αn——切向力Fnt对连接柱销弯曲变形和剪切变形的影响系数;
k0——两块侧板的刚度比较系数;
Lm——侧板1/3段的弧长(mm),;
Ln——两块侧板中心平面间的距离(mm);
E——弹性模量(MPa)。
图11-62 行星架截面展开图
且
式中 βmi——两圆盘形侧板的形状系数(i=1,2);
βn——凸四边形连接柱销的形状系数;
kmi——两侧板的横截面形状系数(i=1,2);
kn——连接柱销的横截面形状系数;
Imi——两侧板相对于ym轴的截面惯性矩(i=1,2)(mm4);
In——连接柱销对于yn轴的截面惯性矩(mm4);
Smi——两侧板的横截面积(i=1,2)(mm2);
Sn——连接柱销截面积(mm2);
lmi——相当悬臂变形的两侧板有效长度(i=1,2)(mm);
ln——相当悬臂变形的连接柱销件的有效长度(mm);
μ——泊松比,μ=0.3。
式中系数由有关图表确定,惯性矩和面积是按照图11-62行星架展开图尺寸计算的。由此,可按式(11-147)计算出行星架两侧板相对位移Δ值,因此,由行星架弹性变形引起输出轴转角为
④曲柄轴的弹性变形。曲柄轴的周向弹性变形由曲柄轴的周向弯曲变形和扭转变形两部分组成,下面分别进行计算。
a.曲柄轴的周向弯曲变形。由本章11.2节关于传动曲柄及行星轮的受力分析,曲柄轴受两摆线轮在轴承孔分布圆周切向力的作用(参见图11-49),引起曲柄有此方向上的弯曲变形,从而导致输出轴产生一定的转角。根据材料力学梁的挠度计算方法,对于变截面的曲柄轴,为了提高计算精度,本文采用有限插分法计算挠度。在突变截面附近选择一些点增加差分点数,如图11-63所示,截面3和截面4处的位移即为所求挠度。
图11-63 曲柄受力变形图
这里,摆线轮给曲柄轴的圆周方向的作用力为
式中 Tc——一片摆线轮承受的转矩(N·mm)。
曲柄轴上的行星轮受到的切向力
式中 T1——输入的转矩(N·mm)。
由于曲柄在截面3和截面4处的位移不同,在计算曲柄变形引起输出轴转角时,应以它们的平均值fm来计算
式中 f3——截面3处的位移(mm);
f4——截面4处的位移(mm)。
由此,由曲柄轴的周向弯曲变形引起输出轴的转角为(www.daowen.com)
图11-64 曲柄受转矩作用图
b.曲柄的扭转变形。如图11-64所示,曲柄受转矩作用时,截面0与截面3之间将产生一定的相对转角ΔφS
式中 G——切变模量(MPa);G=8×104MPa;
T2i——曲柄轴各段转矩(N·mm),而
T——输出轴转矩(N·m);
Jpi——极惯性矩(mm4),;
di——曲柄轴各段轴直径(mm);
lzi——曲柄各段轴长(mm)。
因此,由于曲柄的扭转变形产生的输出轴的转角
由式(11-155)及式(11-157)可得,由于曲柄的变形引起输出轴的转角为
ΔθS=ΔθS1+ΔθS2 (11-158)
⑤轴承的弹性变形
a)轴承的变形与外载荷的关系 轴承内外圈与滚子接触处的总接触变形为
δz=kzQ0.9 (11-159)
式中 Q——滚子所受压力(N);
kz——刚度系数,kz=7.66×10-5/L0.8;
L——滚子有效长度(mm)。
图11-65 轴承受力分析图
图11-65是轴承受径向力作用时的变形情况。在径向力Fzr作用下,轴承内外圈相对位移为δr,为计算方便将最下方的滚子的序号记作0,依次向上编为1,2,…,n等,相邻滚子所对圆心角为φz。在Fzr作用下,各个滚子的受力大小及变形量各不相同,序号为j的滚子受力为Qj,与内外圈接触变形量的和为δj,设轴承游隙为Δr,则根据几何关系,有
式(11-160)称为轴承弹性变形协调条件,由方程组中的不等式可确定受力滚子的最大序号n,n又称为分布系数,2nφz为轴承的负荷分布角。在圆心角2nφz范围内的滚子承受载荷,并产生接触变形。当轴承有预紧量时,Δr为负值,轴承的分布角大于π。除了预紧或游隙影响轴承的负荷分布外,外载荷也影响负荷的分布。在有游隙时,外载的增大会引起分布角的增加。相反,在有预紧量时,外载的增加引起负荷分布角的减小。将式(11-159)代入式(11-160)得
轴承不受外力作用时,各滚子受力一样,与内外圈作用的合力为0。当轴承受外力作用时,各滚子的受力情况发生变化,滚子与内圈作用的合力与外载达到新的平衡状态,其平衡条件为
为了找出轴承内外圈的相对位移量与外载荷的关系,而位移量直接与内力大小有关,内力的大小又与外载、轴承结构、滚子及滚道的变形有关。因此,这种关系的计算不能直接套用公式,必须利用迭代的方式逐步逼近其真实解。
对于滚子轴承,先给出最大负载的初值
Q0=4.6Fzr/z (11-163)
式中 Fzr——实际外载荷(N);
z——滚子数量。
由式(11-159)计算出受载最大的滚子的变形量δ0,由δ0和游隙或预紧Δr根据变形协调条件求出相对位移量δr及各滚子的变形量δj和分布系数n。根据变形量δj和滚子负载Qj的关系,可以求出每一个滚子的负载Qj,再由平衡条件求出在假设滚子最大负荷Q0的情况下的外载荷Fz′r。由于在求初值Q0时没有考虑实际游隙或预紧量,由此求出的载荷F′zr与实际载荷Fzr不同,其误差用下式表示:
F′zr-Fzr/Fzr<ε (11-164)
式中 ε——收敛精度。
若精度不符合要求,需要进行下一轮迭代,按F′zr-Fzr的大小修正Q0的初值
Q*0=Q0±F′zr-Fzr/2 (11-165)
式中当Fz′r>Fzr时取负号,反之Fz′r<Fzr时取正号。Q*0为新一轮迭代的初值。
曲柄支撑轴承为圆锥滚子轴承,目前尚缺乏有效的简易算法进行计算。我们仍利用分析圆柱滚子轴承的方法,对其变形及受力进行合理的变换。βz为圆锥滚子轴线与轴承轴线之间的夹角,由于轴承轴向固定,所以在径向力Fzr作用下,其径向位移δj分为滚子轴向滑动位移δjt和法向位移δjn,其中的法向位移δjn=δi/cosβz是法向力Qj产生的接触变形,由此外载在计算程序中应进行如下的相应变化
b.轴承变形引起输出轴的转角。根据轴承受力变形计算程序,可以算出曲柄支撑轴承及转臂轴承的内外圈相对位移量,从而求出轴承接触变形引起输出轴的转角。
曲柄支撑轴承周向变形引起输出轴的转角Δθb1计算如下:
固定摆线轮,在输出轴上加转矩T,则曲柄支撑轴承受到的圆周方向作用力为
式中 a0——曲柄支撑轴承中心到输出轴中心的距离,其数值等于渐开线行星齿轮传动中心距。
采用上述计算程序,可计算出支撑轴承有转矩T的作用下内外圈产生的相对位移量δr1,则在轴承孔中心圆的周向上产生的转角为
转臂轴承周向变形引起输出轴的转角Δθb2计算如下:
由摆线轮承受的转矩Tc作用在转臂轴承上的周向力FA2可以计算出,当固定输入轴、输出轴时,由转臂轴承变形引起摆线轮在半径为a0的圆周上切向位移量δr2,则引起摆线轮的转角为
随摆线轮的转动,针轮产生相应的转角Δθ5为
式中 iH45′——当输出轴、输入轴固定时,摆线轮相对于针轮的传动比。由此可得,当针轮固定时,由转臂轴承的变形引起输出轴的转角为
式中 i165——当输入轴固定时,输出轴相对于针轮的传动比。
根据以上的分析,曲柄支撑轴承与转臂轴承的变形引起输出轴的转角为
Δθb=Δθb1+Δθb2 (11-172)
以机器人用RV-60AⅡ减速器为例,以上五种弹性变形引起的输出轴转角见表11-31。
表11-31 弹性变形引起输出轴的转角
然后可得RV-60AⅡ减速器的刚度k为
从计算结果看出,在额定转矩作用下,轴承的变形最大,刚度最小,正是基于这一实际情况,所以在本章11.2节中对RV传动整体受力分析时,以此为基础,假定摆线轮与曲柄轴刚度很大,而轴承可看成弹性体。从而总结出将RV传动受力的静不定问题简化成静定问题求解的力分析方法及公式。
由于轴承的刚度很小,难于提高,而机器人传动对RV减速器整体的抗扭刚度又有很高的要求,所以必须尽可能增大其他零件的刚度,以弥补轴承刚度的不足。主要措施是:①摆线轮齿形采用负移距与负等距组合修形,可以保证多齿啮合的共轭齿形,增大了摆线轮与针轮齿啮合的刚度;②针齿不用两支点而是用半埋齿以消除其弯曲变形;③输出轴采用刚性尽可能大的大直径圆盘输出。
2)整体刚度的有限元分析RV传动系统是由渐开线行星传动与摆线针轮行星传动组成的一个完整的封闭式差动机构体系,如果只把某个零件从机构中单独取出来进行有限元计算,那么,机构体系中几个重要计算对象之间的相互影响就很难准确地反映到计算模型中去,计算结果也将失去一些精度。为避免上述缺陷,在计算模型中应将摆线轮、行星架、偏心轴、针轮等同时纳入一块参与计算,这样才能使计算结果更加准确可靠,图11-66为RV传动结构系统的整体有限元计算模型。
传动中的摆线轮以自己的回转轴作平面圆周运动,这个回转运动是通过摆线轮内的偏心轴承实现的。轴承外环即为摆线轮内孔,它的网格划分沿环向是均匀的,在每一个接触点对上用“点对点”的接触单元模拟它们的弹性接触。对轴承滚柱以及内环,则将其凝聚到曲柄点上,用I-DEAS的刚性元及自由度释放技术,可以模拟轴承的“铰接”作用,从而实现了摆线轮与曲柄之间的铰接。
考虑到两片摆线轮是主要传递载荷零件,摆线轮齿部分网格划分密一些,可以真实地反映局部的应力与变形情况,两片摆线轮均采用三角形板单元。摆线轮与各针轮齿啮合情况,采用“点对地”的接触单元模拟它们的弹性接触。单元方向为摆线轮与各针轮齿啮合作用力的方向,即通过摆线轮与针轮啮合节点P。图11-67为摆线轮与针轮的有限元计算模型。
曲柄轴的偏心简化为刚性元以模拟轴的曲柄,这些空间刚性元的尺寸取决于曲柄轴的偏心距的数值。所以,将曲柄轴简化为空间梁元后,在曲柄轴与摆线轮的偏心轴处分
图11-66 RV传动系统的有限元计算模型
图11-67 摆线轮与针轮的有限元计算模型
别安放了若干个上述偏心值为其长度的刚性元。图11-68是行星架与曲柄轴的有限元计算模型。
图11-68 行星架及曲轴的有限元计算模型
选用美国SDRC公司的I-DEAS软件系统,按此模型计算的渐开线齿轮传动、曲柄轴及行星架总体结构扭转刚度与前几节提出的计算模型所得结果进行对比列于表11-32中(以机器人用RV-60AⅡ减速器为例)。
表11-32 扭转变形计算结果
[1]因回转具有方向性,因此称回转矢量,常在机器人上使用。
免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。