Loading [MathJax]/jax/output/HTML-CSS/jax.js
  舰船科学技术  2025, Vol. 47 Issue (2): 101-106    DOI: 10.3404/j.issn.1672-7649.2025.02.017   PDF    
一种三自由度稳定廊桥试验台动力学分析与试验研究
邱建超1, 曹彪1, 唐祖定1, 方一格1, 王生海2, 魏昌淼1     
1. 中国船舶集团有限公司第七〇三研究所无锡分部,江苏 无锡 214151;
2. 大连海事大学 轮机工程学院,辽宁 大连 116026
摘要: 海上廊桥在作业时由于受到风、浪、流等载荷的扰动会产生较大摇摆,极大地降低了作业安全。针对该问题设计一种具有运动补偿能力的三自由度海上廊桥试验台。文章首先建立不规则海浪模型以及海上廊桥的运动学数学模型,基于朗格朗日法建立了该试验台的动力学模型,应用Adams对试验台进行动力学仿真与分析,并搭建试验台进行试验研究,试验结果验证了三自由度稳定廊桥试验台设计合理具有良好的运动补偿效果。
关键词: 海上廊桥     试验台     运动学分析     运动补偿    
Kinetics analysis and experimental research of a three-degree-of-freedom stabilized offshore gangway
QIU Jianchao1, CAO Biao1, TANG Zuding1, FANG Yige1, WANG Shenghai2, WEI Changmiao1     
1. Wuxi Division of the 703 Research Institute of CSSC, Wuxi 214151, China;
2. College of Marine Engineering, Dalian Maritime University, Dalian 116026, China
Abstract: During operation, the offshore gangway will sway greatly due to the disturbance of loads such as wind, waves, and currents, which greatly reduces the safety of the operation. Aiming at this problem, a three-degree-of-freedom offshore gangway test bench with motion compensation capability is designed. The article first establishes the irregular wave model and the kinematics mathematical model of the offshore gangway. A dynamic model of the experimental rig is established based on the Langrangian method, and Adams was used to simulate and analyze the dynamics of the test rig, and set up a test bench for experimental research. The experimental results verify that the design of the three-degree-of-freedom stable gallery bridge test bench is reasonable and has a good motion compensation effect.
Key words: offshore gangway     test platform     kinematic analysis     motion compensation    
0 引 言

海上廊桥主要用于海洋平台与船舶之间人员交换、货物搬运以及维修平台设备等,是一条海上通道。由于海上复杂多变的环境,使得海上廊桥在作业时不可避免的产生较大摇摆,因此研究具有运动补偿能力海上廊桥显得尤为重要。Dai等[1]研制出一种六自由度海上廊桥,并将海上廊桥分为稳定平台以及舷梯两部分,稳定平台由6个液压缸组成,舷梯可进行一定范围的伸缩运动。在稳定平台上装有检测船舶实时运动的传感器,将检测到的物理量转换为电信号传输给控制系统。控制系统根据接收的信号计算出需要补偿的位移,随后控制系统输出一个与船舶运动完全相反的控制信号给平台的6个液压缸,驱动液压缸反方向动作使廊桥整体保持稳定,从而实现运动补偿,这是海上廊桥的雏形。吴彦岐[2]发现虽然stewart六自由度稳定平台可以实现扰动全补偿,但是在应用中仍然面临着设计及控制方面的难题,其首先建立了稳定平台的结构,之后建立了包括位姿正、反解在内的运动学模型,并建立了稳定平台的单刚体动力学模型。杜黎童等[3]利用牛顿—欧拉方程建立了火炮并联平台的动力学模型,基于逆动力学建立轨迹跟踪模型,并借助Matlab和Adams联合仿真对该模型的正确性进行验证。马小腾[4]利用全雅可比及弹性模型得到并联平台的刚度模型,基于此模型对平台姿态误差进行预测,最终通过仿真验证其有效性。许猛[5]基于自抗扰原理搭建线性观测器来对平台的动态不确定性进行预测。郭明轩[6]在继承stewart平台动力学理论基础上,得到平台在运动情况下,基于Newton-Euler矢量力学法的稳定平台动力学方程,求解驱动连杆质量惯性矩带来的连接铰约束力。

