CN108021747A - 一种高铁弓网动态行为动网格不平衡力消除的仿真方法 - Google Patents

一种高铁弓网动态行为动网格不平衡力消除的仿真方法 Download PDF

Info

Publication number
CN108021747A
CN108021747A CN201711239896.9A CN201711239896A CN108021747A CN 108021747 A CN108021747 A CN 108021747A CN 201711239896 A CN201711239896 A CN 201711239896A CN 108021747 A CN108021747 A CN 108021747A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
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.)
Granted
Application number
CN201711239896.9A
Other languages
English (en)
Other versions
CN108021747B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201711239896.9A priority Critical patent/CN108021747B/zh
Publication of CN108021747A publication Critical patent/CN108021747A/zh
Application granted granted Critical
Publication of CN108021747B publication Critical patent/CN108021747B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Current-Collector Devices For Electrically Propelled Vehicles (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，包括以下步骤：步骤1：根据有限元和多体动力学方法构建弓网非线性仿真模型；步骤2：每个仿真步时间t中，采用动网格方法加密接触线上接触点周围的网格；步骤3：通过迭代算法消除动网格中的不平衡力；步骤4：通过Newmark‑β算法求解仿真步时间t时的响应；步骤5：重复步骤2～4，直至受电弓遍历整个接触线完成弓网动力学仿真；本发明可消除采用动网格技术产生的不平衡力，用于非线性有限元求解中，提高求解效率和求解精度。

Description

一种高铁弓网动态行为动网格不平衡力消除的仿真方法
技术领域
本发明涉及高铁弓网动态行为的仿真方法，具体涉及一种高铁弓网动态行为动网格不平衡力消除的仿真方法。
背景技术
高速铁路受电弓和接触网系统负责为高速运行的电力机车输送电能，弓网之间通过滑动接触实现电能的持续传输；为了研究弓网之前复杂的动力学行为，通常采用有限元方法构建接触网的数学模型，通过迭代求解的方法的到弓网交互作用下的动态响应；然而，由于接触网具有大跨距、长距离等特点，弓网的动态仿真的计算量是十分可观对的；为了提升弓网的计算效率，Jimenez-Octavio、Carnicero、Sanchez-Rebollo和Such等人提出了弓网系统的动网格仿真技术，该成果《A moving mesh method to deal with cablestructures subjected to moving loads and its application to the catenary-pantograph dynamic interaction》发表于期刊《Journal of Sound and Vibration》的第349卷216-229页；该方法的核心思想是在弓网动态交互仿真过程中，局部加密接触点附近的接触线网格，起到提高计算精度、提升计算效率的目的；然而，这一方法目前只能用于线性有限元模型中，无法考虑接触线大变形产生的几何非线性；当考虑接触线的几何非线性时，新生成的动网格会产生较大的不平衡力，从而引起较大的误差，甚至无法求解。
发明内容
本发明提供一种可消除采用动网格技术产生的不平衡力的仿真方法，用于非线性有限元求解中，提高求解效率和求解精度的一种高铁弓网动态行为动网格不平衡力消除的仿真方法。
本发明采用的技术方案是：一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，包括以下步骤：
步骤1：根据有限元和多体动力学方法构建弓网非线性仿真模型；
步骤2：每个仿真步时间t中，采用动网格方法加密接触线上接触点周围的网格；
步骤3：通过迭代算法消除动网格中的不平衡力；
步骤4：通过纽马克贝塔算法求解仿真步时间t时的响应；
步骤5：重复步骤2～4，直至受电弓遍历整个接触线完成弓网动力学仿真。
进一步的，所述步骤3通过迭代算法消除动网格中的不平衡力的具体过程如下：
S1：提取动网格区域固定网格节点的位移响应；
S2：根据网格节点的位移初始化动网格节点的位移；
S3：根据弓网非线性仿真模型得到广义单元刚度矩阵和不平衡力矩阵；
S4：根据步骤S1和步骤S2中得到的动、定网格节点位移得到广义局部刚度矩阵和不平衡力向量；
S5：对步骤S3和步骤S4中得到的矩阵和向量施加局部约束条件；
S6：计算动网格位移增量，更新动网格位移响应；
S7：检验收敛性，若收敛则转入步骤4，若不收敛则转入S2。
进一步的，所述步骤1具体过程如下：
S21：采用非线性索单元离散接触网的接触线和承力索，力学平衡表达式：
式中：F1、F2、F3、F4、F5和F6均为两个端点x、y、z轴方向上的端点力；lx、ly和lz分别为索单元两个端点在x轴方向、y轴方向和z轴方向上的相对距离；L0为索单元初始长度；w为单位长度的索单元的自重；E为索单元杨氏模量，A为索单元的截面积；
S22：采用非线性杆单元模拟接触网吊弦的非线性伸缩行为，其力学平衡表达式为：
式中：Fg1、Fg2、Fg3为x、y、z轴方向上的端点力，lgx、lgy和lgz分半为两端点在x轴方向、y轴方向和z轴方向上的相对距离；Lg0为杆单元的初始长度，Eg为杆单元杨氏模量，Ag为杆单元的截面积；
S23：根据接触网系统的拓扑结构，通过有限元方法形成接触网的刚度矩阵Kc和集中质量矩阵Mc；
S24：根据步骤S21～S23得到的式(1)～式(6)和接触网的刚度矩阵Kc和集中质量矩阵Mc构建接触网的力学模型；
S25：基于多体动力学方法构建受电弓模型；
根据多体动力学方法，生成受电弓的集中质量模型：
式中：m1、m2、m3分别为弓头、上框架和下框架的质量；c1、c2、c3分别为弓头、上框架和下框架的阻尼；k1、k2、k3分别为弓头、上框架和下框架的质量；z1、z2、z3分别为弓头、上框架和下框架的位移；f0和fc分别为静态抬升力和接触力；
S26：通过罚函数模拟弓网耦合接触；
通过罚函数模拟弓网之间的耦合交互作用，表达式如下：
fc＝ks(z1-zc)(8)
式中：，zc为接触线接触点处的抬升量，ks为接触刚度；
对接触网和受电弓系统耦合，形成弓网有限元动力学方程：
式中：M、C和K分别为整体质量、阻尼和刚度矩阵，U为整体位移向量，F为外载荷向量；和分别为加速度和速度向量。
进一步的，所述步骤2的具体过程如下：
根据仿真步时间t和车速v更新受电弓在接触线上的接触点：
xc＝vt (10)
式中：xc为受电弓在接触线上的接触点x轴坐标；
定义最小网格长度lg，变化比率为α，动网格单侧区间长度为Lc；动网格节点横坐标为：
xm＝[xc-lgαn,xc-lgαn-1…xc-lgα1,xc-lg,xc,xc+lg,xc…+lgα1…xc+lgαn-1,xc+lgαn]T(11)
式中：xm为动网格节点横坐标向量，n为单侧动网格节点数，满足下述关系：
lgαn≤Lc＜lgαn+1 (12)。
进一步的，所述步骤3具体过程如下：
S31：提取动网格在xc-Lc到xc+Lc区间内x、y、z轴方向上原有定网格位移向量xs、ys、zs；
S32：根据xs、ys、zs以及动网格的位置向量xm，基于线性插值法生成动网格的初始位移向量ym和zm；
S33：根据弓网非线性仿真模型得到广义单元刚度矩阵和不平衡力矩阵；
分别对式(1)、(2)和(3)两端积分，得到：
式中：和为相应的单元柔度矩阵，dL0为初始位移增量，dF为单元不平衡力增量。
柔度矩阵求逆得到：
其中：和为响应的刚度矩阵，dXe为单元节点位移增量向量，dFe为端点力增量向量。
S34：根据动、定网格节点位移生成广义局部刚度矩阵和不平衡力向量；
根据动网格节点位移向量xm、ym、zm和定网格节点位移向量xs、ys、zs，生成每个网格的单元刚度矩阵dFe、和
根据有限元方法生成整体刚度矩阵和
整体静力学方程为：
dF＝KTdX+KGdL0 (17)
式中：dX为动网格增量,dL0为初始长度增量；
S34：施加局部约束条件；
在所有网格节点x方向施加虚约束，即：消除矩阵KT和向量dX对应的列向量；
保证每个单元的总初始长度之和不变，即施加约束
在定网格点x、y、z方向上施加固定约束，即：消除矩阵KT和向量dX对应的行、列向量；
S35：计算动网格位移增量，更新动网格位移相应；
根据式(17)计算动网格增量和初始长度增量；更新动网格位移向量：xm、ym、zm；
S36：判断是否满足max[dX dL0]＜d，如果满足则转入步骤4，如果不满足则返回步骤S32。
进一步的，所述步骤4具体过程如下：
将动网格节点位移向量xm、ym、zm引入到弓网整体位移向量U中，重新生成整体质量矩阵M，整体刚度矩阵K和整体阻尼矩阵C，通过Newmark-β法求解该时刻的位移响应，并进入下一时刻t+△t进行新的迭代计算。
本发明的有益效果是：
(1)本发明通过对仿真步中划分的新的移动网格进行迭代计算可消除新网格的不平衡力；
(2)本发明的仿真技术用于弓网动力学的非线性有限元仿真中，充分保证接触线的几何非线性；
(3)本发明的仿真技术用于弓网动力学的非线性有限元仿真中，可提高其求解效率和求解精度。
附图说明
图1为本发明中接触线上的动网格划分示意图。
图2为本发明中动网格各节点上的不平衡力。
图3为本发明中经迭代算法消除后的动网格节点不平衡力。
图4为本发明动网格接触力计算结果。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
如图1-3所示，一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，包括以下步骤：
步骤1：根据有限元和多体动力学方法构建弓网非线性仿真模型；
其具体过程如下：
S21：采用非线性索单元离散接触网的接触线和承力索，力学平衡表达式：
式中：F1、F2、F3、F4、F5和F6均为两个端点x、y、z轴方向上的端点力；lx、ly和lz分别为索单元两个端点在x轴方向、y轴方向和z轴方向上的相对距离；L0为索单元初始长度；w为单位长度的索单元的自重；E为索单元杨氏模量，A为索单元的截面积；
S22：采用非线性杆单元模拟接触网吊弦的非线性伸缩行为，其力学平衡表达式为：
式中：Fg1、Fg2、Fg3为x、y、z轴方向上的端点力，lgx、lgy和lgz分半为两端点在x轴方向、y轴方向和z轴方向上的相对距离；Lg0为杆单元的初始长度，Eg为杆单元杨氏模量，Ag为杆单元的截面积；
S23：根据接触网系统的拓扑结构，通过有限元方法形成接触网的刚度矩阵Kc和集中质量矩阵Mc；
S24：根据步骤S21～S23得到的式(1)～式(6)和接触网的刚度矩阵Kc和集中质量矩阵Mc构建接触网的力学模型；
S25：基于多体动力学方法构建受电弓模型；
根据多体动力学方法，生成受电弓的集中质量模型：
式中：m1、m2、m3分别为弓头、上框架和下框架的质量；c1、c2、c3分别为弓头、上框架和下框架的阻尼；k1、k2、k3分别为弓头、上框架和下框架的质量；z1、z2、z3分别为弓头、上框架和下框架的位移；f0和fc分别为静态抬升力和接触力；
S26：通过罚函数模拟弓网耦合接触；
通过罚函数模拟弓网之间的耦合交互作用，表达式如下：
fc＝ks(z1-zc) (8)
式中：，zc为接触线接触点处的抬升量，ks为接触刚度；
通过式(8)对接触网和受电弓系统耦合，形成弓网有限元动力学方程：
式中：M、C和K分别为整体质量、阻尼和刚度矩阵，U为整体位移向量，F为外载荷向量；和分别为加速度和速度向量。
步骤2：每个仿真步时间t中，采用动网格方法加密接触线上接触点周围的网格；
根据仿真步时间t和车速v更新受电弓在接触线上的接触点：
xc＝vt (10)
式中：xc为受电弓在接触线上的接触点x轴坐标；
定义最小网格长度lg，变化比率为a，动网格单侧区间长度为Lc；动网格节点横坐标为：
xm＝[xc-lgαn,xc-lgαn-1…xc-lgα1,xc-lg,xc,xc+lg,xc+lgα1…xc+lgαn-1,xc+lgαn]T(11)
式中xm为动网格节点横坐标向量，n为单侧动网格节点数，α1～αn为，式中xc是否为xc；
满足下述关系：
lgαn≤Lc＜lgαn+1 (12)。
步骤3：通过迭代算法消除动网格中的不平衡力；
具体过程如下：
S1：提取动网格区域固定网格节点的位移响应；
根据前仿真步(t-△t)的计算结果，提取动网格在xc-Lc到xc+Lc区间内x、y、z轴方向上原有定网格位移向量xs、ys、zs。
S2：根据网格节点的位移初始化动网格节点的位移；
根据xs、ys、zs以及动网格的位置向量xm，基于线性插值法生成动网格的初始位移向量ym和zm。
S3：根据弓网非线性仿真模型得到广义单元刚度矩阵和不平衡力矩阵；
分别对式(1)、(2)和(3)两端积分，得到：
式中：和为相应的单元柔度矩阵，dL0为初始位移增量，dF为单元不平衡力增量。
柔度矩阵求逆得到：
其中：和为响应的刚度矩阵，dXe为单元节点位移增量向量，dFe为端点力增量向量；
S4：根据步骤S1和步骤S2中得到的动、定网格节点位移得到广义局部刚度矩阵和不平衡力向量；
根据动网格节点位移向量xm、ym、zm和定网格节点位移向量xs、ys、zs，生成每个网格的单元刚度矩阵dFe、和
根据有限元方法生成整体刚度矩阵和
整体静力学方程为：
dF＝KTdX+KGdL0 (17)
式中：dX为动网格增量,dL0为初始长度增量。
S5：对步骤S3和步骤S4中得到的矩阵和向量施加局部约束条件；
1)在所有网格节点x方向施加虚约束，即：消除矩阵KT和向量dX对应的列向量；
2)保证每个单元的总初始长度之和不变，即施加约束
3)在定网格点x、y、z方向上施加固定约束，即：消除矩阵KT和向量dX对应的行、列向量。
S6：计算动网格位移增量，更新动网格位移响应；
根据式(17)计算动网格增量和初始长度增量；更新动网格位移向量：xm、ym、zm。
S7：检验收敛性
判断是否满足max[dX dL0]＜d，如果满足则转入步骤4，如果不满足则返回步骤S32。
步骤4：通过Newmark-β算法求解仿真步时间t时的响应；
将动网格节点位移向量xm、ym、zm引入到弓网整体位移向量U中，重新生成整体质量矩阵M，整体刚度矩阵K和整体阻尼矩阵C，通过纽马克贝塔Newmark-β法求解该时刻的位移响应，并进入下一时刻t+△t进行新的迭代计算。
步骤5：重复步骤2～4，直至受电弓遍历整个接触线完成弓网动力学仿真。
下面以仿真列车运行速度为350km/h的京津线弓网系统动力学行为为例，按照步骤1的方法构建弓网动力学模型；具体模型构建方法可参考Yang Song、Zhigang Liu、Hongrui Wang和Xiaobing Lu等发表在《Vehicle System Dynamics》上第53卷1455-1479页的成果《Nonlinear modelling of high-speed catenary based on analyticalexpressions of cable and truss elements》；在步骤2中当仿真时间步t＝3.3s时，根据式(10)计算得弓网接触点在320.83m处；本实例中取Lc＝6.5m、lg＝0.2m、α＝2，从t-△t时刻到t时刻，接触线上的动网格划分示意图如图1所示；在步骤2中，根据式(11)和式(12)，形成的动网格横坐标向量xm为：
xm＝[-6.2 -3 -1.4 -0.6 -0.2 0 0.2 0.6 1.4 3 6.2]
随后，根据前仿真步(t-△t)的计算结果，提取在动网格区间内原有定网格的位移响应xs、ys、zs；其x方向上的定网格横坐标向量xs为：
xs＝[-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7]
通过线性插值法确定动网格在y和z方向的位移向量ym和zm；此时，动网格的节点不平衡力如图2所示，可以看出，节点1和11处的不平衡力已经超过15N，继续求解会引起较大误差，甚至会导致不收敛。
然后进入步骤3，在动网格区域内，同时考虑定网格向量xs、ys、zs和动网格向量xm、ym、zm，根据式(13)～式(16)生成该区间内每个网格的单元刚度矩阵dFe、和根据有限元方法生成整体刚度矩阵和整体静力学方程可以写成：
其中，[KT KG]为69×91的矩阵，不是方阵，无法直接求解；采用步骤3中的1)和2)虚约束条件，使得该刚度矩阵转化为69×69的方阵；再施加3)中的固定约束条件，进一步减少其自由度，使该刚度矩阵减少为24×24的方阵，从而能够进行求解。
本实施例中设定阈值d＝10-5，通过迭代计算，当步骤4中的收敛条件满足时，动网格节点的不平衡力如图3所示；可以看出，最大的节点不平衡力仅为3×10-4N，通过该算法可以成功地消除动网格节点的不平衡力。
通过步骤5中的迭代求解计算，可以最终输出弓网间的接触压力曲线；图4给出了本发明的计算结果和传统定网格计算结果的比较；可以看出，无论是接触力的波形、波动范围，本发明的计算结果和传统定网格计算结果具有高度的一致性；接触压力统计量如表1所示。
表1.接触压力统计量表
从表1可以看出，本发明的接触力计算结果和传统定网格计算结果相差不大，尤其是接触力标准差、均值等较为重要的统计量具有高度的一致性，证明了本发明的有效性；此外，本发明的计算时间仅需要2200s，相比于传统的定网格计算方法节省了超过40％的时间；体现了本发明的高效性。
本发明仿真方法能够在保证几何非线性的同时有效消除新网格的不平衡力，并且将本发明仿真方法扩展到非线性有限元求解中，可提升求解效率和求解精度。

