摘要
汽车高速行驶时,侧窗附近车内外压差较大且压力脉动剧烈,因此侧窗易发生密封不良现象从而产生泄漏噪声,严重影响车内声环境品质。本文将前侧窗导槽密封条作为研究对象,采用稳态压力渗透法和瞬态动力学方法,分析高速来流下引起密封失效的条件,得出以下结论:采用稳态压力渗透法分析显示,导槽密封失效的临界泄漏静压差处于1
汽车高速行驶时,风噪成为车内噪声的主要来源,而车内乘员感受到的风噪主要来自外形噪声和泄漏噪
侧窗密封系统主要包括导槽密封条(glassrun seal)和内外水切密封条(weatherstrip seal)。由于导槽密封条附近的车内外压差更大,且压力脉动更剧烈,所以更易发生密封不良并形成泄漏路径,其对车内噪声的贡献明显,因此本文将导槽密封条作为研究泄漏噪声的对象。
以某款轿车的侧窗导槽密封条为研究对象,应用非线性有限元方法,考虑导槽密封条的橡胶材料属性、密封结构截面几何及边界条件,首先对其进行二维变形仿真及受力分析,然后探索其在外部定常和非定常压力激励条件下,结构几何变形及不同条件下密封失效从而产生声泄漏的条件。研究结果对于密封条的前期设计及预防泄漏噪声产生具有一定参考价值。
通常导槽密封条由3条不同截面形状的密封条通过接角结构连接而成,如

图1 某款经济型轿车的导槽密封条截面示意图
Fig.1 Cross-section of the glassrun seal of an economy sedan
目前,车身密封系统材料普遍采用为三元乙丙橡胶(EPDM),而EPDM可分为密实胶和海绵胶两大类。本文所研究的密封条由EPDM密实胶挤出成型。对于橡胶材料进行不同类型的分析时,需要考虑其不同的材料特性。当所需的分析类型不受时间和载荷频率影响时,即简化为准静态分析时,只需要考虑橡胶的静态特性。反之,则需考虑和时间相关的橡胶材料特
在准静态分析中,主要考虑橡胶材料的超弹性。橡胶材料作为一种典型的超弹性材料,它的应力与应变之间存在一一对应关系,根据材料的应变能函数可以求得应力。常用的应变能函数的本构模型有Mooney-Rivli
(1) |
式中:和是与温度相关的常数,需要通过试验来确定;和分别是Gauchy-Green变形张量的第一不变量和第二不变量。Mooney-Rivlin模型结构简单,得到较为广泛的应用,本文采用该模型描述材料的超弹性。
在动力学分析中,需要考虑橡胶的时间、温度和能量损耗等相关黏弹性特征,主要包括:蠕变、应力松弛、载荷速率效应、频率相关性以及温度依赖性
当侧窗玻璃插入导槽密封条后,密封条会产生变形并形成内应力。基于长度方向上各截面形变结果一致的假设,将变形问题简化为二维,应用非线性有限元方法,对不同截面处的变形仿真均基于二维平面展开。
由于导槽密封条的左侧(后视镜附近)和上侧区域的车内外压差相对较大,且这两个位置的密封条截面几何形状和边界条件均有所区别,所以需要对两个截面分别进行变形仿真和受力分析。
通过阿基米德测量密度实验,得出本文研究的导槽密封条密度大约在1 300 kg·
材料 | /MPa | /MPa | |
---|---|---|---|
0.50 | 0.615 7 | 0.307 9 | |
0.25 | 0.738 9 | 0.184 7 | |
0.05 | 0.879 6 | 0.044 0 |
完成密度和超弹性的定义后,需要定义材料的黏弹性(viscoelastic)参数。可以根据给出的时间上的松弛实验数据拟合出Prony级数方程,从而用该方程来表示材料的黏弹性特性。
本文在变形仿真阶段设置了两个分析步: 第1步是玻璃插入密封条的准静态变形仿真,通过不断尝试,将该分析步的时间设置为200 s;第二步是黏弹性松弛和蠕变的分析步,该分析步的时间和材料松弛实验数据相一致,时间设置为50 s。
本文的密封条网格类型选取平面应力四边形网格CPS4R。经网格无关性检验后,确定左侧密封条网格数量4 884个,上侧密封条网格数量4 837个。
经过玻璃插入的变形仿真以及黏弹性分析之后,左侧截面的应力分布如

