CN102472691A - 预测轮胎运行噪音的方法 - Google Patents

预测轮胎运行噪音的方法 Download PDF

Info

Publication number
CN102472691A
CN102472691A CN2010800331409A CN201080033140A CN102472691A CN 102472691 A CN102472691 A CN 102472691A CN 2010800331409 A CN2010800331409 A CN 2010800331409A CN 201080033140 A CN201080033140 A CN 201080033140A CN 102472691 A CN102472691 A CN 102472691A
Authority
CN
China
Prior art keywords
pressure head
value
descriptor
road surface
tire
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.)
Granted
Application number
CN2010800331409A
Other languages
English (en)
Other versions
CN102472691B (zh
Inventor
D·罗舒
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.)
Compagnie Generale des Etablissements Michelin SCA
Original Assignee
Michelin Recherche et Technique SA Switzerland
Societe de Technologie Michelin SAS
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 Michelin Recherche et Technique SA Switzerland, Societe de Technologie Michelin SAS filed Critical Michelin Recherche et Technique SA Switzerland
Publication of CN102472691A publication Critical patent/CN102472691A/zh
Application granted granted Critical
Publication of CN102472691B publication Critical patent/CN102472691B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60CVEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
    • B60C99/00Subject matter not provided for in other groups of this subclass
    • B60C99/006Computer aided tyre design or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles
    • G01M17/02Tyres

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Tires In General (AREA)

Abstract

在预测轮胎在路面(12)上的滚动噪音的方法的过程中:检测了路面(12)的实轮廓(16);基于所述实轮廓(16)确定路面(12)的至少一个描述符的至少一个值;使用了用于预测轮胎的滚动噪音B的规律,其中B=g(V1,...,Vn),其中V1,…,Vn为路面(12)的描述符D1,…,Dn的值,其中1≤n≤2,且D1为特征的路面的压头数,且D2为特征的路面的压头尺寸。

Description