本文设计一种三自由度海上廊桥试验台,对其进行运动学及动力学分析,并搭建了试验台,通过试验验证了模型的准确性及试验台的可靠性,为具有运动补偿功能的海上廊桥的设计提供了可行的方案[79]

1 不规则海浪及海上廊桥运动学建模 1.1 不规则海浪运动学建模及仿真

船舶纵摇比横摇的幅值小许多,所以本文以横摇作为激励进行仿真分析[10]。选择5级海况进行仿真研究,5级海况频率增量、频率及有义波高分别为:0.08、0.98及3.2m。其余参数以“育鲲”轮为参考对象,其参数如表1所示。

表 1 “育鲲”轮船舶参数 Tab.1 Parameters of the ship "Yukun"

船舶横摇波倾角公式为:

ae(t)=KBKTω2ngNn=12S(ω2n)Δωcos(ωent+εn)sinμ (1)

式中:ωn为规则波频率;g为重力加速度;ωe为某规则波频率;εn为相位;KBKT为不规则波相位角[11]。分别在波倾角为45°(首斜浪)、90°(正横浪)、135°(尾斜浪)下进行计算分析。用Matlab对船舶横摇波倾角公式进行仿真,仿真曲线如图1所示。

图 1 45°、90°、135°的浪向角 Fig. 1 45°, 90°, 135° wave direction angle
1.2 海上廊桥运动学建模及仿真

该稳定平台为3DOF平台,可实现对船舶横摇、纵摇、垂荡3个自由度运动的补偿[12]

其空间几何模型如图2所示。其中x0y0z0为大地坐标系、x1y1z1是以摇摆平台质心建立坐标系、x2y2z2是下平台的质心建立的坐标系、x3y3z3是以上平台质心建立的坐标系[1314]

图 2 海上廊桥试验台数学模型 Fig. 2 Mathematical model of the test bed of the offshore gangway

根据机器人学[11]的相关理论可得到坐标系On+1相对于坐标系On的旋转矩阵为:

n+1nR=Rx(θn+1x)Ry(θn+1y) (2)

式中:θn+1xθn+1y分别为坐标系On+1相对于坐标系Onx轴、y轴方向旋转的角度。因为上平台的位置取决于摇摆平台的变幅角度,利用坐标变换将T3坐标系在T0坐标系表示。T1T2T3x轴的旋转角度为θ1xθ2xθ3x的齐次变换矩阵为:

01Tx=[10000cosθ1xsinθ1x00sinθ1xcosθ1xk50001] (3)
12T=[100k4010000100001] (4)
23Tx=[10000cosθ2xsinθ2x00sinθ2xcosθ2xk30001] (5)

由式(3)~式(5)可得船舶横摇的坐标变换公式为:

0P=01T12T23T3P (6)

式中:以向量的形式对点的空间位置坐标进行定义,所以0P3P分别为点在T0T3坐标系中的坐标。根据式(3)~式(6),可得到点A1A2A3和点B1B2B3x轴旋转的坐标0pA1x0pA2x0pA3x0pB1x0pB2x0pB3x,进而得到1号、2号、3号电动缸绕x轴转动的长度为L1xL2xL3x的表达式:

L1x=(2)(2k213k1k2+2k22+2k23k1k2cosθ2x2k1k3sinθ2x)/2 (7)
L2x=(2)(2k213k1k2+2k22+2k23+k1k2cosθ2x+2k1k3sinθ2x)/2 (8)
L3x=(k122cosθ2xk1k2+2sinθ2xk1k3+k22+k32) (9)

由波倾角仿真曲线可以得出在波倾角为90°时浪向角最大,因此本文以90°浪向角及5级海况为横摇激励利用Matlab[15]对式(7)~式(9)进行仿真,得到3个电动缸位移随时间变化的曲线如图3所示。

图 3 横摇工况电动缸位移 Fig. 3 Electric cylinder stroke under rolling condition

可以得到电动缸行程振幅较小且较为平缓,合理设计三自由度稳定廊桥试验台可使得试验台做三自由度的复合运动[16]。且电动缸的最大行程为217 mm,为接下来试验台结构设计提供了技术指标。

2 试验台动力学建模及分析 2.1 试验台拉格朗日动力学模型

根据拉格朗日方程[17]

ddt(M(Q,˙Q)˙Q)M(Q,˙Q)Q=F