图2 左侧截面应力分布(t=250 s,单位:MPa)
Fig.2 Stress distribution on the left cross-section (t=250 s, unit: MPa)
经过玻璃插入的变形仿真以及黏弹性分析之后,上侧截面的应力分布如

图3 上侧截面应力分布(t=250 s,单位:MPa)
Fig.3 Stress distribution on the upper cross-section (t=250 s, unit: MPa)
本文研究导槽密封条在稳态压差和时变压差作用下的密封性能。随着车行速度提高,车身表面内外压差增大,车内外稳态压差的最大值一般和成正
(2) |
式中:代表空气密度;U为风速;K为比值。
本文应用雷诺时间平均法Realizable k-ε的湍流模型对研究对象试验车流场进行定常分析,得到压力分布以确定K值大小,从而得到稳态激励下车内外压差最大值和风速的对应关系(考虑侧风时K的取值不同)。
车内外时变压差需通过整车非定常流体计算得到。本文采用大涡模拟方法对车外流体进行非定常计算,得到密封条表面上的时变压力,从而获得时变压差,作为后续非线性动力学分析的激励。
为求得不同工况下导槽密封条内外压差最大值以及对应位置,在导槽密封条和侧窗玻璃接触位置上布置27个监测点。根据压力脉动分布,其中左上部分布置较密,右侧较为稀疏,如

图4 导槽密封条外表面压力测点位置
Fig.4 Location of the pressure measuring points on the outer surface of the glassrun seal
风洞内测试结果显示,车内负压分布平均,0°偏航角时,120 km/h风速下为-100 Pa,160 km/h风速下为-150 Pa。应用雷诺时均法Realizable k-ε的湍流模型对试验车流场进行定常分析,可以得到无偏航角时,不同风速下27个测点上的稳态压差,如

图5 密封条不同位置处的内外压差(0°偏航角)
Fig. 5 Pressure difference between the car interior and exterior at different positions of the seal (yaw angle: 0 deg)
研究偏航角影响时,取偏航角为-20°(驾驶员一侧背风),对风速160 km/h的情况进行整车定常流场计算。27个测点上的稳态压差如

图6 密封条不同位置处的内外压差(-20°偏航角)
Fig. 6 Pressure difference between the car interior and exterior at different positions of the seal (yaw angle: -20 deg)
由
将密封条上侧测点11上的最大压差-936 Pa代入
(3) |
同理,-20°偏航角、160 km/h风速下负压最大位置为测点7,其内外压差为-4 027 Pa,得到导槽密封条内外压差最大值和车速的关系为:
(4) |
通过对整车外流场进行非定常大涡模拟计算,得到密封条表面上的时变压力,从而获得时变压差。

图7 密封条上测点11时变压差(160 km/h风速、0°偏航角)
Fig. 7 Transient pressure difference at point 11 on the sealing strip (wind speed: 160 km/h, yaw angle: 0 deg)

图8 密封条上测点7时变压差(160 km/h风速、-20°偏航角)
Fig. 8 Transient pressure difference at point 7 on the sealing strip (wind speed: 160 km/h, yaw angle: -20 deg)
为验证整车流场非定常计算结果的可靠性,通过风洞试验测试前侧窗表面上的压力脉动级。
试验在同济大学整车气动声学风洞内进行。在试验车前侧窗表面布置17个测点,每组测试采用两个表面传声器,如


图9 表面传声器及测点位置布置
Fig.9 Surface microphones and location of measuring points
由于本文主要关注导槽密封条上的压力变化,因此关注靠近导槽密封结构测点的压力脉动。

