CN107292046B - 一种减振降噪效果的检验方法及装置 - Google Patents
一种减振降噪效果的检验方法及装置 Download PDFInfo
- Publication number
- CN107292046B CN107292046B CN201710535321.5A CN201710535321A CN107292046B CN 107292046 B CN107292046 B CN 107292046B CN 201710535321 A CN201710535321 A CN 201710535321A CN 107292046 B CN107292046 B CN 107292046B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- node
- calculating
- track
- matrix
- 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 84
- 230000009467 reduction Effects 0.000 title claims abstract description 60
- 230000000694 effects Effects 0.000 title claims abstract description 44
- 238000007689 inspection Methods 0.000 title claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims abstract description 85
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 238000012546 transfer Methods 0.000 claims abstract description 46
- 230000003595 spectral effect Effects 0.000 claims abstract description 37
- 238000013016 damping Methods 0.000 claims abstract description 17
- 230000001603 reducing effect Effects 0.000 claims abstract description 5
- 229910000831 Steel Inorganic materials 0.000 claims description 61
- 239000010959 steel Substances 0.000 claims description 61
- 238000004364 calculation method Methods 0.000 claims description 36
- 238000006073 displacement reaction Methods 0.000 claims description 35
- 230000000737 periodic effect Effects 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 17
- 230000001419 dependent effect Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 abstract description 10
- 238000011160 research Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000005096 rolling process Methods 0.000 abstract description 5
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 230000001629 suppression Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 15
- 238000005452 bending Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 8
- 230000003321 amplification Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000003199 nucleic acid amplification method Methods 0.000 description 6
- 230000007423 decrease Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 241001669679 Eleotris Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 235000010931 Mesua ferrea Nutrition 0.000 description 1
- 235000005704 Olneya tesota Nutrition 0.000 description 1
- 240000007909 Prosopis juliflora Species 0.000 description 1
- 235000008198 Prosopis juliflora Nutrition 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000001131 transforming effect Effects 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
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01B—PERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
- E01B19/00—Protection of permanent way against development of dust or against the effect of wind, sun, frost, or corrosion; Means to reduce development of noise
- E01B19/003—Means for reducing the development or propagation of noise
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Railway Tracks (AREA)
Abstract
本发明实施例提供一种减振降噪效果的检验方法及装置,涉及轨道减振技术领域。通过将谱元法和谱传递矩阵法应用到了轨道动力学研究中,采用整体道床轨道作为研究结构,建立平面半轨道模型,利用谱元法计算减振扣件作用下的力传递率。同时,对相同轨道结构建立轨道周期子结构模型,利用谱传递矩阵法计算轨道衰减率,并结合钢轨导纳仿真结果,得到钢轨在单位简谐点激励作用下的声功率级。依据声功率级和力传递率可对扣件的减振效果进行分析,能够为精确、高效和全面的评判和检验扣件的减振降噪效果,对降低铁路轨道的滚动噪声有着重要意义。
Description
技术领域
本发明涉及轨道减振技术领域,具体而言,涉及一种减振降噪效果的检验方法及装置。
背景技术
国内轨道减振技术通常采用隔振原理,即在组成轨道的各个刚性部件之间***低刚度弹性层,隔离振动能量往轨下基础的传递,可减小地面振动以及降低建筑物的二次噪声。轨下的扣件支承刚度的改变将引起振动能量的重分配,被减振措施所隔离的振动能量将反射到上方的***中,增大滚动噪声,影响高架地铁沿线居民的正常生活及乘客乘车舒适度。因此,评判扣件的减振降噪效果,对降低铁路的滚动噪声有着重要意义。
发明内容
有鉴于此,本发明实施例的目的在于提供一种减振降噪效果的检验方法及装置,更为精确、高效和全面的扣件的减振降噪效果。
本发明实施例提供一种减振降噪效果的检验方法,所述方法包括:建立平面半轨道模型,利用谱元法计算力传递率;建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率;依据所述轨道衰减率计算钢轨声功率;依据所述力传递率和所述钢轨声功率检验扣件的减振降噪效果。
本发明实施例还提供一种减振降噪效果的检验装置,所述装置包括:力传递计算模块,用于建立平面半轨道模型,利用谱元法计算力传递率;轨道衰减率计算模块,用于建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率;钢轨声功率计算模块,用于依据所述轨道衰减率计算钢轨声功率;检验模块,用于依据所述力传递率和所述钢轨声功率检验扣件的减振降噪效果。
与现有技术相比,本发明的实施例提供的减振降噪效果的检验方法及装置,通过将谱元法和谱传递矩阵法应用到了轨道动力学研究中,采用整体道床轨道作为研究结构,建立平面半轨道模型,利用谱元法计算减振扣件作用下的力传递率。同时,对相同轨道结构建立轨道周期子结构模型,利用谱传递矩阵法计算轨道衰减率,并结合钢轨导纳仿真结果,得到钢轨在单位简谐点激励作用下的声功率级。依据声功率级和力传递率可对扣件的减振效果进行分析,能够为精确、高效和全面的评判和检验扣件的减振降噪效果,对降低铁路轨道的滚动噪声有着重要意义。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明第一实施例提供的电子设备的方框示意图。
图2为本发明第一实施例提供的减振降噪效果的检验装置的功能模块示意图。
图3为本发明第二实施例提供的减振降噪效果的检验方法的流程图。
图4为本发明第二实施例提供的平面半轨道模型的示意图。
图5为本发明第二实施例提供的谱元法计算模型的示意图。
图6为本发明第二实施例提供的轨道周期子结构模型的示意图。
图7是本发明第二实施例提供的谱传递矩阵法计算模型的示意图。
图8是本发明第二实施例提供的通过仿真实验得到的钢轨位移导纳与频率的关系图。
图9是本发明第二实施例提供的不同扣件刚度下力传递率与频率的关系图。
图10是本发明第二实施例提供的不同扣件刚度下钢轨相对声功率级与频率的关系图。
图标:100-电子设备;101-存储器;102-存储控制器;103-处理器;200-减振降噪效果的检验装置;201-力传递计算模块;202-轨道衰减率计算模块;203-钢轨声功率计算模块;204-检验模块。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本发明的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
第一实施例
请参考图1,是发明第一实施例提供的电子设备100的方框示意图。本发明实施例提供的减振降噪效果的检验装置200可应用于电子设备100中,所述电子设备100可检验计算轨道的力传递率、轨道衰减率、钢轨声功率等参数,以检验扣件的减振降噪的效果。一般的,铁路轨道包括:钢轨、轨枕、扣件和道床,钢轨通过轨枕铺设在道床上,扣件设置在钢轨下,用于隔振。
所述电子设备100可以是,但不限于智能手机、个人电脑(personal computer,PC)、平板电脑、个人数字助理(personal digital assistant,PDA)、移动上网设备(mobileInternet device,MID)等。所述电子设备100的操作***可以是,但不限于,安卓(Android)***、IOS(iPhone operating system)***、Windows phone***、Windows***等。
于本发明实施例中,所述电子设备100还包括存储器101、存储控制器102、处理器103。所述存储器101、存储控制器102、处理器103,各元件相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述减振降噪效果的检验装置200包括至少一个可以软件或固件(firmware)的形式存储于所述存储器101中或固化在所述电子设备100的操作***(operating system,OS)中的软件功能模块。所述处理器103用于执行存储器101中存储的可执行模块,例如,所述减振降噪效果的检验装置200包括的软件功能模块或计算机程序。
处理器103可能是一种集成电路芯片,具有信号的处理能力。上述的处理器103可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器103也可以是任何常规的处理器等。
请参考图2,是本发明第一实施例提供的减振降噪效果的检验装置200的功能模块示意图。所述减振降噪效果的检验装置200通过对谱元法进行推导,对谱传递矩阵法进行了改进,并将两者应用到了轨道动力学研究中。采用整体道床轨道作为研究结构,建立平面半轨道模型,利用谱元法计算减振扣件作用下的力传递率。同时,对相同轨道结构建立轨道周期子结构模型,利用谱传递矩阵法计算轨道衰减率,并结合钢轨导纳仿真结果,得到钢轨在单位简谐点激励作用下的声功率级。依据声功率级和力传递率可对扣件的减振效果进行分析,以便于在实际施工中选取合适刚度的扣件。
所述减振降噪效果的检验装置200包括力传递计算模块201、轨道衰减率计算模块202、钢轨声功率计算模块203以及检验模块204。
所述力传递计算模块201,用于建立平面半轨道模型,利用谱元法计算力传递率。具体的,以两相邻轨枕之间的钢轨作为Timoshenko梁单元;利用Timoshenko梁运动方程,计算Timoshenko梁单元的谱节点荷载和谱节点位移的关系;依据Timoshenko梁单元的谱节点荷载和谱节点位移的关系,计算平面半轨道模型的谱节点荷载与谱节点位移的关系;依据平面半轨道模型的谱节点荷载与谱节点位移的关系计算力传递率。
Timoshenko梁也称为铁木辛柯梁,剪应变在一个给定横截面上是常值。接着引入剪切校正因子来解释这种简化,其值取决于横截面的形状。在横向剪切的存在下,横截面的旋转就由挠曲和横向(平面外)剪变形引起。
所述轨道衰减率计算模块202,用于建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率。
具体的,利用谱传递矩阵法对轨道周期子结构模型进行求解,并建立四节点计算模型。建立关于四节点计算模型的谱节点荷载与谱节点位移的关系:
将变换后的整体谱刚度矩阵划分为3×3子矩阵形式:
其中,
由于四节点的中间节点不受外荷载,因此,可得
变换为输入-输出的形式为:
其中,
t(ω)为谱传递矩阵,为
对于一维周期结构,由Bloch定理得到:
得到所述4×4矩阵特征值问题:
其中,为节点i的谱节点位移向量,为节点i的谱节点荷载向量。
即可得到波矢与频率之间的频散关系;依据所述波矢并通过Δ=-8.686Im(kx)求解所述轨道衰减率,其中,Δ为所述轨道衰减率,kx为一维Bloch波矢。
钢轨声功率计算模块203,用于依据所述轨道衰减率计算钢轨声功率。
具体的,通过公式
|v(x)|=v(0)e-β|x|,
计算钢轨的声功率;
依据钢轨的声功率,计算声功率级:
其中,v(x)是钢轨在x处振动速度的幅值,v(0)为原点处的振动速度幅值,β为衰减系数,ρ0c0为空气中的声特性阻抗,ρ0=1.225kg/m3为空气密度,c0=340m/s为声波在空气中的传播速度,σ是与频率有关的辐射率,为钢轨的固有属性,P为一个截面的周长,W0为基准声功率,W0=10-12W。
第二实施例
请参考图3,是本发明第二实施例提供的减振降噪效果的检验方法的流程图。本实施例的减振降噪效果的检验方法应用于第一实施例的减振降噪效果的检验装置200,所述减振降噪效果的检验方法包括:
步骤S101,建立平面半轨道模型,利用谱元法计算力传递率。
请参考图4,是本发明第二实施例提供的平面半轨道模型的示意图。本实施例中,平面半轨道模型的长度为25m,轨枕间距a采用现场实测值62.5cm。钢轨采用Timoshenko梁模型,由于钢轨在1500Hz频率以上发生明显的截面变形,频率上限为1500Hz。为了减小弹性波反射对计算精度的影响,在钢轨两端建立无反射边界条件。扣件的垂向限位能力用一组拉伸弹簧和扭转弹簧来模拟:垂向刚度为kv,扭转刚度为krot。考虑扣件的阻尼特性,刚度采用复刚度形式:
kcv=kv(1+iηp),其中ηp为轨下垫板阻尼损耗因子。同时,垂向刚度和扭转刚度满足如下的关系表达式:其中,br为轨下垫板沿纵向的长度。
请参考图5,是本发明第二实施例提供的谱元法计算模型的示意图。计算Timoshenko梁单元的谱节点荷载和谱节点位移的关系,以两相邻轨枕之间的钢轨作为Timoshenko梁单元,对于材料均匀一致、等截面、自由振动的Timoshenko梁运动方程如下式所示:
其中,G为剪切弹性模量,A为梁横截面面积,K为剪切矫正因子,v为梁垂向位移,φ为弯曲引起的变形角,ρ为材料密度,E为弹性模量,Y为横截面的惯性矩。内力和位移有如下的关系:
其中,V和M分别为梁内部剪力和弯矩。
Timoshenko梁运动方程表示为谱分量的形式为:
其中,和为谱分量,N为时域采样点数,ωn为圆频率。
经过求解可得,梁单元的边界条件为:
其中,为谱节点位移向量,结合上式可得,
其中,为[4×4]矩阵,
r1=(R1-R2)(1-e1e2)
r2=(R1+R2)(e1-e2)
谱元法Timoshenko梁单元和的形函数分别为和显然该形函数与频率相关。
由公式,
当x=0和L时,
其中,为谱节点荷载向量,而
由于谱元法和材料力学的符号正负规定在梁左端相反,故和前存在负号。因此,谱节点荷载和位移的关系可表示为:
其中,为Timoshenko梁的谱单元刚度矩阵,同样,St(ω)与频率相关。
依据Timoshenko梁单元的谱节点荷载和谱节点位移的关系,计算平面半轨道模型的谱节点荷载与谱节点位移的关系。
采用与有限元相同的坐标变换方式,将单元坐标系下的谱单元刚度矩阵变换到全局坐标系内,再组装到整体谱刚度矩阵中,同时处理模型的约束条件。整体模型的谱节点荷载与位移的关系为:
其中,为整体模型的谱节点荷载向量,为整体模型的谱节点位移向量,S(ω)为整体谱刚度矩阵。
通过对整体模型的谱节点荷载与位移的关系式求解,可计算轨道的力传递率。
本发明实施例中,所述步骤S101可通过力传递计算模块201执行。
步骤S102,建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率。请参考图6,是本发明第二实施例提供的轨道周期子结构模型的示意图。轨道周期子结构模型的长度为a,轨道部件的建模与平面半轨道模型相同,而钢轨两端满足Bloch定理。
利用谱传递矩阵法对轨道周期子结构模型进行求解,建立了谱传递矩阵法计算模型(四节点计算模型),请参考图7,是本发明第二实施例提供的谱传递矩阵法计算模型的示意图。
建立关于四节点计算模型的谱节点荷载与谱节点位移的关系:
将变换后的整体谱刚度矩阵划分为3×3子矩阵形式:
其中,
由于四节点的中间节点不受外荷载,因此,可得变换为输入-输出的形式为:
其中,t(ω)为谱传递矩阵,为
对于一维周期结构,由Bloch定理得到:
得到所述4×4矩阵特征值问题:
其中,为节点i的谱节点位移向量,为节点i的谱节点荷载向量。
通过求解特征值问题,即可得到波矢kx与频率ω之间的频散关系。对于本模型中的钢轨垂向弯曲问题,Bloch波矢的解总以±kx的形式成对出现,分别表示同一种特征波沿相反方向的传播特性,轨道周期子结构模型的解包含两对特征波。
轨道衰减率可通过表达式Δ=-8.686Im(kx)求得,即轨道衰减率为波数虚部的-8.686倍。
本发明实施例中,所述步骤S102可通过轨道衰减率计算模块202执行。
步骤S103,依据所述轨道衰减率计算钢轨声功率。
通过公式
|v(x)|=v(0)e-β|x|,
计算钢轨的声功率;
依据钢轨的声功率,计算声功率级:
其中,v(x)是钢轨在x处振动速度的幅值,v(0)为原点处的振动速度幅值,β为衰减系数,ρ0c0为空气中的声特性阻抗,ρ0=1.225kg/m3为空气密度,c0=340m/s为声波在空气中的传播速度,σ是与频率有关的辐射率,为钢轨的固有属性,P为一个截面的周长,W0为基准声功率,W0=10-12W。上式表示,轨道衰减率越大,钢轨声功率级越低,钢轨原点处的速度幅值越大,钢轨声功率级越高。
本发明实施例中,所述步骤S103可通过声功率计算模块执行。
步骤S104,依据所述力传递率和所述钢轨声功率检验扣件的减振降噪效果。
请参考图8,是通过仿真实验得到的钢轨位移导纳与频率的关系图。钢轨垂向原点位移导纳幅值曲线在96Hz处出现波峰,此时钢轨发生垂向弯曲共振,共振频率记为fr。该频率的大小主要由扣件垂向刚度决定。高于fr频率,位移导纳幅值随着频率的增大而减小。当频率为1053Hz时,跨中钢轨垂向原点位移导纳幅值曲线出现峰值,而扣件处幅值曲线出现波谷,钢轨发生一阶垂向弯曲pinned-pinned共振(频率记为fp),其共振模态对应于扣件处为节点的驻波。钢轨跨中处原点位移导纳幅值与扣件处原点位移导纳幅值基本一致,仅在fp频率处有所不同。
请参考图9,是不同扣件刚度下力传递率与频率的关系图。不同扣件刚度下的振动放大频带分别为:67~110Hz(10kN/mm),91~163Hz(20kN/mm),109~207Hz(30kN/mm),123~245Hz(40kN/mm)。随着扣件刚度的增大,振动放大频带往高频移动,并且宽度加宽。当频率高于振动放大频带范围,扣件的减振性能随着扣件刚度的降低而增大,且扣件刚度越低,相同扣件刚度减少量下的减振效果增幅越大。在300~1000Hz范围内,扣件刚度为10kN/mm的力传递率比扣件刚度为40kN/mm的力传递率低11dB左右。
请参考图10,是不同扣件刚度下钢轨相对声功率级与频率的关系图。在扣件刚度减小后的相应fr频率处,钢轨相对声功率级曲线出现峰值,钢轨声功率级比刚度减小前的声功率级前明显增大,且扣件刚度减小量越多,钢轨相对声功率级曲线峰值越高。这是因为钢轨导纳幅值曲线的峰值所在的频率fr随着扣件刚度的减小而降低,而轨道衰减率也随着扣件刚度的减小而降低,尤其在减小后的频率fr处,轨道衰减率降低明显。当中心频率为80Hz时,扣件刚度为10kN/mm的钢轨声功率级比扣件刚度为40kN/mm的钢轨声功率级高20dB。与此同时,力传递率在频率fr附近大于0dB,扣件对振动有放大的作用,因此,钢轨弯曲共振对减振降噪均为不利。在中心频率200Hz以上,尽管扣件刚度对钢轨导纳幅值的影响较小,但轨道衰减率仍受扣件刚度的影响。随着扣件刚度的减小,钢轨声功率级和扣件减振效果均增大,被隔离的振动能量沿着钢轨纵向传递,并以声能的形式向外辐射。在中心频率315Hz处,扣件刚度为10kN/mm的钢轨声功率级比扣件刚度为40kN/mm的钢轨声功率级高3.7dB。
随着扣件刚度的增大,振动放大频带往高频移动,并且宽度加宽。当频率高于振动放大频带范围,扣件的减振性能随着扣件刚度的降低而增大,且扣件刚度越低,相同扣件刚度减少量下的减振效果增幅越大。在扣件刚度减小后的相应钢轨弯曲共振频率处,钢轨相对声功率级曲线出现峰值,钢轨声功率级比刚度减小前的声功率级前明显增大,且扣件刚度减小量越多,钢轨相对声功率级曲线峰值越高。钢轨弯曲共振对减振降噪均为不利。在中心频率200Hz以上,随着扣件刚度的减小,钢轨声功率级和扣件减振效果均增大,被隔离的振动能量沿着钢轨纵向传递,并以声能的形式向外辐射。
本发明实施例中,所述步骤S104可通过检验模块204执行。
综上所述,本发明实施例提供的减振降噪效果的检验方法及装置,通过将谱元法和谱传递矩阵法应用到了轨道动力学研究中,采用整体道床轨道作为研究结构,建立平面半轨道模型,利用谱元法计算减振扣件作用下的力传递率。同时,对相同轨道结构建立轨道周期子结构模型,利用谱传递矩阵法计算轨道衰减率,并结合钢轨导纳仿真结果,得到钢轨在单位简谐点激励作用下的声功率级。依据声功率级和力传递率可对扣件的减振效果进行分析,能够为精确、高效和全面的评判和检验扣件的减振降噪效果,对降低铁路轨道的滚动噪声有着重要意义。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的***来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (6)
1.一种减振降噪效果的检验方法,其特征在于,所述方法包括:
建立平面半轨道模型,利用谱元法计算力传递率;
建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率;
依据所述轨道衰减率计算钢轨声功率;
依据所述力传递率和所述钢轨声功率检验扣件的减振降噪效果;
其中,所述建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率的步骤具体包括:
利用谱传递矩阵法对轨道周期子结构模型进行求解,并建立四节点计算模型;
利用谱元法求解计算四节点计算模型,得到4×4矩阵特征值问题;
求所述4×4矩阵特征值问题,即可得到波矢与频率之间的频散关系;
依据所述波矢并通过Δ=-8.686Im(kx)求解所述轨道衰减率,其中,Δ为所述轨道衰减率,kx为一维Bloch波矢,Im(kx)为一维Bloch波矢的虚部;
所述利用谱元法求解计算四节点计算模型,得到4×4矩阵特征值问题的步骤具体包括:
建立关于四节点计算模型的谱节点荷载与谱节点位移的关系:
将变换后的整体谱刚度矩阵划分为3×3子矩阵形式:
其中,
由于四节点的中间节点不受外荷载,因此,可得
变换为输入-输出的形式为:
其中,
t(ω)为谱传递矩阵,为
对于一维周期结构,由Bloch定理得到:
得到所述4×4矩阵特征值问题:
其中,为节点i的谱节点位移向量,为节点i的谱节点荷载向量,S为整体谱刚度矩阵中的矩阵元素,a为轨道周期子结构模型的长度,I为单元矩阵。
2.根据权利要求1所述的减振降噪效果的检验方法,其特征在于,所述建立平面半轨道模型,利用谱元法计算力传递率的步骤具体包括:
以两相邻轨枕之间的钢轨作为Timoshenko梁单元;
利用Timoshenko梁运动方程,计算Timoshenko梁单元的谱节点荷载和谱节点位移的关系;
依据Timoshenko梁单元的谱节点荷载和谱节点位移的关系,计算平面半轨道模型的谱节点荷载与谱节点位移的关系;
依据平面半轨道模型的谱节点荷载与谱节点位移的关系计算力传递率。
3.根据权利要求1所述的减振降噪效果的检验方法,其特征在于,所述方法包括:
通过公式
|v(x)|=v(0)e-β|x|,
计算钢轨的声功率;
依据钢轨的声功率,计算声功率级:
其中,v(x)是钢轨在x处振动速度的幅值,v(0)为原点处的振动速度幅值,β为衰减系数,ρ0c0为空气中的声特性阻抗,ρ0=1.225kg/m3为空气密度,c0=340m/s为声波在空气中的传播速度,σ是与频率有关的辐射率,为钢轨的固有属性,P为一个截面的周长,W0为基准声功率,W0=10-12W,W为钢轨的声功率,Lw为声功率级。
4.一种减振降噪效果的检验装置,其特征在于,所述装置包括:
力传递计算模块,用于建立平面半轨道模型,利用谱元法计算力传递率;
轨道衰减率计算模块,用于建立轨道周期子结构模型,利用谱传递矩阵法和谱元法计算轨道衰减率;
钢轨声功率计算模块,用于依据所述轨道衰减率计算钢轨声功率;
检验模块,用于依据所述力传递率和所述钢轨声功率检验扣件的减振降噪效果;
所述轨道衰减率计算模块具体用于:
利用谱传递矩阵法对轨道周期子结构模型进行求解,并建立四节点计算模型;
利用谱元法求解计算四节点计算模型,得到4×4矩阵特征值问题;
求所述4×4矩阵特征值问题,即可得到波矢与频率之间的频散关系;
依据所述波矢并通过Δ=-8.686Im(kx)求解所述轨道衰减率,其中,Δ为所述轨道衰减率,kx为一维Bloch波矢,Im(kx)为一维Bloch波矢的虚部;
所述轨道衰减率计算模块还用于:
建立关于四节点计算模型的谱节点荷载与谱节点位移的关系:
将变换后的整体谱刚度矩阵划分为3×3子矩阵形式:
其中,
由于四节点的中间节点不受外荷载,因此,可得
变换为输入-输出的形式为:
其中,
t(ω)为谱传递矩阵,为
对于一维周期结构,由Bloch定理得到:
得到所述4×4矩阵特征值问题:
其中,为节点i的谱节点位移向量,为节点i的谱节点荷载向量,S为整体谱刚度矩阵中的矩阵元素,a为轨道周期子结构模型的长度,I为单元矩阵。
5.根据权利要求4所述的减振降噪效果的检验装置,其特征在于,所述力传递计算模块具体用于:
以两相邻轨枕之间的钢轨作为Timoshenko梁单元;
利用Timoshenko梁运动方程,计算Timoshenko梁单元的谱节点荷载和谱节点位移的关系;
依据Timoshenko梁单元的谱节点荷载和谱节点位移的关系,计算平面半轨道模型的谱节点荷载与谱节点位移的关系;
依据平面半轨道模型的谱节点荷载与谱节点位移的关系计算力传递率。
6.根据权利要求4所述的减振降噪效果的检验装置,其特征在于,所述钢轨声功率计算模块具体用于:
通过公式
|v(x)|=v(0)e-β|x|,
计算钢轨的声功率;
依据钢轨的声功率,计算声功率级:
其中,v(x)是钢轨在x处振动速度的幅值,v(0)为原点处的振动速度幅值,β为衰减系数,ρ0c0为空气中的声特性阻抗,ρ0=1.225kg/m3为空气密度,c0=340m/s为声波在空气中的传播速度,σ是与频率有关的辐射率,为钢轨的固有属性,P为一个截面的周长,W0为基准声功率,W0=10-12W,W为钢轨的声功率,Lw为声功率级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710535321.5A CN107292046B (zh) | 2017-07-03 | 2017-07-03 | 一种减振降噪效果的检验方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710535321.5A CN107292046B (zh) | 2017-07-03 | 2017-07-03 | 一种减振降噪效果的检验方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292046A CN107292046A (zh) | 2017-10-24 |
CN107292046B true CN107292046B (zh) | 2019-12-03 |
Family
ID=60099379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710535321.5A Active CN107292046B (zh) | 2017-07-03 | 2017-07-03 | 一种减振降噪效果的检验方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292046B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108318125B (zh) * | 2018-01-26 | 2021-03-16 | 华东交通大学 | 一种基于钢轨竖向振动特性检测道砟刚度的方法 |
AT522937B1 (de) * | 2019-09-04 | 2022-01-15 | Verfahren zur Erkennung technischer Irregularitäten von Eisenbahnfahrzeugen aus der Analyse von Schall- und Schwingungsdaten | |
CN110904744B (zh) * | 2019-11-22 | 2021-09-10 | 华东交通大学 | 一种改良钢轨的吸振***及改良吸振性能的方法 |
CN111979850B (zh) * | 2020-08-18 | 2022-04-26 | 华东交通大学 | 一种控制铁路过渡段振动的扣件布置方法和*** |
CN112580210A (zh) * | 2020-12-22 | 2021-03-30 | 中国铁路设计集团有限公司 | 一种一维周期垫层减振道床隔振频段调控设计方法 |
CN112560320A (zh) * | 2020-12-24 | 2021-03-26 | 华东交通大学 | 一种轨道结构弯曲振动带隙的计算方法 |
CN112948987A (zh) * | 2021-01-15 | 2021-06-11 | 西南交通大学 | 一种减振轨道声振性能和服役性能的综合评价方法 |
CN113111448A (zh) * | 2021-03-09 | 2021-07-13 | 广州地铁设计研究院股份有限公司 | 一种减振集能轨枕的设计方法 |
CN117194847B (zh) * | 2023-11-06 | 2024-02-09 | 国合通用(青岛)测试评价有限公司 | 一种减振钢轨的减振效果评估***及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0593608A1 (de) * | 1991-07-11 | 1994-04-27 | Jenbacher Energiesysteme Aktiengesellschaft | Lärmschutzvorrichtung für eisenbahngleise |
CN101697175A (zh) * | 2009-10-26 | 2010-04-21 | 华东交通大学 | 一种轨道交通噪声仿真预测方法 |
JP2010209540A (ja) * | 2009-03-09 | 2010-09-24 | Ihi Corp | 車両走行時の騒音低減方法及び装置 |
CN103456310A (zh) * | 2013-08-28 | 2013-12-18 | 大连理工大学 | 一种基于谱估计的瞬态噪声抑制方法 |
CN104036087A (zh) * | 2014-06-24 | 2014-09-10 | 同济大学 | 基于功率流-边界元模型的高架轨道交通振动噪声仿真预测方法 |
CN106446443A (zh) * | 2016-10-12 | 2017-02-22 | 西南交通大学 | 一种轨道扣件***共振频率的识别方法及装置 |
-
2017
- 2017-07-03 CN CN201710535321.5A patent/CN107292046B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0593608A1 (de) * | 1991-07-11 | 1994-04-27 | Jenbacher Energiesysteme Aktiengesellschaft | Lärmschutzvorrichtung für eisenbahngleise |
JP2010209540A (ja) * | 2009-03-09 | 2010-09-24 | Ihi Corp | 車両走行時の騒音低減方法及び装置 |
CN101697175A (zh) * | 2009-10-26 | 2010-04-21 | 华东交通大学 | 一种轨道交通噪声仿真预测方法 |
CN103456310A (zh) * | 2013-08-28 | 2013-12-18 | 大连理工大学 | 一种基于谱估计的瞬态噪声抑制方法 |
CN104036087A (zh) * | 2014-06-24 | 2014-09-10 | 同济大学 | 基于功率流-边界元模型的高架轨道交通振动噪声仿真预测方法 |
CN106446443A (zh) * | 2016-10-12 | 2017-02-22 | 西南交通大学 | 一种轨道扣件***共振频率的识别方法及装置 |
Non-Patent Citations (2)
Title |
---|
Research on the elastic wave band gaps of curved beam of phononic;Liu Shaogang 等;《Physica B》;20150131;第82-91页 * |
周期结构的振动行为与隔振性能研究;吴志静;《中国博士学位论文全文数据库 基础科学辑》;20170215;第A004-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107292046A (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107292046B (zh) | 一种减振降噪效果的检验方法及装置 | |
Thompson et al. | Track dynamic behaviour at high frequencies. Part 1: theoretical models and laboratory measurements | |
Ji et al. | Seismic damage detection of a full-scale shaking table test structure | |
Hong et al. | Reconstruction of dynamic displacement and velocity from measured accelerations using the variational statement of an inverse problem | |
Kouroussis et al. | On the interest of integrating vehicle dynamics for the ground propagation of vibrations: the case of urban railway traffic | |
Kouroussis et al. | Influence of some vehicle and track parameters on the environmental vibrations induced by railway traffic | |
Zhong et al. | Dynamic responses of prestressed bridge and vehicle through bridge–vehicle interaction analysis | |
Andersson et al. | Prediction of rail corrugation generated by three-dimensional wheel–rail interaction | |
Xu et al. | Extended applications of track irregularity probabilistic model and vehicle–slab track coupled model on dynamics of railway systems | |
He et al. | Experimental and three-dimensional finite element method studies on pounding responses of bridge structures subjected to spatially varying ground motions | |
Shen et al. | An inerter-based electromagnetic damper for civil structures: modeling, testing, and seismic performance | |
Vassilopoulou et al. | Nonlinear dynamic behavior of saddle-form cable nets under uniform harmonic load | |
JP2019189009A (ja) | 構造物疲労強度推定方法及びシステム | |
Shaban et al. | Vehicle effects on seismic response of a simple‐span bridge during shake tests | |
Koroma et al. | Influence of preload and nonlinearity of railpads on vibration of railway tracks under stationary and moving harmonic loads | |
Pu et al. | Model updating–based damage detection of a concrete beam utilizing experimental damped frequency response functions | |
Ebrahimzadeh Hassanabadi et al. | On the vibration of a thin rectangular plate carrying a moving oscillator | |
Kostovasilis et al. | A semi-analytical beam model for the vibration of railway tracks | |
Gras et al. | On a coupling between the Finite Element (FE) and the Wave Finite Element (WFE) method to study the effect of a local heterogeneity within a railway track | |
Li et al. | Influence of static preload on vibration reduction effect of floating slab tracks | |
Dou et al. | Effect of viscoelastic-plastic dynamic properties of rail pads on curved rail dynamic characteristics based on the modified SEM–SM hybrid method | |
Recuero et al. | Dynamics of the coupled railway vehicle–flexible track system with irregularities using a multibody approach with moving modes | |
Ronasi et al. | Inverse identification of wheel–rail contact forces based on observation of wheel disc strains: an evaluation of three numerical algorithms | |
Luo et al. | Dynamic simulation studies and experiments on rubber structures used in rail vehicles | |
Kaewunruen et al. | Application of experimental modal testing for estimating dynamic properties of structural components |
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 |