式中:M(Q,˙Q)=DSDS分别为系统的总动能及总势能;F为广义力。

2.1.1 上平台的动能势能

上平台只能沿竖直运动,所以上平台只有沿竖直方向运动的动能及转动动能。其动能为:

Ed=12mzVZ2 (10)

式中:VzZ方向速度;mz为其质量。其转动时的动能为:

Es=12ωhTIcωh (11)

式中:wh为角速度,ωh=[100cosα]Ic为上平台在动坐标系中的转动惯量矩阵。

综合上述两式可得上平台的动能为:

E1=12˙QTM1(Q)˙Q (12)

式中:

M1(Q)=[mz000Ix000Iycos2α+Izcos2α]

上平台的势能为:

N1=mzgzz (13)
2.1.2 三个伺服电动缸的动能势能

根据文献[18]得3个电动缸的动能势能。3个电动缸的动能为:

Edo=12˙QTMdo(Q)˙Q (14)

式中:m1m2分别为电动缸和杆的质量。

Mdo(Q)=(m1+m2)(JT2(DJT1KJ1)J2)
J1=diag(uT1,uT2,uT3)

式中:ui=BiAili代表第i个连杆的方向。

J2=[I3RzRyP(i)RsxA1RzP(j)RyRsxA1I2RzRyP(i)RsxA2RzP(j)RyRsxA2I1RzRyP(i)RsxA3RzP(j)RyRsxA3]
P(i)=[000001010]P(j)=[001000100]

I3[001]TE为对角矩阵。

K=diag(k1,k1,k1,k2,k2,k2,k3,k3,k3)
ki=[m1l1im2l2ili(m1+m2)][m1l1im2l2ili(m1+m2) + m2m1+m2]

3个伺服电动缸的势能为:

Nsh=(m1+m2)g3i=1[m1l1im2l2ili(m1+m2)+m2m1+m2]×(zpxBsinβ+yBcosβ+zBcosβ) (15)

式中:zp为其质心竖直方向位移;xByBzB为其各铰点在动坐标系下的坐标。

2.1.3 伸缩立柱(从动支链)的动能和势能

因为从动支链只有竖直方向位移,所以其势能为:

Nc=mcg(zcl0) (16)

式中:zc为竖直方向位移;l0为从动支链顶部到其质心的距离。

2.2 试验台的广义力

广义力F[17]X=[zαβ]T上的虚拟力和力矩,设试验台的广义力为F=[F1F2F3]T,试验台各从动支链的力为f=[f1f2f3]T

从动支链的速度与上平台之间的关系[18]

l=JX

式中:X=[zαβ]J为一阶系数,J=J1J2

˙l=J˙X,则虚位移为σl=Jσx。其中,σlσx分别为从动支链、上平台的虚位移矢量。

各驱动支链所做的虚功为:

w1=fTσl=f1σl1+f2σl2+f3σl3 (17)

式中:f为驱动力向量。

运动平台所做的虚功为:

w2=FTσX=Fxσz+Fασα+Fβσβ (18)

根据虚功原理可得:

FTσX=fTσlFTσX=fTJσXF=JTf (19)
2.3 试验台动力学模型

根据所得动能与势能结合拉格朗日方程:

M1¨X+Vm˙X+G=F+Fext (20)

式中:M1为广义惯性力,M1=M1(Q)+Mdo(Q)Vn为广义哥氏力及离心力,Vn=Vnh + VdoG为广义重力,G=Gh+Gdo

Gh=[mz00]T (21)
Gdo=(m1+m2)g3i=1[(m1l1m2l2(m1+m2)l2i)liX(k)]sAiz+(m1+m2)g3i=1[(m1l1m2l2(m1+m2)li)+m2m1+m2]sAizX(k) (22)

式中:Fext为舷梯和外界作用在平台上的力和力矩;VnhVdo分别为主、从动支链产生的哥氏力和离心力。

根据文献[19]可将Vmh求得:

Vnh=[0000Vms22Vms230Vms32Vms33] (23)

式中:Vms22=0Vms23=12(IyIz)sin(2α)˙βVms32=0Vms33=(IYIZ)sin(2α)˙αMdo(Q) = (m1+m2)×(JT2(EJT1KJ1)J2)

