Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

破舱倾覆船体扳正过程数值模拟

潘德位 林成新 周兆欣 孙玉强 刘志杰

薛小荣. 公路工程项目经济评价指标计算方法[J]. 交通运输工程学报, 2008, 8(4): 58-60.
引用本文: 潘德位, 林成新, 周兆欣, 孙玉强, 刘志杰. 破舱倾覆船体扳正过程数值模拟[J]. 交通运输工程学报, 2017, 17(5): 102-112.
XUE Xiao-rong. Calculation method of economic evaluation indices for highway engineering project[J]. Journal of Traffic and Transportation Engineering, 2008, 8(4): 58-60.
Citation: PAN De-wei, LIN Cheng-xin, ZHOU Zhao-xin, SUN Yu-qiang, LIU Zhi-jie. Numerical simulation of righting process for damaged-capsized hull[J]. Journal of Traffic and Transportation Engineering, 2017, 17(5): 102-112.

破舱倾覆船体扳正过程数值模拟

基金项目: 

交通运输部建设科技项目 2013328225080

高等学校博士学科点专项科研基金项目 20122125120013

国家科技支撑计划项目 2014BAK05B06

详细信息
    作者简介:

    潘德位(1986-), 男, 辽宁大连人, 山东交通学院讲师, 工学博士, 从事船舶打捞研究

    通讯作者:

    林成新(1963-), 男, 山东栖霞人, 大连海事大学教授, 工学博士

  • 中图分类号: U698.6

Numerical simulation of righting process for damaged-capsized hull

