非平稳随机循环工况离合器接合优化及跟踪控制

李晓祥, 王安麟, 樊旭灿, 吴景培

(同济大学 机械与能源工程学院, 上海 201804)

离合器作为车辆传动系统的重要部件，在各类型的变速器上得到了广泛的应用[1-4]. 对离合器接合过程进行品质控制可降低车辆冲击度，减少接合滑磨功，减轻对传动零部件冲击，使车辆使用寿命和驾乘舒适性得到提高[5-7]. 国内外对离合器接合过程品质有不同的评价参数，研究多集中在接合过程轨迹最优和执行机构跟踪效果两个方面[8-9]. Haj-Fraj、彭建鑫、胡宏伟等[10-12]分别利用动态规划法、遗传算法和接合曲线修正算法等多种方法实现离合器接合过程最优控制；陈清洪等[13]以滑磨功最小为目标、冲击度为约束进行接合过程的最优控制；Wang、单外平等[14-15]通过试验与数学模型相结合方法模拟离合器接合轨迹对接合过程品质的影响；离合器接合轨迹研究过程中，从动盘端多加载定值信号，较难反映实际载荷的非平稳性. 因离合器执行机构及离合器自身具有较强非线性特性，多采用具有鲁棒性的控制器实现执行机构跟踪[16-18]. 滑模控制方法算法简单，能克服系统的不确定性, 对外界噪声干扰和未建模动态参数摄动具有鲁棒性，对非线性系统具有良好的控制效果[19-22]. Li Junqiu、赵韩、孙承顺等[23-25]采用滑模控制方法实现离合器接合过程中车速、油压及其位移的鲁棒跟踪；Li、赵治国等[26-27]采用滑模控制和不确定观测器相结合的方法实现离合器接合过程扭矩估计和油压轨迹跟踪. 考虑系统内外部扰动等因素，以冲击度和滑磨功为约束条件的离合器接合轨迹及跟踪相结合控制的研究较少.

本文在变速器换挡规律确定的前提下，以装载机V型重载后退换挡中离合器接合过程为研究对象，通过建立冲击度和滑磨功相结合的二次型性能指标泛函和依据Pontryagin极大值原理等方法，实现其时变控制律下的可变扭矩最优轨迹，并采用滑模方法进行跟踪控制设计. 通过对离合器冲击度、滑磨功和跟踪精度等数据分析，验证了设计方法的有效性，改善了离合器接合过程品质.

1 离合器接合过程数学模型

在换挡规律确定前提下，离合器接合过程的研究主要集中在主从动摩擦片开始接触直至两者转速差为零的滑磨时间段，该过程是车辆冲击度和滑磨功产生的阶段[7]. 离合器接合过程系统结构原理如图1所示.

①离合器主动端(动力端)模型； ②接合过程动力学模型； ③执行机构模型； ④压力控制模型； ⑤离合器从动端(载荷端)

图1 接合过程系统结构原理简图

Fig.1 Schematic diagram of structural principle of engagement process system

1.1 离合器主从动端模型参数化构建

以装载机V型重载后退换挡中变速箱输入输出端实测载荷和转速数据为基础进行参数化模型构建，有利于后续台架或车辆试验前的仿真和加载. 部分元器件安装和测试采集的数据曲线如图2所示.

(a)部分元器件安装图 (b)测试数据曲线图2 部分元器件安装图和测试曲线Fig.2 Part component installation diagram and test curves

对0.8～1.5 s内多组变速器输出端扭矩随机试验样本数据的平稳性进行轮次数检验，可知该时间段信号非平稳，其重构数学模型为[28]

xt=dt+st.

(1)

式中：xt为某一时间历程非平稳随机载荷，dt为xt中的确定性分量，st为xt中平稳随机分量.

利用变速器输出扭矩所构建的参数化数学模型均值数据与文献[6,10]等所得数据误差在1%以内，说明该模型在一定程度上可用于系统仿真和试验加载. 变速器输入端扭矩和转速数据去除确定性分量后，两者符合均值为零正态分布，亦可用式(1)进行参数化模型表达.

1.2 接合过程动力学模型

为更好地研究离合器接合过程的动态特性，对模型进行简化并作如下假设：1)离合器接合滑磨过程采用冷却油冷却，忽略摩擦片热变形；2)忽略零部件间隙和轴横向振动的影响；3)忽略轴承、轴承座和齿轮啮合的弹性变形；4)忽略系统中弹性元件的惯性和惯性元件的弹性. 根据上述假设，离合器接合过程滑磨阶段模型如图3所示.

动力学方程为

(2)

式中：Jp为离合器主动盘之前(包括主动盘)的动力系折算到主动盘的等效转动惯量，Tp为动力系折算到离合器主动盘的等效动力矩，ωp为离合器主动盘转速，Tcl为离合器主从动盘传递的扭矩，Jl为离合器从动盘之后(包括从动盘)的传动系折算到从动盘的等效转动惯量，Tl为负载端扭矩折算到从动盘的等效阻力矩，ωl为离合器从动盘转速，cp和cl为等效的旋转黏性阻尼系数.

