CN107038339A - 一种城市轨道交通环境振动的计算方法 - Google Patents

一种城市轨道交通环境振动的计算方法 Download PDF

Info

Publication number
CN107038339A
CN107038339A CN201710206784.7A CN201710206784A CN107038339A CN 107038339 A CN107038339 A CN 107038339A CN 201710206784 A CN201710206784 A CN 201710206784A CN 107038339 A CN107038339 A CN 107038339A
Authority
CN
China
Prior art keywords
vibration
acceleration
track traffic
frequency
urban track
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.)
Granted
Application number
CN201710206784.7A
Other languages
English (en)
Other versions
CN107038339B (zh
Inventor
施卫星
谢志行
拜立岗
章仪青
刘念
王梁坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN201710206784.7A priority Critical patent/CN107038339B/zh
Publication of CN107038339A publication Critical patent/CN107038339A/zh
Application granted granted Critical
Publication of CN107038339B publication Critical patent/CN107038339B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Abstract

本发明涉及一种城市轨道交通环境振动的计算方法,用于对应的城市轨道交通环境振动进行统一计算,所述方法包括下列步骤:测量轨道交通环境的实际振动加速度时程a(j),通过傅里叶变换,计算得到1/3倍频程不同中心频率处的加速度有效值ari;根据加速度有效值ari,进行振级计算,得到城市轨道交通环境的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值;将分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值进行可视化展示。与现有技术相比,本发明具有简单易行、显示直观以及对比明确等优点。

Description

一种城市轨道交通环境振动的计算方法
技术领域
本发明涉及轨道交通领域,尤其是涉及一种城市轨道交通环境振动的计算方法。
背景技术
随着中国城市经济迅猛发展,人口不断增多,交通问题在城市规划中越来越受到重视,城市轨道交通以其安全性、准时性、运送量大等优点逐渐成为解决城市交通问题的重要方法。但是,随着城市轨道交通路网逐渐完善,其引发的环境振动问题也越来越受到人们的重视。在国际上,环境振动已经被列为7大环境公害之一,并且已经出台了一系列规范对环境振动进行限定。国际上现行相关规范主要包括:国际标准ISO 2631-1:1997、ISO 2631-2:2003,日本的JIS Z 8735-1981,德国的DIN 4150-1、DIN 4150-2和DIN 4150-3以及英国的BS 6472-1:2008等。我国结合上述规范以及实际国情也颁布了相关规范。
国家环境保护局最早于1989年批准实施了GB 10070-88《城市区域环境振动标准》(以下简称GB 10070-88)以及GB 10071-88《城市区域环境振动测量方法》(以下简称GB10071-88),该规范规定了城市区域环境振动的标准值及其适用地带范围和监测方法[8,9]。此后,随着行业的发展以及国家对于环境振动的日渐重视,相关的标准呈现出细分化的趋势。而针对于交通振动对于建筑物的影响,主要使用GB 10070-88、JBJ 16-2000J61-2000、JGJ/T 170-2009、GB 50868-2013规范进行评价,由于规范繁多,计算流程复杂,给工程人员对于轨道交通的振动计算造成了很大的困扰。
发明内容
本发明的目的是针对上述问题提供一种城市轨道交通环境振动的计算方法。
本发明的目的可以通过以下技术方案来实现:
一种城市轨道交通环境振动的计算方法,用于对应的城市轨道交通环境振动进行统一计算,所述方法包括下列步骤:
1)测量轨道交通环境的实际振动加速度时程a(j),通过傅里叶变换,计算得到1/3倍频程不同中心频率处的加速度有效值ari
2)根据步骤1)得到的加速度有效值ari,进行振级计算,得到城市轨道交通环境的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值;
3)将步骤2)中得到的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值进行可视化展示。
所述步骤1)具体为:
11)测量轨道交通环境的实际振动加速度时程a(j),其中j=1,2,3,...,N;
12)通过快速傅里叶变换,将步骤11)测量的实际振动加速度时程a(j)变换至频域内,计算频域内各频率点处的幅值|Ai(j)|;
13)根据步骤12)中得到的幅值|Ai(j)|,计算得到1/3倍频程不同中心频率处的加速度有效值ari
所述步骤13)具体为:
131)判断步骤1)测量的实际振动加速度时程a(j)是否为离散值,若是则进入步骤132),若否则进入步骤133);
132)直接根据频域内各频率点处的幅值|Ai(j)|进行计算,得到1/3倍频程不同中心频率处的加速度有效值ari
133)对步骤12)中得到的幅值进行傅里叶逆变换,得到各频带的振动加速度时程ai(j),根据ai(j)的值计算得到1/3倍频程不同中心频率处的加速度有效值ari
所述步骤132)具体为:
其中,Aimax为Ai(j)在1/3倍频程所对应频段的幅值谱的最大幅值,f为相应的频率。
所述根据ai(j)的值计算得到1/3倍频程不同中心频率处的加速度有效值ari具体为:
所述步骤2)具体为:
21)根据步骤1)得到的加速度有效值ari,计算得到振动计权加速度值awi
22)根据步骤21)计算得到的振动计权加速度值awi,计算得到分频加速度振级VAL或分频计权加速度振级VAL';
23)根据步骤22)得到的分频加速度振级VAL,计算得到Z振级VLz
24)根据步骤21)计算得到的振动计权加速度值awi,计算得到竖向四次方振动剂量值VDVz
所述振动计权加速度值awi具体为:
awi=ari·Wik
其中,Wik为1/3倍频程第i段对应的加速度计权因子。
所述分频加速度振级VAL具体为:
所述分频计权加速度振级VAL'具体为:
其中,a0为基准加速度值。
所述Z振级VLz具体为:
其中,Wi为1/3倍频程第i段对应的分贝计权因子。
所述竖向四次方振动剂量值VDVz具体为:
与现有技术相比,本发明具有以下有益效果:
(1)本发明通过将国家各个规范下进行城市轨道交通环境振动评判所需要的参数进行系列化的统一计算,只需测量轨道交通环境的实际振动加速度时程,即可得到与各个规范对应的振级参数,计算程序规范化且计算方式简单,与现有技术每次都要进行独立计算相比,便于工程人员一次化对城市轨道交通环境的振动进行评判,简单方便,实用性能强。
(2)本发明提出的计算方法,是按照计算的最简顺序来实现的,比如说Z振级的求取需要分频加速度振级之后才可进行,因而这种计算方法是最为简便且不会造成重复计算的方法,大大节省了工作人员的时间和计算量。
(3)通过傅里叶变换,将时域上的计算转变成频域进行对比计算,因而可以在频域范围内对振动的成分进行查看,对比明确,更加便于对振动程度的判断。
(4)本发明最终将得到的结果进行可视化展示,结果显示直观,便于工程人员进行观察和监测。
附图说明
图1为本发明的方法流程图;
图2为实施例中六层居民楼的方位图;
图3为加速度时程图;
图4为频谱图;
图5为计算结果与限值对比图;
图6为具体的计算过程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明提出的一种城市轨道交通环境振动的计算方法,用于统一对与国家各个规范相对应的城市轨道交通环境振动进行计算,该方法包括下列步骤:
1)测量轨道交通环境的实际振动加速度时程a(j),通过傅里叶变换,计算得到1/3倍频程不同中心频率处的加速度有效值ari
11)测量轨道交通环境的实际振动加速度时程a(j),其中j=1,2,3,...,N;
12)通过快速傅里叶变换,将步骤11)测量的实际振动加速度时程a(j)变换至频域内,计算频域内各频率点处的幅值|Ai(j)|;
13)根据步骤12)中得到的幅值|Ai(j)|,计算得到1/3倍频程不同中心频率处的加速度有效值ari
131)判断步骤1)测量的实际振动加速度时程a(j)是否为离散值,若是则进入步骤132),若否则进入步骤133);
132)直接根据频域内各频率点处的幅值|Ai(j)|进行计算,得到1/3倍频程不同中心频率处的加速度有效值ari
其中,Aimax为Ai(j)在1/3倍频程所对应频段的幅值谱的最大幅值,f为相应的频率;
133)对步骤12)中得到的幅值进行傅里叶逆变换,得到各频带的振动加速度时程ai(j),根据ai(j)的值计算得到1/3倍频程不同中心频率处的加速度有效值ari
2)根据步骤1)得到的加速度有效值ari,进行振级计算,得到与国家各个规范相对应的城市轨道交通环境的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值:
21)根据步骤1)得到的加速度有效值ari,计算得到振动计权加速度值awi
awi=ari·Wik
其中,Wik为1/3倍频程第i段对应的加速度计权因子;
22)根据步骤21)计算得到的振动计权加速度值awi,计算得到分频加速度振级VAL或分频计权加速度振级VAL':
分频计权加速度振级VAL'具体为:
其中,a0为基准加速度值;
23)根据步骤22)得到的分频加速度振级VAL,计算得到Z振级VLz
其中,Wi为1/3倍频程第i段对应的分贝计权因子;
24)根据步骤21)计算得到的振动计权加速度值awi,计算得到竖向四次方振动剂量值VDVz
3)将步骤2)中得到的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值与国家标准的规范值进行对比,并进行可视化展示。
根据上述步骤进行具体的计算,如图6所示,过程如下:
参照图2所示实际工程,经过现场测量得到如图3所示的工程加速度时程图,经过一次Fourier变换得到如图4所示的频谱图,具体的计算过程如下:
(1)时程加速度有效值计算:在城市轨道交通振动评价指标计算中,首先需要得到振动加速度的有效值a,单位m/s2为,其计算公式如(1.1)所示。对于解析信号而言,采用(1.1)式的中间表达式,而对于离散采样数据,则采用后一个表达式(1.2);其中T为振动数据的计权时间,N为采样数据的点数。
(2)计权加速度有效值计算:结合1/3倍频程不同中心频率处的加速度有效值,按GB/T 13441.1-2007或者ISO 2631-1:1985中的计权因子Wk、Wd以及W、Wh修正后得到计权加速度有效值awi,如式(1.3)所示。
(3)振级计算:如式(1.4)所示,单位为dB。其中:a0为基准加速度值,规范取值为10-6m/s2;在GB/T 50355-2005中的限值是分频铅垂向振动加速度级,由于采用的铅垂向计权因子与GB10070-88中相同,所以可以推导出各1/3倍频程中心频率处的“Z向分频振动加速度级”与“Z振级”间的转换关系如式(1.5)所示。式中La,i为1/3倍频程中心频率处的“振动加速度级”;为GB/T 50355-2005条文说明1.0.5条给出的铅垂向计权因子(即ISO 2631-1:1985之中的铅垂向计权因子W),而(La,i-Wi)即为“分频振级”。在计算水平向振级VLxy时,只需将铅垂向计权因子换成水平向计权因子即可。同时由于交通振动的间歇性和长期性,以计权加速度均方根值为基础的基本评价方法可能会低估振动对人体舒适性的影响,故在GB50868-2013中附加了竖向法,即竖向四次方振动剂量值,计算公式如式(1.6)所示。
上述计算公式如下:
式(1.1),式中ai(j)为经过快速Fourier变换和Fourier逆变换得到的各频带的振动加速度时程,ari为1/3倍频程第i频段的加速度有效值。
式(1.2),式中Aimax为1/3倍频程第i频段的幅值谱的最大幅值,f为相应的频率。
式(1.3),awi=ariWik,式中Wik为1/3倍频程第i频段对应的加速度计权因子,awi为振动计权加速度值。
式(1.4),式中a0为基准加速度值,一般取用a0=10-6m/s2,当a为1/3倍频程第i频段的加速度有效值ari时,得到的VAL为分频加速度振级La,i,当a为振动计权加速度值awi时,得到的VAL为分频计权加速度振级。
式(1.5),式中La,i为分频加速度振级,Wi为1/3倍频程第i频段对应的分贝计权因子,VLz为Z振级。
式(1.6),式中awi为振动计权加速度值,VDVz为竖向四次方振动剂量值。
运用编制的MATLAB程序进行上述计算过程,可以直接得到本工程的与规范限值对比的振级,结果如图5所示,经计算本工程的可知分频最大振级VLmax和铅垂向振级VLz均超过了相关规范之中的限值,会对居民生活造成影响,需要进行减振降噪工作,并且由图3频谱图可以知道,40Hz的振动是本次轨道交通振动的主要振动成分,在进行工程减振时需要主要控制在40Hz左右的振动。通过MATLAB编程,利用傅里叶变换,将时域之中的加速度转变成频域进行对比计算,只需要输入加速度既可以得到最终振级,同时还可以在频域范围内对振动的成分进行查看,具有简单易行、对比明确的优点。

