自动同步离合器恒阻尼控制方法

2020-07-27蒋德松王荣杰王亦春周春良

哈尔滨工程大学学报 2020年4期

蒋德松，王荣杰，王亦春，周春良

(1．集美大学轮机工程学院,福建厦门 361021；2．福建省船舶与海洋工程重点实验室,福建厦门 361021；3．哈尔滨工程大学动力与能源工程学院，黑龙江哈尔滨 150001)

双原动机动力装置是由2种型式或不同功率的推进装置组成，通常情况下，在低速工况下工作时，由其中一台原动机单独工作；而在高速直至全工况时，由另一台高功率的原动机单独工作或两个原动机共同在线工作。双原动机驱动负荷的动力装置比较适宜于负荷变化范围大的动力系统。如大型工程机械，舰船推进等，也可以用于移动应急电站，以提高供电的可靠性[1-2]。自动同步离合器作为双原动机动力装置中的常用的传动部件,其性能直接影响着系统的整体性能，自动同步离合器在接入过程前期，要求阻尼油腔无阻尼力或者产生较小的阻尼力，这样有利于提高切换过程的机动性能。在自动同步离合器的接入过程中，中间传动部件处于零阻尼状态,角加速度为正值，中间传动部件与动力输出部件之间就会产生过大的转速差进而在二者之间产生剧烈的结合碰撞，在接入运动的后期，即后5 mm，液压阻尼油腔产生较大的阻尼力，有效的缓解了中间传动部件与动力输出部件之间的结合碰撞力，从而避免损坏离合器；另一方面，对处于脱离运动状态的自动同步离合器来说，为防止很小的转速差就产生脱开运动，在脱开动作的前段，液压阻尼油腔产生的驻退阻尼可以防止离合器产生频繁的脱开和接入动作，现有的条件下，通过人工更换阻尼螺栓这个配件可以调整阻尼力的大小[3-6]。

本文在双原动机(以柴油机和燃气轮机双原动机为例)及传动部件的仿真平台上通过改变柴油机的控制规律(调速器的PID控制参数)，在恒阻尼条件下就改变工作模式切换过程的时间参数进行可行性分析，并结合实验结果进行比较研究。

1 原动机数学模型

增压中冷柴油机的数学模型由以下6个部分组成[7]：

1)柴油机本体模型为:

(1)

式中：q为供油量波动系数；φ为柴油机曲轴角速度波动系数；ρ为增压压力变化；αn为耗能参数；k为系统结构系数；T为加速时间；θ1为增压压力增益；θ2为负荷装置调整增益。

2)增压器数学模型为：

(2)

式中：ξ为无量纲的涡轮入口压力变化；θq和θρ为无量纲放大系数，表示循环供油量和增压压力的波动对涡轮增压器的影响参数；φK为无量纲的增压器转子角速度改变；TK为速度调整时间；KK为涡轮增压器的自平衡因子。

3)柴油机进气数学模型描述为：

kBρ=φK-θBφ

(3)

式中：kB为进气管自平衡无量纲因子；θB为发动机曲轴角速度无量纲增益。

4)柴油机排气单元数学模型为：

krξ=φ+θrρ-θGq

(4)

式中：kr为排气管自平衡无量纲系数；θr为增压压力放大系数；θG为循环供油量增益。

5)供油系统数学模型为：

kqq=x

(5)

式中：x为无量纲的曲轴角速度变化；kq为供油位置参数。

综合上述部件模型得到二阶线性非齐次微分方程：

(6)

式中系数由柴油机及增压器的静态特性来确定。在Matlab/Simulink环境中得到柴油机的数学模型[8-9]。模型表达了柴油机及增压器各部件仿真模块的连接以及相互之间的参数联系。

6)调速器数学模型描述为：转速测量值和设定值形成的转速差信号作为比例积分微分调节器的输入参数，经功率放大模块驱动执行单元为：

(7)

式中：Kp为增益；Ti为积分环节时间参数；Td为微分环节时间参数。

电磁驱动机构是执行单元主要组成部分，其作用是用来把控制器控制的输出信号线性变换为油门位置。其工作原理是驱动机构和弹簧弹力的力学平衡，执行器单元的数学模型经拉普拉斯变换后表示为：

(8)

式中b1、b2、b3为结构参数。

调速器的仿真模型如图1所示，模型中调速器根据转速差及负荷修正，通过PID算法来控制齿杆的位置，从而控制油门的大小。图中设定转速和初始转速由控制单元产生。

图1 电子调速器仿真模型Fig.1 Simulation model of electronic governor

由于联合动力装置实验台燃气轮机采用的是定转速调节系统，被调量只有1个，即燃油流量Gf。鉴于该燃机的特点，转速的调节范围小，在额定负荷下，变速齿轮箱输出转速为1 500 r/min。经多次实验得出燃机具有4%的静态下垂特性。通过对燃气轮机转速调节器的改造，使燃气轮机机可以在小范围内调速。对于调节后的燃机转速和功率的关系，可以认为是线性的，尤其在驱动模式转换过程中，燃机在较低负荷的运行状态更是如此。其静态下垂特性可表示为：

ng=ng0-60Ng/220.5

(9)

式中：ng、ng0分别为加负荷后转速和无负荷时转速；Ng为功率。