图3 接合过程模型Fig.3 Engagement process model

1.3 执行机构模型

此阶段离合器活塞停止运动，钢片与摩擦片间隙消除，作如下假设：1)离合器腔压力油温度和黏度不变；2)离合器腔无泄漏；3)忽略摩擦片的弹性变形. 离合器执行机构模型如图4所示.

图4 离合器执行机构模型Fig.4 Clutch actuator model

离合器受力平衡方程为

Fs=Fp+Fd-Fcl-Fg-Fseal.

(3)

式中：Fs为离合器复位力；Fp和Fd分别为液压油作用于离合器活塞的静压力和动压力，Fd=ρAp2(αωp)2/(4π)；Fcl为离合器摩擦片对活塞反作用力，Fcl=Tcl/(μ(t)ZR)；Fg为齿圈花键摩擦力；Fseal为离合器活塞密封阻力；ρ为液压油密度；α为液压油旋转速度滞后系数；Ap为离合器活塞面积；μ(t)为摩擦副动摩擦因数；Z为摩擦副数量；R为摩擦力作用等效半径.

离合器腔流量压力动态方程为

(4)

式中：Qc为流入离合器腔流量，xmax为活塞运动最大位移，V0为离合器腔和油道初始总体积，E为流体体积弹性模量，pc为离合器腔控制压力.

1.4 压力控制模型

压力控制通过PWM控制的高速电磁开关阀来实现，在忽略油道液阻力、管路和阀自身泄漏基础上，进入离合器腔控制流量方程为

(5)

由以上数学模型知，离合器接合过程的系统之间有较强的非线性时变特性，可采用鲁棒性较强的非线性控制理论来实现其控制. 选取的典型工况下离合器接合过程控制设计方法见图5.

图5 接合过程控制设计流程Fig.5 Engagement process control design

2 离合器接合轨迹优化

2.1 状态空间方程模型

将滑磨功和冲击度作为离合器接合过程品质的评价指标[7]，其数学模型为

(6)

式中：W为离合器接合过程滑磨功，t0为接合开始时刻，tf为主从动盘转速差为零时刻，j为冲击度，i0为变速器至驱动轮传动比，rw为车轮半径.

考虑滑磨功和冲击度式(6)并结合式(2)，取状态向量x=(ωp,ωp-ωl,Tcl)T，控制变量为u，建立的状态方程为

(7)

式中：

2.2 接合轨迹最优

根据离合器接合过程的特点，本文最优控制问题可归结为终端时刻tf有限，终端状态x2(tf)有约束的最优控制问题[29].

二次型性能指标泛函如式(8)所示，其中对冲击度的评价定义为接合过程冲击度平方总和.

(8)

式中：

R为控制加权矩阵，β为权重系数，γ为系数(i0Jl/rw)2.

构造的Hamiltonian函数为

H(x,u,λ)=0.5(xTQx+βu2)+λ(t)T(Ax(t)+Bu(t)+D(t)).

式中λ(t)为待定的拉格朗日乘子向量函数.

初始条件和终端约束条件：

x(t0)=[ωp0,ωp0-ωl0,Tcl0]T,N[x(tf),tf]=x2(tf)=0.

由极值条件，得

(9)

同时由于

∂2H/∂u2=β>0，

根据Pontryagin极大值原理，可知式(9)中确定的u*(t)是泛函取得极小值的最优控制.

同时，规范方程为

横截条件为

λ(tf)=[0,ν,0]T.

(10)

式中ν为待求值λ2(tf).

干扰矩阵D(t)的存在，使最优控制u*(t)无法通过求解Riccati方程中P(t)来获得λ(t)[29].λ(t)通过如下方程组构建

(11)

式中：P(tf)为3阶零矩阵，G(tf)=[0,1,0]T用于求解v，H(tf)=[0,0,0]T用于补偿向量D(t)，F(tf)=[0,1,0]，m(tf)=0，n(tf)=0.

联立方程(9)、(10)和(11)及相关方程的微分，在通过Riccati微分方程求解P(t)基础上得G(t)、H(t)、F(t)、m(t)和n(t)，待定值v求解后得最优控制：

2.3 仿真及分析

在2.2节分析基础上，确定tf为1.1 s和权重系数β为1～100范围后，根据表1中参数对离合器接合过程进行仿真及分析. 在式(7)对u与Tcl的定义和Tcl≤Tclmax约束下，得离合器接合过程可变扭矩最优轨迹Tcl*(t).

表1 离合器接合过程相关参数Tab.1 Parameters of clutch engagement process

在冲击度小于我国标准17.64 m/s3下使J极小以确定β，图6为不同权重系数β下滑磨功和冲击度曲线. 图7为离合器主从动盘接合过程转速和可变扭矩最优轨迹曲线，得tf为1.09 s，β为1，冲击度最大值为15.05 m/s3，比标准低14.68%.

(a)滑磨功变化曲线 (b)冲击度变化曲线图6 权重系数β对评价指标的影响

