使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

岸桥前大梁结构参数和运行参数的同步优化  PDF

  • 秦仙蓉
  • 詹澎明
  • 王玉龙
  • 郝婼兰
  • 孙远韬
  • 张氢
同济大学机械与能源工程学院 上海,200092

中图分类号: TH21U653.921

最近更新:2023-04-28

DOI:10.16450/j.cnki.issn.1004-6801.2023.02.024

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

对岸桥前大梁的结构参数和运行参数进行了同步优化。首先,建立前大梁系统简化力学模型,对比整机有限元模型与简化模型的有限元模型得到前3阶模态频率,验证了简化模型的有效性;其次,为了提高优化过程的计算效率,利用假设模态法及拉格朗日法得到简化模型的动力学方程;然后,以Runge‑Kutta法计算得到简化模型的吊重摆角最大值、岸桥前大梁中点处的最大挠度及小车总运行时间为优化目标,采用多目标遗传算法对岸桥前大梁系统的结构参数及小车运行参数进行同步优化,得到Pareto解集;最后,根据系统性能的具体需求从解集中选择合理的解,指导岸桥前大梁系统的结构设计以及小车运行参数的设定,为岸桥防摇提供了一个新的思路。

引 言

随着集装箱运输船舶向大型化的方向发展,在装卸过程中,岸桥小车运行速度和主起升速度都有了大幅提高。岸桥速度的提高增大了吊具及其所吊集装箱(以下称“吊重”)的摆动,增加了吊具定位的难度。负载或空载小车在岸桥前大梁上运动时,梁体会产生不小的弹性挠度。如何在不减小吊重的情况下减小梁体的挠度,是一个值得研究的问题。此外,当小车在梁上运动时,吊重会产生摆动,且摆动角度与小车运行参数和岸桥结构参数均有关系。因此,如果在岸桥设计伊始,就对结构参数及小车运行参数等进行优化,将得到更加可靠且性能更优的系统。

国内外学者在起重机械结构优化方面进行了相关研究。Abid

1对桥式起重机主梁内部的横筋和纵筋的数量与位置进行规划,对主梁的几何参数进行优化设计,实现了起重机的轻量化。Savkovic2运用拉格朗日乘子法确定了桥式起重机梁截面几何参数,在满足主梁的刚度和稳定性约束条件的情况下,得到了更轻的结构。Savković3以起重机主梁最轻为优化目标,以主梁的强度、刚度、动刚度和稳定性为约束条件,采用仿生算法进行了优化计算,取得了良好的优化效果。张慧4以岸桥金属结构为研究对象,运用有限元分析方法,以结构的强度、刚度及模态频率为约束,以重量最轻为目标函数,实现了岸桥结构的轻量化设计。计三有5结合人工神经网络,对起重机臂架结构进行了优化,改善了结构的应力状态。张文正6以大型岸桥的强度、刚度和第1阶整机模态频率为约束条件,减轻了岸桥结构的重量。夏浩军7运用改进的遗传算法,以截面积最小为优化目标对起重机主梁截面几何尺寸进行了优化设计,优化结果满足结构强度、刚度和稳定性的要求。但是,上述研究的关注点主要放在桥式或门式起重机主梁的结构优化设计上,通常以结构的质量和应力状态为优化目标,而以岸桥的动力学响应和吊重的摆角为优化目标的研究较少。对于实际运行的岸桥而言,小车运行过程中结构的动响应和吊重的摆动不可忽略,且吊重的摆动与结构的振动是相互耦合8,因此有必要对岸桥前大梁的挠度和吊重的摆角进行优化,这既能提高吊重的定位速度和精度,提高装卸及运输的效率,也能减小因吊重摆动幅度过大或结构振动而导致的安全事故发生的几率。

笔者同时考虑岸桥前大梁系统的结构和运行参数,对岸桥结构和小车运行等参数进行同步优化,既减小岸桥前大梁的挠度响应,又减小吊重的摆动,同时确保小车总运行时间不会太大。首先,建立岸桥前大梁系统简化力学模型;其次,结合假设模态法和拉格朗日方程建立运动方程并利用Runge‑Kutta法求解;最后,以岸桥关键参数为优化目标,采用多目标遗传算法同步优化结构参数与运行参数,根据具体需求从Pareto解集中提取选好的解,对岸桥结构设计与实际运营提供理论支持。

1 岸桥前大梁‑小车‑吊重系统简化建模

