CN116384162A - 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备 - Google Patents

基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备 Download PDF

Info

Publication number
CN116384162A
CN116384162A CN202310650699.5A CN202310650699A CN116384162A CN 116384162 A CN116384162 A CN 116384162A CN 202310650699 A CN202310650699 A CN 202310650699A CN 116384162 A CN116384162 A CN 116384162A
Authority
CN
China
Prior art keywords
matrix
virtual spring
energy
complex
rail
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
CN202310650699.5A
Other languages
English (en)
Other versions
CN116384162B (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.)
East China Jiaotong University
Original Assignee
East China 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 East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN202310650699.5A priority Critical patent/CN116384162B/zh
Publication of CN116384162A publication Critical patent/CN116384162A/zh
Application granted granted Critical
Publication of CN116384162B publication Critical patent/CN116384162B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备，基于能量法框架，利用虚拟弹簧模拟周期边界条件，将边界条件转化为虚拟弹簧的弹性势能；对总能量泛函变分得到周期性轨道结构的运动特征方程；将虚拟弹簧的刚度矩阵展开，并将其含有波数的项和不含波数的项解耦；将周期性轨道结构的运动特征方程进行降价处理，将其转化为线性的特征值矩阵；通过扫描频率，并对线性的特征值矩阵变分得到复波数，进而得到轨道结构的复能带。本发明提供了一种新的周期轨道结构的复能带计算方法，该方法的方法简单，思路清晰，计算效率高且具有能量法在计算耦合结构中的优势。

Description