Fig.6 Influence of weighting coefficientsβon evaluation index

(a)转速最优轨迹曲线 (b)扭矩最优轨迹曲线图7 最优轨迹曲线Fig.7 Optimal trajectory curves

在离合器摩擦片从动盘阻力矩幅值±5%Tl的变化集合下进一步分析最优轨迹的特性，由图8和表2可知，最优轨迹在数学模型增加误差幅度状态下的关键参数变化程度对其性能影响较小.

图8 ±5%Tl变化下最优轨迹曲线Fig.8 Optimal trajectory curves under ±5%Tl

表2 ±5%Tl下关键参数变化百分比Tab.2 Percentage changes of key parameters under ±5%Tl %

考虑变速器较多挡位数量和摩擦片用扭矩传感器可靠性(实际车辆中难以将其安装使用[5])，无法通过检测该扭矩信号实现闭环控制. 为实现离合器接合过程中最优轨迹的精准跟踪，采用对系统参数摄动和外干扰具有完全不变性的滑模控制方法.

3 基于接合轨迹的滑模控制

针对非线性特性执行机构参数摄动下最优轨迹的跟踪精度影响接合过程品质的问题，将难检测的可变扭矩转化为可通过压力传感器检测的离合器压力跟踪问题. 与采用车速反馈跟踪相比，该方法能实现离合器换挡过程中充油阶段的检测，有利于整个换挡过程品质的控制设计[30-32]. 通过式(3)中离合器受力平衡，可实现可变扭矩最优轨迹向压力最优轨迹的转化.

(12)

3.1 滑模控制相关理论

选取的切换面为

(13)

式中c为滑模面系数.

当切换面s=0时，其滑动模态方程为

在c>0保证系统在切换面稳定前提下，通过Hurwitz确定c范围[22].

选取的指数趋近律为

(14)

式中：ε为增益系数，k为趋近律，-εsign(s)项有在一定时间内全局到达切换面的作用，k项保证了趋近时间不过长.

3.2 滑模控制器设计

定义状态变量y2=pc(t)，控制输入为τ，联立式(4)和(5)，建立的跟踪系统状态方程为

式中y1为状态变量y2的积分值，

选择离合器压力最优轨迹pc*(t)作为参考输入r2，引入偏差变量跟踪误差e，建立微分方程为

(15)

联立式(13)和(14)，并将式(15)带入，得

基于Lyapunov稳定性理论，选取函数V=0.5 s2，进行稳定性分析. 在切换面s≠0时，所设计的控制τ能使控制系统状态在有限时间内收敛至切换面s=0，之后进入滑动模态后

为避免或削弱滑动模态中的抖振，采用饱和函数sat(s)代替符号函数sign(s).

3.3 仿真验证及分析

根据前述离合器接合过程动力学模型及滑模控制器，利用Matlab/Simulink建立离合器压力最优轨迹跟踪仿真模型. 取c=30,ε=2，k=2 000, 边界层厚度为0.001，设置不同初始点，仿真结果见图9. 其中(a)初始值为离合器位移xmax时压力跟踪曲线，(b)初始值为0，初始值的不同使得图(b)收敛时间延迟0.004 s，两者最大跟踪偏差分别在0.001%和0.280%以内，设计的控制器使离合器接合压力在t0～tf内有效地实现了对参考输入r2的精准跟踪.

(a)初始值xmax轨迹跟踪曲线 (b)初始值0轨迹跟踪曲线图9 初始值对轨迹跟踪影响Fig.9 Influence of initial values on trajectory tracking

图10为传统PID控制方法与滑模控制方法在图9(b)基础上的跟踪效果对比图. 由图10可知：所采用的滑膜控制方法在开始0.01 s内跟踪偏差比传统PID控制方法小40%～130%，该方法在0.01 s后跟踪偏差小于0.1%，总体跟踪效果相对更好.

(a)轨迹跟踪对比曲线 (b)跟踪误差对比曲线图10 PID和SMC跟踪效果对比Fig.10 Comparison of tracking effects between PID and SMC

为进一步分析具有非线性特性离合器执行机构的参数摄动对所设计控制器性能的影响，对式(12)中摩擦因数μ(t)±10%摄动变化，由图11可知，所设计的控制器能实现精准跟踪，最大跟踪偏差和收敛延迟时间分别在0.3%和0.05 s以内.

图11 ±10% μ(t)下轨迹跟踪曲线Fig.11 Trajectory tracking curves under ±10% μ(t)

4 结 论

1)基于最优化理论的非平稳随机循环工况下，离合器接合轨迹在一定程度上能有效地降低冲击度和滑磨功，利于提高车辆换挡接合品质和驾乘舒适性.

2)采用指数趋近律滑模控制方法实现离合器腔压力偏差0.3%内精准鲁棒跟踪，并克服了系统不确定性及外界摄动影响.

3)非平稳随机循环工况下离合器接合轨迹优化及跟踪控制的方法为解决某土方机械变速器离合器接合压力跟踪偏差过大导致的摩擦片烧结故障问题提供了理论性的工程借鉴.