岸桥的结构组成如图1所示,主要由门框系统、门框连接系统、拉杆系统、梯形架及主梁等部分组成。由图可以看出,岸桥前大梁处于靠海侧,尺寸大,质量重,其最外端处于自由状态,中间仅由2组拉杆约束,且在岸桥工作时空载或负重小车在前大梁上频繁作业,其动力学特性是研究岸桥动力学特性时的重点和难点。

图1  岸桥结构组成示意图

Fig.1  Structure composition of quayside container crane

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

9‑10。然而,当空载或负载小车在岸桥主梁上运动时,吊重会产生摆动,且小车运行参数和岸桥结构参数均会影响吊重摆动的角度,吊重摆角过大会影响吊具的定位精度,影响装卸效率。因此,有必要建立岸桥前大梁‑小车‑吊重耦合系统的简化力学模型,同时考虑前大梁、小车和吊重之间的相互耦合作用。此外,为了提高求解优化问题的计算效率,笔者利用数值方法求解耦合振动响应。

在Zrnić

9‑10建立的岸桥前大梁系统简化力学模型的基础上,考虑吊重的摆动,建立岸桥前大梁‑小车‑吊重耦合系统的简化力学模型,如图2所示。假设梁为Euler‑Bernoulli梁,梁的结构阻尼对振动影响较11‑12,将钢丝绳简化为刚性无质量13,吊重简化为一个通过刚性钢丝绳悬吊在小车上摆动的集中质量。

图2  岸桥前大梁系统简化力学模型

Fig.2  Simplified mechanical model of front girder system of quayside container crane

图2中:M1M2分别为内外侧拉杆的质量;k1k2分别为除前大梁之外的结构等效刚度;ρ为梁的单位长度质量;E为弹性模量;I为截面惯性矩;L1L2分别为前大梁内侧和外侧拉杆到前大梁与海侧门框的连接位置(图中A点)的距离;L为前大梁总长;Mc为小车质量;ml为吊重质量;l为钢丝绳的长度。前大梁与岸桥其他结构在A点以铰接形式连接,这是因为岸桥前大梁需要根据工作需求进行俯仰操作。

简化模型中2个等效弹簧的刚度值k1k2可以通过整机有限元模型结合柔度法得

14。利用ANSYS建立的整机有限元模型,其中海陆侧门框、门框联系横梁、前后大梁、后撑管和梯形架采用三维梁单元BEAM189模拟,前后大梁的拉杆系统采用杆单元LINK10模拟,其余附属结构采用质量点单元MASS21模拟,建立的岸桥整机有限元模型如图3所示。模型的其余参数可以通过设计图纸得到,最终的简化力学模型的参数初始取值如表1所示。

图3  岸桥整机有限元模型

Fig.3  Finite element model of the QCC

表1  岸桥前大梁系统简化力学模型参数初始取值
Tab.1  Parameter values initial of simplified mechanical model for front girder system of bridge
E/PaI/m4A/m2ρ/(kg·m-1)M1/kgM2/kgL1/m
2.1×1011 0.089 25 0.138 8 1 634.4 4 606 13 380 29
M2/kg L/m k1/(N·m-1) k2/(N·m-1) Mc/kg ml/kg l/m
13 380 4.58×106 1×105 4.58×106 1×105 6.5×104 10

为了验证简化力学模型的有效性和模型参数设置的合理性,利用简化力学模型的参数建立有限元模型进行模态分析,与整机模型的模态分析结果进行对比,整机有限元模型(finite element method,简称FEM)和简化模型前3阶弯曲频率及模态置信因子(modal assurance criterion,简称MAC)如表2所示。由表2可知,2种模型的前3阶频率的最大误差为4.8%,其前3阶模态MAC

15均趋近于1,最小值为93.12%,说明模态相关性较高。因此,图2所示的岸桥前大梁系统中的前大梁简化力学模型及其参数设置是合理的,可以用于求解前大梁中点处的挠度和吊重摆角。

表2  整机FEM模型和简化模型前3阶弯曲频率及MAC值
Tab.2  The first three frequencies and MAC of the whole FEM and the simplified model

频率

阶次

整机模型

频率/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

为了减少重复有限元计算带来的优化计算效率不高的问题,利用假设模态法和拉格朗日

10建立简化力学模型的运动微分方程。对于图2所示的弹性梁,梁上xt时刻的挠度可以表示为

yx,t=i=1Nψixqit (1)

其中:ψix为梁的第i阶假设模态函数,ψix=siniπx/2Lqit为对应的广义坐标;N为广义坐标个数。

