CN104573274B - 车辆荷载下基于位移时程面积的结构有限元模型修正方法 - Google Patents
车辆荷载下基于位移时程面积的结构有限元模型修正方法 Download PDFInfo
- Publication number
- CN104573274B CN104573274B CN201510042275.6A CN201510042275A CN104573274B CN 104573274 B CN104573274 B CN 104573274B CN 201510042275 A CN201510042275 A CN 201510042275A CN 104573274 B CN104573274 B CN 104573274B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mover
- msubsup
- 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.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种车辆荷载下基于位移时程面积的结构有限元模型修正方法,基于位移响应时程面积的目标函数和结构的局部刚度一一对应,因此可以选取局部刚度和边界条件作为修正参数,目标函数和修正变量意义明确,有助于提高计算效率并易于收敛,该方法不需要桥梁交通中断,不需要提取结构的动力特性(即不需要进行时频转换),因此可避免现有方法的不足。
Description
技术领域
本发明涉及一种有限元模型修正方法,具体是指一种车辆荷载下基于位移时程面积的结构有限元模型修正方法。
背景技术
随着新材料和施工技术的飞速发展,我国在桥梁建设方面取得的成就举世瞩目,有限元软件的快速发展使大跨复杂桥梁结构的精细化模拟、特殊荷载下(地震荷载、冲击荷载、撞击荷载、爆炸荷载等)的结构响应预测成为可能,一个比较精确的有限元模型有如下好处:(1)可以对结构进行定量化评估;(2)对结构精确分析(疲劳分析等)提供依据;(3)预测结构的响应;(4)对结构控制提供依据。
然而,在进行有限元建模时,由于材料、边界条件等参数的不确定性,导致建立的有限元模型和实际结构往往有较大差距。在此背景下,有限元模型修正技术应运而生,最常用的是基于灵敏度分析的有限元模型修正方法,其方法的原理是通过迭代优化算法来使实际测试数据和有限元模型计算的数据组合而成的目标函数最小化,进而进行参数识别(结构刚度和边界条件),使得修正后的有限元模型和实际结构尽可能一致。
目前,桥梁有限元模型修正方法大体分为两种,第一,基于静态响应(位移、应变等)的有限元模型修正方法;第二,基于振动特性的有限元模型修正方法。基于静态响应的有限元模型修正方法原理是:在桥梁结构上施加静态车辆荷载,利用结构的应变或挠度构建目标函数。基于振动特性的有限元模型修正方法原理是分析振动数据进而提取出桥梁结构的动力特性构建目标函数(基于自振频率的目标函数、基于位移模态的目标函数、基于模态柔度的目标函数、基于模态应变能的目标函数),这些方法存在的缺点如下:
(1)基于静态响应的有限元模型修正方法需要中断桥梁交通,在实际操作中很难做到这一点。
(2)基于振动特性的有限元模型修正方法需要利用傅里叶变换把时域振动数据转换为频域数据,会产生时频转换误差,进而影响模型修正的精度;
(3)由于桥梁的实际激励大小及形式是未知的,在提取桥梁动力特性的过程中,通常假定环境激励(车辆、风、地脉动)满足白噪声的理想分布,但是由于实际激励的复杂性,往往这个假定不能成立,导致动力特性参数的识别值和真实值存在误差,进而影响模型修正的精度。
发明内容
本发明的目的是克服现有技术的不足,提供一种车辆荷载下基于位移时程面积的结构有限元模型修正方法。
本发明采用的技术方案为:一种车辆荷载下基于位移时程面积的结构有限元模型修正方法,包括以下步骤:
步骤1:在目标结构的关键区域准分布式布置位移传感器,测试已知移动荷载作用下各个测试点的位移响应时程dj(t);
对于桥梁结构,跨度为l,梁高度为H,截面j-1,j,j+1,j+2沿桥梁长度方向的坐标分别为xj-1,xj,xj+1,xj+2,假设截面j-1,j,j+1,j+2等间隔,间距为L,移动荷载经过桥梁的整个过程中,截面j-1,j,j+1,j+2处的竖向位移响应为分别为dj-1(x),dj(x),dj+1(x),dj+2(x)(各截面随着移动荷载位移x变化所对应的位移),与之对应的截面j-1,j,j+1,j+2处的位移响应时程为dj-1(t),dj(t),dj+1(t),dj+2(t)(各截面随着时间t变化所对应的位移),其中移动荷载的参数如下:共有n个轴,轴重分别为P1,P2…Pi,Pn,速度为v;
假设结构符合欧拉梁假定,则j截面和j+1截面之间单元的底部平均应变εj,j+1(t)表达为
移动荷载作用下,截面j和j+1之间单元底部的平均应变表达为
式中di(i=1~n)为移动荷载的第i个轴距第1个轴之间的距离,其中d1=0,fj,j+1(x)为截面j和j+1之间单元底部的平均应变影响线,x为第1个轴距左边支座的距离;
把公式(2)左右部分分别沿着结构长度方向积分可以得到
其中为截面j和j+1之间单元底部的平均应变影响线与x轴围成的面积,只和结构的局部刚度相关,是结构的本质属性,和外部荷载等参数无关,且可以表达为
其中是与位置、距离、中和轴高度相关的函数,其中为截面j和j+1之间的平均刚度;
公式(3)左边进一步表示为
其中v为移动荷载的速度,t0为第一个轴刚进入结构的时刻,tn为最后一个轴,即第n个轴,刚离开结构的时刻,为截面j和j+1之间单元底部的平均应变时程的面积,其中横坐标为时间,纵坐标为应变;
把公式(1)和公式(4)代入公式(5),得到
步骤2:计算各个测试点的位移响应时程面积,代入公式(10)计算实测的位移响应函数比值向量St;
令位移响应函数
其中为截面j的位移时程的面积,其中横坐标为时间,纵坐标为位移;
同理,参考点的位移响应函数表示为
则目标位移响应函数相对参考位移响应函数比值为
从公式(9)看出,位移响应函数比值只和局部刚度以及位置参数相关,位移响应函数比值向量
步骤3:建立初始有限元模型,利用公式(10)计算移动荷载作用下位移响应函数比值向量Sa,其中荷载大小及位置同实测中使用的移动荷载;
步骤4:代入公式(11),计算位移响应函数比值的置信准则MAC,若相关性好,则不需要修正;如相关性不好,建立目标函数
定义位移响应函数比值的置信准则
其中Sa和St分别为结构在已知的移动荷载作用下位移响应函数比值向量的有限元模型计算值和实测值,如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全相关,则MAC=1.0;如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全不相关,则MAC=0;
建立目标函数
步骤5:选择单元局部刚度和边界条件作为修正变量,利用一阶优化算法(或其他优化算法),进行公式(13)求解,当迭代误差小于设定值,则终止计算;
至此,有限元模型修正过程就转化为有约束条件下的优化求解过程,即利用优化算法,通过不断迭代结构的设计参数,使目标函数最小化。
其中g和k是设计参数,分别表示单元平均刚度和边界条件(转动弹簧刚度或竖向弹簧刚度),在迭代过程中设计参数都根据实际情况选取上下限,让迭代过程更快收敛,优化算法一般选择一阶优化算法;
步骤6:把目标函数最优值所对应的结构局部刚度值及边界条件(转动弹簧刚度或竖向弹簧刚度)代入初始有限元模型,则可以得到修正后的有限元模型。
本发明提出利用车辆荷载下的位移响应时程面积构建目标函数进行结构有限元模型修正,理论推导表明,基于位移响应时程面积的目标函数和结构的局部刚度一一对应,因此可以选取局部刚度和边界条件作为修正参数,目标函数和修正变量意义明确,有助于提高计算效率并易于收敛,该方法不需要桥梁交通中断,不需要提取结构的动力特性(即不需要进行时频转换),因此可避免现有方法的不足。
本发明的有益效果:
(1)该方法利用位移传感器响应构建目标函数进行有限元模型修正,位移传感器是目前最常用的传感器之一,和其他传感器相比,位移传感器的精度和可靠性可以得到保证;
(2)相比于其他基于振动特性的有限元模型修正方法,该方法不需要经过傅立叶变换把时域数据转换成频域数据,因此可以避免时频转换误差;
(3)相比于传统的静态有限元模型修正方法,该方法利用移动荷载下的响应进行有限元模型修正,不需要中断交通;
(4)该方法构建的目标函数和结构的局部刚度具有一一对应的关系,目标函数和修正变量意义明确,有助于提高计算效率并易于收敛。
附图说明
图1为桥梁在移动荷载作用下示意图;
图2为移动荷载下结构有限元模型修正流程。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。
如图2所示,一种车辆荷载下基于位移时程面积的结构有限元模型修正方法,包括以下步骤:
步骤1:在目标结构的关键区域准分布式布置位移传感器,测试已知移动荷载作用下各个测试点的位移响应时程dj(t);
对于桥梁结构,跨度为l,梁高度为H,截面j-1,j,j+1,j+2沿桥梁长度方向的坐标分别为xj-1,xj,xj+1,xj+2,假设截面j-1,j,j+1,j+2等间隔,间距为L,移动荷载经过桥梁的整个过程中,截面j-1,j,j+1,j+2处的竖向位移响应为分别为dj-1(x),dj(x),dj+1(x),dj+2(x)(各截面随着移动荷载位移x变化所对应的位移),与之对应的截面j-1,j,j+1,j+2处的位移响应时程为dj-1(t),dj(t),dj+1(t),dj+2(t)(各截面随着时间t变化所对应的位移),其中移动荷载的参数如下:共有n个轴,轴重分别为P1,P2…Pi,Pn,速度为v;
假设结构符合欧拉梁假定,则j截面和j+1截面之间单元的底部平均应变εj,j+1(t)表达为
移动荷载作用下,截面j和j+1之间单元底部的平均应变表达为
式中di(i=1~n)为移动荷载的第i个轴距第1个轴之间的距离,其中d1=0,fj,j+1(x)为截面j和j+1之间单元底部的平均应变影响线,x为第1个轴距左边支座的距离;
把公式(2)左右部分分别沿着结构长度方向积分可以得到
其中为截面j和j+1之间单元底部的平均应变影响线与x轴围成的面积,只和结构的局部刚度相关,是结构的本质属性,和外部荷载等参数无关,且可以表达为
其中是与位置、距离、中和轴高度相关的函数,其中为截面j和j+1之间的平均刚度;
公式(3)左边进一步表示为
其中v为移动荷载的速度,t0为第一个轴刚进入结构的时刻,tn为最后一个轴,即第n个轴,刚离开结构的时刻,为截面j和j+1之间单元底部的平均应变时程的面积,其中横坐标为时间,纵坐标为应变;
把公式(1)和公式(4)代入公式(5),得到
步骤2:计算各个测试点的位移响应时程面积,代入公式(10)计算实测的位移响应函数比值向量St;
令位移响应函数
其中为截面j的位移时程的面积,其中横坐标为时间,纵坐标为位移;
同理,参考点的位移响应函数表示为
则目标位移响应函数相对参考位移响应函数比值为
从公式(9)看出,位移响应函数比值只和局部刚度以及位置参数相关,位移响应函数比值向量
步骤3:建立初始有限元模型,利用公式(10)计算移动荷载作用下位移响应函数比值向量Sa,其中荷载大小及位置同实测中使用的移动荷载;
步骤4:代入公式(11),计算位移响应函数比值的置信准则MAC,若相关性好,则不需要修正;如相关性不好,建立目标函数
定义位移响应函数比值的置信准则
其中Sa和St分别为结构在已知的移动荷载作用下位移响应函数比值向量的有限元模型计算值和实测值,如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全相关,则MAC=1.0;如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全不相关,则MAC=0;
建立目标函数
步骤5:选择单元局部刚度和边界条件作为修正变量,利用一阶优化算法(或其他优化算法),进行公式(13)求解,当迭代误差小于设定值,则终止计算;
至此,有限元模型修正过程就转化为有约束条件下的优化求解过程,即利用优化算法,通过不断迭代结构的设计参数,使目标函数最小化。
其中g和k是设计参数,分别表示单元平均刚度和边界条件(转动弹簧刚度或竖向弹簧刚度),在迭代过程中设计参数都根据实际情况选取上下限,让迭代过程更快收敛,优化算法一般选择一阶优化算法;
步骤6:把目标函数最优值所对应的结构局部刚度值及边界条件(转动弹簧刚度或竖向弹簧刚度)代入初始有限元模型,则可以得到修正后的有限元模型。
应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (1)
1.一种车辆荷载下基于位移时程面积的结构有限元模型修正方法,其特征在于:包括以下步骤:
步骤1:在目标结构的关键区域准分布式布置位移传感器,测试已知移动荷载作用下各个测试点的位移响应时程dj(t);
对于桥梁结构,跨度为l,梁高度为H,截面j-1,j,j+1,j+2沿桥梁长度方向的坐标分别为xj-1,xj,xj+1,xj+2,假设截面j-1,j,j+1,j+2等间隔,间距为L,移动荷载经过桥梁的整个过程中,截面j-1,j,j+1,j+2处的竖向位移响应分别为dj-1(x),dj(x),dj+1(x),dj+2(x),其为各截面随着移动荷载位移x变化所对应的位移,与之对应的截面j-1,j,j+1,j+2处的位移响应时程为dj-1(t),dj(t),dj+1(t),dj+2(t),其为各截面随着时间t变化所对应的位移,其中移动荷载的参数如下:共有n个轴,轴重分别为P1,P2…Pi…Pn,速度为v;
假设结构符合欧拉梁假定,则j截面和j+1截面之间单元的底部平均应变εj,j+1(t)表达为
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>H</mi>
<mrow>
<mn>4</mn>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
移动荷载作用下,截面j和j+1之间单元底部的平均应变表达为
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中di为移动荷载的第i个轴距第1个轴之间的距离,其中i=1~n,d1=0,fj,j+1(x)为截面j和j+1之间单元底部的平均应变影响线,x为第1个轴距左边支座的距离;
把公式(2)左右部分分别沿着结构长度方向积分可以得到
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mi>l</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
</mrow>
</msubsup>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mi>l</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>l</mi>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中为截面j和j+1之间单元底部的平均应变影响线与x轴围成的面积,只和结构的截面j和j+1之间的平均刚度相关,是结构的本质属性,和外部荷载参数无关,且可以表达为
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>l</mi>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mover>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中是与位置、距离、中和轴高度相关的函数,其中为截面j和j+1之间的平均刚度;
公式(3)左边进一步表示为
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mi>l</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
</mrow>
</msubsup>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mi>v</mi>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
</msubsup>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>l</mi>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中v为移动荷载的速度,t0为第一个轴刚进入结构的时刻,tn为最后一个轴,即第n个轴,刚离开结构的时刻,为截面j和j+1之间单元底部的平均应变时程的面积,其中横坐标为时间,纵坐标为应变;
把公式(1)和公式(4)代入公式(5),得到
<mrow>
<mi>v</mi>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
</msubsup>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mi>H</mi>
</mrow>
<mrow>
<mn>4</mn>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mover>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤2:计算各个测试点的位移响应时程面积,代入公式(10)计算实测的位移响应函数比值向量St;
令目标位移响应函数
其中为截面j的位移时程的面积,其中横坐标为时间,纵坐标为位移;
同理,参考位移响应函数表示为
<mrow>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>v</mi>
<mi>H</mi>
<msub>
<mover>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
则目标位移响应函数相对参考位移响应函数比值为
<mrow>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>L</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
从公式(9)看出,位移响应函数比值只和局部刚度以及位置参数相关,位移响应函数比值向量
<mrow>
<mi>S</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>{</mo>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
<mo>...</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>{</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<mo>...</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>...</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>}</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3:建立初始有限元模型,利用公式(10)计算移动荷载作用下位移响应函数比值向量Sa,其中荷载大小及位置同实测中使用的移动荷载;
步骤4:代入公式(11),计算位移响应函数比值的置信准则MAC,若相关性好,则不需要修正;如相关性不好,建立目标函数
定义位移响应函数比值的置信准则
<mrow>
<mi>M</mi>
<mi>A</mi>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msup>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Sa和St分别为结构在已知的移动荷载作用下位移响应函数比值向量的有限元模型计算值和实测值,如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全相关,则MAC=1.0;如果实测位移响应函数比值向量和分析位移响应函数比值向量两者完全不相关,则MAC=0;
建立目标函数
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
2
步骤5:选择单元局部刚度和边界条件作为修正变量,利用一阶优化算法,进行公式(13)求解,当迭代误差小于设定值,则终止计算;
至此,有限元模型修正过程就转化为有约束条件下的优化求解过程,即利用优化算法,通过不断迭代结构的设计参数,使目标函数最小化;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>min</mi>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>g</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
其中g和h是设计参数,分别表示单元平均刚度和边界条件;
步骤6:把目标函数最优值所对应的结构局部刚度值及边界条件代入初始有限元模型,则可以得到修正后的有限元模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510042275.6A CN104573274B (zh) | 2015-01-27 | 2015-01-27 | 车辆荷载下基于位移时程面积的结构有限元模型修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510042275.6A CN104573274B (zh) | 2015-01-27 | 2015-01-27 | 车辆荷载下基于位移时程面积的结构有限元模型修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104573274A CN104573274A (zh) | 2015-04-29 |
CN104573274B true CN104573274B (zh) | 2017-11-24 |
Family
ID=53089327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510042275.6A Expired - Fee Related CN104573274B (zh) | 2015-01-27 | 2015-01-27 | 车辆荷载下基于位移时程面积的结构有限元模型修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104573274B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930571B (zh) * | 2016-04-15 | 2019-02-05 | 东南大学 | 基于单位温度响应监测值的大跨钢桥有限元模型修正方法 |
CN105956218B (zh) * | 2016-04-15 | 2019-05-31 | 东南大学 | 基于非均匀温度响应监测值的钢桥有限元模型修正方法 |
CN105956216B (zh) * | 2016-04-15 | 2019-05-31 | 东南大学 | 大跨钢桥基于均匀温度响应监测值的有限元模型修正方法 |
CN106706239B (zh) * | 2016-11-30 | 2020-02-07 | 山西省交通建设工程质量检测中心(有限公司) | 一种桥梁快速荷载试验测试方法 |
CN108009345B (zh) * | 2017-11-29 | 2021-07-20 | 大连交通大学 | 基于等效移动载荷法的结构拓扑优化方法 |
CN110077494B (zh) * | 2019-04-09 | 2020-08-21 | 爱驰汽车有限公司 | 一种车梁位置移动方法及装置 |
CN110779611B (zh) * | 2019-05-24 | 2020-10-02 | 南京航空航天大学 | 一种超声手术刀刀杆纵向振动频率的校准方法及系统 |
CN112362149B (zh) * | 2020-09-21 | 2022-01-18 | 中铁第四勘察设计院集团有限公司 | 基于竖向位移影响面加载动态识别车辆轴重方法和系统 |
CN113392451B (zh) * | 2021-06-09 | 2022-05-17 | 哈尔滨工业大学 | 基于车-桥梁耦合作用力修正的桥梁模型更新方法、系统、存储介质及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927428A (zh) * | 2014-05-09 | 2014-07-16 | 清华大学 | 一种考虑多因素影响的锥齿轮错位量有限元计算方法 |
CN104090977A (zh) * | 2014-07-22 | 2014-10-08 | 东南大学 | 一种桥面移动车辆载荷的随机识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0414992D0 (en) * | 2004-07-02 | 2004-08-04 | Engenuity Ltd | Crush modelling |
-
2015
- 2015-01-27 CN CN201510042275.6A patent/CN104573274B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927428A (zh) * | 2014-05-09 | 2014-07-16 | 清华大学 | 一种考虑多因素影响的锥齿轮错位量有限元计算方法 |
CN104090977A (zh) * | 2014-07-22 | 2014-10-08 | 东南大学 | 一种桥面移动车辆载荷的随机识别方法 |
Non-Patent Citations (2)
Title |
---|
Detailed load rating analyses of bridge populations using nonlinear finite element models and artificial neural networks;Hasançebi O 等;《Computers & Structures》;20131231(第128期);第48-63页 * |
随机实测车辆荷载下大跨径斜拉桥钢箱梁的动力响应特征;杨沐野 等;《河海大学学报自然科学版》;20140930;第42卷(第5期);第427-432页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104573274A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104573274B (zh) | 车辆荷载下基于位移时程面积的结构有限元模型修正方法 | |
CN101539477B (zh) | 一种深水顶张式立管涡激振动与疲劳分析的方法 | |
CN102735386B (zh) | 考虑弯曲刚度的斜拉索索力数值计算方法 | |
Gu et al. | Direct identification of flutter derivatives and aerodynamic admittances of bridge decks | |
CN104866666A (zh) | 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法 | |
CN105260568B (zh) | 基于离散型卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN107678282B (zh) | 考虑未知动力学和外部干扰的mems陀螺智能控制方法 | |
CN103400035B (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN110017929B (zh) | 基于子结构灵敏度分析的船撞桥荷载与损伤同步识别方法 | |
CN104992002B (zh) | 一种面向智能蒙皮天线的应变传感器布局方法 | |
CN106200378B (zh) | 鲁棒非脆弱h∞结构振动控制方法 | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
CN104807414A (zh) | 基于分布式光纤传感技术的地铁隧道沉降变形监测方法 | |
CN109444350A (zh) | 一种基于无人机的大气污染物监测传感器的布局方法 | |
CN104568493A (zh) | 车辆荷载下基于位移时程面积的结构快速损伤识别方法 | |
CN107356523B (zh) | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 | |
CN106503332B (zh) | 一种引入非线性弹簧的微型桩基-土动力响应求解方法 | |
CN105403344A (zh) | 管道实时应力的获取方法 | |
CN104408295A (zh) | 一种大跨桥梁下部结构风-浪耦合作用荷载数值模拟方法 | |
CN102566446B (zh) | 基于线性模型组的无人直升机全包线数学模型构建方法 | |
CN105117539A (zh) | 风力机叶片模态频率及其双峰高斯分布拟合方法 | |
CN104794332B (zh) | 一种高层建筑风致响应分析模型的不确定性分析方法 | |
CN110717216A (zh) | 不规则波下带柔性气囊直升机横摇响应预报方法 | |
CN109398020A (zh) | 一种基于非线性模型的车辆液电耦合式isd悬架的预测控制方法 | |
CN104155043B (zh) | 一种动力定位系统外界环境力测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171124 Termination date: 20210127 |