轮胎与路面的接触应力分布及应力集中的检测方法

摘要

本发明提供了一种轮胎与路面的接触应力分布及应力集中的检测方法，包括：获取压力敏感薄膜的密度显示图像；根据所述密度显示图像获取接触应力数据；计算轮胎与路面的有效接触面积；将轮胎与路面的接触应力划分成多个应力区间；计算各个所述应力区间的接触应力分布面积；计算各个所述应力区间的接触应力分布面积占所述有效接触面积的比例，获取所述应力区间内的接触应力分布累积概率；根据所述接触应力数据和所述接触应力分布累积概率计算韦伯分布函数中的位置参数、形状参数和尺度参数，获得应力分布模型；获取所述接触应力分布累积概率的特定分位值作为应力特征值；计算所述应力特征值与等效应力圆均布载荷的比值，得到应力集中因子。

说明书

技术领域

本发明涉及道路工程技术领域，特别是涉及一种轮胎与路面的接触应力分布及应力集中的检测方法。

背景技术

轮胎接地压力历来都是轮胎制造行业和交通部门所关注的焦点。初期，交通荷载相对较轻，加之轮胎尺寸及材料相对较软等原因，轮胎接地面积呈近似圆形，路面结构力学分析时，为简化计算模型，常常将轮胎接地荷载简化为均布的当量圆形荷载，且被广泛采用。随着整个社会的进步，无论是轮胎制造行业还是路面结构设计者都注意到，重载交通下，大尺寸轮胎、高轮压(超过0.7MPa)，轮胎花纹的改进甚至于轮胎材料的改善，都导致轮胎接地面积区、压力都明显区别于早期的认识。

近期的不少研究，将轮胎与路面接触界面假设为光滑平面，开展了轮胎接地面积、压力分布等内容的虚拟实验研究，认为接触面不呈圆状，接地压力分布不均匀的情况远超想象。国内相关研究显示，某载重车的轮胎在相同胎压下，变化荷载，接地面积的折减率横向花纹在85％左右，纵向花纹在80％左右，等等。然而，真实路面的宏观和微凸体构造对胎/路接触界面的应力分布有直接影响，尤其是应力集中现象，而路面简化为光滑平面显然不合乎实际情况。此外，受技术手段制约，轮胎与实际路面作用的接触压力分布规律及其不均匀程度，无法有效获取，更无法合理、准确地量化评价。目前的接触应力分布的描述仍然停留在接触界面应力代表值的宏观趋势描述和接地轮廓的几何分析上，而应力集中效应主要采用的是应力云图的主观判定方法。因此，如何提供一种轮胎与路面的接触应力分布及应力集中的检测方法以准确的反映轮胎与路面的接触应力分布规律及应力集中效应是本领域技术人员亟需解决的问题。

发明内容

本发明的目的是提供一种轮胎与路面的接触应力分布及应力集中的检测方法，以获取客观、合理且准确的轮胎与路面的接触应力分布规律及应力集中效应，并为轮胎制造行业和道路工程提供重要的设计依据。

为实现上述目的，本发明提供了一种轮胎与路面的接触应力分布及应力集中的检测方法，包括：

获取实验后的轮胎与路面之间的压力敏感薄膜的密度显示图像；

根据所述密度显示图像获取接触应力数据；

根据所述接触应力数据计算轮胎与路面的有效接触面积；

根据所述压力敏感薄膜的量程，将轮胎与路面的接触应力划分成多个应力区间；

根据所述应力区间和所述接触应力数据，计算各个所述应力区间的接触应力分布面积；

计算各个所述应力区间的接触应力分布面积占所述有效接触面积的比例，获取所述应力区间内的接触应力分布累积概率；

根据所述接触应力数据和所述接触应力分布累积概率计算韦伯分布函数中的位置参数、形状参数和尺度参数，获得应力分布模型；

获取所述接触应力分布累积概率的特定分位值作为应力特征值；

