留言板

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

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

基于动力效应的船体远场冰载荷测量与识别

何帅康 陈晓东 孔帅 季顺迎

何帅康, 陈晓东, 孔帅, 等. 基于动力效应的船体远场冰载荷测量与识别[J]. 中国舰船研究, 2021, 16(X): 1–10 doi:  10.19693/j.issn.1673-3185.02065
引用本文: 何帅康, 陈晓东, 孔帅, 等. 基于动力效应的船体远场冰载荷测量与识别[J]. 中国舰船研究, 2021, 16(X): 1–10 doi:  10.19693/j.issn.1673-3185.02065
HE S K, CHEN X D, KONG S, et al. Measurement and identification of ice loads on hull structures in far field based on dynamic effects[J]. Chinese Journal of Ship Research, 2021, 16(0): 1–10 doi:  10.19693/j.issn.1673-3185.02065
Citation: HE S K, CHEN X D, KONG S, et al. Measurement and identification of ice loads on hull structures in far field based on dynamic effects[J]. Chinese Journal of Ship Research, 2021, 16(0): 1–10 doi:  10.19693/j.issn.1673-3185.02065

基于动力效应的船体远场冰载荷测量与识别

doi: 10.19693/j.issn.1673-3185.02065
基金项目: 国家重点研发计划资助项目(2017YFE0111400, 2016YFC1401505);国家自然科学基金资助项目(51639004, 41576179)
详细信息
    作者简介:

    何帅康,男,1997年生,硕士生。研究方向:船舶结构冰载荷测量与识别。E-mail: 408936364@mail.dlut.edu.cn

    陈晓东,男,1987年生,博士,工程师。研究方向:海冰力学与热力学。E-mail: chenxiaodong@dlut.edu.cn

    孔帅,男,1990年生,博士,工程师。研究方向:船舶结构冰载荷的监测反演和数值研究。E-mail: 764116421@qq.com

    季顺迎,男,1972年生,教授,博士,博士生导师。研究方向:极地船舶与海洋工程。E-mail: jisy@dlut.edu.cn

    通讯作者:

    季顺迎

  • 中图分类号: U663.2

