摘要
下水驳船是一种高效的船舶产品及大型海洋结构物的下水和运载工具,在造船业兴旺的今天,可大大减轻船厂对干船坞和船台的使用需求。随着造船模式的发展和建造船舶吨位的增加,下水产品重量也在不断增加。在使用下水驳船作为产品下水工具时,如何保证整个作业过程平稳、安全和快速地完成,是一个非常值得研究的问题。本文针对下水驳船使用过程中各调载水舱调载方案的实时确定,上驳产品对驳船施加载荷的计算方法以及待调载水舱阀门的打开时刻确定等几个方面作了深入的研究。
下水驳船在整个作业过程中,需要始终保持特定的、适宜的浮态,而驳船浮态的调整是通过调节各压载水舱的水量来实现的。何时进行调载,需要调节哪些压载水舱以及调节多少水量,这都是影响整个下水过程的关键因素。通常采用传统的人工手动方式进行调节,但其主要依靠操作人员经验,而且操作繁琐,效率低下。本文首先依照力学平衡原理,分析了驳船的受力特点,建立了力学平衡方程组模型。针对该模型,以总调水量最小为原则,采用牛顿迭代方法求解,并引入了差商代替导数的思想求解雅可比迭代矩阵;针对水舱按行列布置的较有规律的驳船,研究了四种特殊的有针对性的模型求解方法,并通过算例分析比较了这四种方法的特点,给出了方法选择的依据和准则。接着,结合上驳时间要求尽量短的特点,以力学平衡方程和调水舱水量的限制为约束,调水量绝对值之和为目标函数,建立了驳船实时配载优化模型。针对优化模型的非线性规划特点,对模型进行了归一化处理,分别采用了惩罚函数法和乘子法进行求解。通过工程实例计算分析,证明了所建模型和对应解法的正确性和适用性,同时也比较得出了两种模型及求解方法的优缺点。对于舱室布置对称性明显的驳船,适于采用力学平衡方程组模型求解,其计算速度快,结果稳定。而对于舱室布置任意性较强的驳船,则宜于采用优化模型求解,尽管计算速度稍慢,但其适应性较强,配载结果也较为合理。
在下水驳船移运产品的上驳过程中,需要实时确定上驳产品对驳船所施加的压力载荷,从而为配载模型的计算提供必要和准确合理的已知数据。而该部分载荷的计算方法与产品上驳的方式密切相关。本文主要分析了滑道和气囊两种移运方式下的该载荷计算方法。对于滑道移运方式,分别采用了对已上驳部分产品重量载荷直接积分方法、弹性基础梁法和弹性支座的有限元法进行了计算分析。实例计算结果表明,对于支承滑道覆盖整个产品上驳长度的情形,可采用产品重量载荷直接积分方法进行计算,该计算方法具有简单和快速的优点,同时也不会与实际重量产生较大的误差;而对于产品相对于其支承滑道存在较大悬臂端的情形,则宜采用弹性基础梁法进行计算,其公式虽然复杂,但计算结果精度较高,且与实际情形较为接近;有限元法则适用于各种情形,并且计算结果具有很高的精度,但这往往又涉及到其他有限元软件的参与,不能满足上驳过程的实时计算要求。对于气囊移运方式,驳船承受产品的压力载荷既与驳船上支承产品的各气囊实时压力有关,又与驳船上支承产品的气囊根数及其位置密切相关。首先,本文依照力平衡原理,通过适当的假设,给出了支承产品气囊压力的实时计算方法。接着,依照运动学原理,分析了气囊和产品的相对运动关系,给出了不同的产品上驳长度时驳船上支承产品的气囊根数和其对应位置坐标的计算公式。最后,通过相应气囊的合力计算和坐标转换,得到了驳船所承受的气囊压力载荷计算公式。通过实例计算分析表明,在气囊移运方式下,驳船所承受产品的载荷具有离散性、突变性和周期性的特点,一般不宜采用滑道移运的方式进行近似计算。
当配载方案确定后,如何控制阀门开关时刻来保证整个调水过程中驳船处于最佳的浮态,对于产品上驳作业过程具有重要的意义。通常,阀门开关操作存在两种不同的作业方式。在两种作业方式下,阀门开关时刻都会影响调水过程中驳船的浮态。而阀门开关时刻的确定主要与产品的上驳速度、调水舱室的阀门流量和潮位变化密切相关。通过分析调水过程特点,本文选取适当的驳船浮态表征函数作为目标函数,阀门调水时间范围作为约束条件,建立了作业过程中调水舱室阀门开关时刻的优化模型,分别研究了惩罚函数法、网格法和模拟退火法在模型求解中的应用。通过实例计算分析表明,所建的优化模型是正确和有效的,采用的解法均可使驳船在调水过程中保持较好的浮态。但采用网格法的计算量随调水舱室数量或调载时间等分份数的增加而急剧增长,导致求解时间不可接受;惩罚函数法虽求解时间较短,但调水过程中驳船浮态相对差些;而采用模拟退火法所得到的驳船浮态最好,只是计算时间稍长。因此,对于计算时间要求较短而浮态限制较松的情形可采用惩罚函数法;反之,则宜于采用模拟退火法或网格法。
The launch barge is a highly-efficient tool to launch and carry ship products and it can reduce the demand to the dry dock and building berth of a shipyard greatly for today's prosperous shipbuilding industry. The weight of the launched products increases a lot with the development of shipbuilding mode and increasing of the ship tonnage. To ensure the whole launching process completed stably, safely and quickly is worthy of investigation for the launch barge. In this paper, the determination of real-time stowage planning, the calculating method about the loads of ship products and the choice about valve open time of ballast tanks are studied deeply in the launching process respectively.
Special and proper floating conditions of the barge need to be remained in the whole launching process all the while. And the floatation is adjusted by changing the water volume of each ballast tank. When, where and how much water needs to be adjusted are the key factors influencing the whole launching process. Usually this process is completed manually which is of low efficiency and complex operation. In this paper, firstly by analyzing mechenics characteristics of the barge, a static equilibrium equations model is established. And it is solved by Newton iteration method on the least adjusted water criteria in which the idea that derivative is replaced with difference quotient is involved in calculating each item of Jacobi matrix. For those barges which subvision has obvious orderly ranks, four special solving methods are studied to the model, and through an example's analysis, these four methods' characteristics are shown and concluded. Then the real-time stowage optimization model is established combined with the characteristic of load-out time as short as possible. It is constraint of static equilibrium equations and the upper and lower limits of the water quantity, and its objective function is determined as the absolute value of the adjusting water. Penalty function method and multiplier method are adopted as to the nonlinear characteristic of the optimization model and normalization processing is applied meantime. The results of the practical examples prove the correctness and applicability of the two models and corresponding solving method. And at the same time the advantages and disadvantages between them are pointed out. That is, the former model is suitable to the case that the barge subvision is symmetrically, and it is of fast calculation speed and stable results, while the latter model is better to the barge subvision with strong arbitrariness. Although the calculation speed is slower slightly, the optimization model has good adaptability and more reasonable load-out results.
The pressure exerted onto the barge by the transferring object needs to be determined real-timely in the load-out process in order to provide necessary data for calculation of the stowage model, and its calculation is closely related to the transferring mode. This paper mainly discussed the calculation method under two modes transferring by slideway and gasbags. As to the slideway mode, directional integral method to weight load of the project on the barge, elastic foundation beam method and finite element method with elastic support are analyzed respectively. The examples show that directional integral method is adapted to the case that the support slideway covers the whole object transferring length because of its simplicity, speedness and result close to the practice, and elastic foundation beam method having advantage of high precision should be adopted if there is a large cantilever end between the object and supporting slideway although its formula are very complex. The finite element method is suitable for all cases and of high precision, but it doesn't satisfy the demand of real-time calculation because of participation of other software. To the gasbag mode, the pressure is related to both real-time pressure of each gasbag and number and position of those gasbags on the barge supporting the object closely. Firstly, the real-time calculation of each gasbag pressure is derived according to static equilibrium principle with some proper assumptions. Secondly, the relative movement relation between the gasbag and object is analyzed in terms of kinematics principle, and the calculation formulas of gasbags number and corresponding coordinate are given for different transferring length. Finally, the calculation formula of pressure onto the barge from the gasbags can be obtained through the resultant force calculation and coordinate transformation of corresponding gasbags. The example results demonstrate that the loads onto the barge from the object have three characteristics of discreteness, mutagenicity and periodicity under the gasbag transferring mode. So the slideway mode method shouldn't be used to calculate approximately in this case.
How to control the valve open time is of significance to ensure the barge being desirable floatation in the load-out process after the stowage planning determined. Usually, there are two different working modes in the operation of the valves, and the barge floatation in the process will be influenced by the valve open time for both modes. The determination of valve open time is closely related to the object movement velocity, the flux of the ballast tank and the variation in tidal level. Analyzing the characteristics of the water adjusting process, a proper function for characterizing the barge floatation is chosen as objective function, the range of adjusting time is regarded as constraint condition, thus an optimization model on valve open time for ballast tanks of the launch barge can be built, and penalty function method, grid method and simulated annealing method are adopted to solve this model. The example results demonstrate that the optimization model is effective and all the methods are well to maintain good barge floatation in the adjusting process. But the calculated amount of grid method will increase rapidly with the increasing of the number of ballast tanks adjusted or the number of adjusting time's division, and it maybe lead to an unacceptable solution time. Although penalty function method has advantage of short solution time, the barge floatation is not good enough relatively during the adjusting process. While the barge floatation will be the best if simulated annealing method is adopted and just the computing time is slightly longer. So it can draw the conclusion that penalty function method can be used under the case of computing time demanded shortly and floatation limited loosely, and vice versa, grid method or simulated annealing should be adopted.
引文
[1]李沁溢.国外先进船体建造技术现状和发展趋向.造船技术,2007,5:9-12,15.
[2]造船工程师.偶工作的单位[EB/OL].http://www.shipfriends.net/bbs,2007-09-08.
[3]梅孝恒,周延东,樊之夏.带有间断不对称式滑靴的大型结构物拖拉滑移装船技术.中国海上油气(工程),1997,9(2):1-7.
[4]王照宁.集装箱船舶自动配载系统管理模型的建立与研究.大连海运学院学报,1994,20(2):70-84.
[5]杨星,季永青.集装箱船配载优化方法研究.武汉理工大学学报,2002,2:223-226.
[6]张维英.集装箱船全航线配载智能优化研究:(博士学位论文).大连:大连理工大学,2005.
[7]车鉴.集装箱船舶配载计划优化研究:(硕士学位论文).上海:上海交通大学,2008.
[8]杜嘉立,杨盐生,张义军.基于Visual C++散货船自动配载仪的开发.大连海事大学学报,2004,30(4):9-12.
[9]杨军.散货船自动配载系统的研究与实现:(硕士学位论文).武汉:武汉理工大学,2004.
[10]杨久青.计算机辅助油轮配载.大连海运学院学报,1992,18(4):401-408.
[11]郭昌捷,马骏.油船配载控制与剩余强度研究.大连理工大学学报,1995,35(6):867-872.
[12]周丰.船舶配载的计算机算法改进.武汉交通职业学院学报,2004,6(1):77-79.
[13]成鹏飞,刘慧芬.计算机优化物流环节在船舶配载上的应用.矿业研究与开发,2002,22(5):41-44.
[14]刘勇,徐言民.船舶配积载模糊优化专家系统.船海工程,2002,150(6):33-35.
[15]简兆权,丁承民,褚启勤.基于神经网络专家系统的货船配载系统的设计与实现.交通与计算机,1997,15(1):5-9.
[16]孙静,朱绍庐,刘世宁.船舶最佳纵倾及压载微机系统.大连海运学院学报,1988,14(2):48-58.
[17]邱文昌,吴善刚.对船舶最佳压载方案编制问题的研究.上海海运学院学报,2001,22(2):4-8.
[18]王鸿鹏.一种优化编制船舶装载计划的新方法.集美大学学报(自然科学版),2002,7(1):38-44.
[19]张俊,王忠煜,曹志英.滚装船舶车辆优化配载系统得研究与实现.交通与计算机,2000,18(增刊):43-46.
[20]宋振洪.客滚船舶车辆的动态装载算法.计算机工程,2004,30(增刊):522-523.
[21]张明霞,夏益美,林焰,肖进军.全回转起重船作业过程稳性仿真.船舶,2007,36(4):137-141.
[22]段兴锋.基于型值表的船舶装载仪的研究:(硕士学位论文).大连:大连海事大学,2008.
[23]吴祥虎,程远胜,刘均.船舶纵向下水配载优化研究.中国舰船研究,2006,1(2):23-25,34.
[24]林焰,蔡尊德.AutoCAD2002环境下驳船自动配载研究.大连理工大学学报,2005,45(5):676-681.
[25]杨立,林焰,纪卓尚.基于总纵弯矩的驳船配载优化设计.船舶,2004,4:24-27.
[26]张权松.船用动载横倾平衡系统.柴油机,2005,27(3):36-37.
[27]胡思纯,黄大乐.沉箱出运过程的浮船坞调水技术及量化分析.华南港工,2008,110(1):45-51.
[28]赵龙星,陈建,陈维军,杨博良.特殊用途浮船坞作业控制系统研究.上海船舶运输科学研究所学报,2006,29(1):26-28.
[29]宗诚,陈远猷,陈逢源.浮船坞的实时控制.船舶设计通讯,1997,95(3):36-42.
[30]陈谦,吴宏.浮船坞沉浮的智能控制.船舶,2007,1:43-44.
[31]Shields J J.Containership stowage:a computer-aided preplanning system.Marine Technology,1984,21(4):370-383.
[32]Martin G L,Randhawa S U,McDowell E D.Computerised Containership Load Planning:A Methodology and Evaluation.Computers&Industrial Engineering,1988,14(4):429-440.
[33]Botter R C,Brinati M.A.Stowage Container Planning:A Model for Getting an Optimal Solution.Computer Applications in the Automation of Shipyard Operation and Ship Design,IFIP Transactions B(Application in Tech.),1992:217-229.
[34]Wilson I D,Roach P A.The deep-sea container-ships stowage problem modeling and automating the human planning processin:Proceedings of the Second International ICSC Symposium on Intelligent Industrial Automation,ICSC,Canada/Switzerland,1997,129-135.
[35]John J.D.Expert System Applications to Ocean Shipping-A Status Report.Marine Technology,1990,27(5):265-284.
[36]Dubrovsky O,Levitin G,Penn M.A Genetic Algorithm with a Compact Solution Encoding for the Container Ship Stowage Problem.Journal of Heuristics,2002,8:585-599.
[37]Sato K,Itoh H,Awashima Y.Expert System for Oil Tanker Loading/Unloading operation planning.Computer Applications in the Automation of Shipyard Operation and Ship design,IFIP.North Holland,1992.
[38]Subrata K.Chakrabarti.Handbook of offshore engineering.Illinois:Elsevier,2005.
[39]Leland G.Cole.The Treatment Of Ship Bilge,Ballast And Dock Rain-water Run-off By An Electrox Process.In:IEEE Ocean Engineering,1996,19(8):72-74.
[40]Hsien Hua,Lee Pei-Wen Wang.Analytical solution on the surge motion of tension-leg twin platform structural systems.In:IEEE Ocean Engineering,2000,27:393-415.
[41]Ore Arup.Solution for the shallows.In:Offshore engineering,1998,1:42-43.
[42]Neil Ferguson,Hugo Zarate,Toyoharu Kitani,Osamu Inokoshi,Satoshi Masuda.An Analytical Study and Systematic Monitoring Procedure Developed for the Load-Out Operation of the North Rankin Jacket 'A'.In:Proceedings of 15th annual offshore technology conference(held in Houston,Texas),1983:497-504.
[43]E.S.Piter,K.A.Digre,R.M.Tabone.Bullwinkle Loadout Analysis.In:Proceedings of 21st annual offshore technology conference(held in Houston,Texas),1989:469-478.
[44]J.E.Hofferber,N.A.Macleod,T.G.McGarragh,J.F.Thibeaux.Loadout,Transportation,and Installation of the Harmony and Heritage Jackets.In:Proceedings of 23rd annual offshore technology conference(held in Houston,Texas),1991:371-380.
[45]Y.T.Yang,B.N.Park,S.S.Ha.Development of load-out design methodology and numerical strength evaluation for on-ground-build floating storage and offloading system.Ocean Engineering,2005,32:986-1014.
[46]A.Kurniawan,Guowei Ma.Optimization of ballast plan in launch jacket load-out.Structural and Multidisciplinary Optimization,2008,online first.
[47]刘亚东,高国庆.气囊搬运技术在码头大型沉箱施工的应用.福建建筑,2005,94(4):73-74.
[48]冯章杰,陆克发.气囊“滚筒”拖运技术在港口工程中的应用.珠江现代建设,2006,131(3):7-11.
[49]韩邦峰,滕涛.半潜驳不搭岸出运沉箱应用技术.水运工程,2007,4:93-95.
[50]桑运水,王允,郭刚.大型平台滑移装船和浮装就位技术.中国造船,2003,44(增):129-134.
[51]朱崇诚.大型海上平台滑移装船有限元分析:(硕士学位论文).天津:天津大学,2004.
[52]幺子云.大型结构物滑移装船技术的研究:(硕士学位论文).天津:天津大学,2004.
[53]翟晓岗.大型结构物滑移装船过程的优化控制:(硕士学位论文).天津:天津大学,2005.
[54]刘桂涛.基于力反馈的大型结构物滑移装船计算机检测系统研究:(硕士学位论文).天津:天津大学,2005.
[55]周延兴.大型结构物滑移装船监控系统及实验研究:(硕士学位论文).天津:天津大学,2006.
[56]张伟.大型结构物滑移装船过程仿真研究:(硕士学位论文).天津:天津大学,2006.
[57]盛振邦,杨尚荣,陈雪深.船舶静力学.上海:上海交通大学出版社,1991.
[58]王志荣,钱培英.船体及舱室处理新方法在破舱稳性计算中的应用.船舶工程,1988,5:9-14.
[59]朱照辉.基于数据库的船舶破舱稳性计算确定性方法研究:(硕士学位论文).大连:大连理工大学,1999.
[60]杨博良,钱培英“客船分舱与稳性等效规则”分析及电算方法.中国造船,1983,1:45-55.
[61]王杰德,杨永谦等.船体强度与结构设计.北京:国防工业出版社,1995.
[62]奥特加,莱因博尔特.多元非线性方程组迭代解法.北京:科学出版社,1983.
[63]孙昕.关于投影梯度法的一些新的研究成果:(硕士学位论文).大连:辽宁师范大学,2007.
[64]夏洪胜,盛昭瀚,徐南荣.多目标决策问题的交互式简约梯度法.系统工程,1992,10(1):36-41.
[65]霍水泉,李人厚.多目标规划的惩罚函数法.西安交通大学学报,1993,27(3):75-80,90.
[66]陶海,沈祥福.斜拉桥索力优化的强次可行序列二次规划法.力学学报,2006,38(3):381-384.
[67]王向慧,连志春,徐志英,唐云岚.基于Pareto最优概念的多目标进化算法研究.计算机工程与应用,2008,44(27):58-61.
[68]李香平,张红阳.模拟退火算法原理及改进.软件导刊,2008,7(4):47-48.
[69]岳晓辉.基于禁忌搜索算法的蛋白质结构预测的研究:(硕士学位论文).大连:大连理工大学,2006.
[70]唐焕文,秦学志.实用最优化方法.大连:大连理工大学出版社,2004.
[71]孙靖民.机械优化设计.北京:机械工业出版社,2000.
[72]Bertsekas D P.Multiplier Method:A Survey.Automatica,1976,12(2):133-145.
[73]史明仁.一个使用乘子法的最优化FORTRAN通用程序OPTIMA.北京工业大学学报,1984,10(3):51-60.
[74]朱谦阳.航海气象与海洋学.大连:大连海事大学出版社,1999.
[75]程远胜,曾广武.船舶坐墩反力的分析计算.华中理工大学学报,1994,22(10):68-72.
[76]程远胜,曾广武.非概率不确定性及其对船舶坐墩配墩优化的影响.工程力学,2003,20(3):129-133.
[77]张寿富,朱小龙.船舶布墩理论与实践.江苏船舶,1995,3:19-30.
[78]陈铁云.船舶结构力学.上海:上海交通大学出版社,1991.
[79](苏)Ю.A.斯曼斯基.船舶结构力学手册 第三卷.孙海涛译.上海:上海科学技术出版社,1980.
[80]卓家寿,赵宁.不连续介质静、动力分析的刚体-弹簧元法.河海大学学报,1993,21(5):34-43.
[81]刘玉擎.有限弹簧元在桥梁结构分析中的应用研究.桥梁建设,2005,6:19-22.
[82]金尚年.理论力学.北京:高等教育出版社,2002.
[83]周瑛,吴国忠,曾广杰,余飞鸿.用模拟退火算法设计光学有限脉冲响应滤波器.光学学报,2003,23(8):1000-1004.