摘要
对岸桥前大梁的结构参数和运行参数进行了同步优化。首先,建立前大梁系统简化力学模型,对比整机有限元模型与简化模型的有限元模型得到前3阶模态频率,验证了简化模型的有效性;其次,为了提高优化过程的计算效率,利用假设模态法及拉格朗日法得到简化模型的动力学方程;然后,以Runge‑Kutta法计算得到简化模型的吊重摆角最大值、岸桥前大梁中点处的最大挠度及小车总运行时间为优化目标,采用多目标遗传算法对岸桥前大梁系统的结构参数及小车运行参数进行同步优化,得到Pareto解集;最后,根据系统性能的具体需求从解集中选择合理的解,指导岸桥前大梁系统的结构设计以及小车运行参数的设定,为岸桥防摇提供了一个新的思路。
随着集装箱运输船舶向大型化的方向发展,在装卸过程中,岸桥小车运行速度和主起升速度都有了大幅提高。岸桥速度的提高增大了吊具及其所吊集装箱(以下称“吊重”)的摆动,增加了吊具定位的难度。负载或空载小车在岸桥前大梁上运动时,梁体会产生不小的弹性挠度。如何在不减小吊重的情况下减小梁体的挠度,是一个值得研究的问题。此外,当小车在梁上运动时,吊重会产生摆动,且摆动角度与小车运行参数和岸桥结构参数均有关系。因此,如果在岸桥设计伊始,就对结构参数及小车运行参数等进行优化,将得到更加可靠且性能更优的系统。
国内外学者在起重机械结构优化方面进行了相关研究。Abid
笔者同时考虑岸桥前大梁系统的结构和运行参数,对岸桥结构和小车运行等参数进行同步优化,既减小岸桥前大梁的挠度响应,又减小吊重的摆动,同时确保小车总运行时间不会太大。首先,建立岸桥前大梁系统简化力学模型;其次,结合假设模态法和拉格朗日方程建立运动方程并利用Runge‑Kutta法求解;最后,以岸桥关键参数为优化目标,采用多目标遗传算法同步优化结构参数与运行参数,根据具体需求从Pareto解集中提取选好的解,对岸桥结构设计与实际运营提供理论支持。
岸桥的结构组成如

图1 岸桥结构组成示意图
Fig.1 Structure composition of quayside container crane
岸桥前大梁系统结构复杂,很难直接对其进行准确的理论分析和求解,可通过建立有效的简化力学模型将复杂的问题简单化。以往的岸桥动力学研究在建立结构简化模型的过程中往往忽略吊重的摆动,将小车和吊重简化为移动载荷,研究移动载荷在梁体上运动时系统的动力学特
在Zrnić

图2 岸桥前大梁系统简化力学模型
Fig.2 Simplified mechanical model of front girder system of quayside container crane
简化模型中2个等效弹簧的刚度值k1和k2可以通过整机有限元模型结合柔度法得