More Information
Article Text (Baidu Translation)
  • 摘要: 考虑了破舱倾覆船体浮性和稳性, 研究了船体在扳正过程中空间位置和受力状态; 采用欧拉旋转变换方法建立了船体空间力学平衡方程, 根据船舶静力学原理, 得到了破舱倾覆船体稳性和扳正力数学模型; 根据伯努利定理计算了破舱进水量及其对船体重心和浮心位置的影响; 利用GHS软件模拟了破舱倾覆船体的扳正过程, 求解了其最大扳正力和进水量, 计算了船体纵向6个位置的剪力、弯矩和扭矩。计算结果表明: 在最初扳正时, 破舱进水导致倾覆船体扳正力矩降低了130.312 MN·m, 说明破舱进水降低了倾覆船体的稳性, 可以减小最初扳正力, 降低了扳正难度; 在扳正后期时, 破舱进水产生的倾斜力矩最大值为163.594 MN·m, 说明破舱进水降低了船体的稳性, 提高了扳正难度, 仍需要施加较大的扳正力平衡船体; 船体纵向强度分布会随着扳正力和破舱进水量的变化而改变, 多点扳正船体的最大扳正力小于单点最大扳正力的40%, 最大扭矩小于单点扭矩的50%;方案1~4的最大进水量分别为6 269.76、6 781.01、5 830.76、6 653.33t, 因此, 合理布置扳正点的位置, 单点扳正(方案1~3) 的进水量小于多点扳正(方案4)。

     

  • 公路工程项目经济评价是公路工程经济分析的核心内容, 目的在于确保决策的正确性和科学性, 避免或最大限度地减少公路工程项目投资的风险, 了解公路建设方案的经济效果, 最大限度地提高公路工程项目投资的综合经济效益[1-2]。公路工程项目投资决策正确与否在很大程度上取决于公路工程项目方案的评价结果。评价方法的研究始终是一个重点。目前, 确定条件下的经济评价方法已经相对成熟, 研究的重点集中在不确定条件下的经济评价。不确定条件主要指的是随机条件和模糊条件, 随机条件下的经济评价是在确定随机分布的情况下进行[3]。由于项目的一次性和独特性以及企业较少有自己的项目数据库, 不易获取方案现金流量的概率分布, 然而项目评估者却较善于进行模糊性的估计[4], 当方案的现金流量或参数为模糊变量时, 就很有必要探讨模糊条件下的方案评价方法, 目前对这一问题研究较少。评价方法主要涉及评价指标的计算实现问题, 基于此, 本文建立了模糊条件下公路工程项目经济评价指标体系, 利用模糊模拟技术实现指标值的计算。

    ξ是从可能性空间{ϕ, ψ(ϕ), P}到实数集R的一个函数, 则称ξ为一个模糊变量[5-6], 其中: ϕ为非空集; ψ(ϕ)为ϕ的幂集; P为可能性测度。常用的模糊变量有梯形模糊变量和三角形模糊变量。梯形模糊变量由4元组(r1, r2, r3, r4)完全决定, r1 < r2 < r3 < r4; 三角形模糊变量由3元组(r1, r2, r3)完全决定, r1 < r2 < r3ξα水平集为

    ξα={ξ(θ)θUnknown node type: i,Ρ(θ)α}

    ξ的隶属度函数为

    μ(x)=Ρ[θUnknown node type: iξ(θ)=x]xR

    模糊变量的运算[5-6]如下。

    (1) 同一可能空间上的模糊运算。设

    fRnR

    是1个函数, ξ1ξ2、…、ξn是可能性空间{ϕ, ψ(ϕ), P}上的模糊变量, 则

    ξ=f(ξ1,ξ2,,ξn)

    是1个模糊变量, 定义为

    ξ(θ)=f[ξ1(θ),ξ2(θ),,ξn(θ)]θUnknown node type: i

    (2) 不同可能性空间上的模糊运算。

    fRnR是1个函数, 并且ξi是可能性空间{ϕi, ψ(ϕi), Pi}, i=1, 2, …, n上的模糊变量, 则

    ξ=f(ξ1,ξ2,,ξn)

    是乘积可能性空间{ϕ, ψ(ϕ), P}上的模糊变量, 定义为

    ξ(θ1, θ2, …, θn)=f[ξ1(θ1), ξ2(θ2), …, ξn(θn)] θiϕ

    模糊运算后隶属度可由Zadeh扩展原理获得

    μ(x)=supx1,,xi,,xn{min1inμi(xi)x=f(x1,,xi,,xn)}

    模糊变量的比较常由模糊变量的可信性测度、关键值测度和期望值来完成[5-6]

    模糊条件下公路工程项目各方案现金流量等基础数据不确定, 无法比较大小, 无法按确定条件下的指标体系进行经济评价。模糊条件下的项目经济评价就需要寻求模糊理论的支持, 模糊理论中的模糊期望值、可信性值和关键值理论可以解决这一问题, 原因在于模糊变量可通过其期望值、可信值和关键值比较大小, 因此, 模糊条件下的经济评价就转化为计算方案的评价指标的期望值、可信值和关键值。对应于这3个计算值, 在模糊条件下, 方案经济评价时可采用的形式就有3种: 一是求某一指标的模糊期望值; 二是考察指标值大于(或小于)某一指定值的可信性; 三是计算指标值达到某一可信性的最大指标值(关键值)。对应确定条件下的经济评价指标, 模糊条件条件下的经济评价指标体系就是模糊期望值(图 1)、可信值和关键值指标体系。按照计算模糊期望值、可信值和关键值时, 是否考虑资金时间价值, 可将指标体系又可分为静态和动态评价指标。

    图  1  模糊期望值指标体系
    Figure  1.  Index system of fuzzy expectation value

    刘宝碇等提出了模糊模拟技术, 是对模糊系统模型进行抽样试验的一项技术[7-10]f(ξ)的期望值E[f(ξ)]、可信性值C{f(ξ)≤r}(r为任意指定的实数, 如要求净现值大于等于0, 则r为0)和满足C{f(ξ)ˉf}a的最大ˉf在模拟次数充分多的情况下, 可通过模糊模拟求得。其中: a为任意指定可信值; f¯为相应指标值。模拟时, 分别从ϕ中均匀产生θk, 使得

    Ρ(θk)ε

    式中: ε为充分小的正数。并记

    υk=Ρ(θk)k=1,2,n

    该过程等价于分别从ξ1ξ2、…、ξnε水平集中均匀地产生u1u2、…、un, 因此, 有

    υk=μ1(u1k)μ2(u2k)μn(unk)

    式中: μiξi的隶属度函数。

    步骤1:分别从ϕ中均匀产生θk, 使得P(θk)≥ε

    步骤2:置υkP(θk)。

    步骤3:通过估计公式返回L(r)

    L(r)=12{max1kn{vkf[ξ(θk)]r}+min1kn{1-vkf[ξ(θk)]<r}}

    步骤1:均匀产生θk, 使得P(θk)≥ε, 并记υkP(θk)。

    步骤2:找到满足L(r)≥a的最大值r

    步骤3:返回r

    模糊变量的期望值为

    E[f(ξ)]=0+C{f(ξ)r}dr--0C{f(ξ)<r}dr

    其模糊模拟的过程如下。

    步骤1:置变量e为0

    步骤2:均匀产生θk, 使得P(θk)≥ε, 并记υkP(θk)。

    步骤3:置

    a=f[ξ(θ1)]f[ξ(θn)]b=f[ξ(θ1)]f[ξ(θn)]

    步骤4:从[a, b]中均匀产生r

    步骤5:如果r≥0, 那么

    ee+C{f(ξ)r}

    如果r < 0, 那么

    ee-C{f(ξ)r}

    步骤6:重复步骤4、5共n次。

    步骤7:最后得

    E[f(ξ)]=a0+b0+e(b-a)/n

    某公路工程投资项目初始投资估计为三角模糊变量(130, 140, 150), 单位为万元; 投资后当年即可获得正常收益, 年净收益为梯形模糊变量(27, 28, 29, 30), 单位为万元; 寿命期为三角形模糊变量(14, 15, 16), 单位为年; 基准收益率为三角形模糊变量(10%, 11%, 12%)。投资、年净收益、寿命期和基准收益率均为模糊变量, 分别设为ξ1ξ2ξ3ξ4, 则净现值为

    f(ξ)=f(ξ1,ξ2,ξ3,ξ4)=-ξ1+ξ2(1+ξ4)ξ3-1ξ4(1+ξ4)ξ3

    也为模糊变量。采用前面介绍的模糊模拟技术, 通过5 000次的模糊模拟计算, 求得E[f(ξ)]为657 794.7元, C{f(ξ)≥0}为100%, 满足

    C{f(ξ)f¯}=0.9

    的最大f¯值为373 229.1元。工程项目的E[f(ξ)] > 0, 说明方案可行; 方案的可信性C{f(ξ)≥0}=100%, 说明方案一定可行; 方案有90%的把握能获利373 229.1元。

    本文建立了模糊条件下的经济评价指标体系, 提出了评价指标值计算方法, 并利用指标值对模糊条件下的经济项目进行了评价。评价结果客观, 易于解释, 对历史数据要求不高。但是案例中的三角模糊数是在评估人员直接给出的情况下获得, 未来还需要进一步分析在有一定的历史数据时, 如何通过模糊统计分析来获取方案数据的模糊统计估计值。

  • 图  1  坐标系

    Figure  1.  Coordinate systems

    图  2  破舱进水

    Figure  2.  Damaged cabin bilging

    图  3  双体船

    Figure  3.  Catamaran

    图  4  船体和舱室

    Figure  4.  Hull and cabins

    图  5  倾覆船体

    Figure  5.  Capsized hull

    图  6  倾覆船体的稳性力矩曲线

    Figure  6.  Righting moment curves of capsized hull

    图  7  倾覆船体扳正过程中进水量

    Figure  7.  Flooding quantities of capsized hull during uprighting

    图  8  倾覆船体剪力曲线

    Figure  8.  Shear curves of capsized hull

    图  9  倾覆船体弯矩曲线

    Figure  9.  Bending moment curves of capsized hull

    图  10  倾覆船体扭矩曲线

    Figure  10.  Torque curves of capsized hull

    表  1  船体主尺度

    Table  1.   Principal dimensions of hull

    下载: 导出CSV
  • [1] BAČKALOV I, BULIAN G, ROSĚN A, et al. Improvement of ship stability and safety in intact condition through operational measures: challenges and opportunities[J]. Ocean Engineering, 2016, 120: 353-361. doi: 10.1016/j.oceaneng.2016.02.011
    [2] BAČKALOV I. Impact of contemporary ship stability regulations on safety of shallow-draught inland container vessels[J]. Safety Science, 2015, 72: 105-115. doi: 10.1016/j.ssci.2014.09.001
    [3] BOULOUGOURIS E, PAPANIKOLAOU A. Risk-based design of naval combatants[J]. Ocean Engineering, 2013, 65: 49-61. doi: 10.1016/j.oceaneng.2013.02.014
    [4] ELEFTHERIA E, APOSTOLOS P, MARKOS V. Statistical analysis of ship accidents and review of safety level[J]. Safety Science, 2016, 85: 282-292. doi: 10.1016/j.ssci.2016.02.001
    [5] GOERLANDT F, KUJALA P. Traffic simulation based ship collision probability modeling[J]. Reliability Engineering and System Safety, 2011, 96 (1): 91-107. doi: 10.1016/j.ress.2010.09.003
    [6] 赵晓非. 船舶分舱和破舱稳性计算[J]. 大连理工大学学报, 1988, 28 (1): 87-92. https://www.cnki.com.cn/Article/CJFDTOTAL-DLLG198803017.htm

    ZHAO Xiao-fei. The calculation of subdivision and damage stability of ship[J]. Journal of Dalian University of Technology, 1988, 28 (1): 87-92. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-DLLG198803017.htm
    [7] 林焰, 李铁骊, 纪卓尚. 破损船舶自由浮态计算[J]. 大连理工大学学报, 2001, 14 (1): 85-89. doi: 10.3321/j.issn:1000-8608.2001.01.019

    LIN Yan, LI Tie-li, JI Zhuo-shang. Ship damage floating calculation[J]. Journal of Dalian University of Technology, 2001, 14 (1): 85-89. (in Chinese). doi: 10.3321/j.issn:1000-8608.2001.01.019
    [8] 李玉明. 运木船破舱浮性、稳性、强度及破舱进水模拟研究[D]. 大连: 大连海事大学, 2002.

    LI Yu-ming. Study on the buoyancy, stability, longitudinal strength of lumber carrier and simulation of damaged flooding[D]. Dalian: Dalian Maritime University, 2002. (in Chinese).
    [9] RUPONEN P. Adaptive time step in simulation of progressive flooding[J]. Ocean Engineering, 2014, 78: 35-44. doi: 10.1016/j.oceaneng.2013.12.014
    [10] RUPONEN P. Progressive flooding of a damaged passenger ship[D]. Helsinki: Helsinki University of Technology, 2007.
    [11] VERMEER H, VREDEVELDT A W, JOURNÉE J M J. Mathematical modelling of motions and damaged stability of ro-ro ships in the intermediate stages of flooding[C]//Trove. Fifth International Conference on Stability of Ships and Ocean Structures. Canberra: Trove, 1994: 1-9.
    [12] VREDEVELDT A W, JOURNÉE J M J. Roll motions of ships due to sudden water ingress, calculations and experiments[C]//British Maritime Technology. International Conference on Ro-Ro Safety and Vulnerability the Way Ahead. London: British Maritime Technology, 1991: 1-17.
    [13] DOMEH V D K, SOBEY A J, HUDSON D A. A preliminary experimental investigation into the influence of compartment permeability on damaged ship response in waves[J]. Applied Ocean Research, 2015, 52: 27-36. doi: 10.1016/j.apor.2015.05.001
    [14] MANDERBACKA T, RUPONEN P. The impact of the inflow momentum on the transient roll response of a damaged ship[J]. Ocean Engineering, 2016, 120: 346-352. doi: 10.1016/j.oceaneng.2016.02.012
    [15] RODRIGUES J M, SOARES C G. A generalized adaptive mesh pressure integration technique applied to progressive flooding of floating bodies in still water[J]. Ocean Engineering, 2015, 110: 140-151. doi: 10.1016/j.oceaneng.2015.10.002
    [16] YAO T. Hull girder strength[J]. Marine Structures, 2003, 16 (1): 1-13. doi: 10.1016/S0951-8339(02)00052-7
    [17] IVANOV L D. On the relationship between maximum still water shear forces, bending moments, and radii of gyration of the total ship's weight and buoyancy forces[J]. Ships and Offshore Structures, 2007, 2 (1): 39-47. doi: 10.1533/saos.2006.0147
    [18] ŽAJA D, ZAMARIN A, HADJINA M. Longitudinal strength of a container ship[J]. Engineering Review, 2007, 27 (1): 55-66.
    [19] KHAN I A, DAS P K. Reliability analysis of intact and damaged ships considering combined vertical and horizontal bending moments[J]. Ships and Offshore Structures, 2008, 3 (4): 371-384. doi: 10.1080/17445300802369950
    [20] PAIK J K, THAYAMBALLI A K, PEDERSEN P T, et al. Ultimate strength of ship hulls under torsion[J]. Ocean Engineering, 2001, 28 (8): 1097-1133. doi: 10.1016/S0029-8018(01)00015-4
    [21] MOHAMMED E A, BENSON S D, HIRDARIS S E, et al. Design safety margin of a 10 000TEU container ship through ultimate hull girder load combination analysis[J]. Marine Structures, 2016, 46: 78-101. doi: 10.1016/j.marstruc.2015.12.003
    [22] ZHAO Meng-xin. Discussion on uprighting big-angle tilted sunken ships underwater in salvage operation[J]. China Ocean Engineering, 1989, 3 (2): 217-227.
    [23] 赵孟信. 打捞工程中大倾度沉船水下扶正的探讨[J]. 海洋工程, 1988 (1): 11-18. https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC198801001.htm

    ZHAO Meng-xin. The research of uprighting a sunken ship inclined at a big angle in salvage operation[J]. The Ocean Engineering, 1988 (1): 11-18. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HYGC198801001.htm
    [24] 赵孟信. "长征"轮打捞中关于扳正力的分析[J]. 航海技术, 1995 (2): 16-18, 20. https://www.cnki.com.cn/Article/CJFDTOTAL-HHJS199502005.htm

    ZHAO Meng-xin. Righting force analysis in Changzheng salvaging[J]. Marine Technology, 1995 (2): 16-18, 20. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HHJS199502005.htm
    [25] DROBYSHEVSKI Y. A note on uprighting of a ship floating upside-down[J]. Ocean Engineering, 2004, 31 (11/12): 1447-1467.
    [26] ANTONINI A, LAMBERTI A, CECCARELLI G. What could happen if the parbuckling of Costa Concordia had failed: analytical and CFD-based investigation of possible generated wave[C]//LYNETT. Proceedings of 34th Conference on Coastal Engineering. Reston: Coastal Engineering Proceedings, 2014: 55-60.
    [27] DRUCE V. Costa Concordia's seafloor legacy[J]. New Scientist, 2013, 219 (2935): 6.
    [28] 刘汉明, 顾泽月. 救捞作业场辅助设计软件的设计与实现[C]//中国航海学会. 救捞专业委员会2004年学术交流会论文集. 北京: 中国航海学会, 2004: 65-70.

    LIU Han-ming, GU Ze-yue. The design and realization of the auxiliary design software of layout for salvage boat[C]//China Institute of Navigation. Proceedings of Rescue and Salvage Conference 2004. Beijing: China Institute of Navigation, 2004: 65-70. (in Chinese).
    [29] 潘德位. 倾覆船舶扳正过程中的受力分析与计算[D]. 大连: 大连海事大学, 2017.

    PAN De-wei. Force analysis and calculation for righting a capsized ship[D]. Dalian: Dalian Maritime University, 2017. (in Chinese).
    [30] 马振伟. 静水中船舶破损稳性计算研究[D]. 武汉: 武汉交通科技大学, 1999.

    MA Zhen-wei. Calculation of damaged ship on still water[D]. Wuhan: Wuhan Transportation Technology University, 1999. (in Chinese).
    [31] PÉREZ F, CLEMENTE J A. Constrained design of simple ship hulls with B-spline surfaces[J]. Computer-Aided Design, 2011, 43 (12): 1829-1840. doi: 10.1016/j.cad.2011.07.008
    [32] PÉREZ-ARROABS F. Parametric generation of planing hulls[J]. Ocean Engineering, 2014, 81: 89-104.
    [33] CUI Hao, TURAN O, SAYER P. Learning-based ship design optimization approach[J]. Computer-Aided Design, 2012, 44 (3): 186-195.
  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  733
  • HTML全文浏览量:  186
  • PDF下载量:  392
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-06-10
  • 刊出日期:  2017-10-25

目录

/

返回文章
返回