发明名称 一种路面及桥面铺装材料阻尼特性参数的测试方法
摘要 一种路面及桥面铺装材料阻尼特性参数的测试方法。其主要目的在于精确地测出路面或桥面铺装材料的阻尼比。其实现步骤为:根据实测数据计算确定阻尼测试时需考虑的车辆荷载激振的时间频率范围,取铺装材料长梁,并用有限元方法计算得到铺装材料长梁在车辆荷载激振时间频率范围内的所有模态振型,及其各自对应的无阻尼自振频率。对计算得到的各阶模态振型进行测试,获取各阶模态振型对应的激振信号和响应信号。根据所获得的激振信号和响应信号计算得到各阶模态振型的频响函数,并绘制各阶模态振型的频响函数幅值与频率的关系曲线,确定各阶模态振型的阻尼自振频率,从而计算出各阶模态振型对应的阻尼比。
申请公布号 CN102621225A 申请公布日期 2012.08.01
申请号 CN201210096866.8 申请日期 2012.03.30
申请人 东南大学 发明人 张磊;钱振东;王文炜;涂永明;张晓春;毕建新;何铁军;张宁;郭建华;夏井新;闵召辉;王华荣;朱浩然
分类号 G01N29/12(2006.01)I 主分类号 G01N29/12(2006.01)I
代理机构 南京天翼专利代理有限责任公司 32112 代理人 汤志武
主权项 1.一种路面及桥面铺装材料阻尼特性参数的测试方法,其特征在于,包括如下步骤:步骤1:试验前参数的确定步骤1.1:通过路面或桥面的不平度实测数据,即沿行车方向不同位置处行车道表面的高程,计算得到空间频率的范围,具体步骤如下:沿行车方向,对1km路面或桥面进行不平度测试,得到高程序列:{D(m)},m=0,1,2,…N-1,N为所测高程总点数,对该高程序列进行离散时间序列傅里叶变换,得到:<maths num="0001"><![CDATA[<math><mrow><msub><mi>a</mi><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>D</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>jk</mi><msub><mi>O</mi><mn>0</mn></msub><mi>m</mi></mrow></msup></mrow></math>]]></maths>其中,<img file="FDA0000148749190000012.GIF" wi="198" he="105" />N为所测高程总点数,a<sub>k</sub>为离散时间序列傅里叶变换系数,k=0,1,2,…N-1,对求得的离散时间序列傅里叶变换系数a<sub>k</sub>,k=0,1,2,…N-1,进行归一化,也就是依次将离散时间序列傅里叶变换系数a<sub>k</sub>除以a<sub>k</sub>序列中的最大值,得到a<sub>k</sub>,选择大于1%的所有a<sub>k</sub>,假设其中下标的最小值为k<sub>min</sub>,下标的最大值为k<sub>max</sub>,则空间频率n的范围为<img file="FDA0000148749190000013.GIF" wi="393" he="113" />单位为m<sup>-1</sup>,步骤1.2:确定通行车辆的行车速度范围根据道路或桥梁收费站或车辆监测点的车流数据,选择每日平均交通量中所占比例最大的车型作为代表性车型,通过车速监测点的监控数据,获得代表性车型的行车速度v的范围v<sub>min</sub>~v<sub>max</sub>,即为通行车辆的行车速度范围,步骤1.3:确定阻尼测试时需要考虑的车辆荷载激振的时间频率范围根据步骤1.1得到的空间频率n的范围<img file="FDA0000148749190000014.GIF" wi="350" he="113" />和步骤1.2得到的通行车辆行车速度v的范围v<sub>min</sub>~v<sub>max</sub>,确定车辆荷载激振的时间频率范围,所述的车辆荷载激振的时间频率范围可以根据如下公式计算确定:f=nv其中,n为空间频率,v为行车速度,f为时间频率,所述的车辆荷载激振的时间频率范围用圆频率表示为:Ω=2πf=2πnv其中,Ω为圆频率,步骤1.4:从待测路面或桥面上取一根长宽比大于5且小于20的铺装材料长梁作为试件,铺装材料长梁横截面的宽度和高度均大于铺装材料中集料规格尺寸的<img file="FDA0000148749190000015.GIF" wi="84" he="57" />倍且不大于100mm,步骤1.5:在步骤1.3中确定的阻尼测试时需要考虑的车辆荷载激振的时间频率范围内,采用有限元方法计算得到铺装材料长梁在车辆荷载激振的时间频率范围内的N′阶模态振型,所述N′阶模态振型包括第1阶模态振型,第2阶模态振型,第3阶模态振型,……,第N′阶模态振型,N′为模态振型阶数,第1阶模态振型,第2阶模态振型,第3阶模态振型,……,第N′阶模态振型所对应的模态振型向量依次分别为:{Ψ′}<sub>1</sub>,{Ψ′}<sub>2</sub>,{Ψ′}<sub>3</sub>,……,{Ψ′}<sub>N′</sub>,并同时计算出每一阶模态振型对应的无阻尼自振频率Ω<sub>1</sub>,Ω<sub>2</sub>,Ω<sub>3</sub>,……,Ω<sub>N′</sub>,有限元方法计算得到铺装材料长梁的N′阶模态振型的具体步骤如下:步骤1.5.1:定义单元类型为梁单元,并进行梁单元实常数的设置,具体包括铺装材料长梁横截面的面积、惯性矩、高度,其中铺装材料长梁横截面的高度为h,面积为A=bh,惯性矩为<img file="FDA0000148749190000021.GIF" wi="210" he="113" />b为铺装材料长梁的宽度,铺装材料长梁的长度为L′,并定义铺装材料长梁的材料属性,包括弹性模量和密度,其中,按《公路工程沥青及沥青混合料试验规程》(JTG E20-2011)中T0715-2011沥青混合料弯曲试验,测定其弹性模量,并按《公路工程沥青及沥青混合料试验规程》(JTG E20-2011)中T0706-2011压实沥青混合料密度试验,即水中重法,测定其密度,根据铺装材料长梁的尺寸进行建模并划分网格,网格划分数为M′,具体如下:令x′为<img file="FDA0000148749190000022.GIF" wi="53" he="107" />取整,网格划分数M′在如下范围取值[x′,4x′],步骤1.5.2:将铺装材料长梁两端的边界条件设为简支,并定义分析类型为模态分析,进行分析参数的设置,步骤1.5.3:进行求解和结果分析,步骤2:阻尼测试步骤2.1:采用压电材料传感器,作为激振力发射端,又兼作为响应信号接收端,分别同时连接到信号发生器和存储器,步骤2.2:令r=1,r为所测试的模态振型的阶次,步骤2.3:根据第r阶模态振型,在铺装材料长梁上布设r个压电材料传感器,其中,r≤N′,所述压电材料传感器布设的具体位置为:<maths num="0002"><![CDATA[<math><mfrac><mrow><mrow><mo>(</mo><mn>2</mn><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>L</mi></mrow><mrow><mn>2</mn><mi>r</mi></mrow></mfrac></math>]]></maths>其中,i为压电材料传感器的编号,i=1,2,…,r,r为所测试的模态振型的阶次,L为铺装材料长梁的跨径,令其等于铺装材料长梁的长度L′,同时分别设置第i个压电材料传感器的初始相位<img file="FDA0000148749190000024.GIF" wi="63" he="38" />i=1,2,…,r,当i为奇数时初始相位<img file="FDA0000148749190000025.GIF" wi="34" he="38" />设为0,当i为偶数时初始相位<img file="FDA0000148749190000026.GIF" wi="34" he="38" />设为π,步骤2.4:确定第r阶模态振型的频率扫描范围,具体过程如下:由步骤1.5中获得的第r阶模态振型的无阻尼自振频率Ω<sub>r</sub>,根据所测试的铺装材料长梁的材料种类选取其可能的阻尼比ξ<sup>ke</sup>,对一般道路工程而言,结构阻尼比通常取0~14%,由下式:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&omega;</mi><mi>r</mi><mi>ke</mi></msubsup><mo>=</mo><msqrt><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><msup><mi>&xi;</mi><mi>ke</mi></msup><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msub><mi>&Omega;</mi><mi>r</mi></msub></mrow></math>]]></maths>计算得到第r阶模态振型可能的阻尼自振频率<img file="FDA0000148749190000028.GIF" wi="88" he="57" />其中<img file="FDA0000148749190000029.GIF" wi="67" he="57" />为第r阶模态振型可能的阻尼自振频率,并确定频率扫描范围,所述频率扫描范围为:<img file="FDA00001487491900000210.GIF" wi="341" he="65" />其中,<img file="FDA00001487491900000211.GIF" wi="107" he="64" />为第r阶模态振型频率扫描范围的下限,<img file="FDA00001487491900000212.GIF" wi="107" he="64" />取<img file="FDA00001487491900000213.GIF" wi="455" he="57" /><img file="FDA00001487491900000214.GIF" wi="107" he="65" />为第r阶模态振型频率扫描范围的上限,<img file="FDA0000148749190000031.GIF" wi="107" he="64" />取<img file="FDA0000148749190000032.GIF" wi="473" he="57" />步骤2.5:选取第r阶模态振型,在步骤2.4确定的第r阶模态振型的频率扫描范围<img file="FDA0000148749190000033.GIF" wi="341" he="64" />内,以<img file="FDA0000148749190000034.GIF" wi="96" he="108" />为步长进行频率扫描,令<img file="FDA0000148749190000035.GIF" wi="249" he="125" /><img file="FDA0000148749190000036.GIF" wi="235" he="127" /><img file="FDA0000148749190000037.GIF" wi="244" he="120" />T为扫描时间,一般在2s~10s范围内取值,由公式<img file="FDA0000148749190000038.GIF" wi="598" he="108" />可得第r阶模态振型第i个压电材料传感器在在t<sub>q</sub>时刻的激振力序列{y<sub>qri</sub>},其中y<sub>qri</sub>为第r阶模态振型下第i个压电材料传感器在t<sub>q</sub>时刻的激振力,t<sub>q</sub>=qΔt,<img file="FDA0000148749190000039.GIF" wi="160" he="108" />S为第r阶模态振型的频率扫描的激振力激励的时间点数,q=1,2,3,…,S,<img file="FDA00001487491900000310.GIF" wi="35" he="38" />为第r阶模态振型下第i个压电材料传感器的初始相位,由信号发生器分别产生对应于第i个压电材料传感器的激振信号{y<sub>qri</sub>},i=1,2,…,r,并同时将激振信号发送给各对应的压电材料传感器,由各个压电材料传感器产生激振力,激起铺装材料长梁的第r阶模态振型,同时各压电材料传感器又兼作为响应信号接收端,接收铺装材料长梁的振动响应信号,并将接收到的铺装材料长梁的振动响应信号发送至存储器,步骤2.6:如果r+1>N′时,则停止振型测试,进入步骤3,否则,令r=r+1,返回步骤2.3,继续进行振型测试,步骤3:阻尼比确定方法步骤3.1:令r=1,r为所测试的模态振型阶次,步骤3.2:根据步骤2.5所获得的第r阶模态振型的激振信号和响应信号获取第r阶模态振型的频响函数,当忽略噪声的影响时,进入步骤3.2.1,当考虑噪声的影响时,进入步骤3.2.2,步骤3.2.1:在存储器中选取r个压电材料传感器中的任意一个压电材料传感器的激振信号和响应信号,并对选取的激振信号和响应信号进行傅里叶变换,根据激振信号的傅里叶变换和响应信号的傅里叶变换,根据以下公式求得在忽略噪声的影响时第r阶模态振型的频响函数H<sub>r</sub>(ω),<maths num="0004"><![CDATA[<math><mrow><msub><mi>H</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mi>ri</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>F</mi><mi>ri</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中,ω为振动圆频率,也可采用时间频率<img file="FDA00001487491900000312.GIF" wi="165" he="97" />作为自变量,即用ω=2πf′代入上述公式,i为压电材料传感器的编号,i=1,2,…,r,X<sub>ri</sub>(ω)为第r阶模态振型下第i个压电材料传感器采集到的响应信号的傅里叶变换,F<sub>ri</sub>(ω)为第r阶模态振型下第i个压电材料传感器激振信号的傅里叶变换,步骤3.2.2:在存储器中选取r个压电材料传感器中的任意一个压电材料传感器的激振信号和响应信号,根据选取的激振信号和响应信号,求得第r阶模态振型的激振信号与响应信号之间的互功率谱、激振信号的自功率谱和响应信号的自功率谱,根据以下公式求得在考虑噪声的影响时第r阶模态振型的频响函数H′<sub>r</sub>(ω),<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>H</mi><mi>r</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>X</mi><mi>ri</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msubsup><mi>X</mi><mi>ri</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>F</mi><mi>ri</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msubsup><mi>X</mi><mi>ri</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><msub><mi>G</mi><mi>XXri</mi></msub><msub><mi>G</mi><mi>FXri</mi></msub></mfrac></mrow></math>]]></maths>其中,ω为振动圆频率,也可采用时间频率<img file="FDA0000148749190000041.GIF" wi="165" he="97" />作为自变量,即用ω=2πf′代入上述公式,i为压电材料传感器的编号,i=1,2,…,r,X<sub>ri</sub>(ω)为第r阶模态振型下第i个压电材料传感器采集到的响应信号的傅里叶变换,F<sub>ri</sub>(ω)为第r阶模态振型下第i个压电材料传感器激振信号的傅里叶变换,*为共轭符号,G<sub>FXri</sub>为第r阶模态振型下第i个压电材料传感器的激振信号与响应信号之间的互功率谱,G<sub>XXri</sub>为第r阶模态振型下第i个压电材料传感器的响应信号的自功率谱,步骤3.3:根据步骤3.2求得的第r阶模态振型的频响函数,绘制第r阶模态振型的频响函数幅值与频率的关系曲线,关系曲线中幅值最大值对应的频率即为第r阶模态振型的阻尼自振频率ω<sub>r</sub>,其中r为所测试的模态振型的阶次,步骤3.4:根据步骤1.5确定的第r阶模态振型的无阻尼自振频率Ω<sub>r</sub>,以及步骤3.3求得的第r阶模态振型的阻尼自振频率ω<sub>r</sub>,根据公式:<maths num="0006"><![CDATA[<math><mrow><msub><mi>&omega;</mi><mi>r</mi></msub><mo>=</mo><msqrt><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msub><mi>&Omega;</mi><mi>r</mi></msub></mrow></math>]]></maths>就可以求得第r阶模态振型对应的阻尼比ξ<sub>r</sub>,步骤3.5:如果r+1>N′时,则停止阻尼比计算,否则,令r=r+1,返回步骤3.2,继续进行阻尼比计算。
地址 210096 江苏省南京市四牌楼2号