图10 试验测点11压力脉动级频谱图
Fig.10 Frequency spectra of the pressure fluctuation level on the measuring point 11
采用稳态压力渗透法计算密封界面的临界泄漏稳态压差,如

图11 弹性体密封失效的判断方法
Fig. 11 Failure assessment method of an elastomer seal
随着压力P的不断增加,接触压力也会增加,但增加程度没有流体压力大。当流体压力超过接触压力时,弹性体和硬质材料失去原有的接触区域,形成泄漏路
采用静态渗流分析的有限元变形仿真方法,对左侧截面内外两个唇边分别进行临界泄漏静压差分析。由于车内压力高于车外压力,因此静态压力作用下,密封失效的过程应该是内唇边和外唇边依次发生失效。在变形仿真结果的基础上,再进行静态流体渗透分析。不断增加内外压差时,11 000 Pa时内侧唇边先失效,内外唇边均发生失效的临界泄漏静压差为95 000 Pa。利用

图12 左侧截面两唇边与玻璃接触宽度时变曲线
Fig. 12 The transient contact width between the glass and the two lips of the left cross-section

图13 左侧截面两唇边失效示意图
Fig. 13 Failure illustration of the two lips of the left cross-section
同样,对上侧截面进行临界泄漏压差分析。不断增大内外压力差P,当P增大到7 000 Pa时,内侧唇边和玻璃内侧失去接触;当P增加到12 000 Pa时,上唇边也与玻璃失去接触;当P增加到58 000 Pa时,外唇边失去接触,3个唇边均和侧窗玻璃失去接触,形成一个连通车内车外的气吸路径。由

图14 上侧截面三唇边和玻璃接触宽度时变曲线
Fig. 14 The transient contact width between the glass and the three lips of the upper cross-section

图15 上侧截面三唇边失效示意图
Fig. 15 Failure illustration of the three lips of the upper cross-section
导槽密封条通常采用的安装方法为:密封条嵌入侧窗的导轨钢槽之中,并利用导槽密封条的突出部分卡在钢槽的凸台部分,由此将密封条固定在钢槽之内。观察分析结果的动态过程,发现施加静态渗透流体后,导槽密封条的两唇边会围绕钢槽两端向外旋转。将导槽密封条和钢槽导轨产生接触的外边全部固定,模拟导槽密封条被完全卡紧的状态,在这种边界条件下进行变形、松弛以及临界泄漏静压差的分析。固定方式的改变对两截面所有唇边的临界泄漏压差的影响如
位置 | 嵌入安装时临界 压差/Pa | 固定外边时临界 压差/Pa |
---|---|---|
内唇边 | 11 000 | 20 000 |
外唇边 | 95 000 | 120 000 |
位置 | 嵌入安装时临界 压差/Pa | 固定外边时临界 压差/Pa |
---|---|---|
内唇边 | 7 000 | 6 000 |
上唇边 | 12 000 | 15 000 |
外唇边 | 58 000 | 70 000 |
材料 | 左侧临界泄漏压差/Pa | 上侧临界泄漏压差/ Pa |
---|---|---|
M1 | 95 000 | 58 000 |
M2 | 95 000 | 45 000 |
M3 | 97 000 | 40 000 |
由于导槽密封条嵌入安装在钢槽导轨内,再与玻璃配合,所以玻璃外边缘可能存在偏离设计位置的问题。
位移/mm | 左侧外唇边临界 泄漏压差/Pa | 上侧外唇边临界 泄漏压差/Pa |
---|---|---|
0 | 95 000 | 58 000 |
0.2 | 80 000 | 52 000 |
0.5 | 55 000 | 75 000 |
0.7 | 40 000 | 55 000 |
1.0 | 1 000 | 50 000 |
1.2 | 已失去接触 | 49 000 |
1.5 | 已失去接触 | 40 000 |
1.7 | 已失去接触 | 35 000 |

