复杂边界下可变长度轴向移动弦线系统侧向振动获取方法技术方案

本发明专利技术公开了一种复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，包括如下步骤：步骤一：根据动力学特性和力学原理，建立可变长度轴向移动弦线系统在质量

【技术实现步骤摘要】
复杂边界下可变长度轴向移动弦线系统侧向振动获取方法


[0001]本专利技术属于机械振动测量
，具体的为一种复杂边界下可变长度轴向移动弦线系统侧向振动获取方法。

技术介绍

[0002]轴向移动系统具有阻尼低、刚度小、质量轻、柔韧性好和能耗低等优点，被广泛地应用于现代工程领域中,如航空航天、海洋、电子、土木、机械和汽车等工程等。噪声和振动伴随着这些设备的运转，尤其是横向振动对这些设备的功能和安全造成了很大的影响。
[0003]传统的研究技术是基于哈密顿原理建立的偏微分运动方程以及基于拉格朗日方程建立的有限单元动力学方程，利用数值计算方法，如伽辽金法、龙格库塔法、纽马克法以及时变状态空间方程等方法求解以上方程，以获得轴向移动绳索设备的横向振动响应。在进行理论分析时，解析解是更为有力的方法，但当前一般的解析解技术仅仅限于一个周期内的计算，如公开号为CN109614745A的中国专利提出的一种获取混合边界条件下轴向移动绳索设备横向振动的方法，该方法提出了一种获取轴向移动绳索一个振动周期T中各阶段的反射波响应方程的方法，但该方法无法在任意长时间进行计算。

技术实现思路

