2. 大连海事大学 轮机工程学院,辽宁 大连 116026
2. College of Marine Engineering, Dalian Maritime University, Dalian 116026, China
海上廊桥主要用于海洋平台与船舶之间人员交换、货物搬运以及维修平台设备等,是一条海上通道。由于海上复杂多变的环境,使得海上廊桥在作业时不可避免的产生较大摇摆,因此研究具有运动补偿能力海上廊桥显得尤为重要。Dai等[1]研制出一种六自由度海上廊桥,并将海上廊桥分为稳定平台以及舷梯两部分,稳定平台由6个液压缸组成,舷梯可进行一定范围的伸缩运动。在稳定平台上装有检测船舶实时运动的传感器,将检测到的物理量转换为电信号传输给控制系统。控制系统根据接收的信号计算出需要补偿的位移,随后控制系统输出一个与船舶运动完全相反的控制信号给平台的6个液压缸,驱动液压缸反方向动作使廊桥整体保持稳定,从而实现运动补偿,这是海上廊桥的雏形。吴彦岐[2]发现虽然stewart六自由度稳定平台可以实现扰动全补偿,但是在应用中仍然面临着设计及控制方面的难题,其首先建立了稳定平台的结构,之后建立了包括位姿正、反解在内的运动学模型,并建立了稳定平台的单刚体动力学模型。杜黎童等[3]利用牛顿—欧拉方程建立了火炮并联平台的动力学模型,基于逆动力学建立轨迹跟踪模型,并借助Matlab和Adams联合仿真对该模型的正确性进行验证。马小腾[4]利用全雅可比及弹性模型得到并联平台的刚度模型,基于此模型对平台姿态误差进行预测,最终通过仿真验证其有效性。许猛[5]基于自抗扰原理搭建线性观测器来对平台的动态不确定性进行预测。郭明轩[6]在继承stewart平台动力学理论基础上,得到平台在运动情况下,基于Newton-Euler矢量力学法的稳定平台动力学方程,求解驱动连杆质量惯性矩带来的连接铰约束力。
本文设计一种三自由度海上廊桥试验台,对其进行运动学及动力学分析,并搭建了试验台,通过试验验证了模型的准确性及试验台的可靠性,为具有运动补偿功能的海上廊桥的设计提供了可行的方案[7 − 9]。
1 不规则海浪及海上廊桥运动学建模 1.1 不规则海浪运动学建模及仿真船舶纵摇比横摇的幅值小许多,所以本文以横摇作为激励进行仿真分析[10]。选择5级海况进行仿真研究,5级海况频率增量、频率及有义波高分别为:0.08、0.98及3.2m。其余参数以“育鲲”轮为参考对象,其参数如表1所示。
![]() |
表 1 “育鲲”轮船舶参数 Tab.1 Parameters of the ship "Yukun" |
船舶横摇波倾角公式为:
ae(t)=KBKTω2ngN∑n=1√2S(ω2n)Δωcos(ωent+εn)sinμ。 | (1) |
式中:
![]() |
图 1 45°、90°、135°的浪向角 Fig. 1 45°, 90°, 135° wave direction angle |
该稳定平台为3DOF平台,可实现对船舶横摇、纵摇、垂荡3个自由度运动的补偿[12]。
其空间几何模型如图2所示。其中
![]() |
图 2 海上廊桥试验台数学模型 Fig. 2 Mathematical model of the test bed of the offshore gangway |
根据机器人学[11]的相关理论可得到坐标系
n+1nR=Rx(θn+1x)Ry(θn+1y)。 | (2) |
式中:
01Tx=[10000cosθ1x−sinθ1x00sinθ1xcosθ1xk50001], | (3) |
12T=[100k4010000100001], | (4) |
23Tx=[10000cosθ2x−sinθ2x00sinθ2xcosθ2xk30001]。 | (5) |
由式(3)~式(5)可得船舶横摇的坐标变换公式为:
0P=01T12T23T3P。 | (6) |
式中:以向量的形式对点的空间位置坐标进行定义,所以
L1x=√(√2)(2k21−3k1k2+2k22+2k23−k1k2cosθ2x−2k1k3sinθ2x)/2, | (7) |
L2x=√(√2)(2k21−3k1k2+2k22+2k23+k1k2cosθ2x+2k1k3sinθ2x)/2, | (8) |
L3x=√(k12−2cosθ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。 |
式中:
上平台只能沿竖直运动,所以上平台只有沿竖直方向运动的动能及转动动能。其动能为:
Ed=12mzVZ2。 | (10) |
式中:
Es=12ωhTIcωh。 | (11) |
式中:
综合上述两式可得上平台的动能为:
E1=12˙QTM1(Q)˙Q。 | (12) |
式中:
M1(Q)=[mz000Ix000Iycos2α+Izcos2α]。 |
上平台的势能为:
N1=mzgzz。 | (13) |
根据文献[18]得3个电动缸的动能势能。3个电动缸的动能为:
Edo=12˙QTMdo(Q)˙Q。 | (14) |
式中:
Mdo(Q)=(m1+m2)(JT2(D−JT1KJ1)J2), |
J1=diag(uT1,uT2,uT3)。 |
式中:
J2=[I3RzRyP(i)RsxA1RzP(j)RyRsxA1I2RzRyP(i)RsxA2RzP(j)RyRsxA2I1RzRyP(i)RsxA3RzP(j)RyRsxA3]; |
P(i)=[00000−1010];P(j)=[001000−100]; |
K=diag(k1,k1,k1,k2,k2,k2,k3,k3,k3), |
ki=[m1l1i−m2l2ili(m1+m2)][m1l1i−m2l2ili(m1+m2) + m2m1+m2]。 |
3个伺服电动缸的势能为:
Nsh=(m1+m2)g3∑i=1[m1l1i−m2l2ili(m1+m2)+m2m1+m2]×(zp−xB′sinβ+yB′cosβ+zB′cosβ)。 | (15) |
式中:
因为从动支链只有竖直方向位移,所以其势能为:
Nc=mcg(zc−l0)。 | (16) |
式中:
广义力
从动支链的速度与上平台之间的关系[18]:
l=JX。 |
式中:
由
各驱动支链所做的虚功为:
w1=fT⋅σl=f1⋅σl1+f2⋅σl2+f3⋅σl3。 | (17) |
式中:
运动平台所做的虚功为:
w2=FT⋅σX=Fx⋅σz+Fα⋅σα+Fβ⋅σβ。 | (18) |
根据虚功原理可得:
FT⋅σX=fT⋅σl,FT⋅σX=fT⋅JσX,F=JT⋅f。 | (19) |
根据所得动能与势能结合拉格朗日方程:
M1¨X+Vm˙X+G=F+Fext。 | (20) |
式中:
Gh=[mz00]T, | (21) |
Gdo=(m1+m2)g3∑i=1[−(m1l1−m2l2(m1+m2)l2i)∂li∂X(k)]sAiz+(m1+m2)g3∑i=1[−(m1l1−m2l2(m1+m2)li)+m2m1+m2]∂sAiz∂X(k)。 | (22) |
式中:
根据文献[19]可将
Vnh=[0000Vms22Vms230Vms32Vms33]。 | (23) |
式中:
将在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 |
本文的测控系统整体方案如图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.
|