汽车悬架减振器能量特性测试系统及计量方法
【专利摘要】本发明创造公开一种汽车悬架减振器能量特性测试系统及计量方法，该系统包括架体，设置在架体上可产生正弦运动的曲柄连杆滑块机构，用于固定减振器的夹具和下安装支座，速度传感器，及用于测量阻尼力的力传感器；曲柄连杆滑块机构用于实现减振器的拉伸与压缩行程往复运动，具有曲柄，曲轴连杆可绕着曲柄的另一端点进行周向运动；曲轴连杆的顶端与下安装支座连接从而带动减振器随着曲轴连杆作上下运动，滑块在曲轴连杆的带动下进行上下往复运动。该计量方法基于汽车悬架减振器能量特性测试系统将计算机图像识别技术应用于减振器示功图面积的求解与标定，并提出一种减振器吸收能量定量分析方法，为减振器的性能检测提供新的评价方法及量化标准。
【专利说明】
汽车悬架减振器能量特性测试系统及计量方法
技术领域
[0001] 本发明创造涉及汽车悬架减振器领域，具体涉及一种汽车悬架减振器能量特性测 试系统及计量方法。
【背景技术】
[0002] 减振器是汽车悬架中主要的减振元件，与弹性元件并联安装，其性能直接关系到 汽车的舒适性和安全性。汽车悬架是保证汽车行驶平顺性与操作稳定性的关键总成，减振 器作为其中的核心器件承担着非常大的职责。因此，无论对于出厂前的检测还是日后的维 修保养，汽车减振器的性能和失效模式的研究是十分必要的。
[0003] 减振器在生产及维修时，都需要经过测试来明确其性能。减振器的性能主要是通 过减振器的外特性来表征的。通过减振器示功机可以得到减振器的示功图和速度特性图。 可以根据示功图曲线的形状来判断减振器产品是否存在缺陷，即分析减振器是否失效。
[0004] 目前，已经可以通过减振器示功机来测试最大拉伸阻力和最大压缩阻力这两个参 数合格与否，其识别方式为"门檻法"，超过了所设的"门滥值"就不合格。该识别方法的应用 过于单一，只能识别出通过最大拉伸阻力、压缩阻力表征的缺陷类型。一旦遇到最大拉伸阻 力、压缩阻力都在合格范围内，却存在着其他缺陷类型的示功图，现有检测方式就无法自动 识别。在这种情况下，仍然需要依靠人工在线的方式逐一查验，根据检测人员的主观经验进 行分类识别。
[0005] 可以看出依靠人工识别的方式不利于减振器的大批量生产与检测，当前减振器示 功机所配置的自动识别功能也较为粗糙。而且减振器在复原、压缩行程中所测试的多数点 的阻尼特性无法得到体现，且减振器在一个往返行程内所衰减车辆振动的能量多少往往也 被忽略。
【发明内容】
[0006] 本发明创造目的是提供一种汽车悬架减振器能量特性测试系统及计量方法，该系 统用于减振器能量特性的测试，采用曲柄连杆滑块机构对减振器施加激励;采用速度传感 器测试减振器的位移或速度;采用力传感器测试减振器的动静态阻尼力和冲击载荷。通过 测试既可以获得减振器外能量特性，又可以提高对减振器是否存在性能畸变的识别能力， 提高曲柄连杆滑块机构、减振器的装夹及其机构是否存在间隙的识别能力。汽车悬架减振 器能量特性计量方法基于汽车悬架减振器能量特性测试系统将计算机图像识别技术应用 于减振器示功图面积的求解与标定，并提出一种减振器吸收能量定量分析方法，为减振器 的性能检测提供新的评价方法及量化标准。
[0007] 本发明创造采用的技术方案为：
[0008] 汽车悬架减振器能量特性测试系统，包括架体，设置在架体上可产生正弦运动的 曲柄连杆滑块机构，测试时用于固定减振器上吊环的夹具和用于固定减振器下吊环的下安 装支座，速度传感器，及用于测量阻尼力的力传感器;所述的曲柄连杆滑块机构用于实现减 振器的拉伸与压缩行程往复运动，具有曲柄，曲柄的一端固定在架体上，曲柄的另一端与曲 轴连杆底端相连接，所述曲轴连杆可绕着曲柄的另一端点进行周向运动；曲轴连杆的顶端 与下安装支座连接从而带动减振器随着曲轴连杆作上下运动，曲轴连杆上还设有与之配合 的滑块，所述滑块在曲轴连杆的带动下进行上下往复运动;所述速度传感器设置在曲轴连 杆上用于测试由曲柄转动产生的速度。
[0009] 所述的汽车悬架减振器能量特性测试系统，所述力传感器固连于架体的上支座 上，并与夹具串联;所述力传感器包括应变片式力传感器和压电式力传感器，所述应变片式 力传感器和压电式力传感器串联。
[0010] 所述的汽车悬架减振器能量特性测试系统，还包括减振器侧向力装置，所述减振 器侧向力装置与减振器固定连接，具有测力环，和固装在架体上的用于支撑钢索的滚轮，所 述侧力环的一端固定连接在减振器的中间部位，侧力环的另一端与钢索连接，钢索通过滚 轮连接砝码;减振器随下面的曲柄连杆滑块机构作上下运动时，测力环通过钢索带动砝码 随减振器作小范围上下运动。
[0011] 汽车悬架减振器能量特性计量方法，应用上述汽车悬架减振器能量特性测试系 统，包括以下步骤：
[0012] 1)建立减振器数学模型，包括建立减振器压缩行程和复原行程数学模型；
[0013] 2)建立减振器示功图面积数学模型，基于汽车悬架减振器能量特性测试系统，在 激振行程和频率下，减振器的压缩行程和复原行程产生示功曲线，该示功曲线组成的图形 即减振器示功图，该示功曲线所围成的面积即减振器示功图面积;基于减振器数学模型建 立减振器示功图面积数学模型；
[0014] 3)示功图面积能量标定，
[0015] 3.1)基于Matleb平台利用计算机图像识别技术对减振器示功图进行图像处理，包 括图像读取、图像灰度化及图像二值化，检测减振器示功图像素点个数并求解减振器示功 图面积；
[0016] 3.2)用LabVIEW对减振器示功图面积数学模型编程，并利用减振器在汽车悬架减 振器能量特性测试系统中所采集的数个采样点的速度值及阻尼力值数据，计算求解减振器 在一个工作行程中吸收的总能量。
[0017]所述的汽车悬架减振器能量特性计量方法，步骤1)建立减振器数学模型是根据减 振器工作机理，分别建立减振器的各行程开阀前后的分段数学模型；
[0018] 工作行程中活塞所受阻尼力为：
[0019] F = Pi(Ah-Ag)-P2 · Ah+fi (1)
[0020] 油液从上腔流入下腔的流量为Q:
[0021] Q= (Ah-Ag)vi (2)
[0022] 贮油腔油液补偿流入工作缸下腔的流量为Qc:
[0023] Qc = Ag · νι (3)
[0024] 式中：F为减振器受到的阻尼力N;Ah为活塞面积m2;Ag为活塞杆面积复原行 程活塞的运动速度m/^Pi为工作缸上腔的压强Mpa;P2工作缸下腔的压强Mpadi为活塞所受 摩擦力N;
[0025] 1.1)建立复原行程数学模型
[0026] 1.1.1)复原行程开阀前数学模型
[0027] 假设减振器正常工作时油封密闭且无油液泄漏现象，减振器油只从工作缸的上腔 经过活塞上复原阀进入工作缸下腔，根据流体串并联原理，
[0028]油液通过复原阀组件流入下腔的流量&为：
[0030]贮油腔油液通过补偿流入工作缸下腔的流量为Qcl:
[0032 ] 根据油液连续性方程由Q = Qi，Qc = Qcl得：
[0034]则开阀前活塞所受阻尼力^可表示为：
[0036] 式中：Cd为油液流量系数;Am为阀体节流孔面积m2 ;P为油液密度Kg/m3; An为补偿阀 节流孔面积m2; P3为贮油腔压强Mpa;
[0037] 1.1.2)复原行程开阀后数学模型
[0038] 随着复原行程中活塞速度的增加，工作缸上腔的压力也逐渐增加，当上腔油液压 力升高并且大于复原阀组件内阀片的承受压力时，复原阀片变形产生开阀状态，此时减振 器进入开阀后的阻尼状态；
[0039] 开阀后油液由工作缸上腔经过活塞复原阀上的常通孔和阀片发生形变形成的环 状缝隙进入下腔产生的流量为Q 2:
[0041]式中，A2为阀片发生形变产生的缝隙面积m2,根据油液连续性方程由〇 = 〇2得：
[0043] 则开阀后活塞所受阻尼力^可表示为：
[0045] 1.2)压缩行程数学模型的建立
[0046] 减振器压缩行程的数学模型与复原行程的数学模型建立方法同步骤1.1)，根据活 塞运动速度及阻尼力形成方式分别建立开阀前与开阀后减振器数学模型。
[0047]所述的汽车悬架减振器能量特性计量方法，步骤2)减振器示功图面积分为上下两 部分，分别为复原行程产生的示功曲线围成的上部分和减振器压缩行程产生的示功曲线围 成的下部分，可视为微小位移ds和瞬时阻尼力h所构成矩形单元的面积累加之和；
[0048]将减振器单个复原行程所用时间t均分为η等份，由于示功机的激振形式为正弦激 励信号，由速度传感器可测曲柄转动产生速度特性，即减振器活塞速度特性为：
[0049] v = 23TfRcos(23rft) =3rfDcos(23rft) (11)
[0050] 式中：f为输入信号的激振频率Hz ;R为曲柄长度mm;D为减振器总行程mm;
[0051] 减振器活塞位移特性由其速度特性表达式积分得到为：
[0053] 式中，S为减振器活塞位移，
[0054] 由减振器速度传感器可测得第i个时间点与第i+Ι个时间点的活塞运行速度，其表 达式分别为：
[0055] vi = 3TfDcos(23rfti) (13)
[0056] vi+i = 3TfDcos(23rfti+i) (14)
[0057] 对减振器速度表达式求积分可得到第i个时间点与第i + 1个时间点之间的微小位 移ds表达式为
[0059]给定活塞一定的激振频率，在汽车悬架减振器能量特性测试系统上可采集到两个 相对应的数据，即活塞运动速度Vl及减振器瞬时阻尼力Fi表达式在步骤1)数学模型中已 推导;将第i点与第i+Ι点所对应的瞬时阻尼力的均值作为矩形单元的高，单个矩形的面积 可表示为：
[0060] Aii = Fi〇ds (16)
[0061] 式中:Ah为矩形单元的面积，N · mm; 为平均阻尼力，N;
[0062] 因此，复原行程示功曲线围成的面积为：
[0064] 同理可求压缩行程其示功曲线围成的面积Ay，因此该减振器示功图面积可表示 为：
[0065] A=Ai+Ay (18)。
[0066] 所述的汽车悬架减振器能量特性计量方法，步骤3.1)具体包括如下步骤：
[0067] 3.1.1)基于Canny边缘检测算法对减振器示功曲线进行边缘检测，检测该示功曲 线是否平滑，Canny算法可对目标曲线平滑处理，其方法是采用高斯函数：
[0069] 式中，〇为高斯滤波器平滑程度控制参数；
[0070] 3.1.2)平滑后的曲线使用一阶微分算子，获得各像素点处的梯度幅值和方向，点 (i，j)处2个方向的偏导数G x( i，j)和Gy( i，j)分别为：
[0071] Gx(i,j) = [I(i,j+l)-I(i,j)+I(i+l,j+l)-I(i+l,j)]/2 (20)
[0072] Gy(i,j) = [I(i,j)-I(i+l,j)+I(i,j+l)-I(i+l,j+l)]/2 (21);
[0073] 3.1.3)通过细化各像素点的梯度幅值图像，可精确识别示功曲线边缘线，非极大 值抑制即只需保留幅值的局部极大值；
[0074] 3.1.4)使用Canny算子从候选边缘点中检测和连接得到最终的边缘。
[0075] 所述的汽车悬架减振器能量特性计量方法，步骤3.1.3)具体为Canny算子在梯度 幅值图像G中以点（i，j)为中心3X3的邻域内进行插值运算，若点（i，j)处的梯度幅值G(i， j)小于相邻的插值，那么该像素点（i，j)为非边缘点，相反则可视为候选边缘点，从而得到 候选边缘曲线图像。
[0076] 所述的汽车悬架减振器能量特性计量方法，步骤3.1.4)具体为根据图像特点选取 两个合理的阈值，分别为高阈值和低阈值，然后在候选边缘曲线图像中对任一像素点（i，j) 进行检测，如果被检测像素点（i，j)梯度幅值G(i，j)大于所选取的高阈值，可将其视为边缘 点;相反如果其值小于低阈值，那么判定其为非边缘点;若梯度幅值介于高低阈值之间的像 素点，可先将其视为疑似边缘点，接着根据边缘的连通性对其进行判断，若该像素点的邻接 像素中存在边缘点，则认为该像素点也为边缘点，否则为非边缘点。
[0077]本发明创造具有以下有益效果：
[0078]本发明创造汽车悬架减振器能量特性测试系统包括可产生正弦运动的曲柄连杆 滑块机构、测试时固定减振器上下吊环的夹具与下安装支座、测量阻尼力及速度位移的传 感器；曲柄连杆滑块机构能实现减振器拉伸与压缩行程往复运动，曲柄一端固定，连杆可绕 着曲柄另一端点进行周向运动，带动滑块进行上下往复运动。
[0079]力传感器与夹具串联;包括应变片式力传感器和压电式力传感器，两传感器串联 在一起，利用应变片式及压电式力传感器自身特点，对减振器性能测试时产生的静态力及 冲击载荷进行有效测量。应变片式力传感器用来实现减振器静态力、准静态力测试；由于压 电式传感器的自身特点，压电晶体受压后其表面产生电荷，在动态交变应力作用下，电荷可 以供给测量电路一定的电压或电流，因此适合减振器动态力的测试，并可预测减振器空行 程及示功曲线畸变现象。设置在曲柄连杆机构两端的速度/位移传感器用来识别压电式力 传感器力值的方向，可判别其力值产生于减振器运动的哪一拉伸或压缩行程，可有效测量 出该系统的各部件间的工作间隙，以及长期工作后各机构的磨损间隙，自身对系统进行维 护。
[0080] 振器侧向力装置中砝码可以根据整车需求侧向力的大小，可以任意更换，砝码的 大小会影响减振器内部油液摩擦。当减振器随下面的曲轴连杆上下运动的时候，钢索也会 随减振器小范围运动。该侧向力装置可更好的模拟减振器在整车上的工作状态。
【附图说明】
[0081] 图1为减振器示功图面积 [0082]图2为一有畸变减振器不功图
[0083]图3为图2所示图像二值化后减振器示功图 [0084] 图4 LabVIEW程序框图 [0085]图5能量计算结果图
[0086] 图6为实施例汽车悬架减振器能量特性测试系统示意图。
[0087] 图7为实施例汽车悬架减振器能量特性测试系统中侧向力装置工作状态示意图。 [0088] 其中：1-上支座，2-应变片式力传感器，3-压电式力传感器，4-夹具，5-减振器，6-下安装支座，7-滑块，8-连杆，9-速度传感器，10-曲柄，11-侧力环，12-钢索，13-滚轮，14-砝 码。
【具体实施方式】
[0089] 实施例汽车悬架减振器能量特性测试系统，
[0090] 如图6和7所示包括架体，设置在架体上可产生正弦运动的曲柄连杆滑块机构，测 试时用于固定减振器5上吊环的夹具4和用于固定减振器5下吊环的下安装支座6,速度传感 器9,及用于测量阻尼力的力传感器;所述的曲柄连杆滑块机构用于实现减振器5的拉伸与 压缩行程往复运动，具有曲柄10,曲柄10的一端固定在架体上，曲柄10的另一端与曲轴连杆 8底端相连接，所述曲轴连杆8可绕着曲柄10的另一端点进行周向运动；曲轴连杆8的顶端与 下安装支座6连接从而带动减振器5随着曲轴连杆8作上下运动，曲轴连杆8上还设有与之配 合的滑块7,所述滑块7在曲轴连杆8的带动下进行上下往复运动;所述速度传感器9设置在 曲轴连杆8上用于测试由曲柄10转动产生的速度。所述力传感器固连于架体的上支座1上， 并与夹具4串联;所述力传感器包括应变片式力传感器2和压电式力传感器3,所述应变片式 力传感器2和压电式力传感器3串联。
[0091] 该测试系统，还包括减振器侧向力装置，所述减振器侧向力装置与减振器5固定连 接，具有测力环11，和固装在架体上的用于支撑钢索12的滚轮13,所述侧力环11的一端固定 连接在减振器5的中间部位，侧力环11的另一端与钢索12连接，钢索12通过滚轮13连接砝码 14;减振器随下面的曲柄连杆滑块机构作上下运动时，测力环11通过钢索12带动砝码14随 减振器5作小范围上下运动。
[0092] 工作原理，该系统加力于减振器5的活塞杆，使其往复运动。测取减振器示功特性 曲线特性采用正弦激励方式，使减振器活塞按简谐规律运动。
[0093] 该系统测试为:减振器运动机械能=阻尼力X位移;减振器消耗功率=阻尼力X 速度。
[0094] 实施例汽车悬架减振器能量特性计量方法
[0095] 基于上述汽车悬架减振器能量特性测试系统，包括以下步骤：
[0096] 1)建立减振器数学模型:建立减振器压缩行程和复原行程数学模型；
[0097] 根据减振器的工作机理，由于活塞杆相对工作缸运动速度的不同，活塞总成内的 复原阀片和底阀总成内的压缩阀片会有不同的开度，会对活塞往复运动时产生阻尼力有不 同的影响，因此需要分别建立减振器的各行程开阀前后的分段数学模型。
[0098]工作行程中活塞所受阻尼力为：
[0099] F = Pi(Ah-Ag)-P2 · Ah+fi (1)
[0100] 油液从上腔流入下腔的流量为Q:
[0101] Q= (Ah-Ag)vi (2)
[0102] 贮油腔油液补偿流入工作缸下腔的流量为Q。：
[0103] Qc = Ag · νι (3)
[0104] 式中:F为减振器受到的阻尼力，N;Ah为活塞面积，!112^为活塞杆面积，m2;V1为复原 行程活塞的运动速度，工作缸上腔的压强，Mpa;P 2工作缸下腔的压强，MpaJi为活 塞所受摩擦力，N。
[0105] 1.1)复原行程数学模型的建立
[0106] 1.1.1)开阀前数学模型
[0107] 假设减振器正常工作时油封密闭且无油液泄漏现象，减振器油只从工作缸的上腔 经过活塞上复原阀进入工作缸下腔，根据流体串并联原理，
[0108] 油液通过复原阀组件流入下腔的流量可表示为&:
[0110]贮油腔油液通过补偿流入工作缸下腔的流量为Qu:
[0112] (4)、（5)式中：Cd为油液流量系数;Am为阀体节流孔面积，m2;P为油液密度，Kg/m 3;An 为补偿阀节流孔面积，m2; P3为贮油腔压强，Mpa。
[0113]根据油液连续性方程得：
[0117] 1.1.2)开阀后数学模型
[0118] 随着复原行程中活塞速度的增加，工作缸上腔的压力也逐渐增加，当上腔油液压 力升高并且大于复原阀组件内阀片的承受压力时，复原阀片变形产生开阀状态，此时减振 器进入开阀后的阻尼状态。
[0119] 开阀后油液由工作缸上腔经过活塞复原阀上的常通孔和阀片发生形变形成的环 状缝隙进入下腔产生的流量为Q 2:
[0121]式中，A2为阀片发生形变产生的缝隙面积，m2。根据油液连续性方程得：
[0123] 则开阀后活塞所受阻尼力^可表示为：
[0125] 1.2)压缩行程数学模型的建立
[0126] 减振器压缩行程的数学模型与复原行程的数学模型建立方法相似同步骤1.1)，根 据活塞运动速度及阻尼力形成方式分别建立开阀前与开阀后减振器数学模型，具体为：
[0127] 1.2.1)压缩行程开阀前数学模型
[0128] 压缩阀开阀前，活塞处于低速运动状态，此时油液通过活塞压缩阀上的常通孔和 流通阀阀片环流间隙产生阻尼力。压缩阀阀片的常通孔节流及流通阀阀片环流间隙可认为 都属于薄壁小孔节流。
[0129] 由下腔通过流通阀阀片环流间隙流入上腔的流量Q3为：
[0131]由下腔通过压缩阀常通孔流入储油腔的流量为Q4:
[0133]根据油液连续性方程得：
[0135] 则阻尼力F3可表示为：
[0137] 式中:Ai为流通阀阀片发生形变产生的缝隙面积m2;A3为活塞上压缩阀常通孔截面 积m 2;
[0138] 1.2.2)压缩行程开阀后数学模型
[0139] 活塞处于高速运动状态时，减振器的压缩阀阀片变形形成环状缝隙节流，此时油 液通过活塞压缩阀上的常通孔和阀片变形形成的环状缝隙进入储油腔产生阻尼力。
[0140] 其流量为Q5:
[0142]根据油液连续性方程由得：
[0144] 则阻尼力F4可表示为：
[0146] 式中，Vy为复原行程活塞的运动速度，A4为压缩阀片发生形变产生的缝隙面积。
[0147] 2)建立减振器示功图面积数学模型
[0148] 所述减振器示功图面积是基于汽车悬架减振器能量特性测试系统，在一定的激振 行程及频率下，减振器的压缩行程和复原行程产生示功曲线，该曲线所围成的面积;基于减 振器数学模型建立该减振器示功图面积数学模型；
[0149] 减振器示功图面积如图1所示A1与A2的总和。其意义是表示减振器经过一个复原 与压缩行程，其产生的阻尼力衰减车身振动所做的功，其值的大小可表征减振器的能量效 能，既减振器衰减振动能量效率及能力。示功图面积分为上下两部分，可视为微小位移ds和 瞬时阻尼力h所构成矩形单元的面积累加之和。
[0150] 将减振器单个复原行程所用时间t均分为η等份，由于示功机的激振形式为正弦激 励信号，由速度传感器可测曲柄转动产生速度特性，即减振器活塞速度特性为：
[0151] v = 23TfRcos(23rft) =3rfDcos(23rft) (18)
[0152] 式中：f为输入信号的激振频率Hz ;R为曲柄长度mm;D为减振器总行程mm;
[0153] 减振器活塞位移特性由其速度特性表达式积分得到为：
[0155] 式中，S为减振器活塞位移，
[0156] 由减振器速度传感器可测得第i个时间点与第i+Ι个时间点的活塞运行速度，其表 达式分别为：
[0157] vi = 3TfDcos(23rfti) (20)
[0158] vi+i = 3TfDcos(23rfti+i) (21)
[0159] 对减振器速度表达式求积分可得到第i个时间点与第i + 1个时间点之间的微小位 移ds表达式为
[0161] 给定活塞一定的激振频率，在汽车悬架减振器能量特性测试系统上可采集到两个 相对应的数据，即活塞运动速度Vl及减振器瞬时阻尼力Fi表达式在步骤1)数学模型中已 推导;将第i点与第i + Ι点所对应的瞬时阻尼力的均值作为矩形单元的高，单个矩形的面积 可表示为：
[0162] Aii = Fi〇ds (23)
[0163] 式中:Ah为矩形单元的面积： 为平均阻尼力，N;
[0164] 因此，复原行程示功曲线围成的面积为：
[0166] 同理可求压缩行程其示功曲线围成的面积Ay，因此该减振器示功图面积可表示 为：
[0167] A=Ai+Ay (25)。
[0168] 3)示功图面积能量标定；
[0169] 如图2所示有畸变减振器示功图，其中近似椭圆形白色曲线为示功曲线轮廓线，即 目标曲线;黑色部分为背景区域。
[0170] 3.1)对图2所示减振器示功图进行图像处理，其主要技术路线包括图像读取、图像 灰度化及图像二值化；
[0171] 把减振器示功图看作具有不同灰度级的两类区域，利用目标曲线（白色)与其背景 在灰度特性上的差异，通过选取的阈值可界定图像中任一像素点的所属区域，其二值化后 图像如图3所示。通常，图像处理的输入与输出均为图像，图片格式无变化;而图像识别的不 同之处在于，其输出则是关于图片类别或结构分析数据。
[0172 ] 3.1.1)基于Canny边缘检测算法对目标曲线进行边缘检测，检测曲线是否平滑；
[0173]基于Canny边缘检测算法对目标曲线进行边缘检测，用范函求导方法推导出高斯 函数的一阶导数即为最优边缘检测算子的最佳近似。Canny算法可对目标曲线平滑处理，其 方法是采用高斯函数：
[0175] 式中，σ为高斯滤波器平滑程度控制参数。
[0176] 3.1.2)平滑的曲线使用一阶微分算子，获得各像素点处的梯度幅值和方向；
[0177] 平滑后示功曲线I各像素点处的梯度幅值的计算方法是使用一阶微分算子，这样 可获得各像素点处的梯度幅值和方向。点（i，j)处2个方向的偏导数G x(i，j)和Gy(i，j)分别 为：
[0178] Gx(i,j) = [I(i,j+l)-I(i,j)+I(i+l,j+l)-I(i+l,j)]/2 (27)
[0179] Gy(i,j) = [I(i,j)-I(i+l,j)+I(i,j+l)-I(i+l,j+l)]/2 (28)
[0180] 3.1.3)通过细化各像素点的的梯度幅值图像，可以精确识别示功曲线边缘线，非 极大值抑制即只需保留幅值的局部极大值；
[0181] 通过细化各像素点的的梯度幅值图像，可以精确识别示功曲线边缘线，非极大值 抑制即只需保留幅值的局部极大值[5]Xanny算子在梯度幅值图像G中以点（i，j)为中心3X 3的邻域内进行插值运算，若点（i，j)处的梯度幅值G(i，j)小于相邻的插值，那么该像素点 (i，j)为非边缘点，相反则可视为候选边缘点，从而得到候选边缘曲线图像。
[0182] 3.1.4)使用Canny算子从候选边缘点中检测和连接得到最终的边缘；
[0183] Canny算子采用双阈值法从候选边缘点中检测和连接得到最终的边缘。首先根据 图像特点选取两个合理的阈值，其中设定高阈值为240，低阈值为120，然后在候选边缘曲线 图像中对任一像素点（i，j)进行检测。如果被检测像素点（i，j)梯度幅值G(i，j)大于所选取 的高阈值240，可将其视为边缘点；相反如果其值小于低阈值120，那么判定其为非边缘点。 若梯度幅值介于高低阈值之间的像素点，可先将其视为疑似边缘点，接着根据边缘的连通 性对其进行判断。若该像素点的邻接像素中存在边缘点，则认为该像素点也为边缘点，否则 为非边缘点。
[0184] 使用Matleb平台编程计算结果求得示功曲线轮廓内包含142554个像素点，像素面 积为257。
[0185] 3.2)用LabVIEW对减振器示功图面积数学模型编程，之后导入减振器在汽车悬架 减振器能量特性测试系统中所采集的2000个采样点的位移及力值数据，便可求得减振器所 吸收的能量；
[0186] 具体为，依据示功图面积微单元累加法，使用LabVIEW编程并计算减振器上下往复 一个周期所吸收的能量总和。其程序框图如图4所示。
[0187] 导入减振器在汽车悬架减振器能量特性测试系统中所采集的2000个采样点的位 移及力值数据，便可求得减振器所吸收的能量。如图5所示为减振器能量计算结果，测试得 到其复原阻尼力为1406N，压缩阻尼力为-942N，使用LabVIEW平台计算减振器上下往复一个 周期衰减的能量总和为69.613J，因此单位像素面积所表示的能量为0.27J。
[0188] 上述方法可应用于减振器性能分析并完善其性能评价方法，可用来判断减振器畸 变及检验减振器能量吸收能力。本发明首先建立了汽车悬架减振器在复原与压缩行程的数 学模型及定义减振器示功图面积，然后基于Matleb平台将计算机图像识别技术应用于减振 器示功图像素点个数的检测与面积的求解;随后使用LabVIEW平台计算求解减振器在一个 工作行程中吸收总能量，从而定量分析出减振器示功图面积与所吸收能量的关系，为减振 器的性能检测提供新的评价方法及量化标准。
【主权项】
1. 汽车悬架减振器能量特性测试系统，其特征在于，包括架体，设置在架体上可产生正 弦运动的曲柄连杆滑块机构，测试时用于固定减振器(5)上吊环的夹具(4)和用于固定减振 器(5)下吊环的下安装支座(6)，速度传感器(9)，及用于测量阻尼力的力传感器;所述的曲 柄连杆滑块机构用于实现减振器(5)的拉伸与压缩行程往复运动，具有曲柄（10)，曲柄（10) 的一端固定在架体上，曲柄（10)的另一端与曲轴连杆(8)底端相连接，所述曲轴连杆(8)可 绕着曲柄（10)的另一端点进行周向运动；曲轴连杆(8)的顶端与下安装支座(6)连接从而带 动减振器(5)随着曲轴连杆(8)作上下运动，曲轴连杆(8)上还设有与之配合的滑块(7)，所 述滑块(7)在曲轴连杆(8)的带动下进行上下往复运动;所述速度传感器(9)设置在曲轴连 杆(8)上用于测试由曲柄(10)转动产生的速度。2. 如权利要求1所述的汽车悬架减振器能量特性测试系统，其特征在于，所述力传感器 固连于架体的上支座（1)上，并与夹具(4)串联;所述力传感器包括应变片式力传感器(2)和 压电式力传感器(3)，所述应变片式力传感器(2)和压电式力传感器(3)串联。3. 如权利要求1所述的汽车悬架减振器能量特性测试系统，其特征在于，还包括减振器 侧向力装置，所述减振器侧向力装置与减振器(5)固定连接，具有测力环（11)，和固装在架 体上的用于支撑钢索（12)的滚轮（13)，所述侧力环（11)的一端固定连接在减振器(5)的中 间部位，侧力环（11)的另一端与钢索（12)连接，钢索（12)通过滚轮(13)连接砝码（14);减振 器随下面的曲柄连杆滑块机构作上下运动时，测力环（11)通过钢索（12)带动砝码（14)随减 振器(5)作小范围上下运动。4. 汽车悬架减振器能量特性计量方法，应用权要求1-3所述的汽车悬架减振器能量特 性测试系统，其特征在于，包括以下步骤： 1) 建立减振器数学模型，包括建立减振器压缩行程和复原行程数学模型； 2) 建立减振器示功图面积数学模型，基于汽车悬架减振器能量特性测试系统，在激振 行程和频率下，减振器的压缩行程和复原行程产生示功曲线，该示功曲线组成的图形即减 振器示功图，该示功曲线所围成的面积即减振器示功图面积;基于减振器数学模型建立减 振器示功图面积数学模型； 3) 示功图面积能量标定， 3.1) 基于Matleb平台利用计算机图像识别技术对减振器示功图进行图像处理，包括图 像读取、图像灰度化及图像二值化，检测减振器示功图像素点个数并求解减振器示功图面 积； 3.2) 用LabVIEW对减振器示功图面积数学模型编程，并利用减振器在汽车悬架减振器 能量特性测试系统中所采集的数个采样点的速度值及阻尼力值数据，计算求解减振器在一 个工作行程中吸收的总能量。5. 如权利要求4所述的汽车悬架减振器能量特性计量方法，其特征在于，步骤1)建立减 振器数学模型是根据减振器工作机理，分别建立减振器的各行程开阀前后的分段数学模 型； 工作行程中活塞所受阻尼力为： F = Pi(Ah-Ag)-P2 · Ah+fi (1) 油液从上腔流入下腔的流量为Q: Q= (Ah-Ag)vi (2) 贮油腔油液补偿流入工作缸下腔的流量为Q。： Qc=Ag · vi (3) 式中：F为减振器受到的阻尼力N;Ah为活塞面积m2;Ag为活塞杆面积m 2;V1为复原行程活 塞的运动速度m/sTl·为工作缸上腔的压强Mpa;P2工作缸下腔的压强MpajiS活塞所受摩擦 力N; 1.1) 建立复原行程数学模型 1.1.1) 复原行程开阀前数学模型 假设减振器正常工作时油封密闭且无油液泄漏现象，减振器油只从工作缸的上腔经过 活塞上复原阀进入工作缸下腔，根据流体串并联原理， 油液通过复原阀组件流入下腔的流量&为：贮油腔油液通过补偿流入工作缸下腔的流量为Qcl: v y根据油液连续性方程由Q=Qi，Qc=Qci得：则开阀前活塞所受阻尼力^可表示为：式中：Cd为油液流量系数;Am为阀体节流孔面积m2 ;P为油液密度Kg/m3; An为补偿阀节流 孔面积m2; P3为贮油腔压强Mpa; 1.1.2) 复原行程开阀后数学模型 随着复原行程中活塞速度的增加，工作缸上腔的压力也逐渐增加，当上腔油液压力升 高并且大于复原阀组件内阀片的承受压力时，复原阀片变形产生开阀状态，此时减振器进 入开阀后的阻尼状态； 开阀后油液由工作缸上腔经过活塞复原阀上的常通孔和阀片发生形变形成的环状缝 隙进入下腔产生的流量为⑴·式中，A2为阀片发生形变产生的缝隙面积m2,根据油液连续性方程由〇=〇2得：则开阀后活塞所受阻尼力F2可表示为：1.2)压缩行程数学模型的建立 减振器压缩行程的数学模型与复原行程的数学模型建立方法同步骤1.1)，根据活塞运 动速度及阻尼力形成方式分别建立开阀前与开阀后减振器数学模型。6. 如权利要求5所述的汽车悬架减振器能量特性计量方法，其特征在于，步骤2)减振器 示功图面积分为上下两部分，分别为复原行程产生的示功曲线围成的上部分和减振器压缩 行程产生的示功曲线围成的下部分，可视为微小位移ds和瞬时阻尼力Fi所构成矩形单元的 面积累加之和； 将减振器单个复原行程所用时间t均分为η等份，由于示功机的激振形式为正弦激励信 号，由速度传感器可测曲柄转动产生速度特性，即减振器活塞速度特性为： v = 23TfRcos(23rft) = 3TfDcos(23rft) (11) 式中:f为输入信号的激振频率Hz ;R为曲柄长度mm;D为减振器总行程mm; 减振器活塞位移特性由其速度特性表达式积分得到为：式中，S为减振器活塞位移， 由减振器速度传感器可测得第i个时间点与第i+Ι个时间点的活塞运行速度，其表达式 分别为： Vi = 3TfDcos(23rfti) (13)vi+i = 3TfDcos(23rfti+i) (14) 对减振器速度表达式求积分可得到第i个时间点与第i+1个时间点之间的微小位移ds 表达式为给定活塞一定的激振频率，在汽车悬架减振器能量特性测试系统上可采集到两个相对 应的数据，即活塞运动速度Vi及减振器瞬时阻尼力FnFi表达式在步骤1)数学模型中已推 导;将第i点与第i + Ι点所对应的瞬时阻尼力的均值作为矩形单元的高，单个矩形的面积可 表示为： An = Fi〇ds (16) 式中:An为矩形单元的面积，为平均阻尼力，N; 因此，复原行程示功曲线围成的面积为：同理可求压缩行程其示功曲线围成的面积Ay，因此该减振器示功图面积可表示为： A=Ai+Ay (18) 〇7. 如权利要求4所述的汽车悬架减振器能量特性计量方法，其特征在于，步骤3.1)具体 包括如下步骤： 3.1.1)基于Canny边缘检测算法对减振器示功曲线进行边缘检测，检测该示功曲线是 否平滑，Canny算法可对目标曲线平滑处理，其方法是采用高斯函数：式中，σ为高斯滤波器平滑程度控制参数； 3.1.2) 平滑后的曲线使用一阶微分算子，获得各像素点处的梯度幅值和方向，点（i，j) 处2个方向的偏导数Gx( i，j)和Gy(i，j)分别为： Gx(i,j) = [I(i,j+l)-I(i,j)+I(i+l,j+l)-I(i+l,j)]/2 (20) Gy(i,j) = [I(i,j)-I(i+l,j)+I(i,j+l)-I(i+l,j+l)]/2 (21); 3.1.3) 通过细化各像素点的梯度幅值图像，可精确识别示功曲线边缘线，非极大值抑 制即只需保留幅值的局部极大值； 3.1.4) 使用Canny算子从候选边缘点中检测和连接得到最终的边缘。8. 如权利要求7所述的汽车悬架减振器能量特性计量方法，其特征在于，步骤3.1.3)具 体为Canny算子在梯度幅值图像G中以点（i，j)为中心3 X 3的邻域内进行插值运算，若点（i， j)处的梯度幅值G(i，j)小于相邻的插值，那么该像素点（i，j)为非边缘点，相反则可视为候 选边缘点，从而得到候选边缘曲线图像。9. 如权利要求7所述的汽车悬架减振器能量特性计量方法，其特征在于，步骤3.1.4)具 体为根据图像特点选取两个合理的阈值，分别为高阈值和低阈值，然后在候选边缘曲线图 像中对任一像素点（i，j)进行检测，如果被检测像素点（i，j)梯度幅值G(i，j)大于所选取的 高阈值，可将其视为边缘点；相反如果其值小于低阈值，那么判定其为非边缘点；若梯度幅 值介于高低阈值之间的像素点，可先将其视为疑似边缘点，接着根据边缘的连通性对其进 行判断，若该像素点的邻接像素中存在边缘点，则认为该像素点也为边缘点，否则为非边缘 点。
【文档编号】G01M17/04GK105865811SQ201610236338
【公开日】2016年8月17日
【申请日】2016年4月14日
【发明人】王天利, 田雪, 曾庆东, 王雪, 杨旭光
【申请人】辽宁工业大学