2.4 试验台动力学仿真分析

将在Solidworks中建立的三维模型导入到Adams中,之后对其进行几何前处理、添加质量属性、添加运动副以及驱动等操作,得到海上廊桥的试验台模型。将Matlab仿真得到数据进行处理,导入Adams中进行动力学仿真,对其进行动力学仿真和曲线后处理,结果如图4所示。

图 4 海上廊桥位移、速度曲线 Fig. 4 Displacement and velocity curve of offshore gangway

可知,试验台与给定的横摇激励误差很小,验证了该模型的准确性,且位移最大为0.043 m,速度最大为0.017 m/s。其瞬时速度非常小,试验台整体具有很强的运动稳定性。

3 试验研究 3.1 海上廊桥试验台搭建

为了检验廊桥试验台设计的合理性,搭建试验台如图5所示。试验台的部分试验参数如表2所示。该试验台分为3DOF并联稳定平台和伸缩舷梯两部分。稳定平台由上、下平台及3个电动缸构成。能够进行3DOF(横摇、纵摇及升沉)运动。下平台与模拟平台固连,试验时与模拟平台位姿一致。上平台与伸缩舷梯固定支架固连,实现变幅动作。上、下平台安装角度传感器对平台及模拟平台的角度进行实时检测。下平台上安装有伸缩立柱,防止试验台塌陷。

图 5 海上廊桥试验台 Fig. 5 Offshore gangway test bed

表 2 部分试验设备参数 Tab.2 Some test equipment parameters
3.2 海上廊桥试验台运动控制系统设计

本文的测控系统整体方案如图6所示。

图 6 测控系统整体方案 Fig. 6 Measurement and control system overall scheme

该测控系统利用工控机将硬件与虚拟仪器连接,基于LabVIEW编写上位机程序,该试验台数采系统使用数据采集卡(USB-4704)对数据进行采集,使用运动控制卡(PCI-1285E)进行控制,通过选用具有RS485通信功能的伺服电机驱动器实现对电机转速、转矩等参数的实时测量。基于子VI模块对PCI-1285E进行编程,最终通过“PC+运动控制卡”完成试验台测控系统的搭建。

3.3 试验

本节试验内容分为两部分:第一部分测量试验台上、下平台的角度,分析该试验台减摇性能,试验测得的船舶运动模拟平台倾角与上平台的倾角曲线如图7所示。第二部分通过通讯协议(Modbus)读取电机实际转矩,将试验给定工况的角度和周期代入海上廊桥试验台数学模型,使用Adams对得到的执行机构运动参数进行仿真,再将仿真结果与试验数据进行对比。结果如图7图8所示。

图 7 横摇5°,周期20 s时倾角对比图 Fig. 7 Roll 5°, period 20 s angle comparison diagram

图 8 横摇5°、周期20 s时3个电动缸试验与仿真转矩 Fig. 8 Test and simulation torque

图7可知,下平台在该激励下,角度最大、最小幅值分别为5.8°、−6.11°;在该激励下上平台角度最大、最小幅值分别为2.54°、−2.65°。运动补偿下效果达到57%。由图8可知,在该激励下,1号电动缸试验测得电机转矩的最大、最小幅值分别为0.38、0.18 N·m;2号电动缸试验测得电机转矩的最大、最小幅值分别为0.59、0 N·m;3号电动缸试验测得电机转矩的最大幅值分别为0.25、0.11 N·m。从3个电动缸的最大值与最小值可知,即使在横摇为5°的情况下,3个电动缸的振幅并不大而且均是周期性变化,电动缸的受力环境较好,由试验测得的3个电动缸电机转矩数据和Matlab与Adams联合仿真得到的数据基本一致,误差处于范围内,证明仿真建模的正确性和三自由度稳定廊桥装置设计合理。

4 结 语

1)在平台随海况做相应运动时,电动缸也同时在运动,而且曲线的振幅较小且变化较为平缓,三自由度稳定廊桥试验台设计合理。

2)试验验证了本文设计的试验台具有良好的运动补偿效果,在横摇角度为5°时可达57%。

3)仿真得到的电动缸电机转矩与试验测得的转矩误差在20%以内,表明模型仿真结果的准确性。

