1.一种基于流固耦合的高速铁路轮轨区域振动及噪声分析方法,其特征在于,包括:根据高速铁路的车轮、钢轨、轮轨的振动导纳和联合粗糙度计算出轮轨相互作用力,根据所述轮轨相互作用力基于流固耦合构建车轮—轨道振动噪声联合预测模型,利用所述车轮—轨道振动噪声联合预测模型预测得到车轮、钢轨、轨道板的振动噪声;构建轮轨区域流场的流体域模型,利用所述流体域模型得到流场内的空气动压力,根据所述空气动压力仿真得到轮轨区域的空气气动噪声;根据所述车轮、钢轨、轨道板的振动噪声和轮轨区域的空气气动噪声,基于流固耦合理论对轮轨区域的近场噪声和远场噪声进行声源贡献分析;所述的根据所述空气动压力仿真得到轮轨区域的空气气动噪声,包括:构建气动噪声的空气域预测模型,将频域下的流场内的空气动压力导入所述气动噪声的空气域预测模型,仿真得到轮轨区域的空气气动噪声;所述气动噪声的空气域预测模型的构建包括:选择湍流模型,流体材料默认选择空气,在之前设置的进口、出口、底部、上部、边界选择相应的流动面并进行了网格的划分,其中进口边界条件类型选择速度进口,速度依照实际行车速度进行选取,出口边类型选择压力出口,在建模的对称面上,在其类型选择设置为对称边界,而其他面上的边界条件类型均设置为流体不可通过的边界面,设置完毕后选择基于压力法的压力修正算法进行求解,计算步长根据最高分析频率倒数的1/2进行选取,计算结束后得到流场内的空气动压力下的时域结果;将所述空气动压力下的时域结果通过FFT变换得到频域下的流场内的空气动压;空气域中各点i的频域下的气动压力p i 与各点的气动噪声N i 间的相互关系如下式:
2.根据权利要求1所述的方法,其特征在于,所述车轮、钢轨和轮轨的振动导纳的计算方法包括:通过分别在钢轨及车轮上的某一点j点施加作用力F j ,利用谐响应分析法计算车轮及钢轨在加载点处的位移响应D j ,则车轮及钢轨加载点处的振动导纳 分别通过下式进行计算:轮轨接触导纳α contact 轮轨接触刚度K contact 的倒数,即:α contact =1/K contact (3)。
3.根据权利要求1所述的方法,其特征在于,所述车轮—轨道振动噪声联合预测模型采用有限元分析软件选择无砟轨道进行建模,钢轨采用三维实体单元模拟,扣件系统采用弹簧阻尼单元模拟,轨道板采用三维实体单元模拟,所述车轮—轨道振动噪声联合预测模型中包含一根钢轨和半块轨道板,并且通过在轨道板中间部位施加对称约束以及在下部轨道板结构结点进行全约束来模拟边界条件。
4.根据权利要求3所述的方法,其特征在于,所述的构建轮轨区域流场的流体域模型,包括:在混凝土底座以下基础上的一端的对称面的最下方轴上建立半径4.25m的流体域模型,运用布尔运算将流体域内的车轮体、钢轨、轨道板、底座和底座下基础全部舍去,只留下纯流体域模型,将流体域模型的外表面导出并进行网格划分,将外表面网格划分为进口面、出口面、对称侧面、顶面和底面,所述进口面、出口面、对称侧面的最小网格尺寸为0.01m,网格类型取三角,所述顶面和底面的最小网格尺寸为0.5m,网格类型取三角,钢轨、车轮的空气接触面的网格划分的最小网格尺寸为0.01m,网格类型取三角,基于所述流体域模型的外表面的网格划分和钢轨、车轮的空气接触面的网格划分进行3D体网格的划分。
5.根据权利要求1所述的方法,其特征在于,所述的根据所述车轮、钢轨、轨道板的振动噪声和轮轨区域的空气气动噪声,基于流固耦合理论对轮轨区域的近场噪声和远场噪声进行声源贡献分析,包括:近场噪声通过提取空气域中的各节点处的空气压力值进行换算后得到,远场噪声基于声学无限元及数值插值的算法,根据近场噪声的计算结果获取空气域最外侧表面生成的无限元网格单元j处对应的空气压力p j 进行计算,通过设置无限元封闭面形成的几何体的质心位置得到空气域以外区域内各点i距离质心位置的距离r,通过选择适合的无限元边界插值阶次n,根据下式计算得到空气域以外i处对应的空气压力p i :
6.根据权利要求2所述的方法,其特征在于,所述的根据所述车轮、钢轨、轨道板的振动噪声和轮轨区域的空气气动噪声,基于流固耦合理论对轮轨区域的近场噪声和远场噪声进行声源贡献分析,包括:对于轮轨区域的近场噪声分析,每个列车运行速度下共放置8个观测点用于观测列车的噪声变化值,每个观测点设置一个传声器,8个观测点在线路旁摆放位置包括;测点1、2是根据《GB/T 5111-2011》对于声学的轨道机车车辆发射噪声测量传声器位置放置的要求,将传声器置于距离钢轨7.5m、距轨顶面以上1.2m±0.2m,观测点1、2是环境评价指标测量点,对于两点分别采用不同列车运行速度下的A计权声压级进行噪声评价,测点3、4、5、6这四个点是横向距两轨中心轴线1m处的处于同一横截面上的观测点,在竖直垂向高度上,测点3距离距钢轨底面高0.976m,测点4距离距钢轨底面高0.576m,测点5距离距钢轨底面高0.176m即在轨顶面的高度上,测点6距离距钢轨底面高0.088m即在轨腰面的高度上,测点7、8这两个观测点是纵向钢轨底面高0.176m即在轨顶面的高度上的,而观测点5也处于该平面上,观测点7水平距离在距两轨中心轴线7.5m处,与环评测点1同纵面;观测点8水平距离在距两轨中心轴线1.5m处,与环评测点2同纵面;观测点5、7、8组成一组点用于观测不同列车运行速度下横向截面的测点噪声的变化情况。