Claims (10)

1.一种城市轨道交通环境振动的计算方法,用于对应的城市轨道交通环境振动进行统一计算,其特征在于,所述方法包括下列步骤:
1)测量轨道交通环境的实际振动加速度时程a(j),通过傅里叶变换,计算得到1/3倍频程不同中心频率处的加速度有效值ari
2)根据步骤1)得到的加速度有效值ari,进行振级计算,得到城市轨道交通环境的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值;
3)将步骤2)中得到的分频加速度振级、Z振级、分频计权加速度振级和竖向四次方振动剂量值进行可视化展示。
2.根据权利要求1所述的城市轨道交通环境振动的计算方法,其特征在于,所述步骤1)具体为:
11)测量轨道交通环境的实际振动加速度时程a(j),其中j=1,2,3,...,N;
12)通过快速傅里叶变换,将步骤11)测量的实际振动加速度时程a(j)变换至频域内,计算频域内各频率点处的幅值|Ai(j)|;
13)根据步骤12)中得到的幅值|Ai(j)|,计算得到1/3倍频程不同中心频率处的加速度有效值ari
3.根据权利要求2所述的城市轨道交通环境振动的计算方法,其特征在于,所述步骤13)具体为:
131)判断步骤1)测量的实际振动加速度时程a(j)是否为离散值,若是则进入步骤132),若否则进入步骤133);
132)直接根据频域内各频率点处的幅值|Ai(j)|进行计算,得到1/3倍频程不同中心频率处的加速度有效值ari
133)对步骤12)中得到的幅值进行傅里叶逆变换,得到各频带的振动加速度时程ai(j),根据ai(j)的值计算得到1/3倍频程不同中心频率处的加速度有效值ari
4.根据权利要求2所述的城市轨道交通环境振动的计算方法,其特征在于,所述步骤132)具体为:
其中,Aimax为Ai(j)在1/3倍频程所对应频段的幅值谱的最大幅值,f为相应的频率。
5.根据权利要求2所述的城市轨道交通环境振动的计算方法,其特征在于,所述根据ai(j)的值计算得到1/3倍频程不同中心频率处的加速度有效值ari具体为:
6.根据权利要求1所述的城市轨道交通环境振动的计算方法,其特征在于,所述步骤2)具体为:
21)根据步骤1)得到的加速度有效值ari,计算得到振动计权加速度值awi
22)根据步骤21)计算得到的振动计权加速度值awi,计算得到分频加速度振级VAL或分频计权加速度振级VAL';
23)根据步骤22)得到的分频加速度振级VAL,计算得到Z振级VLz
24)根据步骤21)计算得到的振动计权加速度值awi,计算得到竖向四次方振动剂量值VDVz
7.根据权利要求6所述的城市轨道交通环境振动的计算方法,其特征在于,所述振动计权加速度值awi具体为:
awi=ari·Wik
其中,Wik为1/3倍频程第i段对应的加速度计权因子。
8.根据权利要求6所述的城市轨道交通环境振动的计算方法,其特征在于,所述分频加速度振级VAL具体为:
所述分频计权加速度振级VAL'具体为:
其中,a0为基准加速度值。
9.根据权利要求6所述的城市轨道交通环境振动的计算方法,其特征在于,所述Z振级VLz具体为:
其中,Wi为1/3倍频程第i段对应的分贝计权因子。
10.根据权利要求6所述的城市轨道交通环境振动的计算方法,其特征在于,所述竖向四次方振动剂量值VDVz具体为:
2
CN201710206784.7A 2017-03-31 2017-03-31 一种城市轨道交通环境振动的计算方法 Expired - Fee Related CN107038339B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710206784.7A CN107038339B (zh) 2017-03-31 2017-03-31 一种城市轨道交通环境振动的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710206784.7A CN107038339B (zh) 2017-03-31 2017-03-31 一种城市轨道交通环境振动的计算方法

