本发明公开一种液力变矩器解锁和滑差工况扭转振动隔振率计算方法，以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程；求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。本发明提供的方法，能够预测液力变矩器在解锁和滑差状态的扭转振动隔振率，为控制发动机传递到变速器的扭转振动和变速器齿轮敲击声提供理论依据。
1.液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：包括如下步骤，S1：简化液力变矩器在解锁和滑差状态时的动力传动系统旋转部件，依据各旋转部件的设计参数，计算旋转部件的等效惯量、扭转刚度和阻尼；
S2：通过台架实验得到液力变矩器的特性曲线，依据特性曲线计算得到液力变矩器泵轮、涡轮的动态扭矩；
S3：通过台架实验得到锁止离合器摩擦片的摩擦系数，计算锁止离合器的摩擦扭矩；
S4：以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程；
其中，θ＝{θp,θc,θtb,θtm,θw,θv}为角位移矢量；
J＝diag(Jp,Jc,Jtb,Jtm,Jw,Jv) (8)其中Jp是发动机曲轴、驱动盘和液力变矩器泵轮的等效惯量；Jc是锁止离合器的等效惯量；Jtb是液力变矩器涡轮的等效惯量；Jtm是变速器齿轮的等效惯量；Jw是轮胎的等效惯量；
其中，Kc,Cc分别是锁止离合器的扭转刚度和阻尼；Kis,Cis分别是变速器输入轴的扭转刚度和阻尼；Kds,Cds分别是驱动轴的等效扭转刚度和阻尼；Kw,Cw分别是轮胎的等效扭转刚度和阻尼；
T＝{Te‑Tp‑Tf Tf Tt 0 0 ‑Tv} (11)Te为发动机扭矩，Tv为整车等效阻力矩，Tp为泵轮扭矩，Tt为涡轮扭矩，Tf为解锁和滑差状态的锁止离合器摩擦扭矩；
S5：求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。
2.如权利要求1所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：S5：通过台架实验采集得到发动机激励扭矩，通过实车道路实验采集得到车辆等效阻力矩，求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。
3.如权利要求2所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：S1：计算以下旋转部件的等效惯量：
1)发动机曲轴(1)、驱动盘(2)和液力变矩器泵轮(4)的等效惯量Jp；
2)变速器输入轴(6)的扭转刚度和阻尼Kis,Cis；
3)驱动轴(8)的等效扭转刚度和阻尼Kds,Cds；
4.如权利要求3所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：S2：通过台架实验得到液力变矩器的特性曲线，依据特性曲线计算得到液力变矩器泵轮、涡轮的动态扭矩；
解锁和滑差状态时液力变矩器泵轮、涡轮的动态扭矩特性可表达为：式中：Tp‑泵轮扭矩，Tt‑涡轮扭矩，wp‑泵轮角速度，wt‑涡轮角速度；a0,a1,a2,b0,b1,b2为待定系数，由液力变矩器的特性曲线可求得，而液力变矩器的特性曲线由单体台架实验得到，包括变矩比、泵轮扭矩容量系数和效率。
5.如权利要求4所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：S3中:通过台架实验得到锁止离合器摩擦片的摩擦系数，计算锁止离合器的摩擦扭矩；
基于Karnopp摩擦模型计算解锁和滑差状态的锁止离合器摩擦扭矩，计算表达式为：式中：Ns‑摩擦面数目，R‑摩擦片有效半径，Fn‑法向压力，‑相对滑摩速度，μ为动摩擦系数，可由下式求得：式中，μ0‑初始动摩擦系数，ms‑摩擦系数变化梯度，均由台架实验采集得到。
6.如权利要求5所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，其特征在于：S5中：通过台架实验采集得到发动机激励扭矩，通过实车道路实验采集得到车辆等效阻力矩；车辆加速过程的发动机扭矩Te可用下式表达：Te＝Tavg+T2rd×chirp(t,f1,t1,f2) (12)式中，Tavg‑发动机平均扭矩，T2rd‑发动机二阶波动扭矩，chirp(t,f1,t1,f2)表示从0到t1时间段内，扭矩波动频率由f1扫频到f2；发动机平均扭矩和发动机二阶波动扭矩均由发动机台架实验采集得到；
车辆加速过程的整车等效阻力矩Tv可用下式表达：式中，θv‑轮胎角速度，A,B,D为待定系数，通过实车道路实验确定。
7.一种可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1‑6任一项所述的液力变矩器解锁和滑差工况扭转振动隔振率计算方法的步骤。
液力变矩器解锁和滑差工况扭转振动隔振率计算方法和可读
存储介质
技术领域
[0001] 本发明属于汽车技术领域，具体涉及液力变矩器解锁和滑差工况扭转振动隔振率计算方法。
背景技术
[0002] 随着汽车产业的不断发展，人们对车内驾乘舒适性要求越来越高，而车内噪声振动是舒适性的一个重要指标。当发动机传递到变速器的扭转振动过大时，会导致变速器齿轮产生敲击声，降低车内声品质，容易引起消费者抱怨。
[0003] 液力变矩器是AT自动变速器的重要组成部件，它在传递发动机扭矩的同时，衰减发动机产生的扭转振动，以抑制变速器齿轮敲击噪声。扭转振动隔振率是液力变矩器的重要性能指标之一。在动力传动系统NVH性能调教和机车匹配过程中，能准确预测扭转振动隔振率尤为重要。
[0004] 液力变矩器有三种工作状态，分别是闭锁、解锁和滑差。目前，针对液力变矩器扭转振动隔振率的分析主要集中于闭锁状态，而缺乏预测解锁和滑差状态的扭转振动隔振率的有效方法。因此，迫切需要建立能够预测液力变矩器工作在解锁和滑差状态的扭转振动隔振率的数值建模，并确定计算方法。
发明内容
[0005] 本发明公开的液力变矩器解锁和滑差工况扭转振动隔振率计算方法及可读存储介质，克服目前存在的技术问题。
[0006] 本发明公开的液力变矩器解锁和滑差工况扭转振动隔振率计算方法，包括如下步骤，
[0007] S1：简化液力变矩器在解锁和滑差状态时的动力传动系统旋转部件，依据各旋转部件的设计参数，计算旋转部件的等效惯量、扭转刚度和阻尼；
[0008] S2：通过台架实验得到液力变矩器的特性曲线，依据特性曲线计算得到液力变矩器泵轮、涡轮的动态扭矩；
[0009] S3:通过台架实验得到锁止离合器摩擦片的摩擦系数，计算锁止离合器的摩擦扭矩；
[0010] S4：以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程；
[0011] S5：求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。
[0012] 进一步地，
[0013] S5：通过台架实验采集得到发动机激励扭矩，通过实车道路实验采集得到车辆等效阻力矩，求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。
[0014] 进一步地，
[0015] S1：计算以下旋转部件的等效惯量：
[0016] 1)发动机曲轴、驱动盘和液力变矩器泵轮的等效惯量Jp；
[0017] 2)锁止离合器的等效惯量Jc；
[0018] 3)液力变矩器涡轮的等效惯量Jtb；
[0019] 4)变速器齿轮的等效惯量Jtm；
[0020] 5)轮胎的等效惯量Jw；
[0021] 6)整车的等效惯量Jv；
[0022] 计算以下旋转部件的扭转刚度和阻尼：
[0023] 1)锁止离合器的扭转刚度和阻尼Kc,Cc；
[0024] 2)变速器输入轴的扭转刚度和阻尼Kis,Cis；
[0025] 3)驱动轴的等效扭转刚度和阻尼Kds,Cds；
[0026] 4)轮胎的等效扭转刚度和阻尼Kw,Cw。
[0027] 进一步地，
[0028] S2：通过台架实验得到液力变矩器的特性曲线，依据特性曲线计算得到液力变矩器泵轮、涡轮的动态扭矩；
[0029] 解锁和滑差状态时液力变矩器泵轮、涡轮的动态扭矩特性可表达为：
[0030]
[0031] 式中：Tp‑泵轮扭矩，Tt‑涡轮扭矩，wp‑泵轮角速度，wt‑涡轮角速度；a0,a1,a2,b0,b1,b2为待定系数，由液力变矩器的特性曲线可求得，而液力变矩器的特性曲线由单体台架实验得到，包括变矩比、泵轮扭矩容量系数和效率。
[0032] 进一步地，
[0033] S3中:通过台架实验得到锁止离合器摩擦片的摩擦系数，计算锁止离合器的摩擦扭矩；
[0034] 基于Karnopp摩擦模型计算解锁和滑差状态的锁止离合器摩擦扭矩，计算表达式为：
[0035]
[0036] 式中：Ns‑摩擦面数目，R‑摩擦片有效半径，Fn‑法向压力， 相对滑摩速度，μ为动摩擦系数，可由下式求得：
[0037]
[0038] 式中，μ0‑初始动摩擦系数，ms‑摩擦系数变化梯度，均由台架实验采集得到。
[0039] 进一步地，
[0040] S4中：以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程；
[0041] 传动系统的系统动力学控制方程为：
[0042]
[0043] 其中， 为角位移矢量；
[0044] 惯量矩阵J为：
[0045] J＝diag(Jp,Jc,Jtb,Jtm,Jw,Jv) (8)
[0046] 阻尼矩阵C为：
[0047]
[0048] 刚度矩阵K为：
[0049]
[0050] 系统激励矩阵T为：
[0051] T＝{Te‑Tp‑Tf Tf Tt 0 0 ‑Tv}T (11)
[0052] Te为发动机扭矩，Tv为整车等效阻力矩。
[0053] 进一步地，
[0054] S5中：通过台架实验采集得到发动机激励扭矩，通过实车道路实验采集得到车辆等效阻力矩；车辆加速过程的发动机扭矩Te可用下式表达：
[0055] Te＝Tavg+T2rd×chirp(t,f1,t1,f2) (12)
[0056] 式中，Tavg‑发动机平均扭矩，T2rd‑发动机二阶波动扭矩，chirp(t,f1,t1,f2)表示从0到t1时间段内，扭矩波动频率由f1扫频到f2；发动机平均扭矩和发动机二阶波动扭矩均由发动机台架实验采集得到；
[0057] 车辆加速过程的整车等效阻力矩Tv可用下式表达：
[0058]
[0059] 式中，θv‑轮胎角速度，A,B,D为待定系数，通过实车道路实验确定。
[0060] 本发明有益技术效果为：提供的数值建模方法，以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程；求解系统动力学方程，计算液力变矩器在解锁和滑差状态的扭转振动隔振率。能够预测液力变矩器在解锁和滑差状态的扭转振动隔振率，为控制发动机传递到变速器的扭转振动和变速器齿轮敲击声提供理论依据。
附图说明
[0061] 图1为AT前驱车传动系统结构组成图；
[0062] 图2为液力变矩器原始特性曲线图；
[0063] 图3为液力变矩器解锁和滑差状态的传动系统非线性集中参数模型图；
[0064] 图4为随时间变化的发动机和变速器角速度求解结果图；
[0065] 图5为随转速变化的发动机和变速器角加速度求解结果图；
[0066] 图6为随转速变化的液力变矩器扭转振动隔振率求解结果图。
具体实施方式
[0067] 下面结合附图对本发明做详细说明。
[0068] 如图1所示为AT前驱车传动系统的一个实施例，系统包括：发动机1，驱动盘2，锁止离合器3，泵轮4，涡轮5，变速器输入轴6，变速器齿轮系7，驱动轴8和轮胎9。动力传递路线为：从发动机输出到驱动盘，从驱动盘输出到泵轮，再到涡轮，由涡轮到变速器输入轴，延续至变速器齿轮系到驱动轴，再到轮胎。其中，锁止离合器摩擦片也会传递部分扭矩。
[0069] 以下实施例是液力变矩器解锁和滑差工况扭转隔振率计算方法，包括以下步骤：
[0070] S1：简化液力变矩器在解锁和滑差状态时的动力传动系统旋转部件，依据各旋转部件的设计参数，计算旋转部件的等效惯量、扭转刚度和阻尼。
[0071] 具体需要计算以下旋转部件的等效惯量：
[0072] (1)发动机曲轴1、驱动盘2和液力变矩器泵轮4的等效惯量—Jp；
[0073] (2)锁止离合器3的等效惯量—Jc；
[0074] (3)液力变矩器涡轮5的等效惯量—Jtb；
[0075] (4)变速器齿轮7的等效惯量—Jtm；
[0076] (5)轮胎9的等效惯量—Jw；
[0077] (6)整车的等效惯量—Jv。
[0078] 还要计算以下旋转部件的扭转刚度和阻尼：
[0079] (1)锁止离合器3的扭转刚度和阻尼—Kc,Cc；
[0080] (2)变速器输入轴6的扭转刚度和阻尼—Kis,Cis；
[0081] (3)驱动轴8的等效扭转刚度和阻尼—Kds,Cds；
[0082] (4)轮胎9的等效扭转刚度和阻尼—Kw,Cw。
[0083] S2：通过台架实验得到液力变矩器的特性曲线，依据特性曲线计算得到液力变矩器泵轮、涡轮的动态扭矩。
[0084] 具体地，解锁和滑差状态的液力变矩器泵轮、涡轮的动态扭矩特性可表达为：
[0085]
[0086] 式中：Tp‑泵轮扭矩，Tt‑涡轮扭矩，wp‑泵轮角速度，wt‑涡轮角速度。a0,a1,a2,b0,b1,b2为待定系数，由液力变矩器的特性曲线可求得。
[0087] 其中液力变矩器的特性曲线由单体台架实验得到，如图2所示，包括：
[0088] 变矩比k为：
[0089] k＝Tt/Tp (2)
[0090] 泵轮扭矩容量系数Cp为：
[0091] Cp＝Tp/np2 (3)
[0092] 效率η为：
[0093] η＝Ttnt/Tpnp (4)
[0094] 式中：np‑泵轮转速，nt‑涡轮转速。
[0095] 通过液力变矩器的特性曲线可得到一系列关于泵轮扭矩、涡轮扭矩、泵轮角速度和涡轮角速度的数值，代入式(1)即可计算求得a0,a1,a2,b0,b1,b2，继而得到液力变矩器解锁和滑差状态的泵轮、涡轮动态扭矩。
[0096] S3：通过台架实验得到锁止离合器摩擦片的摩擦系数，计算锁止离合器的摩擦扭矩。
[0097] 基于Karnopp摩擦模型计算解锁和滑差状态的锁止离合器摩擦扭矩，计算表达式为：
[0098]
[0099] 式中：Ns‑摩擦面数目，R‑摩擦片有效半径，Fn‑法向压力， 相对滑摩速度，μ为动摩擦系数，可由下式求得：
[0100]
[0101] 式中，μ0‑初始动摩擦系数，ms‑摩擦系数变化梯度，均由台架实验采集得到。
[0102] S4：以液力变矩器泵轮、涡轮的动态扭矩以及锁止离合器摩擦片的摩擦扭矩替代液力变矩器解锁和滑差状态时实际传递的扭矩，建立系统动力学控制方程。
[0103] 当液力变矩器处于解锁和滑差状态时，发动机扭矩主要通过两条路径传递到变速器，一是通过液力传递，二是通过锁止离合器摩擦传递。其中，液力传递的扭矩可用泵轮、涡轮的动态扭矩表示，锁止离合器摩擦传递的扭矩可用摩擦扭矩表示。需要注意的是，当液力变矩器处于完全断开状态时，锁止离合器摩擦扭矩Tf为0。
[0104] 液力变矩器处于解锁和滑差状态的传动系统可以等效为包含6个自由度的系统模型，如图3所示。除液力变矩器外，其余转子之间采用刚度和阻尼连接。
[0105] 传动系统的系统动力学控制方程为：
[0106]
[0107] 其中，θ＝{θp,θc,θtb,θtm,θw,θv}T为角位移矢量，上标“T”表示矩阵转置。
[0108] 惯量矩阵J为：
[0109] J＝diag(Jp,Jc,Jtb,Jtm,Jw,Jv) (8)
[0110] 阻尼矩阵C为：
[0111]
[0112] 刚度矩阵K为：
[0113]
[0114] 系统激励矩阵T为：
[0115] T＝{Te‑Tp‑Tf Tf Tt 0 0 ‑Tv}T (11)
[0116] Te为发动机扭矩，Tv为整车等效阻力矩。
[0117] S5：通过台架实验采集得到发动机激励扭矩，通过实车道路实验采集得到车辆等效阻力矩，求解系统动力学方程，根据求解方程得到的解，得到车辆加速过程中随时间变化的发动机和变速器角速度(如图4所示)，进一步得到随转速变化的发动机和变速器角加速度(如图5所示)，再计算液力变矩器在解锁和滑差状态的扭转振动隔振率(如图6所示)。
[0118] 车辆加速过程的发动机扭矩Te可用下式表达：
[0119] Te＝Tavg+T2rd×chirp(t,f1,t1,f2) (12)
[0120] 式中，Tavg‑发动机平均扭矩，T2rd‑发动机二阶波动扭矩，chirp(t,f1,t1,f2)表示从0到t1时间段内，扭矩波动频率由f1扫频到f2。发动机平均扭矩和发动机二阶波动扭矩均由发动机台架实验采集得到。
[0121] 车辆加速过程的整车等效阻力矩Tv可用下式表达：
[0122]
[0123] 式中，θv‑轮胎角速度，A,B,D为待定系数，通过实车道路实验确定。
[0124] 在本实施例中，还提供了一种可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的熟悉目的地退出导航的方法的步骤。
[0125] 本领域普通技术人员应知，所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read‑Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