计算所述应力特征值与等效应力圆均布载荷的比值，得到应力集中因子。

可选的，所述根据所述密度显示图像获取接触应力数据，具体包括：

扫描所述压力敏感薄膜的密度显示图像；

构建压力敏感薄膜空间坐标系；

获取所述压力敏感薄膜空间坐标系中各个坐标点(x，y)对应的幅值f(x，y)，得到所述坐标点(x，y)的接触应力值；

在横坐标x取值范围0～M和纵坐标y取值范围0～N内获取多个接触应力值，确定接触应力数据为接触应力矩阵：

可选的，所述根据所述接触应力数据计算轮胎与路面的有效接触面积，具体包括：

获取所述幅值f(x，y)非零的坐标点；

根据所述非零的坐标点，确定多个对应的非零区域；

计算多个所述非零区域的面积，得到轮胎与路面的有效接触面积。

可选的，多个所述应力区间分别为0MPa～0.1MPa，0.1MPa～0.2MPa，……，(k-1)/10MPa～k/10MPa，其中k为应力区间的个数。

可选的，所述根据所述接触应力数据和所述接触应力分布累积概率计算韦伯分布函数中的位置参数、形状参数和尺度参数，获得应力分布模型，具体包括：

确定韦伯分布函数为：其中，a为位置参数；c为形状参数；1/b为尺度参数，属于疏散性指标；

根据所述接触应力数据和所述接触应力分布累积概率采用数值解法计算所述韦伯分布函数中的位置参数a、形状参数c和尺度参数1/b。

可选的，所述根据所述接触应力数据和所述接触应力分布累积概率采用数值解法计算所述韦伯分布函数中的位置参数a、形状参数c和尺度参数1/b，具体包括：

获取接触应力数值x；

按由小到大的顺序排列所述接触应力值x，获得接触应力值序列{x_i}，其中i＝1,2……，n，且x₁＜x₂＜…＜x_n，x_i表示每个坐标点的接触应力值，n为有效接触面积范围内的接触应力值数量；

根据所述接触应力值序列{x_i}，按照公式i＝1,2,…,n和公式i＝1,2,…,n计算，获得序列{ξ_i}；其中，G_n(x_i)为接触应力分布累积概率，所述接触应力分布累积概率为不小于给定接触应力值x_i的接触应力值出现概率；

根据序列{ξ_i}和按特定步长给定的设定位置参数值a_m按公式η_i＝ln(x_i-a),i＝1,2,…,n计算序列{η_i}_m，并将设定位置参数值a_m，序列{η_i}_m和序列{ξ_i}代入公式计算残差平方和Q_m；

计算相邻两个所述设定位置参数值a_m对应的残差平方和Q_m之差的绝对值；

判断所述绝对值是否小于给定的精度ε；

若是，确定所述给定的设定位置参数值a_m为位置参数a；

根据所述位置参数a，利用最小二乘法计算形状参数c和尺度参数1/b。

可选的，所述特定分位值为95％的分位值。

可选的，所述等效应力圆均布载荷为《公路沥青路面设计规范》规定的等效应力圆均布荷载0.7MPa。

可选的，所述轮胎与路面之间的压力敏感薄膜的密度显示图像包括同一单轮负荷的轮胎与多个实验路面之间的压力敏感薄膜的密度显示图像；所述实验路面包括无构造光滑路面、细型沥青混凝土路面和粗型沥青混凝土路面。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

