船舶浮态计算的一种优化方法.pdfVIP

第 32 卷 　第 2 期 大 连 海 事 大 学 学 报 Vol. 32 　No. 2 2006 年 5 月 Journal of Dalian Maritime University May 　2006 文章编号 (2006) 船舶浮态计算的一种优化方法 孙承猛 ,刘寅东 (大连理工大学 船舶工程学院 ,辽宁 大连 　116024) 摘要 :求船舶空间浮态就是计算船舶在某种给定载况下的三个浮态参数 ,它归结为求解一个隐式非线性方 程组. 采用传统的逐次线性化方法迭代求解方程组计算较繁 、工作量大 ,且不易编程. 为此探讨了优化方法在 船舶空间浮态计算中的应用. 利用最小二乘法的思想 ,将原问题转化为相应的数学优化模型 ,采用惩罚函数法 对其求解. 该方法不需要计算船舶水线面面积 、漂心和惯性矩等要素 ,计算结果准确 、可靠 ,且易于编程实现. 通过算例证实了该方法可行 、适用. 关键词 :船舶静力学 ;船舶浮态 ;优化算法 ;最小二乘思想 U 661. 22 　　　　　　　　　　　　　 A 中图分类号 : 文献标识码 : 式方程组, 通 常采用逐次线性化方法迭代求 0 　引　言 解[ 1] . 该方法公式推导比较繁琐, 而且每次迭代 求船舶空间浮态就是计算船舶在某种给定载 都需要计算含有水线面面积 、漂心、惯性矩 、排水 ( ) 况下的浮态参数 本文指符拉索夫角参数 . 从本 体积及浮心等多项要素的Jacobi 矩阵, 不但计算 质上讲 ,船舶的浮态参数取决于船舶平衡方程组 工作量大, 而且不易准确求出各水线面要素, 不利 的求解. 它主要包括三个独立方程 :一个力平衡方 于编程. 利用非线性方程组的求解可化为求解一 [1 ] 程和两个力矩平衡方程. 其形式可表述为 个非线性规划[ 2] 的思想, 建立相应的优化模型, f 1 ( X) = V - V 0 = 0 采用惩罚函数法进行求解. 在每次迭代过程中, 只 f ( X) = m + m ·tan θ= 0 ( ) 2 x z xy 1 需计算水线面以下的排水体积及浮心坐标, 不需 f 3 ( X) = my z + m xy ·tan φ = 0 要计算水线面上的诸要素, 从而回避了这些要素 θ φ T 式中: X = [ dm , tan , tan ] ; 的求解准确性问题, 减少了计算量, 也易于编程实 my z = V ·X b - V 0 ·X g 现. m = V ·Y - V ·Y ( ) x z b 0 g 2 m xy = V ·Zb - V 0 ·Zg 1 　优化模型及求解方法 θφ 其中: d m 、、 分别为浮态参数, 即平均吃水 、横 1. 1 　优化模型