4)在船舶运动模拟平台进行试验时,即使在横摇为5°的情况下,3个电动缸扭矩峰值分别为0.38、0.59、0.25 N·m,3个电动缸的扭矩最小值分别为0.18、0、0.11 N·m,3个电动缸振幅不大,而且均是周期性变化,电动缸的受力环境较好。

综上所述,本文设计的三自由度稳定廊桥试验台能够适应工况变化,具有良好的运动补偿效果,设计方案具有良好的可行性和实用性。

参考文献
[1]
DAI X, SONG S, XU W, et al. Modal space neural network compensation control for Gough-Stewart robot with uncertain load[J]. Neurocomputing, 2021, 449: 245-257. DOI:10.1016/j.neucom.2021.03.119
[2]
吴彦岐. 基于船舶运动预测的六自由度稳定平台[D]. 哈尔滨:哈尔滨工业大学, 2022.
[3]
杜黎童, 荆学东, 王旭飞, 等. 基于并联稳定平台的火炮跟踪仿真研究[J]. 火力与指挥控制, 2023, 48(2): 154-160. DOI:10.3969/j.issn.1002-0640.2023.02.024
[4]
马小腾. 舰载雷达稳定平台构型设计及结构优化与误差补偿研究[D]. 秦皇岛: 燕山大学, 2022.
[5]
许猛. 并联三自由度船载稳定平台稳定控制研究[D]. 大连: 大连海事大学, 2020.
[6]
郭明轩. 海上廊桥试验台的设计与减摇试验研究[D]. 大连: 大连海事大学, 2021.
[7]
戴元风. 舰载并联稳定平台的姿态控制研究[D]. 武汉: 华中科技大学, 2021.
[8]
奚瑞. 基于三自由度稳定平台的海上廊桥结构设计研究[D]. 大连: 大连海事大学, 2020.
[9]
吉阳, 王虎, 陈海泉, 等. 船用起重机减摇装置液压系统设计与试验研究[J]. 合肥工业大学学报(自然科学版), 2019, 42(8): 1041-1046+1076.
[10]
邱建超, 牛安琪, 仇伟晗, 等. 3_UPU/PU_RRP海上稳定廊桥运动学建模与仿真[J]. 舰船科学技术, 2023, 45(1): 78-82. DOI:10.3404/j.issn.1672-7649.2023.01.014
[11]
JOHN J C, 贠超. 机器人学导论[M]. 北京: 机械工业出版社, 2006.
[12]
QIU J, QIU W, NIU A, et al. Modeling and analysis of offshore gangway under dynamic load[J]. Journal of Marine Science and Engineering, 2023, 11(1): 77. DOI:10.3390/jmse11010077
[13]
NIU A, WANG S, SUN Y, et al. Dynamic modeling and analysis of a novel offshore gangway with 3UPU/UP-RRP series-parallel hybrid structure[J]. Ocean Engineering, 2022, 266: 113122. DOI:10.1016/j.oceaneng.2022.113122
[14]
顾永凤, 邱广庭. 基于ADAMS的船用浮式起重机补偿装置平台运动学分析与研究[J]. 广州航海学院学报, 2018, 26(3): 9-13+23. DOI:10.3969/j.issn.1009-8526.2018.03.003
[15]
邱建超, 陈海泉, 仇伟晗, 等. 船载串并混联海上稳定廊桥动力学建模与分析[J]. 中国造船, 2023, 64(1): 146-160. DOI:10.3969/j.issn.1000-4882.2023.01.013
[16]
孔令富, 张世辉, 肖文辉, 等. 基于牛顿—欧拉方法的6-PUS并联机构刚体动力学模型[J]. 机器人, 2004, 4(5): 395-399. DOI:10.3321/j.issn:1002-0446.2004.05.003
[17]
张青云, 赵新华, 刘凉, 等. 空间刚柔耦合并联机器人动力学建模及仿真[J]. 机械设计, 2020, 37(4): 61-66.
[18]
丁祝顺. 载荷稳定平台高精度稳定控制技术研究[D]. 长沙: 国防科技大学, 2021.
[19]
王磊. 基于减摇陀螺的海上引桥平台稳定性研究[D]. 镇江: 江苏科技大学, 2021.
一种三自由度稳定廊桥试验台动力学分析与试验研究
邱建超, 曹彪, 唐祖定, 方一格,