CN115292972B - 地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 - Google Patents
地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 Download PDFInfo
- Publication number
- CN115292972B CN115292972B CN202211224770.5A CN202211224770A CN115292972B CN 115292972 B CN115292972 B CN 115292972B CN 202211224770 A CN202211224770 A CN 202211224770A CN 115292972 B CN115292972 B CN 115292972B
- Authority
- CN
- China
- Prior art keywords
- digital twin
- simulation calculation
- structure subsystem
- track
- digital
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Bridges Or Land Bridges (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,包括以下步骤:建立上部结构子***数字孪生仿真计算模型;建立下部结构子***数字孪生仿真计算模型;建立轨道不平顺数字孪生仿真计算模块;得到在轮轨接触面的耦合动力关系;建立地震波数字孪生仿真计算模块;得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系;建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块;得到高速列车运行数字孪生安全评价指标。本发明能够针对地震作用下高速铁路桥上行车数字孪生精细协同仿真进行数字化构建,解决地震作用下高速铁路桥上行车物理模型与数据驱动模型精准映射建模、设计、分析、诊断、决策的底层关键技术问题。
Description
技术领域
本发明属于交通运输业桥梁工程技术领域,具体涉及一种地震作用下高速铁路桥上行车数字孪生精细协同仿真方法。
背景技术
截至2021年底,全国高铁营运里程已达4万公里,到2025年,以“八纵八横”为主骨架的高速铁路网将覆盖50万人口以上城市,高铁营运里程将达到5万公里,桥梁结构在大多数高铁线路中占比超过50%。另一方面,我国地震具有活动频度高且具有强度大、震源浅、分布广的特点,西部地区断层和地震带广泛分布,地下构造活动较为活跃。这就意味着地震时高速铁路桥梁上行车的概率大大增加,地震时桥梁的振动对行车安全有显著的影响。
近年来地震时高速铁路桥梁上列车运行安全问题和动力响应特性引起了铁路工程界的广泛关注,并围绕地车辆模型和桥梁模型、轮对和桥梁相互作用模型和地震输入方式等方面开展了深入的研究。但多数研究停留于使用商业软件Simpack和Opensees进行联合仿真,无法做到精细化的数字孪生仿真计算。
针对上述问题,很有必要研究出一套自主可控的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,用以***解决准确模拟真实环境中的地震作用下高速铁路桥上行车孪生模型精细协同仿真问题。
发明内容
本发明为解决现有技术存在的问题而提出,其目的是提供一种地震作用下高速铁路桥上行车数字孪生精细协同仿真方法。
本发明的技术方案是:地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,包括以下步骤:
A.建立上部结构子***数字孪生仿真计算模型;
B.建立下部结构子***数字孪生仿真计算模型;
C.建立轨道不平顺数字孪生仿真计算模块;
D.基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系;
E.建立地震波数字孪生仿真计算模块;
F.基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系;
G. 建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块;
H.得到高速列车运行数字孪生安全评价指标。
更进一步的,步骤A建立上部结构子***数字孪生仿真计算模型,具体过程如下:
首先,高速列车根据拉格朗日方程原理建立振动运动方程;
然后,无砟轨道中的扣件、轨枕、道床板采用参振质量的形式进行模拟;
再后,无砟轨道中的钢轨采用带刚臂单元模拟;
最后,桥梁梁体采用带刚臂单元联合变截面空间梁单元模拟。
更进一步的,步骤B建立下部结构子***数字孪生仿真计算模型,具体过程如下:
首先,支座采用支座模型进行模拟;
然后,桥墩采用弹塑性纤维梁柱单元模拟,并对桥墩的横截面进行网格划分。
最后,建立包含弹塑性纤维梁柱单元的下部结构子***动力方程。
桥墩采用弹塑性纤维梁柱单元模拟过程中,还包括以下过程:
首先,桥墩混凝土和桥墩钢筋分别赋予混凝土纤维和钢筋纤维对应的本构模型;
然后,桥墩混凝土纤维采用Kent-Scott-Park混凝土本构模型进行模拟;
最后,桥墩钢筋纤维采用修正梅内戈托与平托模型进行模拟。
更进一步的,步骤C建立轨道不平顺数字孪生仿真计算模块,具体过程如下:
将各类无砟轨道不平顺功率谱密度函数转换成随机轨道不平顺时程曲线,从而使无砟轨道不平顺作为车桥耦合分析***的自激励;
将计算得到的无砟轨道不平顺时程曲线作用于步骤A建立的上部结构子***数字孪生仿真计算模型中,是上部结构子***数字孪生仿真计算模型中高速列车与无砟轨道之间相互作用的激励源。
更进一步的,步骤D基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系,具体过程如下:
首先,以上部结构子***中的高速列车轮和下部结构子***中的无砟轨道之间的接触为界面;
然后,高速列车和无砟轨道在轮轨接触点通过位移协调和轮轨相互作用力相互联系,接触点求解过程中需代入步骤C计算得到的无砟轨道随机不平顺时程曲线;
再后,轮轨相互作用力中的轮轨法向力采用非线性赫兹接触理论,轮轨相互作用力中的轮轨切向蠕滑力采用Kaller理论;
最后,将求解得到的轮轨相互作用力代入步骤A建立的上部结构子***数字孪生仿真计算模型中。
更进一步的,步骤E建立地震波数字孪生仿真计算模块,具体过程如下:
首先,可采用人工合成地震波的方法,建立地震波数字孪生仿真计算模块;
然后,可采用天然地震波调整法,建立地震波数字孪生仿真计算模块。
最后,将地震波数字孪生仿真计算模块计算得到的随机地震波代入步骤B建立的下部结构子***数字孪生仿真计算模型中。
更进一步的,步骤F基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系,具体过程如下:
首先,以步骤A建立的上部结构子***与步骤B建立的下部结构子***之间以支座处为边界;
然后,运用支座处力的平衡和位移协调条件,根据上一时程上部结构子***位移量,计算出上部结构子***作用于边界处的荷载;
再后,步骤B建立的下部结构子***接受上部结构子***传递的荷载;
最后,将传递的荷载代入步骤B建立的动力方程求解对应的动力位移响应,并将该位移响应传递回步骤A建立的上部结构子***作为下一个时程的位移量,如此循环直至结束。
更进一步的,步骤G建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块,具体过程如下:
首先,利用Newton-Raphson迭代算法对步骤B所涉及的桥墩弹塑性纤维梁柱单元进行仿真分析;
然后,分别计算出桥墩在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态;
最后,分别计算出梁体在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态。
更进一步的,步骤H得到高速列车运行数字孪生安全评价指标,具体过程如下:
首先,得到地震情况下步骤A建立的上部结构子***中的车桥耦合动力响应的计算结果;
然后,计算出高速列车运行安全评价指标。
本发明的有益效果如下:
本发明通过分别建立上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块、列车和轨道在轮轨接触面的耦合动力关系、地震波数字孪生仿真计算模块、上部结构子***与下部结构子***之间的数字孪生耦合动力关系、弹塑性纤维梁柱单元数字孪生地震仿真分析模块、高速列车运行数字孪生安全评价指标,进而建立了一套完全自主可控的完整的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法技术。
本发明能够针对交通运输领域的地震作用下高速铁路桥上行车数字孪生精细协同仿真进行数字化构建,解决地震作用下高速铁路桥上行车物理模型与数据驱动模型精准映射建模、设计、分析、诊断、决策的底层关键技术问题。
附图说明
图1为本发明中数字孪生精细协同仿真的模块示意图;
图2为本发明中上部结构子***与下部结构子***的示意图;
图3为本发明中箱梁梁体横截面的示意图;
图4为本发明中箱梁梁体横截面的任意截面输入界面;
图5为本发明中桥墩横截面纤维单元网格划分示意图;
图6为本发明中Kent-Scott-Park混凝土本构模型示意图;
图7为本发明中钢筋纤维的修正梅内戈托与平托模型示意图;
图8为本发明中轮轨接触示意图;
图9为本发明中简支箱梁截面尺寸和Simpack模型图;
图10为本发明中桥墩横截面尺寸及下部结构子***的桥墩实物、支座实物和网格剖分示意图。
具体实施方式
以下,参照附图和实施例对本发明进行详细说明:
如图1至图10所示,地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,包括以下步骤:
A.建立上部结构子***数字孪生仿真计算模型;
B.建立下部结构子***数字孪生仿真计算模型;
C.建立轨道不平顺数字孪生仿真计算模块;
D.基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系;
E.建立地震波数字孪生仿真计算模块;
F.基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系;
G. 建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块;
H.得到高速列车运行数字孪生安全评价指标。
步骤A建立上部结构子***数字孪生仿真计算模型,具体过程如下:
首先,高速列车根据拉格朗日方程原理建立振动运动方程;
然后,无砟轨道中的扣件、轨枕、道床板采用参振质量的形式进行模拟;
再后,无砟轨道中的钢轨采用带刚臂单元模拟;
最后,桥梁梁体采用带刚臂单元联合变截面空间梁单元模拟。
步骤B建立下部结构子***数字孪生仿真计算模型,具体过程如下:
首先,支座采用支座模型进行模拟;
然后,桥墩采用弹塑性纤维梁柱单元模拟,并对桥墩的横截面进行网格划分,如图5。
最后,建立包含弹塑性纤维梁柱单元的下部结构子***动力方程。
桥墩采用弹塑性纤维梁柱单元模拟过程中,还包括以下过程:
首先,桥墩混凝土和桥墩钢筋分别赋予混凝土纤维和钢筋纤维对应的本构模型;
然后,桥墩混凝土纤维采用Kent-Scott-Park混凝土本构模型进行模拟,如图6;
最后,桥墩钢筋纤维采用修正梅内戈托与平托模型进行模拟,如图7。
步骤C建立轨道不平顺数字孪生仿真计算模块,具体过程如下:
将各类无砟轨道不平顺功率谱密度函数转换成随机轨道不平顺时程曲线,从而使无砟轨道不平顺作为车桥耦合分析***的自激励。其中的无砟轨道不平顺功率谱密度函数可选自铁路规范《高速铁路无砟轨道不平顺谱》(TB/T 3352-2014)中所描述的轨距不平顺、水平不平顺、轨向不平顺、高低不平顺等对应的无砟轨道功率谱密度函数,也可来自高铁无砟轨道不平顺实测数据转化的不平顺功率谱密度函数。
将计算得到的无砟轨道不平顺时程曲线作用于步骤A建立的上部结构子***数字孪生仿真计算模型中,是上部结构子***数字孪生仿真计算模型中高速列车与无砟轨道之间相互作用的激励源。
步骤D基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系,具体过程如下:
首先,以上部结构子***中的高速列车轮和下部结构子***中的无砟轨道之间的接触为界面,如图8;
然后,高速列车和无砟轨道在轮轨接触点通过位移协调和轮轨相互作用力相互联系,接触点求解过程中需代入步骤C计算得到的无砟轨道随机不平顺时程曲线。
再后,轮轨相互作用力中的轮轨法向力采用非线性赫兹接触理论,轮轨相互作用力中的轮轨切向蠕滑力采用Kaller理论。
最后,将求解得到的轮轨相互作用力代入步骤A建立的上部结构子***数字孪生仿真计算模型中。
步骤E建立地震波数字孪生仿真计算模块,具体过程如下:
首先,可采用人工合成地震波的方法,建立地震波数字孪生仿真计算模块;
然后,可采用天然地震波调整法,建立地震波数字孪生仿真计算模块。
最后,将地震波数字孪生仿真计算模块计算得到的随机地震波代入步骤B建立的下部结构子***数字孪生仿真计算模型中。
步骤F基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系,具体过程如下:
首先,以步骤A建立的上部结构子***与步骤B建立的下部结构子***之间以支座处为边界;
然后,运用支座处力的平衡和位移协调条件,根据上一时程上部结构子***位移量,计算出上部结构子***作用于边界处的荷载;
再后,步骤B建立的下部结构子***接受上部结构子***传递的荷载;
最后,将传递的荷载代入步骤B建立的动力方程求解对应的动力位移响应,并将该位移响应传递回步骤A建立的上部结构子***作为下一个时程的位移量,如此循环直至结束。
步骤G建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块,具体过程如下:
首先,利用Newton-Raphson迭代算法对步骤B所涉及的桥墩弹塑性纤维梁柱单元进行仿真分析;
然后,分别计算出桥墩在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态;
最后,分别计算出梁体在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态。
步骤H得到高速列车运行数字孪生安全评价指标,具体过程如下:
首先,得到地震情况下步骤A建立的上部结构子***中的车桥耦合动力响应的计算结果;
然后,计算出高速列车运行安全评价指标。
具体的,步骤A建立上部结构子***数字孪生仿真计算模型中,上部结构子***包括高速列车、无砟轨道和桥梁梁,根据孪生对象的实际物理模型尺寸建模。
具体的,步骤B建立下部结构子***数字孪生仿真计算模型中,下部结构子***由桥墩和支座组成,在建立孪生模型时,支座可采用6方向刚度弹性支座模拟,也可采用其他支座模型模拟。
具体的,步骤C建立轨道不平顺数字孪生仿真计算模块中,采用三角级数法来作为不平顺数字孪生仿真计算模块的主算法。
具体的,步骤F基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系中,运用在边界处的位移协调条件,上部结构子***实时接受下部子结构传递的位移实现支座的强迫振动。
具体的,步骤H得到高速列车运行数字孪生安全评价指标,安全评价指标包括脱轨系数、轮重减载率、横向力、车体加速度指标、舒适度指标、平稳性指标。
具体的,步骤A中根据拉格朗日方程建立高速列车模型振动方程的过程如下,以42自由度振动方程为例:
拉格朗日方程的表达式为:
车体的自由度:
转向架的自由度:
轮对的自由度:
基于式(2)~(4)的自由度,分别计算出车辆的总动能、总势能、阻尼耗散能,然后代入式(1)后,即可得到42自由度列车车辆的振动方程。
具体的,钢轨和梁体的截面尺寸可采用任意截面的输入方式进行实际物理模型尺寸建模。
更为具体的,自定义的任意截面输入规则,约定规则如下:
首先,各边按逆时针方向输入,用分号";"隔开;
然后,每条边的信息包含5个参量,中间用逗号","分开;
再后,5个参量的物理意义如下:
①截面号,>0外轮廓,<0内轮廓
②起始端点x坐标
③起始端点z坐标
④半径r,=0直线段,>0圆弧半径
⑤圆弧标示,=0整圆,>0劣弧,<0优弧。
最后,根据上述输入的钢轨和梁体的具体物理模型尺寸,采用带刚臂单元联合变截面空间梁单元模拟。
如附图3所示的箱梁横截面,采用上述任意截面输入规则,具体数据可见表1,基于任意截面输入规则的自主可控数字孪生仿真平台箱梁横截面输入界面如附图4所示。
表1 箱梁横截面输入数据
然后根据上述输入的钢轨和梁体的具体物理模型尺寸,采用带刚臂单元联合变截面空间梁单元模拟。
具体的,步骤B中对于桥墩混凝土的Kent-Scott-Park混凝土本构模型,其应力-应变曲线方程可写为:
具体的,对于钢筋纤维的修正梅内戈托与平托本构模型,其应力-应变曲线方程可写为:
其中:
具体的,步骤C中采用三角级数法作为不平顺数字孪生仿真计算模块的主算法,具体计算公式为:
具体的,步骤D中轮轨法向力采用非线性赫兹接触理论,计算公式如下:
轮轨切向蠕滑力采用Kaller理论,与蠕滑率之间的关系可表示为:
具体的,步骤E中采用人工合成地震波的方法及天然地震波调整法两种方法作为地震波数字孪生仿真计算模块的主要方法,从而形成地震波数字孪生仿真计算模块。
常采用三角级数法人工合成地震波,其计算公式如下:
天然地震波调整法,主要对强度、频率、持时三方面进行调整,天然地震波则采用“两频段双指标法”、“两频段双指标+70%能量持时法”进行选波。
具体的,步骤F中上部结构子***与下部结构子***之间以支座处为边界,上部结构子***实时传递支座处的作用力,运用支座处力的平衡条件和向外插值,下部结构子***接受上部结构子***传递的荷载,下部结构子***可根据动力方程求解响应。可用下列表达式表示:
又一实施例
将利用本发明的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法对某单线高速铁路线路7跨标准跨32m简支箱梁的C50预应力混凝土箱梁地震作用下的行车安全性进行对比研究。
对比研究采用商业软件Simpack与Opensees联合仿真的形式进行。
附图9为简支箱梁的截面尺寸和在商业软件Simpack中建立的高速列车、无砟轨道和桥梁梁体计算模型示意图。
附图10为桥墩横截面尺寸示意图及下部结构子***的桥墩实物、支座实物和网格剖分示意图。
墩高12.5m,桥墩混凝土等级为C35,该部分模型采用商业软件Opensees进行建模仿真。混凝土采用商业软件Opensees中的Concrete02模型,该模型可以考虑箍筋的约束作用对混凝土强度与延性的影响和混凝土的受拉力学性能。钢筋采用Steel02模型,桥梁阻尼比取5%。
采用德国低干扰谱生成轨道不平顺,将轨道高低不平顺和水平不平顺作为***自激励,选用Kern County波作为地震激励。高速列车运行速度分别为300km/h,对地震波进行调幅,研究高速列车在不同地震强度分别为0.035g、0.12g和0.25g的地震作用下的行车安全性。
提取商业有限元联合仿真计算结果,行车安全相关的指标脱轨系数最大值为0.9583(0.25g的地震作用下)及横向轮轨力最大值为52.972KN(动车值)。本发明数字孪生仿真计算结果的脱轨系数最大值为0.9602(0.25g的地震作用下)及横向轮轨力最大值为53.005KN(动车值),最大误差在千分之二以内。
商业有限元联合仿真计算软件的计算结果验证了本发明所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法的准确性。
本发明通过分别建立上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块、列车和轨道在轮轨接触面的耦合动力关系、地震波数字孪生仿真计算模块、上部结构子***与下部结构子***之间的数字孪生耦合动力关系、弹塑性纤维梁柱单元数字孪生地震仿真分析模块、高速列车运行数字孪生安全评价指标,进而建立了一套完全自主可控的完整的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法技术。
本发明能够针对交通运输领域的地震作用下高速铁路桥上行车数字孪生精细协同仿真进行数字化构建,解决地震作用下高速铁路桥上行车物理模型与数据驱动模型精准映射建模、设计、分析、诊断、决策的底层关键技术问题。
Claims (8)
1.地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:包括以下步骤:
(A)建立上部结构子***数字孪生仿真计算模型;
(B)建立下部结构子***数字孪生仿真计算模型;
(C)建立轨道不平顺数字孪生仿真计算模块;
(D)基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系;
(E)建立地震波数字孪生仿真计算模块;
(F)基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系;
(G)针对下部结构子***数字孪生仿真计算模型,建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块;
(H)得到高速列车运行数字孪生安全评价指标;
步骤(A)建立上部结构子***数字孪生仿真计算模型,具体过程如下:
首先,高速列车根据拉格朗日方程原理建立振动运动方程;
然后,无砟轨道中的扣件、轨枕、道床板采用参振质量的形式进行模拟;
再后,无砟轨道中的钢轨采用带刚臂单元模拟;
最后,桥梁梁体采用基于任意截面输入格式的带刚臂单元联合变截面空间梁单元模拟;
步骤(G)建立弹塑性纤维梁柱单元数字孪生地震仿真分析模块,具体过程如下:
首先,利用Newton-Raphson迭代算法对步骤(B)所涉及的桥墩弹塑性纤维梁柱单元进行仿真分析;
然后,分别计算出桥墩在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态;
最后,分别计算出梁体在指定地震波情况下的动力响应及对应的各部位钢筋、混凝土应力应变指标,判别是否进入延性状态;
自定义的任意截面输入规则,约定规则如下:
首先,各边按逆时针方向输入,用分号";"隔开;
然后,每条边的信息包含5个参量,中间用逗号","分开;
再后,5个参量的物理意义如下:
①截面号,>0外轮廓,<0内轮廓
②起始端点x坐标
③起始端点z坐标
④半径r,=0直线段,>0圆弧半径
⑤圆弧标示,=0整圆,>0劣弧,<0优弧;
步骤C中采用三角级数法作为不平顺数字孪生仿真计算模块的主算法,具体计算公式为:
对于钢筋纤维的修正梅内戈托与平托本构模型,其应力-应变曲线方程可写为:
为刚度折减率;
常采用三角级数法人工合成地震波,其计算公式如下:
天然地震波调整法,主要对强度、频率和持时三方面进行调整,天然地震波则采用“两频段双指标法”、“两频段双指标+70%能量持时法”中的一种进行选波。
2.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(B)建立下部结构子***数字孪生仿真计算模型,具体过程如下:
首先,支座采用支座模型进行模拟;
然后,桥墩采用弹塑性纤维梁柱单元模拟,并对桥墩的横截面进行网格划分;
最后,建立包含弹塑性梁柱单元的下部结构子***动力方程。
3.根据权利要求2所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:桥墩采用弹塑性纤维梁柱单元模拟过程中,还包括以下过程:
首先,桥墩混凝土和桥墩钢筋分别赋予混凝土纤维和钢筋纤维对应的本构模型;
然后,桥墩混凝土纤维采用Kent-Scott-Park混凝土本构模型进行模拟;
最后,桥墩钢筋纤维采用修正梅内戈托与平托模型进行模拟。
4.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(C)建立轨道不平顺数字孪生仿真计算模块,具体过程如下:
将各类无砟轨道不平顺功率谱密度函数转换成随机轨道不平顺时程曲线,从而使无砟轨道不平顺作为车桥耦合分析***的自激励;
将计算得到的无砟轨道不平顺时程曲线作用于步骤(A)建立的上部结构子***数字孪生仿真计算模型中,是上部结构子***数字孪生仿真计算模型中高速列车与无砟轨道之间相互作用的激励源。
5.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(D)基于上部结构子***数字孪生仿真计算模型、轨道不平顺数字孪生仿真计算模块,得到列车和轨道在轮轨接触面的耦合动力关系,具体过程如下:
首先,以上部结构子***中的高速列车轮和下部结构子***中的无砟轨道之间的接触为界面;
然后,高速列车和无砟轨道在轮轨接触点通过位移协调和轮轨相互作用力相互联系,接触点求解过程中需代入步骤(C)计算得到的无砟轨道随机不平顺时程曲线;
再后,轮轨相互作用力中的轮轨法向力采用非线性赫兹接触理论,轮轨相互作用力中的轮轨切向蠕滑力采用Kaller理论;
最后,将求解得到的轮轨相互作用力代入步骤(A)建立的上部结构子***数字孪生仿真计算模型中。
6.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(E)建立地震波数字孪生仿真计算模块,具体过程如下:
首先,可采用人工合成地震波的方法,建立地震波数字孪生仿真计算模块;
然后,可采用天然地震波调整法,建立建立地震波数字孪生仿真计算模块;
最后,将地震波数字孪生仿真计算模块计算得到的随机地震波代入步骤(B)建立的下部结构子***数字孪生仿真计算模型中。
7.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(F)基于上部结构子***数字孪生仿真计算模型、下部结构子***数字孪生仿真计算模型、得到上部结构子***与下部结构子***之间的数字孪生耦合动力关系,具体过程如下:
首先,以步骤(A)建立的上部结构子***与步骤(B)建立的下部结构子***之间以支座处为边界;
然后,运用支座处力的平衡和位移协调条件,根据上一时程上部结构子***位移量,计算出上部结构子***作用于边界处的荷载;
再后,步骤(B)建立的下部结构子***接受上部结构子***传递的荷载;
最后,将传递的荷载代入步骤(B)建立的动力方程求解对应的动力位移响应,并将该位移响应传递回步骤(A)建立的上部结构子***作为下一个时程的位移量,如此循环直至结束。
8.根据权利要求1所述的地震作用下高速铁路桥上行车数字孪生精细协同仿真方法,其特征在于:步骤(H)得到高速列车运行数字孪生安全评价指标,具体过程如下:
首先,得到地震情况下步骤(A)建立的上部结构子***中的车桥耦合动力响应的计算结果;
然后,计算出高速列车运行安全评价指标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211224770.5A CN115292972B (zh) | 2022-10-09 | 2022-10-09 | 地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211224770.5A CN115292972B (zh) | 2022-10-09 | 2022-10-09 | 地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115292972A CN115292972A (zh) | 2022-11-04 |
CN115292972B true CN115292972B (zh) | 2023-01-24 |
Family
ID=83834178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211224770.5A Active CN115292972B (zh) | 2022-10-09 | 2022-10-09 | 地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115292972B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117610306B (zh) * | 2023-12-15 | 2024-06-18 | 大连海事大学 | 一种移动荷载作用下简支梁桥数字孪生构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114638046A (zh) * | 2022-05-12 | 2022-06-17 | 中国铁路设计集团有限公司 | 一种铁路桥墩数字孪生变截面仿真计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110147623B (zh) * | 2019-05-23 | 2023-04-18 | 南京金蓝智慧城市规划设计有限公司 | 高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法 |
CN111338300B (zh) * | 2020-02-27 | 2020-11-06 | 广东工业大学 | 生产线基于数字孪生的物理仿真方法及其*** |
-
2022
- 2022-10-09 CN CN202211224770.5A patent/CN115292972B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114638046A (zh) * | 2022-05-12 | 2022-06-17 | 中国铁路设计集团有限公司 | 一种铁路桥墩数字孪生变截面仿真计算方法 |
Non-Patent Citations (1)
Title |
---|
强地震下高速铁路桥上行车精细化模拟及行车安全性分析;国巍 等;《工程力学》;20180630;正文第1-4章 * |
Also Published As
Publication number | Publication date |
---|---|
CN115292972A (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lei et al. | Dynamic response analyses of vehicle and track coupled system on track transition of conventional high speed railway | |
Lei et al. | Analyses of dynamic response of vehicle and track coupling system with random irregularity of track vertical profile | |
CN106650067B (zh) | 一种输出中高频响应的车-轨-桥耦合振动频域分析方法 | |
Peixer et al. | Running safety evaluation of a train moving over a high-speed railway viaduct under different track conditions | |
Yang et al. | An iterative interacting method for dynamic analysis of the maglev train–guideway/foundation–soil system | |
CN115292972B (zh) | 地震作用下高速铁路桥上行车数字孪生精细协同仿真方法 | |
Cao et al. | A semi-analytical/FEM model for predicting ground vibrations induced by high-speed train through continuous girder bridge | |
Guo et al. | Running safety analysis of a train on the Tsing Ma Bridge under turbulent winds | |
CN106951668A (zh) | 一种铁路大跨度斜拉桥的刚度控制方法及斜拉桥 | |
Karoumi | Response of cable-stayed and suspension bridges to moving vehicles: Analysis methods and practical modeling techniques | |
Zhang et al. | Train-induced vibration and noise radiation of a prestressed concrete box-girder | |
Sun et al. | Prediction and mitigation analysis of ground vibration caused by running high-speed trains on rigid-frame viaducts | |
Jiang et al. | Running safety and seismic optimization of a fault-crossing simply-supported girder bridge for high-speed railways based on a train-track-bridge coupling system | |
CN112364426A (zh) | 一种基于行车安全及动力响应的铁路桥墩伤损评定方法、***、终端设备及可读存储介质 | |
Fang et al. | Dynamic behavior of railway Vehicle-Ballasted track system with unsupported sleepers based on the hybrid DEM-MBD method | |
Zhang et al. | Effects of fundamental factors on coupled vibration of wind-rail vehicle-bridge system for long-span cable-stayed bridge | |
Xia et al. | Application of train-bridge-interaction analysis to bridge design of high-speed railways in China | |
Li et al. | Influences of soil-structure interaction on coupled vibration of train-bridge system: theoretical and experimental study | |
Yi | Dynamic analysis of high-speed railway alignment: theory and practice | |
Correia et al. | Dynamic analysis of rail track for high speed trains. 2D approach | |
CN115795628A (zh) | 基于联合仿真的盾构下穿高铁跨线桥动力响应建模分析方法 | |
Ahlbeck | Effects of track dynamic impedance on vehicle-track interactions | |
Tang et al. | Dynamic responses of train under deformation caused by underpass shield tunnel construction based on stochastic field theory | |
Hoorpah | Dynamic calculations of high-speed railway bridges in France–some case studies | |
Liu et al. | Retracted: Running Safety Analysis of Long-Span Railway Cable-Stayed Bridge Under Seismic Excitations |
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 |