门式起重机虚拟样机建模与动力学仿真

门式起重机虚拟样机建模与动力学仿真摘要：研究了门式起重机动力学虚拟样机的建模方法，根据刚性体部件、柔性体部件和钢丝绳的特点，分别给出了不同的建模方法。通过对起升工况和大车运行工况进行动力学仿真，得到了整机的位移时间历程、动应力以及各连接件之间的作用力等，较真实地反映了门式起重机在作业过程中的动态特性。起重机在起动和制动的过程中，承受着强烈的冲击振动,这种冲击振动产生的动载荷在起重机运行过程中的影响是必须加以考虑的。长期以来，在起重机设计中都是将动态问题简化为静态问题处理，一些国家和国际起重机协会的起重机设计规范均采用一个动载系数花来考虑这种动态载荷。其最大的缺陷是不能够较为准确地反映起重机实际工况的动态性能。近年来，国内外不少研究者根据起重机的实际结构及受载特点，对起重机结构和机构的振动进行了大量的动态仿真计算和研究分析,使动态仿真方法逐步应用于起重机的设计计算分析中。本文主要探讨将虚拟样机技术运用于起重量为250t的门式起重机的动态设计，就门式起重机的起升工况和大车运行工况进行仿真研究，探索建立门式起重机动力学虚拟样机的方法。门式起重机的虚拟样机的建模:门式起重机的整机结构形式如图1所示。主要由小车、门架结构、大车运行机构、副机房、司机室、测风装置和栏杆梯子平台等组成。门式起重机的门架结构的主梁和支腿均采用箱型焊接结构，主梁和支腿之间通过法兰盘连接。小车位于门架主梁的轨道上,主起升机构和小车运行机构放置于小车上。副机房固定在门架主梁的伸出端，其作用是控制副起升。1.副机房2.小车3.测风装置4,主吊钩5.司机室6.大车运行机构7,栏杆梯子平台8.副吊钩9.门架结构刚柔耦合门式起重机模型:本文建立的门式起重机模型是刚柔耦合的虚拟样机模型。采用Pro/E建立门式起重机中的刚体部件，然后导入ADAMS。柔性体部件的建立(见图2)是采用ANSYS有限元分析计算软件，有限元模型一般采用建立刚性节点的方法实现柔性体与外部构件之间的连接。如图2中的10个外部节点即为在ADAMS中的联接点,其中，1〜4处装配大车车轮;5〜9处装配小车系统;7、10处装配副机房。ANSYS中的柔性体通过mnf.中性文件导入ADAMS中。图2门架结构的有限元模型钢丝绳是起重机的主要挠性构件。钢丝绳在ADAMS中的建模方法是用一段段的圆柱体通过轴套力(Bushing)连接来模拟钢丝绳的。这种方法较真实地反映了钢丝绳的拉伸、弯曲等力学性能，但分段很多，计算规模较大。在很多情况下，钢丝绳弹性对系统影响很小或只有拉伸的影响，因此可以简单地用一个弹簧阻尼器来模拟。本文采用的就是这种方法。本模型中的钢丝绳由3段圆柱体结合Joint生成，上两段圆柱体之间施加了一个移动副，对其施加运动函数，用于模拟钢丝绳的起升运动;下两段圆柱体之间施加了一个弹簧力，赋予刚度和阻尼属性，模拟钢丝绳的动力学性能。当刚性体部件、柔性体部件和钢丝绳等都导入ADAMS后，添加相应的约束和载荷，就生成了门式起重机的整机模型，见图3,通过对整机进行仿真，即可输出仿真分析结果。

