CN114218827B - 一种抑制轮胎空腔共振噪音的参数设计方法 - Google Patents
一种抑制轮胎空腔共振噪音的参数设计方法 Download PDFInfo
- Publication number
- CN114218827B CN114218827B CN202111495374.1A CN202111495374A CN114218827B CN 114218827 B CN114218827 B CN 114218827B CN 202111495374 A CN202111495374 A CN 202111495374A CN 114218827 B CN114218827 B CN 114218827B
- Authority
- CN
- China
- Prior art keywords
- boundary
- tire
- representing
- matrix
- tire cavity
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000013461 design Methods 0.000 title claims abstract description 13
- 230000002401 inhibitory effect Effects 0.000 title claims abstract description 13
- 239000011358 absorbing material Substances 0.000 claims abstract description 80
- 239000012530 fluid Substances 0.000 claims abstract description 33
- 239000000463 material Substances 0.000 claims abstract description 15
- 238000010521 absorption reaction Methods 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 81
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000004613 tight binding model Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 239000011148 porous material Substances 0.000 abstract 1
- 238000002474 experimental method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010013082 Discomfort Diseases 0.000 description 1
- 206010043268 Tension Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000011796 hollow space material Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000004814 polyurethane Substances 0.000 description 1
- 229920002635 polyurethane Polymers 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60C—VEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
- B60C99/00—Subject matter not provided for in other groups of this subclass
- B60C99/006—Computer aided tyre design or simulation
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明公开了一种抑制轮胎空腔共振噪音的参数设计方法,是将抑制轮胎空腔共振噪音在轮胎内表面贴附的吸声材料等效成满足Delany‑Bazley‑Miki理论模型的流体;再构建基于多域边界元法构建轮胎空腔声模态分析的声学模型;然后使用围道积分法求解非线性特征值问题的特征值,计算内部贴附吸声材料的轮胎空腔特征频率;最后基于轮胎空腔特征频率虚部来表征抑制轮胎空腔共振噪音效果,从而设计出抑制轮胎空腔共振噪音的多孔材料的最佳流阻率、厚度和宽度。本发明能计算内部贴附吸声材料的轮胎空腔频率,为选取抑制轮胎空腔共振噪音的吸声材料参数提供理论依据。
Description
技术领域
本发明是一种抑制轮胎空腔共振噪音的参数设计方法,为轮胎内部贴附吸声材料的参数设计提供依据。
背景技术
在汽车的行驶过程中,凹凸不平的路面对轮胎表面产生激励,轮胎表面发生振动,进而使得胎内空气产生振动。当空气振动达到共振频率时将导致轮胎空腔共振,因而产生较大的空腔共振噪音,这就是轮胎空腔噪声的产生机理;轮胎空腔共振频率与轮胎空腔的尺寸密切相关,一般在150Hz-250Hz的范围内,这个频率范围,噪声穿透力强,递减慢,直达人耳骨,使人的交感神经紧张等不适,对汽车的驾乘舒适性具有重要影响。
为了抑制轮胎在行驶过程种的空腔噪音,在轮胎内部贴附吸声材料是目前减小轮胎空腔共振的一种有效手段,其中是利用吸音材料的吸音效果抑制轮胎空腔噪音;但吸声材料会影影响轮胎的动态性能。当吸声材料尺寸太大会影响车轮重量以及材料的安装,增加油耗和成本;流阻选择不佳,影响抑制轮胎空腔噪声效果。现有技术中,如何选取吸声材料的最佳流阻以及尺寸目前尚无有效方法。
发明内容
本发明是为避免上述现有技术所存在的不足,提供了一种抑制轮胎空腔共振噪音的参数设计方法,以期能给出吸声材料流阻和尺寸对轮胎空腔共振频率的影响规律,从而为吸声材料的最佳流阻以及合适尺寸提供依据。
本发明为解决技术问题采用如下技术方案:
本发明一种抑制轮胎空腔共振噪音的参数设计方法,所述轮胎空腔是由轮胎内表面与轮辋共同形成的呈环空腔,在轮胎内表面上安装吸声材料层,并沿轮胎周向整圈铺设有环形吸声材料层,且所述环形吸声材料层的横断面呈矩形;其特点是:所述参数设计方法是按如下步骤进行:
步骤1、用Delany-Bazley-Miki等效流体模型表征轮胎内表面上贴附的环形吸声材料层;
步骤2、构建所述吸声材料层的轮胎三维模型,并对所述轮胎三维模型的表面划分单元网格;
步骤3、基于多域边界元法构建轮胎空腔声模态分析的声学模型,并基于所述声学模型,利用式(1)构建多域边界元特征方程:
式(1)中:A为多域边界元系数矩阵;λ表示内部贴附吸声材料轮胎的特征值;表示特征值λ对应的特征向量;
步骤4、使用围道积分法计算内部贴附吸声材料轮胎的特征频率λ及其对应的特征向量基于特征频率λ虚部来表征轮胎空腔共振噪音的抑制程度,从而以特征频率λ虚部绝对值最大对应的环形吸声材料层的流阻率σ,厚度C与宽度L作为能抑制轮胎空腔共振噪音的环形吸声材料层的最佳流阻率、厚度和宽度。
本发明所述的抑制轮胎空腔共振噪音的参数设计方法的特点也在于,所述步骤1是如下过程进行:
将所述环形吸声材料层等效为满足所述Delany-Bazley-Miki等效流体模型的流体,并利用式(2)和式(3)得到等效流体域Ω2中的复波数k2和复特征阻抗Z2:
k2=k1(1+0.122X-0.618+i0.180X-0.618) (2)
Z2=Z1(1+0.079X-0.632+i0.120X-0.632) (3)
式(2)和式(3)中:k1为在空气域Ω1中的波数,X表示中间量,且X=ρck1/2πσ,ρ为空气密度,c为空气中声速,σ为所述流体的流阻率,i为虚数单位,Z1为空气特征流阻,Z1=ρc。
步骤3中所述式(1)中的多域边界元系数矩阵A由式(4)表征:
式(4)中:表示在边界ΓA与边界ΓA所划分的网格相关的系数矩阵,其中,表示复数域,NA为边界ΓA进行网格划分后的三角形单元的个数;/>和/>表示边界ΓA与边界ΓB所划分的网格相关的系数矩阵,NB为边界ΓB进行网格划分后的三角形单元的个数;/>表示边界ΓB与边界ΓA所划分的网格相关的系数矩阵;/>和/>表示边界ΓB与边界ΓB所划分的网格相关的系数矩阵;/>和表示边界ΓD与边界ΓC所划分的网格相关的系数矩阵,ND为边界ΓD进行网格划分后的三角形单元的个数,NC为边界ΓC进行网格划分后的三角形单元的个数;
表示在边界ΓD与边界ΓD所划分的网格相关的系数矩阵;/>和/>表示边界ΓC与边界ΓC所划分的网格相关的系数矩阵;/>表示边界ΓC与边界ΓD所划分的网格相关的系数矩阵;ΓA表示轮胎的轮辋边界,ΓB和ΓC表示轮胎空腔中空气域Ω1与等效流体域Ω2的边界,ΓB表示空气域Ω1侧的边界,ΓC表示等效流体域Ω2侧的边界,边界ΓB和ΓC是共用边界,且NB=NC,但法向相反,ΓD表示轮胎空腔的内表面;矩阵Nd是求解多域边界元问题的总自由度,且Nd=NA+NB+NC+ND;并有:
式(5)-式(16)中:δlm为克罗内克δ函数,l≠m时,δlm=0,l=m时,δlm=1;G1(k1,xl,y)表示在空气域Ω1中的边界点xl与y之间的格林函数,G2(k2,xl,y)表示在等效流体域Ω2中的边界点xl与y之间的格林函数;ΓAA表示边界ΓA与边界ΓA,ΓAB表示边界ΓA与边界ΓB,ΓBA表示边界ΓA与边界ΓB,ΓBB表示边界ΓB与边界ΓB,ΓDC表示边界ΓD与边界ΓC,ΓDD表示边界ΓD与边界ΓD,ΓCC表示边界ΓC与边界ΓC,ΓCD表示边界ΓC与边界ΓD;n1(y)表示在空气域Ω1中的边界y处的外法向,n1(y)表示在等效流体域Ω2中的边界y处的外法向。
所述步骤4的围道积分法是按如下步骤进行:
步骤a、利用式(17)和式(18)建立两个Hankel矩阵,即矩阵和矩阵/>
式(17)和式(18)中,K是正整数,第l个矩阵为/>其表达式为:
式(19)中:C为复平面内的闭合积分约当曲线,T(z)为与多域边界元系数矩阵相关的中间矩阵,并且T(z)=PA(z)-1Q,矩阵和矩阵/>是两个随机非零矩阵,L是正整数且大于在复平面内所围成的区域C内包含的特征值的最大重数,且K×L大于闭合积分约当曲线C构成的区域内的总特征值数;
步骤b、利用式(20)对Hankel矩阵进行特征值分解:
式(20)中,矩阵U,是两个酉矩阵,矩阵VH是矩阵V的共轭转置,Σ是非负实对角矩阵且对角元素为σ1,σ2…σK×L;
步骤c、定义δ为正阈值,依次判断对角元素σ1,σ2…σK×L是否小于δ×σ1,若小于,则将相应对角元素置“0”;否则,保留相应对角元素;
步骤d、令Hankel矩阵求解Hankel矩阵/>的特征值γj及其特征向量ηj,并将特征值γj并转换为特征频率λj,λj=γjc/2π;
步骤e、利用式(9)计算内部贴附吸声材料的特征频率λj对应的真实的特征向量
式(21)中,矩阵S=[S0,S1,…,Sl,…,SK-1],Sl表示第l个矩阵,并有:
与已有技术相比,本发明有益效果体现在:
1、本发明方法分析了不同环形吸声材料层的流阻率σ,厚度C与宽度L对轮胎空腔共振频率虚部的结果,从而给出了环形吸声材料层的流阻和尺寸对轮胎空腔共振频率的影响规律,为环形吸声材料层的最佳流阻以及合适尺寸提供依据。
2、本发明方法将吸声材料等效为满足Delany-Bazley-Miki理论模型的流体的方式,从而使构复杂的环形吸声材料层简化成可用参数描述的流体模型,并可以使用本发明的方法计算不同环形吸声材料层的流阻率σ,厚度C与宽度L对轮胎空腔共振共振频率。
3、本发明方法使用多域边界元法建立内表面贴附的吸声材料的轮胎声学模型,边界条件可以实现轮胎结构与空腔内部声场的耦合,可以准确表征实际的内表面贴附的吸声材料的轮胎,使得本发明方法更加适用于轮胎实际模型。
附图说明
图1为本发明中轮胎结构示意图;
图2为本发明中的含吸声材料的轮胎横截面的尺寸示意图;
图3为本发明所建立的轮胎内部粘贴吸声材料的声学模型示意图;
图4为本发明方法中轮胎特征频率随吸声材料流阻率σ变化曲线;
图5为本发明方法中轮胎特征频率实部随吸声材料厚度C变化曲线;
图6为本发明方法中轮胎特征频率虚部随吸声材料厚度C变化曲线;
图7为本发明方法中轮胎特征频率实部随吸声材料宽度L变化曲线;
图8为本发明方法中轮胎特征频率虚部随吸声材料宽度L变化曲线;
图9为轮胎内部有无贴附吸声材料传递函数对比图。
图中标号:1吸声材料层,2a轮辋,2b轮辐,3充气轮胎,4胎面内表面,5轮胎空腔。
具体实施方式
本实施例中,参见图1和图2,本实施例中,安装在胎面内表面4上的吸声材料层1,吸声材料层1沿轮胎周向整圈铺设,横断面呈矩形,矩形断面的吸声材料层1的流阻率σ,厚度为C,宽度为L。轮胎空腔5是由轮胎内表面4与轮辋2a共同形成的呈环空腔,在轮胎内表面4上安装吸声材料层,并沿轮胎周向整圈铺设有环形吸声材料层1,且环形吸声材料层的横断面呈矩形;本实施例中,充气轮胎是指由轮胎的胎面内表面4与轮辋2a共同形成有轮胎空腔5的充气轮胎3;充气轮胎3装设在轮辋2a上,由轮辐2b支撑;通过轮辋2a将充气轮胎3的空腔进行密闭的方式,将充气轮胎3安装在车轮上。
如图1和图2所示,一种抑制轮胎空腔共振噪音的参数设计方法是按如下步骤进行:
步骤1、用Delany-Bazley-Miki等效流体模型表征轮胎内表面4上贴附的环形吸声材料层1;
将环形吸声材料层1等效为满足Delany-Bazley-Miki等效流体模型的流体,并利用式(1)和式(2)得到等效流体域Ω2中的复波数k2和复特征阻抗Z2:
k2=k1(1+0.122X-0.618+i0.180X-0.618) (1)
Z2=Z1(1+0.079X-0.632+i0.120X-0.632) (2)
式(1)和式(2)中:k1为在空气域Ω1中的波数,X表示中间量,且X=ρck1/2πσ,ρ为空气密度,c为空气中声速,σ为流体的流阻率,i为虚数单位,Z1为空气特征流阻,Z1=ρc。
步骤2、构建吸声材料层的轮胎三维模型,并对轮胎三维模型的表面划分单元网格;
步骤3、基于多域边界元法构建轮胎空腔声模态分析的声学模型,并基于声学模型,利用式(3)构建多域边界元特征方程:
式(3)中:A为多域边界元系数矩阵;λ表示内部贴附吸声材料轮胎的特征值;表示特征值λ对应的特征向量;式(3)中的多域边界元系数矩阵A由式(4)表征:
式(4)中:表示在边界ΓA与边界ΓA所划分的网格相关的系数矩阵,其中,/>表示复数域,NA为边界ΓA进行网格划分后的三角形单元的个数;/>和表示边界ΓA与边界ΓB所划分的网格相关的系数矩阵,NB为边界ΓB进行网格划分后的三角形单元的个数;/>表示边界ΓB与边界ΓA所划分的网格相关的系数矩阵;/>和/>表示边界ΓB与边界ΓB所划分的网格相关的系数矩阵;和/>表示边界ΓD与边界ΓC所划分的网格相关的系数矩阵,ND为边界ΓD进行网格划分后的三角形单元的个数,NC为边界ΓC进行网格划分后的三角形单元的个数;/>表示在边界ΓD与边界ΓD所划分的网格相关的系数矩阵;/>和/>表示边界ΓC与边界ΓC所划分的网格相关的系数矩阵;/>表示边界ΓC与边界ΓD所划分的网格相关的系数矩阵;ΓA表示轮胎的轮辋边界,ΓB和ΓC表示轮胎空腔中空气域Ω1与等效流体域Ω2的边界,ΓB表示空气域Ω1侧的边界,ΓC表示等效流体域Ω2侧的边界,边界ΓB和ΓC是共用边界,且NB=NC,但法向相反,ΓD表示轮胎空腔的内表面;矩阵/>Nd是求解多域边界元问题的总自由度,且Nd=NA+NB+NC+ND;并有:
式(5)-式(16)中:δlm为克罗内克δ函数,l≠m时,δlm=0,l=m时,δlm=1;G1(k1,xl,y)表示在空气域Ω1中的边界点xl与y之间的格林函数,G2(k2,xl,y)表示在等效流体域Ω2中的边界点xl与y之间的格林函数;ΓAA表示边界ΓA与边界ΓA,ΓAB表示边界ΓA与边界ΓB,ΓBA表示边界ΓA与边界ΓB,ΓBB表示边界ΓB与边界ΓB,ΓDC表示边界ΓD与边界ΓC,ΓDD表示边界ΓD与边界ΓD,ΓCC表示边界ΓC与边界ΓC,ΓCD表示边界ΓC与边界ΓD;n1(y)表示在空气域Ω1中的边界y处的外法向,n1(y)表示在等效流体域Ω2中的边界y处的外法向。
本实施例中,步骤3式(4)中矩阵A是按如下方式由基于多域边界元法的声场模型推导获得:
如图3所示,假设将轮胎内部空腔划分两个子区域Ω1和Ω2,Ω1代表空气子域,Ω2代表满足Delany-Bazley-Miki理论模型的等效流体子域,对应的是在轮胎内表面上贴附的环形吸声材料层(1)。ΓA,ΓD表示轮胎内部空腔结构两侧的实际边界,边界ΓB和ΓC是共用的等效流体域(Ω2)与空气域(Ω1)的边界。对每个子区域分别采用传统的单一区域边界元法:
式(4-1)中,下标i是子区域的索引,当i=1时表示的是空气域(Ω1),当i=2时表示的是等效流体域(Ω2);是在i子域中场点x和y之间的格林函数,r是场点x和y之间的距离,i是虚数单位,c(x)系数由点x处的几何特征决定,如果点x在结构体内,c(x)=0;如果点x在结构体外,即声场空间中,c(x)=1;如果点x在结构体表面,且表面光滑,则c(x)=0.5。qi(x)是pi(x)沿场点x的法向导数。在共用边界面ΓB,ΓC的声压与法向振速保持连续:
在轮胎内部空腔结构两侧的实际边界的实际边界ΓA,ΓD施加刚性边界条件在分界面ΓB,ΓC的应用该连续性条件即可得到式(4-4)。
式(4-4)中:表示在边界ΓA与边界ΓA所划分的网格相关的系数矩阵,/>表示边界ΓB与边界ΓA所划分的网格相关的系数矩阵;/>表示在边界ΓD与边界ΓD所划分的网格相关的系数矩阵;/>表示边界ΓC与边界ΓD所划分的网格相关的系数矩阵;/>为边界ΓA上声压组成的矩阵;/>为边界ΓB与边界ΓC上声压组成的矩阵;/>为上法向振速组成的矩阵;/>为边界ΓD上声压组成的矩阵;/>为边界ΓA上法向振速组成的矩阵;/>为边界ΓD上法向振速组成的矩阵。并有:
步骤4、使用围道积分法计算内部贴附吸声材料轮胎的特征频率λ及其对应的特征向量基于特征频率λ虚部来表征轮胎空腔5共振噪音的抑制程度,从而以特征频率λ虚部绝对值最大对应的环形吸声材料层1的流阻率σ,厚度C与宽度L作为能抑制轮胎空腔共振噪音的环形吸声材料层1的最佳流阻率、厚度和宽度。
具体实施中,围道积分法是按如下步骤进行:
步骤a、利用式(17)和式(18)建立两个Hankel矩阵,即矩阵和矩阵/>
式(17)和式(18)中,K是正整数,第l个矩阵为/>其表达式为:
式(19)中:C为复平面内的闭合积分约当曲线,T(z)为与多域边界元系数矩阵相关的中间矩阵,并且T(z)=PA(z)-1Q,矩阵和矩阵/>是两个随机非零矩阵,L是正整数且大于在复平面内所围成的区域C内包含的特征值的最大重数,且K×L大于闭合积分约当曲线C构成的区域内的总特征值数;
步骤b、利用式(20)对Hankel矩阵进行特征值分解:
式(20)中,矩阵U,是两个酉矩阵,矩阵VH是矩阵V的共轭转置,Σ是非负实对角矩阵且对角元素为σ1,σ2…σK×L;
步骤c、定义δ为正阈值,依次判断对角元素σ1,σ2…σK×L是否小于δ×σ1,若小于,则将相应对角元素置“0”;否则,保留相应对角元素;
步骤d、令Hankel矩阵求解Hankel矩阵/>的特征值γj及其特征向量ηj,并将特征值γj并转换为特征频率λj,λj=γjc/2π;
步骤e、利用式(9)计算内部贴附环形吸声材料层1的轮胎空腔5的特征频率λj对应的真实的特征向量
式(21)中,矩阵S=[S0,S1,…,Sl,…,SK-1],Sl表示第l个矩阵,并有:
通过本方法求得的内部贴附吸声材料的轮胎空腔特征频率为fe复频率,复频率fe的实部即为轮胎空腔共振频率,复频率fe的虚部绝对值的大小可以体现出轮胎空腔共振噪音的大小。复频率fe的虚部绝对值越大,轮胎空腔共振噪音越小;复频率fe的虚部绝对值越小,轮胎空腔共振噪音越大。
实施例:
以某一款型号为205/55R16的轮胎为案例,建立三维模型并划分网格,空气的密度为ρ=1.2kg/m3,空气中的声速c=343m/s。图4至图8,Re(fe)表示复频率的实部,Im(fe)表示复频率的虚部。
根据具体实施方式可知,改变环形吸声材料层1的流阻率σ,厚度C以及宽度L,可绘制出轮胎空腔特征频率fe变化的曲线:当环形吸声材料层1的厚度C=0.03m,宽度L=0.125m时,将环形吸声材料层1的流阻率σ从50Ns/m4逐渐增大至30000Ns/m4,轮胎空腔特征频率fe随环形吸声材料层1的流阻率σ的曲线如图4所示,其中点A为曲线中波谷,对应的流阻率为σ=1850Ns/m4,轮胎空腔特征频率fe=213.9-8.95i Hz。可以看出轮胎空腔特征频率fe的实部随流阻率σ增大而一直下降;同时也可以看出在流阻率在约σ=1850Ns/m4之前,轮胎空腔特征频率fe的虚部绝对值随流阻率σ增大变得越来越大;当流阻率在σ=1850Ns/m4之后,轮胎空腔共振频率fe的虚部绝对值随流阻率σ增大越来越小。
当环形吸声材料层1的流阻率σ=29100Ns/m4,宽度L=0.125m时,将环形吸声材料层1的厚度C从0.01m逐渐增大至0.06m,空腔特征频率fe的实部随环形吸声材料层1的厚度C的曲线如图5所示,空腔特征频率fe的虚部随吸声材料的厚度C的曲线如图6所示。可以从图5看出,轮胎空腔特征频率fe的实部厚度C的增大变得越来越小;可以从图6看出,轮胎空腔特征频率fe虚部的绝对自厚度C的增大变得越来越大。
当环形吸声材料层1的流阻率σ=29100Ns/m4,厚度C=0.03m时,将环形吸声材料层1的宽度L从0.0125m逐渐增大至0.175m,空腔特征频率fe实部随吸声材料的宽度L的曲线如图7所示,空腔特征频率fe虚部随环形吸声材料层1的宽度L的曲线如图8所示。可以从图7看出,轮胎空腔特征频率fe的实部宽度L的增大变得越来越小;可以从图8看出,轮胎空腔特征频率fe虚部的绝对自宽度L的增大变得越来越大。
根据由本发明方法得到的设计曲线选择环形吸声材料层1的参数和尺寸,将其粘贴在真实轮胎内,粘贴时应紧密贴合在轮胎内表面,并相对于轮胎子午线横截面中的赤道对称布置,使吸声材料层关于轮胎子午线横截面对称分布,参见图2。
为验证本发明方法,以型号为205/55R16的轮胎为试验对象,环形吸声材料层1为聚氨酯海绵,该环形吸声材料层1经测试得其流阻率σ=29100Ns/m4,横截面厚度C=0.03m,宽度L=0.125m,并符合Delany-Bazley-Miki吸声材料理论模型。一个轮胎内表面贴有两条相同的吸声材料,单条长度为0.91m,刚好紧密贴合在轮胎内表面4,为保证测试条件的相同,所有实验均用同一款轮胎和轮辋开展。
图9所示为内侧贴附环形吸声材料层1的轮胎与内侧没有贴附环形吸声材料层1的轮胎的振动传递率的对比图,曲线L1为内侧没有贴附环形吸声材料层1的轮胎的传递率曲线,曲线L2为内侧贴附环形吸声材料层1的轮胎的传递率曲线,由图9中曲线L1和L2可知,内侧贴附环形吸声材料层1使得轮胎振动传递率峰值在150Hz-250Hz空腔共振频段实现了大幅度降低,在轮胎内侧贴附环形吸声材料层1可以衰减轮胎空腔共振噪声,并且具有显著效果。因此,通过理论验证实验证实了本发明对内部贴附环形吸声材料层1的轮胎空腔共振的有效性,在有效降低轮胎空腔共振噪声的同时,也为选取抑制轮胎空腔共振噪音的环形吸声材料层1参数提供理论依据。
Claims (4)
1.一种抑制轮胎空腔共振噪音的参数设计方法,所述轮胎空腔(5)是由轮胎内表面(4)与轮辋(2a)共同形成的呈环空腔,在轮胎内表面(4)上安装吸声材料层,并沿轮胎周向整圈铺设有环形吸声材料层(1),且所述环形吸声材料层的横断面呈矩形;其特征是:所述参数设计方法是按如下步骤进行:
步骤1、用Delany-Bazley-Miki等效流体模型表征轮胎内表面(4)上贴附的环形吸声材料层(1);
步骤2、构建所述吸声材料层的轮胎三维模型,并对所述轮胎三维模型的表面划分单元网格;
步骤3、基于多域边界元法构建轮胎空腔声模态分析的声学模型,并基于所述声学模型,利用式(1)构建多域边界元特征方程:
式(1)中:A为多域边界元系数矩阵;λ表示内部贴附吸声材料轮胎的特征值;表示特征值λ对应的特征向量;
步骤4、使用围道积分法计算内部贴附吸声材料轮胎的特征频率λ及其对应的特征向量基于特征频率λ虚部来表征轮胎空腔(5)共振噪音的抑制程度,从而以特征频率λ虚部绝对值最大对应的环形吸声材料层(1)的流阻率σ,厚度C与宽度L作为能抑制轮胎空腔共振噪音的环形吸声材料层(1)的最佳流阻率、厚度和宽度。
2.根据权利要求1所述的抑制轮胎空腔共振噪音的参数设计方法,其特征在于,所述步骤1是如下过程进行:
将所述环形吸声材料层(1)等效为满足所述Delany-Bazley-Miki等效流体模型的流体,并利用式(2)和式(3)得到等效流体域Ω2中的复波数k2和复特征阻抗Z2:
k2=k1(1+0.122X-0.618+i0.180X-0.618) (2)
Z2=Z1(1+0.079X-0.632+i0.120X-0.632) (3)
式(2)和式(3)中:k1为在空气域Ω1中的波数,X表示中间量,且X=ρck1/2πσ,ρ为空气密度,c为空气中声速,σ为所述流体的流阻率,i为虚数单位,Z1为空气特征流阻,Z1=ρc。
3.根据权利要求1所述的抑制轮胎空腔共振噪音的参数设计方法,其特征是:步骤3中所述式(1)中的多域边界元系数矩阵A由式(4)表征:
式(4)中:表示在边界ΓA与边界ΓA所划分的网格相关的系数矩阵,其中,/>表示复数域,NA为边界ΓA进行网格划分后的三角形单元的个数;/>和/>表示边界ΓA与边界ΓB所划分的网格相关的系数矩阵,NB为边界ΓB进行网格划分后的三角形单元的个数;/>表示边界ΓB与边界ΓA所划分的网格相关的系数矩阵;/>和表示边界ΓB与边界ΓB所划分的网格相关的系数矩阵;/>和/>表示边界ΓD与边界ΓC所划分的网格相关的系数矩阵,ND为边界ΓD进行网格划分后的三角形单元的个数,NC为边界ΓC进行网格划分后的三角形单元的个数;
表示在边界ΓD与边界ΓD所划分的网格相关的系数矩阵;/>和表示边界ΓC与边界ΓC所划分的网格相关的系数矩阵;/>表示边界ΓC与边界ΓD所划分的网格相关的系数矩阵;ΓA表示轮胎的轮辋边界,ΓB和ΓC表示轮胎空腔中空气域Ω1与等效流体域Ω2的边界,ΓB表示空气域Ω1侧的边界,ΓC表示等效流体域Ω2侧的边界,边界ΓB和ΓC是共用边界,且NB=NC,但法向相反,ΓD表示轮胎空腔的内表面;矩阵Nd是求解多域边界元问题的总自由度,且Nd=NA+NB+NC+ND;并有:
式(5)-式(16)中:δlm为克罗内克δ函数,l≠m时,δlm=0,l=m时,δlm=1;G1(k1,xl,y)表示在空气域Ω1中的边界点xl与y之间的格林函数,G2(k2,xl,y)表示在等效流体域Ω2中的边界点xl与y之间的格林函数;ΓAA表示边界ΓA与边界ΓA,ΓAB表示边界ΓA与边界ΓB,ΓBA表示边界ΓA与边界ΓB,ΓBB表示边界ΓB与边界ΓB,ΓDC表示边界ΓD与边界ΓC,ΓDD表示边界ΓD与边界ΓD,ΓCC表示边界ΓC与边界ΓC,ΓCD表示边界ΓC与边界ΓD;n1(y)表示在空气域Ω1中的边界y处的外法向,n1(y)表示在等效流体域Ω2中的边界y处的外法向。
4.根据权利要求3所述的抑制轮胎空腔共振噪音的参数设计方法,其特征是,所述步骤4的围道积分法是按如下步骤进行:
步骤a、利用式(17)和式(18)建立两个Hankel矩阵,即矩阵和矩阵/>
式(17)和式(18)中,K是正整数,第l个矩阵为/>其表达式为:
式(19)中:C为复平面内的闭合积分约当曲线,T(z)为与多域边界元系数矩阵相关的中间矩阵,并且T(z)=PA(z)-1Q,矩阵和矩阵/>是两个随机非零矩阵,L是正整数且大于在复平面内所围成的区域C内包含的特征值的最大重数,且K×L大于闭合积分约当曲线C构成的区域内的总特征值数;
步骤b、利用式(20)对Hankel矩阵进行特征值分解:
式(20)中,矩阵U,是两个酉矩阵,矩阵VH是矩阵V的共轭转置,Σ是非负实对角矩阵且对角元素为σ1,σ2...σK×L;
步骤c、定义δ为正阈值,依次判断对角元素σ1,σ2…σK×L是否小于δ×σ1,若小于,则将相应对角元素置“0”;否则,保留相应对角元素;
步骤d、令Hankel矩阵求解Hankel矩阵/>的特征值γj及其特征向量ηj,并将特征值γj并转换为特征频率λj,λj=γjc/2π;
步骤e、利用式(9)计算内部贴附吸声材料的特征频率λj对应的真实的特征向量
式(21)中,矩阵S=[S0,S1,…,Sl,…,SK-1],Sl表示第l个矩阵,并有:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111495374.1A CN114218827B (zh) | 2021-12-09 | 2021-12-09 | 一种抑制轮胎空腔共振噪音的参数设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111495374.1A CN114218827B (zh) | 2021-12-09 | 2021-12-09 | 一种抑制轮胎空腔共振噪音的参数设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114218827A CN114218827A (zh) | 2022-03-22 |
CN114218827B true CN114218827B (zh) | 2024-02-20 |
Family
ID=80700353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111495374.1A Active CN114218827B (zh) | 2021-12-09 | 2021-12-09 | 一种抑制轮胎空腔共振噪音的参数设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114218827B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010015582A1 (fr) * | 2008-08-06 | 2010-02-11 | Societe De Technologie Michelin | Procede de creation d'un modele dynamique vibratoire d'un ensemble pneumatique roue cavite |
CN109649093A (zh) * | 2019-01-23 | 2019-04-19 | 合肥工业大学 | 一种轮胎空腔共振噪音衰减装置的参数设计方法 |
CN109649092A (zh) * | 2019-01-23 | 2019-04-19 | 合肥工业大学 | 充气轮胎空腔共振噪音消减装置的设计方法 |
-
2021
- 2021-12-09 CN CN202111495374.1A patent/CN114218827B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010015582A1 (fr) * | 2008-08-06 | 2010-02-11 | Societe De Technologie Michelin | Procede de creation d'un modele dynamique vibratoire d'un ensemble pneumatique roue cavite |
CN109649093A (zh) * | 2019-01-23 | 2019-04-19 | 合肥工业大学 | 一种轮胎空腔共振噪音衰减装置的参数设计方法 |
CN109649092A (zh) * | 2019-01-23 | 2019-04-19 | 合肥工业大学 | 充气轮胎空腔共振噪音消减装置的设计方法 |
Non-Patent Citations (1)
Title |
---|
基于有限元和边界元的轮胎振动声辐射仿真计算;曾旭;左曙光;吴旭东;倪天心;赵愿玲;;佳木斯大学学报(自然科学版);20120715(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114218827A (zh) | 2022-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109649093B (zh) | 一种轮胎空腔共振噪音衰减装置的参数设计方法 | |
JP4715902B2 (ja) | タイヤの振動特性評価方法及びタイヤの振動特性評価用コンピュータプログラム、並びにタイヤの製造方法 | |
Kindt et al. | Development and validation of a three-dimensional ring-based structural tyre model | |
CN109649092B (zh) | 充气轮胎空腔共振噪音消减装置的设计方法 | |
Yi et al. | Characteristics of sound pressure in the tire cavity arising from acoustic cavity resonance excited by road roughness | |
CN114218827B (zh) | 一种抑制轮胎空腔共振噪音的参数设计方法 | |
JP5613195B2 (ja) | タイヤのシミュレーション方法 | |
CN107531089B (zh) | 车轮 | |
Dare | Generation mechanisms of tire-pavement noise | |
JP4233473B2 (ja) | タイヤの振動特性評価方法及びタイヤの振動特性評価用コンピュータプログラム、並びにタイヤの製造方法 | |
Wan et al. | An approach for assessing the effects of porous materials on controlling the tire cavity resonance noise | |
Ni et al. | Radiated noise from tire/wheel vibration | |
Bao et al. | Design of Helmholtz resonator group in a lightweight aluminum alloy wheel for reducing tire acoustic cavity resonance noise | |
Azizi | Generation mechanisms of tire/road noise | |
US11900908B2 (en) | Wheel with a sound absorbing device | |
Krylov | Simplified analytical models for prediction of vehicle interior noise | |
Bao et al. | Numerical and experimental investigations of multiple resonators on reducing tyre cavity resonance noise | |
Kindt et al. | Three-dimensional ring model for the prediction of the tire structural dynamic behavior | |
Nakajima et al. | Tire noise | |
JP2023533308A (ja) | 車両の通過騒音を低減する方法 | |
Eisenblaetter | Experimental investigation of air related tyre/road noise mechanisms | |
JPH01115704A (ja) | 低騒音タイヤ車輪 | |
Cinkraut | Transfer Path Analysis of a Passenger Car | |
Marshall | Tire noise and vibration | |
Díaz et al. | Numerical prediction of the dynamic behaviour of rolling tyres |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |