[发明专利]精确获取变长度绳索设备横向振动的方法

摘要

本发明公开了一种精确获取变长度绳索设备横向振动的方法，其轴向移动设备为轴向移动的绳索设备，绳索设备的运动工况为绳索长度呈线性变化，其特征是按照绳索设备中行波运动规律将振动周期分为三个阶段；结合绳索设备中行波移动规律和边界条件分别获得三个阶段的绳索设备的边界入射波和边界反射波；将边界入射波和边界反射波进行叠加，分别获得所述绳索设备在三个阶段的横向位移。本发明适用于绳索长度呈线性变化的移动绳索设备的多种复杂边界条件以及多种速度工况，获得的振动位移响应精确，能满足检验轴向移动绳索横向振动多种数值计算方法的可行性和有效性的需要。

主权项

1.一种精确获取变长度绳索设备横向振动的方法，轴向移动设备为轴向移动的绳索设备，所述绳索设备的运动工况为绳索长度呈线性变化，其特征是所述方法按如下过程进行：确定绳索设备的运动模型，获得所述绳索设备的运动方程；确定绳索设备的运动初始条件和边界条件；确定绳索设备的振动周期T；按照所述绳索设备中行波运动规律将所述振动周期T分为三个阶段：分别是[tn‑1,tan]，[tan,tbn]和[tbn,tn]，其中tn‑1式(3)中：ω_t为ω对t的一阶偏导数；函数为固定坐标系中绳索设备上不同位置的初始横向位移；函数ψ(x)为固定坐标系中绳索设备上不同位置的初始速度；l(t)为绳索设备在时刻t的长度；设置边界条件为式(4)所示的两端固定的形式：步骤3：针对绳索设备的运动工况为绳索长度呈线性变化，有式(5)如下：l(t)＝l0+vt             (5)式(5)中l0为绳索的初始长度；第n个周期的时长Tn为：式(5)中，ln‑1为第n‑1个周期结束时的绳索长度；利用式(5)和式(6)计算获得各振动周期的时长：T1，T2...Tn，n为自然数；步骤4：按照行波移动变化规律将振动周期分为三个阶段：定义tn为系统经过n次循环回到初始状态的时间，则有：Tn＝tn‑tn‑1                (8)当v>0，绳索从左向右移动，即绳索伸长，则在周期Tn内的三个阶段分别是[tn‑1,tan]，[tan,tbn]和[tbn,tn]，其中tn‑1当v<0，绳索从右向左移动，即绳索长度缩短，则在三个阶段中tan和tbn值相互交换；步骤5：结合运动初始条件和边界条件获取绳索设备在第一阶段[tn‑1,tan]的横向振动位移：第一阶段中，左移行波G分为行波G1和行波G2，右移行波F分为行波F1和行波F2；其中，行波G2是行波F1在右边界x＝l(t)处的反射波，行波F2是行波G1在左边界x＝0的处的反射波；结合运动初始条件和式(2)分别得到行波F1和行波G1的表达式如式(10)和式(11)：其中，ξ为积分变量；K为积分常数；结合边界条件和式(2)得到行波F2、行波G2与行波F1、行波G1的关系如式(12)：由式(10)、式(11)和式(12)得到式(13)和式(14)：则，第一阶段[tn‑1,tan]右移行波F和左移行波G分别如式(15)、式(16)所表征：利用式(17)获得第一阶段移动绳索的横向振动位移ω(x,t)ω(x，t)＝F(x‑vrt)+G(x+vlt)，0≤x≤l(t)、tn‑1≤t≤tan           (17)步骤6：结合初始条件和边界条件获取移动绳索在第二阶段[tan,tbn]的横向振动位移：在第二阶段中，对于v>0，左移行波G分为行波G1、行波G2和行波G3，右移行波F即为行波F2；行波F2为入射波，G3是F2在右边界x＝l(t)处的反射波；由边界条件和式(2)得到G3和F2的关系如式(18)：G3(x)＝‑F2(2ln‑1‑x)                  (18)结合式(13)和式(18)，获得行波G3(x+vlt)如式(19)：对于v<0，左移行波G即为行波G2，右移行波F分为行波F1、行波F2和行波F3；行波F3是行波G2在左边界x＝0处的反射波；利用边界条件和式(2)得到F3和G2的关系如式(20)：结合式(14)和式(20)，获得行波F3(x‑vrt)如式(21)：则，第二阶段[tan,tbn]右移行波F和左移行波G分别如式(22)和式(23)：利用式(24)获得第二阶段移动绳索的横向振动位移ω(x,t)：ω(x，t)＝F(x‑vrt)+G(x+vlt)，0≤x≤l(t)、tan≤t≤tbn          (24)步骤7：结合初始条件和边界条件获取移动绳索在第三阶段[tbn,tn]的横向振动位移：在第三阶段中，左移行波G分为行波G2和行波G3，右移行波F分为行波F2和行波F3，行波G3是行波F2在右边界x＝l(t)处的反射波，行波F3是行波G2在左边界x＝0处的反射波；右移行波F和左移行波G分别如式(25)和式(26)：利用式(27)获得第三阶段移动绳索的横向振动位移ω(x,t)为：ω(x，t)＝F(x‑vrt)+G(x+vlt)，0≤x≤l(t)、tbn≤t≤tn       (27)。