图3门式起重机虚拟样机模型起升工况的动力学仿真分析整个起升过程，在起升机构刚起动的瞬间，绳索系统是松弛的，因此起升机构运动的第1阶段是在空转，这时松弛的绳索被收紧;当绳索开始受力时,这一阶段结束，这时ml具有一定的速度，而m0、m2还处于静止状态。从绳索受力开始直至滑轮组的弹性张力等于荷载重力Q为止,这是运动的第2阶段。在这一阶段中，m2还处于静止状态，而ml和m0则在滑轮组弹性力和门架主梁弹性力的作用下产生振动，这一阶段钢丝绳的弹性力将从零增大到Q［1］。第3阶段从吊重离地瞬间开始［2、3］。它的计算模型简化为图4。图4中m0为小车系统的推算质量;k0为主梁的刚度系数；ml为电动机转子和机构中所有运动部分的推算质量；m2为吊重的质量;k为钢丝绳的刚度系数；P为电动机的驱动力；s0、s1、s2分别表示质量m0、ml、m2的位移。ml、m2的位移。//////.7?图4起升机构的计算模型由拉格朗日方程［4］可以得到m0、ml、m2的3个运动微分方求解该方程组可以计算出该三质量二自由度的完整解。但考虑到m0的动变形与静变形相差很少，故不考虑m0的运动方程，系统简化为二质量一自由度系统。在吊重离地的瞬间，ml的速度为v,m2的速度为零,可解得钢丝绳（及传动系统）的力FF=0+kf况-5))=kAsillt+8?一=(P-Q)-Q, +也2 -为： 一

离地起升工况中钢丝绳的受力情况可以用以下方程组表示：［| ,--..・,・妙■ * -F=ks（t） tC［Qftx）F=Q+kA'sin仙f+为+ (P-Q)tC[tift>)mi+州？F=Q 强伉“fl?2z t 厂F=Q+ ' (-Pz-Q)(\-co^)t)tt[4h)=Q 强3,R其中，［0,t1）为离地之前的时间；s（t）为离地之前钢丝绳随时间的伸长量；［t1,t2）为从离地开始到平稳上升之间的时间；［t2,t3）为平稳上升的时间；［t3,t4）为从制动开始到静止的时间。该方程组的曲线形式见图5。Q+U图5Q+U图5吊重离地起升工况的动载荷理论曲线2.1.2起升制动工况计算下面讨论起升机构在上升平稳后制动工况时系统的动载荷计算。起重机在这种工况时，受到起升机构中电动机的制动力Pz的作用。同样可解得钢丝绳（包括传动系统）的力F—Q k=Q- ~~为： '：-'<2.2起升工况的仿真分析[5]

Pz-Q)(1-cost^t)对于门式起重机的离地起升工况，在ADAMS所建虚拟样机中实现方法是，首先吊重静止放在地面上，地面与吊重之间施加接触力，初始接触力等于吊重重力，即钢丝绳尚未提供拉力，然后拉动起升钢丝绳上端将其提起，平稳上升后，对其制动。根据实际情况给钢丝绳上端施加描述起升运动的函数。2.2.1模型的相关参数设置根据电机选型后的起制动时间验算，得到启动过程为0.98s,制动过程为1.01s。起升运动函数的加速度曲线见图6。图6钢丝绳的起升加速度曲线钢丝绳的刚度则可以通过虎克定律计算得到。2.2.2起升工况仿真结果在仿真起始时首先做静平衡分析。然后模拟起升运动的运动副拉动钢丝绳上端将吊重拉离地面，得到平稳上升中的变形云图和应力云图，可知起升工况中门架的最大应力为129.06MPa。门架的变形曲线图和应力曲线图见图7、图8。

CtsneTranslstionaiDeiormefttonY-1-1-2-2-3-37门架振动的变形曲线CtsneTranslstionaiDeiormefttonY-1-1-2-2-3-35工8门架振动的最大应力曲线50E*006图9中的实线为吊重起升过程中钢丝绳的受力情况，虚线为起升瞬间地面对吊重的接触力变0E*006化曲线。图9钢丝绳张力与地面接触力曲线由图9可以看出，在1.38s时,吊重离开地面，此时钢丝绳承受较大的冲击力,其最大拉力为1.33x10*。在5s时，起升过程开始制动，此时钢丝绳的张力有一个较大的波动。