Measurement and identification of ice loads on hull structures in far field based on dynamic effects

  • 摘要:   目的  提出一种基于动力效应的冰载荷反演识别方法。  方法  针对船−冰相互作用过程和结构响应特性之间的卷积关系,使用格林函数建立描述冰载荷的卷积积分方程。针对“天恩”号极地运输船的北极航行,安装应变传感器测量船艏区域的冰激响应。对于求解矩阵的病态性与测量中的随机噪声,引入正则化方法解决识别计算中的不适定问题,并获取动冰载荷时程的数值近似稳定解。依托“天恩”号极地运输船的北极航行实测数据,通过船体在海冰作用下的应变响应对冰载荷进行反演识别。  结果  结果表明,格林函数能够有效识别远场冰载荷的时间历程。  结论  载荷识别的有效性受应变信号强度的影响较为明显。此外,在不同海冰破坏模式下,浮冰尺寸和船−冰接触面积是影响冰载荷特征的关键因素。
  • 图  1  “天恩”号极地运输船的航行路线及沿途冰况

    Figure  1.  The arctic route of ice class vessel Tian En and ice conditions

    图  2  “天恩”号极地运输船的监测区域及测点信息

    Figure  2.  The monitoring area and gauge map of ice class vessel Tian En

    图  3  脉冲函数表示的载荷时程曲线

    Figure  3.  The load history represented by the impulse function

    图  4  冰载荷识别的流程图

    Figure  4.  The flow chart of ice load identification

    图  5  “天恩”号测试区域的有限元模型

    Figure  5.  The FE model of monitored area of Tian En

    图  6  数值计算的测点、加载区域示意图以及单位冲击载荷下的格林函数响应

    Figure  6.  The sketches of loading area, gauging matrix and the response of Green's function under specific surface loading in the numerical simulation

    图  7  数值算例的计算结果

    Figure  7.  The results of the numerical simulation

    图  8  “天恩”号船体与浮冰的相互作用

    Figure  8.  The interaction between Tian En and the ice floe

    图  9  不同工况下的肋骨剪切应变差值测量结果

    Figure  9.  The measured results of frame shear strain difference under various navigation conditions

    图  10  不同工况的测点及子加载区域设置

    Figure  10.  The monitored area and gauge map under different conditions

    图  11  不同工况下的冰载荷识别结果

    Figure  11.  The identified ice load under different conditions

    图  12  不同工况下的最大冰载荷时程曲线与测点的最大应变响应时程曲线

    Figure  12.  The time-history curves of maximum ice load and maximum strain response of measuring point under different conditions

    图  13  不同工况下的总冰载荷时程曲线

    Figure  13.  The time history curve of global ice load under different conditions

    表  1  “天恩”号极地运输船船体主要参数

    Table  1.   Main dimensions of ice class vessel Tian En

    船长/m型宽/m型深/m设计吃水/m设计航速/kn载重量/t
    189.9928.515.81114.837 124
    下载: 导出CSV

    表  2  “天恩”号极地运输船的有限元模型主要参数

    Table  2.   Key parameters of the FE model of ice class vessel Tian En

    参数数值
    弹性模量E/GPa209
    泊松比λ0.3
    密度ρ/(kg·m−3)7 850
    下载: 导出CSV

    表  3  计算工况主要参数

    Table  3.   The main parameters of navigation conditions

    编号日期UTC
    时间
    航速
    /(m·s–1)
    吃水
    /m
    冰状冰厚
    /m
    海冰
    破坏模式
    工况12019-08-0212:025.146.5小冰盘1.2局部挤压、
    整体劈裂
    工况22019-08-033:084.246.5冰块1.5翻转滑移
    下载: 导出CSV
  • 高世龙, 刘加钊, 张晓. 冰上丝绸之路背景下商船北极航行的经济性评析[J]. 对外经贸实务, 2018(1): 26–29. doi:  10.3969/j.issn.1003-5559.2018.01.006

    GAO S L, LIU J Z, ZHANG X. The economic evaluation of Arctic shipping under the background of the Polar Silk Road[J]. Practice in Foreign Economic Relations and Trade, 2018(1): 26–29 (in Chinese). doi:  10.3969/j.issn.1003-5559.2018.01.006
    王超, 封振, 李兴, 等. 航行于碎冰区船舶冰阻力与冰响应探析[J]. 中国舰船研究, 2018, 13(1): 73–78. doi:  10.3969/j.issn.1673-3185.2018.01.011

    WANG C, FENG Z, LI X, et al. Analysis on ice resistance and ice response of ships sailing in brash ice[J]. Chinese Journal of Ship Research, 2018, 13(1): 73–78 (in Chinese). doi:  10.3969/j.issn.1673-3185.2018.01.011
    蔡伟, 朱凌. 冰碰载荷下船体板弹塑性动力响应分析[J]. 中国舰船研究, 2020, 15(6): 155–161. doi:  10.19693/j.issn.1673-3185.01804

    CAI W, ZHU L. Elastic-plastic responses of ship plates under ice impact loading[J]. Chinese Journal of Ship Research, 2020, 15(6): 155–161 (in Chinese). doi:  10.19693/j.issn.1673-3185.01804
    刘璐, 胡冰, 季顺迎. 破冰船引航下极地船舶结构冰荷载的离散元分析[J]. 水利水运工程学报, 2020(3): 11–18.

    LIU L, HU B, JI S Y. Discrete element analysis of ice loads on polar ships under pilotage of icebreaker[J]. Hydro-Science and Engineering, 2020(3): 11–18 (in Chinese).
    韩端锋, 乔岳, 薛彦卓, 等. 冰区航行船舶冰阻力研究方法综述[J]. 船舶力学, 2017, 21(8): 1041–1054. doi:  10.3969/j.issn.1007-7294.2017.08.014

    HAN D F, QIAO Y, XUE Y Z, et al. A review of ice resistance research methods for ice-going ships[J]. Journal of Ship Mechanics, 2017, 21(8): 1041–1054 (in Chinese). doi:  10.3969/j.issn.1007-7294.2017.08.014
    LEIRA B, BØRSHEIM L, ESPELAND Ø, et al. Ice-load estimation for a ship hull based on continuous response monitoring[J]. Proceedings of the Institution of Mechanical Engineers, Part M: Journal of Engineering for the Maritime Environment, 2009, 223(4): 529–540. doi:  10.1243/14750902JEME141
    SUYUTHI A, LEIRA B J, RISKA K. A generalized probabilistic model of ice load peaks on ship hulls in broken-ice fields[J]. Cold Regions Science and Technology, 2014, 97: 7–20. doi:  10.1016/j.coldregions.2013.09.012
    SUOMINEN M, KARHUNEN J, BEKKER A, et al. Full-scale measurements on-board PSRV S. A. Agulhas II in the Baltic Sea[C]//POAC '13: Proceedings of the 22nd International Conference on Port and Ocean Engineering under Arctic Conditions. Espoo, Finland: POAC, 2013.
    SUOMINEN M, KUJALA P, ROMANOFF J, et al. Influence of load length on short-term ice load statistics in full-scale[J]. Marine Structures, 2017, 52: 153–172. doi:  10.1016/j.marstruc.2016.12.006
    KOTILAINEN M, VANHATALO J, SUOMINEN M, et al. Predicting ice-induced load amplitudes on ship bow conditional on ice thickness and ship speed in the Baltic Sea[J]. Cold Regions Science and Technology, 2017, 135: 116–126. doi:  10.1016/j.coldregions.2016.12.006
    FREDERKING R, JOHNSTON M. Comparison of local ice pressures on the CCGS Terry Fox with other data[C]//The Eighteenth International Offshore and Polar Engineering Conference. Vancouver, Canada: ISOPE, 2008.
    LEE T K, LEE J H, KIM H, et al. Field measurement of local ice pressures on the ARAON in the Beaufort Sea[J]. International Journal of Naval Architecture and Ocean Engineering, 2014, 6(4): 788–799.
    KWON Y H, LEE T K, CHOI K. A study on measurements of local ice pressure for ice breaking research vessel "ARAON" at the Amundsen Sea[J]. International Journal of Naval Architecture and Ocean Engineering, 2015, 7(3): 490–499. doi:  10.1515/ijnaoe-2015-0035
    UTO S, OKA S, MURAKAMI C, et al. Ice load exerted on the hull of icebreaker PM Teshio in the south sea of okhotsk[C]//POAC '05: Proceedings of the 18th International Conference on Port and Ocean Engineering under Arctic Conditions. Potsdam, NY, USA: POAC, 2005.
    季顺迎, 雷瑞波, 李春花, 等. “雪龙”号科考船在冰区航行的船体振动测量研究[J]. 极地研究, 2017, 29(4): 427–435.

    JI S Y, LEI L B, LI C H, et al. Measurement of ice-induced local vibration of R/V Xuelong icebreaker during its navigation in ice-covered fields[J]. Chinese Journal of Polar Research, 2017, 29(4): 427–435 (in Chinese).
    王键伟, 段庆林, 季顺迎. 冰区航行中船舶结构冰载荷的现场测量与反演方法研究进展[J].力学进展, 2020, 50: 202003.

    WANG J W, DUAN Q L, JI S Y. Research progress of field measurements and inversion methods of ice loads on ship structure during ice navigation[J]. Advances in Mechanics, 2020, 50: 202003(in Chinese).
    刘瀛昊, 佟福山, 高良田. 基于原型测量的极地航行船舶船体冰载荷分析[J]. 振动与冲击, 2017, 36(7): 226–233.

    LIU Y H, TONG F S, GAO L T. Ice-induced load analysis for hull of an ice-going vessel based on full-scale measurement[J]. Journal of Vibration and Shock, 2017, 36(7): 226–233 (in Chinese).
    孔帅, 崔洪宇, 季顺迎. 船体结构冰载荷反演方法及试验验证[J]. 中国机械工程, 2020, 31(3): 281–288. doi:  10.3969/j.issn.1004-132X.2020.03.005

    KONG S, CUI H Y, JI S Y. Ice load identification method of ship structures and experimental verification[J]. China Mechanical Engineering, 2020, 31(3): 281–288 (in Chinese). doi:  10.3969/j.issn.1004-132X.2020.03.005
    崔洪宇, 胡大士, 孔帅, 等. 基于正则化方法的雪龙号破冰船冰载荷反演的研究[J]. 中国造船, 2020, 61(1): 109–119. doi:  10.3969/j.issn.1000-4882.2020.01.011

    CUI H Y, HU D S, KONG S, et al. Study on inversion of ice load for Xue Long icebreaker based on regularization method[J]. Shipbuilding of China, 2020, 61(1): 109–119 (in Chinese). doi:  10.3969/j.issn.1000-4882.2020.01.011
    吴坤铭. 应变片布置及组桥方法在工程中的应用研究[J]. 安徽建筑大学学报, 2019, 27(3): 1–6.

    WU K M. Application research of strain gage layout and bridging method in engineering[J]. Journal of Anhui Institute of Architecture & Industry, 2019, 27(3): 1–6 (in Chinese).
    刘杰. 动态载荷识别的计算反求技术研究[D]. 长沙: 湖南大学, 2011.

    LIU J. Research on computational inverse techniques in dynamic load identification[D]. Changsha: Hunan University, 2011 (in Chinese).
    韩旭, 刘杰, 李伟杰, 等. 时域内多源动态载荷的一种计算反求技术[J]. 力学学报, 2009, 41(4): 595–602. doi:  10.3321/j.issn:0459-1879.2009.04.017

    HAN X, LIU J, LI W J, et al. A computational inverse technique for reconstruction of multisource loads in time domain[J]. Chinese Journal of Theoretical and Applied Mechanics, 2009, 41(4): 595–602 (in Chinese). doi:  10.3321/j.issn:0459-1879.2009.04.017
    孙兴盛, 刘杰, 丁飞, 等. 基于矩阵摄动的随机结构动态载荷识别技术[J]. 机械工程学报, 2014, 50(13): 148–156. doi:  10.3901/JME.2014.13.148

    SUN X S, LIU J, DING F, et al. Identification method of dynamic loads for stochastic structures based on matrix perturbation theory[J]. Journal of Mechanical Engineering, 2014, 50(13): 148–156 (in Chinese). doi:  10.3901/JME.2014.13.148
    王福昌, 胡顺田. 求解病态线性方程组的模拟退火算法[J]. 大学数学, 2009, 25(4): 69–72. doi:  10.3969/j.issn.1672-1454.2009.04.014

    WANG F C, HU S T. Simulated annealing algorithm for solving III-conditioned linear systems[J]. College Mathematics, 2009, 25(4): 69–72. doi:  10.3969/j.issn.1672-1454.2009.04.014
    卢立勤, 乔百杰, 张兴武, 等. 共轭梯度最小二乘迭代正则化算法在冲击载荷识别中的应用[J]. 振动与冲击, 2016, 35(22): 176–182, 196.

    LU L Q, QIAO B J, ZHANG X W, et al. Application of conjugate gradient least squares iteration regularization algorithm in impact load identification[J]. Journal of Vibration and Shock, 2016, 35(22): 176–182, 196 (in Chinese).
  • [1] 蔡伟, 朱凌2020: 冰碰载荷下船体板弹塑性动力响应分析 , 中国舰船研究, 15, 155-161. doi: 10.19693/j.issn.1673-3185.01804
    [2] 李想, 李红霞, 黄一2020: 核电平台连接机构设计与运动响应分析 , 中国舰船研究, 15, 152-161. doi: 10.19693/j.issn.1673-3185.01786
    [3] 刘佳仑, 杨帆, 马枫, 严新平2020: 智能船舶航行功能测试验证的方法体系 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01780
    [4] 乔丹, 马宁, 顾解忡2019: 不同堆垛模式下集装箱船风载荷特性研究 , 中国舰船研究, 14, 105-115. doi: 10.19693/j.issn.1673-3185.01410
    [5] 王超, 韩康, 汪春辉, 李兴2019: 冰区航行船舶推进器特殊性分析 , 中国舰船研究, 14, 1-7. doi: 10.19693/j.issn.1673-3185.01340
    [6] 杨鹏, 寇冠元, 朱学康, 刘成义, 高长华2019: 浮体浅水波浪载荷数值计算方法研究 , 中国舰船研究, 14, 19-26. doi: 10.19693/j.issn.1673-3185.01346
    [7] 李德聪, 段宏, 吴国民, 周心桃, 杨雄辉2018: 船内爆炸载荷特性及对舰船结构毁伤研究综述 , 中国舰船研究, 13, 7-16. doi: 10.3969/j.issn.1673-3185.2018.01.002
    [8] 叶礼裕, 王超, 郭春雨, 常欣2018: 潜艇破冰上浮近场动力学模型 , 中国舰船研究, 13, 51-59. doi: 10.3969/j.issn.1673-3185.2018.02.007
    [9] 杨鹏, 叶梦熊, 李金成, 何景异, 杜兆伟2018: 对称格林函数简化计算的数学推导 , 中国舰船研究, 13, 22-28. doi: 10.19693/j.issn.1673-3185.01345
    [10] 常欣, 李鹏, 王超, 叶礼裕2018: 冰载荷和纵倾角对螺旋桨强度的影响 , 中国舰船研究, 13, 7-15,32. doi: 10.19693/j.issn.1673-3185.01037
    [11] 李井煜, 卢晓平2017: 基于Rankine源和Kelvin源格林函数求解兴波阻力的复合算法 , 中国舰船研究, 12, 1-5. doi: 10.3969/j.issn.1673-3185.2017.06.001
    [12] 龚榆峰, 张正艺, 刘敬喜, 董问, 解德2017: 基于ABAQUS的海冰单元开发及冰载荷直接计算法 , 中国舰船研究, 12, 75-85. doi: 10.3969/j.issn.1673-3185.2017.03.011
    [13] 邓琦, 毛筱菲, 吴铭浩2016: 主体瘦长度对三体船耐波性和波浪载荷的影响 , 中国舰船研究, 11, 8-14. doi: 10.3969/j.issn.1673-3185.2016.06.002
    [14] 吕保达, 刘敬喜, 解德, 龚榆峰2015: 基于UMAT的冰—结构相互作用数值仿真 , 中国舰船研究, 10, 39-45. doi: 10.3969/j.issn.1673-3185.2015.01.006
    [15] 陈明明, 王志东, 杨爽, 吴贺贺2012: 两船在波浪载荷作用下的运动响应分析 , 中国舰船研究, 7, 24-28. doi: 10.3969/j.issn.1673-3185.2012.02.005
    [16] 王晓强, 李陈峰, 任慧龙2012: 破损舰船运动与波浪载荷预报方法 , 中国舰船研究, 7, 30-35. doi: 10.3969/j.issn.1673-3185.2012.04.006
    [17] 孙谦, 陈志坚2011: 基于实船航行试验的舰船总振动阻尼识别方法 , 中国舰船研究, 6, 34-37,44. doi: 10.3969/j.issn.1673-3185.2011.06.007
    [18] 黄晓琼, 陈立, 杨雄辉, 周心桃2009: 三体船连接桥结构波浪载荷研究 , 中国舰船研究, 4, 42-46. doi: 10.3969/j.issn.1673-3185.2009.04.009
    [19] 黄晓琼, 陈立, 杨雄辉, 周心桃2008: 三体船波浪诱导载荷计算研究 , 中国舰船研究, 3, 16-20. doi: 10.3969/j.issn.1673-3185.2008.05.004
    [20] 尹群, 仇栋熠2007: 水下爆炸载荷下舰船设备抗冲击生命力模糊综合评判方法研究 , 中国舰船研究, 2, 19-23. doi: 10.3969/j.issn.1673-3185.2007.05.005
  • 加载中
