CN113742838B - 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质 - Google Patents

一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质 Download PDF

Info

Publication number
CN113742838B
CN113742838B CN202110789761.XA CN202110789761A CN113742838B CN 113742838 B CN113742838 B CN 113742838B CN 202110789761 A CN202110789761 A CN 202110789761A CN 113742838 B CN113742838 B CN 113742838B
Authority
CN
China
Prior art keywords
test data
tire
longitudinal force
test
abscissa
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110789761.XA
Other languages
English (en)
Other versions
CN113742838A (zh
Inventor
夏丹华
吕剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongce Rubber Group Co Ltd
Original Assignee
Zhongce Rubber Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhongce Rubber Group Co Ltd filed Critical Zhongce Rubber Group Co Ltd
Priority to CN202110789761.XA priority Critical patent/CN113742838B/zh
Publication of CN113742838A publication Critical patent/CN113742838A/zh
Application granted granted Critical
Publication of CN113742838B publication Critical patent/CN113742838B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Operations Research (AREA)
  • Tires In General (AREA)

Abstract

本发明属于轿车轮胎力学特性数据处理领域，具体涉及一种针对瞬态复合工况下的轮胎纵向力关于滑移率试验数据进行分区拟合的方法、智能设备和计算机可读载体介质。本发明将滑移率测量范围分为七个区域，对每个区间的纵向力关于滑移率的试验数据进行三次多项式拟合，且约束每个节点处的函数值和一阶导数值相等。通过分区拟合，可将加载段和卸载段的试验数据进行高精度拟合成光滑的曲线，用以计算轮胎特征值、对标分析或用于轮胎模型辨识。具有准确性高、重复性好、效率高的特点。

Description

一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载 体介质
技术领域
本发明属于轿车轮胎力学特性数据处理领域，具体涉及一种针对瞬态复合工况（侧倾、侧偏及纵滑共同作用）下的轮胎纵向力关于滑移率试验数据进行分区拟合的方法、智能设备和计算机可读载体介质。
背景技术
轮胎力和力矩测试数据主要用于对标分析、品控及动力学建模，是轮胎产品性能研发的重要技术手段。轮胎力和力矩测试方法分为稳态测试法和瞬态测试法，分别用于测评轮胎稳态性能和瞬态性能。
为了消除轮胎不均匀性的影响，采用瞬态测试法试验时通常设置较高的数据采样频率，如60km/h条件下常设置为250Hz以上。因此，这种情况下采用瞬态测试法测得的轮胎纵向力关于滑移率的加载与卸载两部分的试验数据可能达到几千个数据点甚至更多。因此，开发有效的数据处理方法对于轮胎瞬态力学性能的准确计算与评估至关重要。
目前，瞬态试验数据的处理方法主要有多项式、光滑样条或魔术公式的整体拟合。采样多项式对试验数据进行整体拟合的精度较差，光滑样条拟合结果的光滑度和精度易受试验数据影响，魔术公式的拟合结果比前两者好，但是精度仍有待提升。
综上所述，为准确地对瞬态复合工况的纵向力关于滑移率试验数据进行拟合处理，需开发一种准确性高、重复性好的数据处理方法。
发明内容
为了解决上述的技术问题，本发明的第一个目的是提供一种瞬态复合工况的轮胎纵向力分区拟合方法，该方法通过分区拟合，可将加载段和卸载段的试验数据进行高精度拟合成光滑的曲线，用以计算轮胎特征值、对标分析或用于轮胎模型辨识，具有准确性高、重复性好、效率高的特点。
为了实现上述的目的，本发明采用了以下的技术方案：
一种瞬态复合工况的轮胎纵向力分区拟合方法，所述的复合工况为侧倾、侧偏及纵滑共同作用，该方法包括以下的步骤：
一、准备试验所需的轮胎作为试验胎；
二、将试验胎安装至合适的轮辋上并充气调整至所需的试验气压，根据试验要求进行静置停放合适的时间；
三、将充气停放后的轮胎安装至具有瞬态侧倾、侧偏及纵滑测试功能的轮胎六分力试验机上，设置试验所需的轮胎气压、路面速度、垂向载荷、外倾角、侧偏角与滑移率，进行分区测试或者整体扫掠测试；
四、提取瞬态复合工况下的轮胎纵向力关于滑移率的试验数据，并根据垂向载荷的命令值与反馈值，对纵向力数据进行修正；
五、按照以下方式对试验数据进行处理：
1）以滑移率为横坐标X轴，纵向力为纵坐标Y轴，建立直角坐标系；
2）在直角坐标系中绘制纵向力关于滑移率的试验数据散点图；
3）在横坐标中设置连接点；
4）根据所设置的连接点，将纵向力关于滑移率的试验数据散点图进行分区：
5）分别对每个区的纵向力关于滑移率试验数据进行单独的三次多项式拟合，并约束节点处的函数值相等和一阶导数值相等；
六、根据实际需求设置横坐标侧偏角序列；
七、将拟合得到的多项式系数分别回代至对应的三次多项式中，获得拟合后的纵向力关于滑移率的数据结果。
优选，所述步骤四根据F Xtarget =(F Ztarget /F Zmeasured )*F Xmeasured 对试验数据进行修正，其中F Ztarget 为目标垂向载荷，F Zmeasured 为测得的垂向载荷，F Xmeasured 为测得的纵向力，F Xtarget 为修正后的纵向力。
优选，所述步骤2）如果是采用分区测试，则需对试验数据进行拼接处理。
本发明中步骤4）中的分区可以采用5-10个分区，优选是7个分区；再优选，所述步骤3）连接点依次设为-15、-7、-1.5、1.5、7、15，滑移率测量范围分为七个区域，步骤4）分区方法如下：
4.1）将横坐标SR min 至-15对应的试验数据标记为一区，对应的三次多项式方程记为F x1 =a1 κ 1 3+b1 κ 1 2+c1 κ 1 +d1；
4.2）将横坐标-15至-7对应的试验数据标记为二区，对应的三次多项式方程记为F x2 =a2 κ 2 3+b2 κ 2 2+c2 κ 2 +d2；
4.3）将横坐标-7至-1.5对应的试验数据标记为三区，对应的三次多项式方程记为F x3 =a3 κ 3 3+b3 κ 3 2+c3 κ 3 +d3；
4.4）将横坐标-1.5至+1.5对应的试验数据标记为四区，对应的三次多项式方程记为F x4 =a4 κ 4 3+b4 κ 4 2+c4 κ 4 +d4；
4.5）将横坐标+1.5至+7对应的试验数据标记为五区，对应的三次多项式方程记为F x5 =a5 κ 5 3+b5 κ 5 2+c5 κ 5 +d5；
4.6）将横坐标+7至+15对应的试验数据标记为六区，对应的三次多项式方程记为F x6 =a6 κ 6 3+b6 κ 6 2+c6 κ 6 +d6；
4.7）将横坐标+15至SR max 对应的试验数据标记为七区，对应的三次多项式方程记为F x7 =a7 κ 7 3+b7 κ 7 2+c7 κ 7 +d7。
再优选，所述步骤5）中三次多项式拟合如下：F x1 (-15) = F x2 (-15)且dF x1 (-15) = dF x2 (-15)、F x2 (-7) = F x3 (-7)且dF x2 (-7) = dF x3 (-7)、、F x3 (-1.5) = F x4 (-1.5)且dF x3 (- 1.5) = dF x4 (-1.5)、F x4 (1.5) = F x5 (1.5)且dF x4 (1.5) = dF x5 (1.5)、F x5 (7) = F x6 (7)且dF x5 (7) = dF x6 (7)、F x6 (15) = F x7 (15)且dF x6 (15) = dF x7 (15)。
优选，所述Samin为--25~-40区间，SAmax为25~40区间。
优选，所述步骤六-SR min ~SR max 的步进设为0.1。
进一步，本发明还公开了所述方法或的数据结果用以计算轮胎特征值、对标分析和/或用于轮胎模型辨识。
进一步，本发明还公开了一种智能设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现所述方法中步骤五至步骤七。
进一步，本发明还公开了一种存储程序指令的非暂时性计算机可读载体介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被所述处理器执行时实现所述方法中步骤五至步骤七。
本发明的有益效果
（1）采用本发明的拟合结果准确性高、可重复性好；
（2）本发明的分区拟合方法，各个区间的数据拟合互不干扰；
（3）本发明的数据处理方法效率高。
附图说明
图1为 瞬态复合工况下的轮胎纵向力关于滑移率试验数据处理流程图。
图2为 瞬态复合工况下的轮胎纵向力关于滑移率试验数据分区示意图。
图3、图4为 瞬态复合工况下的轮胎纵向力关于滑移率试验数据处理结果。
具体实施方式
采用本发明对225/60R18半钢子午线轮胎进行瞬态复合工况的力学特性测量，以下是具体的实施例过程。
1、抽取1条225/60R18规格的半钢子午线轮胎作为试验胎；
2、将试验轮胎安装至6.5J×18的轮辋上并充气调整至250kPa的气压，然后静置于实验室环境（室温为24±3℃）下停放3小时；
3、将充气停放后的试验胎安装至MTS flat-trac CT型号试验机上，设置合适的胎压、路面速度、垂向载荷、侧倾角、侧偏角及滑移率（SR min 取-30°、SR max 取+30°）对轮胎进行预实验；
4、预实验结束后，根据实际需要采用分区测试或整体扫掠测试，提取瞬态复合工况下的轮胎纵向力关于滑移率的试验数据；
5、根据F Xtarget =(F Ztarget /F Zmeasured )*F Xmeasured 对试验数据进行修正，其中F Ztarget 为目标垂向载荷，F Zmeasured 为测得的垂向载荷，F Xmeasured 为测得的纵向力，F Xtarget 为修正后的纵向力；
6、按照以下方式对试验数据进行处理：
6.1 对试验数据进行分区，将横坐标的滑移率值落在-30至-15范围的纵向力关于滑移率数据划分为一区，将-15至-7对应的试验数据划分为二区，将-7至-1.5对应的试验数据划分为三区，将-1.5至1.5对应的试验数据划分为四区，将1.5至7对应的试验数据划分为五区，将7至15对应的试验数据划分为六区，将15至30对应的试验数据划分为七区；
6.2 分别对一区至七区的纵向力关于滑移率的试验数据进行单独的三次多项式拟合，记F xi = ai κ i 3+bi κ i 2+ci κ i +di，其中i=1:7；
6.3 拟合过程中须保证：F x1 (-15) = F x2 (-15)且dF x1 (-15) = dF x2 (-15)、F x2 (-7) = F x3 (-7)且dF x2 (-7) = dF x3 (-7)、F x3 (-1.5) = F x4 (-1.5)且dF x3 (-1.5) = dF x4 (-1.5)、F x4 (1.5) = F x5 (1.5)且dF x4 (1.5) = dF x5 (1.5)、F x5 (7) = F x6 (7)且dF x5 (7) = dF x6 (7)、F x6 (15) = F x7 (15)且dF x6 (15) = dF x7 (15)；
6.4 拟合得到7组三次多项式的系数值，即ai、bi、ci、di（i=1:7）；
6.5 设置滑移率序列为-30:0.1:30，即步长为0.1°；
6.6 将得到的7组三次多项式系数分别回代至对应的三次多项式中，获得纵向力关于滑移率的曲线数据结果。以上为对本发明实施例的描述，通过对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的。本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施列，而是要符合与本文所公开的原理和新颖点相一致的最宽的范围。

Claims (7)

1.一种瞬态复合工况的轮胎纵向力分区拟合方法，所述的复合工况为侧倾、侧偏及纵滑共同作用，其特征在于，该方法包括以下的步骤：
一、准备试验所需的轮胎作为试验胎；
二、将试验胎安装至合适的轮辋上并充气调整至所需的试验气压，根据试验要求进行静置停放合适的时间；
三、将充气停放后的轮胎安装至具有瞬态侧倾、侧偏及纵滑测试功能的轮胎六分力试验机上，设置试验所需的轮胎气压、路面速度、垂向载荷、外倾角、侧偏角与滑移率，进行分区测试或者整体扫掠测试；
四、提取瞬态复合工况下的轮胎纵向力关于滑移率的试验数据，并根据垂向载荷的命令值与反馈值，对纵向力数据进行修正；
根据F Xtarget =(F Ztarget /F Zmeasured )×F Xmeasured 对试验数据进行修正，其中F Ztarget 为目标垂向载荷，F Zmeasured 为测得的垂向载荷，F Xmeasured 为测得的纵向力，F Xtarget 为修正后的纵向力;
五、按照以下方式对试验数据进行处理：
1）以滑移率为横坐标X轴，纵向力为纵坐标Y轴，建立直角坐标系；
2）在直角坐标系中绘制纵向力关于滑移率的试验数据散点图；
3）在横坐标中设置连接点；连接点依次设为-15、-7、-1.5、1.5、7、15，滑移率测量范围分为七个区域;
4）根据所设置的连接点，将纵向力关于滑移率的试验数据散点图进行分区：
4.1）将横坐标SR min 至-15对应的试验数据标记为一区，对应的三次多项式方程记为F x1 =a1 κ 1 3+b1 κ 1 2+c1 κ 1 +d1；
4.2）将横坐标-15至-7对应的试验数据标记为二区，对应的三次多项式方程记为F x2 =a2 κ 2 3+b2 κ 2 2+c2 κ 2 +d2；
4.3）将横坐标-7至-1.5对应的试验数据标记为三区，对应的三次多项式方程记为F x3 =a3 κ 3 3+b3 κ 3 2+c3 κ 3 +d3；
4.4）将横坐标-1.5至+1.5对应的试验数据标记为四区，对应的三次多项式方程记为F x4 =a4 κ 4 3+b4 κ 4 2+c4 κ 4 +d4；
4.5）将横坐标+1.5至+7对应的试验数据标记为五区，对应的三次多项式方程记为F x5 =a5 κ 5 3+b5 κ 5 2+c5 κ 5 +d5；
4.6）将横坐标+7至+15对应的试验数据标记为六区，对应的三次多项式方程记为F x6 =a6 κ 6 3+b6 κ 6 2+c6 κ 6 +d6；
4.7）将横坐标+15至SR max 对应的试验数据标记为七区，对应的三次多项式方程记为F x7 =a7 κ 7 3+b7 κ 7 2+c7 κ 7 +d7；
5）分别对每个区的纵向力关于滑移率试验数据进行单独的三次多项式拟合，并约束节点处的函数值相等和一阶导数值相等；
三次多项式拟合如下：F x1 (-15) = F x2 (-15)且dF x1 (-15) = dF x2 (-15)、F x2 (-7) = F x3 (- 7)且dF x2 (-7) = dF x3 (-7)、F x3 (-1.5) = F x4 (-1.5)且dF x3 (-1.5) = dF x4 (-1.5)、F x4 (1.5) = F x5 (1.5)且dF x4 (1.5) = dF x5 (1.5)、F x5 (7) = F x6 (7)且dF x5 (7) = dF x6 (7)、F x6 (15) = F x7 (15)且dF x6 (15) = dF x7 (15);
六、根据实际需求设置横坐标侧偏角序列；
七、将拟合得到的多项式系数分别回代至对应的三次多项式中，获得拟合后的纵向力关于滑移率的数据结果。
2.根据权利要求1所述的方法，其特征在于，步骤2）如果是采用分区测试，则需对试验数据进行拼接处理。
3.根据权利要求1所述的方法，其特征在于，SR min 为-25~-40区间，SR max 为25~40区间。
4.根据权利要求1或3所述的方法，其特征在于，步骤六-SR min ~SR max 的步进设为0.1。
5.权利要求1-4任意一项权利要求所述方法的数据结果用以计算轮胎特征值、对标分析和/或用于轮胎模型辨识。
6.一种智能设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述计算机程序被所述处理器执行时实现权利要求1-4任意一项权利要求所述方法中步骤五至步骤七。
7.一种存储程序指令的非暂时性计算机可读载体介质，其特征在于，所述计算机可读载体介质上存储计算机程序，所述计算机程序被处理器执行时实现权利要求1-4任意一项权利要求所述方法中步骤五至步骤七。
CN202110789761.XA 2021-07-13 2021-07-13 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质 Active CN113742838B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110789761.XA CN113742838B (zh) 2021-07-13 2021-07-13 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110789761.XA CN113742838B (zh) 2021-07-13 2021-07-13 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质

Publications (2)

Publication Number Publication Date
CN113742838A CN113742838A (zh) 2021-12-03
CN113742838B true CN113742838B (zh) 2023-09-26

Family

ID=78728607

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110789761.XA Active CN113742838B (zh) 2021-07-13 2021-07-13 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质

Country Status (1)

Country Link
CN (1) CN113742838B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107117072A (zh) * 2017-05-06 2017-09-01 哈尔滨工业大学 轮毂/轮边电机驱动电动汽车横摆稳定控制中的期望横摆率设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1840552B1 (en) * 2006-03-31 2012-05-02 The Yokohama Rubber Co., Ltd. Tire transient response data calculating method, tire designing method, and vehicle motion predicting method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107117072A (zh) * 2017-05-06 2017-09-01 哈尔滨工业大学 轮毂/轮边电机驱动电动汽车横摆稳定控制中的期望横摆率设计方法

Also Published As

Publication number Publication date
CN113742838A (zh) 2021-12-03

Similar Documents

Publication Publication Date Title
EP2655103B1 (en) Method and system for estimating the inflation pressure of a tire
US20120197548A1 (en) Tire rubber index calculating method, device, and computer-readable storage medium
EP3109588B1 (en) Method for measuring tread radius of tire, and deivce for measuirng tread radius used therefor
CN113742838B (zh) 一种瞬态复合工况轮胎纵向力分区拟合方法、设备和可读载体介质
KR101829698B1 (ko) 자동차 타이어의 롤링 반경을 시뮬레이션하는 방법
CN113761472B (zh) 一种瞬态纯纵滑工况轮胎纵向力分区拟合方法、设备和可读载体介质
CN113761471B (zh) 一种瞬态复合工况轮胎侧向力分区拟合方法、设备和可读载体介质
CN113761474B (zh) 一种瞬态复合工况轮胎回正力矩分区拟合方法、设备和可读载体介质
JP5515779B2 (ja) タイヤ接地面が路面から受ける物理量を予測する方法、タイヤ摩耗の予測方法、タイヤ摩耗の予測装置、および、プログラム
CN113761473B (zh) 瞬态纯侧偏工况轮胎回正力矩分区拟合方法、设备和可读载体介质
CN113742837B (zh) 一种瞬态纯侧偏工况轮胎侧向力分区拟合方法、设备和可读载体介质
JP2016528099A5 (zh)
EP2979075B1 (en) Tire uniformity improvement using estimates based on convolution/deconvolution with measured lateral force variation
CN115219247A (zh) 轮胎侧偏刚度的测试及辨识方法
US10377194B2 (en) Method of modeling tire performance
CN109708903A (zh) 一种轮胎滚动阻力系数测量方法、系统和装置
CN109996669B (zh) 通过使用微波监视固化状态
KR101829695B1 (ko) 사전 결정된 속도의 주행 조건에서 타이어를 모델링하는 방법
JP2004037286A (ja) 室内タイヤ耐久試験方法
CN113553656B (zh) 轮胎复合滑移特性测量及数据处理方法、设备和计算机可读载体介质
CN113553657B (zh) 轮胎侧偏特性的测量及数据处理方法、设备和计算机可读载体介质
CN113553658B (zh) 轮胎纵滑特性测量及数据处理方法、设备和计算机可读载体介质
KR101934690B1 (ko) 콘볼루션/디콘볼루션에 기초한 추정을 사용한 타이어 균일성 개선
JP2017150897A (ja) コーナリング性能評価方法
CN109668577B (zh) 一种用于移动机器人的轮子半径自动校准方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 310008 No. 1, No. 1 Street, Qiantang District, Hangzhou, Zhejiang

Applicant after: Zhongce Rubber Group Co.,Ltd.

Address before: 310008 No. 1, Baiyang street, Qiantang District, Hangzhou City, Zhejiang Province

Applicant before: ZHONGCE RUBBER GROUP Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant