文章信息
- 周建容, 陈春俊, 艾永军, 李新
- ZHOU Jianrong, CHEN Chunjun, AI Yongjun, LI Xin
- 高速列车明线交会流致振动耦合响应分析
- Coupled response analysis of flow-induced vibration of high-speed trains passing each other
- 中国测试, 2019, 45(5): 134-139, 150
- CHINA MEASUREMENT & TEST, 2019, 45(5): 134-139, 150
- http://dx.doi.org/10.11857/j.issn.1674-5124.2018070077
-
文章历史
- 收稿日期: 2018-07-19
- 收到修改稿日期: 2018-08-19
高速列车内部噪声传播路径包括空气声传播与固体声传播,目前高速列车多为空调列车,密封性较好,故车内噪声主要由固体声传播引起。据研究表明,在低马赫数情况下,脉动压力是气动噪声之源[1],因此研究车外气压通过二次固体传声对车内的影响机理具有重要意义。明线会车是高速列车行驶中较常见的一种工况,高速列车明线交会时,列车与空气挤压形成会车压力波[2],不仅会引起车体的横向振动,而且车厢作为典型的弹性薄壁腔体结构,会车压力波还会引起列车车壁与车厢内流场产生耦合作用,引起车厢内气压变化以及噪声增加。因此研究会车时产生的气压冲击通过车体的流致振动耦合作用对车内气压变化的影响,对探究车内噪声传播机理以及减振降噪具有重要意义。
在弹性薄壁腔体结构的声振效应方面,文献[3]基于汽车声振耦合模型,以车身振动数据为边界条件,获取了车内声学响应。文献[4]建立了装载机驾驶室以及声腔的有限元模型,通过声振耦合模型进行了频率响应分析和降噪设计。目前,声振耦合效应方面的研究主要集中在汽车、装载机、航天、板-腔系统等方面,能运用其进行噪声的预测、低噪声设计以及结构的优化设计,但是在高速列车上基于声振耦合模型研究较少。文献[5]对高速列车车厢结构进行了声振耦合模型的谐响应特性分析,但列车建模简化较多,且主要分析了列车模型的模态特性,对气压加载特性分析较少。高速列车车厢结构复杂,数值分析模型建立困难,为此可利用夹芯板等效的方法建立高速列车车身结构的等效模型,夹芯板等效最初是由Libove等[6]学者对波纹夹芯板进行了研究,后面也有不少学者投入到夹芯板等效的研究中,但是等效方法不统一,等效结果验证较困难,且基本都只考虑了位移变形等效,未进行模态等效,等效模型建立不完整。
以某线路试验高速动车组中间车厢为研究对象,分别对铝合金中空型材的车身与蜂窝板结构的车门建立列车等效模型,基于流致振动耦合有限元方法,分析车厢耦合系统的模态频率特性;通过列车交会侧壁面传感器实测的会车压力波信号,对车厢耦合系统进行气压冲击加载,在流致振动耦合作用下分析车内气压响应情况,研究气压传递规律,为列车室内的噪声研究以及减振降噪提供参考。
1 流致振动耦合分析原理随着高速列车的提速,车外气动作用加剧,车外流体的激扰将引起列车振动,而高速列车车厢作为弹性薄壁腔体结构,车体的振动又将与车内流场产生耦合作用。在封闭的结构流场中,流场会对结构振动产生影响,此时结构有限元方程[7]为
${{ M}_{\rm{s}}}\mathop {{\rm{ }}s}\limits^{ \bullet \bullet } + {{ K}_{\rm{s}}}s = {F_{\rm{s}}} + {{ A}^{\rm{T}}}p$ | (1) |
式中:
在考虑结构-流场相互耦合作用时,结构振动也会对室内流场产生影响,此时流场的有限元方程[7]为
${{ M}_{\rm{p}}}\mathop {{\rm{ }}p}\limits^{ \bullet \bullet } + {{ K}_{\rm{p}}}p + {\rho _{\rm{p}}}{ A}\mathop {{\rm{ }}s}\limits^{ \bullet \bullet } = 0$ | (2) |
式中:
求解结构-流场的耦合问题时,需要耦合求解离散后的流场有限元方程和结构有限元方程,为
$\left[ {\begin{array}{*{20}{c}} {{{ M}_{\rm{s}}}}&0 \\ {{\rho _{\rm{p}}}{ A}}&{{{ M}_{\rm{p}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\mathop s\limits^{ \bullet \bullet } } \\ {\mathop p\limits^{ \bullet \bullet } } \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{{ K}_{\rm{s}}}}&{ - {{ A}^{\rm{T}}}} \\ 0&{{{ K}_{\rm{p}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} s \\ p \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{F_{\rm{s}}}} \\ 0 \end{array}} \right]$ | (3) |
高速列车车体承载结构是以大型中空铝合金挤压型材焊接而成的双壳结构,可以减小车体质量,并且提高整体刚度性能,但该结构在数值分析中建模困难、计算量大,难以直接求解,故可通过建立等厚度的各向异性板材来等效三角夹芯板的中间夹芯层。
三角夹芯板结构见图1(a),取夹芯层的一斜板单元AB,见图1(b),以推导y方向的弹性常数为例[8],假设等效后的模型在y方向处于均匀拉伸状态,y方向的等效应力在x方向上保持不变,取x方向的一个单元长度分析,若AB板受到力
由图1受力分析有:
$M = H\frac{{{h_{\rm{c}}}}}{2},\;\;H = {\sigma _{{y}}}{h_{\rm{c}}}$ | (4) |
式中
可得板AB的弯曲挠度为
${\delta _{ m}} = \frac{{M{l^2}}}{{6E{I_{\rm c}}}} = \frac{{H{l^3}\sin \varphi }}{{E{t_{\rm c}}^3}}$ | (5) |
式中:E——夹芯板的弹性模量,Pa;
板AB的拉伸长度为
${\delta _{ n}} = \frac{{Hl\cos \varphi }}{{E{t_{\rm c}}}}$ | (6) |
则y方向的变形为
$\delta _y^H = {\delta _{ m}}\sin \varphi + {\delta _{ n}}\cos \varphi = \frac{{Hl}}{{E{t_{\rm c}}^3}}({h_{\rm c}}^2 + {\cos ^2}\varphi {t_{\rm c}}^2)$ | (7) |
可得y方向的等效弹性模量为
${E_y} = \frac{{{\sigma _y}}}{{{\varepsilon _y}}} = \frac{{{\sigma _y}}}{{\frac{{\delta _y^H}}{p}}} = \frac{{E{t_{\rm c}}^3\cos \varphi }}{{{h_{\rm c}}^3\left(1 + {{\cos }^2}\varphi {{\left(\displaystyle\frac{{{t_{\rm c}}}}{{{h_{\rm c}}}}\right)}^2}\right)}}$ | (8) |
z方向的变形为
$\delta _z^H = {\delta _{ m}}\cos \varphi - {\delta _{ n}}\sin \varphi = \frac{{Hl}}{{E{t_{\rm c}}^3}}({l^2} - {t_{\rm c}}^2)\cos \varphi \sin \varphi $ | (9) |
则泊松比为
${\upsilon _{yz}} = - \frac{{{\varepsilon _z}}}{{{\varepsilon _y}}} = \frac{{{{\cos }^2}\varphi ({l^2} - {t_{\rm c}}^2)}}{{{h_{\rm c}}^2 + {{\cos }^2}\varphi {t_{\rm c}}^2}}$ | (10) |
同理,可推导其他等效弹性常数。
根据等效前后质量守恒,等效密度为
$\rho = \frac{{{\rho _{\rm c}}{t_{\rm c}}}}{{p\sin \varphi }}$ | (11) |
式中
高速列车的门板为铝合金蜂窝夹层结构,该结构强度高,隔热隔振性能好,但也会导致列车建模、求解困难。通过将蜂窝结构的上下板与中间夹层分开,可将蜂窝夹层结构等效为一均质等厚度的正交异性板[9],从而减少计算量。
蜂窝结构单元见图2(a),以推导x方向的等效弹性常数为例,假设等效后的模型在x方向受到均质拉伸作用,x方向受到的应力为
由图2受力分析有:
$M = \frac{{P{l_1}\sin \theta }}{2},\;\;P = {\sigma _1}(h + {l_1}\sin \theta )b$ | (12) |
式中b为蜂窝芯杆件的厚度,m。
AB杆件弯曲挠度为
${\delta _1} = \frac{{P{l_1}^3\sin \theta }}{{12{E_{\rm s}}I}}$ | (13) |
式中:
杆件AB的伸长量为
${\delta _2} = \frac{{{\sigma _1}{l_1}}}{{{E_{\rm s}}}} = \frac{{P{l_1}\cos \theta }}{{{E_{\rm s}}bt}}$ | (14) |
则x方向的变形位移为
${\varepsilon _x} = \frac{{{\delta _1}\sin \theta + {\delta _2}\cos \theta }}{{{l_1}\cos \theta }} = \frac{{P{{\sin }^2}\theta {l_1}^3}}{{{E_{\rm s}}b{l_1}\cos \theta {t^3}}}\left(1 + {\cot ^2}\theta \frac{{{t^2}}}{{{l_1}^2}}\right)$ | (15) |
y方向的变形位移为
${\varepsilon _y} = \frac{{{\delta _2}\sin \theta - {\delta _1}\cos \theta }}{{h + {l_1}\sin \theta }} = - \frac{{P\sin \theta \cos \theta {l_1}^3}}{{{E_{\rm s}}b(h + {l_1}\sin \theta ){t^3}}}\left(1 - \frac{{{t^2}}}{{{l_1}^2}}\right)$ | (16) |
则泊松比为
${v_{xy}} = \frac{{{l_1}{{\cos }^2}\theta ({l_1}^2 - {t^2})}}{{(h\sin \theta + {l_1}{{\sin }^2}\theta )({l_1}^2 + {t^2}{{\cot }^2}\theta )}}$ | (17) |
x方向的弹性模量为
${E_{{\rm{c}}x}} = {E_{\rm s}}\frac{{{t^3}\cos \theta }}{{(h{{\sin }^2}\theta + {l_1}{{\sin }^3}\theta )({l_1}^2 + {{\cos }^2}\theta {t^2})}}$ | (18) |
同理可推导其他的弹性常数。
根据等效前后质量守恒,等效密度
$\rho = \frac{2}{{\sqrt 3 }}\frac{t}{{{l_1}}}{\rho _{\rm s}}$ | (19) |
式中
运用ANSYS有限元软件分别对原三角夹芯板和等效板的变形位移和模态进行求解计算,原夹芯板模型选用壳单元shell63分析,由层合板等效的等效模型采用solid185和shell63单元分析,两种模型求解得到的变形位移见图3,前4阶模态频率见表1。
参数 | 第一阶 | 第二阶 | 第三阶 | 第四阶 |
原始模型/Hz | 46.12 | 87.09 | 106.90 | 152.96 |
等效模型/Hz | 46.22 | 87.28 | 106.40 | 152.50 |
误差/% | 0.22 | 0.22 | 0.47 | 0.3 |
图3中原始模型与等效模型最大位移变形分别为0.625 mm与0.617 mm,相差1.28%,两种模型变形基本一致。表1中两种模型的模态频率最大相差0.47%,模态特性也表现一致,表明建立的车身等效模型正确。
2.3.2 车门蜂窝板等效模型验证分别运用ANSYS建立蜂窝结构的原始模型与等效模型,位移变形见图4,模态频率见表2。
参数 | 第一阶 | 第二阶 | 第三阶 | 第四阶 |
原始模型/Hz | 4 368.6 | 7 656.1 | 8 349.2 | 11 887 |
等效模型/Hz | 4 270.5 | 7 399.6 | 8 113.2 | 11 581 |
误差/% | 2.24 | 3.35 | 2.82 | 2.57 |
图4中原始模型与等效模型最大位移分别为0.271 m、0.270 m,相差0.37%,表2中前4阶模态频率最大相差3.35%,表明蜂窝结构的等效模型与原始模型表现基本一致,可以将该等效模型运用到列车建模中。
2.4 列车流致振动模型以某线路试验高速动车组中间车厢为研究对象,依次建立车厢结构、流场以及车厢结构-流场耦合系统有限元模型,模型见图5。车厢车体主要由侧墙、车顶、底架、端墙4部分组成,侧墙与车顶以铝合金中空挤压型材为主,侧墙上安装有门、窗,车门结构为铝合金蜂窝结构,车窗为多层复合结构。底架由牵引梁、枕梁、边梁和波纹地板等组成。车体侧墙和车门分别由铝合金中空型材等效模型和铝合金蜂窝结构等效模型建立,端墙结构由铝合金中空型材和梁组成,室内流场介质为空气,采用Fluid30单元,在车体与流场接触面之间建立结构-流场耦合面。
3 列车交会时耦合系统响应分析 3.1 车厢耦合系统模态分析
车厢车体壁面振动与车厢内流场会产生耦合作用,结构振动将引起车厢内流场流动,而流场又反作用于结构,两者相互制约,使得结构振动与车厢噪声被放大或抑制。分别对车厢结构、流场以及结构-流场耦合系统进行模态分析,选取前50阶模态频率进行对比,表3列举了3种模型的部分模态频率,流场的第1阶模态以及耦合系统的前7阶模态均为刚性模态。
车厢结构 | 流场 | 耦合系统 | |||||
阶数 | 频率/Hz | 阶数 | 频率/Hz | 阶数 | 频率/Hz | ||
1 | 14.907 | 2 | 6.959 | 8 | 5.978 | ||
2 | 24.999 | 3 | 13.918 | 9 | 12.026 | ||
3 | 30.714 | 4 | 20.878 | 10 | 16.498 | ||
4 | 31.993 | 5 | 27.840 | 11 | 18.153 | ||
5 | 32.490 | 6 | 34.804 | 12 | 21.947 |
3.2 会车压力波加载
高速列车明线交会时,车体与空气之间相互挤压,使列车在车厢表面产生一个移动的压力波,对车体结构与车内的舒适性均会造成很大的影响。选取某线路试验高速列车在300 km/h速度级下会车压力波作为车厢流致振动耦合模型气压加载数据。实测气压数据滤除高频干扰及气压基准,以均布载荷的方式加载至列车会车一侧,并在列车底架枕梁上施加固定约束,列车加载数据见图6。
3.3 耦合系统响应分析
选取车厢内部不同位置测点观察列车的振动特性以及车内压力的变化趋势,测点选取见图7,A点为会车侧壁板上测点,B点为会车侧车窗中点,C点为会车侧车门中点,D点为端墙上测点,E点为车厢中部测点。
不同位置的车身变形随着加载时间变化历程见图8,所选测点中车窗变形最大,车身侧壁测点与车门测点较为接近,端墙上测点变形最小。各测点的位移响应与压力波载荷趋势一致,且频率特征点相吻合。
车身结构振动引起车内气压变化,参见声压级定义,将气压信号从时域转化到频域内表示,气压转化为气压级为
$SPL = 10\lg {(p/{p_{\rm ref}})^2}$ | (20) |
式中:
观察车内气压分布情况,取车窗B的邻近测点、车厢中部测点E作为观测点,所得气压级见图9。
从气压级图中可看出,车内气压主要集中在中低频段,其中低频较大,越往高频移动,气压越小;气压级在6.1 Hz、14.67 Hz频率点处较大,6.1 Hz与耦合系统第8阶模态频率相吻合,14.67 Hz与车厢结构的第1阶模态频率相吻合,表明车厢内气压响应与耦合系统、车厢结构的固有特性有关。
4 线路试验验证选取A点侧壁处测点的车内仿真压力与该位置线路实测值进行比较,频谱图见图10,由于实际列车运行时列车空调系统、门窗等部位的泄露等,均会影响车内气压变化,导致车内气压实测频谱比仿真频谱分布更广泛,但仿真值的频谱特征点能与对应的实测频谱一致。
通过经验模态分解方法使实测气压信号依据自身成分特性进行自适应分解,得到各本征模态层,见图11,其频率依次由高向低排列。
利用相关性系数法分析各个本征模态层与仿真压力信号的相关性:
$r = \frac{{\sum\limits_{i = 1}^n {{X_j}(i)Y(i)} }}{{\sqrt {\sum\limits_{i = 1}^n {{X_j}^2(i)} } \sqrt {\sum\limits_{i = 1}^n {{Y^2}(i)} } }}$ | (21) |
式中:
相关性系数
5 结束语
本文将多重等效方法运用到高速列车的模型建立上,可以简化列车模型,减小数值模拟计算量,为列车模型的求解提供便利。建立了车厢流致振动耦合模型,分别对车厢、流场以及耦合系统进行了模态分析。以实测会车压力波信号对车厢流致振动耦合模型进行气压冲击加载,观察到车窗位移比车门与车身的变形大,端墙上的变形最小,振动位移频率分布与加载的会车压力波频率相吻合,表明车体的结构振动特性主要与列车加载载荷有关。车内气压级主要集中在低频段,高频段气压级相对较小,气压级在6.1 Hz、14.67 Hz处较大,分别与耦合系统第一阶非刚性模态频率和结构的第一阶模态频率相吻合,表明车内气压响应特性与耦合系统固有特性及车体结构固有特性有关。以实测车内压力信号与仿真信号进行对比验证,仿真信号频域基本包含在实测信号内,利用经验模态分解方法分解得到实测信号的各个本征模态层,结合相关性系数定义,实测信号第3层数据与仿真信号相关性较高,说明实测压力信号中包含有仿真信号的成分,表明所建车厢流致振动模型能基本反映车厢振动及室内气压特性,建立的模型可靠性较高。
[1] |
王亚南, 陈春俊, 何洪阳. 高速列车脉动压力的大涡模拟及小波分解[J].
机械设计与制造, 2015(8): 86-88.
DOI:10.3969/j.issn.1001-3997.2015.08.023 |
[2] |
RICCO P, BARON A, MOLTENI P. Nature of pressure waves induced by a high-speed train travelling through a tunnel[J].
Journal of Wind Engineering and Industrial Aerodynamics, 2007, 95(8): 781-808.
DOI:10.1016/j.jweia.2007.01.008 |
[3] |
侯献军, 郭金, 杜松泽, 等. 基于声模态和板件贡献分析的车身降噪研究[J].
汽车技术, 2018(5): 41-45.
|
[4] |
张俊红, 李忠鹏, 毕凤荣, 等. 基于板件贡献分析的装载机驾驶室低噪声设计[J].
振动.测试与诊断, 2016, 36(3): 568-574.
|
[5] |
徐凯, 李跃明. 高速列车车厢结构声-振耦合响应数值分析[J].
计算机辅助工程, 2011(3): 42-48.
DOI:10.3969/j.issn.1006-0871.2011.03.008 |
[6] |
LIBOVE C, BATDORF S. Elastic constants for corrugated-core sandwith plates[R].NASA TN, 1951.
|
[7] |
袁培佩.乘用车声固耦合有限元模态分析[D]. 西安: 长安大学, 2014.
|
[8] |
王青伟, 赵才其. 三角形桁架夹芯层等效弹性常数研究和夹芯板参数优化设计[J].
特种结构, 2010, 27(5): 61-66.
DOI:10.3969/j.issn.1001-3598.2010.05.015 |
[9] |
富明慧, 尹久仁. 蜂窝芯层的等效弹性参数[J].
力学学报, 1999(1): 113-118.
DOI:10.3321/j.issn:0459-1879.1999.01.014 |