[0004]有鉴于此，本专利技术的目的在于提供一种复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，旨在获取轴向移动弦线系统在可变长度和复杂边界条件下的侧向振动。
[0005]为达到上述目的，本专利技术提供如下技术方案：
[0006]一种复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，包括如下步骤：
[0007]步骤一：根据动力学特性和力学原理，建立可变长度轴向移动弦线系统在质量
‑
阻尼
‑
弹簧边界和Dirichlet边界条件下的动力学模型；
[0008]可变长度轴向移动弦线系统的弦线的两端中，第一端为激励输入端，第二端为等效集中质量端；以弦线的第一端为原点建立固定坐标系OXY，可变长度轴向移动弦线系统在质量
‑
阻尼
‑
弹簧边界和Dirichlet边界条件下的动力学模型为：
[0009][0010]其中：ρ，T和V分别为弦线的线密度，静态张力和恒定移动速度；x表示弦线上的点在X方向上与原点距离；w(x，t)表示弦线上与原点距离为x的点的横向位移；l0为弦线的初始长度；l(t)为弦线的实时长度；t表示时间；
[0011]边界条件为：
[0012][0013]其中，m
e
，c
e
和k
e
分别为等效集中质量，等效阻尼和等效刚度；e(t)表示第一端输入的位移激励；表示横向位移对时间的一阶导数；表示横向位移对时间的二阶导数；w
′
表
示横向位移对空间的一阶导数；w
″
表示横向位移对空间的二阶导数；表示横向位移对时间、空间的二阶混合导数；
[0014]步骤二：将计算域转换为特征线域，在特征线域的端点建立弹性波反射方程，所述弹性波反射方程包括在Dirichlet边界条件下的波反射方程和在质量
‑
阻尼
‑
弹簧的复杂边界条件下的波反射方程；
[0015]步骤三：采用Duhamel积分进行求解，得到波反射方程解析解；
[0016]步骤四：根据小步长内波形线性假设，采用快速递推的方法，得到特征线域内波反射方程半解析解；
[0017]步骤五：根据左右行波的叠加方法，将特征线域的计算结果叠加得到原始时空域，得到可变长度轴向移动弦线系统的振动结果。
[0018]进一步，所述步骤二中，将计算域转换为特征线域，将动力学模型的解w表示为两个行波叠加的形式：
[0019]w＝F(x
‑
V
r
t)+G(x+V
l
t)
[0020]V
r
＝c十V
[0021]V
l
＝c
‑
V
[0022]其中：V
r
和V
l
表示固定坐标系中弦线朝第二端移动的第二行波的速度和朝第一端移动的第一行波的速度；表示弦线在静止条件下的波速；F(x
‑
V
r
t)为速度为V
r
的第二行波；G(x+V
l
t)为速度为V
l
的第一行波；
[0023]确定弦线振动的初始横向位移和速度条件：
[0024][0025]其中：函数φ(x)为固定坐标系中弦线上不同位置的初始横向位移；函数ψ(x)为固定坐标系中弦线上不同位置的初始速度；
[0026]得到第一行波和第二行波在初始时刻的表达式：
[0027][0028]其中，C为积分常数；
[0029]在x
‑
t域中使用第一横波和第二横波表示弦线中任意时间任意点的横向位移：
[0030]w(x，t)＝F(ξ)+G(η)
[0031][0032]其中：x表示在固定坐标系中点的位置；t表示任意时间；ξ表示引入的第二行波参
数，也即特征线域内的第二坐标；η表示引入的第一行波参数，也即特征线域内的第一坐标；
[0033]引入第一行波参数ξ和第二行波参数η，在ξ
‑
η域中表示行波，将求解动力学方程问题转化为求解边界上的第一行波和第二行波问题，其中：
[0034]第一行波参数ξ为：
[0035]ξ＝
‑
V
r
t
[0036]第二行波参数η为：
[0037]η＝l0+ct
[0038]其中：c表示弦线在静止条件下的波速。
[0039]进一步，所述步骤二中，在Dirichlet边界条件下的波反射方程为：
[0040]F(
‑
V
r
t)+G(V
l
t)＝e(t)
[0041]其中：F(
‑
V
r
t)表示第一端点的第二行波，G(V
l
t)表示第一端点的第一行波；
[0042]第一行波与第二反射行波之间的关系为：
[0043][0044]其中：F(ξ)表示第二行波；F
′
(ξ)表示第二行波对第二行波参数的一阶导数；F
″
(ξ)表示第二行波对第二行波参数的二阶导数；表示第一行波；表示第一行波对第一行波参数的一阶导数之后代入表示第一行波对第一行波参数的二阶导数之后代入为第一端点的位移激励在特征线域内的表示；Y
′
(ξ)表示第一端点的位移激励在特征线域内对第二行波参数的第一阶导数；Y
″
(ξ)表示第一端点的位移激励在特征线域内对第二行波参数的第二阶导数。
[0045]进一步，所述步骤二中，在质量
‑
阻尼
‑
弹簧的复杂边界条件下的波反射方程为本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，其特征在于：包括如下步骤：步骤一：根据动力学特性和力学原理，建立可变长度轴向移动弦线系统在质量
‑
阻尼
‑
弹簧边界和Dirichlet边界条件下的动力学模型；可变长度轴向移动弦线系统的弦线的两端中，第一端为激励输入端，第二端为等效集中质量端；以弦线的第一端为原点建立固定坐标系OXY，可变长度轴向移动弦线系统在质量
‑
阻尼
‑
弹簧边界和Dirichlet边界条件下的动力学模型为：其中：ρ，T和V分别为弦线的线密度，静态张力和恒定移动速度；x表示弦线上的点在X方向上与原点距离；w(x，t)表示弦线上与原点距离为x的点的横向位移；l0为弦线的初始长度；l(t)为弦线的实时长度；t表示时间；边界条件为：其中，m
e
，c
e
和k
e
分别为等效集中质量，等效阻尼和等效刚度；e(t)表示第一端输入的位移激励；表示横向位移对时间的一阶导数；表示横向位移对时间的二阶导数；w
′
表示横向位移对空间的一阶导数；w
″
表示横向位移对空间的二阶导数；表示横向位移对时间、空间的二阶混合导数；步骤二：将计算域转换为特征线域，在特征线域的端点建立弹性波反射方程，所述弹性波反射方程包括在Dirichlet边界条件下的波反射方程和在质量
‑
阻尼
‑
弹簧的复杂边界条件下的波反射方程；步骤三：采用Duhamel积分进行求解，得到波反射方程解析解；步骤四：根据小步长内波形线性假设，采用快速递推的方法，得到特征线域内波反射方程半解析解；步骤五：根据左右行波的叠加方法，将特征线域的计算结果叠加得到原始时空域，得到可变长度轴向移动弦线系统的振动结果。2.根据权利要求1所述的复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，其特征在于：所述步骤二中，将计算域转换为特征线域，将动力学模型的解w表示为两个行波叠加的形式：w＝F(x
‑
V
r
t)+G(x+V
l
t)V
r
＝c+VV
l
＝c
‑
V其中：V
r
和V
l
表示固定坐标系中弦线朝第二端移动的第二行波的速度和朝第一端移动的第一行波的速度；表示弦线在静止条件下的波速；F(x
‑
V
r
t)为速度为V
r
的第二行波；G(x+V
l
t)为速度为V
l
的第一行波；确定弦线振动的初始横向位移和速度条件：
其中：函数φ(x)为固定坐标系中弦线上不同位置的初始横向位移；函数ψ(x)为固定坐标系中弦线上不同位置的初始速度；得到第一行波和第二行波在初始时刻的表达式：其中，C为积分常数；在x
‑
t域中使用第一横波和第二横波表示弦线中任意时间任意点的横向位移：w(x，t)＝F(ξ)+G(η)其中：x表示在固定坐标系中点的位置；t表示任意时间；ξ表示引入的第二行波参数，也即特征线域内的第二坐标；η表示引入的第一行波参数，也即特征线域内的第一坐标；在ξ—η域中表示行波，将求解动力学方程问题转化为求解边界上的第一行波和第二行波问题，其中：第一行波参数ξ为：ξ＝
‑
V
r
t第二行波参数η为：η＝l0+ct其中：c表示弦线在静止条件下的波速。3.根据权利要求2所述的复杂边界下可变长度轴向移动弦线系统侧向振动获取方法，其特征在于：所述步骤二中，在Dirichlet边界条件下的波反射方程为：F(
‑
V
r
t)+G(V
l
t)＝e(t)其中：F(
‑
V
r
t)表示第一端点的第二行波，G(V
l
t)表示第一端点的第一行波；；第一行波与第二反射行波之间的关系为：
其中：F(ξ)表示第二行波；F
′
(ξ)表示第二行波对第二行波参数的一阶导数；F
″
(ξ)表示第二行波对第二行波参数的二阶导数；表示第一行波；表示第一行波对第一行波参数的一阶导数之后代入波对第一行波参数的一阶导数之后代入表示第一行波对第一行波参数的二阶导数之后代入Y(ξ)为第一端点的位移激励在特...

【专利技术属性】
技术研发人员：皮阳军，张建，彭越，张鑫宇，蒲华燕，罗均，
申请(专利权)人：重庆大学，
类型：发明
国别省市：