图(13) / 表 (3)
计量
  • 文章访问数:  10
  • HTML全文浏览量:  5
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-10
  • 修回日期:  2020-10-19
  • 网络出版日期:  2020-12-10

基于动力效应的船体远场冰载荷测量与识别

doi: 10.19693/j.issn.1673-3185.02065
    基金项目:  国家重点研发计划资助项目(2017YFE0111400, 2016YFC1401505);国家自然科学基金资助项目(51639004, 41576179)
    作者简介:

    何帅康,男,1997年生,硕士生。研究方向:船舶结构冰载荷测量与识别。E-mail: 408936364@mail.dlut.edu.cn

    陈晓东,男,1987年生,博士,工程师。研究方向:海冰力学与热力学。E-mail: chenxiaodong@dlut.edu.cn

    孔帅,男,1990年生,博士,工程师。研究方向:船舶结构冰载荷的监测反演和数值研究。E-mail: 764116421@qq.com

    季顺迎,男,1972年生,教授,博士,博士生导师。研究方向:极地船舶与海洋工程。E-mail: jisy@dlut.edu.cn

    通讯作者: 季顺迎
  • 中图分类号: U663.2

摘要:   目的  提出一种基于动力效应的冰载荷反演识别方法。  方法  针对船−冰相互作用过程和结构响应特性之间的卷积关系,使用格林函数建立描述冰载荷的卷积积分方程。针对“天恩”号极地运输船的北极航行,安装应变传感器测量船艏区域的冰激响应。对于求解矩阵的病态性与测量中的随机噪声,引入正则化方法解决识别计算中的不适定问题,并获取动冰载荷时程的数值近似稳定解。依托“天恩”号极地运输船的北极航行实测数据,通过船体在海冰作用下的应变响应对冰载荷进行反演识别。  结果  结果表明,格林函数能够有效识别远场冰载荷的时间历程。  结论  载荷识别的有效性受应变信号强度的影响较为明显。此外,在不同海冰破坏模式下,浮冰尺寸和船−冰接触面积是影响冰载荷特征的关键因素。