经验算N取6时,假设模态法可以很好地拟合岸桥前大梁系统的前3阶振动响应。

系统动能可以表示为

T=Tb+Tc+Tl+TM1+TM2 (2)

其中:TbTcTlTM1TM2分别为梁、小车、吊重、质量块M1和质量块M2的动能。

系统势能可以表示为

U=Ub+Uc+Ul+Uk1+Uk2 (3)

其中:UbUcUlUk1Uk2分别为梁、小车、吊重、弹簧k1和弹簧k2的势能。

利用拉格朗日方程建立系统的运动微分方程,分别对广义坐标qii=1, 2, …, 6)和θ应用拉格朗日方程,可得系统的运动微分方程为

i=160Lρψixψjxdx+M1ψiL1ψjL1+M2ψiL2ψjL2+Mc+mlψixcψjxcq¨i(t)+2Mc+mlx˙ci=16ψi'xcψjxcq˙it+i=16EI0Lψixψjxdx+k1ψiL1ψjL1+k2ψiL2ψjL2+Mc+mlx¨cψi'xcψjxc+x˙c2ψixcψjxcqi(t)-mllsinθψjxcθ¨-mllcosθψjxcθ˙2-
Mc+mlgψjxc=0                  (j=1,2,,6) (4)
-mllsinθi=16ψixcq¨i(t)-2mllsinθx˙ci=16ψi'xcq˙it-mllsinθx¨ci=16ψi'xc+x˙c2i=16ψixcqit+
mll2θ¨-mllx˙csinθ θ˙+mllx¨ccosθ+mlglsinθ=0                                                                                                      (5)

式(4)式(5)可知,该系统为7元2阶非线性时变微分方程组,难以求得解析解,因此笔者利用Runge‑Kutta积分方

16对微分方程进行求解,得到梁中点处的挠度响应与吊重摆角响应。

2 岸桥前大梁系统优化问题的数学建模

在岸桥前大梁系统简化力学模型中,内外侧拉杆的质量M1M2、等效刚度k1k2、梁单位长度质量ρ、前大梁内外侧拉杆到A点的距离L1L2对系统响应有重要影响,所以选取这7个结构参数作为优化设计变量,但在实际设计中往往根据经验选取。

为了提高生产率,通常希望小车能够以最优的运行参数运行,即采用Bang‑Bang法使小车以最大加减速度和最大速度运行。小车的加速度为

x¨t=amax          (0tta)         0               tattf-ta-amax      tf-tattf (6)

其中:amax为最大加(减)速度;ta为加(减)速时间;tf为总运行时间。

加(减)速时间ta和总运行时间tf

ta=vmax/amaxtf=vmax/amax+xd/vmax (7)

其中:vmax为最大运行速度;xd为期望的小车最大位移。

Bang‑Bang法的加速度及速度曲线如图3所示。假设小车以图4所示的运动规律运行,则需确定xdamaxvmax这3个独立的运行参数,由于小车最大运行位移xd不可改变(本研究中取65 m),只需确定amaxvmax这2个参数即可确定小车的总运行时间,因此选取最大加(减)速度和最大运行速度作为优化设计变量。

图4  Bang-Bang加速度和速度曲线

Fig.4  Bang-Bang acceleration and velocity curve

选取岸桥前大梁系统的7个结构参数和2个小车运行参数共9个参数作为优化设计变量,各设计变量的含义、初始值及其取值范围如表3所示。其中,由于L1L2的取值对梁的动力学特征影响较大,并且实际设计过程中的可变范围较小,所以取其变化范围为初值±5%,其余5个结构参数的变化范围取为初值±20%。此外,本研究中的岸桥,其小车运行额定速度为4 m/s,属于中高速的岸桥,根据《起重机设计规范

17选取其加(减)速度值的范围。

表3  优化设计变量含义、初始值及取值范围
Tab.3  Meaning, initial value and range of optimal design variables
设计变量含义初始值取值范围
x1 L1/m 29 [27.55, 30.45]
x2 L2/m 59 [56.05, 61.95]
x3 k1/(N·m-1) 2.49×107 [1.99×107, 2.99×107]
x4 k2/(N·m-1) 4.58×106 [3.66×106, 5.50×106]
x5 M1/kg 4 606 [3 685, 5 527]
x6 M2/kg 13 380 [10 704, 16 056]
x7 ρ/kg·m-1 1 634.4 [1 308, 1 961]
x8 vmax/(m·s-1) 4.0 [2.5, 4.0]
x9 amax/(m·s-2) 0.6 [0.5, 0.8]