切换过程距离燃机运行的危险边界较远，因此考虑其功率、转速的动态变化过程，对其他如压气机出口温度和压力、涡轮排气温度等在仿真中不作计算。建立燃气轮机转速控制模型，如图2所示。

图2 燃气轮机转速控制模型Fig.2 Gas turbine speed control model

2 基于FLANN-PSO的传动装置非线性动态系统辨识

自动同步离合器仿真模型中，输入参数为主动部件转速，从动部件的转速以及主动部件的力矩，输出则为发动机转速信号和自动同步离合器状态参数。当主动部件和从动输出部件之间的转速差大于一定值时自动同步离合器接入，反之自动同步离合器脱开。离合器的完全接入或脱开之间存在一个滞后环节，同时，模型中对离合器的负载扭矩区间作了限值，设计的扭矩峰值为工作值的1.1倍，最小扭矩工作值的0.9倍，这2个参数可以调整[10-12]。为了简化辨识模型和减少运算工作量，本文基于函数链神经网络(functional link artificial neural networks，FLANN)和粒子群优化算法(particle swarm optimization，PSO)，构建一种非线性动态系统自适应辨识方法，辨识过程中，具有三角函数特征的神经网络用来逼近系统，而粒子群算法主要用于训练，获得与函数神经网络中各个三角函数相对应的权值，以达到系统辨识的目的。非线性动态装置自适应辨识模型如文献[13]所示[13]。

系统辨识的模型在(k+1)时刻的输出为：

(10)

式中：f和g分别表示函数神经网络结构的前馈和反馈部分，其权值由粒子群算法来获取，输入x(k)经非线性扩展为：

φ(k)=[φ0(k),φ1(k),φ2(k),…,φ2n-1(k),φ2n(k)]T

(11)

其中：

(12)

由此可看出，扩展后三角函数有n项，还有1项(首项)为单位输入，由此可得到输入向量含有(2n+

1)个分量。设式(11)和式(12)中在第k时刻输入向量的权值为：

W1(k)=[w10(k),w11(k),w12(k),…,

w1,2n-1(k),w1,2n(k)]T

(13)

式(10)中第1部分的输出估计值为：

(14)

式(10)中第2部分输出估计值也由式(12)得到：

(15)

其中：

W2(k)=[w20(k),w21(k),w22(k),…,

w2,2n-1(k),w2,2n(k)]T

(16)

φ(k)=[φ0(k),φ1(k),φ2(k),…,

φ2n-1(k),φ2n(k)]T

(17)

即：

(18)

图3所示模型的输出估计值为：

图3 非线性动态装置自适应模型Fig.3 Adaptive model of nonlinear dynamic device

(19)

粒子群算法是一种基于种群的随机优化技术[14]，在优化求解的时候，首先寻求空间内随机初始化粒子，每个粒子都有自己的位置和速度参数，同时含有与被优化函数相关的适应值。在迭代过程中通过参数匹配的程序来实现自我更新，并评价新的适应函数值，适应函数值可由下式计算得到：

(20)

3 仿真计算及实验结果分析

柴油机和燃气轮机携带各自的自动同步离合器通过传动系统与电涡流测功器相联接，为减少自动同步离合器的结合碰撞带来的实验风险，电涡流测功器把负荷设置在较低的区间。部分实验结果如图4所示，该实验结果是在200 N·m负荷不同阻尼孔0、2.5以及4 mm阻尼孔条件下获取的。

图4 动力切换实验在200 N·m负荷不同阻尼孔条件下的时间参数Fig.4 Time parameters of power switching experiment with different damper hole condition under load 200 N·m

图4中3组数据分别显示了200 N·m负荷不同阻尼孔条件下的时间参数，每组数据有2部分。上部分显示动力由燃气轮机到柴油机切换柴油机端自动同步离合器接入过程曲线；下部分显示动力由燃气轮机到柴油机切换燃气轮机端自动同步离合器脱开过程曲线。接入时间的试验结果分别为0.37、0.29和0.26 s。

自动同步离合器的配件中准备了3种阻尼螺栓即0、2.5以及4 mm这3种规格，仿真计算中模型采用0 mm的阻尼螺栓，在200 N·m负荷条件下调节柴油机的控制器即图1电子调速器模型的PID控制模块的参数(用调节方程描述，记为调节方程1,调节方程2,调节方程3)图5分别显示出了相同工况在3种不同调节方程控制下的切换时间参数中的接入时间。柴油机端分别为0.38、0.28和0.26 s(经过四舍五入得到)，与实验结果基本吻合。

图5 不同调节方程控制作用下的接入时间Fig.5 Meshing time controlled by different governing equations

从仿真结果可以看出，以负荷200 N·m负荷为例，为了免除了更换阻尼螺栓的工作，2.5和4 mm阻尼孔条件下的接入时间,可以通过改变调节方程在0 mm阻尼孔条件下获取。同时有理由推断在系统稳定以及传动系统强度能承受接入碰撞力的前提下，接入时间可以在相应的时间段内连续可控。大量的实验结果显示阻尼孔从0～4 mm在高负荷下,接入时间的改变非常微小。如图6所示分别显示了0、2.5以及4 mm阻尼孔条件下接入时间随负荷的改变情况。