塔式多绳摩擦提升机主振动的计算机模拟.docx

塔式多绳摩擦提升机主振动的计算机模拟 0 动态动力特性仿真研究 为了研究提升机的整体振动，可以通过传递矩阵法简单地计算各阶的固有频率和相应的主振源，了解整体系统的振动特性。这为提升机动态设计提供了强有力的基础和可行的方法，但在某些初始条件下，系统的响应无法深入理解，也无法理解特定一阶主振动下的振动就是主振动。在这项工作中，我们调查了基于系统的振动，并得出了基于系统的振动计算公式。使用计算机以动态分析每个群落的主振动，为深入理解主振动提供了一种直观、可靠的方法，并为系统在任何初始条件下的响应提供了一种可行的方法。 1 计算公共公式的导出 1.1 刚性铰控制链状结构 对于提升机整体振动的研究,一般将提升机离散化为一链状结构来研究,如图1所示.在此设此链状结构具有20个质点,其中1,6,16,20点均为刚性铰支点.取y为广义坐标来表示各集中质量点的位移值,横坐标值x确定各集中质量点在主轴上的位置. 1.2 结构系统的动力传递 系统经离散化后各集中质量点的质量矩阵 mi(i=1,2,…20)依次为各集中质量点的质量. 由传递矩阵法可求出此链状结构系统的各阶固有频率及主振型为: 1.3 正则坐标的运动方程 由模态分析法可求出系统的模态质量矩阵 [M]=[A]T[m][A]=???????M1M2?M20???????. 根据下式可以将上述模态矩阵[A]正则化为正则模态矩阵[AN]: [AN]=[A][α], 其中, 矩阵中αi(i=1,2,…,20)为正则化因子. 对应于正则模态矩阵的正则坐标为{YN},且正则坐标与原广义坐标{Y}的关系为 {Y}=[AN]{YN}. (1) 当采用正则坐标时,系统的运动方程式则呈下述形式: 对于正定系统,可求出各正则坐标的一般解为 YNi=Aicos?ωnit+Bisin?ωnit, (i=1,2,…,20), (2) 待定常数Ai,Bi可以由初始时刻t=0时各正则坐标及其速度的初始值YNi,Y˙Ni表示.则(1)式可改写为 YNi=YNi0cos?ωnit+Y˙Ni0ωNisinωnit,(i=1,2,…,20), (3) 式中YNi0,Y˙Ni0分别表示各点正则坐标初始值及正则坐标表示的速度的初始值.可由下式求出: {YN}t=0=[AN]T[m]{Y}t=0, (4) {Y˙N}t=0=[AN]T[m]{Y˙}t=0, (5) 将(4)(5)两式代入(3)式,并将计算结果代入(1)式,得 在此,我们取初始条件为 其中{Y(r)N},{Y˙(r)N}分别为正则坐标表示的第r主振动下各集中质量点的位移和速度. 这样,激发的就是第r阶主振型下的振动,即第r主振动.则可得下式: 为便于计算机运算,我们将上式进一步改写为 由(6)式可以求出第r阶主振型下各集中质量点在t时刻的位移值,从而可以确定各点在该时刻的位置. 2 质量点位置的确定 由上可知,在求出系统的各阶主振型和各阶固有频率后,我们就可以求出系统的正则模态矩阵,进而根据公式(6)计算出在第r主振动下各集中质量点在某一时刻t的位移值.这样,就可以确定各集中质量点在该时刻的位置.根据数值计算方法,我们可以用分段均差插值法将该时刻整个链状系统的振动曲线求出.继而可以求出整个周期T(T=2πωnr)内每相隔一定时间间隔Δt各点的位移值.这样,对于整个系统在一个周期内的振动情况,我们就可以以动画的形式将其模拟出来.通过选择r的值,可以控制程序对哪一阶主振型下系统的振动即哪一阶主振动进行模拟. 2.1 程序脚本 下面给出了本文所用计算机程序的流程(图2): 2.2 阶主振型下的整体振动 该程序在读入已知数据后,能够根据用户的要求计算出任一阶主振型下系统各集中质量点在任意时刻的位移值,从而以动画的形式模拟出在该阶主振型下系统的整体振动,即与该阶主振型相对应的主振动.让我们能直观、生动的看到系统在任一阶主振型下的振动情况.给我们了解系统的整体振动特性提供了直观的感性认识,也有助于我们分析系统的实际振动情况. 3 计算机算法验证 模拟取得了较为理想的效果.将同一阶主振型下的计算机模拟结果与相应的振型图进行比较,可看到通过两种方式所得到的结果完全吻合,都能反映出系统的整体振动特性.图3中a为计算机模拟的一阶主振动在各集中质量点的位移达到最大值时的图像,b是由传递矩阵法求出的一阶振型图.图3所示结果显示,主振动中各集中质量点的位移比值基本保持振型图中所显示的位移比值.这也证明了本论文结果的正确性. 不足之处在于,采用的计算机算法存在着不可避免的误差.在最后一个集中质量点(第20点)左边存在由于算法问题而造成的较大误差.另外,这里所研究的主振动就是系统对特定初始条件(各点初始位移的比值与某一阶振型相符)的响应.而对系统