根据式(4)式(5)可知,前大梁、小车以及吊重三者之间存在动力学耦合关系。为了获得岸桥前大梁系统的综合最优特性,需要同时优化岸桥前大梁结构参数和小车运行参数,既减小岸桥前大梁的挠度响应,又减小吊重的摆动,同时确保小车总运行时间不会太大。因此,笔者提出以小车总运行时间f1x)、吊重最大摆角响应f2x)及梁中点处的最大挠度响应f3x)为优化设计目标函数。本优化问题为多目标非线性约束优化问题,优化数学模型为

findxminf1x,f2x,f3xs.t.xLxxU

其中:x=x1,x2,,x9T,为选定的9个设计变量;xLxU分别为各设计变量的上、下界值;f1x)为小车总运行时间,即式(7)中的tff2x),f3x)分别为吊重摆角响应最大值和梁中点处挠度响应的最大值,其无明确解析表达式,需通过Runge‑Kutta法对式(4)式(5)进行求解得到。

3 基于遗传算法的多目标优化

采用多目标遗传算法对上述优化模型进行求解,优化过程中,多目标遗传算法的种群大小为100,进化代数为200,交叉概率为0.8,变异概率为0.2。

图5为优化后归一化目标函数值的Pareto解集及其投影,图中共有60组解,每组解中3个目标函数所占的权重各不相同,从中提取f1x),f2x)及f3x)分别最小的解,即图5中的解1、解2和解3。此外,采用权重系数法从Pareto解集中选取一个综合性能较好的解,定义综合指标F

  

  

图5  Pareto解集及其投影图

Fig.5  Pareto set and projection

F=w1f1xf1*x+w2f2xf2*x+w3f3xf3*x (8)

其中:w1w2w3分别为各目标函数的权重,w1+w2+w3=1;f1*xf2*xf3*x分别为各目标函数的初始值。

为了使3个目标函数的值均较小,初选w1=w2=w3=1/3,选取使得F最小的解,如图5中的解4。优化后设计变量值和目标函数值如表4所示。为了更清楚地表现4组解的差别,分别求解4组参数下的梁中点处的挠度响应曲线和吊重摆角响应曲线,优化前后结果如图67所示。图7中4组参数计算得到的吊重摆角响应趋势、响应幅值与文献[

18]中类似工况下的实验结果基本一致,说明了图7中的各组解吊重响应的合理性,也验证了系统运动微分方程的正确性。

表4  优化后设计变量值和目标函数值
Tab.4  Optimized design variable values and objective function values
符号x1x2x3x4x5x6x7x8x9f1(x)f2(x)f3(x)
含义 L1/m L2/m k1/(N·m-1) k2/(N·m-1) M1/kg M2/kg ρ/(kg·m-1) vmax/(m·s-1) amax/(m·s-2) tf/s θmax/(°) ymax/mm
解1 30.28 61.66 2.49×107 5.24×106 5 237 14 116 1 687 4.0 0.8 21.3 14.73 69.3
解2 30.45 59.00 2.96×107 4.31×106 4 882 13 380 1 667 3.5 0.5 25.7 5.79 71.4
解3 30.45 59.00 2.94×107 5.45×106 4 836 13 380 1 700 2.7 0.6 28.1 12.32 63.1
解4 30.16 59.00 2.86×107 5.40×106 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

表4和图67中可以看出:解1的小车总运行时间最小而吊重的摆角响应最大;解2的吊重摆角响应值最小而梁中点处的挠度响应最大;解3的梁中点挠度响应值最小而小车的运行时间最长。根据权重系数法得到的最优解4的3个目标函数值均较小,即在解4对应的结构参数和小车运行参数设置下,岸桥前大梁系统的小车运行时间、梁中点处的挠度响应及摆角响应均较小。在实际设计过程中,也可以根据系统性能的具体需求,改变式(8)中各目标函数的权重值,从Pareto解集中选择合理解,以指导岸桥前大梁系统的结构设计以及小车运行参数的设定。

4 结 论

1) 建立了考虑吊重摆角的岸桥前大梁系统简化模型,并与整机有限元模型进行比较,证明了简化模型的有效性和参数设置的合理性。

2) 利用假设模态法和拉格朗日方程得到了系统的运动微分方程,并利用数值方法进行求解,提高了优化效率。

3) 对岸桥前大梁系统的结构参数和运行参数进行了同步优化,利用权重系数法从Pareto解集中选取得到最优参数组合,优化后系统的挠度响应、摆角响应和小车运行时间均有了一定程度的减小,为岸桥防摇提供了一种新的思路。

