基于曲柄连杆传动的舱口盖启闭装置多体动力学仿真.pdf

舱口盖启闭装置是保证舰船安全的重要分系统，其运动控制技术是实现船体水密和装甲防护等功能的关键。本文基于多体运动学和动力学理论分析，建立曲柄连杆结构舱口盖启闭装置的动力学仿真分析模型，进而结合各种海况下的经验方程对舱口盖进行受力分析，计算得到了装置中曲柄连杆机构各个部件在开关盖过程中的支反力、转矩以及输出功率等参数，为曲柄连杆传动模式的舱口盖启闭装置结构优化及电机选型提供理论依据。关键词：启闭装置；曲柄连杆机构；多体动力学中图分类号：U663.1文章编号：1 6 7 2-7 6 49(2 0 2 3)1 6-0 0 32-0 5Multi-body dynamics simulation of hatch cover opening and closing deviceAbstract:Hatch cover opening and closing device is an important subsystem to ensure the safety of ships,and of thisdevice,the motion control technology of opening and closing mechanism is the key to achieve the watertight and armor pro-tection of the ship.Based on the multi-body kinematics and dynamics theory analysis,this paper establishes the dynamicssimulation analysis model of the hatch cover opening and closing device with crank link structure,and then analyzes theforce on the hatch cover by combining the empirical equations under various sea conditions.The dynamic parameters such asthe reaction force,torque and output power of each component of the crank link mechanism in the device in the process ofopening and closing the hatch cover are calculated,It provides a theoretical basis for the structure optimization of the ware-house lid hoist and the motor selection of the crank connecting rod drive mode.Key words:the opening and closing device;crank-rocker linkage;multi-body dynamics0引言舱口盖启闭装置是舰船装备的主要组成部分1-2 。舱口盖用于弹库的水密防护，并为筒（箱）弹提供装甲防护，启闭装置功能则是依据导弹的状态执行舱口盖的开启到位或关闭、锁定等指令动作。此外，在导弹意外点火时，启闭装置需能应急打开舱口盖，为意外点火导弹提供飞离通道。因此，舱口盖启闭装置的运动控制对导弹成功发射和舰船安全性至关重要。目前，启闭装置的传动模式主要分为液压传动3、螺杆传动4 以及曲柄连杆传动5-6 。液压传动系统受季节气候及海况影响较大，在极端情况下液压系统存在脉冲冲击、液压泄漏等不利因素，影响舱口盖系统的收稿日期：2 0 2 2-1 2-2 5作者简介：李红（1 9 9 0），女，博士，工程师，研究方向为水面发射装置结构设计。文献标识码：Abased on crank-linkage transmissionLI Hong,LIU Jun-zhou,WANG Bei,HE Chao-xun,CHEN Jin-cao(The 713 Research Institute of CSSC,Zhengzhou 450015,China)doi:10.3404/j.issn.1672-7649.2023.16.007安全性。螺杆传动装置容易受到腐蚀。相比而言，曲柄连杆结构作为一种常见的机械式传动机构，可以方便实现复杂运动，传递较大的动力，同时便于加工生产成本低，可用于舱口盖的启闭装置中7-8 。殷炜棋等9 采用Simscape对曲柄连杆机构运动轨迹进行了仿真研究，张林仙等1 0-1 2 对内燃机曲柄连杆机构进行了动力学分析。然而，海上作业工况包括了舰艇的摇摆、风的作用等不确定因素的影响，作用于舱口盖的负载转矩随开启角度的变化而变化，导致曲柄连杆启闭机构的变化规律难以描述1 3-1 4。因此，本文首先对舱口盖在不同工况下的受力进行分析，并将多种工况下的受力进行叠加，通过舱口盖启闭装置的第45卷三维实体仿真模型来分析得到该多体系统的极限运动过程，从而获得了开关盖过程中启闭机构各构件的运动参数极限值。利用机构的运动参数极限值可以为舱口盖启闭装置性能优化设计以及机构电机选型提供理论指导1 5。1多体系统运动及动力学分析1.1运动学分析假设某多体系统由n个构件组成，机械系统的仿真模型将构件与构件采用运动副联接，并采用广义坐标为代数方程，则其运动学约束方程可表示为：K()=(q),(),.,(g)T=0假设系统中所有运动副的自由度为m，为使系统具有确定运动，要使系统实际自由度为0，需要为系统施加m-n个驱动约束方程。同时，驱动约束通常认为是系统广义坐标系和时间的函数，表示为：P(q,)=P(q,1),d(q,t),.-(q,)T=.由于驱动约束在其集合内部及其与运动学约束合集中是独立和相容的，此时，驱动系统将作确定运动。由式（1）系统运动学约束和式（2）驱动约束方程组合成系统所受的全部约束，表示为：ak(a.=0.d(q,t)=P(q,t)式（3）为m个广义坐标的非线性方程组，其构成了系统的位置方程。对式（3）求导可得速度约束方程为：d(q,9,t)=g(q,t)g+(q,t)=g(q,)g-V=0。进一步对速度约束方程（4）求导可得加速度方程：d(q,q,q,t)=q(q,t)q+(q(q,t)gq+2gr(q,t)q+(q,t)=0。式中：为雅可比矩阵，如果 的维数为m，q 的维数为n，么q为mxn矩阵。由式（1）式（5）可知通过求解系统的约束方程，可获得零自由度系统的位置、速度、加速度。1.2动力学分析对于刚体系统，假设刚体B用质心笛卡尔坐标和反映刚体方位的欧拉角作为广义坐标可表示为：q=xyz4e$李红，等：基于曲柄连杆传动的舱口盖启闭装置多体动力学仿真singsing0B=singsinp0-singcosO10(1)则构件在知心笛卡尔坐标系下的角速度可表达为：W=B。(10)系统在质心笛卡尔坐标系下的的动能可表示为：RTMR+TB.JBy.T=2(2)式中：M为构件的质量矩阵；J为构件在质心坐标系下的惯量阵。考虑约束方程的动力学方程为：d(oT)nQ+i=1(3)式中：T为系统广义坐标表达的动能；Q为广义坐标下的广义力，最后一项为约束方程的雅可比矩阵与拉格朗日乘子，其表达了广义坐标的约束反力。进一步引入广义动量：8T(4)P=因此，质心笛卡尔坐标系下的移动方向和转动方向动量分别表示为：(5)PR0qR)OTBTJB。设广义力为：Q=HTF。式中，H为外力的坐标转换矩阵。(6)F=f(q,u,t)。33令：R=x y z,=$T。定义一个欧拉坐标系，该坐标系的3个单位矢量分别为3个欧拉转动的轴，该坐标系到刚体B质心坐标的转换矩阵为：(9)(11)(12)(13)8TMR,(7)(8)(14)(15)(16)(17)34因此可建立质心笛卡尔坐标系下的动力学方程-微分-代数方程：oTP+QTa+HTF=0.可知，当外力已知条件下，联合约束方程和微分-代数方程便可求解多体系统的相关运动参数。1.3仿真模型考虑到冲击惯性的影响，启闭装置采用实体建模，模型所用材料均采用低合金高强度结构钢，密度为7 8 7 0 kg/m3，利用ADAMS动力学仿真软件添加相应的运动副。为便于说明机构间的相互关系以及后续的受力分析，将实体模型简化为二维结构原理图（见图1）。图中，AB为曲柄，BC为连杆，CD为摇杆，DE为舱口盖。舱口盖通过转臂与摇杆固连在同一转轴上。其中，运动副A为主动曲柄与电机轴之间的转动副，数量为1；运动副B、C 为主动曲柄与连杆、连杆与摇杆之间的转动副，左右各1 个；运动副D为转臂轴与盖体之间的转动副，数量为1；E为舱口盖。E口图1 舱口盖启闭装置二维结构原理图Fig.1 Schematic diagram of two dimensional structure of hatchcover hoist1.4模型受力分析启闭机构开关盖过程中，受到舰艇的摇摆、风向作用等多种不确定因素的影响。另外，作用于舱口盖的负载转矩随开启角度变化而变化。因此，模型受力分析时兼顾不同的受力工况，同时进行如下简化：1）假设作用在舱口盖上的过载系数取模型改进舰在6 级海况下的最大摇摆载荷。2）假设盖体在开关过程中除了开盖结构传递的动力外，其他外力包括舰艇的摇摆载荷、风阻等阻力。此情况下，舱口盖受力最为恶劣，从而可确保盖体解锁启闭机构在6 级海况及任何风向时均正常工作。舰船科学技术3）假设舰船航行所产生的风始终平行于舰甲板，相对风速为舰的航速；自然风作为阻力始终垂直作用于舱口盖，风速为6 级海况下的自然风速。(18)根据上述假设，开盖与关盖过程舱口盖在其局部笛卡尔坐标系下的受力图如图2 所示。HW2FwlFL(a)开盖过程F2FwFw2D图2 复杂海况下舱口盖受理分析图Fig.2Analysis of hatch cover acceptance undercomplex sea conditions图中，L为舱口盖转轴到质心距离；为舱口盖开启角度；Fx为x方向摇摆过载力；Fz1、F z 2 为z方向摇摆过载力；Fw1为舰船航速引起的风载；Fw2为自然风引起的风载。(19)BFz=nzmgFw1(2)=CiknC2v1(2)A。式中：nx和nz分别为x、z 方向摇摆过载系数；Ci、C2为风速（气动阻力）系数；kn为风压高度系数；A为迎风面积；V1(2)为风速；V1为舰船航行时产生的相对风速；V2为6 级海况下的风速。将上述各个方向局部坐标系下的力进行向量运算，得到在系统基本坐标系下外力矩阵为：HTF=Fx Fy Fz JT,式中，Fx、F y、F z 分别为启闭机构在其整体迪卡尔坐标系下X、Y、Z三方向的合力。在笛卡尔坐标系下，将式（2 2）代人式（1 8）即第45卷甲板面FXL甲板面(b)关盖过程Fx=nxmg(20)(21)(22)第45卷可求得机构中各个部件的相关动力学参数。1.5结果及讨论舱口盖开启过程瞬时角度和角速度如图3所示。可以看出，舱口盖开启过程是非匀速进行的，整个过程可分为解锁-开启-锁死3个阶段。舱口盖开启最大角度为9 7.0 6，开盖总时间约为2.4s。100+舱口盖开启角度开启角速度7550$25-00.0图3舱口盖开启角度与角速度-时间曲线Fig.3 Opening angle and angular velocity curve of hatch cover atdifferent times图4为曲柄-摇杆传动比-时间曲线。可以看出，该连杆机构在解锁和锁死过程的传动比趋于无穷大值，而传动比最小值则存在于开启过程。整个开盖过程中，舱口盖、曲柄及摇杆主要运动参数汇总如表1 所示。另外，由于受力差异，开关盖过程中舱口盖负载转矩也存在明显差距。图5表明，开盖过程中相同时100一一曲柄-摇杆传动比80%/160402000.0图4曲柄-摇杆传动比-时间曲线Fig.4 Transmission ratio curve of crank rocker at different times表1 各系数及其数值汇总表Tab.1Summary of coefficients and their values参数nx数值0.530210.65801.20.6131.00.9063注：Ao为舱口盖面积，nz“+”表示摇摆载荷垂直于甲板并指向下方；“二”表示摇摆载荷垂直于甲板并指向上方。李红，等：基于曲柄连杆传动的舱口盖启闭装置多体动力学仿真8060402000.51.0时间/s0.51.0时间/snzCiC2kAo/m2V/ms-1v2/ms-114.424.4352000开盖过程1600关盖过程120080040000.0图5开关盖过程中舱口盖负载转矩-时间变化曲线Fig.5 Load torque curve of hatch cover at different times刻盖体负载力矩明显高于关盖过程。0.7 5s时，开盖过程盖体负载力矩达到最大值1 7 51 Nm，而关盖过程盖体负载力矩则在0.7 0 s时达到最大值8 56 Nm。机构中1.52.01.52.00.52.54组铰链归一化后的支反力随时间的变化曲线如图6所示。可以看出，0.2 5s时，铰链A、B（C）支反力达到最大值；0.2 42 s时，铰链D支反力达到最大值。图7 和图8 分别为舱口盖开盖和关盖过程曲柄转矩（归一化）和曲柄功率-角度变化曲线。可以看出，开盖过程中曲柄转矩与瞬时功率远高于关盖过程；开盖角度为34.7 4时，曲柄转矩和瞬时功率达到最大值。关盖角度为7 1.48 时，曲柄转矩和瞬时功率达到最大值。开关盖过程中，舱口盖功率最大值及其对应时间汇总如表2 所示。54260.53关盖过程2.500.0图6 开关盖过程中舱口盖启闭机构中4组铰链支反力（归一化）-时间变化曲线Fig.6 Support reaction(normalized)curve of four groups of hinges in hatch cover hoisting mechanism at different times2结语1）通过多体系统的运动学和动力学理论分析可知，对舱口盖启闭机构在开关过程的受力分析是保证1.0时间/s开盖过程+铰链A-铰链B/C铰链D1.01.50.51.0时间/s1.51.52.02.02.02.52.52.53620001500100050000图7 开关盖过程曲柄转矩（归一化）随舱口盖角度Fig.7 Curve of crank torque(normalized)changing withhatch opening angle1.50.50.00图8 开关盖过程曲柄功率随舱口盖角度变化曲线Fig.8Curve of crank power changing with hatch opening angle表2 开关盖过程中曲柄连杆结构相关动力学参数汇总Tab.2 Dynamic paramaters at the limit position of the hatch coverhoisiting mechanism during the process from opening and项目/参数最大值对应时间/s开启角度/()97.06角速度/():s-58.78瞬时传动比80瞬时功率/W1476仿真计算结果准确的前提条件。2）舱口盖启闭装置仿真计算结果表明，舱口盖启闭装置用电机的最大转矩与功率都出现于开盖过程中，其最大功率值可作为电机选型的依据。3）为确保舱口盖启闭装置可以在极端条件下正常运行，在对舱口盖受力分析时将各种可能存在的载荷舰船科学技术取最大值并线性叠加，因此电机最大功率计算结果比开盖过程实际偏高。关盖过程参考文献：1 ITO K,KUROIWA T,AKAGI S.Design optimization ofhatch cover opening and closing mechanism of ships.Applica-tion of minimax non-linear optimization methodJj.Transac-福tions of the Japan Society of Mechanical Engineers,1985,一51(462):460468.2040角度/)）变化曲线开盖过程关盖过程12040角度/）closing最小值对应时间/s2.400.221.3000.07或2.401.331.017882第45卷608060801001000.070.07或2.381.240.617【2 王忠强，秦荣俊,支山亚.基于质量优化目标的舱口盖结构设计及强度校核J.船舶,2 0 2 2,(1):7 0-7 7.3曲志云.9 2 2 型打捞救生船的船舱口盖启闭装置.船舶，1992,(2):41-42.【4吴韵侠，胡修池，高卫红.引黄涵闸螺杆式启闭机的维护技术研究与应用J.水利水电技术,2 0 1 1,42(2):7 4-7 5.5 BATES K T,FALKINGHAM P L.Estimating maximum biteperformance in Tyrannosaurus rex using multi-body dynamics.Biology Letters,2012,8(4):660-664.【6】吴楠，廖日东，张保成，等.柴油机曲柄连杆机构多体动力学仿真分析J.内燃机工程,2 0 0 5,2 6(5)：6 9-7 3.【7】刘永亮，任克亮，马旭轮，等.平面四杆机构在舰载垂直发射装置舱口盖系统中的应用分析.装备环境工程，2 0 1 9,16(5):1317.8 SINGH R,CHAUDHARY H,SINGH A K.Defect-free opti-mal synthesis of crank-rocker linkage using nature-inspired op-timization algorithmsJ.Mechanism and Machine Theory,2017,116:105122.9殷炜棋，杨鲁义，朱国仁.基于Simscape的曲柄连杆机构运动轨迹仿真研究.工程与试验,2 0 2 2,6 2(2):1 8-2 1.10 张林仙，张生保.基于多体系统仿真的内燃机曲柄连杆机构动力学分析.装备制造技术,2 0 0 6(4):1 4-1 6.11 杜爱民，杜玉彪，楼狄明，等.发动机曲柄连杆机构动力学仿真分析平台研究J.计算机仿真，2 0 1 0,2 7(1 1):2 7 8-2 8 3.12叶奇，孙鹏文，张敏.发动机曲柄连杆机构多体动力学建模的若干问题J.机电工程，2 0 0 7,2 4(1 2):53-55.13 享郭亦平，杜春旺，李明，等.恶劣海况下船舶航向控制仿真及应用研究.舰船科学技术,2 0 0 8,30(1):1 42-1 52.GUO Y P,DU C W,LI M,et al.Simulation and research ofapplication of ship yaw control in the bad sea conditionJ.Ship Science and Technology,2008,30(1):142-152.14 LEPIDI M,GATTULLI V.Non-linear interactions in the flexi-ble multi-body dynamics of cable-supported bridge cross-sec-tions JJ.International Journal of Non-Linear Mechanics,2016(80):1428.15潘爱华.起重机起升机构电动机的选用.起重运输机械，2006(7):31-34.