图3 岸桥整机有限元模型
Fig.3 Finite element model of the QCC
E/Pa | I/ | A/ | M1/kg | M2/kg | L1/m | |
---|---|---|---|---|---|---|
0.089 25 | 0.138 8 | 1 634.4 | 4 606 | 13 380 | 29 | |
M2/kg | L/m |
k1/(N· |
k2/(N· | Mc/kg | ml/kg | l/m |
13 380 | 10 |
为了验证简化力学模型的有效性和模型参数设置的合理性,利用简化力学模型的参数建立有限元模型进行模态分析,与整机模型的模态分析结果进行对比,整机有限元模型(finite element method,简称FEM)和简化模型前3阶弯曲频率及模态置信因子(modal assurance criterion,简称MAC)如
频率 阶次 | 整机模型 频率/Hz | 简化模型 频率/Hz | 误差/% | MAC值/% |
---|---|---|---|---|
1 | 1.47 | 1.48 | 0.68 | 98.82 |
2 | 3.58 | 3.41 | -4.70 | 93.12 |
3 | 5.01 | 5.25 | 4.80 | 96.42 |
为了减少重复有限元计算带来的优化计算效率不高的问题,利用假设模态法和拉格朗日
(1) |
其中:为梁的第i阶假设模态函数,;为对应的广义坐标;N为广义坐标个数。
经验算N取6时,假设模态法可以很好地拟合岸桥前大梁系统的前3阶振动响应。
系统动能可以表示为
(2) |
其中:Tb,Tc,Tl,TM1和TM2分别为梁、小车、吊重、质量块M1和质量块M2的动能。
系统势能可以表示为
(3) |
其中:Ub,Uc,Ul,Uk1和Uk2分别为梁、小车、吊重、弹簧k1和弹簧k2的势能。
利用拉格朗日方程建立系统的运动微分方程,分别对广义坐标qi(i=1, 2, …, 6)和θ应用拉格朗日方程,可得系统的运动微分方程为
(4) |
(5) |
由
在岸桥前大梁系统简化力学模型中,内外侧拉杆的质量M1和M2、等效刚度k1和k2、梁单位长度质量ρ、前大梁内外侧拉杆到A点的距离L1和L2对系统响应有重要影响,所以选取这7个结构参数作为优化设计变量,但在实际设计中往往根据经验选取。
为了提高生产率,通常希望小车能够以最优的运行参数运行,即采用Bang‑Bang法使小车以最大加减速度和最大速度运行。小车的加速度为
(6) |
其中:amax为最大加(减)速度;ta为加(减)速时间;tf为总运行时间。
加(减)速时间ta和总运行时间tf为
(7) |
其中:vmax为最大运行速度;xd为期望的小车最大位移。
Bang‑Bang法的加速度及速度曲线如

图4 Bang-Bang加速度和速度曲线
Fig.4 Bang-Bang acceleration and velocity curve
选取岸桥前大梁系统的7个结构参数和2个小车运行参数共9个参数作为优化设计变量,各设计变量的含义、初始值及其取值范围如
设计变量 | 含义 | 初始值 | 取值范围 |
---|---|---|---|
x1 | L1/m | 29 | [27.55, 30.45] |
x2 | L2/m | 59 | [56.05, 61.95] |
x3 |
k1/(N· |
2.49×1 |
[1.99×1 |
x4 |
k2/(N· |
4.58×1 |
[3.66×1 |
x5 | M1/kg | 4 606 | [3 685, 5 527] |
x6 | M2/kg | 13 380 | [10 704, 16 056] |
x7 | 1 634.4 | [1 308, 1 961] | |
x8 |
vmax/(m· | 4.0 | [2.5, 4.0] |
x9 |
amax/(m· | 0.6 | [0.5, 0.8] |
根据
其中:,为选定的9个设计变量;,分别为各设计变量的上、下界值;f1(x)为小车总运行时间,即
采用多目标遗传算法对上述优化模型进行求解,优化过程中,多目标遗传算法的种群大小为100,进化代数为200,交叉概率为0.8,变异概率为0.2。


图5 Pareto解集及其投影图
Fig.5 Pareto set and projection
(8) |
其中:w1,w2和w3分别为各目标函数的权重,w1+w2+w3=1;,和分别为各目标函数的初始值。
为了使3个目标函数的值均较小,初选w1=w2=w3=1/3,选取使得F最小的解,如
符号 | x1 | x2 | x3 | x4 | x5 | x6 | x7 | x8 | x9 | f1(x) | f2(x) | f3(x) |
---|---|---|---|---|---|---|---|---|---|---|---|---|
含义 | L1/m | L2/m |
k1/(N· |
k2/(N· | M1/kg | M2/kg |
ρ/(kg· |
vmax/(m· |
amax/(m· | tf/s | θmax/(°) | ymax/mm |
解1 | 30.28 | 61.66 |
2.49×1 |
5.24×1 | 5 237 | 14 116 | 1 687 | 4.0 | 0.8 | 21.3 | 14.73 | 69.3 |
解2 | 30.45 | 59.00 |
2.96×1 |
4.31×1 | 4 882 | 13 380 | 1 667 | 3.5 | 0.5 | 25.7 | 5.79 | 71.4 |
解3 | 30.45 | 59.00 |
2.94×1 |
5.45×1 | 4 836 | 13 380 | 1 700 | 2.7 | 0.6 | 28.1 | 12.32 | 63.1 |
解4 | 30.16 | 59.00 |
2.86×1 |
5.40×1 | 4 468 | 12 443 | 1 683 | 3.7 | 0.5 | 24.6 | 6.13 | 64.5 |

图6 优化前后各解前大梁中点挠度响应结果对比
Fig.6 Comparison of midpoint deflection of girder before and after optimization

图7 优化前后各解吊重摆角响应结果对比
Fig.7 Comparison of swinging angle before and after optimization
从
1) 建立了考虑吊重摆角的岸桥前大梁系统简化模型,并与整机有限元模型进行比较,证明了简化模型的有效性和参数设置的合理性。
2) 利用假设模态法和拉格朗日方程得到了系统的运动微分方程,并利用数值方法进行求解,提高了优化效率。
3) 对岸桥前大梁系统的结构参数和运行参数进行了同步优化,利用权重系数法从Pareto解集中选取得到最优参数组合,优化后系统的挠度响应、摆角响应和小车运行时间均有了一定程度的减小,为岸桥防摇提供了一种新的思路。
参考文献
ABID M, AKMAL M H, PARVEZ S. Optimization of box type girder of overhead crane[C]∥13th International Conference on Advanced Design and Manufacture. Surrey, England: Springer-Verlag London Ltd.,2008:609-618. [百度学术]
SAVKOVIC M M, BULATOVIC R R, GASIC M M, et al. Optimization of the box section of the single-grider bridge crane by GRG algorithm according to domestic standards and eurocodes[C]∥International Scientific Conference on Defensive Technologies OTEH. Belgrade Serbia:[s.n.],2016:211-217. [百度学术]
SAVKOVIĆ M M, BULATOVIĆ R R, GAŠIĆ M M, et al. Optimization of the box section of the main girder of the single-girder bridge crane by applying biologically inspired algorithms[J]. Engineering Structures, 2017, 148:452-465. [百度学术]
张慧博. 超大型岸边集装箱起重机金属结构静动态特性分析与轻量化研究[D]. 上海:上海交通大学, 2008. [百度学术]
计三有, 苏运波. 基于神经网络的起重机臂架结构优化[J]. 机械设计与制造, 2008(8):23-25. [百度学术]
JI Sanyou, SU Yunbo. Gantry crane′s arm frame structure parameters optimization based on neural network[J]. Machinery Design & Manufacture, 2008(8):23-25.(in Chinese) [百度学术]
张文正, 蒋伟康. 大型岸桥结构的轻量化技术研究及应用[J]. 机械强度, 2012(4):615-620. [百度学术]
ZHANG Wenzheng, JIANG Weikang. Optimal design for lightweight of structure of quayside container crane[J]. Journal of Mechanical Strength, 2012(4):615-620.(in Chinese) [百度学术]
夏浩军, 邱荣华. 改进遗传算法的起重机主梁优化设计方法[J]. 控制工程, 2017, 24(7):1415-1418. [百度学术]
XIA Haojun, QIU Ronghua. Optimized design method of crane girder based on improved genetic algorithm[J]. Control Engineering of China, 2017, 24(7):1415-1418.(in Chinese) [百度学术]
谢伟平,黄金,周家玲 ,等.重物-桥吊耦合系统振动分析[J].振动与冲击,2015,34(15):127-132. [百度学术]
XIE Weiping, HUANG Jin, ZHOU Jialing, et al. Vibration analysis of a suspension weight-bridge crane coupled system[J]. Journal of Vibration and Shock, 2015, 34(15):127-132.(in Chinese) [百度学术]
ZRNIĆ N D, OGUAMANAM D, BOŠNJAK S. Dynamics and modelling of mega quayside container cranes[J]. Fme Transactions, 2006, 34(4):193-198. [百度学术]
ZRNIĆ N D, HOFFMANN K, BOSNJAK S M. Modelling of dynamic interaction between structure and trolley for mega container cranes[J]. Mathematical and Computer Modelling of Dynamical Systems, 2009, 15(3):295-311. [百度学术]
MICHALTSOS G T . Dynamic behavior of a single-span subjected to loads moving with variable speeds [J]. Journal of Sound & Vibration, 2002, 258(2):359-372. [百度学术]
ABDEL-RAHMAN E M, NAYFEH A H, MASOUD Z N. Dynamics and control of cranes: a review[J]. Journal of Vibration & Control, 2003, 9(7):863-908. [百度学术]
邹军, 陈志坚. 桥式起重机水平运行及抓斗防摆规律研究[J]. 山东大学学报(自然科学版), 1998, 33(4):393-397. [百度学术]
ZOU Jun, CHEN Zhijian. The research about horizontal moving rule of bridge crane and swing rule of grab[J]. Journal of Shandong University (Natural Science), 1998, 33(4):393- 397.(in Chinese) [百度学术]
王阳. 船用起重机防摇装置设计研究[J]. 舰船科学技术, 2013, 35(7):105-108. [百度学术]
WANG Yang. Research and design of an anti-swing mechanism for a shipborne crane[J]. Ship Science and Technology, 2013, 35(7):105-108.(in Chinese) [百度学术]
宗周红, 任伟新. 桥梁有限元模型修正和模型确认[M]. 北京:人民交通出版社, 2012:16. [百度学术]
张策. 机械动力学 [M]. 北京:高等教育出版社, 2008:193-196. [百度学术]
GB/T 3811―2008 起重机设计规范[S]. 北京: 中国国家标准化管理委员会, 2008. [百度学术]
吴庆祥,汪小凯,王孝文,等.起重机前馈防摇控制算法仿真及实验研究[J].武汉理工大学学报,2016,38(6):109-116. [百度学术]
WU Qingxiang, WANG Xiaokai, WANG Xiaowen, et al. Simulation and experimental study of feedforward and anti -swing control algorithm for overhead cranes[J]. Journal of Wuhan University of Technology, 2016, 38(6):109-116.(in Chinese) [百度学术]