参考文献

1

ABID MAKMAL M HPARVEZ S. Optimization of box type girder of overhead crane[C]∥13th International Conference on Advanced Design and Manufacture. Surrey, EnglandSpringer-Verlag London Ltd.2008609-618. [百度学术] 

2

SAVKOVIC M MBULATOVIC R RGASIC M Met 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.]2016211-217. [百度学术] 

3

SAVKOVIĆ M MBULATOVIĆ R RGAŠIĆ M Met al. Optimization of the box section of the main girder of the single-girder bridge crane by applying biologically inspired algorithms[J]. Engineering Structures2017148452-465. [百度学术] 

4

张慧博. 超大型岸边集装箱起重机金属结构静动态特性分析与轻量化研究[D]. 上海上海交通大学2008. [百度学术] 

5

计三有苏运波. 基于神经网络的起重机臂架结构优化[J]. 机械设计与制造20088):23-25. [百度学术] 

JI SanyouSU Yunbo. Gantry crane′s arm frame structure parameters optimization based on neural network[J]. Machinery Design & Manufacture20088):23-25.(in Chinese) [百度学术] 

6

张文正蒋伟康. 大型岸桥结构的轻量化技术研究及应用[J]. 机械强度20124):615-620. [百度学术] 

ZHANG WenzhengJIANG Weikang. Optimal design for lightweight of structure of quayside container crane[J]. Journal of Mechanical Strength20124):615-620.(in Chinese) [百度学术] 

7

夏浩军邱荣华. 改进遗传算法的起重机主梁优化设计方法[J]. 控制工程2017247):1415-1418. [百度学术] 

XIA HaojunQIU Ronghua. Optimized design method of crane girder based on improved genetic algorithm[J]. Control Engineering of China2017247):1415-1418.(in Chinese) [百度学术] 

8

谢伟平黄金周家玲.重物-桥吊耦合系统振动分析[J].振动与冲击20153415):127-132. [百度学术] 

XIE WeipingHUANG JinZHOU Jialinget al. Vibration analysis of a suspension weight-bridge crane coupled system[J]. Journal of Vibration and Shock20153415):127-132(in Chinese) [百度学术] 

9

ZRNIĆ N DOGUAMANAM DBOŠNJAK S. Dynamics and modelling of mega quayside container cranes[J]. Fme Transactions2006344):193-198. [百度学术] 

10

ZRNIĆ N DHOFFMANN KBOSNJAK S M. Modelling of dynamic interaction between structure and trolley for mega container cranes[J]. Mathematical and Computer Modelling of Dynamical Systems2009153):295-311. [百度学术] 

11

MICHALTSOS G T . Dynamic behavior of a single-span subjected to loads moving with variable speeds [J]. Journal of Sound & Vibration20022582):359-372. [百度学术] 

12

ABDEL-RAHMAN E MNAYFEH A HMASOUD Z N. Dynamics and control of cranes: a review[J]. Journal of Vibration & Control200397):863-908. [百度学术] 

13

邹军陈志坚. 桥式起重机水平运行及抓斗防摆规律研究[J]. 山东大学学报(自然科学版)1998334):393-397. [百度学术] 

ZOU JunCHEN Zhijian. The research about horizontal moving rule of bridge crane and swing rule of grab[J]. Journal of Shandong University (Natural Science)1998334):393- 397.(in Chinese) [百度学术] 

14

王阳. 船用起重机防摇装置设计研究[J]. 舰船科学技术2013357):105-108. [百度学术] 

WANG Yang. Research and design of an anti-swing mechanism for a shipborne crane[J]. Ship Science and Technology2013357):105-108.(in Chinese) [百度学术] 

15

宗周红任伟新. 桥梁有限元模型修正和模型确认[M]. 北京人民交通出版社201216. [百度学术] 

16

张策. 机械动力学 [M]. 北京高等教育出版社2008193-196. [百度学术] 

17

GB/T 3811―2008 起重机设计规范[S]. 北京中国国家标准化管理委员会2008. [百度学术] 

18

吴庆祥汪小凯王孝文.起重机前馈防摇控制算法仿真及实验研究[J].武汉理工大学学报2016386):109-116. [百度学术] 

WU QingxiangWANG XiaokaiWANG Xiaowenet alSimulation and experimental study of feedforward and anti -swing control algorithm for overhead cranes[J]. Journal of Wuhan University of Technology2016386):109-116.(in Chinese) [百度学术]