图16 玻璃位移对外唇边临界泄漏压差的影响
Fig. 16 The critical static pressure difference of the outer lips of the seal under different inward offsets of the glass
由以上可知,即使在橡胶硬化或玻璃有一定的位置偏差时,外唇边的临界泄漏静压差的数量级普遍在1
试验车在无偏航及160 km/h风速下,其导槽密封条的上侧截面所承受的内外压差最大,因此上侧截面是无偏航行驶时密封失效的“危险位置”。将整车非定常CFD计算求得测点11位置上随时间变化的车内外压力差作为瞬态动力学分析的激励,其时间步长为0.000 2 s,其0.5 s内测点11位置的压力差随时间变化的情况如

图17 上侧截面时变压力施加位置
Fig. 17 Excitation position of transient pressure on the upper cross-section
当玻璃处于基准设计位置,时变压力激励下外唇边和玻璃之间的接触宽度仅有小幅降低并迅速稳定,如

图18 上侧截面外唇边接触宽度时变曲线(玻璃基准位置,160 km/h风速、0°偏航角)
Fig. 18 The transient contact width of the outer lips of the upper cross-section (glass position: baseline, wind speed: 160km/h, yaw angle: 0 deg)
不断将玻璃向内移动做相同分析,将玻璃向内移动1.7 mm(保持外唇边接触最大偏移量)得到外唇边和玻璃的接触宽度如

图19 上侧截面外唇边接触宽度时变曲线(玻璃位移1.7 mm,160 km/h风速、0°偏航角)
Fig. 19 The transient contact width of the outer lips of the upper cross-section (glass position: 1.7mm inward offset, wind speed: 160km/h, yaw angle: 0deg)
因此可以初步得出结论:对于本文的研究对象而言,在160 km/h风速、无偏航条件下的时变压差作用下,导槽密封条的上侧截面始终保持密封,甚至在玻璃向内偏离设计位置较大时,外唇边依然可以保持密封状态。由此可以推断整条密封条各个位置均可以保持密封状态。
研究对象车辆在-20°偏航角、160 km/h风速下,其导槽密封条的左侧截面所承受的内外压差最大,因此左侧截面是-20°偏航行驶时密封失效的“危险位置”。同样,将测点7位置上随时间变化的车内外压力差作为激励,其时间步长0.000 2 s,取0.5 s内的数据,则测点7位置的压力差随时间变化的情况如

图20 左侧截面时变压力施加位置
Fig. 20 Excitation position of transient pressure on the left cross-section
此时玻璃处于基准设计位置,分析后得到时变压力激励下外唇边和玻璃之间的接触宽度先降低并随后稳定,如

图21 左侧截面外唇边接触宽度时变曲线(玻璃基准位置,160 km/h风速、-20°偏航角)
Fig. 21 The transient contact width of the outer lips of the left cross-section (glass position: baseline, wind speed: 160 km/h, yaw angle: -20 deg)
将玻璃相对于基准设计位置向内偏移0.2 mm,时变压差激励下外唇边和玻璃之间的接触宽度随时间不断变化,且在零和非零之间不断振荡,如