基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备
技术领域
本发明属于交通运输工程技术领域，具体涉及一种基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备。
背景技术
轨道结构是一种典型的周期结构，周期结构具有带隙特性，这一特性使得结构在特定频率范围内的声波/弹性波能够顺利通过，而其它频带则因衰减而无法通过。周期结构的带隙特性可以由实能带和复能带两个方面表征。实能带结构通过给定波数的实数取值，从而求解波的衰减频率，仅能反应无衰减波的传播规律，难以反应衰减效应。与之相反的，复能带结构是给定频率来求解波数。因此不仅能够刻画波的传播规律还能获取波的衰减细节。
传递矩阵法是一种天然的复能带解法。该方法是通过状态向量传递矩阵的特征值分解来计算带隙特性，由于传递矩阵元素中包含动刚度项（与频率有关），却不含波数项，因此可以方便的得到复能带。然而，传递矩阵法的几何适用性较差，主要应用于简单一维模型的复能带，如离散支撑钢轨。
扩展平面波的复能带结构计算流程直观，能够计算二维周期结构的复能带。但该方法的不足在于：对两组元材料差异较大的周期结构的收敛性较差，且难以处理组合结构中各子结构的耦合，难以计算无砟轨道模型。
有限元法的在计算复杂周期结构中具有巨大优势，是一种好用的分析方法，但计算精度依赖于网格划分的合理程度。意味着对模型做微小改动也必须重新划分网格，严重影响几何参数分析和优化时的计算效率。
能量法具有将泛函边值问题转化为极值问题的特点，这有利于耦合问题的求解。这也在日益复杂的周期结构计算中，具有较大的优势。然而，目前利用能量法计算周期性轨道结构的复能带的方法很少。
发明内容
为了解决上述现有技术的问题，本发明利用虚拟弹簧模型，提出基于虚拟弹簧模型的轨道结构复能带计算方法，该方法简单，思路清晰，计算效率高且具有能量法在计算耦合结构中的优势。
为解决上述技术问题，本发明的技术方案是：基于虚拟弹簧模型的轨道结构复能带计算方法，包括下述步骤：
S1：将轨道结构的周期边界条件用虚拟弹簧模型模拟，将周期边界条件转化为虚拟弹簧的弹性势能；根据轨道结构的材料参数得到钢轨和扣件的应变能、以及钢轨和扣件的动能；
S2：对总能量泛函变分得到周期性轨道结构的运动特征方程；
S3：将虚拟弹簧的刚度矩阵展开，并将其含有波数的项和不含波数的项解耦；
S4：将周期性轨道结构的运动特征方程进行降阶处理，将其转化为线性的特征值矩阵；
S5：通过扫描频率，并对线性的特征值矩阵变分得到复波数，进而得到轨道结构的复能带。
进一步优选，步骤S1中，对于周期性离散支撑钢轨，钢轨采用欧拉梁进行描述，扣件采用弹簧描述，扣件下部的轨枕和轨道板简化为刚体，将轨道结构划分为若干个胞元；扣件间距为
Figure SMS_1
，钢轨的垂向位移表示为/>
Figure SMS_2
，扣件刚度为/>
Figure SMS_3
；使用虚拟弹簧模型来模拟周期边界条件，利用虚拟弹簧将相邻胞元的头尾相连，以此模拟周期边界条件。
进一步优选，步骤S1中，根据能量法计算原理，钢轨的垂向位移由位移形函数
Figure SMS_4
和与时间相关的未知系数/>
Figure SMS_5
表示：
Figure SMS_6
，
式中未知系数矩阵
Figure SMS_7
，位移形函数矩阵
Figure SMS_8
, M表示用于模拟位移的位移形函数个数，即截断系数，
Figure SMS_9
为第m个未知系数，/>
Figure SMS_10
为第m个位移形函数，m∈1,2,…,M；T表示转置；根据布洛赫理论，位移和转角边界条件为：
Figure SMS_11
，
其中,
Figure SMS_12
表示轨道单元的起始位置的钢轨垂向位移，/>
Figure SMS_13
表示轨道单元的结束位置的钢轨垂向位移，/>
Figure SMS_14
表示轨道x方向的波数，i为复数的虚数单位。
进一步优选，步骤S1中，用于模拟周期边界条件的虚拟弹簧的弹性势能表示为：
Figure SMS_15
，
式中，
Figure SMS_16
为虚拟弹簧的弹性势能，/>
Figure SMS_17
表示虚拟弹簧的刚度矩阵，上标H表示共轭转置，/>
Figure SMS_18
为虚拟的位移弹簧的刚度，/>
Figure SMS_19
为虚拟的转角弹簧的刚度；
钢轨和扣件的应变能和动能表示为：
Figure SMS_20
，
Figure SMS_21
，
式中，
Figure SMS_24
表示钢轨和扣件的应变能，/>
Figure SMS_25
表示钢轨和扣件的动能；E，I，ρ和A分别表示钢轨的弹性模量、截面转动惯量、密度和横截面积；/>
Figure SMS_27
和/>
Figure SMS_23
分别为轨道结构一个胞元的刚度矩阵和质量矩阵，/>
Figure SMS_26
表示/>
Figure SMS_28
对时间的导数，/>
Figure SMS_29
表示未知系数矩阵/>
Figure SMS_22
对时间的导数。
进一步优选，步骤S2中，轨道结构的总能量泛函表示为：
Figure SMS_30
，根据拉格朗日方程/>
Figure SMS_31
，得到其运动特征方程：/>
Figure SMS_32
，/>
Figure SMS_33
为频率。
进一步优选，步骤S3中，将虚拟弹簧的刚度矩阵展开为：
Figure SMS_34
，
式中，
Figure SMS_36
表示位移形函数矩阵/>
Figure SMS_40
对x方向的一阶导数；/>
Figure SMS_44
为钢轨x方向坐标为0时的位移形函数矩阵；/>
Figure SMS_38
为钢轨x方向坐标为l时的位移形函数矩阵；/>
Figure SMS_42
为/>
Figure SMS_46
的共轭转置；/>
Figure SMS_48
为/>
Figure SMS_35
对x方向的一阶导数，/>
Figure SMS_39
为/>
Figure SMS_43
的共轭转置；/>
Figure SMS_47
表示与波数相关的量，/>
Figure SMS_37
，/>
Figure SMS_41
；而后将/>
Figure SMS_45
中含有波数的项和不含波数的项分离：
Figure SMS_49
，
式中，
Figure SMS_51
分别为/>
Figure SMS_55
展开后系数为/>
Figure SMS_57
，/>
Figure SMS_52
和1的矩阵；将/>
Figure SMS_54
展开后代入到运动特征方程/>
Figure SMS_56
，将已知矩阵：/>
Figure SMS_58
合并形成过渡矩阵
Figure SMS_50
；将/>
Figure SMS_53
解耦后的运动特征方程表示为：
Figure SMS_59
。
进一步优选，步骤S4中，定义一个降阶向量
Figure SMS_60
，并令/>
Figure SMS_61
，即：/>
Figure SMS_62
，将其代入方程/>
Figure SMS_63
中，得到降阶后的运动特征方程：
Figure SMS_64
。
进一步优选，步骤S4中，将
Figure SMS_65
代入降阶后的运动特征方程，形成线性的特征值矩阵：
Figure SMS_66
，
对特征值矩阵分解，求出
Figure SMS_67
，进而求出波数/>
Figure SMS_68
；输入不同的频率/>
Figure SMS_69
，便可得到对应频率下的波数/>
Figure SMS_70
。
本发明提供了一种非易失性计算机存储介质，计算机存储介质存储有计算机可执行指令，该计算机可执行指令可执行上述的基于虚拟弹簧模型的轨道结构复能带计算方法。
本发明还提供一种计算机程序产品，计算机程序产品包括存储在非易失性计算机存储介质上的计算机程序，计算机程序包括程序指令，当程序指令被计算机执行时，使计算机执行上述的基于虚拟弹簧模型的轨道结构复能带计算方法。
本发明提供一种电子设备，包括：至少一个处理器，以及与所述至少一个处理器通信连接的存储器，其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行基于虚拟弹簧模型的轨道结构复能带计算方法。
本发明的有益效果：本发明基于能量法的框架，利用虚拟弹簧模型模拟周期边界条件，将位移形函数与边界条件解耦，能够简单快捷的表示轨道结构胞元的位移场。通过对模拟周期边界条件的虚拟弹簧的刚度矩阵处理，得到能够变分求解的特征值方程，通过给定频率并对上述特征值方程变分便可得到复频散进而得到复能带结构。本发明为研究轨道结构复能带提供了一种新的计算方法。该方法计算过程简单，思路清晰，计算效率高且具有能量法在计算耦合结构中的优势。
附图说明
图1是基于虚拟弹簧模型的轨道结构复能带计算方法的步骤流程图；
图2是轨道结构的结构示意图；
图3是轨道结构计算模型图；
图4为本发明提出的方法得到的复能带的实数部分示意图；
图5为传递矩阵法得到的复能带的实数部分示意图；
图6为本发明提出的方法得到的复能带的虚数部分示意图；
图7为传递矩阵法得到的复能带的虚数部分示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例（离散支撑钢轨模型），对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，为本实施例提供的一种基于虚拟弹簧模型的轨道结构复能带计算方法，包括下述步骤S1至步骤S5，以下详细说明各个步骤的实现方式。
S1：将轨道结构的周期边界条件用虚拟弹簧模型模拟，将周期边界条件转化为虚拟弹簧的弹性势能；根据轨道结构的材料参数得到钢轨和扣件的应变能、以及钢轨和扣件的动能。
图2所示为周期性离散支撑钢轨，钢轨采用欧拉梁进行描述，扣件采用弹簧描述，扣件下部的轨枕和轨道板简化为刚体。轨道结构划分为若干个胞元，虚线框出的部分即为轨道结构的一个胞元，扣件间距为
Figure SMS_71
。钢轨的垂向位移表示为/>
Figure SMS_72
，扣件刚度为/>
Figure SMS_73
。本方法使用虚拟弹簧模型来模拟周期边界条件，轨道结构计算模型如图3所示。利用虚拟弹簧将相邻胞元的头尾相连，以此模拟周期边界条件。
根据能量法计算原理，钢轨的垂向位移由位移形函数
Figure SMS_74
和与时间相关的未知系数/>
Figure SMS_75
表示：
Figure SMS_76
，
式中未知系数矩阵
Figure SMS_77
, 位移形函数矩阵/>
Figure SMS_78
, M表示用于模拟位移的位移形函数个数，即截断系数，/>
Figure SMS_79
为第m个未知系数，/>
Figure SMS_80
为第m个位移形函数，m∈1,2,…,M；T表示转置；根据布洛赫(Bloch)理论，周期单元需要满足位移，转角，弯矩和剪力的周期边界条件，在能量法的运算框架下，可以简化为位移和转角边界条件为：
Figure SMS_81
，
其中,
Figure SMS_82
表示轨道单元的起始位置的钢轨垂向位移，/>
Figure SMS_83
表示轨道单元的结束位置的钢轨垂向位移，/>
Figure SMS_84
表示轨道x方向的波数，i为复数的虚数单位。
用于模拟周期边界条件的虚拟弹簧的弹性势能表示为：
Figure SMS_85
，
式中，
Figure SMS_86
为虚拟弹簧的弹性势能，/>
Figure SMS_87
表示虚拟弹簧的刚度矩阵，上标H表示共轭转置，/>
Figure SMS_88
为虚拟的位移弹簧的刚度，/>
Figure SMS_89
为虚拟的转角弹簧的刚度；
钢轨和扣件的应变能和动能表示为：
Figure SMS_90
，
Figure SMS_91
，
式中，
Figure SMS_94
表示钢轨和扣件的应变能，/>
Figure SMS_96
表示钢轨和扣件的动能；E，I，ρ和A分别表示钢轨的弹性模量、截面转动惯量、密度和横截面积；/>
Figure SMS_98
和/>
Figure SMS_93
分别为轨道结构一个胞元的刚度矩阵和质量矩阵，/>
Figure SMS_95
表示/>
Figure SMS_97
对时间的导数，/>
Figure SMS_99
表示未知系数矩阵/>
Figure SMS_92
对时间的导数。
S2：对总能量泛函变分得到周期性轨道结构的运动特征方程。
根据能量法的计算原理，轨道结构的总能量泛函表示为：
Figure SMS_100
，根据拉格朗日方程/>
Figure SMS_101
，得到其运动特征方程：/>
Figure SMS_102
，/>
Figure SMS_103
为频率。
对于实能带问题，扫描第一不可约布里渊区,
Figure SMS_104
，并运动特征方程进行特征值分解，便可得到轨道结构的实能带。与实能带算法相反，复能带算法是通过给定频率
Figure SMS_105
，来求解波数/>
Figure SMS_106
。然而，由于未知的波数/>
Figure SMS_107
出现在周期边界约束中而无法直接求解，我们考虑该运动特征方程转化为线性特征值问题。
S3：将虚拟弹簧的刚度矩阵展开，并将其含有波数的项和不含波数的项解耦。
为得到线性特征值方法，首先需要将波数从矩阵中分离出来。得益于虚拟弹簧模型，波数仅存在于虚拟弹簧的弹性势能中。因此，仅需将虚拟弹簧的刚度矩阵
Figure SMS_108
中的波数
Figure SMS_109
解耦。将虚拟弹簧的刚度矩阵展开为：
Figure SMS_110
，
式中，
Figure SMS_113
表示位移形函数矩阵/>
Figure SMS_118
对x方向的一阶导数；/>
Figure SMS_123
为钢轨x方向坐标为0时的位移形函数矩阵；/>
Figure SMS_112
为钢轨x方向坐标为l时的位移形函数矩阵；/>
Figure SMS_116
为/>
Figure SMS_120
的共轭转置；/>
Figure SMS_124
为/>
Figure SMS_111
对x方向的一阶导数，/>
Figure SMS_115
为/>
Figure SMS_119
的共轭转置；/>
Figure SMS_122
表示与波数相关的量，/>
Figure SMS_114
，/>
Figure SMS_117
；而后将/>
Figure SMS_121
中含有波数的项和不含波数的项分离：
Figure SMS_125
，
式中，
Figure SMS_127
分别为/>
Figure SMS_130
展开后系数为/>
Figure SMS_132
，/>
Figure SMS_128
和1的矩阵；将/>
Figure SMS_131
展开后代入到运动特征方程/>
Figure SMS_133
，将已知矩阵：/>
Figure SMS_134
合并形成过渡矩阵
Figure SMS_126
；将/>
Figure SMS_129
解耦后的运动特征方程表示为：
Figure SMS_135
。
S4：将周期性轨道结构的运动特征方程进行降阶处理，将其转化为线性的特征值矩阵。
将
Figure SMS_136
解耦后的运动特征方程的/>
Figure SMS_137
的次数设置为2。那么，需要降低λ的次数，这在数学上称为降阶。定义一个降阶向量/>
Figure SMS_138
，并令/>
Figure SMS_139
，即：/>
Figure SMS_140
，将其代入方程
Figure SMS_141
中，得到降阶后的运动特征方程：
Figure SMS_142
，
进一步的，将
Figure SMS_143
代入降阶后的运动特征方程，形成线性的特征值矩阵：
Figure SMS_144
，
对特征值矩阵分解，求出
Figure SMS_145
，进而求出波数/>
Figure SMS_146
；输入不同的频率/>
Figure SMS_147
，便可得到对应频率下的波数/>
Figure SMS_148
。
S5：通过扫描频率，并对线性的特征值矩阵变分得到复波数，进而得到轨道结构的复能带。
本实施例中，仅考虑钢轨的垂向振动，周期性离散支撑钢轨的材料参数如表1所示。为验证本实施例准确性，取相同材料参数，采用本发明提出的方法与传递矩阵法计算，获取轨道复能带结果如图4-图7所示，两者结果符合良好。如图4和图5所示，实线框部分为轨道结构的带隙频率区域，对比两种方法得到的轨道结构复能带的实数部分，可知采用本方法计算得到在129Hz以内和1340.5Hz~1353Hz存在两条带隙频段；由传递矩阵法计算得到的带隙频段为129 Hz以内和1342~1355 Hz。如图6和图7所示，对比两种方法得到的轨道结构复能带的虚数部分，可知129Hz以内带隙的最大衰减的波数为1.1，而1340.5Hz~1353Hz的带隙的最大衰减波数为1.1e-3。
Figure SMS_149
本实施例还提供了一种非易失性计算机存储介质，计算机存储介质存储有计算机可执行指令，该计算机可执行指令可执行上述任意实施例中的基于虚拟弹簧模型的轨道结构复能带计算方法。
本实施例还提供一种计算机程序产品，计算机程序产品包括存储在非易失性计算机存储介质上的计算机程序，计算机程序包括程序指令，当程序指令被计算机执行时，使计算机执行上述实施例的基于虚拟弹簧模型的轨道结构复能带计算方法。
本实施例还提供一种电子设备，包括：至少一个处理器，以及与所述至少一个处理器通信连接的存储器，其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行基于虚拟弹簧模型的轨道结构复能带计算方法。
上面所述的实施例仅仅是对本发明的实施方式进行描述，并非对本发明的构思和范围进行限定。在不脱离本发明设计构思的前提下，本领域普通人员对本发明的技术方案做出的各种变型和改进，均应落入到本发明的保护范围，本发明请求保护的技术内容，已经全部记载在权利要求书中。

Claims (10)

1.基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，包括下述步骤：
S1：将轨道结构的周期边界条件用虚拟弹簧模型模拟，将周期边界条件转化为虚拟弹簧的弹性势能；根据轨道结构的材料参数得到钢轨和扣件的应变能、以及钢轨和扣件的动能；
S2：对总能量泛函变分得到周期性轨道结构的运动特征方程；
S3：将虚拟弹簧的刚度矩阵展开，并将其含有波数的项和不含波数的项解耦；
S4：将周期性轨道结构的运动特征方程进行降阶处理，将其转化为线性的特征值矩阵；
S5：通过扫描频率，并对线性的特征值矩阵变分得到复波数，进而得到轨道结构的复能带。
2.根据权利要求1所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S1中，对于周期性离散支撑钢轨，钢轨采用欧拉梁进行描述，扣件采用弹簧描述，扣件下部的轨枕和轨道板简化为刚体，将轨道结构划分为若干个胞元；扣件间距为
Figure QLYQS_1
，钢轨的垂向位移表示为/>
Figure QLYQS_2
，扣件刚度为/>
Figure QLYQS_3
；使用虚拟弹簧模型来模拟周期边界条件，利用虚拟弹簧将相邻胞元的头尾相连，以此模拟周期边界条件。
3.根据权利要求2所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S1中，根据能量法计算原理，钢轨的垂向位移由位移形函数
Figure QLYQS_4
和与时间相关的未知系数/>
Figure QLYQS_5
表示：
Figure QLYQS_6
，
式中未知系数矩阵
Figure QLYQS_7
, 位移形函数矩阵/>
Figure QLYQS_8
,M表示用于模拟位移的位移形函数个数，即截断系数，/>
Figure QLYQS_9
为第m个未知系数，/>
Figure QLYQS_10
为第m个位移形函数，m∈1,2,…,M；T表示转置；根据布洛赫理论，位移和转角边界条件为：
Figure QLYQS_11
，
其中,
Figure QLYQS_12
表示轨道单元的起始位置的钢轨垂向位移，/>
Figure QLYQS_13
表示轨道单元的结束位置的钢轨垂向位移，/>
Figure QLYQS_14
表示轨道x方向的波数，i为复数的虚数单位。
4.根据权利要求3所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S1中，用于模拟周期边界条件的虚拟弹簧的弹性势能表示为：
Figure QLYQS_15
，
式中，
Figure QLYQS_16
为虚拟弹簧的弹性势能，/>
Figure QLYQS_17
表示虚拟弹簧的刚度矩阵，上标H表示共轭转置，/>
Figure QLYQS_18
为虚拟的位移弹簧的刚度，/>
Figure QLYQS_19
为虚拟的转角弹簧的刚度；
钢轨和扣件的应变能和动能表示为：
Figure QLYQS_20
，
Figure QLYQS_21
，
式中，
Figure QLYQS_23
表示钢轨和扣件的应变能，/>
Figure QLYQS_25
表示钢轨和扣件的动能；E，I，ρ和A分别表示钢轨的弹性模量、截面转动惯量、密度和横截面积；/>
Figure QLYQS_27
和/>
Figure QLYQS_24
分别为轨道结构一个胞元的刚度矩阵和质量矩阵，/>
Figure QLYQS_26
表示/>
Figure QLYQS_28
对时间的导数，/>
Figure QLYQS_29
表示未知系数矩阵/>
Figure QLYQS_22
对时间的导数。
5.根据权利要求4所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S2中，轨道结构的总能量泛函表示为：
Figure QLYQS_30
，根据拉格朗日方程/>
Figure QLYQS_31
，得到其运动特征方程：/>
Figure QLYQS_32
，/>
Figure QLYQS_33
为频率。
6.根据权利要求5所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S3中，将虚拟弹簧的刚度矩阵展开为：
Figure QLYQS_34
，
式中，
Figure QLYQS_36
表示位移形函数矩阵/>
Figure QLYQS_39
对x方向的一阶导数；/>
Figure QLYQS_43
为钢轨x方向坐标为0时的位移形函数矩阵；/>
Figure QLYQS_37
为钢轨x方向坐标为l时的位移形函数矩阵；/>
Figure QLYQS_41
为/>
Figure QLYQS_45
的共轭转置；
Figure QLYQS_48
为/>
Figure QLYQS_35
对x方向的一阶导数，/>
Figure QLYQS_40
为/>
Figure QLYQS_44
的共轭转置；/>
Figure QLYQS_47
表示与波数相关的量，
Figure QLYQS_38
，/>
Figure QLYQS_42
；而后将/>
Figure QLYQS_46
中含有波数的项和不含波数的项分离：
Figure QLYQS_49
，
式中，
Figure QLYQS_51
分别为/>
Figure QLYQS_54
展开后系数为/>
Figure QLYQS_56
，/>
Figure QLYQS_52
和1的矩阵；将/>
Figure QLYQS_55
展开后代入到运动特征方程/>
Figure QLYQS_57
，将已知矩阵：/>
Figure QLYQS_58
合并形成过渡矩阵/>
Figure QLYQS_50
；将/>
Figure QLYQS_53
解耦后的运动特征方程表示为：
Figure QLYQS_59
。
7.根据权利要求6所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S4中，定义一个降阶向量
Figure QLYQS_60
，并令/>
Figure QLYQS_61
，即：/>
Figure QLYQS_62
，将其代入方程/>
Figure QLYQS_63
中，得到降阶后的运动特征方程：
Figure QLYQS_64
。
8.根据权利要求7所述的基于虚拟弹簧模型的轨道结构复能带计算方法，其特征在于，步骤S4中，将
Figure QLYQS_65
代入降阶后的运动特征方程，形成线性的特征值矩阵：
Figure QLYQS_66
，
对特征值矩阵分解，求出
Figure QLYQS_67
，进而求出波数/>
Figure QLYQS_68
；输入不同的频率/>
Figure QLYQS_69
，便可得到对应频率下的波数/>
Figure QLYQS_70
。
9.一种计算机程序产品，计算机程序产品包括存储在非易失性计算机存储介质上的计算机程序，计算机程序包括程序指令，其特征在于，当程序指令被计算机执行时，使计算机执行权利要求1所述的基于虚拟弹簧模型的轨道结构复能带计算方法。
10.电子设备，包括：至少一个处理器，以及与所述至少一个处理器通信连接的存储器，其中，所述存储器存储有可被所述至少一个处理器执行的指令，其特征在于，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1所述的基于虚拟弹簧模型的轨道结构复能带计算方法。
CN202310650699.5A 2023-06-05 2023-06-05 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备 Active CN116384162B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310650699.5A CN116384162B (zh) 2023-06-05 2023-06-05 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310650699.5A CN116384162B (zh) 2023-06-05 2023-06-05 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备

Publications (2)

Publication Number Publication Date
CN116384162A true CN116384162A (zh) 2023-07-04
CN116384162B CN116384162B (zh) 2023-08-08

Family

ID=86979106

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310650699.5A Active CN116384162B (zh) 2023-06-05 2023-06-05 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备

Country Status (1)

Country Link
CN (1) CN116384162B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060144129A1 (en) * 2004-11-05 2006-07-06 Board Of Regents Of University Of Nebraska Method and apparatus for noncontact relative rail displacement, track modulus and stiffness measurement by a moving rail vehicle
CN105205035A (zh) * 2015-10-13 2015-12-30 哈尔滨工程大学 一种非均匀弹性约束边界条件矩形板结构面内振动分析方法
CN106758738A (zh) * 2016-12-20 2017-05-31 武汉理工大学 弹性转动约束边界钢‑混组合梁高腹板纵向加劲肋设置的方法
CN107808048A (zh) * 2017-10-26 2018-03-16 哈尔滨工程大学 一种均匀截面梁结构的振动分析方法
CN107818209A (zh) * 2017-10-26 2018-03-20 哈尔滨工程大学 一种弹性板结构的振动分析方法
CN112560320A (zh) * 2020-12-24 2021-03-26 华东交通大学 一种轨道结构弯曲振动带隙的计算方法
CN112613114A (zh) * 2020-11-25 2021-04-06 上海航天控制技术研究所 一种含摩擦边界的板结构模态求解方法
US20210334423A1 (en) * 2020-04-27 2021-10-28 Suzhou University of Science and Technology Method for calculation of natural frequency of multi-segment continuous beam
US20220147666A1 (en) * 2020-11-07 2022-05-12 Southwest Jiaotong University Method for determining propagation characteristics of guided waves of variable cross-section rail of turnout
US20220284152A1 (en) * 2020-02-25 2022-09-08 Shandong University Peridynamics method and system for tunnel rock mass failure water inrush catastrophe simulation

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060144129A1 (en) * 2004-11-05 2006-07-06 Board Of Regents Of University Of Nebraska Method and apparatus for noncontact relative rail displacement, track modulus and stiffness measurement by a moving rail vehicle
CN105205035A (zh) * 2015-10-13 2015-12-30 哈尔滨工程大学 一种非均匀弹性约束边界条件矩形板结构面内振动分析方法
CN106758738A (zh) * 2016-12-20 2017-05-31 武汉理工大学 弹性转动约束边界钢‑混组合梁高腹板纵向加劲肋设置的方法
CN107808048A (zh) * 2017-10-26 2018-03-16 哈尔滨工程大学 一种均匀截面梁结构的振动分析方法
CN107818209A (zh) * 2017-10-26 2018-03-20 哈尔滨工程大学 一种弹性板结构的振动分析方法
US20220284152A1 (en) * 2020-02-25 2022-09-08 Shandong University Peridynamics method and system for tunnel rock mass failure water inrush catastrophe simulation
US20210334423A1 (en) * 2020-04-27 2021-10-28 Suzhou University of Science and Technology Method for calculation of natural frequency of multi-segment continuous beam
US20220147666A1 (en) * 2020-11-07 2022-05-12 Southwest Jiaotong University Method for determining propagation characteristics of guided waves of variable cross-section rail of turnout
CN112613114A (zh) * 2020-11-25 2021-04-06 上海航天控制技术研究所 一种含摩擦边界的板结构模态求解方法
CN112560320A (zh) * 2020-12-24 2021-03-26 华东交通大学 一种轨道结构弯曲振动带隙的计算方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
冯青松等: "周期性钢弹簧浮置板轨道垂向振动带隙特性研究", 《中国铁道科学》, vol. 43, no. 1, pages 17 - 28 *
冯青松等: "基于人工弹簧模型的周期结构带隙计算方法研究", 《力学学报》, vol. 53, no. 6, pages 1684 - 1697 *
周渤;石先杰;: "连续多跨梁结构振动特性分析", 机械设计与制造, no. 08 *
张俊;李天匀;朱翔;: "多开口矩形板自由振动特性分析", 振动与冲击, no. 14 *
张俊;李天匀;朱翔;郭文杰;陈繁;: "中心典型形状开口的矩形薄板自由振动特性分析", 中国舰船研究, no. 02 *
盛曦;赵才友;王平;陈俊豪;魏晓;: "整体道床轨道扣件刚度对钢轨声功率特性的影响", 西南交通大学学报, no. 05 *

Also Published As

Publication number Publication date
CN116384162B (zh) 2023-08-08

Similar Documents

Publication Publication Date Title
Fu et al. Hybrid FEM–SBM solver for structural vibration induced underwater acoustic radiation in shallow marine environment
Ni et al. Using polynomial chaos expansion for uncertainty and sensitivity analysis of bridge structures
Liu et al. Analytical solutions for forced long waves on a sloping beach
Zhao et al. Seismic fragility analysis of nuclear power plants considering structural parameter uncertainty
Ritto et al. A new measure of efficiency for model reduction: application to a vibroimpact system
Zuo et al. Numerical and experimental investigations on the vibration band-gap properties of periodic rigid frame structures
Shen et al. Construction of peridynamic beam and shell models on the basis of the micro-beam bond obtained via interpolation method
CN112949121A (zh) 一种导波传播特性的求解方法及系统
Ma et al. Nonlinear waves in a floating thin elastic plate, predicted by a coupled SPH and FEM simulation and by an analytical solution
Gatti et al. Investigation of the earthquake ground motion coherence in heterogeneous non-linear soil deposits
Karmakar et al. Conditional simulation of non-Gaussian wind velocity profiles: Application to buffeting response of Vincent Thomas suspension bridge
Huang et al. A Bayesian probabilistic approach for damage identification in plate structures using responses at vibration nodes
Li et al. A nonmode-shape-based model updating method for offshore structures using extracted components from measured accelerations
Schneider et al. Bayesian parameter updating in linear structural dynamics with frequency transformed data using rational surrogate models
CN116384162B (zh) 基于虚拟弹簧模型的轨道结构复能带计算方法及电子设备
Muratov et al. Grid-characteristic method on unstructured tetrahedral meshes
CN111400898A (zh) 基于主模态法和应变的阵列天线振动变形预测方法及设备
Skrinar Elastic beam finite element with an arbitrary number of transverse cracks
Chen et al. Partitioned analysis of soil‐structure interaction for Nuclear Island Buildings
Zuo et al. Wave propagation of laminated composite plates via GPU-based wavelet finite element method
Van Rij et al. Structural loads analysis for wave energy converters
Golovko et al. Solution of axisymmetric dynamic problems for cylindrical shells on an elastic foundation.
Cheng et al. Wave energy conversion by an array of oscillating water columns deployed along a long-flexible floating breakwater
Tabatabaie SASSI FE Program for Seismic Response Analysis of Nuclear Containment Structures
Hoover et al. Immersed boundary thin shell analysis using 3D B-Spline background mesh

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