预测轮胎运行噪音的方法
技术领域
本发明涉及轮胎的技术领域,更特别地涉及轮胎在道路路面上的滚动噪音的技术领域。 
背景技术
从现有技术中已知了预测轮胎在道路路面上的滚动噪音的方法。在该方法中,获得路面的三维实轮廓。所述实轮廓包括由掺入基体(例如柏油)的鹅卵石或砾石碎屑所形成的压头。 
为了简化所获得的实轮廓的使用,确定实压头的实描述符例如包括压头的平均高度、压头的平均密度等。 
接着,使用预测轮胎的滚动噪音的规律。所述规律取决于被采用为与滚动噪音相关的许多描述符。该规律由用于预测噪音的经验规律与所确定的描述符之间的关系所决定。因此所述规律具有基本上相关的性质。 
如今,所述规律所依赖的大量描述符不允许快速计算所预测的噪音。此外,某些描述符基于复杂的计算表,或基于求解相对漫长的数字方程体系而获得。 
最后,所采用的描述符不总是与描述路面和轮胎之间的物理相互作用的物理现实相关。因此这些描述符被不必要地进行计算,而不会改进所预测噪音的准确性。此外,未清楚建立噪音和路面的描述符之间的关系。 
所有这些缺点导致使用在测量和计算时间方面冗长,且在所预测的噪音的准确性方面平平的方法。 
发明内容
本发明的目的是确定一个或多个描述符的类型,以及与用于准确且简单预测轮胎的滚动噪音相关和必需的一个或多个描述符。 
为此目的,本发明的主题为一种预测轮胎在路面上的滚动噪音的方法,其中: 
-绘制路面的实轮廓, 
-基于所述实轮廓确定路面的至少一个描述符的至少一个值, 
-使用B=g(V1,...,Vn)类型的用于预测轮胎的滚动噪音的规律,其中V1,…,Vn为路面的描述符D1,…,Dn的值,且其中1≤n≤2,且 
-D1为路面压头的特征数目,且 
-D2为路面压头的特征尺寸。 
本发明人已出乎意料地发现,单独考虑路面的两个描述符D1和D2就与准确预测滚动噪音有关。因此,数目类型的描述符D1和尺寸类型的描述符D2使得以简单方式预测滚动噪音成为可能。此外,本发明人已发现路面压头的特征数目相比于特征尺寸具有更显著的影响,使得在n=1的情况中,单独的压头的特征数目有可能预测滚动噪音。 
这种预测规律对于路面设计特别有利。实际上,不同于在铺设路面的道路上例如通过声音传感器实验测量滚动噪音,有可能基于仅一米或数米的道路样品来预测噪音。这提供了道路建设成本的节约和实验成本的节约。 
在给定条件下滚动的给定轮胎的滚动噪音基本上取决于路面的粗糙度。该粗糙度导致轮胎的振动,所述轮胎的振动产生噪音。路面也显示吸收和反射特性。滚动噪音随这些特性而变化。在已考虑解决物理相互作用的效应的物理方法下,根据本发明的方法可以分开这些特性的每一个的贡献,且道路设计者有可能确定导致有或多或少噪音的路面的一个或多个特性。相反,噪音测量是综合性的,无法区分噪音起源的物理机理。单独的测量不提供特性和滚动噪音之间的清楚关系,且不允许控制路面设计。 
在第一实施方案中,噪音预测规律具有如下类型: 
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo
其中R为压头的特征尺寸值,N为压头的特征数目值,C为轮胎的滚动速度,E为轮胎橡胶的刚度模量,v为轮胎橡胶的泊松比,且Ro、No、Co、Eo和Bo为参照值。k属于区间[1.2-2]。 
因此有可能预测在任何路面上滚动的任何轮胎的滚动噪音。 
在第二实施方案中,噪音预测规律具有如下类型: 
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ( C Co ) k ] + Bo '
其中R为压头的特征尺寸值,N为压头的特征数目值,C为轮胎的滚动速度,且Ro、No、Co和Bo为参照值。k属于区间[1.2-2]。 
对于给定轮胎,有可能使用根据第二实施方案的规律来相互比较数个路面。实际上,该规律不包含与轮胎相关的任何项。 
在第三实施方案中,噪音预测规律具有如下类型: 
B = 20 log [ ( No N ) 3 5 ( C Co ) k ] + Bo ' '
其中N为压头的特征数目值,C为轮胎的滚动速度,且No、Co和Bo″为参照值。k属于区间[1.2-2]。 
在第四实施方案中,噪音预测规律具有如下类型: 
B = 20 log [ ( No N ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo ' ' '
其中N为压头的特征数目值,C为轮胎的滚动速度,E为轮胎橡胶的刚度模量,v为轮胎橡胶的泊松比,且No、Co、Eo和Bo″′为参照值。k属于区间[1.2-2]。 
在已考虑涉及特征尺寸的项的幂(其与涉及速度和涉及路面压头的特征数目的项的幂相关)下,有可能获得对滚动噪音的良好近似,而无需通过第三和第四实施方案的规律来确定路面压头的特征尺寸。 
在第五实施方案中,噪音预测规律具有如下类型: 
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ] + Bo ' ' ' '
其中R为压头的特征尺寸值,N为压头的特征数目值,且Bo″″为参照值。 
根据所述方法的一个任选特性: 
-基于实轮廓,路面的至少一个描述符的至少一个变化规律确定为随实轮廓的可变参数而变化, 
-对于每个变化规律,在参数的所谓的相互作用值下确定每个描述符的所谓的相互作用值,以及 
-基于至少一个描述符的相互作用值来预测滚动噪音值。 
通过确定变化规律并通过寻找每个描述符的相互作用值,因此有可能正确描述路面和轮胎之间的物理相互作用,并因此正确描述噪音。 
任选地,将实轮廓划分为数个层,每个层对应于实轮廓的一个高度: 
-每个描述符的变化规律确定为随实轮廓的高度而变化,以及 
-对于每个变化规律,在高度的相互作用值下确定每个虚描述符的相互作用值。 
可变参数为高度z。高度的相互作用值对应于轮廓高度,在该高度下的轮胎挤压如实表示轮胎在路面上的实挤压。 
任选地: 
-一个或多个描述符包括至少一个路面的实压头的所谓的实描述符; 
-确定每个实描述符的相互作用值; 
-基于至少一个实描述符的相互作用值来预测滚动噪音值。 
根据所述方法的其他任选特性: 
-实描述符为路面的每单位面积的实压头数目,压头的特征数目值为路面的每单位面积的实压头数目的相互作用值。 
-相互作用值为路面的每单位面积的实压头数目的变化规律的最大值。 
与现有技术噪音预测方法(其中必须计算数十个描述符)不同,容易地计算每单位面积路面的压头数目。 
任选地: 
-一个或多个描述符包括用于模拟实压头的虚压头的至少一个所谓的虚描述符; 
-在一方面,基于轮胎和路面之间的相互作用规律来确定每个虚描述符的相互作用值,在另一方面,基于虚压头的每个变化规律来确定每个虚描述符的相互作用值,以及 
-基于至少一个虚描述符的相互作用值来预测滚动噪音值。 
根据所述方法的另一任选特征,虚压头由每个可变参数值的相同球体表示,虚描述符为球体的半径,特征尺寸值为球体半径的相互作用值。 
对于相互作用值,每个实压头通过虚压头进行模拟,所述虚压头表达轮胎和路面之间的相互作用的物理现实而不是如现有技术中的几何现实。因此,待操作的数据量更小。实际上,对于每个参数值,尽管虚压头模拟的实压头各不相同,但虚压头均相同。 
本发明的主题也是一种计算机程序,一种计算机程序,其包括代码指令,当该计算机程序在计算机上执行时,所述代码指令能够控制例如以上限定的方法的步骤的执行。 
本发明的主题也是一种记录数据的媒介,其已记录的形式包括例如上文所定义的程序。 
本发明的主题也是使诸如以上限定的程序能够在远程通信网络上获得,其目的是为了它的下载。 
最后,本发明的主题为一种制造路面的方法,其包括根据例如上文所定义的方法的预测步骤。 
附图说明
通过阅读如下说明书并同时参照附图可更好地理解本发明,所述说明书仅以非限制性的实例给出: 
-图1表示道路路面的实轮廓和已过滤的路面的实轮廓; 
-图2显示了将实轮廓划分为多个层的步骤; 
-图3表示实轮廓的实压头数目的变化规律随高度z的变化; 
-图4显示了通过实压头确定轮胎的实变形的步骤; 
-图5显示了确定轮胎中的由实压头压入的实体积的步骤; 
-图6显示了通过根据本发明的方法获得的虚轮廓; 
-图7表示轮胎和路面之间的模拟的相互作用; 
-图8显示了用于确定轮胎的虚变形和由轮胎中的虚压头压入的虚体积所参照的球体的透镜部分; 
-图9和10表示虚描述符的变化规律随高度z的变化。 
在图中,已表示出对应于路面的常规纵向(X)、横向(Y)和竖 直(Z)取向的互相正交的轴X、Y、Z。 
具体实施方式
已在图1和2中显示了旨在与轮胎相互作用的道路路面12。路面12包括其中掺入不同尺寸的砾石碎屑的基体,所述基体包含柏油、沥青和/或水泥。这些砾石碎屑(也称为粗砂)在路面12的表面上形成压头14。压头14旨在与轮胎接触。 
路面12根据下文所述的预测方法而制得。 
优选利用激光轮廓仪绘制路面12的样品的二维的所谓实轮廓16。在此情况中,实轮廓16的深度在长度L的路面上进行测量。通过获得二维轮廓16,对于绘制三维轮廓,收集时间和待实施的硬件减少。实轮廓16在平面X,Z中在图1中由虚线表示,并显示出波峰18和波谷20。认为每个实轮廓16表示路面的厚度为dy的切片。 
作为变型,绘制数个二维轮廓16以构建三维样品。 
轮廓16通过截止滤波器22进行过滤,所述截止滤波器22由图1中的折线表示,并称为波峰线。波峰线22选择并连接显示高度z的波峰18,使得在轮廓16上滚动的轮胎变得仅与所选的波峰22接触。最后,通过将所选的波峰22的高度重新初始化至零高度而调整轮廓16。然后获得经调整的轮廓24,其由图1中的实线表示。 
在图2中已经显示了轮廓24。将轮廓24离散化为随轮廓24的可变参数而变化。在此情况中,可变参数为轮廓24的高度z。然后沿Z轴将轮廓24划分为数个层26,例如一百个层,每个层对应于一个高度z。所述层基本上为平面,并平行于X,Y平面。 
定义了压头14的所谓的实描述符NR、γR和VR。指数R表明这些描述符基于实轮廓24测得。也定义了高度z的剖面,其对应于垂直于Z轴的水平面并在高度z延伸。对于每个实描述符NR、γR和VR,现在将计算随可变参数z变化的变化规律。 
参照图2和3,NR表示每单位面积的实压头14的数目。对于每个高度z,计算由高度z的每个剖面所截取的每单位面积的轮廓24的压头14的数目。然后变化规律NR(z)确定为随可变参数z而变化。接着,基于变化规律NR(z)确定所谓的相互作用值NR,E。相互作用值NR,E为路 面12的每单位面积的实压头14的数目的最大值。在图3中,对于zR,E=17mm,每平方米为NR,E=115个压头。 
在三维轮廓的情况中,对于每个高度z确定与轮胎接触的实压头14的数目。对于路面的同一个切片dy而言,如果波峰18在所考虑的片dy中不共面,则与轮胎接触的实压头14的数目变化。 
在二维轮廓的情况中,假设对于每个高度z,实压头14的数目的最大值为实压头14的数目的实值。实际上,局部极大值通常表明在所研究的切片dy中沿横向方向存在波峰18。 
参照图4,γR表示由每个压头14所导致的轮胎的实局部变形。对于每个高度z,计算由γz,i=cotan(αi)定义的每个变形γz,i,其中αi为在高度z存在的每个压头14i的顶点的半角。接着,对于每个高度z,计算变形γz,i的平均值。因此变化规律γR(z)确定为随着可变参数z而变化。 
参照图5,VR表示在轮胎中由实压头14压入的每单位面积的实体积。对于每个高度z,计算由 
Figure BDA0000132790940000071
定义的在高度z的由每个压头14i压入的每个体积ΔVz,i,其中Lcz,i为在高度z的每个压头i的弦长,hm,i为在长度Lcz,i上每个压头14i的平均高度。接着通过加入体积ΔVz,i(z)来计算每个高度z的VR(z)。因此确定变化规律VR(z)随可变参数z而变化,并由 
Figure BDA0000132790940000072
所定义,其中L表示轮廓24的长度。 
路面12通过在图6中表示的虚轮廓28进行模拟。虚轮廓28包括用于模拟实压头14的虚压头30。定义虚压头30的所谓虚描述符。每个虚描述符根据变化规律而随着可变参数z进行变化。 
这些虚描述符包括主要虚描述符av、Rv、hv和次要虚描述符γV、VV。指数V表明所述描述符为虚的。次要虚描述符γV、VV和实描述符γR、VR描述了轮胎分别与虚压头30和实压头14之间的相互作用的类似物理量。以类似于实轮廓24的方式,定义虚轮廓28的高度z的剖面,所述高度z的剖面对应于垂直于Z轴的水平面并在高度z延伸。 
虚压头30由球体32表示。作为变型,虚压头30可由适合于待模拟的路面的类型和适合于物理相互作用的椎体、圆柱条或其他形状表示。 
轮胎和半径Rv的虚球体32之间的相互作用已在图7中表示。对 于施加于轮胎的负载Z,球体32产生深度为2hv的压痕,并包括与具有最大半径av和高度hv的轮胎接触的透镜状部分34。轮胎和球体32之间的相互作用遵循Hertz相互作用定律: 
a v = 3 ( 1 - v 2 ) . Z . R v 4 . E
其中v和E分别为轮胎橡胶的泊松比和刚度模量。 
半径为Rv的球体32的透镜状部分34在图8中表示。对于每个次级虚描述符γv、Vv,现在将计算每个变化规律随主要虚描述符av、Rv、hv的至少一个的变化,在本文随av和Rv的变化。 
γv表示由每个虚压头30(在本文由半球部分34)导致的轮胎的虚变形。实变形基于在两个维度上的实压头14的轮廓16进行计算。虚变形基于在三个维度上的球体32进行计算。因此,根据穿过每个球体的X,Z平面组(就像球体32沿虚轮廓28的随机设置),认为虚变形对应于由每个球体32所导致的变形的平均值。因此,基于二维实轮廓,所述方法有可能确定虚轮廓的三维描述符。 
在每个剖面X,Z中,平均变形对应于由中间接触点施加的变形,所述中间接触点即坐标为x=(av.cos θ)/2,y=(av.sin θ)/2的点。在中间接触点的斜率由如下关系给出: 
dz dx = - a v 2 . cos θ . [ R v 2 - a v 2 4 . ( 1 - 3 . sin 2 θ ] - 1 2
为了获得在剖面X,Z组上的轮胎的虚变形γV,我们计算: 
γ m , v = 2 π ∫ 0 π / 2 - a v 2 . cos θ . [ R v 2 - a v 2 4 . ( 1 - 3 . sin 2 θ ] - 1 2 dθ
因此,变化规律γv(z)确定为随av(z)和Rv(z)而变化,因此随z而变化: 
γ m , v ( z ) = a v π ( R v 2 - a v 2 4 ) 1 2 + 1 8 π [ a v ( R v 2 - a v 2 4 ) 1 2 ] 3
VV表示由轮胎中的虚压头30(在本文由半球部分34)压入的虚体 积。 
因此,变化规律VV(z)确定为以如下方式随av(z)和Rv(z)而变化,因此随z而变化: 
V V ( z ) = π . R v 3 . ( 2 3 + cos 3 β 3 - cos β ) 其中 cos β = R v - h v a v
接着,对于可变参数z的每个值,基于每个实描述符γR(z)、VR(z)的每个变化规律和每个次要虚描述符γv(z)、Vv(z)的每个变化规律,确定每个主要虚描述符av(z)、Rv(z)的值。在此情况中,对于可变参数z的每个值,在一方面在实压入体积VR(z)和虚压入体积Vv(z)之间建立等式,在另一方面在实变形γR(z)和虚变形γv(z)之间建立等式。因此每个主要虚描述符av(z)、Rv(z)的变化规律确定为随可变参数z而变化。然后可谈及具有施加的变形和施加的体积的模型。 
将会注意到,对于可变参数z的每个值,虚球体32是相同的,并均由相同的主要描述符av(z)、Rv(z)进行描述。 
具有施加的变形和施加的体积的主要虚描述符av(z)、Rv(z)的变化规律分别在图9A和9B中表示。 
然后,在一方面基于上文指出的轮胎和路面之间的相互作用的Hertz规律,在另一方面基于主要虚描述符av(z)、Rv(z)的每个变化规律,确定可变参数的所谓的相互作用值(本文为高度z的相互作用值zE)。在示出的实例中,数字求解出由这三个规律形成的系统,获得zE=0.19mm。对于每个变化规律,也确定在相互作用值zE下的每个描述符av(z)、Rv(z)的相互作用值。获得av(zE)=1.85mm和Rv(zE)=5.14mm。 
接着,预测在轮胎和路面之间的物理相互作用。在此情况中,基于用于预测B=g(V1,...,Vn)类型的噪音B的规律,预测轮胎在路面12上的滚动噪音B,其中V1,...,Vn为路面的描述符D1,...,Dn的值。在第一实施方案中,n=2,D1为路面12的压头的特征数目,且D2为路面12的压头的特征尺寸。在示出的实例中,压头的特征数目为每单位面积路面12的实压头12的数目的相互作用值NR,E,压头的特征尺寸为虚球体32的半径的相互作用值Rv(zE)。 
在该实施方案中,基于如下定义的滚动噪音B预测规律来预测噪音B: 
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo (等式1) 
其中C为轮胎的滚动速度,v为橡胶的泊松比,且E为橡胶的刚度模量,Bo为在具有参照虚描述符Ro和No的参照路面上以参照速度Co滚动的具有刚度模量Eo的参照轮胎的参照滚动噪音。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另一变型中,k=2。Bo取决于进行测量的位置:在带有轮胎的车辆内部,在铺设路面的轨道边缘等。 
在第二实施方案中,n、D1和D2与第一实施方案相同。基于如下定义的滚动噪音B预测规律来预测噪音B: 
B = 20 log [ ( R v ( z E ) Ro ) 1 5 ( No N R , E ) 3 5 ( C Co ) k ] + Bo ' (等式2) 
当使用参照轮胎比较数个不同的路面时,橡胶的刚度模量E和泊松比v是已知的,并包含于参照常数Bo′中。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另一变型中,k=2。 
在第三实施方案中,在已考虑对于某些类型的路面描述符Rv相对于其他描述符具有可忽略的影响下,基于如下定义的滚动噪音B预测规律来预测给定轮胎的滚动噪音B: 
B = 20 log [ ( No N R , E ) 3 5 ( C Co ) k ] + Bo ' ' (等式3) 
在该实施方案中,n=1,路面压头的特征数目D1为仅有的所用的描述符。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另一变型中,k=2。通常观察到在通过等式2和3所预测的噪音B之间的小差异。 
在第四实施方案中,基于如下定义的滚动噪音B预测规律来预测任意轮胎的滚动噪音B: 
B = 20 log [ ( No N R , E ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo ' ' ' (等式4) 
在该实施方案中,n=1,路面压头的特征数目D1为仅有的所用的描述符。k属于区间[1.2-2]。优选地,k=1.6。作为变型,k=1.2。在另 一变型中,k=2。 
在第五实施方案中,基于如下定义的滚动噪音B预测规律来预测给定轮胎的滚动噪音B: 
B = 20 log [ ( R v ( z E ) Ro ) 1 5 ( No N R , E ) 3 5 ] + Bo ' ' ' ' (等式5) 
仅考虑路面的特性的该关系可以获得良好的预测。实际上,所述规律考虑Rv(zE),尽管其显示比NR,E更低的数学权数,但其可以间接描述用于制造路面的方法的特性,如岩石性质、破碎模式的影响。该关系也可以了解用于制造路面的方法的特性的更改如何更改了滚动噪音。 
预测给定雕刻式轮胎的滚动噪音B的实例。
利用在各种路面T1-T4上的雕刻式轮胎进行试验。已根据本发明的方法确定每个这些路面的描述符NR,E和Rv(zE)。利用激光轮廓仪获得每个路面T1-T4的1米样品的实轮廓,所述激光轮廓仪的空间分辨率沿X轴为dx=0.1mm,沿Z轴为dz=0.01mm。使用第二实施方案的预测规律(等式2),其中,值Bo′=45dB。该Bo′值对应于对客运车辆轮胎在车辆外部通常测得的试验平均值。 
给定车辆在每个路面T1-T4上以60km/h行驶。然后在三个不同位置测量噪音:接近轮胎、在车辆内部和当车辆经过时在车辆外部(也称为coast-by噪音)。结果在下文表1中给出。 
Figure BDA0000132790940000112
表1 
对于接近轮胎测得的噪音,获得如下关系B测量=0.73×B预测+52,相关系数R2=0.8。 
对于在车辆外部轮胎测得的噪音,获得如下关系B测量=1.26×B预测-21.40,相关系数R2=0.96。 
对于在车辆内部轮胎测得的噪音,获得如下关系B测量=1.23×B预测-22.10,相关系数R2=0.97。 
注意到,可通过根据测量位置校正Bo′的值而更准确地估计实验值B测量。 
无论测量位置是哪里,预测模型和进行的测量之间的关系是非常令人满意的,并有可能基于一个或多个路面的样品比较并预测由给定轮胎的滚动所产生的噪音B。 
预测给定无花纹轮胎的滚动噪音B的实例。
为了表征道路的粗糙度而不与轮胎的刻纹相互作用,使用无花纹轮胎。在显示道路网的路面的多样性的各种路面R1-R6上进行试验。已根据本发明的方法确定每个这些路面的描述符NR,E和Rv(zE)。利用激光轮廓仪获得给定路面的1.7m样品的实轮廓,所述激光轮廓仪的空间分辨率沿X轴为dx=0.4mm,沿Z轴为dz=0.01mm。使用第二和第三实施方案的预测规律(等式2和3),值Bo′=Bo″=45dB。 
给定车辆在每个路面R1-R6上以60km/h行驶。进行接近轮胎的数个噪音实验测量,并平均每个路面R1-R6的实验测量值。结果在下文表2中给出。 
Figure BDA0000132790940000121
表2 
在一方面,注意到路面组R1-R6的滚动噪音B的预测是优良的,由此确定所选择描述符的相关性以及用于非常不同的路面的滚动噪音 B预测规律的适合性。 
此外,注意到,利用根据第三实施方案的规律所预测的噪音B的值大于利用根据第二实施方案的规律所预测的噪音B的值。在该实施例的架构内,相比于根据第二实施方案的规律,根据第三实施方案的规律可以获得更接近于实验值的噪音B的预测值。 
最后,注意到预测的噪音B的精确度大于之前实施例的预测的噪音B的精确度。 
根据本发明的方法的全部或部分可通过代码指令实施,当所述方法的步骤在计算机上执行时,所述代码指令能够控制所述方法的步骤的执行。所述指令可以从用于记录数据(所述数据例如具有硬盘或闪存、CD或DVD类型)的媒介上所记录的计算机程序发出。可以为这种程序做好准备使其可以利用,其目的是为了在远程通信网络(如英特网或无线网)上进行下载。因此,可以经由该网络将程序的升级发送至与所述网络连接的计算机。 

Claims (17)

1.一种预测轮胎在路面(12)上的滚动噪音(B)的方法,其特征在于:
-绘制路面(12)的实轮廓(16),
-基于所述实轮廓(16)确定路面(12)的至少一个描述符(av、Rv、NR)的至少一个值(av(zE)、Rv(zE)、NR,E),
-使用B=g(V1,...,Vn)类型的用于预测轮胎的滚动噪音B的规律,其中V1,…,Vn为路面(12)的描述符D1,…,Dn的值,1≤n≤2,且
-D1为路面压头的特征数目(NR),且
-D2为路面压头的特征尺寸(Rv)。
2.根据权利要求1所述的方法,其中所述噪音预测规律具有如下类型:
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo
其中R为压头的特征尺寸值(Rv(zE)),N为压头的特征数目值(NR,E),C为轮胎的滚动速度,E为轮胎橡胶的刚度模量,v为轮胎橡胶的泊松比,k属于区间[1.2-2],且Ro、No、Co、Eo和Bo为参照值。
3.根据权利要求1所述的方法,其中所述噪音预测规律具有如下类型:
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ( C Co ) k ] + Bo '
其中R为压头的特征尺寸值(Rv(zE)),N为压头的特征数目值(NR,E),C为轮胎的滚动速度,k属于区间[1.2-2],且Ro、No、Co和Bo为参照值。
4.根据权利要求1所述的方法,其中所述噪音预测规律具有如下类型:
B = 20 log [ ( No N ) 3 5 ( C Co ) k ] + Bo ' '
其中N为压头的特征数目值(NR,E),C为轮胎的滚动速度,k属于区间[1.2-2],且No、Co和Bo″为参照值。
5.根据权利要求1所述的方法,其中所述噪音预测规律具有如下类型:
B = 20 log [ ( No N ) 3 5 ( C Co ) k ( 0.75 1 - v 2 E Eo ) 1 5 ] + Bo ' ' '
其中N为压头的特征数目值(NR,E),C为轮胎的滚动速度,E为轮胎橡胶的刚度模量,v为轮胎橡胶的泊松比,k属于区间[1.2-2],且No、Co、Eo和Bo″′为参照值。
6.根据权利要求1所述的方法,其中所述噪音预测规律具有如下类型:
B = 20 log [ ( R Ro ) 1 5 ( No N ) 3 5 ] + Bo ' ' ' '
其中R为压头的特征尺寸值,N为压头的特征数目值,且Bo″″为参照值。
7.根据前述权利要求任一项所述的方法,其中:
-基于实轮廓(16),路面(12)的至少一个描述符(av、Rv、NR)的至少一个变化规律(av(z)、Rv(z)、NR(z))被确定为随实轮廓(16)的可变参数(z)而变化,
-对于每个变化规律,在参数的所谓的相互作用值(zR,E、zE)下确定每个描述符(av、Rv、NR)的所谓的相互作用值(av(zE)、Rv(zE)、NR,E),以及
-基于至少一个描述符(Rv、NR)的相互作用值(Rv(zE)、NR,E)来预测滚动噪音值(B)。
8.根据权利要求7所述的方法,其中将实轮廓(16)划分为数个层,每个层对应于实轮廓(16)的高度(z):
-每个描述符(av、Rv、NR)的变化规律(av(z)、Rv(z)、NR(z))确定为随实轮廓(16)的高度(z)变化,以及
-对于每个变化规律(av(z)、Rv(z)、NR(z)),在高度(z)的相互作用值(zE)下确定每个虚描述符的相互作用值(av(zE)、Rv(zE)、NR,E)。
9.根据前述权利要求任一项所述的方法,其中:
-一个或多个描述符(av、Rv、NR)包括至少一个路面(12)的实压头(24)的所谓的实描述符(NR);
-确定每个实描述符(NR)的相互作用值(NR,E);
-基于至少一个实描述符(NR)的相互作用值(NR,E)来预测滚动噪音值(B)。
10.根据权利要求9所述的方法,其中所述实描述符为路面的每单位面积的实压头数目(NR),压头的特征数目值为路面的的每单位面积实压头数目(NR)的相互作用值(NR,E)。
11.根据权利要求10所述的方法,其中相互作用值(NR,E)为路面的每单位面积的实压头数目(NR)的变化规律(NR(z))的最大值。
12.根据权利要求9至11任一项所述的方法,其中:
-一个或多个所述描述符(av、Rv、NR)包括用于模拟实压头(14)的虚压头(24)的至少一个所谓的虚描述符(av、Rv);
-在一方面,基于轮胎和路面(12)之间的相互作用规律来确定每个虚描述符(av、Rv)的相互作用值(av(zE)、Rv(zE)),在另一方面,基于虚压头(24)的每个变化规律(av(z)、Rv(z))来确定每个虚描述符(av、Rv)的相互作用值(av(zE)、Rv(zE)),以及
-基于至少一个虚描述符(Rv)的相互作用值(Rv(zE))来预测滚动噪音值(B)。
13.根据权利要求12所述的方法,其中对于每个可变参数(z)值而言,所述虚压头(30)由相同的球体(32)表示,所述虚描述符为球体(32)的半径(Rv),所述特征尺寸值为球体(32)半径的相互作用值Rv(zE)。
14.一种计算机程序,其特征在于它包括代码指令,当该计算机程序在计算机上执行时,所述代码指令能够控制根据前述权利要求任一项所述的方法的步骤的执行。
15.一种用于记录数据的媒介,其以记录的形式包括根据权利要求14所述的程序。
16.使根据权利要求14所述的程序在远程通信网络上可用,其目的是它的的下载。
17.一种制造路面的方法,其特征在于它包括根据权利要求1至13任一项所述的方法的预测步骤。
CN201080033140.9A 2009-07-28 2010-07-19 预测轮胎运行噪音的方法 Active CN102472691B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR0955276A FR2948764B1 (fr) 2009-07-28 2009-07-28 Procede de prevision d'un bruit de roulement d'un pneumatique
FR0955276 2009-07-28
PCT/FR2010/051515 WO2011012793A1 (fr) 2009-07-28 2010-07-19 Procède de prévision d'un bruit de roulement d'un pneumatique

Publications (2)

Publication Number Publication Date
CN102472691A true CN102472691A (zh) 2012-05-23
CN102472691B CN102472691B (zh) 2016-01-20

Family

ID=41786204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201080033140.9A Active CN102472691B (zh) 2009-07-28 2010-07-19 预测轮胎运行噪音的方法

Country Status (6)

Country Link
US (1) US9097576B2 (zh)
EP (1) EP2459982B1 (zh)
JP (1) JP5654593B2 (zh)
CN (1) CN102472691B (zh)
FR (1) FR2948764B1 (zh)
WO (1) WO2011012793A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108614935A (zh) * 2018-04-24 2018-10-02 哈尔滨工大泰铭科技有限公司 一种快速实现轮胎花纹节距噪声预测的方法
CN115782480A (zh) * 2022-11-25 2023-03-14 安徽佳通乘用子午线轮胎有限公司 一种滚动轮胎结构噪声时域预测方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2948765B1 (fr) * 2009-07-28 2013-10-18 Michelin Soc Tech Procede de prevision d'un effet physique d'interaction entre un pneumatique et un revetement routier
FR2948764B1 (fr) * 2009-07-28 2011-08-26 Michelin Soc Tech Procede de prevision d'un bruit de roulement d'un pneumatique
US8793004B2 (en) * 2011-06-15 2014-07-29 Caterpillar Inc. Virtual sensor system and method for generating output parameters
FR3090103A3 (fr) 2018-12-14 2020-06-19 Michelin & Cie Modèle de prédiction de l’usure et de la fin de vie d’un pneumatique
WO2020120923A1 (fr) 2018-12-14 2020-06-18 Compagnie Generale Des Etablissements Michelin Modèle de prédiction de l'usure et de la fin de vie d'un pneumatique

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002090264A (ja) * 2000-09-21 2002-03-27 Yokohama Rubber Co Ltd:The タイヤ車外騒音予測方法およびタイヤ車外騒音予測プログラムを記録した記録媒体
CN1486419A (zh) * 2001-01-12 2004-03-31 �����ּ�����˾ 用于预测车内舒适度的方法
JP2007230458A (ja) * 2006-03-02 2007-09-13 Toyo Tire & Rubber Co Ltd タイヤからの放射騒音シミュレーション方法

Family Cites Families (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3864967A (en) * 1973-08-20 1975-02-11 Electronic Associates Tire load force variation testing system
US4727501A (en) * 1987-02-27 1988-02-23 The Uniroyal Goodrich Tire Company Method of simulating tire tread noise and apparatus
US5065618A (en) * 1990-11-14 1991-11-19 Hodges Transportation Inc. Method and apparatus for determining terrain surface profiles
US5880362A (en) 1995-09-06 1999-03-09 Engineering Technology Associates, Inc. Method and system for simulating vehicle and roadway interaction
US5610330A (en) * 1996-01-16 1997-03-11 Ford Motor Company Effective road profile control method for a spindle-coupled road simulator
JP3595411B2 (ja) * 1996-04-22 2004-12-02 株式会社ブリヂストン 騒音シミュレート方法及びシミュレータ
JP3318235B2 (ja) * 1997-06-16 2002-08-26 株式会社ブリヂストン 車内騒音の予測方法及び装置
DE69804723T2 (de) * 1997-07-18 2002-11-21 Kabushiki Kaisha Toyota Chuo Kenkyusho, Aichi Radzustands- Abschätzungsvorrichtung
JP3331310B2 (ja) * 1997-09-25 2002-10-07 富士重工業株式会社 路面摩擦係数検出装置
US6077302A (en) * 1998-02-12 2000-06-20 Egs, Inc. System and method for analyzing and designing vibration isolators
EP1489532A3 (en) * 1998-09-07 2006-03-08 Bridgestone Corporation Method of estimating tire performance
JP4137271B2 (ja) 1999-02-23 2008-08-20 株式会社ブリヂストン タイヤの振動・騒音シミュレーション方法
SE514309C2 (sv) * 1999-05-28 2001-02-05 Ericsson Telefon Ab L M Förfarande för att bestämma ytstruktur
JP2001165647A (ja) * 1999-12-06 2001-06-22 Toyota Central Res & Dev Lab Inc 路面騒音評価装置
JP3650342B2 (ja) * 2001-05-28 2005-05-18 住友ゴム工業株式会社 タイヤ・ホイール性能のシミュレーション方法及び装置
JP3927080B2 (ja) * 2002-06-12 2007-06-06 住友ゴム工業株式会社 タイヤのシミュレーション方法
JP4087193B2 (ja) * 2002-08-23 2008-05-21 株式会社ブリヂストン 騒音予測装置及び方法
JP4177671B2 (ja) * 2003-01-10 2008-11-05 住友ゴム工業株式会社 タイヤ情報の提供方法及び提供装置
JP3940093B2 (ja) * 2003-04-28 2007-07-04 住友ゴム工業株式会社 タイヤのシミュレーション方法
JP4529562B2 (ja) * 2004-07-06 2010-08-25 横浜ゴム株式会社 接触特性の評価方法及び接触状態の評価用コンピュータプログラム
JP4608306B2 (ja) * 2004-12-21 2011-01-12 住友ゴム工業株式会社 タイヤのシミュレーション方法
US7533010B2 (en) * 2005-09-01 2009-05-12 The Goodyear Tire & Rubber Company Method and system for tread pattern noise optimization
US8180605B1 (en) * 2005-11-23 2012-05-15 Livermore Software Technology Corporation Methods and systems for creating a smooth contact-impact interface in finite element analysis
EP1798648B1 (en) * 2005-12-13 2017-03-15 Sumitomo Rubber Industries, Ltd. Method for pneumatic tire simulation
JP4745845B2 (ja) * 2006-02-02 2011-08-10 東洋ゴム工業株式会社 タイヤからの放射騒音シミュレーション方法
JP4913418B2 (ja) * 2006-02-09 2012-04-11 東洋ゴム工業株式会社 タイヤからの放射騒音シミュレーション方法
DE102006041867B4 (de) * 2006-09-06 2008-12-04 Continental Automotive Gmbh Verfahren und Vorrichtung zur Überwachung des Rauschens eines Sensors
US7950273B2 (en) * 2007-11-09 2011-05-31 GM Global Technology Operations LLC Engine misfire and rough road detection systems and methods
JP4792049B2 (ja) * 2008-01-09 2011-10-12 住友ゴム工業株式会社 タイヤのノイズ性能のシミュレーション方法及びタイヤの製造方法
JP4621271B2 (ja) * 2008-07-29 2011-01-26 住友ゴム工業株式会社 タイヤのシミュレーション方法
JP5243145B2 (ja) * 2008-08-08 2013-07-24 株式会社ブリヂストン タイヤ放射音解析方法及びタイヤ放射音解析プログラム
JP2010230641A (ja) * 2009-03-06 2010-10-14 Toyota Motor Corp タイヤの有限要素モデルを用いたタイヤ作用力解析方法とそれを用いたタイヤ振動またはタイヤ騒音の解析方法
FR2948765B1 (fr) * 2009-07-28 2013-10-18 Michelin Soc Tech Procede de prevision d'un effet physique d'interaction entre un pneumatique et un revetement routier
FR2948764B1 (fr) * 2009-07-28 2011-08-26 Michelin Soc Tech Procede de prevision d'un bruit de roulement d'un pneumatique
US8555698B2 (en) * 2011-01-26 2013-10-15 Bridgestone Americas Tire Operations, Llc Engineered surfaces for laboratory tread wear testing of tires
US20120273306A1 (en) * 2011-04-28 2012-11-01 John Alan Pangrazio Activity Reporting System
WO2013090830A1 (en) * 2011-12-16 2013-06-20 University Of Southern California Autonomous pavement condition assessment
EP2801835B1 (en) * 2011-12-26 2019-07-24 The University of Tokyo Measurement method and measurement device
JP5662971B2 (ja) * 2012-07-11 2015-02-04 住友ゴム工業株式会社 タイヤのシミュレーション方法
KR20140121226A (ko) * 2013-04-05 2014-10-15 현대자동차주식회사 노면 프로파일과 로드 노이즈의 상관성 정량화 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002090264A (ja) * 2000-09-21 2002-03-27 Yokohama Rubber Co Ltd:The タイヤ車外騒音予測方法およびタイヤ車外騒音予測プログラムを記録した記録媒体
CN1486419A (zh) * 2001-01-12 2004-03-31 �����ּ�����˾ 用于预测车内舒适度的方法
JP2007230458A (ja) * 2006-03-02 2007-09-13 Toyo Tire & Rubber Co Ltd タイヤからの放射騒音シミュレーション方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GERT HEINRICH等: "Rubber friction, tread deformation and tire traction", 《WEAR》 *
J. CESBRON等: "Numerical and experimental study of multi-contact on an elastic half-space", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108614935A (zh) * 2018-04-24 2018-10-02 哈尔滨工大泰铭科技有限公司 一种快速实现轮胎花纹节距噪声预测的方法
CN115782480A (zh) * 2022-11-25 2023-03-14 安徽佳通乘用子午线轮胎有限公司 一种滚动轮胎结构噪声时域预测方法
CN115782480B (zh) * 2022-11-25 2024-08-16 安徽佳通乘用子午线轮胎有限公司 一种滚动轮胎结构噪声时域预测方法

Also Published As

Publication number Publication date
CN102472691B (zh) 2016-01-20
JP2013500491A (ja) 2013-01-07
EP2459982A1 (fr) 2012-06-06
FR2948764A1 (fr) 2011-02-04
JP5654593B2 (ja) 2015-01-14
FR2948764B1 (fr) 2011-08-26
EP2459982B1 (fr) 2014-05-07
US20120158383A1 (en) 2012-06-21
WO2011012793A1 (fr) 2011-02-03
US9097576B2 (en) 2015-08-04

Similar Documents

Publication Publication Date Title
CN102472691A (zh) 预测轮胎运行噪音的方法
Elhattab et al. Drive-by bridge damage monitoring using Bridge Displacement Profile Difference
Haigermoser et al. Road and track irregularities: measurement, assessment and simulation
Chen et al. Damage detection of highway bridges based on long-gauge strain response under stochastic traffic flow
Zhan et al. Bridge identification and damage detection using contact point response difference of moving vehicle
Cantero et al. Railway infrastructure damage detection using wavelet transformed acceleration response of traversing vehicle
He et al. Novel virtual simply supported beam method for detecting the speed and axles of moving vehicles on bridges
He et al. Mass‐normalized mode shape identification method for bridge structures using parking vehicle‐induced frequency change
CN112964551B (zh) 基于应变和频率的路面承载力动态测试方法
Tan et al. Extracting mode shapes from drive-by measurements to detect global and local damage in bridges
Mei et al. Bridge mode shape identification using moving vehicles at traffic speeds through non‐parametric sparse matrix completion
Ye et al. Investigating the pavement vibration response for roadway service condition evaluation
Keenahan et al. Determination of road profile using multiple passing vehicle measurements
CN103132411A (zh) 一种高速铁路轨道静态平顺性分析调整方法
CN107228772A (zh) 一种盾构隧道结构损伤评估方法
CN104932392B (zh) 碾压监控方法、装置及系统
CN116484510B (zh) 动力学行为分析方法、装置、计算机设备和存储介质
Zhang et al. Reference-free damage identification method for highway continuous girder bridges based on long-gauge fibre Bragg grating strain sensors
CN102472690B (zh) 用于预测轮胎和道路表面之间的物理相互作用效应的方法
Nouri et al. Probabilistic estimation of dynamic impact factor for masonry arch bridges using health monitoring data and new finite element method
Fu et al. Importance of examining particle movements in modelling ballast bed via discrete element method
Li et al. New approach for estimating international roughness index based on the inverse pseudo excitation method
CN105319155A (zh) 一种考虑粗糙表面特征的岩石毛管力确定方法
Gong et al. Analysis of tyre-pavement contact behaviour of heavy truck in full-scale test
Wu et al. Fast calibration for vibration-based pavement roughness measurement based on model updating of vehicle dynamics

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: COMPAGNIE GENERAL DES ETABLISSEMENTS MICHELIN

Free format text: FORMER OWNER: SOCIETE DE TECHNOLOGIE MICHELIN

Effective date: 20121026

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20121026

Address after: French Clermont Ferrand

Applicant after: Compagnie General Des Etablissements Michelin

Applicant after: Michelin Research & Technology Co., Ltd.

Address before: French Clermont Ferrand

Applicant before: Michelin Co., Ltd.

Applicant before: Michelin Research & Technology Co., Ltd.

C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20170928

Address after: French Clermont Ferrand

Patentee after: Compagnie General Des Etablissements Michelin

Address before: French Clermont Ferrand

Co-patentee before: Michelin Research & Technology Co., Ltd.

Patentee before: Compagnie General Des Etablissements Michelin