Publications (2)

Publication Number Publication Date
CN107038339A true CN107038339A (zh) 2017-08-11
CN107038339B CN107038339B (zh) 2020-11-27

Family

ID=59534291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710206784.7A Expired - Fee Related CN107038339B (zh) 2017-03-31 2017-03-31 一种城市轨道交通环境振动的计算方法

Country Status (1)

Country Link
CN (1) CN107038339B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109100101A (zh) * 2018-07-05 2018-12-28 北京市劳动保护科学研究所 一种模拟随机振动源强的方法
CN112698383A (zh) * 2020-12-29 2021-04-23 浙江工业大学 一种城市轨道交通引起的环境振动预测方法及系统
CN113386822A (zh) * 2021-06-01 2021-09-14 唐智科技湖南发展有限公司 一种车轮多边失圆综合诊断方法、装置及设备
CN113962127A (zh) * 2021-10-22 2022-01-21 北京市科学技术研究院城市安全与环境科学研究所 一种竖向波屏障减振效率快速估算方法
CN116109786A (zh) * 2023-02-16 2023-05-12 中铁四院集团南宁勘察设计院有限公司 一种城市轨道交通振动地图的构建方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5366081B2 (ja) * 2009-03-05 2013-12-11 地方独立行政法人大阪府立産業技術総合研究所 振動生成方法および振動生成装置
CN103616192A (zh) * 2013-12-06 2014-03-05 中联重科股份有限公司 一种挖掘机的振动舒适度评价方法及系统
CN104180894A (zh) * 2014-08-28 2014-12-03 华东交通大学 铁路环境振动实时监测分析系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5366081B2 (ja) * 2009-03-05 2013-12-11 地方独立行政法人大阪府立産業技術総合研究所 振動生成方法および振動生成装置
CN103616192A (zh) * 2013-12-06 2014-03-05 中联重科股份有限公司 一种挖掘机的振动舒适度评价方法及系统
CN104180894A (zh) * 2014-08-28 2014-12-03 华东交通大学 铁路环境振动实时监测分析系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
凌育洪 等: "地铁引起的振动对框架结构的影响及隔振研究_以某教学楼为例", 《振动与冲击》 *
曾泽民: "地铁车辆段列车运行引发振动与噪声效应的现场试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
王兆亮: "城轨引起近地表振动的传播规律研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109100101A (zh) * 2018-07-05 2018-12-28 北京市劳动保护科学研究所 一种模拟随机振动源强的方法
CN112698383A (zh) * 2020-12-29 2021-04-23 浙江工业大学 一种城市轨道交通引起的环境振动预测方法及系统
CN113386822A (zh) * 2021-06-01 2021-09-14 唐智科技湖南发展有限公司 一种车轮多边失圆综合诊断方法、装置及设备
CN113962127A (zh) * 2021-10-22 2022-01-21 北京市科学技术研究院城市安全与环境科学研究所 一种竖向波屏障减振效率快速估算方法
CN116109786A (zh) * 2023-02-16 2023-05-12 中铁四院集团南宁勘察设计院有限公司 一种城市轨道交通振动地图的构建方法