Claims (6)

1.一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，其特征在于，包括以下步骤：
步骤1：根据有限元和多体动力学方法构建弓网非线性仿真模型；
步骤2：每个仿真步时间t中，采用动网格方法加密接触线上接触点周围的网格；
步骤3：通过迭代算法消除动网格中的不平衡力；
步骤4：通过纽马克贝塔算法求解仿真步时间t时的响应；
步骤5：重复步骤2～4，直至受电弓遍历整个接触线完成弓网动力学仿真。
2.根据权利要求1所述的一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，其特征在于，所述步骤3通过迭代算法消除动网格中的不平衡力的具体过程如下：
S1：提取动网格区域固定网格节点的位移响应；
S2：根据网格节点的位移初始化动网格节点的位移；
S3：根据弓网非线性仿真模型得到广义单元刚度矩阵和不平衡力矩阵；
S4：根据步骤S1和步骤S2中得到的动、定网格节点位移得到广义局部刚度矩阵和不平衡力向量；
S5：对步骤S3和步骤S4中得到的矩阵和向量施加局部约束条件；
S6：计算动网格位移增量，更新动网格位移响应；
S7：检验收敛性，若收敛则转入步骤4，若不收敛则转入S2。
3.根据权利要求1所述的一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，其特征在于，所述步骤1具体过程如下：
S21：采用非线性索单元离散接触网的接触线和承力索，力学平衡表达式：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>l</mi> <mi>x</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>F</mi> <mn>1</mn> </msub> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>A</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <msub> <mi>F</mi> <mn>1</mn> </msub> <mi>w</mi> </mfrac> <mo>{</mo> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>4</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>5</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>6</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>+</mo> <msub> <mi>F</mi> <mn>6</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>l</mi> <mi>y</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>F</mi> <mn>2</mn> </msub> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>A</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>F</mi> <mn>2</mn> </msub> <mi>w</mi> </mfrac> <mo>{</mo> <mi>l</mi> <mi>n</mi> <mo>(</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>wL</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>wL</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> <mo>&amp;rsqb;</mo> <mo>-</mo> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>l</mi> <mi>z</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>F</mi> <mn>3</mn> </msub> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>A</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msubsup> <mi>wL</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>A</mi> </mrow> </mfrac> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mi>w</mi> </mfrac> <mo>&amp;lsqb;</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>4</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>5</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>6</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中：F1、F2、F3、F4、F5和F6均为两个端点x、y、z轴方向上的端点力；lx、ly和lz分别为索单元两个端点在x轴方向、y轴方向和z轴方向上的相对距离；L0为索单元初始长度；w为单位长度的索单元的自重；E为索单元杨氏模量，A为索单元的截面积；
S22：采用非线性杆单元模拟接触网吊弦的非线性伸缩行为，其力学平衡表达式为：
<mrow> <msub> <mi>F</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>g</mi> </msub> <msub> <mi>A</mi> <mi>g</mi> </msub> </mrow> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> </mfrac> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>F</mi> <mrow> <mi>g</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>g</mi> </msub> <msub> <mi>A</mi> <mi>g</mi> </msub> </mrow> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> </mfrac> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>F</mi> <mrow> <mi>g</mi> <mn>3</mn> </mrow> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>g</mi> </msub> <msub> <mi>A</mi> <mi>g</mi> </msub> </mrow> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> </mfrac> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mi>z</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>L</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式中：Fg1、Fg2、Fg3为x、y、z轴方向上的端点力，lgx、lgy和lgz分半为两端点在x轴方向、y轴方向和z轴方向上的相对距离；Lg0为杆单元的初始长度，Eg为杆单元杨氏模量，Ag为杆单元的截面积；
S23：根据接触网系统的拓扑结构，通过有限元方法形成接触网的刚度矩阵Kc和集中质量矩阵Mc；
S24：根据步骤S21～S23得到的式(1)～式(6)和接触网的刚度矩阵Kc和集中质量矩阵Mc构建接触网的力学模型；
S25：基于多体动力学方法构建受电弓模型；
根据多体动力学方法，生成受电弓的集中质量模型：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>m</mi> <mn>1</mn> </msub> <mfrac> <mrow> <msup> <mi>d</mi> <mn>2</mn> </msup> <msub> <mi>z</mi> <mn>1</mn> </msub> </mrow> <mrow> <msup> <mi>dt</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>f</mi> <mi>c</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>m</mi> <mn>2</mn> </msub> <mfrac> <mrow> <msup> <mi>d</mi> <mn>2</mn> </msup> <msub> <mi>z</mi> <mn>2</mn> </msub> </mrow> <mrow> <msup> <mi>dt</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>3</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>m</mi> <mn>3</mn> </msub> <mfrac> <mrow> <msup> <mi>d</mi> <mn>2</mn> </msup> <msub> <mi>z</mi> <mn>3</mn> </msub> </mrow> <mrow> <msup> <mi>dt</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>3</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <mfrac> <mrow> <msub> <mi>dz</mi> <mn>3</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>k</mi> <mn>3</mn> </msub> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>=</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中：m1、m2、m3分别为弓头、上框架和下框架的质量；c1、c2、c3分别为弓头、上框架和下框架的阻尼；k1、k2、k3分别为弓头、上框架和下框架的质量；z1、z2、z3分别为弓头、上框架和下框架的位移；f0和fc分别为静态抬升力和接触力；
S26：通过罚函数模拟弓网耦合接触；
通过罚函数模拟弓网之间的耦合交互作用，表达式如下：
fc＝ks(z1-zc) (8)
式中：，zc为接触线接触点处的抬升量，ks为接触刚度；
对接触网和受电弓系统耦合，形成弓网有限元动力学方程：
<mrow> <mi>M</mi> <mover> <mi>U</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>C</mi> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>K</mi> <mi>U</mi> <mo>=</mo> <mi>F</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式中：M、C和K分别为整体质量、阻尼和刚度矩阵，U为整体位移向量，F为外载荷向量；和分别为加速度和速度向量。
4.根据权利要求3所述的一种用于高铁弓网动态行为动网格不平衡力消除的仿真方法，其特征在于，所述步骤2的具体过程如下：
根据仿真步时间t和车速v更新受电弓在接触线上的接触点：
xc＝vt (10)
式中：xc为受电弓在接触线上的接触点x轴坐标；
定义最小网格长度lg，变化比率为α，动网格单侧区间长度为Lc；动网格节点横坐标为：
xm＝[xc-lgαn,xc-lgαn-1…xc-lgα1,xc-lg,xc,xc+lg,xc+lgα1…xc+lgαn-1,xc+lgαn]T (11)
式中：xm为动网格节点横坐标向量，n为单侧动网格节点数，满足下述关系：
lgαn≤Lc＜lgαn+1 (12)。
5.根据权利要求4所述的一种用于高铁弓网动态行为动网格仿真的不平衡力消除方法，其特征在于，所述步骤3具体过程如下：
S31：提取动网格在xc-Lc到xc+Lc区间内x、y、z轴方向上原有定网格位移向量xs、ys、zs；
S32：根据xs、ys、zs以及动网格的位置向量xm，基于线性插值法生成动网格的初始位移向量ym和zm；
S33：根据弓网非线性仿真模型得到广义单元刚度矩阵和不平衡力矩阵；
分别对式(1)、(2)和(3)两端积分，得到：
<mrow> <mtable> <mtr> <mtd> <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>dl</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dl</mi> <mi>y</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dl</mi> <mi>z</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>x</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>x</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>x</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>y</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>y</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>y</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>z</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>z</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>z</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>F</mi> <mn>3</mn> </msub> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>x</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>y</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>l</mi> <mi>z</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>L</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>dL</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msubsup> <mi>F</mi> <mi>C</mi> <mi>e</mi> </msubsup> <mi>d</mi> <mi>F</mi> <mo>+</mo> <msubsup> <mi>F</mi> <mi>G</mi> <mi>e</mi> </msubsup> <msub> <mi>dL</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
式中：和为相应的单元柔度矩阵，dL0为初始位移增量，dF为单元不平衡力增量。
柔度矩阵求逆得到：
<mrow> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msubsup> <mi>F</mi> <mi>C</mi> <mi>e</mi> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msup> <mi>dF</mi> <mi>e</mi> </msup> <mo>=</mo> <msubsup> <mi>K</mi> <mi>T</mi> <mi>e</mi> </msubsup> <msup> <mi>dX</mi> <mi>e</mi> </msup> <mo>+</mo> <msubsup> <mi>K</mi> <mi>G</mi> <mi>e</mi> </msubsup> <msub> <mi>dL</mi> <mn>0</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
其中：和为响应的刚度矩阵，dXe为单元节点位移增量向量，dFe为端点力增量向量。
<mrow> <mtable> <mtr> <mtd> <mrow> <msup> <mi>dF</mi> <mi>e</mi> </msup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>d</mi> <msub> <mi>F</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mi>d</mi> <msub> <mi>F</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>4</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>5</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dF</mi> <mn>6</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <msup> <mi>dX</mi> <mi>e</mi> </msup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>d</mi> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dy</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dz</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dx</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dy</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>dz</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <msubsup> <mi>K</mi> <mi>T</mi> <mi>e</mi> </msubsup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> </mrow> </mtd> <mtd> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>K</mi> <mi>G</mi> <mi>e</mi> </msubsup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> <msubsup> <mi>F</mi> <mi>G</mi> <mi>e</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>w</mi> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>+</mo> <msubsup> <mi>K</mi> <mi>C</mi> <mi>e</mi> </msubsup> <msubsup> <mi>F</mi> <mi>G</mi> <mi>e</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
S34：根据动、定网格节点位移生成广义局部刚度矩阵和不平衡力向量；
根据动网格节点位移向量xm、ym、zm和定网格节点位移向量xs、ys、zs，生成每个网格的单元刚度矩阵dFe、和
根据有限元方法生成整体刚度矩阵和
整体静力学方程为：
dF＝KTdX+KGdL0 (17)
式中：dX为动网格增量,dL0为初始长度增量；
S34：施加局部约束条件；
在所有网格节点x方向施加虚约束，即：消除矩阵KT和向量dX对应的列向量；
保证每个单元的总初始长度之和不变，即施加约束
在定网格点x、y、z方向上施加固定约束，即：消除矩阵KT和向量dX对应的行、列向量；
S35：计算动网格位移增量，更新动网格位移相应；
根据式(17)计算动网格增量和初始长度增量；更新动网格位移向量：xm、ym、zm；
S36：判断是否满足max[dX dL0]＜d，如果满足则转入步骤4，如果不满足则返回步骤S32。
6.根据权利要求5所述的一种用于高铁弓网动态行为动网格仿真的不平衡力消除方法，其特征在于，所述步骤4具体过程如下：
将动网格节点位移向量xm、ym、zm引入到弓网整体位移向量U中，重新生成整体质量矩阵M，整体刚度矩阵K和整体阻尼矩阵C，通过纽马克贝塔算法Newmark-β法求解该时刻的位移响应，并进入下一时刻t+△t进行新的迭代计算。
CN201711239896.9A 2017-11-30 2017-11-30 一种高铁弓网动态行为动网格不平衡力消除的仿真方法 Expired - Fee Related CN108021747B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711239896.9A CN108021747B (zh) 2017-11-30 2017-11-30 一种高铁弓网动态行为动网格不平衡力消除的仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711239896.9A CN108021747B (zh) 2017-11-30 2017-11-30 一种高铁弓网动态行为动网格不平衡力消除的仿真方法

Publications (2)

Publication Number Publication Date
CN108021747A true CN108021747A (zh) 2018-05-11
CN108021747B CN108021747B (zh) 2021-04-09

Family

ID=62077568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711239896.9A Expired - Fee Related CN108021747B (zh) 2017-11-30 2017-11-30 一种高铁弓网动态行为动网格不平衡力消除的仿真方法

Country Status (1)

Country Link
CN (1) CN108021747B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110362786A (zh) * 2019-06-14 2019-10-22 西南交通大学 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法
CN113191057A (zh) * 2021-05-12 2021-07-30 四川长虹空调有限公司 不平衡力方向的确定方法及装置
CN113255195A (zh) * 2021-07-02 2021-08-13 成都国佳电气工程有限公司 一种考虑实测分布式阻尼的接触网有限元建模方法
CN114757079A (zh) * 2022-04-29 2022-07-15 中铁第一勘察设计院集团有限公司 接触网棘轮卡滞故障的仿真方法及系统
CN114936488A (zh) * 2022-04-29 2022-08-23 中铁第一勘察设计院集团有限公司 接触网异物侵扰故障的仿真方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1696356A1 (en) * 2005-02-24 2006-08-30 Siemens Aktiengesellschaft Flow acoustic simulation with the Lattice-Boltzmann method
CN105631128A (zh) * 2015-12-29 2016-06-01 西南交通大学 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1696356A1 (en) * 2005-02-24 2006-08-30 Siemens Aktiengesellschaft Flow acoustic simulation with the Lattice-Boltzmann method
CN105631128A (zh) * 2015-12-29 2016-06-01 西南交通大学 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110362786A (zh) * 2019-06-14 2019-10-22 西南交通大学 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法
CN110362786B (zh) * 2019-06-14 2023-03-14 西南交通大学 一种高速铁路弓网耦合系统动力学响应分析的解耦建模方法
CN113191057A (zh) * 2021-05-12 2021-07-30 四川长虹空调有限公司 不平衡力方向的确定方法及装置
CN113191057B (zh) * 2021-05-12 2022-04-12 四川长虹空调有限公司 不平衡力方向的确定方法及装置
CN113255195A (zh) * 2021-07-02 2021-08-13 成都国佳电气工程有限公司 一种考虑实测分布式阻尼的接触网有限元建模方法
CN113255195B (zh) * 2021-07-02 2021-10-08 成都国佳电气工程有限公司 一种考虑实测分布式阻尼的接触网有限元建模方法
CN114757079A (zh) * 2022-04-29 2022-07-15 中铁第一勘察设计院集团有限公司 接触网棘轮卡滞故障的仿真方法及系统
CN114936488A (zh) * 2022-04-29 2022-08-23 中铁第一勘察设计院集团有限公司 接触网异物侵扰故障的仿真方法及系统

Also Published As

Publication number Publication date
CN108021747B (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
CN108021747A (zh) 一种高铁弓网动态行为动网格不平衡力消除的仿真方法
Majka et al. Effects of speed, load and damping on the dynamic response of railway bridges and vehicles
Zhang et al. Study on dynamics of coupled systems in high-speed trains
CN105631128B (zh) 一种高速铁路弓-网-车-轨垂向耦合大系统动力学建模仿真方法
Zhu et al. Vibration of elevator cables with small bending stiffness
Van Do et al. Dynamic responses of Euler–Bernoulli beam subjected to moving vehicles using isogeometric approach
Li et al. Nonstationary Random Vibration Performance of Train‐Bridge Coupling System with Vertical Track Irregularity
Gu Resonance in long-span railway bridges carrying TGV trains
Sadeghi et al. An efficient algorithm for nonlinear analysis of vehicle/track interaction problems
CN107194120A (zh) 一种基于有限质点法的覆冰输电线路找形方法
CN104992025A (zh) 一种基于科瑟拉弹性杆模型的线缆的建模方法及装置
CN109398020A (zh) 一种基于非线性模型的车辆液电耦合式isd悬架的预测控制方法
CN114218699A (zh) 基于绝对节点坐标及附加约束法的刚性接触网建模方法
Zhang et al. Nonlinear dynamic analysis of a stochastic delay wheelset system
Qiao et al. Dynamic analysis of an integrated train–bridge–foundation–soil system by the substructure method
JP2019194747A (ja) 鉄道複線桁のねじり振動評価方法
Ding et al. Status and progress of kansei engineering design method
Zeng et al. Three-dimensional rail-bridge coupling element of unequal lengths for analyzing train-track-bridge interaction system
CN103926947A (zh) 大跨度索桥结构非线性系统半主动振动控制方法
Hung et al. Predicting dynamic responses of frame structures subjected to stochastic wind loads using temporal surrogate model
CN107092730A (zh) 适用于显式分析的三维无限元人工边界建立方法
CN114970134A (zh) 含惯容减振系统的野外作业平台动力学建模及参数优化设计方法
Pan et al. Curved Ballasted Track–Vehicle Dynamic Interaction: Effects of Curve Radius and Track Structural Nonlinearity
Liu et al. Estimation of wheelset natural vibration characteristics based on transfer matrix method with various elastic beam models
Liu et al. A non-iterative integration scheme enriching the solution to the coupled maglev vehicle–bridge system

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: 20210409

Termination date: 20211130

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