工程科学与技术   2021, Vol. 53 Issue (2): 187-194
滚轮平盘牵引式无级变速器自适应变速机构动力学分析
李超1, 姚进2, 李华2, 胡波1     
1. 成都理工大学 核技术与自动化学院,四川 成都 610059;
2. 四川大学 机械工程学院,四川 成都 610065
基金项目: 四川省科技厅应用基础研究计划项目(2012JY0085)
摘要: 无级变速器(CVT)是一种理想的机械传动方案,随着现代工业和车辆技术的发展,希望CVT的变速机构对于负载改变能够实现自动/自适应变速,针对典型的同轴移动类牵引式CVT—滚轮平盘牵引式CVT,提出了一种基于螺纹–弹簧力封闭的自适应变速机构,并对其进行动力学研究。首先,给出具有螺纹–弹簧式自适应变速机构的滚轮平盘牵引式CVT结构模型;然后,分别对该CVT进行了运动学分析并给出运动学参数计算模型,根据赫兹接触理论建立其牵引点接触模型,采用分离体分析法建立了其动力学模型;最后,重点分析了变速过程中滚轮与平盘之间侧滑力的变化规律,分析出滚轮自适应负载变化的变速运动规律,研究了弹簧刚度系数对变速过程的影响,并对自适应变速过程的滚轮偏移规律进行了仿真实验验证。结果表明:自适应变速过程中,当侧滑率系数大于0.55时,侧滑力大小基本恒定;自适应变速过程为一振动过程,且变速机构中弹簧的刚度系数和负载力矩会直接影响变速响应时间和变速稳定时滚轮侧滑量大小;在同样大小的负载作用下,不同的初始负载条件会影响到变速振动的振幅,但对于变速稳定时滚轮的偏移量无较大影响;自适应变速过程的滚轮偏移规律的仿真结果与计算结果趋势吻合,但变速稳定时间有明显区别,仿真结果更接近实际情况。
关键词: 滚轮平盘    牵引式无级变速器    自适应变速    动力学分析    侧滑力    
Dynamic Analysis of Adaptive Shifting Mechanism of Roller Flat Disk Traction-type Continuously Variable Transmission
LI Chao1, YAO Jin2, LI Hua2, HU Bo1     
1. College of Nuclear Technol. and Automation Eng., Chengdu Univ. of Technol., Chengdu 610059, China;
2. School of Mechanical Eng., Sichuan Univ., Chengdu 610065, China
Abstract: The continuously variable transmission (CVT) is an ideal solution of mechanical transmission. As the development of modern industry and vehicle technology, the CVT is expected to realize the adaptive shifting when the load changes. For the roller flat disk CVT, which is a kind of traction CVT that the output shaft is coaxial with the shifting direction of the shifting component, a screw-spring adaptive shifting mechanism (ASM) of load torque was proposed and the dynamic analysis was conducted. Firstly, the structure model of roller flat disk CVT with screw-spring ASM was described. Then, the kinematic model was established by using analytical method, the contact model was established with Hertz contact theory, and the dynamic model was conducted through the separation analysis method, respectively. Finally, the change regular of side-slip force between roller and flat disk, as well as the transformation law of roller during adaptive shifting of different load torque were mainly analyzed, and the effect of balance spring stiffness coefficient on shifting process was also studied. Simulation experiment was also carried out to verify the transformation law of roller during the adaptive shifting process. The results showed that: In the process of adaptive shifting, the side-slip force is basically constant when the side-slip coefficient is higher than 0.55; The adaptive shifting process is actually a vibration process, and the stiffness coefficient and the load variety could directly affect the dynamic response characteristics; When the initial load conditions are different, the same output load only affects the amplitude of the vibration, but there is no significant effect on the offset of the roller in the speed stabilization state. The simulation results of the transformation law during adaptive shifting process are consistent with the trend of the calculated results, but there is a clear difference between the stable time. The simulation results are closer to the actual situation.
Key words: roller flat disk    traction continuously variable transmission    adaptive shifting    dynamic analysis    side slip force    