图22 左侧截面外唇边接触宽度时变曲线(玻璃位移0.2 mm,160km/h风速、-20°偏航角)
Fig. 22 The transient contact width of the outer lips of the left cross-section (glass position: 0.2 mm inward offset, wind speed: 160 km/h, yaw angle: -20 deg)
(1) 在稳态流场中,无侧风时的导槽密封条上最大压差位置在上侧,-20°偏航角时左侧截面为压差最大位置。
(2) 静态分析中,左侧截面的临界泄漏压差为95 000Pa,对应-20°偏航角下的临界车速为778 km/h;上侧截面的临界泄漏压差为58 000 Pa,对应-20°偏航角下的临界车速为607 km/h。可见静态分析下导槽密封条在实际驾驶情况中极难失效。此外边界条件、本构模型中参数的比值、材料的硬度和玻璃的相对位置都会对唇边的临界泄漏静压差产生影响。
(3) 非线性动力学分析中,在无偏航及风速160 km/h情况下,导槽密封条的上侧截面是密封失效的危险位置。在0°偏航角、160 km/h风速带来的时变压差作用下,整个密封条可保持密封状态。
(4) 非线性动力学分析中,在-20°偏航角、风速160 km/h情况下,导槽密封条的左侧截面是密封失效的危险位置。该工况的时变压差作用下,玻璃处于基准设计位置时导槽密封条依然保持密封状态;而玻璃相对基准位置向内移动0.2 mm后,左侧截面处产生密封失效现象。由此可见左侧位置在-20°偏航角、风速160km/h情况下比较容易发生失效。
(5) 对密封条失效的临界条件(失效压力)分析时,采用静态分析和瞬态动力学分析得到的结果相差较大。由于采用静态分析的结果与现实不符,应采用瞬态动力学分析方法,才可能对密封条发生失效并产生泄漏噪声的临界条件进行更准确地预估。
参考文献
GEORGE A R. Automobile aerodynamic noise[J]. SAE Transactions, 1990, 99:434. [百度学术]
AHMED S R. A survey of automobile aeroacoustic activities in Germany[J]. SAE Transactions. 1995, 104:1100. [百度学术]
BUCHHEIM R, DOBRZYNSKI W, MANKAU H, et al. Vehicle interior noise related to external aerodynamics[J]. International Journal of Vehicle Design, 1982, 3(4): 398. [百度学术]
LOREA A, CASTELLUCCIO V, COSTELLI A, et al. A wind-tunnel method for evaluating the aerodynamic noise of cars[J]. SAE Transactions, 1986, 95:1096. [百度学术]
GUR Y, MORMAN K N, SINGH N. Analysis of door and glass run seal systems for aspiration[C]//Proceedings of the 1997 Noise and Vibration Conference-P-309. Detroit: SAE, 1997: 971902. [百度学术]
JUNG W W, OH S J. The influence of vehicle elements to aspiration wind noise[C]//International Congress & Exposition. Detroit: SAE, 1995: 950624. [百度学术]
SAHA P, MYERS R D. Importance of sealants for interior noise control of automobiles[C]//International Congress & Exposition. Detroit: SAE, 1992: 920412. [百度学术]
GOTO T, KAWAI K. Performance improvement in leak noise reduction[C]//SAE 2002 World Congress & Exhibition: Detroit: SAE, 2002: 2002-01-1238. [百度学术]
贺银芝, 杨志刚, 王毅刚. 汽车车身密封对车内气动噪声影响的机理及试验研究[J]. 汽车工程, 2012, 34(8):692. [百度学术]
HE Yinzhi, YANG Zhigang, WANG Yigang. Mechanism and experimental research on the effect of automobile body seal on aerodynamic noise in vehicles[J]. Automotive Engineering, 2012, 34(8): 692. [百度学术]
赵建才, 姚振强, 王伟. 密封条结构的非线性分析与优化设计[J]. 合成橡胶工业, 2005,28(6):421. [百度学术]
ZHAO Jiancai, YAO Zhenqiang, WANG Wei. Non-linear analysis and optimization design of seal structure[J]. Synthetic Rubber Industry, 2005, 28(6): 421. [百度学术]
MOONEY M. A theory of large elastic deformation[J]. Journal of Applied Physics, 2004, 11(9): 582. [百度学术]
RIVLIN R S, SAUNDERS D W. Large elastic deformations of isotropic materials VII experiments on the deformation of rubber[J]. Philosophical Transactions of the Royal Society of London, 1951, 243(865): 251. [百度学术]
左亮. 机车车辆中常用橡胶件的有限元分析[D]. 成都:西南交通大学, 2008. [百度学术]
ZUO liang. Finite element analysis of rubber parts used in locomotives and vehicles[D]. Chengdu: Southwest Jiaotong University, 2008. [百度学术]
LIU Q, WANG Z, LOU Y, et al. Elastic leak of a seal[J]. Extreme Mechanics Letters, 2014, 1: 54. [百度学术]