CN115828379B - 一种风-汽车-列车-桥梁耦合振动分析方法 - Google Patents
一种风-汽车-列车-桥梁耦合振动分析方法 Download PDFInfo
- Publication number
- CN115828379B CN115828379B CN202211434913.5A CN202211434913A CN115828379B CN 115828379 B CN115828379 B CN 115828379B CN 202211434913 A CN202211434913 A CN 202211434913A CN 115828379 B CN115828379 B CN 115828379B
- Authority
- CN
- China
- Prior art keywords
- subsystem
- train
- bridge
- automobile
- tsr
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 84
- 230000008878 coupling Effects 0.000 title claims abstract description 63
- 238000010168 coupling process Methods 0.000 title claims abstract description 63
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 63
- 230000005284 excitation Effects 0.000 claims abstract description 23
- 238000006073 displacement reaction Methods 0.000 claims description 49
- 230000001133 acceleration Effects 0.000 claims description 47
- 239000011159 matrix material Substances 0.000 claims description 30
- 230000003993 interaction Effects 0.000 claims description 19
- 230000004044 response Effects 0.000 claims description 19
- 238000013016 damping Methods 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000000354 decomposition reaction Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 4
- 230000008846 dynamic interplay Effects 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000001808 coupling effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- 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
- Bridges Or Land Bridges (AREA)
Abstract
本发明公开了一种风‑汽车‑列车‑桥梁耦合振动分析方法,包括:建立系统模型,模拟外部激励,定义时间步长,利用域分解技术将复杂的风‑汽车‑列车‑桥梁系统划分为列车、汽车和桥梁三个子系统,并根据第i个时间步长中的i能否整除对应子系统时间步长与单位的比值对三个子系统各自独立进行动力特性分析,获得三个子系统对应的振动分析结果,并耦合成为风‑汽车‑列车‑桥梁耦合振动分析的结果。可扩展的多时间步算法为汽车子系统、列车子系统和桥梁子系统分别分配了不同的时间步长,极大地提高了求解风‑汽车‑列车‑桥梁系统动态交互问题的计算效率,弥补了现有技术在风‑汽车‑列车‑桥梁系统分析方面的空缺。
Description
技术领域
本发明涉及桥梁安全技术领域,具体为一种风-汽车-列车-桥梁耦合振动分析方法。
背景技术
公铁两用桥是一种可同时容纳同一层或两层不同的公路和铁路的桥梁。随着中国经济的不断发展和对公路、铁路运输需求的不断增长,在过去的二十年中,许多铁路公路结合的桥梁相继建成,如天兴洲长江大桥、苏通长江大桥、平潭海峡大桥、五峰山大桥和规划中的常泰长江大桥。由于公铁两用桥通常建于风多发地区,因此在风、汽车和列车的联合荷载作用下,确保桥梁的安全性和功能是非常重要的。
然而,目前公铁两用桥的处理非常具有挑战性,因为它需要充分考虑桥梁、汽车和列车三个子系统之间相互作用,这比传统公路桥或铁路桥要复杂得多。事实上,目前的研究大多集中在风荷载作用下的车-桥或列车-桥系统,而很少有研究关注风-汽车-列车-桥的相互作用,并且求解风-汽车-列车-桥梁系统计算成本高。
发明内容
针对现有技术中的上述不足,本发明提供的风-汽车-列车-桥梁耦合振动分析方法解决了桥梁安全技术领域风-汽车-列车-桥系统研究方向空缺的问题。
为了达到上述发明目的,本发明采用的技术方案为:提供一种风-汽车-列车-桥梁耦合振动分析方法,其包括以下步骤:
S1、建立汽车、列车和桥梁子系统的数值模型;
S2、模拟风-汽车-列车-桥梁系统的动态激励源,将其作为风-汽车-列车-桥梁系统的外部激励输入;
S3、根据汽车、列车和桥梁子系统的数值模型和风-汽车-列车-桥梁系统的外部激励输入构建风-汽车-列车-桥梁耦合振动分析系统;
S4、将风-汽车-列车-桥梁耦合振动系统划分为桥梁子系统、汽车子系统和列车子系统,定义汽车子系统的时间步长Δtr、列车子系统的时间步长Δtt、桥梁子系统的时间步长Δtb,单位时间步长dt,并分别计算汽车、列车和桥梁子系统的时间步长与dt的比值TSRr、TSRt和TSRb;
S5、计算第i个时间的步长ti,并判断i是否整除TSRt:
若是,则进入步骤S6;
若否,则进入步骤S8;
S6、进行列车子系统第i个时间步的动力响应分析,获得列车子系统耦合振动分析结果;
S7、根据列车子系统耦合震振动分析结果,判断是否满足列车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S8;
若否,则返回步骤S6;
S8、判断i是否整除TSRr:
若是,则进入步骤S9;
若否,则进入步骤S11;
S9、进行汽车子系统第i个时间步的动力响应分析,获得汽车子系统耦合振动分析结果;
S10、根据汽车子系统耦合震振动分析结果,判断是否满足汽车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S11;
若否,则返回步骤S9;
S11、判断i是否整除TSRb:
若是,则进入步骤S12;
若否,则将i的值加1,并返回步骤S5;
S12、进行桥梁子系统第i个时间步的动力响应分析,并假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变,获得桥梁子系统耦合振动分析结果;
S13、判断是否所有汽车和列车都驶离桥梁:
若是,则将列车子系统耦合振动分析结果、汽车子系统耦合振动分析结果和桥梁子系统耦合振动分析结果作为风-汽车-列车-桥梁耦合振动分析结果,并结束分析;
若否,则将i的值加1,并返回步骤S5。
进一步地:所述步骤S3中,所述风-汽车-列车-桥梁耦合振动分析系统由风-汽车-列车-桥梁耦合振动方程表示,其表达式为:
其中,Mr为汽车子系统的质量矩阵,Mt为列车子系统的质量矩阵,Mb为桥梁子系统的质量矩阵,Kr为汽车子系统的刚度矩阵,Kt为列车子系统的刚度矩阵,Kb为桥梁子系统的刚度矩阵,Cr为汽车子系统的阻尼矩阵,Ct为列车子系统的阻尼矩阵,Cb为桥梁子系统的阻尼矩阵,Xr为汽车子系统的位移向量,Xt为列车子系统的位移向量,Xb为桥梁子系统的位移向量,为汽车子系统的速度向量,/>为列车子系统的速度向量,/>为桥梁子系统的速度向量,/>为汽车子系统的加速度向量,/>为列车子系统的加速度向量,/>为桥梁子系统的加速度向量,Frb和Fbr分别表示汽车作用在桥梁上的力向量和桥梁作用在汽车上的力向量;Ftb和Fbt分别表示列车作用在桥梁上的力向量和桥梁作用在列车上的力向量;Fbw表示风作用在桥梁的力向量;Frw和Ftw分别表示风作用在汽车上的力向量和风作用在列车上的力向量。
进一步地:所述Frb和Frw表示为:
Frw=FrwST(Xr)+FrwBU(Xr)
其中,FrbR(.)表示路面粗糙度引起的作用在汽车上的激励力;Frbp(.)表示桥梁变形引起的作用在汽车上的附加力;FrwST(.)表示风荷载作用在汽车上的静风力;FrwBU(.)表示风荷载作用在汽车上的抖振风力。
进一步地:所述Ftb和Ftw表示为:
Ftw=FtwST(Xt)+FtwBU(Xt)
其中,FtbN(.)表示作用在列车上的轮轨法向力;FtbRH(.)表示作用在列车上的蠕滑力;FtwST(.)表示风荷载作用在列车上的静风力;FtwBU(.)表示风荷载作用在列车上的抖振风力。
进一步地:所述Fbr、Fbt和Fbw表示为:
Fbw=FbwST(Xb)+FbwBU(Xb)+FbwSE(Xb)
其中,FbrR(.)表示路面粗糙度引起的汽车作用在桥面上的激励力;FbrP(.)表示桥梁变形引起的汽车作用在桥面上的附加力;FbtN(.)表示列车作用在桥梁上的轮轨法向力;FbtRH(.)表示列车作用在桥梁上的蠕滑力;FbwST(.)表示风荷载作用在桥梁结构上的静风力;FbwBU(.)表示风荷载作用在桥梁结构上的抖振风力;FbwSE(.)表示作用在桥梁结构上的自激风力。
进一步地:所述步骤S4中,所述汽车子系统、列车子系统和桥梁子系统的时间步长表示为:
Δtr=TSRrdt
Δtt=TSRtdt
Δtb=TSRbdt
其中,Δtr、Δtt和Δtb分别为汽车子系统、列车子系统和桥梁子系统的时间步长;dt是单位时间步长;TSRr、TSRt和TSRb定义为汽车、列车和桥梁子系统的时间步长与dt的比值,且TSRt≤TSRr≤TSRb。
上述进一步方案的有益效果为:建立的风-汽车-列车-桥梁耦合振动系统科学合理、概念清晰,符合大跨桥梁在内外部激励下桥梁子系统、汽车子系统和列车子系统间的耦合作用状况。
进一步地:所述步骤S6中,列车子系统第i个时间步的动力响应分析包括以下分步骤:
S61、获取列车在桥梁上的当前位置;
S62、计算列车在第i个时间步下风作用在列车上的力和列车作用在桥上的力
S63、根据和/>计算列车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为列车子系统在第i个时间步下的位移,/>为列车子系统在第i个时间步下的速度,/>为列车子系统在第i个时间步下的加速度。
S64、将和/>在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间列车子系统的加速度、速度和位移,并将其作为列车子系统耦合振动分析结果,其表达式为:
其中,为列车子系统在第i-TSRt个时间步下的位移,/>为列车子系统在第i-TSRt个时间步下的速度,/>为列车子系统在第i-TSRt个时间步下的加速度,为列车子系统在第i-TSRt+m个时间步下的位移,/>为列车子系统在第i-TSRt+m个时间步下的速度,/>为列车子系统在第i-TSRt+m个时间步下的加速度,m=1,2,…TSRt-1。
进一步地:所述步骤S9中,汽车子系统第i个时间步的动力响应分析包括以下分步骤:
S91、获取汽车在桥梁上的当前位置;
S62、计算汽车在第i个时间步下风作用在列车上的力和列车作用在桥上的力
S93、根据和/>计算汽车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为汽车子系统在第i个时间步下的位移,/>为汽车子系统在第i个时间步下的速度,/>为汽车子系统在第i个时间步下的加速度。
S94、将和/>在区间[(i-TSRr)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,并将其作为汽车子系统耦合振动分析结果,其表达式为:
其中,为汽车子系统在第i-TSRr个时间步下的位移,/>为汽车子系统在第i-TSRr个时间步下的速度,/>为汽车子系统在第i-TSRr个时间步下的加速度,为汽车子系统在第i-TSRr+m个时间步下的位移,/>为汽车子系统在第i-TSRr+m个时间步下的速度,/>为汽车子系统在第i-TSRr+m个时间步下的加速度,m=1,2,…TSRr-1。
进一步地:所述步骤S12中,桥梁子系统第i个时间步的动力响应分析包括以下分步骤:
S121、计算第i个时间步下桥梁作用在列车上的力和桥梁作用在汽车上的力
S122、计算第i个时间步下风作用在桥梁上的力
S123、根据和/>计算桥梁子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为桥梁子系统在第i个时间步下的位移,/>为桥梁子系统在第i个时间步下的速度,/>为桥梁子系统在第i个时间步下的加速度。
S124、通过在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,其表达式为:
其中,为桥梁子系统在第i-TSRb个时间步下的位移,/>为桥梁子系统在第i-TSRb个时间步下的速度,/>为桥梁子系统在第i-TSRb个时间步下的加速度,为桥梁子系统在第i-TSRb+m个时间步下的位移,/>为桥梁子系统在第i-TSRb+m个时间步下的速度,/>为桥梁子系统在第i-TSRb+m个时间步下的加速度,m=1,2,…TSRb-1。
进一步地:所述步骤S12中,假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变的表达式为:
其中,和/>均与振动相关,且在每个时间步开始时未知。
上述进一步方案的有益效果为:利用域分解技术将复杂的风-汽车-列车-桥梁系统划分为三个子系统,并为汽车子系统、列车子系统和桥梁子系统分别分配了不同的时间步长,极大地提高了求解风-汽车-列车-桥梁系统的计算效率。
本发明的有益效果为:
1.建立了一个风-汽车-列车-桥梁动力耦合振动系统,为桥梁安全技术领域提供了一种全新的分析方法;
2.建立的风-汽车-列车-桥梁耦合振动系统科学合理、概念清晰,符合大跨桥梁在内外部激励下桥梁子系统、汽车子系统和列车子系统间的耦合作用状况;
3.基于可扩展的多时间步方法利用域分解技术将复杂的风-汽车-列车-桥梁系统划分为三个子系统,即低频率的桥梁子系统和高频率的汽车子系统和列车子系统,并为汽车子系统、列车子系统和桥梁子系统分别分配了不同的时间步长,极大地提高了求解风-汽车-列车-桥梁系统的计算效率。
附图说明
图1为本发明提供的基于可扩展的多时间步算法的风-汽车-列车-桥梁分析方法的流程图。
图2为汽车、列车和大跨桥梁的各子系统划分方式示意图。
图3为可扩展的多时间步方法示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,在本发明的一个实施例中,包括以下步骤:
S1、建立汽车、列车和桥梁子系统的数值模型;
S2、模拟风-汽车-列车-桥梁系统的动态激励源,将其作为风-汽车-列车-桥梁系统的外部激励输入;
S3、根据汽车、列车和桥梁子系统的数值模型和风-汽车-列车-桥梁系统的外部激励输入构建风-汽车-列车-桥梁耦合振动分析系统;
S4、将风-汽车-列车-桥梁耦合振动系统划分为桥梁子系统、汽车子系统和列车子系统,定义汽车子系统的时间步长Δtr、列车子系统的时间步长Δtt、桥梁子系统的时间步长Δtb,单位时间步长dt,并分别计算汽车、列车和桥梁子系统的时间步长与dt的比值TSRr、TSRt和TSRb;
S5、计算第i个时间的步长ti,并判断i是否整除TSRt:
若是,则进入步骤S6;
若否,则进入步骤S8;
S6、进行列车子系统第i个时间步的动力响应分析,获得列车子系统耦合振动分析结果;
S7、根据列车子系统耦合震振动分析结果,判断是否满足列车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S8;
若否,则返回步骤S6;
S8、判断i是否整除TSRr:
若是,则进入步骤S9;
若否,则进入步骤S11;
S9、进行汽车子系统第i个时间步的动力响应分析,获得汽车子系统耦合振动分析结果;
S10、根据汽车子系统耦合震振动分析结果,判断是否满足汽车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S11;
若否,则返回步骤S9;
S11、判断i是否整除TSRb:
若是,则进入步骤S12;
若否,则将i的值加1,并返回步骤S5;
S12、进行桥梁子系统第i个时间步的动力响应分析,并假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变,获得桥梁子系统耦合振动分析结果;
S13、判断是否所有汽车和列车都驶离桥梁:
若是,则将列车子系统耦合振动分析结果、汽车子系统耦合振动分析结果和桥梁子系统耦合振动分析结果作为风-汽车-列车-桥梁耦合振动分析结果,并结束分析;
若否,则将i的值加1,并返回步骤S5。
具体地,本实施例中步骤S1中汽车的数值模型包括若干个刚体、阻尼器、弹簧、悬挂系统和轮轴,将汽车模型模拟为一个质量-弹簧-阻尼系统;列车的数值模型包括若干个刚体、阻尼器、弹簧和悬挂系统,将列车模型模拟为一个质量-弹簧-阻尼系统;桥梁的数值模型通过有限元法建立,主梁、桥塔、桥墩和基础采用三维梁单元模拟,斜拉索采用空间杆单元模拟,附属结构及二期恒载通过施加质量单元来模拟;相互作用中的风和路面粗糙度等激励通过谐波合成法来模拟。
如图2所示,在本发明的一个实施例中,所述风-汽车-列车-桥梁耦合振动分析系统包括:汽车-桥相互作用模型、列车-桥相互作用模型、风-桥相互作用模型、风-汽车相互作用模和风-列车相互作用模型;
所述汽车-桥相互作用模型包括桥面和轮胎之间的横向接触力和竖向接触力,所述横向接触力分为路面粗糙度引起的激励力和桥梁变形引起的附加力;列车-桥相互作用模型包括轮轨法向力和轮轨蠕滑力;风-桥相互作用模型包括作用在桥梁结构上的静风力,自激风力和抖振风力;风-汽车相互作用模型包括静风力和抖振风力;风-列车相互模型作用包括静风力和抖振风力。
本实施例步骤S3中,所述风-汽车-列车-桥梁耦合振动分析系统由风-汽车-列车-桥梁耦合振动方程表示,其表达式为:
其中,Mr为汽车子系统的质量矩阵,Mt为列车子系统的质量矩阵,Mb为桥梁子系统的质量矩阵,Kr为汽车子系统的刚度矩阵,Kt为列车子系统的刚度矩阵,Kb为桥梁子系统的刚度矩阵,Cr为汽车子系统的阻尼矩阵,Ct为列车子系统的阻尼矩阵,Cb为桥梁子系统的阻尼矩阵,Xr为汽车子系统的位移向量,Xt为列车子系统的位移向量,Xb为桥梁子系统的位移向量,为汽车子系统的速度向量,/>为列车子系统的速度向量,/>为桥梁子系统的速度向量,/>为汽车子系统的加速度向量,/>为列车子系统的加速度向量,/>为桥梁子系统的加速度向量,Frb和Fbr分别表示汽车作用在桥梁上的力向量和桥梁作用在汽车上的力向量;Ftb和Fbt分别表示列车作用在桥梁上的力向量和桥梁作用在列车上的力向量;Fbw表示风作用在桥梁的力向量;Frw和Ftw分别表示风作用在汽车上的力向量和风作用在列车上的力向量。
所述Frb和Frw表示为:
Frw=FrwST(Xr)+FrwBU(Xr)
其中,FrbR(.)表示路面粗糙度引起的作用在汽车上的激励力;Frbp(.)表示桥梁变形引起的作用在汽车上的附加力;FrwST(.)表示风荷载作用在汽车上的静风力;FrwBU(.)表示风荷载作用在汽车上的抖振风力。
所述Ftb和Ftw表示为:
Ftw=FtwST(Xt)+FtwBU(Xt)
其中,FtbN(.)表示作用在列车上的轮轨法向力;FtbRH(.)表示作用在列车上的蠕滑力;FtwST(.)表示风荷载作用在列车上的静风力;FtwBU(.)表示风荷载作用在列车上的抖振风力。
所述Fbr、Fbt和Fbw表示为:
Fbw=FbwST(Xb)+FbwBU(Xb)+FbwSE(Xb)
其中,FbrR(.)表示路面粗糙度引起的汽车作用在桥面上的激励力;FbrP(.)表示桥梁变形引起的汽车作用在桥面上的附加力;FbtN(.)表示列车作用在桥梁上的轮轨法向力;FbtRH(.)表示列车作用在桥梁上的蠕滑力;FbwST(.)表示风荷载作用在桥梁结构上的静风力;FbwBU(.)表示风荷载作用在桥梁结构上的抖振风力;FbwSE(.)表示作用在桥梁结构上的自激风力。
如图3所示,在本发明的一个实施例中,基于利用域分解技术将复杂的风-汽车-列车-桥梁系统划分为三个子系统,即低频率的桥梁子系统和高频率的汽车子系统和列车子系统,为汽车子系统、列车子系统和桥梁子系统分别分配了不同的时间步长,极大地提高了求解风-汽车-列车-桥梁系统的计算效率。
所述可拓展的多时间步算法的求解过程为:
在进行风-汽车-列车-桥梁耦合系统动力响应分析之前,首先定义汽车子系统的时间步长Δtr、列车子系统的时间步长Δtt和桥梁子系统的时间步长Δtb:
Δtr=TSRrdt
Δtt=TSRtdt
Δtb=TSRbdt
其中,Δtr、Δtt和Δtb分别为汽车子系统、列车子系统和桥梁子系统的时间步长;dt是单位时间步长;TSRr、TSRt和TSRb定义为汽车、列车和桥梁子系统的时间步长与dt的比值,且TSRt≤TSRr≤TSRb。
进一步地,在本发明的一个实施例步骤S6中,列车子系统第i个时间步的动力响应分析包括以下分步骤:
S61、获取列车在桥梁上的当前位置;
S62、计算列车在第i个时间步下风作用在列车上的力和列车作用在桥上的力
S63、根据和/>计算列车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为列车子系统在第i个时间步下的位移,/>为列车子系统在第i个时间步下的速度,/>为列车子系统在第i个时间步下的加速度。
S64、将和/>在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间列车子系统的加速度、速度和位移,并将其作为列车子系统耦合振动分析结果,其表达式为:
其中,为列车子系统在第i-TSRt个时间步下的位移,/>为列车子系统在第i-TSRt个时间步下的速度,/>为列车子系统在第i-TSRt个时间步下的加速度,为列车子系统在第i-TSRt+m个时间步下的位移,/>为列车子系统在第i-TSRt+m个时间步下的速度,/>为列车子系统在第i-TSRt+m个时间步下的加速度,m=1,2,…TSRt-1。
进一步地,在本发明的一个实施例步骤S9中,汽车子系统第i个时间步的动力响应分析包括以下分步骤:
S91、获取汽车在桥梁上的当前位置;
S62、计算汽车在第i个时间步下风作用在列车上的力和列车作用在桥上的力
S93、根据和/>计算汽车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为汽车子系统在第i个时间步下的位移,/>为汽车子系统在第i个时间步下的速度,/>为汽车子系统在第i个时间步下的加速度。
S94、将和/>在区间[(i-TSRr)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,并将其作为汽车子系统耦合振动分析结果,其表达式为:
其中,为汽车子系统在第i-TSRr个时间步下的位移,/>为汽车子系统在第i-TSRr个时间步下的速度,/>为汽车子系统在第i-TSRr个时间步下的加速度,为汽车子系统在第i-TSRr+m个时间步下的位移,/>为汽车子系统在第i-TSRr+m个时间步下的速度,/>为汽车子系统在第i-TSRr+m个时间步下的加速度,m=1,2,…TSRr-1。
进一步地,在本发明的一个实施例步骤S12中,桥梁子系统第i个时间步的动力响应分析包括以下分步骤:
S121、计算第i个时间步下桥梁作用在列车上的力和桥梁作用在汽车上的力/>
S122、计算第i个时间步下风作用在桥梁上的力
S123、根据和/>计算桥梁子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为桥梁子系统在第i个时间步下的位移,/>为桥梁子系统在第i个时间步下的速度,/>为桥梁子系统在第i个时间步下的加速度。
S124、通过在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,其表达式为:
其中,为桥梁子系统在第i-TSRb个时间步下的位移,/>为桥梁子系统在第i-TSRb个时间步下的速度,/>为桥梁子系统在第i-TSRb个时间步下的加速度,为桥梁子系统在第i-TSRb+m个时间步下的位移,/>为桥梁子系统在第i-TSRb+m个时间步下的速度,/>为桥梁子系统在第i-TSRb+m个时间步下的加速度,m=1,2,…TSRb-1。
在本发明实施例步骤S12中,假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变的表达式为:
其中,和/>均与振动相关,且在每个时间步开始时未知。
风-汽车-列车-桥梁耦合系统基于三个子系统之间位移协调和相互作用力的耦合关系,使用分离迭代法独立求解。
在本发明的描述中,需要理解的是,术语“中心”、“厚度”、“上”、“下”、“水平”、“顶”、“底”、“内”、“外”、“径向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性或隐含指明的技术特征的数量。因此,限定由“第一”、“第二”、“第三”的特征可以明示或隐含地包括一个或者更多个该特征。
本发明建立了一个风-汽车-列车-桥梁动力耦合振动系统,其系统科学合理、概念清晰,为桥梁安全技术领域提供了一种全新的分析方法,并为汽车子系统、列车子系统和桥梁子系统分别分配了不同的时间步长,极大地提高了求解风-汽车-列车-桥梁系统的计算效率。
Claims (9)
1.一种风-汽车-列车-桥梁耦合振动分析方法,其特征在于,包括以下步骤:
S1、建立汽车、列车和桥梁子系统的数值模型;
S2、模拟风-汽车-列车-桥梁系统的动态激励源,将其作为风-汽车-列车-桥梁系统的外部激励输入;
S3、根据汽车、列车和桥梁子系统的数值模型和风-汽车-列车-桥梁系统的外部激励输入构建风-汽车-列车-桥梁耦合振动分析系统;
S4、将风-汽车-列车-桥梁耦合振动系统划分为桥梁子系统、汽车子系统和列车子系统,定义汽车子系统的时间步长Δtr、列车子系统的时间步长Δtt、桥梁子系统的时间步长Δtb,单位时间步长dt,并分别计算汽车、列车和桥梁子系统的时间步长与dt的比值TSRr、TSRt和TSRb;
S5、计算第i个时间的步长ti,并判断i是否整除TSRt:
若是,则进入步骤S6;
若否,则进入步骤S8;
S6、进行列车子系统第i个时间步的动力响应分析,获得列车子系统耦合振动分析结果;
所述步骤S6中,列车子系统第i个时间步的动力响应分析包括以下分步骤:
S61、获取列车在桥梁上的当前位置;
S62、计算列车在第i个时间步下风作用在列车上的力和列车作用在桥上的力/>
S63、根据和/>计算列车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为列车子系统在第i个时间步下的位移,/>为列车子系统在第i个时间步下的速度,/>为列车子系统在第i个时间步下的加速度;
S64、将和/>在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间列车子系统的加速度、速度和位移,并将其作为列车子系统耦合振动分析结果,其表达式为:
其中,为列车子系统在第i-TSRt个时间步下的位移,/>为列车子系统在第i-TSRt个时间步下的速度,/>为列车子系统在第i-TSRt个时间步下的加速度,/>为列车子系统在第i-TSRt+m个时间步下的位移,/>为列车子系统在第i-TSRt+m个时间步下的速度,/>为列车子系统在第i-TSRt+m个时间步下的加速度,m=1,2,…TSRt-1;Mt为列车子系统的质量矩阵,Kt为列车子系统的刚度矩阵,Ct为列车子系统的阻尼矩阵;
S7、根据列车子系统耦合震振动分析结果,判断是否满足列车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S8;
若否,则返回步骤S6;
S8、判断i是否整除TSRr:
若是,则进入步骤S9;
若否,则进入步骤S11;
S9、进行汽车子系统第i个时间步的动力响应分析,获得汽车子系统耦合振动分析结果;
S10、根据汽车子系统耦合震振动分析结果,判断是否满足汽车子系统和桥梁子系统之间的收敛准则:
若是,则进入步骤S11;
若否,则返回步骤S9;
S11、判断i是否整除TSRb:
若是,则进入步骤S12;
若否,则将i的值加1,并返回步骤S5;
S12、进行桥梁子系统第i个时间步的动力响应分析,并假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变,获得桥梁子系统耦合振动分析结果;
S13、判断是否所有汽车和列车都驶离桥梁:
若是,则将列车子系统耦合振动分析结果、汽车子系统耦合振动分析结果和桥梁子系统耦合振动分析结果作为风-汽车-列车-桥梁耦合振动分析结果,并结束分析;
若否,则将i的值加1,并返回步骤S5。
2.根据权利要求1所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述步骤S3中,所述风-汽车-列车-桥梁耦合振动分析系统由风-汽车-列车-桥梁耦合振动方程表示,其表达式为:
其中,Mr为汽车子系统的质量矩阵,Mt为列车子系统的质量矩阵,Mb为桥梁子系统的质量矩阵,Kr为汽车子系统的刚度矩阵,Kt为列车子系统的刚度矩阵,Kb为桥梁子系统的刚度矩阵,Cr为汽车子系统的阻尼矩阵,Ct为列车子系统的阻尼矩阵,Cb为桥梁子系统的阻尼矩阵,Xr为汽车子系统的位移向量,Xt为列车子系统的位移向量,Xb为桥梁子系统的位移向量,为汽车子系统的速度向量,/>为列车子系统的速度向量,/>为桥梁子系统的速度向量,为汽车子系统的加速度向量,/>为列车子系统的加速度向量,/>为桥梁子系统的加速度向量,Frb和Fbr分别表示汽车作用在桥梁上的力向量和桥梁作用在汽车上的力向量;Ftb和Fbt分别表示列车作用在桥梁上的力向量和桥梁作用在列车上的力向量;Fbw表示风作用在桥梁的力向量;Frw和Ftw分别表示风作用在汽车上的力向量和风作用在列车上的力向量。
3.根据权利要求2所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述Frb和Frw表示为:
Frw=FrwST(Xr)+FrwBU(Xr)
其中,FrbR(.)表示路面粗糙度引起的作用在汽车上的激励力;Frbp(.)表示桥梁变形引起的作用在汽车上的附加力;FrwST(.)表示风荷载作用在汽车上的静风力;FrwBU(.)表示风荷载作用在汽车上的抖振风力。
4.根据权利要求2所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述Ftb和Ftw表示为:
Ftw=FtwST(Xt)+FtwBU(Xt)
其中,FtbN(.)表示作用在列车上的轮轨法向力;FtbRH(.)表示作用在列车上的蠕滑力;FtwST(.)表示风荷载作用在列车上的静风力;FtwBU(.)表示风荷载作用在列车上的抖振风力。
5.根据权利要求2所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述Fbr、Fbt和Fbw表示为:
Fbw=FbwST(Xb)+FbwBU(Xb)+FbwSE(Xb)
其中,FbrR(.)表示路面粗糙度引起的汽车作用在桥面上的激励力;FbrP(.)表示桥梁变形引起的汽车作用在桥面上的附加力;FbtN(.)表示列车作用在桥梁上的轮轨法向力;FbtRH(.)表示列车作用在桥梁上的蠕滑力;FbwST(.)表示风荷载作用在桥梁结构上的静风力;FbwBU(.)表示风荷载作用在桥梁结构上的抖振风力;FbwSE(.)表示作用在桥梁结构上的自激风力。
6.根据权利要求1所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述步骤S4中,所述汽车子系统、列车子系统和桥梁子系统的时间步长表示为:
Δtr=TSRrdt
Δtt=TSRtdt
Δtb=TSRbdt
其中,Δtr、Δtt和Δtb分别为汽车子系统、列车子系统和桥梁子系统的时间步长;dt是单位时间步长;TSRr、TSRt和TSRb定义为汽车、列车和桥梁子系统的时间步长与dt的比值,且TSRt≤TSRr≤TSRb。
7.根据权利要求2所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述步骤S9中,汽车子系统第i个时间步的动力响应分析包括以下分步骤:
S91、获取汽车在桥梁上的当前位置;
S62、计算汽车在第i个时间步下风作用在列车上的力和列车作用在桥上的力/>
S93、根据和/>计算汽车子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为汽车子系统在第i个时间步下的位移,/>为汽车子系统在第i个时间步下的速度,/>为汽车子系统在第i个时间步下的加速度;
S94、将和/>在区间[(i-TSRr)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,并将其作为汽车子系统耦合振动分析结果,其表达式为:
其中,为汽车子系统在第i-TSRr个时间步下的位移,/>为汽车子系统在第i-TSRr个时间步下的速度,/>为汽车子系统在第i-TSRr个时间步下的加速度,/>为汽车子系统在第i-TSRr+m个时间步下的位移,/>为汽车子系统在第i-TSRr+m个时间步下的速度,/>为汽车子系统在第i-TSRr+m个时间步下的加速度,m=1,2,…TSRr-1。
8.根据权利要求2所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述步骤S12中,桥梁子系统第i个时间步的动力响应分析包括以下分步骤:
S121、计算第i个时间步下桥梁作用在列车上的力和桥梁作用在汽车上的力/>
S122、计算第i个时间步下风作用在桥梁上的力
S123、根据和/>计算桥梁子系统在ti=i*dt时刻的振动方程,其表达式为:
其中,为桥梁子系统在第i个时间步下的位移,/>为桥梁子系统在第i个时间步下的速度,/>为桥梁子系统在第i个时间步下的加速度;
S124、通过在区间[(i-TSRt)×dt,i×dt]进行线性插值,得到第(i-TSRt+1)和(i-1)个时间步之间汽车子系统的加速度、速度和位移,其表达式为:
其中,为桥梁子系统在第i-TSRb个时间步下的位移,/>为桥梁子系统在第i-TSRb个时间步下的速度,/>为桥梁子系统在第i-TSRb个时间步下的加速度,/>为桥梁子系统在第i-TSRb+m个时间步下的位移,/>为桥梁子系统在第i-TSRb+m个时间步下的速度,/>为桥梁子系统在第i-TSRb+m个时间步下的加速度,m=1,2,…TSRb-1。
9.根据权利要求8所述的风-汽车-列车-桥梁耦合振动分析方法,其特征在于:所述步骤S12中,假设汽车子系统、列车子系统和桥梁子系统之间的相互作用力在当前时间步Δtb内保持不变的表达式为:
其中,和/>均与振动相关,且在每个时间步开始时未知。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211434913.5A CN115828379B (zh) | 2022-11-16 | 2022-11-16 | 一种风-汽车-列车-桥梁耦合振动分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211434913.5A CN115828379B (zh) | 2022-11-16 | 2022-11-16 | 一种风-汽车-列车-桥梁耦合振动分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115828379A CN115828379A (zh) | 2023-03-21 |
CN115828379B true CN115828379B (zh) | 2023-10-31 |
Family
ID=85528515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211434913.5A Active CN115828379B (zh) | 2022-11-16 | 2022-11-16 | 一种风-汽车-列车-桥梁耦合振动分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115828379B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657117A (zh) * | 2017-09-26 | 2018-02-02 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种公铁两用桥梁车‑桥‑风浪流耦合振动分析方法 |
CN108268711A (zh) * | 2018-01-04 | 2018-07-10 | 嘉兴学院 | 一种风车轨桥耦合模型及桥上抗风行车准则制定方法 |
CN108763716A (zh) * | 2018-05-22 | 2018-11-06 | 西南交通大学 | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 |
CN110673743A (zh) * | 2019-08-21 | 2020-01-10 | 佳都新太科技股份有限公司 | 一种基于三维空间场景的虚实交互系统及方法 |
CN111898304A (zh) * | 2020-08-06 | 2020-11-06 | 西南交通大学 | 风车流桥耦合振动分析方法及系统 |
WO2022021934A1 (zh) * | 2020-07-27 | 2022-02-03 | 中铁大桥勘测设计院集团有限公司 | 一种大跨度铁路桥梁行车安全防风设计方法及减风装置 |
CN114117651A (zh) * | 2022-01-26 | 2022-03-01 | 湖南大学 | 风-车-桥全动态耦合分析方法及系统 |
-
2022
- 2022-11-16 CN CN202211434913.5A patent/CN115828379B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107657117A (zh) * | 2017-09-26 | 2018-02-02 | 中交公路长大桥建设国家工程研究中心有限公司 | 一种公铁两用桥梁车‑桥‑风浪流耦合振动分析方法 |
CN108268711A (zh) * | 2018-01-04 | 2018-07-10 | 嘉兴学院 | 一种风车轨桥耦合模型及桥上抗风行车准则制定方法 |
CN108763716A (zh) * | 2018-05-22 | 2018-11-06 | 西南交通大学 | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 |
CN110673743A (zh) * | 2019-08-21 | 2020-01-10 | 佳都新太科技股份有限公司 | 一种基于三维空间场景的虚实交互系统及方法 |
WO2022021934A1 (zh) * | 2020-07-27 | 2022-02-03 | 中铁大桥勘测设计院集团有限公司 | 一种大跨度铁路桥梁行车安全防风设计方法及减风装置 |
CN111898304A (zh) * | 2020-08-06 | 2020-11-06 | 西南交通大学 | 风车流桥耦合振动分析方法及系统 |
CN114117651A (zh) * | 2022-01-26 | 2022-03-01 | 湖南大学 | 风-车-桥全动态耦合分析方法及系统 |
Non-Patent Citations (2)
Title |
---|
大跨度公铁平层桥梁列车-汽车-桥梁动力作用研究;任文渊;中国优秀硕士学位论文全文数据库 工程科技II辑(第第5期期);全文 * |
大跨度公铁平层桥梁列车-汽车-桥梁动力作用研究;任文渊;中国优秀硕士学位论文全文数据库(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115828379A (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106650067B (zh) | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 | |
CN108763716B (zh) | 一种基于振型叠加法的过桥列车倾覆危险性评价方法 | |
Xia et al. | Dynamic analysis of train–bridge system and its application in steel girder reinforcement | |
CN111079323B (zh) | 基于人-车-路耦合振动模型的动力响应预测方法及系统 | |
Guo et al. | Running safety analysis of a train on the Tsing Ma Bridge under turbulent winds | |
Karoumi | Response of cable-stayed and suspension bridges to moving vehicles: Analysis methods and practical modeling techniques | |
CN107273665A (zh) | 磁悬浮车辆动态限界获取方法 | |
Zeng et al. | Hunting stability of high-speed railway vehicles under steady aerodynamic loads | |
CN110717296A (zh) | 高速铁路乘客舒适度评价方法、装置、存储介质及设备 | |
CN115828379B (zh) | 一种风-汽车-列车-桥梁耦合振动分析方法 | |
Zhao et al. | Influence of fluid viscous damper on the dynamic response of suspension bridge under random traffic load | |
Jin et al. | Fragility analysis for vehicle derailment on railway bridges under earthquakes | |
Yin et al. | Three-dimensional vibrations of a suspension bridge under stochastic traffic flows and road roughness | |
CN102367644B (zh) | 钢构高架路主结构 | |
CN113076594B (zh) | 浮车型低地板有轨电车限界计算方法、系统、设备及介质 | |
Gong et al. | Computer simulation of dynamic interactions between vehicle and long span box girder bridges | |
Zhu et al. | Dynamic analysis of suspension-type monorail long-span cable-stayed bridge in a wind-vehicle-bridge system | |
CN113378265B (zh) | 基于多体动力学的桥梁动力性能的评估方法 | |
CN107451305A (zh) | 一种列车‑轨道‑结构耦合系统异步长高效动力分析方法 | |
Kozek et al. | Vibration damping of a flexible car body structure using piezo-stack actuators | |
Xia et al. | Dynamic responses of railway suspension bridges under moving trains | |
Wang et al. | Nonlinear dynamic response analysis of a long-span suspension bridge under running train and turbulent wind | |
Li et al. | Comparison of Vibration and Noise Characteristics of Urban Rail Transit Bridges with Box-Girder and U-Shaped Sections | |
CN113177339B (zh) | 一种地震-风-浪-车-桥耦合振动分析方法 | |
Sun et al. | Simulation Analysis of a Dual-Purpose Intelligent Mobile Platform for Highway and Railway |
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 |