本发明提供的轮胎与路面的接触应力分布及应力集中的检测方法利用压力敏感薄膜的密度显示图像获取轮胎与路面的接触应力数据，并根据该接触应力数据获取应力分布模型以及应力集中因子，从而通过应力分布模型来反映轮胎与路面的接触应力分布规律，通过应力集中因子反映轮胎与路面的接触应力集中效应。由于本发明的检测方法中是通过数学模型和数值分别反映接触应力的分布规律和集中效应，能够更直观、更客观的反映轮胎与路面的接触应力分布规律和集中效应；并且，本发明是直接从设置于轮胎与路面之间的压力敏感薄膜上获取接触应力数据，那么在检测不同路面与轮胎的接触应力时，分别获取不同路面与轮胎之间的压力敏感薄膜即可，检测方法简单，可操作性强，而且压力敏感薄膜不会损坏路面或轮胎的构造，也就能准确的反映轮胎与路面的接触应力，因此，也提高了本检测方法的准确性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提供的轮胎与路面的接触应力分布及应力集中的检测方法的流程图；

图2为轮胎与无构造路面之间的压力敏感薄膜的密度显示图像；

图3为图2数值化后的图像；

图4为轮胎与细型沥青混凝土路面之间的压力敏感薄膜的密度显示图像；

图5为图4数值化后的图像；

图6为轮胎与无构造路面的接触应力分布累积概率与应力分布模型的拟合图；

图7为实验测得的轮胎与无构造路面的接触应力分布累积概率与应力分布模型计算的接触应力分布累积概率的曲线图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的目的是提供一种轮胎与路面的接触应力分布及应力集中的检测方法，以获取客观、合理且准确的轮胎与路面的接触应力分布规律及应力集中效应，并为轮胎制造行业和道路工程提供重要的设计依据。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1为本发明提供的轮胎与路面的接触应力分布及应力集中的检测方法的流程图。如图1所示，

轮胎与路面的接触应力分布及应力集中的检测方法，包括：

步骤101：获取实验后的轮胎与路面之间的压力敏感薄膜的密度显示图像；

步骤102：根据所述密度显示图像获取接触应力数据；

步骤103：根据所述接触应力数据计算轮胎与路面的有效接触面积；

步骤104：根据所述压力敏感薄膜的量程，将轮胎与路面的接触应力划分成多个应力区间；具体的，多个所述应力区间分别为0MPa～0.1MPa，0.1MPa～0.2MPa，……，(k-1)/10MPa～k/10MPa，其中k为应力区间的个数；

步骤105：根据所述应力区间和所述接触应力数据，计算各个所述应力区间的接触应力分布面积；

步骤106：计算各个所述应力区间的接触应力分布面积占所述有效接触面积的比例，获取所述应力区间内的接触应力分布累积概率；

步骤107：根据所述接触应力数据和所述接触应力分布累积概率计算韦伯分布函数中的位置参数、形状参数和尺度参数，获得应力分布模型；

步骤108：获取所述接触应力分布累积概率的特定分位值作为应力特征值；其中，特定分位值可以为95％的分位值；

步骤109：计算所述应力特征值与等效应力圆均布载荷的比值，得到应力集中因子；其中，等效应力圆均布载荷可以为《公路沥青路面设计规范》规定的等效应力圆均布荷载0.7MPa；或是轮胎与无构造路面之间的平均基础压力值。

本实施例通过应力分布模型来反映轮胎与路面的接触应力分布规律，通过应力集中因子反映轮胎与路面的接触应力集中效应，即通过数学模型和数值分别反映接触应力的分布规律和集中效应，该方法能够更直观、更客观的反映轮胎与路面的接触应力分布规律和集中效应；并且，本发明是直接从设置于轮胎与路面之间的压力敏感薄膜上获取接触应力数据，那么在检测不同路面与轮胎的接触应力时，分别获取不同路面与轮胎之间的压力敏感薄膜即可，检测方法简单，可操作性强，而且压力敏感薄膜不会损坏路面或轮胎的构造，也就能准确的反映轮胎与路面的接触应力，因此，也提高了本检测方法的准确性。

具体的，对于步骤101，获取实验后的轮胎与路面之间的压力敏感薄膜的密度显示图像的方法，具体包括以下步骤：

选取实验路面；

在轮胎与实验路面之间放置至少3个量程的压力敏感薄膜；

施加静态单轮负荷至少2min后，卸压；