无级变速器(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中各参数说明如下:力封闭弹簧选用刚度系数 ${k_1}$ 的弹簧;输入轴扭矩 ${T_{{\rm{in}}}}$ ,输入转速 ${\omega _{{\rm{in}}}}$ $P(x,y)$ 为稳定状态是滚轮与平盘的接触点;平盘直径 ${D_1}$ ,工作半径 ${r_1}$ ;滚轮直径 ${D_2}$ ,工作半径 ${r_2}$ ;滚轮边缘圆弧半径 ${r_r}$ ,输出轴扭矩 $T_{\rm{out}} $ ,输出转速 $\omega _{\rm{out}} $

2 自适应变速机构动力学建模 2.1 运动学模型

1)理论传动比

假设接触点 $P(x,y)$ 为平盘径向外缘时,力封闭弹簧压缩量 $x = 0$ ,则该CVT的传动比可表示为:

$i = \frac{{{\omega _1}}}{{{\omega _2}}}=\frac{{{r_2}}}{{{r_1}}} = \frac{{{R_2}}}{{{R_1} - x}}$ (1)

式中, $x$ 为力封闭弹簧压缩量, $\omega_1 $ 为平盘角速度, $\omega _2 $ 为滚轮角速度。

2)自旋角速度

图2自旋角速度示意图。平盘角速度为 ${\omega _{1}}$ ,滚轮角速度为 ${\omega _{2}}$ ,输出轴角速度为 ${\omega _{4}}$ ,当处于稳定平衡状态时: ${\omega _{2}}={\omega _{4}}$

图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]。而这个微小的速度差产生的本质原因在于牵引传动接触点传递扭矩的实质是弹流润滑油膜的剪切作用。

定义滑移率[3-4]

$ {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所示,平盘与滚轮之间的加载力由加载凸轮机构提供,加载凸轮通过中间球体与端面接触产生加载力。假设 ${N_{\rm{c}}}$ 为端面凸轮与球体之间的压力, ${F_{{\rm{tr}}}}$ 为牵引力, ${F_{{\rm{ca}}}}$ 为输入轴端由凸轮产生的轴向力, ${F_{{\rm{ct}}}}$ 为输入轴端由凸轮产生的切向力, ${R_{\rm{c}}}$ 为凸轮升程段对应的中径, $\lambda $ 为凸轮升程段对应的倾斜角, ${L_{\rm{c}}}$ 为凸轮升程。

根据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)

式中, ${\;\mu _{\rm{c}}}$ 为凸轮和球体接触摩擦系数。

由式(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)

如果忽略摩擦,即 ${\;\mu _{\rm{c}}} = 0$ ,则轴向力 ${F_{{\rm{ca}}}}$ 可简化为:

${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)

其中 ${F_{{\rm{tr}}}} = \dfrac{{{T_{{\rm{in}}}}}}{{{R_1} - x}}$ ,则可得:

$ {\mu _{\rm{e}}} = \frac{{{L_{\rm{c}}}}}{{2\text{π}({R_1}{\rm{ - }}x{\rm{)}}}} $ (14)
2.2 接触模型

根据赫兹接触理论,当平盘与滚轮之间有法向加载力,则接触点处应为椭圆[11-12]。椭圆的方向和尺寸与接触点的法向力、平盘与滚轮的材料特性和结构尺寸均相关,其接触应力在接触区内按照椭球体规律分布,表达如下:

$ p={p}_{\rm{H}}\left(1-\frac{{x}^{2}}{{b}^{2}}-\frac{{y}^{2}}{{a}^{2}}\right) $ (15)

最大赫兹接触应力 ${p_{\rm{H}}}$ 为:

${p_{\rm{H}}} = \frac{{3N}}{{2\text{π}ab}}$ (16)

式中, $a$ 为接触椭圆的长轴, $b$ 为接触椭圆的短轴,N为法向载荷。

在Hamrock等[13-14]对于椭圆接触问题等温弹流润滑的研究中,采用了Harris提出的迭代算法,该方法可以较高精度地确定ab的值,如下:

$ 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^{'}}}$ 为当量弹性模量,可由以下公式求得:

$ {{E^{'}}}=\frac{{{\rm{1 - }}{\nu _{1}}^{2}}}{{{E_1}}} + \frac{{{\rm{1 - }}{\nu _{2}}^{2}}}{{{E_2}}} $ (20)

式中, ${\nu _{1}}$ 为平盘材料泊松比, $ {\nu }_{2} $ 为滚轮材料泊松比, $ {E}_{1} $ 为平盘材料弹性模量, ${E_2}$ 为滚轮材料弹性模量。

根据Brewe等[15]从迭代方法中采用数据处理回归出的近似公式来分别计算 $k$ $\varepsilon $ 值,可得:

$ {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)
2.3 动力学模型

将各部件分离,进行受力分析,为了简化方程数量,将固定连接部件视为整体部件,如图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)

式中, ${F_{{\rm{ib}}}}$ 为输入端轴承推力, ${F_{{\rm{ca}}}}$ 为端面凸轮轴向推力, ${T_{{\rm{in}}}}$ 为输入扭矩, ${F_{{\rm{ct}}}}$ 为端面凸轮周向推力, ${R_{\rm{c}}}$ 为端面凸轮作用半径, ${J_{{\rm{in}}}}$ 为输入端转动惯量, $\omega $ 为输入角速度。

图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)

式中, $N$ 为接触点法向压力, ${m_1}$ 为平盘组件质量, ${F_{{\rm{tr}}}}$ 为接触点牵引力, ${r_1}$ 为平盘接触半径, ${J_1}$ 为平盘组件转动惯量, ${\omega _{1}}$ 为平盘组件角速度。

图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)

式中, ${F_{{\rm{sa}}}}$ 为螺旋传动轴向推力, ${F_{{\rm{sf}}}}$ 为侧滑阻力, ${F_{{{\rm{k}}_1}}}$ 为弹簧弹力, ${m_2}$ 为滚轮质量, ${r_2}$ 为滚轮接触半径, ${T_{\rm{s}}}$ 为螺旋传动扭矩, ${J_2}$ 为滚轮转动惯量, ${\omega _2}$ 为滚轮角速度。

图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{ob}}}}$ 为输出端轴承推力, ${T_{{\rm{Load}}}}$ 为负载扭矩, ${J_{{\rm{out}}}}$ 为输出轴转动惯量, ${\omega _{{\rm{out}}}}$ 为输出角速度。

根据螺纹传动原理,有:

${F_{{\rm{sa}}}}=\frac{{{2}{T_{{\rm{Load}}}}}}{{{D_{\rm{s}}}\tan (\varphi + \rho )}}$ (31)

${F_{{{\rm{k}}_1}}}={k_1}x$ 和式(31)代入式(27),得:

$ {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}}}}$ 为变速时滚轮在平盘上径向移动时侧滑阻力,根据Tevaarwerk和Johnson的牵引传动理论,变速过程中忽略自旋影响时,侧滑力可以由以下公式求得[10,16]

$ {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{s}}}$ 为侧滑率系数(侧滑方向), ${C_{\rm{c}}}$ 为滑移率系数(牵引传动方向)。

在整体动力学研究中,应将该侧滑力 ${F_{{\rm{sf}}}}$ 等效为黏性阻尼,则等效阻尼系数 ${C_{{\rm{sf}}}}$ 计算如下:

$ {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)

式中, ${D_{\rm{s}}}$ 为螺旋传动的中径, $\varphi $ 为螺旋传动的螺旋升角, $\;\rho $ 为螺旋传动的摩擦角。

于是,公式(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)

假设初始条件为: $x(0) = 0,x'(0) = 0$ ,则得出:

$ \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} $
3 计算结果与分析 3.1 参数设定

对滚轮平盘牵引式CVT自适应变速机构进行动力学分析时,采用的参数如表1所示。根据设定的参数,可以求出CVT接触模型参数如下:接触椭圆长半轴a为1.9 mm,接触椭圆短半轴b为1.6 mm,最大赫兹接触压力 ${p_{\rm{H}}}$ 为1.548 kN/mm2

表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时,即使滑移率系数 ${C_{\rm{s}}}$ 在0.015~0.10之间变化,侧滑力的大小也趋于恒定值;而变速过程中,滑移率系数 ${C_{\rm{s}}}$ 在0.015左右,侧滑率系数接近于1,则侧滑力大小将保持在420 N左右。

图4 侧滑力与侧滑率系数之间的关系曲线 Fig. 4 Relationship between the side-slip force and the side-slip rate

力封闭弹簧的劲度系数 $ {k}_{1} $ 是自适应变速机构设计的重要参数,其大小对于自适应变速的动力学响应有着重要影响。一般选取原则是能够快速稳定响应自适应扭矩变速的前提下取较小的值以减小能耗。

图5为自适应变速动力学响应曲线,图5(a)与(b)为输出负载由空载突变至有不同负载下,不同 ${k_1}$ 对于动力学响应的影响曲线,而图5(c)为初始负载条件不同时的动力学响应。从图5中可知:输出负载 ${T_{{\rm{Load}}}} = 5\;{\rm{ N}} \cdot {\rm{m}}$ 时,自适应变速响应在15 ms之内基本稳定,而 ${T_{{\rm{Load}}}}= 10\;{\rm{ N}} \cdot {\rm{m}}$ 时,对应的响应稳定时间接近30 ms;而同一负载增量下, ${k_1}$ 越大,响应稳定越快,侧滑稳定时滚轮移动量越小;初始负载条件不同时,仅影响到变速振动的振幅,但对于变速稳定状态时滚轮的偏移量无较大影响。

图5 自适应变速动力学响应曲线 Fig. 5 Dynamic response of adaptive variable speed changing mechanism

此外,根据计算统计出 ${k_1}$ 每增加50 N/mm下动力响应稳定时 $x$ 的偏移量如图6所示,偏移量的减小量如图7所示。

图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)中可见,当 ${k_1}$ 越大时,偏移量越小,且偏移减小也越不明显。

图8 ${T_{{\rm{Load}}}}$ 变化与变速偏移量的关系。从图8可知, ${T_{{\rm{Load}}}}$ 改变量与变速稳定时滚轮偏移量之间近似为线性递增关系,其线性斜率随弹簧劲度系数 ${k_1}$ 的增加而减小。

图8 变速稳定时滚轮的偏移量与TLoad关系 Fig. 8 Relationship between the offset of the roller and the TLoad when the speed is stable

图910 ${k_1}$ ${T_{{\rm{Load}}}}$ 对传动比的影响。从图910中可知: ${k_1}$ 对传动比的影响呈负相关性,且输出负载 ${T_{{\rm{Load}}}}$ 变化越大,影响越明显;而 ${T_{{\rm{Load}}}}$ 对传动比的影响呈正相关性,且 $ {k}_{1} $ 越小,影响越明显。

图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,优选 ${k_1} = 750$ N/mm;此条件下,最大变速稳定下滚轮偏移量的理论值为 ${x_{\max }} = 22.82\;{\rm{mm}}$ ,传动比范围为0.8~1.47之间。

3.3 仿真实验

为了验证计算结果的正确性,利用Adams平台进行 ${T_{{\rm{Load}}}}= 10\;{\rm{ N}} \cdot {\rm{m}}$ 时系统的动力学仿真实验,如图11所示为建立的仿真模型。其中:简化了输出轴和滚轮之间的螺纹连接,采用等效恒力施加在滚轮上;简化滚轮与平盘接触处的侧滑力计算,采用库仑摩擦模型等效;简化加载凸轮机构,采用恒力加载方式代替。在仿真条件中,取与计算相同的参数,仿真滚轮的偏移量,结果如图12所示。

图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)自适应负载变速过程是一个振动过程, $ {k}_{1} $ 越大,动力响应稳定越快,侧滑稳定时滚轮移动量x越小,传动比越小;而负载变化 ${T_{{\rm{Load}}}}$ 越大,动力响应达到稳定状态的时间就越长,侧滑稳定时滚轮移动量x越大,传动比越大。

3)负载变化 ${T_{{\rm{Load}}}}$ 与侧滑稳定时滚轮移动量x之间可采用线性关系表示,线性斜率与 ${k_1}$ 相关。

4)初始负载条件不同时,相同输出负载仅影响到变速振动的振幅,但对于变速稳定状态时滚轮的偏移量无较大影响。

5)选择 ${k_1}$ =750 N/mm,由于滚轮平盘牵引式CVT结构限制,设定参数下最大输出扭矩20 N·m,变速稳定下的滚轮最大偏移量的理论值为 ${x_{\max }} = 22.82\;{\rm{mm}}$ ,传动比范围为0.80~1.47之间。

参考文献
[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