English Abstract

何帅康, 陈晓东, 孔帅, 等. 基于动力效应的船体远场冰载荷测量与识别[J]. 中国舰船研究, 2021, 16(X): 1–10 doi:  10.19693/j.issn.1673-3185.02065
引用本文: 何帅康, 陈晓东, 孔帅, 等. 基于动力效应的船体远场冰载荷测量与识别[J]. 中国舰船研究, 2021, 16(X): 1–10 doi:  10.19693/j.issn.1673-3185.02065
HE S K, CHEN X D, KONG S, et al. Measurement and identification of ice loads on hull structures in far field based on dynamic effects[J]. Chinese Journal of Ship Research, 2021, 16(0): 1–10 doi:  10.19693/j.issn.1673-3185.02065
Citation: HE S K, CHEN X D, KONG S, et al. Measurement and identification of ice loads on hull structures in far field based on dynamic effects[J]. Chinese Journal of Ship Research, 2021, 16(0): 1–10 doi:  10.19693/j.issn.1673-3185.02065
    • 全球变暖导致的北冰洋冰雪融化与北极航道的逐步开发,使商船航行于北极地区成为可能。2013年中远海运集团所属“永盛”轮在北极东北航道的首航,拉开了我国商船北极航行的序幕[1]。相比于传统航道,北极航道大幅缩短了东亚与北欧地区的海上运输距离,其商业价值不容忽视。海冰是船舶在冰区航行中所面临的主要威胁之一,船−冰相互作用可能导致船体局部破坏甚至整船倾覆,因此准确、有效地衡量海冰对船舶结构的影响是保障航行安全的重要前提。

      国内外针对冰载荷的研究主要采用理论分析、数值模拟[2]、模型试验[3]和实船走航测试等方法[4]。其中,实船走航测试能得到最为真实可靠的船体冰载荷特性[5]。早期冰载荷的确定采用布置压力盒的方式对冰载荷进行直接测量,而在冰载荷的作用下,压力盒很容易脱落或损坏,因此目前多采用测量船体冰激响应(如船体应变、加速度)的方法,间接识别冰载荷。Leira等[6]和Suyuthi等[7]针对KV Svalbard船在巴伦支海的航行,安装光纤光栅传感器测量了船−冰作用下的肋骨剪切应变;Suominen等[8-9]和Kotilainen等[10]于2012年在波罗的海对PSRV S. A. AgulhasⅡ船进行了全船冰载荷、冰激振动、噪声等的测量;Frederking等[11]于2001年在纽芬兰东北海岸对CCGS Terry Fox船进行了小冰山撞击试验,分析了各海冰参数对总冰力的影响;Lee等[12]和Kwon等[13]于2010~2017年间在楚科奇海对IBRV Araon船进行了多次局部冰压现场试验和破冰性能研究;Uto等[14]针对PM Teshio船在鄂霍次克海南部冰区的航行,进行了破冰模式与冰载荷关系的研究试验。国内的现场测量主要针对“雪龙”号破冰船的极地科考航行,进行了冰激振动加速度、应变信息的测量及分析[15]

      根据船体冰激响应反求冰载荷是冰载荷测量研究中的难点[16]。Leira等[6]根据肋骨剪应变,通过应变与应力之间的关系求出剪应力,并将剪应力沿肋骨积分求得冰载荷;Suominen等[8]考虑到冰载荷在相邻肋骨间的传递,提出基于静力平衡的影响系数矩阵法,并进行了理论推导和数值分析;Kwon等[13]针对IBRV Araon船的肋骨剪应变,分别采用直接计算和影响系数矩阵法进行局部冰压力的识别,验证了影响系数矩阵法的工程可用性。然而,无论是直接计算还是采用影响系数矩阵法,都仅在静力平衡的条件下[17]进行识别计算,而冰载荷作为典型的随机动载荷,其动力响应不应被忽略。

      本文考虑到冰载荷的动力响应,拟采用格林函数法[18-19]建立冰载荷识别模型,对“天恩”号极地运输船于2019年在北极东北航道航行的冰载荷进行反演识别,并采用正则化方法避免病态问题求解过程中的不适定问题,识别出局部冰载荷的时间历程,以此为冰区航行船舶的船体冰载荷监测提供一定参考。

    • “天恩”号是中远海运特种运输股份有限公司针对极地气候打造的首艘冰级轮,于2017年完成建造,2018年取道北极东北航道顺利完成北欧首航,其船体主要参数见表1。“天恩”号的冰级为CCS Ice Class B1级,可自由通行0.8 m厚的当年冰冰区。2019年7月20日,“天恩”号由江苏太仓港出发,取道“冰上丝绸之路”,途经长江水域、中国沿海、日本海、白令海、北冰洋、北海、丹麦海峡、波罗的海于8月17日到达目的地挪威。其中,在8月1日~8月8日期间,“天恩”号在俄罗斯破冰船的引航下顺利穿过了海冰覆盖区域,图1为“天恩”号在破冰船引航下于冰区航行的航线与实况图。

      表 1  “天恩”号极地运输船船体主要参数

      Table 1.  Main dimensions of ice class vessel Tian En

      船长/m型宽/m型深/m设计吃水/m设计航速/kn载重量/t
      189.9928.515.81114.837 124

      图  1  “天恩”号极地运输船的航行路线及沿途冰况

      Figure 1.  The arctic route of ice class vessel Tian En and ice conditions

      船舶在冰区航行时,海冰对船体造成的荷载远大于波浪、风等环境因素,因此在分析时认为船体的结构变形响应主要由海冰与船体碰撞所产生。船舶舷侧肋骨作为主要的承力构件,在船体监测中被重点关注,因此一般将肋骨剪切应变作为船体结构响应监测信号和载荷识别信号。

      由于“天恩”号的船−冰作用区所在位置处为抗横倾水舱,无法安装应变信息采集装置,因此将监测区域设置于远场的水手长储物间左舷侧外板区域(Fr231~Fr242),监测区域和船−冰作用区大体位置如图2(a)所示。“天恩”号船艏区域为进行冰区加强采用混合骨架式建造,图2(b)所示为监测区域内部实景,椭圆框内为测点布置区域,区域内加密的肋骨(每半个肋位处设置一根肋骨)均为T型钢材料,应变传感器均安装在T型材腹板处。

      图  2  “天恩”号极地运输船的监测区域及测点信息

      Figure 2.  The monitoring area and gauge map of ice class vessel Tian En

      为扩大应变监测范围,选择在每个肋位处布置测点,在11根肋骨上分结构上、下位置布置了共24个测点,其中测点1~22的布置如图2(c)所示。剩余2个测点与测点17,18在同一肋位上,但安装位置更靠近于肋骨翼板,如图2(d)所示,用于比较靠近外板与远离外板的应变信号强弱。考虑到电阻应变片传感器具有灵敏度高、尺寸小、粘结牢固、适用条件好等特点,在每根肋骨的上下2个测点处均呈$90^\circ $直角布置且沿船宽方向对称分布2个电阻应变片,应变片的安装方式及桥路连接方式分别如图2(e)图2(f)所示。

      根据应变状态理论,在yz平面(肋骨腹板所在平面)内任意方向的线应变[20]可表示为

      $$\varepsilon \left( \alpha \right) = \frac{{{\varepsilon _y} + {\varepsilon _{\textit{z}}}}}{2} + \frac{{{\varepsilon _y} + {\varepsilon _{\textit{z}}}}}{2}\cos \left( {2\alpha } \right) + \frac{1}{2}{\gamma _{y{\textit{z}}}}{\rm{sin}}\left( {2\alpha } \right)$$ (1)

      式中:$\varepsilon \left( \alpha \right)$表示离y$\alpha $角度方向的线应变;${\varepsilon _y}$${\varepsilon _{\textit{z}}}$分别表示yz方向上的正应变;${\gamma _{y{\textit{z}}}}$表示测点处的剪应变。

      $\alpha =45^\circ $$\alpha = - 45^\circ $代入式(1)得:

      $$\varepsilon \left( {45^\circ } \right) = \frac{1}{2}\left( {{\varepsilon _y} + {\varepsilon _{\textit{z}}} + {\gamma _{y{\textit{z}}}}} \right)$$ (2)
      $$\varepsilon \left( { - 45^\circ } \right) = \frac{1}{2}\left( {{\varepsilon _y} + {\varepsilon _{\textit{z}}} - {\gamma _{y{\textit{z}}}}} \right)$$ (3)

      由以上两式可得

      $${\gamma _{y{\textit{z}}}} = \varepsilon \left( {45^\circ } \right) - \varepsilon \left( { - 45^\circ } \right)$$ (4)

      因此,肋骨平面内任意一点的剪切应变可由该测点处所布置的两片垂直对称分布的应变片测量得出。从而,肋骨上下两测点的剪切应变差值可表示为

      $$\Delta \gamma = {\gamma _{{\rm{up}}}} - {\gamma _{{\rm{low}}}}$$ (5)

      式中:$\Delta \gamma $为两测点间的剪切应变差值,可用于后续的载荷识别计算;${\gamma _{{\rm{up}}}}$为肋骨结构上测点处的剪切应变;${\gamma _{{\rm{low}}}}$为肋骨结构下测点处的剪切应变。

      将两电阻应变片以图2(f)所示形式接入电路中。其中:${R_1}$${R_2}$分别代表$45^\circ $$ - 45^\circ $布置的应变片;${R_{\rm{c}}}$为固定阻值电阻;AC端接电源$E$${U_{BD}}$为输出电压即所测量的电压。应变变化导致$ {R}_{1} $$ {R}_{2} $阻值变化,$ {U}_{BD} $$ \varepsilon \left(45°\right)-\varepsilon \left(-45°\right) $,即${\gamma _{yz}} $,呈线性关系,因此通过测量电路两端电压值可相应转化为肋骨处的剪切应变。

    • 根据线性时不变系统的叠加原理可知,任意动态载荷引起的应变响应可用一系列单位脉冲载荷的响应叠加而成[21-22],如图3所示。因此,作用在船体结构上的冰载荷和冰激应变响应之间可以表示为卷积积分的形式:

      图  3  脉冲函数表示的载荷时程曲线

      Figure 3.  The load history represented by the impulse function

      $$\varepsilon \left( t \right) = \int_0^t {g\left( {t - \tau } \right)f\left( \tau \right){\rm{d}}\tau } $$ (6)

      式中:$\varepsilon \left( t \right)$为冰载荷下的应变响应时间历程,本文采用肋骨结构上、下测点间的剪切应变差值时程${\rm{\Delta }}\gamma $进行计算;$g\left( t \right)$为作用点载荷与对应测点响应之间的格林函数,用于表征结构的动态特性;$f\left( t \right)$为冰载荷时间历程。

      在时域内,设置$\Delta t$为采样时间间隔,m为采样点数,可将式(6)的结构动力响应卷积分关系式离散化[23]

      $$ \left[\begin{array}{c}{\varepsilon }_{1}\\ {\varepsilon }_{2}\\ \vdots \\ {\varepsilon }_{m}\end{array}\right]=\left[\begin{array}{c}{g}_{1}\\ {g}_{2}\\ \vdots \\ {g}_{m}\end{array} \begin{array}{c}0\\ {g}_{1}\\ \vdots \\ {g}_{m-1}\end{array} \begin{array}{c}\cdots \\ \cdots \\ \ddots \\ \cdots \end{array} \begin{array}{c}0\\ 0\\ \vdots \\ {g}_{1}\end{array}\right]\left[\begin{array}{c}{f}_{1}\\ {f}_{2}\\ \vdots \\ {f}_{m}\end{array}\right]\Delta t$$ (7)

      式中:${\varepsilon _i}$$t = i \cdot \Delta t$时刻下实测的结构响应值;${g_i}$$t = i \cdot \Delta t$时刻下的格林函数值;${f_i}$$t = i \cdot \Delta t$时刻下的待识别载荷值。

      冰与船的撞击通常表现出多源载荷作用的形式,而其中单个载荷引起的响应依然可以用式(7)表示,结构总的响应则是各单源载荷引起响应的线性叠加,因此多源载荷下响应与激励的关系可表示为

      $$ \left[\begin{array}{c}{{{\varepsilon}} }_{1}\\ {{{\varepsilon}} }_{2}\\ \vdots \\ {{{\varepsilon}} }_{N}\end{array}\right]=\left[\begin{array}{c}{{{G}}}_{11}\\ {{{G}}}_{21}\\ \vdots \\ {{{G}}}_{{N}1}\end{array} \begin{array}{c}{{{G}}}_{12}\\ {{{G}}}_{22}\\ \vdots \\ {{{G}}}_{{N}2}\end{array} \begin{array}{c}\cdots \\ \cdots \\ \ddots \\ \cdots \end{array} \begin{array}{c}{{{G}}}_{1{M}}\\ {{{G}}}_{2{M}}\\ \vdots \\ {{{G}}}_{{NM}}\end{array}\right]\left[\begin{array}{c}{{{F}}}_{1}\\ {{{F}}}_{2}\\ \vdots \\ {{{F}}}_{M}\end{array}\right]\Delta t$$ (8)

      式中:$M$为待识别载荷子区域的数目;$N$为用于载荷识别的测点数目,需满足条件$N \geqslant M$${{{G}}_{ij}}$$j$区域对应$i$测点的格林函数矩阵;${{{\varepsilon}} _j}$$j$测点的结构响应时程;${{{F}}_i}$$i$区域的载荷时程。

      为方便描述,可将式(7)和式(8)统一写作:

      $${{\varepsilon}} = {{GF}}$$ (9)

      本文结构响应$\varepsilon $采用肋骨上、下测点剪切应变差值$\Delta \gamma $来描述,因此船体结构冰载荷${F_{\rm{ice}}}$识别的正问题可由下式表示:

      $$\Delta \gamma = {{GF}}{_{{\rm{ice}}}}$$ (10)

      远场冰载荷识别中的格林函数矩阵${{G}}$的条件数$cond\left( {{G}} \right)$很大,当输入数据有微小扰动或计算过程中的舍入误差都可能引起输出数据的很大扰动,导致解严重失真[24]。而在数据测量中,由于测量的系统误差和信息采集中不可避免的噪声影响,船体实测结构响应和实际响应存在一定的差异,可能导致方程无解或者有无穷多解,因此需引入正则化方法来克服冰载荷识别问题的病态性。

      正则化的基本思想是给不适定问题添加一个附加条件,求得一组可以接受的解,再从中选出最适合方程的解。本文采用共轭梯度最小二乘迭代算法(CGLS)对冰载荷识别反问题进行求解,其基本思想是应用共轭梯度法求解下式的最小值问题[25]

      $$min||\Delta {{\gamma}} - {{GF}}_{m}|{|_2}$$ (11)

      式中:${{{F}}_{{m}}}$m次迭代后的待识别冰载荷的解。

      具体的冰载荷识别流程图如图4所示。其主要思路是通过有限元分析,分别在每个子加载区域加载单位冲击载荷,求出各加载区域与每个测点响应之间的格林函数矩阵;通过简单处理实测的应变信息找到船体冰激响应数据,构建方程(10)并通过迭代正则化方法解出方程的近似稳定解,即为所识别的冰载荷时间历程。

      图  4  冰载荷识别的流程图

      Figure 4.  The flow chart of ice load identification

    • 为评估基于格林函数和正则化方法的时域反卷积识别方法对“天恩”号部分舷侧结构所受冰载荷的识别能力,建立相应的有限元模型并在模型上进行多源远场载荷识别的数值算例验证。根据“天恩”号的船艏型线图、船艏基本结构图和现场实拍照片,采用ANSYS软件建立船艏区域部分舷侧外板的有限元模型。有限元模型完全采用SHELL181单元进行建立,模型的主要参数和示意图分别如表2图5所示。骨材尺寸信息及板厚通过现场测量得到。模型四周边界条件设置为刚性固定。

      表 2  “天恩”号极地运输船的有限元模型主要参数

      Table 2.  Key parameters of the FE model of ice class vessel Tian En

      参数数值
      弹性模量E/GPa209
      泊松比λ0.3
      密度ρ/(kg·m−3)7 850

      图  5  “天恩”号测试区域的有限元模型

      Figure 5.  The FE model of monitored area of Tian En

      数值算例的测点及加载区域示意图如图6(a)所示。由于强肋骨结构与普通肋骨尺寸、厚度信息不同,4组测点从右至左分别设置在Fr239,Fr237,Fr236和Fr235处的4根普通肋骨上。图中三角形处为肋骨的测点选取位置,同一肋骨的上、下两个测点组成一个测点组,从右至左命名为第1~4个测点组,用于定义肋骨处的剪应变差值。矩形框内为设置的子加载区域,子加载区域水平宽度为船艏区域的肋骨间距,高度设置为水线上下各0.5 m。首先在图中的1号和2号子加载区域上分别施加1 Pa的单位冲击载荷,载荷时间步长设置为0.005 s,得到各测点的响应,其中1号子加载区域到各测点的格林函数响应如图6(b)所示。将各子加载区域对应各测点的响应代入式(8)求得G中的每个元素,组装后求出监测区域对应加载区域的格林函数矩阵G

      图  6  数值计算的测点、加载区域示意图以及单位冲击载荷下的格林函数响应

      Figure 6.  The sketches of loading area, gauging matrix and the response of Green's function under specific surface loading in the numerical simulation

      求出对应的格林函数矩阵后,在子加载区域分别施加已知的载荷时程,并将测点组的剪切应变响应代入式(10)解出载荷时程。在以往的冰载荷实际测量中可以看出冰载荷呈现为类似三角波的形式[10],再考虑到船−冰作用的随机性,分别在1号和2号子加载区域中施加单三角波和多三角波形式的载荷。图7所示为无噪声影响下的多源远场载荷识别算例的结果,图中实线表示施加的均布载荷时程,识别载荷的采样频率设置为50 Hz,从图中可以看出,基于格林函数和正则化方法的冰载荷识别模型能较好地反演出真实载荷时程。

      图  7  数值算例的计算结果

      Figure 7.  The results of the numerical simulation

    • 为验证格林函数和正则化方法对实船冰激响应的识别能力,基于所建立的“天恩”号部分外板有限元模型,对2个典型的船−冰相互作用工况下的船体冰激应变响应进行冰载荷识别计算,并对不同工况下的船体冰载荷特征进行对比分析。

    • 根据冰区航行过程中数字摄像机所采集的图像信息,本次冰区航行中船−冰相互作用下的海冰破坏模式主要包括挤压破碎、劈裂破坏和翻转滑移。其中挤压破碎主要出现在海冰碰撞点处的局部区域,整体的劈裂破坏或者翻转滑移取决于海冰尺寸及船−冰碰撞角度。以下选取2种不同海冰破坏模式下的冰激应变响应进行冰载荷识别分析,2种计算工况的主要参数见表3

      表 3  计算工况主要参数

      Table 3.  The main parameters of navigation conditions

      编号日期UTC
      时间
      航速
      /(m·s–1)
      吃水
      /m
      冰状冰厚
      /m
      海冰
      破坏模式
      工况12019-08-0212:025.146.5小冰盘1.2局部挤压、
      整体劈裂
      工况22019-08-033:084.246.5冰块1.5翻转滑移

      图8所示为2种不同工况下“天恩”号极地运输船与浮冰发生相互作用的视频图像。其中工况1为“天恩”号航行过程中撞上的最大浮冰块,船体左舷与一块小面积的平整冰迎面撞上。其作用过程可以描述为:最初的浮冰与船体碰撞,产生裂纹;裂纹开始扩展且碰撞点处的海冰发生局部挤压现象破碎成大量小块碎冰;同时浮冰的局部裂纹扩展至浮冰整体使之发生劈裂破坏,从而大块浮冰被分割成多块较小的浮冰;最后船−冰不再碰撞,冰载荷事件结束。工况2为一块较大浮冰块与艏部左舷区域相撞,局部的挤压破碎现象不明显,碰撞产生的能量主要转化为浮冰块的动能,使之发生翻转滑移,与船体不再发生碰撞。

      图  8  “天恩”号船体与浮冰的相互作用

      Figure 8.  The interaction between Tian En and the ice floe

      根据船−冰作用图像信息所记录的时间可从应变采集系统中找到对应的冰激应变时程。在识别计算中选取了Fr235,Fr236,Fr237,Fr239,Fr240和Fr241处的6根相同规格的普通肋骨,分析肋骨上、下两个测点之间的剪切应变差值,其时程曲线如图9所示。其中:图9(a)为工况1下的冰激剪切应变差值时程曲线,应变变化周期约为1.8 s,应变峰值约为9 μm/m;图9(b)为工况2下的冰激剪切应变差值时程曲线,应变变化周期约为0.3 s,应变峰值约为1.2 μm/m。从中可以明显发现,冰激应变周期和峰值与海冰破坏模式有很明显的关系,局部挤压状态下的海冰,碰撞后发生小块碎冰堆积现象,与船体持续碰撞,产生的应变响应较高且持续周期较长。而发生翻转滑移现象的小块浮冰,碰撞后会“弹开”,冰载荷事件下仅表现为一次短暂的撞击现象,因此产生的应变响应较低且周期较短。

      图  9  不同工况下的肋骨剪切应变差值测量结果

      Figure 9.  The measured results of frame shear strain difference under various navigation conditions

      以下选取工况1和工况2进行局部冰载荷识别计算。2个工况的测点及子加载区域的选取如图10所示,三角形处为肋骨的传感器布置位置,上、下两个应变片传感器组成一个测点组,矩形框内为子加载区域,共设置了6个测点组和6个子加载区域。6个测点组响应分别是Fr235,Fr236,Fr237,Fr239,Fr240和Fr241处的6根普通肋骨上、下两应变传感器的剪切应变差值;6个子加载区域根据视频信息中船−冰作用点的位置确定,子加载区域水平宽度为船艏区域的肋骨间距,高度设置为水线上下各0.5 m。

      图  10  不同工况的测点及子加载区域设置

      Figure 10.  The monitored area and gauge map under different conditions

    • 在确定了测点及子加载区域后,分别在6个子加载区域上施加1 Pa的单位冲击载荷,载荷时间步长设置为0.005 s。根据各子加载区域对应的测点响应信息求得对应的格林函数矩阵,将矩阵G和实测肋骨剪切应变差值代入式(10)计算,通过CGLS正则化算法解出冰载荷时间历程。工况1、工况2计算出的冰载荷时程曲线如图11所示,2个工况均可定义为一次冰载荷事件,没有明显的二次碰撞现象,其船−冰碰撞区域面积、作用周期及冰载荷峰值均有所不同。对于每个工况,各子加载区域冰载荷的周期相同而幅值不同,均可以等效为三角波的形式。从曲线光滑程度上看,工况1的冰载荷时程曲线较为平滑,信噪比较高,而工况2下的局部冰载荷数值受噪声影响较大,导致冰载荷识别在一定程度上失真。此外,将计算出的各工况下子加载区域最大冰载荷和测点的最大应变响应绘制在一幅图中,如图12所示。从中可以明显看出,船−冰相互作用下的冰载荷的作用周期与船体的应变响应周期一致,对船体局部区域应变的监测可以有效地实现对冰载荷作用周期的监测。

      图  11  不同工况下的冰载荷识别结果

      Figure 11.  The identified ice load under different conditions

      图  12  不同工况下的最大冰载荷时程曲线与测点的最大应变响应时程曲线

      Figure 12.  The time-history curves of maximum ice load and maximum strain response of measuring point under different conditions

    • 船−冰相互作用是一个动态的加载过程,根据计算出的各子区域载荷时程,可以看出冰载荷的有效作用面积往往不局限于一个肋位。工况1的有效作用肋位为3个,工况2的有效作用肋位为2个,且冰载荷随着肋位的变化表现出衰减的趋势。其原因是随着船−冰相互作用过程的进行,浮冰逐渐破碎,导致有效作用在各肋位上的浮冰尺寸总体呈现出减小的趋势。

      识别计算出的冰载荷为均布载荷的形式,实际上是冰压力,而船−冰作用下的冰层压力分布并不是均匀的,水线附近会产生一个明显的高压力区。因此在船体结构监测中,通常采用线载荷(kN/m)的形式描述冰载荷[6]。线载荷和冰压力的关系可由式(12)表示:

      $${F_{\rm{L}}} = {F_{{\rm{ice}}}} \cdot B$$ (12)

      式中:${F_{\rm{L}}}$为线载荷;${F_{{\rm{ice}}}}$为识别计算中求得的冰压力;$B$为划分子区域的单位宽度。

      通过计算转换,可将2个工况下的总冰载荷以线载荷的形式表示,如图13所示。工况1(冰厚1.2 m)对应的最大总冰载荷为323.5 kN/m,工况2(冰厚1.5 m)对应的最大总冰载荷为69.3 kN/m。可见在不同海冰破坏模式下,冰厚已不再是影响船体冰载荷特征的最主要因素,浮冰尺寸和船−冰作用面积成为影响冰载荷大小、周期的关键因素。结合海冰图像信息分析,工况1中有局部海冰挤压破碎堆积现象的发生,在这个过程中海冰与船体保持接触,持续加载,从而导致冰载荷在大小和周期上都大于工况2。此外,通过对比发现,不同工况下冰载荷时程曲线在持续增长阶段的斜率较为相近,即在加载过程中,冰载荷的增长速率不受冰厚、海冰破坏模式等海冰参数的影响,只与海冰破坏强度有关。

      图  13  不同工况下的总冰载荷时程曲线

      Figure 13.  The time history curve of global ice load under different conditions

    • 本文针对“天恩”号极地运输船在2019年7月北极东北航道的航行,安装应变传感器和图像采集装置测量了冰载荷作用下的肋骨剪切应变及海冰参数。基于实测的冰激响应数据,采用格林函数和正则化方法克服了病态方程求解的不适定问题,进行了冰载荷的识别计算。计算选取了2个不同海冰破坏模式的工况,结果表明格林函数和正则化方法能有效识别出船−冰作用下的冰载荷时程,识别效果受冰激响应信号信噪比的影响较大。对识别结果对比分析发现,冰载荷时程表现出类三角波的形式,冰载荷作用周期等同于应变响应周期。对于不同的海冰破坏模式,浮冰尺寸和船−冰接触面积成为影响冰载荷大小、周期、作用面积等特性的关键因素。加载过程中的冰载荷增长速率取决于海冰破坏强度,与海冰参数无关。今后还需要进行更多实测数据的识别分析来发现、完善冰载荷作用规律,并建立海冰参数与冰载荷之间的对应关系,为冰区船舶的结构设计与安全航行提供参考依据。

参考文献 (25)

目录

    /

    返回文章
    返回