2. 四川大学 机械工程学院,四川 成都 610065
2. School of Mechanical Eng., Sichuan Univ., Chengdu 610065, China
无级变速器(continuously variable transmission, CVT)是一种理想的机械传动方案,因为无级变速器可以提供一系列的连续传动比使得车辆发动机始终处于高效点工作,这就极大地提高了整车的效率,减小了车辆排放,达到节能减排的作用[1-4]。CVT具有动力高、油耗低、排放少的特点,据统计,采用CVT的汽车比采用自动变速器(AT)的汽车油耗减少7%~15%,废气排放降低约10%,生产成本降低20%~30%;而采用CVT传动系的混合动力车的油耗可能减少30%,排放可减少45%[5-6]。牵引式CVT具有功率密度高、效率高、传递功率大等特点,将会成为未来无级变速器技术的发展方向,也有望在其他领域推广[7]。因此,牵引式CVT技术在日本、英国、美国等发达国家已经被广泛研究。
牵引式CVT主要由变速机构、加载机构与传动机构3部分组成,实现变速的关键是变速机构。现代工业技术的发展,希望CVT的变速机构对于负载改变能够实现自动/自适应变速;而目前自适应变速还多采用电控方式实现,其工作原理是通过监测负载力/力矩、整车运动速度/加速度等信号,将其反馈给控制器,控制变速机构动作实现变速控制。然而,电控方式需要复杂的控制电路和检测传感器,成本较高。因此,设计出可根据负载而自动变速的自适应变速机构,对于提高牵引式CVT的结构紧凑性,降低控制复杂度等具有较大推动作用。在此背景下,作者将针对最典型的同轴移动类牵引式CVT—滚轮平盘牵引式CVT,提出了一种基于螺纹−弹簧力封闭的自适应变速机构,并对此CVT进行了运动学计算和动力学分析建模,以给牵引式CVT自适应变速研究提供新的思路。
1 自适应变速机构描述如图1所示,将传统的滚轮平盘牵引式CVT[8]的输出轴与滚轮之间通过滚珠螺纹传动连接,在滚轮的另一端设置力封闭弹簧机构,其在某一位置可与螺纹轴向力平衡,滚珠螺纹传动经牵引传动到力封闭弹簧,构成了滚轮平盘牵引式CVT的自适应变速机构。其工作原理为:动力从输入轴输入,经加载凸轮机构传递至平盘,平盘与滚轮之间压紧,通过牵引运动传递动力,将动力从平盘传递至滚轮,滚轮通过滚珠螺纹机构将动力传递至输出轴输出;如忽略接触处的摩擦力,稳定传动时,滚珠螺纹机构轴向推力与弹簧弹力平衡,当输出轴负载发生改变时,滚珠螺纹机构轴向推力改变,平衡打破,滚轮沿着平盘径向移动,直到达到新的平衡状态实现稳定传动;而新平衡状态下,接触点位置改变,实现了变速。
![]() |
1. 平盘;2. 滚轮;3. 输入轴;4. 输出轴;5. 力封闭机构;6. 加载凸轮机构;7. 滚珠螺纹机构。 图1 滚轮平盘牵引式CVT自适应变速机构 Fig. 1 Shifting mechanism model of roller flat disk traction-type CVT |
为了便于计算,图1中各参数说明如下:力封闭弹簧选用刚度系数
1)理论传动比
假设接触点
$i = \frac{{{\omega _1}}}{{{\omega _2}}}=\frac{{{r_2}}}{{{r_1}}} = \frac{{{R_2}}}{{{R_1} - x}}$ | (1) |
式中,
2)自旋角速度
图2自旋角速度示意图。平盘角速度为
![]() |
图2 自旋角速度 Fig. 2 Spin angular velocity |
平盘与滚轮之间的相对角速度为:
$ {\omega }_{21}={\omega }_{2}-{\omega }_{1} $ | (2) |
则自旋角速度可表示为:
${\omega _{{2}{{1}_{{\rm{spin}}}}}} = {\omega _{21}} \cdot \sin \; \alpha = {\omega _1}$ | (3) |
3)滑移率
牵引装置传递扭矩时,其实际的传动比与理论传动比会有所差异,这是因为在牵引传递过程中存在一定的滑移现象。所谓滑移,即为两个滚动元件在接触点表面的速度相对差值[9]。而这个微小的速度差产生的本质原因在于牵引传动接触点传递扭矩的实质是弹流润滑油膜的剪切作用。
$ {C_{\rm{r}}} = \frac{{{V_1} - {V_2}}}{{{V_1}}} = \frac{{\left| {{\omega _1}} \right|{r_1} - \left| {{\omega _2}} \right|{r_2}}}{{\left| {{\omega _1}} \right|{r_1}}} $ | (4) |
当考虑滑移现象时,平盘和滚轮之间的角速度关系可表示为:
$ \frac{{{\omega _2}}}{{{\omega _1}}} = \frac{{{R_2}}}{{{R_1} - x}}(1 - {C_{\rm{r}}}) $ | (5) |
4)法向压力与有效牵引系数
如图1所示,平盘与滚轮之间的加载力由加载凸轮机构提供,加载凸轮通过中间球体与端面接触产生加载力。假设
根据Zhang等[10]提出的端面凸轮法向力计算模型,可得:
${L_{\rm{c}}} = 2\text{π}{R_c}{\rm{tan}}\;\lambda $ | (6) |
${T_{{\rm{in}}}} = {F_{{\rm{ct}}}}{R_{\rm{c}}}$ | (7) |
凸轮结构受力分析:
${F_{{\rm{ct}}}} = {N_{\rm{c}}}(\sin\;\lambda + {\mu _{\rm{c}}}\cos\;\lambda )$ | (8) |
$ {F_{{\rm{ca}}}} = {N_{\rm{c}}}(\cos\;\lambda - {\mu _{\rm{c}}}\sin\; \lambda ) $ | (9) |
式中,
由式(5)~(8)可得:
${F_{{\rm{ca}}}} = \frac{{2\text{π}{T_{{\rm{in}}}}(\cos\;\lambda - {\mu _{\rm{c}}}\sin \;\lambda )\tan\;\lambda }}{{{L_{\rm{c}}}(\sin\;\lambda + {\mu _{\rm{c}}}\cos\;\lambda )}}$ | (10) |
如果忽略摩擦,即
${F_{{\rm{ca}}}} = \frac{{2\text{π}{T_{{\rm{in}}}}}}{{{L_{\rm{c}}}}}$ | (11) |
可得平盘与滚轮之间接触的法向压力为:
$N = {F_{{\rm{ca}}}} = \frac{{2\text{π}{T_{{\rm{in}}}}}}{{{L_{\rm{c}}}}}$ | (12) |
稳定运动状态下,平盘与滚轮之间的牵引力由EHD液体油膜产生,可近似表达为:
$ {F_{{\rm{tr}}}} = {\mu _{\rm{e}}}N = {\mu _{\rm{e}}}\frac{{2\text{π}{T_{{\rm{in}}}}}}{{{L_{\rm{c}}}}} $ | (13) |
其中
$ {\mu _{\rm{e}}} = \frac{{{L_{\rm{c}}}}}{{2\text{π}({R_1}{\rm{ - }}x{\rm{)}}}} $ | (14) |
根据赫兹接触理论,当平盘与滚轮之间有法向加载力,则接触点处应为椭圆[11-12]。椭圆的方向和尺寸与接触点的法向力、平盘与滚轮的材料特性和结构尺寸均相关,其接触应力在接触区内按照椭球体规律分布,表达如下:
$ p={p}_{\rm{H}}\left(1-\frac{{x}^{2}}{{b}^{2}}-\frac{{y}^{2}}{{a}^{2}}\right) $ | (15) |
最大赫兹接触应力
${p_{\rm{H}}} = \frac{{3N}}{{2\text{π}ab}}$ | (16) |
式中,
在Hamrock等[13-14]对于椭圆接触问题等温弹流润滑的研究中,采用了Harris提出的迭代算法,该方法可以较高精度地确定a、b的值,如下:
$ a = {\left(\frac{{6{{{k}}^2}\varepsilon RN}}{{\text{π}{{E^{'}}}}}\right)^{1/3}} $ | (17) |
$b = \frac{a}{k}$ | (18) |
R为当量曲率半径,结合本文中平盘滚轮式CVT结构,可得:
$\frac{1}{R} = \frac{1}{{{R_2}}} + \frac{1}{{{r_{\rm{r}}}}}$ | (19) |
$ {{E^{'}}}=\frac{{{\rm{1 - }}{\nu _{1}}^{2}}}{{{E_1}}} + \frac{{{\rm{1 - }}{\nu _{2}}^{2}}}{{{E_2}}} $ | (20) |
式中,
根据Brewe等[15]从迭代方法中采用数据处理回归出的近似公式来分别计算
$ {k}={1.033\;9\left(\frac{{R}_{y}}{{R}_{x}}\right)}^{0.636} $ | (21) |
$ \varepsilon=1.000\;3+\frac{0.596\;8}{\dfrac{{R}_{y}}{{R}_{x}}} $ | (22) |
将各部件分离,进行受力分析,为了简化方程数量,将固定连接部件视为整体部件,如图3所示。图中省略轴承、螺纹、键等连接结构表达,仅给出分析相关的力的表达。
![]() |
图3 分离体受力模型 Fig. 3 Separate force model of the roller flat disk traction-type CVT |
如图3(a)所示,对输入轴受力分析:
${F_{{\rm{ib}}}} - {F_{{\rm{ca}}}} = 0$ | (23) |
$ {T_{{\rm{in}}}} - {F_{{\rm{ct}}}}{R_{\rm{c}}} = {J_{{\rm{in}}}}{\dot \omega _{{\rm{in}}}} $ | (24) |
式中,
如图3(b)所示,对平盘受力分析:
$ {F_{{\rm{ca}}}} - N = 0 $ | (25) |
$ {F_{{\rm{ct}}}}{R_{\rm{c}}} - {F_{{\rm{tr}}}}{r_1} = {J_1}{\dot \omega _1} $ | (26) |
式中,
如图3(c)所示,对滚轮受力分析:
${F_{{\rm{sa}}}} - {F_{{\rm{sf}}}} - {F_{{{\rm{k}}_1}}} = {m_2} \cdot \ddot x$ | (27) |
$ {F_{{\rm{tr}}}}{r_2} - {T_{\rm{s}}} = {J_2} \cdot {\dot \omega _2} $ | (28) |
式中,
如图3(d)所示,输出轴受力分析:
$ {F_{{\rm{ob}}}} - {F_{{\rm{sa}}}} = 0 $ | (29) |
${T_{\rm{s}}} - {T_{{\rm{Load}}}} = {J_{{\rm{out}}}}{\dot \omega _{{\rm{out}}}}$ | (30) |
式中,
根据螺纹传动原理,有:
${F_{{\rm{sa}}}}=\frac{{{2}{T_{{\rm{Load}}}}}}{{{D_{\rm{s}}}\tan (\varphi + \rho )}}$ | (31) |
将
$ {m_2} \cdot \ddot x + {F_{{\rm{sf}}}} + {k_1}x = \frac{{{2}{T_{{\rm{Load}}}}}}{{{D_{\rm{s}}}\tan (\varphi + \rho )}} $ | (32) |
式中,
$ {F_{{\rm{sf}}}}{\rm{ = (\text{π} }}ab{{\rm{\overline \tau }}_{\rm{c}}}{\rm{)}}{J_{\rm{s}}} = ({\mu _{\max }}N){J_{\rm{s}}} $ | (33) |
$ {J_{\rm{s}}}=\psi \frac{{{C_{\rm{s}}}}}{{\sqrt {{C_{\rm{c}}}^2 + {C_{\rm{s}}}^2} }} $ | (34) |
$\psi =\frac{{1}}{\text{π}}\left[\frac{\text{π}}{2} - \sin\left(\frac{{1 - {s^2}}}{{1 + {s^2}}}\right) + \frac{{2s}}{{1 + {s^2}}}\right]$ | (35) |
$ {s} = \frac{{\text{π}m}}{{4{\mu _{\max }}}}\sqrt {{C_{\rm{c}}}^2 + {C_{\rm{s}}}^2} $ | (36) |
式中,
在整体动力学研究中,应将该侧滑力
$ {C_{{\rm{sf}}}} = \frac{{4{F_{{\rm{sf}}}}}}{{\text{π}\sqrt {\dfrac{{{k_1}}}{{{m_2}}}X} }} $ | (37) |
一般来讲,处于稳定平衡运行状态时侧滑力较小;忽略实际结构中轴承等接触处的摩擦力,则上式中X可采用近似估算法,则有:
$ X=\frac{{{2}{T_{{\rm{Load}}}}}}{{{k_1}{D_{\rm{s}}}\tan (\varphi + \rho )}} $ | (38) |
式中,
于是,公式(32)进一步改写为:
$\begin{aligned}[b] {m_2} \cdot \frac{{{{\rm{d}}^2}x}}{{{\rm{d}}{t^2}}} &+ \frac{{{2}{D_{\rm{s}}}\tan (\varphi + \rho ){F_{{\rm{sf}}}}\sqrt {{k_1}{m_2}} }}{{\text{π}{T_{{\rm{out}}}}}}\frac{{{\rm{d}}x}}{{{\rm{d}}t}} + {k_1}x =\\ & \frac{{{2}{T_{{\rm{Load}}}}}}{{{D_{\rm{s}}}\tan (\varphi + \rho )}} \end{aligned} $ | (39) |
假设初始条件为:
$ \begin{aligned}[b] x =& \frac{{M(B + \sqrt {{B^2} - 4AC} )}}{{2C\sqrt {{B^2} - 4AC} }} \cdot {{\rm{e}}^{ - \frac{{B - \sqrt {{B^2} - 4AC} }}{{2A}}t}} +\\ & \frac{{M(B - \sqrt {{B^2} - 4AC} )}}{{2C\sqrt {{B^2} - 4AC} }} \cdot {{\rm{e}}^{ - \frac{{B + \sqrt {{B^2} - 4AC} }}{{2A}}t}} \end{aligned} $ | (40) |
其中,
$ \begin{aligned}[b] &A = {m_2},B = \frac{{{\rm{2}}{D_{\rm{s}}}\tan (\varphi + \rho ){F_{{\rm{sf}}}}\sqrt {{k_1}{m_2}} }}{{\text{π}{T_{{\rm{out}}}}}},C = {k_1},\\ &D = \frac{{{\rm{2}}{T_{{\rm{Load}}}}}}{{{D_{\rm{s}}}\tan (\varphi + \rho )}}{\text{。}} \end{aligned} $ |
对滚轮平盘牵引式CVT自适应变速机构进行动力学分析时,采用的参数如表1所示。根据设定的参数,可以求出CVT接触模型参数如下:接触椭圆长半轴a为1.9 mm,接触椭圆短半轴b为1.6 mm,最大赫兹接触压力
表1 滚轮平盘牵引式CVT参数设定 Tab. 1 Basic parameters of the roller flat disk traction-type CVT |
![]() |
3.2 计算结果与分析
在实际变速过程中,滑移率系数与侧滑率系数均会发生改变,其中滑移率系数变化较小,而侧滑率系数较大;侧滑率系数对于侧滑力的大小有明显影响,这也就直接影响到了自适应变速过程。
图4为根据Tevaarwerk和Johnson的牵引传动理论[10, 16]计算出的侧滑力与侧滑率系数之间的数据曲线。从图4中可知:当侧滑率系数低于0.10时,侧滑率系数对于侧滑力影响较大;当侧滑率系数高于0.55时,即使滑移率系数
![]() |
图4 侧滑力与侧滑率系数之间的关系曲线 Fig. 4 Relationship between the side-slip force and the side-slip rate |
力封闭弹簧的劲度系数
图5为自适应变速动力学响应曲线,图5(a)与(b)为输出负载由空载突变至有不同负载下,不同
![]() |
图5 自适应变速动力学响应曲线 Fig. 5 Dynamic response of adaptive variable speed changing mechanism |
此外,根据计算统计出
![]() |
图6 变速稳定时滚轮的偏移量 Fig. 6 Offset of the roller with k1 when the speed is stable |
![]() |
图7 变速稳定时滚轮的偏移减小量 Fig. 7 Offset reduction of the roller with k1 when the speed is stable |
从图(6)、(7)中可见,当
图8为
![]() |
图8 变速稳定时滚轮的偏移量与TLoad关系 Fig. 8 Relationship between the offset of the roller and the TLoad when the speed is stable |
图9、10为
![]() |
图9 k1对传动比变化的影响 Fig. 9 Influence of k1 on the change of transmission ratio |
![]() |
图10 TLoad对传动比变化的影响 Fig. 10 Influence of TLoad on the change of transmission ratio |
结合动力学响应与滚轮平盘牵引式CVT的结构特征的限制,可得其最大负载输出扭矩为20 Nm,优选
为了验证计算结果的正确性,利用Adams平台进行
![]() |
图11 仿真模型 Fig. 11 Simulation model |
![]() |
图12 TLoad=10 N·m时的动力学仿真结果 Fig. 12 Dynamic simulation result at TLoad=10 N·m |
对比图12中的仿真结果与图5(b)的结果,可以看出:自适应变速响应趋势吻合;若对其进行坐标统一,变速稳定状态时滚轮的偏移量也非常接近,仿真结果中滚轮偏移为负值,是因为仿真模型中坐标系原点位于平盘的中心轴上。而自适应变速稳定时间有明显区别,仿真模型中滚轮侧滑稳定时间较图5(b)中减少,这与仿真模型中采用了库伦简化模型处理接触处侧滑力有一定关系;而在实际情况下,因为侧滑力和阻尼的存在,微量振动一般很难发生,因此仿真模型更接近实际情况。总之,从对比来看,仿真结果可以有效验证计算结果的准确性。
4 结 论针对滚轮平盘牵引式CVT进行了自适应负载变速机构设计,并建立了相应的计算模型,包括运动学计算模型、接触计算模型和动力学计算模型。在此基础上,根据滚轮平盘牵引式CVT结构特征与设定参数进行了动力学研究与仿真实验,得出以下结论:
1)变速过程中,侧滑率对于侧滑力的影响规律为:侧滑率系数较低时,侧滑力随着侧滑率系数的增长急剧加大;而当侧滑率系数较高时,侧滑力大小基本恒定。
2)自适应负载变速过程是一个振动过程,
3)负载变化
4)初始负载条件不同时,相同输出负载仅影响到变速振动的振幅,但对于变速稳定状态时滚轮的偏移量无较大影响。
5)选择
[1] |
Fuchs R,Hasuda Y,Rothernbuehler Y,et al. Control concepts of continuously variable transmissions (CVT)[J]. JTEKT Engineering Journal, 2006, 1001(E): 24-29. |
[2] |
Shen Chujing,Yuan Shihua,Wei Chao,et al. Transmission efficiency of hydraulic clamp type traction transmission device[J]. Journal of Xi’an Jiaotong University, 2012, 46(7): 32-37. [沈楚敬,苑士华,魏超,等. 液压压紧式牵引传动装置传动效率研究[J]. 西安交通大学学报, 2012, 46(7): 32-37.] |
[3] |
Carbone G,Mangialardi L,Mantriota G. A comparison of the performances of full and half toroidal traction drives[J]. Mechanism and Machine Theory, 2004, 39(9): 921-942. DOI:10.1016/j.mechmachtheory.2004.04.003 |
[4] |
Li Chao,Li Hua,Li Qingtao,et al. Modeling, kinematics and traction performance of no-spin mechanism based on roller-disk type of traction drive continuously variable transmission[J]. Mechanism and Machine Theory, 2019, 133: 278-294. |
[5] |
Liu Changqi,Li Maosheng. Continuously variable transmission fluids technological market analysis in China[J]. Lubrication Engineering, 2009, 34(4): 105-107. [刘长期,李茂生. 国内无级变速器油技术市场分析[J]. 润滑与密封, 2009, 34(4): 105-107. DOI:10.3969/j.issn.0254-0150.2009.04.029] |
[6] |
Cretu O S,Glovnea R P. Constant power continuously variable transmission (CP-CVT):Operating principle and analysis[J]. Journal of Mechanical Design, 2005, 127(1): 114-119. DOI:10.1115/1.1828457 |
[7] |
Tyreas G C Nikolakopoulos P. G. Development and friction estimation of the half-toroidal continuously variable transmission a wind generator application[J]. Simulation Modelling Practice and Theory, 2016(66): 63-80. |
[8] |
周有强.机械无级变速器[M].北京:机械工业出版社,2001.
|
[9] |
Li Qingtao,Liao Min,Wang Shuang. A zero-spin design methodology for transmission components generatrix in traction drive CVTs[J]. Journal of Mechanical Design, 2017, 140(3): 33301. |
[10] |
Zhang Y,Zhang X,Tobler W. A systematic model for the analysis of contact,side slip and traction of toroidal drives[J]. Journal of Mechanical Design, 2000, 122(4): 523-528. DOI:10.1115/1.1290248 |
[11] |
Li Hua,Li Qingtao,Yao Jin. Analysis for contact stress of half toroidal traction continuously variable transmission[J]. Journal of Sichuan University (Engineering Science Editon), 2014(2): 166-170. [李华, 李青涛, 姚进. 半环面牵引式无级变速器接触应力分析[J]. 四川大学学报(工程科学版), 2014(2): 166-170.] |
[12] |
Huerta I M,Mohammadpour M,Rahnejat H.Analysis of lubricated contact in continuously variable transmissions (CVT)[C]//The 3rd Biennial International Conference on Powertrain Modelling and Control (PMC 2016).Loughborough:PMC,2016.
|
[13] |
Hamrock B J,Dowson D. Isothermal elastohydrodynamic lubrication of point contact,Part Ⅰ:Theoretical formation[J]. Journal of Lubrication Technology, 1976, 98(2): 223-228. DOI:10.1115/1.3452801 |
[14] |
Hamrock B J,Dowson D. Isothermal elastohydrodynamic lubrication of point contact,Part Ⅱ:Ellipticity parameter results[J]. Journal of Lubrication Technology, 1976, 98(3): 375-381. DOI:10.1115/1.3452861 |
[15] |
Brewe D E,Hamrock B J. Simplified solution for elliptical-contact deformation between two elastic solids[J]. Journal of Lubrication Technology, 1977, 99(4): 485-487. DOI:10.1115/1.3453245 |
[16] |
Tevaarwerk J L,Johnson K L. The influence of fluid rheology on the performance of traction drives[J]. Journal of Lubrication Technology, 1979, 101(3): 266-273. DOI:10.1115/1.3453346 |