取出压力敏感薄膜并放置于干燥无紫外线环境中3～5h，获得敏度密度显示图像，如图2和4所示。

其中，实验路面可以为无构造路面，如钢板；或细型沥青混凝土路面；或粗型沥青混凝土路面。通过选取不同的实验路面可以获得不同的压力敏感薄膜的密度显示图像，从而在相同施加压力和单轮负荷的情况下获得轮胎与不同路面的接触应力数据。进而实现对不同构造的路面与轮胎接触应力的分布规律和集中效应进行评价。本实施例中，压力敏感薄膜为双片薄膜结构，其厚度为0.09mm×2，这样将压力敏感薄膜放置在轮胎与路面之间，不会改变轮胎与路面的接触特性。该双片薄膜结构主要由两个聚酯片基组成，一个涂有微囊生色物质，而另一个涂有显色物质，轮胎与路面的接触压力水平可以通过色彩变化的密度显示出来，其测量结果精度高，可达0.016mm²，受环境以及测量因素的影响较小。为了保证接触应力的准确，压力敏感薄膜的静载受压时间至少2min。

对于步骤102，根据密度显示图像获取接触应力数据的方法，具体包括以下步骤：

扫描所述压力敏感薄膜的密度显示图像；

构建压力敏感薄膜空间坐标系；

获取所述压力敏感薄膜空间坐标系中各个坐标点(x，y)对应的幅值f(x，y)，得到所述坐标点(x，y)的接触应力值；

在横坐标x取值范围0～M和纵坐标y取值范围0～N内获取多个接触应力值，确定接触应力数据为接触应力矩阵：

其中，在扫描密度显示图像之前需要对扫描设备进行校正，以确保扫描后结果的准确性。经过上述步骤，将密度显示图像数值化后的图分别如图3和5所示，图中，压力敏感薄膜空间坐标系中的每一点的显示的颜色所对于的数值为坐标系右侧的应力颜色比色条的中相同颜色对应的压力数值，即为接触应力值，由该压力敏感薄膜空间坐标系即可获得接触应力值，进而获得接触应力数据。并且该获得方法能够直观的显示出接触应力的变化及分布情况，能够准确的获得接触应力值。

对于步骤103，根据所述接触应力数据计算轮胎与路面的有效接触面积的方法，具体包括如下步骤：

获取所述幅值f(x，y)非零的坐标点；

根据所述非零的坐标点，确定多个对应的非零区域；

计算多个所述非零区域的面积，得到轮胎与路面的有效接触面积。

其中，每一坐标点所占区域面积的压力敏感薄膜的精度，即0.016mm²。因此，在本实施方式中，获取幅值f(x，y)非零的坐标点的个数即可获得有效接触面积。其他实施方式可以利用获取坐标点圈定的区域，在计算该圈定区域的面积即可获得有效接触面积。相比这种圈定区域的方法，本实施方式中的根据压力敏感薄膜的精度确定点的面积，进而求得的有效接触面积更准确。

对于步骤107，根据接触应力数据和接触应力分布累积概率计算韦伯分布函数中的位置参数、形状参数和尺度参数，获得应力分布模型的方法，具体包括如下步骤：

确定韦伯分布函数为：

$F\left(x\right)=\left(\begin{array}{cc}1-\exp \{-{\left(\frac{x-a}{b}\right)}^c\},& x\ge a\\ 0,& x<a\end{array}\right),(b,c>0)---\left(1\right),$

其中，a为位置参数；c为形状参数；1/b为尺度参数，属于疏散性指标；

对式(2)的等式两边做对数变换，得到

$\ln \left[\frac{1}{G\left(x\right)}\right]={\left(\frac{x-a}{b}\right)}^c---\left(2\right),$

其中，G(x)为随机变量接触应力值x的接触应力分布累积概率，即

G(x)＝1-F(x) (3)；

为了方便计算，假设参数ξ，η，和B，令