Also Published As

Publication number Publication date
CN107038339B (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
CN107038339A (zh) 一种城市轨道交通环境振动的计算方法
US20230314272A1 (en) Method for real-time identification, monitoring, and early warning of vortex-induced vibration event of long-span suspension bridge
Kwon et al. Gust-front factor: New framework for wind load effects on structures
Cao et al. Measurement and analysis of vibrations in a residential building constructed on an elevated metro depot
He et al. EMD-based random decrement technique for modal parameter identification of an existing railway bridge
CN106126904A (zh) 多层建筑物的吸收爆破振动能量舒适性评价方法
Zhang et al. Elevator ride comfort monitoring and evaluation using smartphones
Boggs et al. Acceleration indexes for human comfort in tall buildings–Peak or RMS
WO2023087890A1 (zh) 一种基于动应力、振动和oma综合分析判断构架模态共振方法
Dong et al. Vibration source identification of offshore wind turbine structure based on optimized spectral kurtosis and ensemble empirical mode decomposition
Lin et al. Modal identification from nonstationary ambient response data using extended random decrement algorithm
CN104458173B (zh) 钢框架结构突变损伤的识别方法及系统
CN109459130A (zh) 一种基于智能手机的建筑结构舒适度测试系统
Zhu et al. Ride comfort evaluation of stochastic traffic flow crossing long-span suspension bridge experiencing vortex-induced vibration
CN108133098A (zh) 基于fe-sea混合法的桥梁局部振动及结构噪声的预测方法
Dai et al. Study on the damping ratios of reinforced concrete structures from seismic response records
Shakib et al. Evaluation of accidental eccentricity in buildings due to rotational component of earthquake
Yu et al. Wind-induced response of an L-shaped cable support glass curtain wall
Li et al. Dynamic system identification of a high-rise building during Typhoon Kalmaegi
Wang et al. Discrete analytical mode decomposition with automatic bisecting frequency selection for structural dynamic response analysis and modal identification
Su et al. Impulsive resistant optimization design of tuned viscous mass damper (TVMD) based on stability maximization
CN106446386B (zh) 模态能量法中模态间耦合强度的一种界定方法
He et al. Acceleration‐based effective damping ratio of tuned mass damper under seismic excitations
Ge et al. Structural dynamic responses of building structures with non-viscous dampers under Kanai–Tajimi spectrum excitation
CN102270249B (zh) 识别零部件特征频率的方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201127

CF01 Termination of patent right due to non-payment of annual fee