=一L33X10”_=]086=一L33X10”_=]086mg/2250XI0土X9.8/2数冷为：此外，通过仿真还可以得到机构中各部件的速度曲线、加速度曲线及相互之间的作用力等。3大车运行机构动力学仿真大车运行机构在启动和制动时，悬吊着的吊重将会摆动。吊重摆动时，将会对起重机产生附加的水平力。3.1大车运行工况计算理论3.1.1大车运行启动工况计算图10为大车运行机构运行时的示意图，图11所示的三质量二自由度的弹性系统为大车运行机图10大车运行启动工况示意图构的简化系统。图11大车运行启动工况简图计算模型中的参数说明如下：m0为起重机大车的推算质量；ml为大车传动系统的推算质量;m2为吊重的质量;k01为大车运行机构高速部分的刚度系数;l为吊重起升时钢丝绳的长度;k20为吊重摆动的刚度系数,k20=Q/l=m2g/l;P为大车运行机构的电动机的驱动力;W为大车的运行静阻力；s0、s1、s2分别为质量m0、ml、m2的位移。由于高速部分的刚性系数k01比吊重摆动的刚性系数k20大得多。因此可以用一个二质量一自由度的系统来解决吊重的摆动问题，见图12。图12大车运行启动工况计算简图

由拉格朗日方程得到运动微分方程Co+叫iJ7+服)，为■筮，=p-w&-dSm工2~^20『•免■力/=°组:Idr已知s0-s2为吊重与大车之间的相对摆幅，因此可以解得摆动的水平力为:式!(]Nq一+用“相3=I中，3.1.2大车运行制动工况计算门式起重机在大车运行制动时，受到运行机构中电动机的制动力Pz的作用。此处亦将系统简化为二质量一自由度系统，见图13。图此处亦将系统简化为二质量一自由度系统，见图13。图13大车运行制动工况计算简图同样，解得摆动的水平力F=龙°<-%-S)为：ricoscot)大车运行工况中的水平动载荷,可以用以下方程组表TOC\o"1-5"\h\z\o"CurrentDocument"（P-W） 广f= ■co女te［如<f=o $em〃）\o"CurrentDocument"血（-P--W） 厂F= H-coswr;*-［.、&］\o"CurrentDocument"示人 5—的 一’其中，［0,t1）为大车运行起动的时间；［t1,t2）为大车平稳运行的时间；［t2,t3］为大车运行制动的时间。该方程组的曲线形式见图14。图14图14大车运行工况中的水平动载荷理论曲线3.2大车运行工况的仿真分析10.010.0在仿真起始时首先做静平衡分析。对大车和大车轨道之间施加移动副,对该移动副赋予加速度函数。根据电机选型后的起制动时间验算,得到启动过程为3.6s,制动过程为4.3s。由大车运行的变形云图和应力云图可知大车运行过程中门架的最大应力为129.35MPa。门架的竖直方向的变形曲线见图15，图中实线表示前主梁的竖直方向变形，虚线表示后主梁“MAKKER.149.151TranslfltC'r^|_DefcrTTiatioriY3拍 七岫牝— 的竖直方向的变形。图15大车运行的竖直方向变形曲线门架水平方向的变形曲线见图16，由图可知前后主梁的水平方向变形一致。

MARKERjl4导J MARKERjl4导J TrarmtabonwDefDriTrttHjn心- 每00~-01TimiP(f-rc1图16大车运行的水平方向变形曲线大车运行工况的最大应力曲线见图17，图中实线表示前主梁的最大应力曲线，虚线表示后主Dooocooooon-r897fi-c.-4331fl_?才『痴21h2・Z2231wiDooocooooon-r897fi-c.-4331fl_?才『痴21h2・Z2231wi1_1—_114—<111梁的最大应力曲线。图17大车运行的最大应力曲线结论通过对门式起重机的起升工况和大车运行工况进行动力学仿真，得到了整机的动