一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方法
【专利摘要】本发明的目的在于提供一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方法，根据活塞销中心的瞬时坐标和曲柄销中心的瞬时坐标获取连杆的瞬时长度，再利用连杆瞬时长度的变化量和变化率获取连杆的动态内力；根据曲柄销中心处集中质量的动力学控制方程获取曲柄的瞬时长度，再利用曲柄瞬时长度的变化量和变化率来获取曲柄的动态内力；根据主轴颈中心处转动惯量的动力学控制方程获取曲柄的瞬时角位移，再利用曲柄瞬时角位移的相对变化量和相对变化率来获取曲柄的动态内力矩；最后根据曲柄销中心处连杆动态内力在垂直于曲柄方向的分力和曲柄的瞬时长度预测曲柄的驱动扭矩。本发明可以采用常用的数值方法求解，具有计算规模小、预测速度快的优点。
【专利说明】
一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测 方法
技术领域
[0001] 本发明涉及的是一种预测方法，具体地说是一种机械结构驱动扭矩的预测方法。
【背景技术】
[0002] 曲柄连杆机构是内燃机、斯特林发动机、活塞式压缩机等机械中主要的运动机构 之一，由活塞组、连杆组和曲柄等部件组成。曲柄连杆机构的运动特点是曲柄做旋转运动， 活塞做往复运动，连接活塞和曲柄的连杆做平面运动。
[0003] 曲柄的驱动扭矩代表曲柄连杆机构的做功能力。目前常用的曲柄连杆机构动力学 预测方法是质点力系法。在预测时不考虑零部件的自身振动的影响，将零部件直接简化为 集中于运动中心的质点。考虑各质点的受力状态时，先假设曲柄做等角速度转动，然后根据 质点系的运动学关系确定各质点的运动加速度，再根据牛顿第二定律计算各质点的惯性 力，最后考虑气缸压力和惯性力的合力作用。曲柄连杆机构具有较高的运行速度，并且受到 周期性变化的气缸压力和惯性力的作用。在这样的工作条件下，曲柄、连杆等零部件自身可 能发生一定振动。这些现象不仅会影响曲柄连杆机构自身的运行速度和加速度，而且会影 响其负载机构的工作平稳性。
【发明内容】
[0004] 本发明的目的在于提供能够提高预测精度的一种考虑曲柄和连杆振动的曲柄连 杆机构驱动扭矩的预测方法。
[0005] 本发明的目的是这样实现的：
[0006] 本发明一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方法，其特征 是：
[0007] (1)获取曲柄连杆机构中零部件的质量，包括曲柄销质量MpinC、以曲柄销中心线为 对称中心并且附着在曲柄销上面的两个曲柄臂的质量Mac、连杆质量Mr、活塞质量Mpistcm、活 塞销质量MpinP、活塞环质量M Ring;获取曲柄连杆机构中零部件的转动惯量，包括主轴颈转动 惯量Ijc>urnal、以主轴颈中心线为对称中心并且附着在主轴颈上面的两个曲柄臂的转动惯量 IA J、两个曲柄臂在主轴颈与曲柄销之间部分的转动惯量IAW;
[0008] (2)获取曲柄连杆机构中零部件的刚度参数，包括曲柄的弯曲刚度系数Kc、曲柄的 扭转刚度系数kc和连杆的纵向刚度系数如，并分别取临界阻尼系数的2%作为相应的阻尼系 数；
[0009] (3)获取曲柄连杆机构中零部件的几何尺寸参数，包括曲柄的长度1CQ、活塞顶横截 面直径d PQ、连杆小端中心与连杆重心的距离Irl、连杆大端中心与连杆重心的距离1RB、连杆 小端中心与大端中心的原始距离Iro ;
[0010] (4)获取曲柄连杆机构的驱动力源数据列表，即气缸压力fc(t);
[0011] (5)获取曲柄连杆机构的负载运行角速度GH(t);
[0012] (6)获取活塞的初始位移和初始速度，曲柄的初始角位移和初始角速度；
[0013] (7)离散连杆的质量，即将连杆质量向连杆小端中心和大端中心离散；
[0 014 ] (8)建立曲柄连杆机构驱动扭矩的预测模型：活塞销中心处的集中质量Mp包括活 塞的质量Mpiston、活塞销的质量Mpinp、活塞环的质量MRing、连杆小端中心的简化质量Mlittle，曲 柄销中心处的集中质量Me包括曲柄销的质量Mpinc、以曲柄销中心线为对称中心并且附着在 曲柄销上面的两个曲柄臂的质量Mac、连杆大端中心的简化质量MBig，主轴颈中心的转动惯量 1C包括主轴颈的转动惯量Ijcmrnal、以主轴颈中心线为对称中心并且附着在主轴颈上面的两 个曲柄臂的转动惯量IAJ、两个曲柄臂在主轴颈与曲柄销之间部分的转动惯量IAW、集中质量 Me的转动惯量I Crank ;
[0015] (9)输入当前时刻驱动力源数据，即气缸压力fG(t);
[0016] (10)根据活塞销中心处集中质量MP的动力学控制方程获取当前时刻活塞销中心 的瞬时坐标(XP，yP)，利用曲柄的瞬时角位移0c(t)和瞬时长度lc(t)获取当前时刻曲柄销中 心的瞬时坐标(xc，yc);
[0017] (11)根据活塞销中心的瞬时坐标(XP，yP)、曲柄销中心的瞬时坐标(xc，ye)获取连 杆的瞬时长度lR(t)，利用连杆瞬时长度lR(t)的变化量[lRQ-lR(t)]和变化率dl R(t)/dt来获 取连杆的动态内力；
[0018] (12)根据曲柄销中心处集中质量Me的动力学控制方程获取曲柄的瞬时长度lc(t)， 利用曲柄瞬时长度lc(t)的变化量[l CQ-lc(t)]和变化率dlc(t)/dt来获取曲柄的动态内力；
[0019] (13)根据主轴颈中心处转动惯量Ic的动力学控制方程获取曲柄的瞬时角位移0C (t)，利用瞬时角位移0c(t)的相对变化量[0c(t)-0L(t)]和相对变化率[d0 c(t)/dt-coL(t)] 来获取曲柄的动态内力矩；
[0020] (14)预测当前时刻曲柄连杆机构的驱动扭矩TT(t);
[0021] (15)返回步骤(9)，进行下一个时刻的预测过程，直至预测结束。
[0022] 本发明还可以包括：
[0023] 1、获取连杆的瞬时长度lR(t)方法为：
[0025] 利用连杆瞬时长度lR(t)的变化量[lRQ-l R(t)]和变化率dlR(t)/dt来表示连杆的纵 向振动，即通过下式获取连杆的动态内力：
[0027] 其中，KR与CR分别为连杆的纵向刚度系数和阻尼系数;连杆动态内力FR(t)沿连杆 中心线传递到曲柄销中心处，分为互相垂直的二个分力F T(t)和FN(t)，FT(t)的作用方向与 曲柄圆相切，F N(t)的作用方向与曲柄圆垂直，两者的大小分别为： | ^(〇 = F,,(t)sin[0c(t) + 0r{t)]
[0028] ] 、 _ ^ 1 [Fv{/) = Fa,(0cos[^.(0 + ^,(0]
[0029] 其中，0c(t)为曲柄的瞬时角位移，0P(t)为连杆的瞬时摆角。
[0030] 2、获取曲柄的瞬时长度lc(t)方法为：
[0032]其中，Fc (t)为曲柄的动态内力，if (〇为集中质量Me的离心惯性力；
[0033]利用曲柄瞬时长度lc(t)的变化量[lCQ-lc(t)]和变化率dlc(t)/dt来表示曲柄的弯 曲振动，即通过下式获取曲柄的动态内力：
[0035]其中，Kc与Cc分别为曲柄的弯曲刚度系数和阻尼系数。
[0036] 3、获取曲柄的瞬时角位移0c(t)方法为：
[0038]其中，TT(t)为曲柄的驱动扭矩，Tc(t)为曲柄的动态内力矩；
[0039] 利用瞬时角位移0C(t)的相对变化量[0c(t)-0L(t)]和相对变化率[d0 c(t)/dt-?L (t)]来表示曲柄的扭转振动，即通过下式获取曲柄的动态内力矩：
其中，kc与cc分别为曲柄的扭转刚度系 数和阻尼系数，《L(t)为负载的运行角速度，0L(t)为负载在当前时刻的角位移，其大小为 〇L(t) ? t，t为时间。
[0041] 4、当前时刻曲柄连杆机构的驱动扭矩TT(t)的获取方法：
[0042] TT(t) =FT(t)lc(t)。
[0043]本发明的优势在于：
[0044] (1)本发明提供的预测方法考虑了曲柄和连杆振动的影响。其中，曲柄的振动包括 曲柄的弯曲振动和扭转振动，连杆的振动为纵向振动。这相对于以往的不考虑零部件振动 影响的动态扭矩预测方法，本发明提供的预测方法更完备。
[0045] (2)本发明提供的预测方法不需要复杂的数学理论，并且可以采用常用的数值方 法求解，例如有限差分法和四阶龙格库塔法，具有计算规模小、计算速度快的优点。
[0046] (3)按照曲柄的排列相位，根据本发明提供的预测方法可以方便得到由多个曲柄 连杆单元组成的曲柄连杆机构驱动扭矩的预测方法。
【附图说明】
[0047]图1为本发明的流程图；
[0048]图2为本发明的预测模型图。
【具体实施方式】
[0049] 下面结合附图举例对本发明做更详细地描述：
[0050] 结合图1~2,本发明按如下步骤实现：
[0051 ]步骤1 :获取曲柄连杆机构中零部件的质量和转动惯量参数，包括曲柄销质量 MPinC、以曲柄销中心线为对称中心并且附着在曲柄销上面的两个曲柄臂的部分质量Mac、连 杆质量Mr、活塞质量Mpiston、活塞销质量Mpinp、活塞环质量MRing。获取曲柄的转动惯量，包括主 轴颈转动惯量Ijcmrrml、以主轴颈中心线为对称中心并且附着在主轴颈上面的两个曲柄臂部 分的转动惯量IA J、两个曲柄臂在主轴颈与曲柄销之间部分的转动惯量IAW。
[0052]步骤2:获取曲柄连杆机构中零部件的刚度参数和阻尼参数，包括曲柄的弯曲刚度 系数Kc、曲柄的扭转刚度系数kc和连杆的纵向刚度系数KR，并分别取临界阻尼系数的2%作 为相应的阻尼系数。
[0053]步骤3:获取曲柄连杆机构中零部件的几何尺寸参数，包括曲柄的长度1CQ、活塞顶 横截面直径dPQ、连杆小端中心与连杆重心的距离Irl、连杆大端中心与连杆重心的距离Irb、 连杆小端中心与大端中心的原始距离Iro。
[0054]步骤4:获取曲柄连杆机构的驱动力源数据列表，即气缸压力fG(t)。
[0055] 步骤5:获取曲柄连杆机构的负载运行角速度〇u(t)。
[0056] 步骤6:获取活塞的初始位移和初始速度，曲柄的初始角位移和初始角速度。
[0057]步骤7 :离散连杆的质量，即将连杆质量Mr向连杆小端中心和大端中心离散。离散 后的连杆小端质量为]^1；1；16，连杆大端质量为]/^，其大小分别为]/[1^1；1；16=]\ /[[?1[^/1[?()，]/^ = Mr1rl/1ro〇
[0058]步骤8:建立如图2所示的曲柄连杆机构驱动扭矩的预测模型。其中，活塞销中心处 的集中质量MP包括活塞的质量Mpiston、活塞销的质量MpinP、活塞环的质量MRing、连杆小端中心 的简化质量MLittle。曲柄销中心处的集中质量Me包括曲柄销的质量MpinC、以曲柄销中心线为 对称中心并且附着在曲柄销上面的两个曲柄臂的部分质量Mac、连杆大端中心的简化质量 MBig。主轴颈中心的转动惯量1C包括主轴颈的转动惯量I Journal、以主轴颈中心线为对称中心 并且附着在主轴颈上面的两个曲柄臂部分的转动惯量IAJ、两个曲柄臂在主轴颈与曲柄销之 间部分的转动惯量IAW、集中质量Me的转动惯量I Crank。
[0059]步骤9:输入当前时刻驱动力源数据，即气缸压力fG(t)。
[0060]步骤10:根据活塞销中心处集中质量MP的动力学控制方程获取当前时刻活塞销中 心的瞬时坐标(XP，yP)，利用曲柄的瞬时角位移0c(t)和瞬时长度l C(t)获取当前时刻曲柄销 中心的瞬时坐标(xc, ye)。不考虑活塞的横向运动，所以有xp = 0。
[0061 ]步骤11:根据活塞销中心的瞬时坐标(XP，yP)、曲柄销中心的瞬时坐标(Xc，ye)获取 连杆的瞬时长度lR(t)，即：
(1)
[0063] 利用连杆瞬时长度lR(t)的变化量[lRQ-l R(t)]和变化率dlR(t)/dt来表示连杆的纵 向振动。即通过下式获取连杆的动态内力：
(2)
[0065]其中，KR与CR分别为连杆的纵向刚度系数和阻尼系数。连杆动态内力FR(t)沿连杆 中心线传递到曲柄销中心处，分为互相垂直的二个分力FT(t)和FnUKFtU)的作用方向与 曲柄圆相切，称为切向力。F N(t)的作用方向与曲柄圆垂直，称为法向力。两者的大小分别 为： |F, (t) = F,,(〇sin[ft (/) + 0,Xt)\
[0066] 丨厂、，⑴=/^(Ocosh ⑴ + 沒"(〇] (3)
[0067] 其中，0c(t)为曲柄的瞬时角位移，0P(t)为连杆的瞬时摆角。
[0068] 步骤12:根据曲柄销中心处集中质量Me的动力学控制方程获取曲柄的瞬时长度lc ⑴，即
(4)
[0070] 其中，FC (t)为曲柄的动态内力，为集中质量Me的离心惯性力。
[0071] 利用曲柄瞬时长度lc(t)的变化量[lCQ-lC(t)]和变化率dlc(t)/dt来表示曲柄的弯 曲振动，即通过下式获取曲柄的动态内力：
C5)
[0073]其中，Kc与Cc分别为曲柄的弯曲刚度系数和阻尼系数。
[0074]步骤13:根据主轴颈中心处转动惯量Ic的动力学控制方程获取曲柄的瞬时角位移 9c(t)，即
(6)
[0076] 其中，TT(t)为曲柄的驱动扭矩，Tc(t)为曲柄的动态内力矩。
[0077] 利用瞬时角位移0C(t)的相对变化量[0C(t)-0L(t)]和相对变化率[d0 C(t)/dt-c0L (t)]来表示曲柄的扭转振动，即通过下式获取曲柄的动态内力矩：
(7)
[0080] 其中，kc与cc分别为曲柄的扭转刚度系数和阻尼系数，〇u(t)为负载的运行角速 度，t(t)为负载在当前时刻的角位移，其大小为《L(t) ? t，t为时间。
[0081] 步骤14:预测当前时刻曲柄连杆机构的驱动扭矩TT(t)，即
[0082] TT(t)=Fx(t)lc(t) (8)
[0083]其中，Ft (t)为曲柄销中心处连杆动态内力Fr (t)在垂直于曲柄方向的分力。
[0084] 步骤15:返回步骤(9)，进行下一个时刻的预测过程。
【主权项】
1. 一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方法，其特征是： (1) 获取曲柄连杆机构中零部件的质量，包括曲柄销质量MPinC、以曲柄销中心线为对称 中心并且附着在曲柄销上面的两个曲柄臂的质量Mac、连杆质量Mr、活塞质量Mpistcm、活塞销 质量Mpinp、活塞环质量MRing ;获取曲柄连杆机构中零部件的转动惯量，包括主轴颈转动惯量 Ijcmrnal、以主轴颈中心线为对称中心并且附着在主轴颈上面的两个曲柄臂的转动惯量IAJ、 两个曲柄臂在主轴颈与曲柄销之间部分的转动惯量I?; (2) 获取曲柄连杆机构中零部件的刚度参数，包括曲柄的弯曲刚度系数Kc、曲柄的扭转 刚度系数k和连杆的纵向刚度系数KR，并分别取临界阻尼系数的2%作为相应的阻尼系数； (3) 获取曲柄连杆机构中零部件的几何尺寸参数，包括曲柄的长度lco、活塞顶横截面直 径dPQ、连杆小端中心与连杆重心的距离1 RL、连杆大端中心与连杆重心的距离1RB、连杆小端 中心与大端中心的原始距离Iro ; (4) 获取曲柄连杆机构的驱动力源数据列表，即气缸压力fc(t); (5) 获取曲柄连杆机构的负载运行角速度〇u(t); (6) 获取活塞的初始位移和初始速度，曲柄的初始角位移和初始角速度； (7) 离散连杆的质量，即将连杆质量向连杆小端中心和大端中心离散； (8) 建立曲柄连杆机构驱动扭矩的预测模型:活塞销中心处的集中质量MP包括活塞的质 量Mpiston、活塞销的质量MpinP、活塞环的质量MRing、连杆小端中心的简化质量Mlittle，曲柄销中 心处的集中质量Me包括曲柄销的质量MPinC、以曲柄销中心线为对称中心并且附着在曲柄销 上面的两个曲柄臂的质量Mac、连杆大端中心的简化质量M Big，主轴颈中心的转动惯量I c包括 主轴颈的转动惯量Ijcmrnal、以主轴颈中心线为对称中心并且附着在主轴颈上面的两个曲柄 臂的转动惯量Iaj、两个曲柄臂在主轴颈与曲柄销之间部分的转动惯量Iaw、集中质量Me的转 动惯量I Crank; (9) 输入当前时刻驱动力源数据，即气缸压力fG(t); (1〇)根据活塞销中心处集中质量MP的动力学控制方程获取当前时刻活塞销中心的瞬时 坐标(Xp，yp)，利用曲柄的瞬时角位移9c(t)和瞬时长度lc(t)获取当前时刻曲柄销中心的瞬 时坐标(xc，yc); (11) 根据活塞销中心的瞬时坐标(xp，yp)、曲柄销中心的瞬时坐标(xc，yc)获取连杆的瞬 时长度lR(t)，利用连杆瞬时长度l R(t)的变化量[lRQ-lR(t)]和变化率dlR(t)/dt来获取连杆 的动态内力； (12) 根据曲柄销中心处集中质量Me的动力学控制方程获取曲柄的瞬时长度lc(t)，利用 曲柄瞬时长度lc⑴的变化量[l c『lc(t)]和变化率dlc(t)/dt来获取曲柄的动态内力； (13) 根据主轴颈中心处转动惯量Ic的动力学控制方程获取曲柄的瞬时角位移0c(t)，利 用瞬时角位移9c(t)的相对变化量[0 c(t)-0L(t)]和相对变化率[d0c(t)/dt-coL(t)]来获取 曲柄的动态内力矩； (14) 预测当前时刻曲柄连杆机构的驱动扭矩TT(t); (15) 返回步骤(9)，进行下一个时刻的预测过程，直至预测结束。2. 根据权利要求1所述的一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方 法，其特征是： 获取连杆的瞬时长度1r (t)方法为：利用连杆瞬时长度lR(t)的变化量[lRQ-lR(t)]和变化率dlR(t)/dt来表不连杆的纵向振 动，即通过下式获取连杆的动态内力：其中，KR与CR分别为连杆的纵向刚度系数和阻尼系数;连杆动态内力FR(t)沿连杆中心 线传递到曲柄销中心处，分为互相垂直的二个分力FT(t)和FN(t)，FT(t)的作用方向与曲柄 圆相切，F N(t)的作用方向与曲柄圆垂直，两者的大小分别为：其中，0c(t)为曲柄的瞬时角位移，0P(t)为连杆的瞬时摆角。3. 根据权利要求2所述的一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方 法，其特征是： 获取曲柄的瞬时长度lc (t)方法为：其中，Fc (t)为曲柄的动态内力，为集中质量Me的离心惯性力； 利用曲柄瞬时长度lc(t)的变化量[lct)-lc(t)]和变化率dlc(t)/dt来表不曲柄的弯曲振 动，即通过下式获取曲柄的动态内力：其中，Kc与Cc分别为曲柄的弯曲刚度系数和阻尼系数。4. 根据权利要求3所述的一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方 法，其特征是： 获取曲柄的瞬时角位移9c(t)方法为：其中，TT(t)为曲柄的驱动扭矩，Tc(t)为曲柄的动态内力矩； 利用瞬时角位移(t)的相对变化量[0c(t)-0L(t)]和相对变化率[d0c(t)/dt- ω L(t)] 来表示曲柄的扭转振动，即通过下式获取曲柄的动态内力矩：其中，kc与cc分别为曲柄的扭转刚度系数和阻尼系数，〇u(t)为负载的运行角速度， (t)为负载在当前时刻的角位移，其大小为〇H(t) · t，t为时间。5. 根据权利要求4所述的一种考虑曲柄和连杆振动的曲柄连杆机构驱动扭矩的预测方 法，其特征是： 当前时刻曲柄连杆机构的驱动扭矩TT(t)的获取方法： Ττ(?) =FT(t)lc(t) 〇
【文档编号】G06F17/50GK105930592SQ201610265552
【公开日】2016年9月7日
【申请日】2016年4月26日
【发明人】国杰, 张文平, 张新玉, 明平剑, 柳贡民, 曹贻鹏
【申请人】哈尔滨工程大学