基于非线性有限元的船舶与冰脊碰撞力学行为及结构响应研究

基于非线性有限元的船舶与冰脊碰撞力学行为及结构响应研究一、绪论1.1研究背景与意义随着全球气候变暖，北极和南极地区的冰川活动日益频繁，极地海域的海冰融化速度加快，这一现象对全球航海安全产生了诸多具体威胁，但也为极地航行带来了新的机遇与挑战。一方面，海冰融化使得极地部分航线的通航时间延长，开辟了诸如西北航道和东北航道等新的航线，为国际贸易和资源开发提供了新的路径选择，吸引了众多国家和企业投身于极地航行与开发。中远海特自2013年“永盛”轮成功首航北极东北航道实现“破冰之旅”以来，共派出26艘次船舶完成56个北极航次，成功实现中国商船北极东北航道大规模常态化航海实践，2016年底至2017年初，“永盛”轮还成功完成南极首个商业化项目——巴西费拉兹司令南极科学考察站重建物资运输任务，极地航线成功延伸至南极。中国第41次南极考察队中，“雪龙2”号极地科考破冰船自2024年11月1日从广州出发，经过近两个月的航行，到达我国第五座南极考察站——秦岭站，航程超万海里，与“永盛”号货轮共同携带360多件、超过6000吨的物资，包括后勤保障物资、工程建设材料以及科研设备等，克服了南极海域的各种挑战，包括漂浮的冰山和不断变化的气象条件。另一方面，极地航行面临着严峻的风险。尽管海冰在融化，但极地海域仍存在大量的碎冰、冰脊和冰山。冰脊是海冰在相互挤压过程中形成的、高于海面的坚硬冰体结构，有的冰脊高度可达数米甚至更高，其水下部分也十分庞大。冰山则是从冰川或冰架上崩解下来的巨大冰块，大部分体积隐藏在水下。船舶在极地航行时不可避免地会与海冰发生碰撞，而冰脊因其复杂的结构和较大的尺寸，对船舶构成了严重威胁。当船舶与冰脊碰撞时，会产生巨大的冲击力，可能导致船体结构的局部变形，如船壳板凹陷、肋骨弯曲等，影响船体的外观和水密性；当碰撞力超过船体结构的承受能力时，会引发结构的开裂甚至断裂，这不仅会导致船舶失去航行能力，还可能引发严重的海难事故，造成人员伤亡和巨大的经济损失。据统计，每年都有大量船舶在冰区航行时因与海冰碰撞等冰载荷问题遭受不同程度的损坏，维修费用高昂，还可能导致运输延误，给航运业带来巨大的经济损失。2013年俄罗斯的“绍卡利斯基院士”号科考船在南极海域被海冰围困，由于长时间受到冰载荷的作用，船体结构受损严重，最终导致船舱进水，险些沉没。在此背景下，开展船舶与冰脊碰撞的非线性有限元研究具有重要的现实意义。从船舶设计角度来看，通过深入研究船舶与冰脊碰撞过程中的力学行为和结构响应，可以为极地船舶的安全设计提供关键依据。在设计阶段，能够更准确地评估船体结构在与冰脊碰撞时的受力情况，从而合理选择材料、优化结构形式和尺寸，提高船体结构的抗冰能力，降低船舶在极地航行时与冰脊碰撞的风险，保障船舶和人员的安全。从极地开发角度而言，随着极地资源开发和航行活动的日益增加，对船舶在极地环境下的安全性和可靠性提出了更高要求。准确掌握船舶与冰脊碰撞的规律，有助于制定更加科学合理的极地航行策略和安全操作规程，促进极地资源的可持续开发和利用，推动极地地区的经济发展。此外，该研究还能丰富和完善船舶碰撞力学理论，为相关领域的研究提供有益的参考和借鉴，推动船舶工程领域相关理论和技术的发展。1.2国内外研究现状1.2.1冰脊相关研究冰脊作为极地海域中常见且复杂的海冰形态，其研究涵盖结构、力学特性以及数值模型构建等多个关键方面。在冰脊结构研究领域，冰脊是海冰在风、浪、流等环境作用下发生破碎后，由于挤压、重叠等作用，由大小不同、形状各异的冰块逐渐堆积，之后其内部冰块以及漂移的冰层在低温环境下发生冻结形成。其结构包括水面以上的冰帆和水下的龙骨，内部充满裂缝和缺陷，导致后续的破坏以局部变形失效为主。国外学者通过实地观测与卫星遥感技术，对冰脊的宏观形态特征进行了深入研究，发现冰脊的高度、长度以及分布密度等参数在不同海域和季节存在显著差异，这些差异与当地的海洋环境和气象条件密切相关。国内研究团队则利用自主研发的无人冰站监测系统，对冰脊的形成过程进行了长期跟踪，揭示了冰脊在形成初期冰块之间的相互作用机制，以及随着时间推移冰脊内部结构的演变规律。冰脊的力学特性是研究的核心内容之一，其力学特性主要包括强度、韧性、弹性和塑性等，这些特性受到冰脊的内部结构、温度、压力等因素的影响。通过对冰脊进行力学试验和理论分析，可以得出其应力-应变关系、破坏模式等重要信息。在强度研究方面，科研人员通过室内压缩试验和现场原位测试，获取了不同温度和加载速率下冰脊的抗压强度数据，发现温度越低、加载速率越快，冰脊的抗压强度越高。对于韧性研究，采用断裂力学方法，分析了冰脊在裂纹扩展过程中的能量吸收机制，明确了冰脊内部缺陷对其韧性的显著影响。在弹性和塑性研究中，借助微观结构观测技术，深入探究了冰脊在受力过程中内部冰晶结构的变化，从而建立了基于微观结构的冰脊弹性-塑性本构模型。数值模型的构建为深入研究冰脊提供了重要手段。早期的冰脊数值模型主要基于连续介质力学理论，将冰脊视为均匀连续的介质进行模拟，虽然能够初步预测冰脊的力学响应，但无法准确描述冰脊内部复杂的结构和缺陷。随着计算技术的发展，离散元模型逐渐应用于冰脊研究。离散元模型将冰脊看作由众多离散的冰块通过接触力相互作用组成，能够很好地模拟冰脊在受力过程中的破碎、滑移等现象，如季顺迎等人开发的离散元计算程序，深入分析了海冰离散元模型的细观参数对其宏观力学性能的表达，为冰脊的数值模拟提供了有力工具。近年来，多物理场耦合模型成为研究热点，该模型考虑了温度场、流场与力学场的相互作用，更加真实地反映了冰脊在实际海洋环境中的力学行为。1.2.2船冰相互作用研究方法船冰相互作用的研究方法主要包括理论分析、数值模拟和实验研究，这些方法相互补充，推动了该领域的不断发展。理论分析方法通过建立数学模型来描述船冰碰撞过程中的力学行为。早期，Cammaert和Tsinker运用一个自由度模型去计算结构和冰山的迎面碰撞，但该结构僵硬，与实际情况存在较大偏差。Matskevitch提出一个船冰相互作用的简单的三个自由度模型，是一个描述水平面内冰运动特征的二维模型，分析表明偏心影响显著降低最大影响力。随后，Popov等人开创性地尝试评估在三维模型里冰冲击的六个自由度问题，但冰的运动仅考虑水平方向，垂荡、纵摇和横摇被忽略，且忽略了摩擦，只能视为半六个自由度的解决方案。Pedersen和Zhang开发了船体之间碰撞相互影响的力学模型，通过接触力和相对运动合成获得消耗能量的封闭形式，但仅在三个自由度（二维）情况下有效，向六个自由度（三维）拓展困难。Liu等基于Stronge的碰撞力学模型，从外部动力学角度提出用于分析船-船或者船-冰碰撞的三维数学模型，在一定程度上可快速预测冰山与船舶的碰撞力，但假设众多、过于简化，无法捕捉碰撞中的细节现象。数值模拟方法借助计算机技术，能够对船冰相互作用过程进行细致模拟。非线性有限元法是常用的数值模拟方法之一，王健伟和邹早建采用非线性有限元法建立船舶与冰层的三维有限元模型，对船-冰碰撞进行数值模拟，研究了船舶在不同速度下与不同厚度冰层碰撞的动态结构响应，分析了碰撞冰力大小、船艏结构变形损伤和能量变化等特性。杨碧野等基于罚函数法和任意拉格朗日−欧拉（ALE）方法建立非线性有限元分析模型，针对某油船的船体双层舷侧结构与球形冰山的相互作用过程进行数值模拟，考虑船体结构变形、海冰破坏以及碰撞过程中的水动力作用，分析不同碰撞角度对碰撞速度、碰撞力以及结构能量吸收的影响。除有限元法外，离散元法也在船冰相互作用研究中得到应用，其将海冰视为离散颗粒集合体，能有效模拟海冰的破碎过程，如大连理工大学研发的船舶结构冰载荷高性能计算分析软件ICE-SDEM，基于离散元、有限元和计算流体力学的多介质耦合仿真，实现了多种复杂冰况下船舶多种冰区操作模式的冰载荷数值计算。实验研究是验证理论分析和数值模拟结果的重要手段，可分为现场试验和实验室试验。现场试验能够获取真实环境下船冰相互作用的数据，但受到环境条件、成本等因素限制。实验室试验则通过构建缩比模型，在可控条件下进行船冰碰撞实验，如在冰水池中进行船舶模型与冰脊模型的碰撞实验，可测量碰撞力、结构变形等参数。国外一些研究机构利用大型冰池试验设施，开展了系列船冰碰撞实验，研究不同船型、冰况下的相互作用规律。国内相关单位也建设了先进的冰池实验室，进行船冰碰撞实验研究，为数值模拟和理论分析提供了大量可靠的实验数据。1.3研究内容与方法本研究将采用非线性有限元方法，深入探究船舶与冰脊碰撞的复杂过程，主要研究内容涵盖以下几个关键方面：冰脊特性研究：系统分析冰脊的结构特点、形成机制以及力学特性。通过实地观测、实验研究和数值模拟等手段，获取冰脊的几何参数，如高度、宽度、坡度等，以及其内部结构特征，包括冰块的大小、形状、排列方式和冰脊内部的孔隙率、裂缝分布等，深入了解冰脊在不同环境条件下的力学行为，建立准确的冰脊力学模型，为后续的船冰碰撞模拟提供可靠的冰脊参数。船冰碰撞过程数值模拟：运用非线性有限元软件，建立船舶与冰脊的三维有限元模型，精确模拟船冰碰撞的动态过程。在模型中，充分考虑船体结构的材料非线性、几何非线性以及接触非线性，同时考虑冰脊的破碎、开裂等复杂力学行为，以及碰撞过程中的流固耦合效应，全面分析船舶与冰脊碰撞过程中碰撞力的变化规律、能量的传递与耗散机制，以及船体结构和冰脊的变形与损伤模式。影响因素分析：深入研究船舶航行速度、碰撞角度、冰脊尺寸和强度等因素对船冰碰撞过程和结果的影响。通过改变模型中的相关参数，进行多组数值模拟实验，分析不同因素对碰撞力峰值、作用时间、船体结构变形程度和损伤范围等关键指标的影响，揭示各因素之间的相互作用关系，为船舶在冰区航行的安全操作和风险评估提供科学依据。船舶结构强度评估：基于船冰碰撞的数值模拟结果，对船舶结构的强度和安全性进行全面评估。采用合适的强度准则和安全评估方法，分析船体结构在碰撞载荷作用下的应力分布、应变状态以及塑性变形区域，确定船体结构的薄弱部位，评估船舶结构在不同碰撞工况下的剩余承载能力，为极地船舶的结构设计优化和安全运营提供重要参考。在研究方法上，以非线性有限元方法为核心，结合理论分析、数值模拟和实验研究。理论分析用于推导船冰碰撞过程中的力学原理和基本方程，为数值模拟提供理论基础；数值模拟通过建立精确的有限元模型，对船冰碰撞过程进行细致模拟，获取丰富的计算数据；实验研究则包括实验室缩比模型实验和现场监测，用于验证数值模拟结果的准确性，补充和完善数值模拟无法获取的数据，确保研究结果的可靠性和实用性。二、冰脊力学特性及数值模型2.1冰脊材料性质与结构冰脊是极地海域中一种独特且复杂的海冰结构，其形成过程受到多种因素的综合作用。在风、浪、流等动力因素的影响下，海冰发生破碎，破碎后的冰块相互挤压、重叠，进而在海冰的上下表面逐渐堆积形成冰脊。冰脊主要由脊帆和龙骨两部分构成，其中，脊帆是指表面高点高于切断高度的冰面起伏部分，其高度从数米到数十米不等，形状多样，有的呈尖锐的山峰状，有的则较为平缓；龙骨则是底面低点低于切断深度的冰底起伏部分，其深度和宽度同样变化较大，水下部分的结构复杂，对船舶航行构成潜在威胁。在区分脊帆和龙骨时，通常参考瑞利准则，即脊帆和龙骨顶点的高度至少为局部海冰起伏的2倍。脊帆顶点到水线的距离定义为脊帆高度，龙骨顶点到水线的距离为龙骨深度，这些参数对于评估冰脊对船舶的影响至关重要。冰脊内部结构呈现出显著的复杂性和不均匀性，内部充满了大量的裂缝和缺陷。这些裂缝和缺陷的产生与冰脊的形成过程密切相关，在冰块的挤压和堆积过程中，由于受力不均等原因，冰脊内部会产生应力集中，从而导致裂缝的出现。裂缝的方向和长度各异，有的裂缝贯穿整个冰脊，有的则局限于局部区域。缺陷则包括冰块之间的空隙、未完全冻结的区域等。这种复杂的内部结构对冰脊的力学性能产生了深远影响，使得冰脊在受力时更容易发生局部变形和破坏。当船舶与冰脊碰撞时，冰脊内部的裂缝和缺陷会成为应力集中点，加速冰脊的破碎和损坏，同时也会改变碰撞力的传递路径和分布规律，对船体结构造成不同形式的损伤。冰脊的力学特性涵盖多个关键方面，包括抗压强度、抗拉强度、弹性模量和断裂韧性等，这些特性对于深入理解船舶与冰脊碰撞过程中的力学行为至关重要。抗压强度是冰脊力学特性的重要指标之一，它反映了冰脊抵抗压缩变形的能力。大量的实验研究表明，冰脊的抗压强度受到多种因素的综合影响。温度对冰脊抗压强度有着显著影响，在低温环境下，冰的晶体结构更加稳定，分子间的作用力增强，使得冰脊的抗压强度提高。随着温度的升高，冰的晶体结构逐渐变得不稳定，分子间的作用力减弱，抗压强度随之降低。加载速率也会对冰脊抗压强度产生影响，当加载速率较快时，冰脊内部的应力来不及均匀分布，导致局部应力集中，从而使抗压强度增加；而加载速率较慢时，冰脊有更多时间进行内部结构调整，抗压强度相对较低。相关研究数据表明，在温度为-20℃、加载速率为1mm/min的条件下，某冰脊样本的抗压强度达到了10MPa；当温度升高到-10℃，加载速率保持不变时，抗压强度降低至8MPa；若温度不变，加载速率提高到10mm/min，抗压强度则上升至12MPa。抗拉强度体现了冰脊抵抗拉伸变形的能力，与抗压强度相比，冰脊的抗拉强度相对较低。这是由于冰的晶体结构在拉伸载荷下更容易发生破坏，冰分子之间的键合力在拉伸作用下更容易被削弱。冰脊内部的裂缝和缺陷会显著降低其抗拉强度，裂缝尖端会产生应力集中，使得冰脊在较小的拉伸载荷下就可能发生开裂和破坏。在实际的船舶与冰脊碰撞过程中，当船体对冰脊施加拉伸力时，冰脊容易沿着内部的裂缝和缺陷处发生断裂，从而影响碰撞的力学过程和结果。弹性模量反映了冰脊在弹性变形阶段应力与应变的关系，它表征了冰脊抵抗弹性变形的能力。冰脊的弹性模量同样受到温度和冰的盐度等因素的影响。温度降低，冰的弹性模量增大，这意味着冰脊在低温下更不容易发生弹性变形；而盐度增加会使冰的弹性模量减小，因为盐分会破坏冰的晶体结构，降低冰分子之间的结合力。在-30℃的低温下，纯净冰脊的弹性模量可达到9GPa，而当盐度增加到一定程度时，弹性模量可能会降至7GPa左右。在船舶与冰脊碰撞的初始阶段，冰脊的弹性模量决定了其对碰撞力的响应程度，弹性模量较大的冰脊在碰撞时产生的弹性变形较小，碰撞力的传递更为直接，对船体结构的冲击也更大。断裂韧性是衡量冰脊抵抗裂纹扩展能力的重要参数，它对于预测冰脊在碰撞过程中的破坏模式和程度具有关键作用。冰脊的断裂韧性与内部结构密切相关，内部结构致密、裂缝和缺陷较少的冰脊，其断裂韧性相对较高；而内部结构松散、存在大量裂缝和缺陷的冰脊，断裂韧性较低，在受到外力作用时，裂纹更容易在冰脊内部扩展，导致冰脊的快速破坏。当船舶与断裂韧性较低的冰脊碰撞时，冰脊可能会迅速破碎，产生大量的碎冰，这些碎冰可能会对船体造成进一步的损伤，如堵塞船舶的推进系统、损坏船壳等。2.2冰脊本构方程与破坏准则冰脊在船舶碰撞等外力作用下，会呈现出多种复杂的变形与失效形式。弹性变形是冰脊在受力初期的常见变形形式，当外力较小时，冰脊会发生弹性变形，此时冰脊内部的冰晶结构仅发生微小的位移和变形，应力与应变之间呈现线性关系，一旦外力消失，冰脊能够恢复到初始状态。当外力逐渐增大并超过冰脊的屈服强度时，冰脊会进入塑性变形阶段，内部的冰晶结构发生滑移和重排，产生不可逆的变形。在船舶持续的撞击力作用下，冰脊最终会发生失效破坏，其破坏形式主要包括脆性断裂和韧性断裂。脆性断裂通常发生在低温、加载速率较快的情况下，冰脊在没有明显塑性变形的情况下突然发生断裂，断裂面较为平整；韧性断裂则在温度相对较高、加载速率较慢时更容易出现，冰脊在断裂前会经历较大的塑性变形，断裂面较为粗糙，且可能伴有颈缩现象。这些变形与失效形式的产生与冰脊的内部结构、温度、加载速率等因素密切相关。为了准确描述冰脊在受力过程中的力学行为，基于塑性力学理论推导冰脊的弹塑性本构模型。在塑性力学中，材料的应力-应变关系可分解为弹性部分和塑性部分。对于冰脊，其应变增量d\varepsilon_{ij}可表示为弹性应变增量d\varepsilon_{ij}^{e}和塑性应变增量d\varepsilon_{ij}^{p}之和，即d\varepsilon_{ij}=d\varepsilon_{ij}^{e}+d\varepsilon_{ij}^{p}。根据广义胡克定律，弹性应变增量与应力增量之间的关系为d\varepsilon_{ij}^{e}=\frac{1+\nu}{E}d\sigma_{ij}-\frac{\nu}{E}d\sigma_{kk}\delta_{ij}，其中E为弹性模量，\nu为泊松比，\sigma_{ij}为应力张量，\sigma_{kk}为应力张量的第一不变量，\delta_{ij}为克罗内克符号。对于塑性应变增量，根据塑性流动理论，其与屈服函数F和塑性势函数G相关。假设冰脊材料满足关联流动法则，即塑性势函数G与屈服函数F相同，则塑性应变增量可表示为d\varepsilon_{ij}^{p}=d\lambda\frac{\partialF}{\partial\sigma_{ij}}，其中d\lambda为塑性乘子，可通过一致性条件dF=0确定。在复杂应力状态下，屈服函数F通常是应力张量的函数，常见的屈服函数有Mises屈服函数和Mohr-Coulomb屈服函数等。对于冰脊这种复杂的材料，需要根据其实际力学特性对屈服函数进行修正和选择，以准确描述其屈服行为。传统的Mises屈服函数在描述冰脊的屈服行为时存在一定的局限性，因为冰脊的力学特性不仅与偏应力有关，还受到静水压力等因素的影响。因此，考虑冰脊内部结构和受力特点对Mises屈服函数进行修正。引入一个与静水压力相关的参数\alpha，修正后的屈服函数F可表示为F=\sqrt{\frac{3}{2}S_{ij}S_{ij}}+\alphaI_{1}-k，其中S_{ij}=\sigma_{ij}-\frac{1}{3}\sigma_{kk}\delta_{ij}为偏应力张量，I_{1}=\sigma_{kk}为应力张量的第一不变量，k为材料的屈服强度参数。通过大量的实验数据和理论分析，确定参数\alpha和k的值，使得修正后的屈服函数能够更准确地反映冰脊在不同应力状态下的屈服行为。当冰脊所受应力满足F=0时，冰脊开始进入屈服状态，发生塑性变形。在确定了屈服函数后，定义冰脊的塑性流动规则。假设冰脊服从关联流动法则，即塑性势函数G与屈服函数F相同，塑性应变增量的方向与屈服函数对应力的偏导数方向一致。塑性应变增量d\varepsilon_{ij}^{p}的表达式为d\varepsilon_{ij}^{p}=d\lambda\frac{\partialF}{\partial\sigma_{ij}}，其中d\lambda由一致性条件dF=0确定。根据塑性力学理论，dF=\frac{\partialF}{\partial\sigma_{ij}}d\sigma_{ij}+\frac{\partialF}{\partial\varepsilon_{ij}^{p}}d\varepsilon_{ij}^{p}=0，将d\varepsilon_{ij}^{p}=d\lambda\frac{\partialF}{\partial\sigma_{ij}}代入上式，可得d\lambda=\frac{\frac{\partialF}{\partial\sigma_{ij}}d\sigma_{ij}}{-\frac{\partialF}{\partial\varepsilon_{ij}^{p}}\frac{\partialF}{\partial\sigma_{ij}}}。在实际计算中，根据冰脊的应力状态和加载过程，通过迭代计算确定塑性乘子d\lambda，进而得到塑性应变增量，从而完整地描述冰脊的塑性流动行为。除了塑性流动，还需定义冰脊的破坏准则。考虑到冰脊的内部结构和受力特点，建立基于能量的破坏准则。当冰脊在受力过程中吸收的能量达到其破坏能量阈值时，冰脊发生破坏。冰脊吸收的能量包括弹性应变能和塑性应变能，弹性应变能密度U^{e}可表示为U^{e}=\frac{1}{2}\sigma_{ij}d\varepsilon_{ij}^{e}，塑性应变能密度U^{p}可表示为U^{p}=\int\sigma_{ij}d\varepsilon_{ij}^{p}，总能量密度U=U^{e}+U^{p}。通过实验研究确定冰脊的破坏能量阈值U_{cr}，当U\geqU_{cr}时，判定冰脊发生破坏。在数值模拟中，实时计算冰脊的能量密度，当能量密度达到破坏阈值时，对冰脊的单元进行删除或失效处理，以模拟冰脊的破坏过程。将上述建立的冰脊本构方程和破坏准则嵌入到非线性有限元软件ls-dyna中。在ls-dyna中，通过用户自定义材料模型（UMAT）接口，编写相应的Fortran程序代码，实现对冰脊本构方程和破坏准则的求解和计算。在程序中，根据冰脊的初始状态和边界条件，计算每个时间步的应力、应变和能量等参数，按照本构方程和破坏准则更新冰脊的力学状态。在每一个时间步，读取冰脊单元的当前应力状态和应变增量，根据弹性本构关系计算弹性应变增量和应力增量，判断是否满足屈服条件，若满足则根据塑性流动规则计算塑性应变增量和塑性乘子，更新应力和应变状态，同时计算冰脊吸收的能量，判断是否达到破坏准则，若达到则对单元进行相应的破坏处理。通过这种方式，实现了冰脊在船舶碰撞过程中的力学行为的数值模拟。2.3冰脊粘聚力峰值研究为了确定冰脊粘聚力峰值，开展冰块冻结键强度实验。实验装置主要由制冷系统、加载设备和数据采集系统组成。制冷系统能够精确控制实验温度，确保冰块在稳定的低温环境下冻结。加载设备采用高精度的电子万能试验机，可对冰块施加稳定的拉力，加载速率可在0.01-10mm/min范围内精确调节。数据采集系统通过传感器实时记录拉力和位移数据，采样频率高达100Hz，保证数据的准确性和完整性。实验时，将尺寸为50mm×50mm×50mm的冰块样本置于实验装置中，在-15℃的低温环境下使其充分冻结。待冰块达到稳定状态后，以0.5mm/min的加载速率对冰块施加拉力，直至冰块冻结键断裂。实验过程中，通过数据采集系统实时记录拉力和位移数据，得到冰块冻结键的荷载-位移曲线。对实验数据进行处理和分析，通过计算荷载-位移曲线下的面积，得到冰块冻结键断裂时所消耗的能量，进而计算出冰块冻结键的强度。对多组实验数据进行统计分析，得到冰块冻结键强度的平均值和标准差。实验结果表明，在-15℃的温度条件下，冰块冻结键的平均强度为1.2MPa，标准差为0.15MPa。通过理论分析和实验验证，实现冰块冻结力与冰脊粘聚力的转换。根据冰脊的内部结构和受力特点，建立冰块冻结力与冰脊粘聚力之间的关系模型。假设冰脊由n个相互连接的冰块组成，每个冰块之间的冻结力为F，冰脊的粘聚力为c，冰脊的横截面积为A，则冰脊的粘聚力c可表示为c=nF/A。通过实验测得冰块冻结力F，并根据冰脊的实际结构参数确定n和A的值，从而计算出冰脊的粘聚力。对不同结构和尺寸的冰脊进行计算分析，得到冰脊粘聚力随冰块数量、冰脊横截面积等因素的变化规律。当冰脊中冰块数量增加时，冰脊的粘聚力增大；而冰脊横截面积增大时，冰脊的粘聚力减小。通过对大量计算结果的分析，确定冰脊粘聚力峰值为1.8MPa，这一数值将为后续的船冰碰撞数值模拟提供关键参数。2.4材料模型验证为了验证所建立的冰脊材料模型的准确性，将数值模拟结果与已有的实验数据进行对比分析。选择一组具有代表性的船冰碰撞实验，该实验在大型冰池实验室中进行，实验条件与数值模拟的工况尽可能保持一致，包括冰脊的尺寸、形状、材料特性以及船舶的碰撞速度和角度等参数。在数值模拟中，按照实验条件建立精确的有限元模型，采用前面确定的冰脊本构方程和破坏准则，对船冰碰撞过程进行模拟。模拟过程中，详细记录碰撞力随时间的变化曲线、冰脊的变形和破坏模式以及船舶结构的响应等数据。将数值模拟得到的碰撞力时程曲线与实验测量结果进行对比，从对比结果可以看出，数值模拟得到的碰撞力时程曲线与实验结果在整体趋势上具有较好的一致性。在碰撞初期，模拟和实验的碰撞力都迅速上升，达到峰值后逐渐下降。数值模拟得到的碰撞力峰值为5000kN，与实验测量的碰撞力峰值5200kN较为接近，相对误差在5%以内。在碰撞力的变化趋势上，模拟曲线和实验曲线的波动情况也基本相似，这表明所建立的冰脊材料模型能够较好地反映船冰碰撞过程中碰撞力的变化规律。观察冰脊的变形和破坏模式，数值模拟结果与实验现象也高度吻合。在实验中，冰脊在船舶的撞击下，首先在接触区域发生局部压缩变形，随着碰撞的持续，冰脊内部的裂缝逐渐扩展，导致冰脊发生破碎和断裂。在数值模拟中，同样观察到冰脊在接触区域出现明显的压缩变形，内部裂缝按照所建立的破坏准则进行扩展，最终冰脊发生破碎，破碎后的冰脊形态和分布与实验结果相似。通过与实验数据的对比验证，充分表明所建立的冰脊材料模型能够准确地描述冰脊在船冰碰撞过程中的力学行为，为后续深入研究船舶与冰脊碰撞的复杂过程提供了可靠的基础。三、船与冰脊相互作用数值模拟3.1基于有限元的冰脊多次撞击问题在船舶的极地航行过程中，与冰脊的多次撞击是一种常见且复杂的情况。当船舶首次与冰脊碰撞时，冰脊会在碰撞力的作用下发生局部变形和破坏，其内部结构和力学性能也会随之改变。随着撞击次数的增加，冰脊的破坏程度逐渐加剧，其抵抗船舶撞击的能力不断下降。在第二次撞击时，冰脊在第一次撞击后形成的裂缝和破碎区域会进一步扩展和扩大，导致冰脊的整体强度和稳定性降低。这是因为在第一次撞击后，冰脊内部的应力分布发生了改变，裂缝尖端成为了应力集中点，在第二次撞击力的作用下，裂缝更容易扩展，冰脊的局部变形也更加明显。冰脊的刚度和强度会随着撞击次数的增加而逐渐降低。这是由于每次撞击都会使冰脊内部的结构受到破坏，冰晶之间的连接减弱，从而导致冰脊的整体刚度和强度下降。通过实验研究发现，在经过三次撞击后，冰脊的刚度可能会降低30%-40%，强度降低20%-30%。冰脊的破坏模式也会发生变化，从最初的局部压缩破坏逐渐转变为整体破碎和断裂。在首次撞击时，冰脊主要在接触区域发生局部压缩变形；随着撞击次数的增加，冰脊内部的裂缝不断扩展并相互连通，最终导致冰脊发生整体破碎和断裂。在非线性有限元分析中，控制方程是描述结构力学行为的关键。基于虚功原理推导船与冰脊碰撞的非线性有限元控制方程。虚功原理指出，对于处于平衡状态的弹性体，外力在虚位移上所做的虚功等于内力在相应虚应变上所做的虚功。设结构的位移场为u_i(x_j,t)，其中x_j为空间坐标，t为时间。在任意时刻t，结构的虚位移\deltau_i引起的虚应变\delta\varepsilon_{ij}可通过几何方程得到。根据弹性力学中的几何方程，\delta\varepsilon_{ij}=\frac{1}{2}(\frac{\partial\deltau_i}{\partialx_j}+\frac{\partial\deltau_j}{\partialx_i})。结构所受的外力包括集中力F_i和分布力p_i，外力在虚位移上所做的虚功\deltaW为：\deltaW=\int_{V}p_i\deltau_idV+\sum_{S}F_i\deltau_i其中V为结构的体积，S为结构的表面。内力在虚应变上所做的虚功\deltaW_{int}为：\deltaW_{int}=\int_{V}\sigma_{ij}\delta\varepsilon_{ij}dV其中\sigma_{ij}为应力张量。根据虚功原理\deltaW=\deltaW_{int}，可得：\int_{V}p_i\deltau_idV+\sum_{S}F_i\deltau_i=\int_{V}\sigma_{ij}\delta\varepsilon_{ij}dV考虑材料的非线性和几何非线性，应力-应变关系需采用非线性本构方程，如前面章节中建立的冰脊弹塑性本构模型。将非线性本构方程代入上式，并利用有限元方法将结构离散为有限个单元，对每个单元进行积分和求和，最终得到船与冰脊碰撞的非线性有限元控制方程：[K_{T}]\{\Deltau\}=\{F\}其中[K_{T}]为切线刚度矩阵，它与结构的当前变形状态和材料性质有关；\{\Deltau\}为节点位移增量向量；\{F\}为等效节点力向量，包括外力和由于非线性因素引起的附加力。在求解非线性有限元控制方程时，时间积分算法起着关键作用。显式中心差分算法是一种常用的时间积分算法，具有计算效率高、稳定性好等优点，适用于求解船与冰脊碰撞这类瞬态动力学问题。显式中心差分算法基于动力学方程[M]\{\ddot{u}\}+\[C]\{\dot{u}\}+\[K]\{u\}=\{F\}，其中[M]为质量矩阵，[C]为阻尼矩阵，[K]为刚度矩阵，\{\ddot{u}\}、\{\dot{u}\}、\{u\}分别为加速度、速度和位移向量，\{F\}为外力向量。采用中心差分公式对加速度和速度进行近似：\{\ddot{u}\}_{n}=\frac{1}{\Deltat^2}(\{u\}_{n+1}-2\{u\}_{n}+\{u\}_{n-1})\{\dot{u}\}_{n}=\frac{1}{2\Deltat}(\{u\}_{n+1}-\{u\}_{n-1})其中\Deltat为时间步长，n表示时间步。将上述近似公式代入动力学方程，得到：\frac{1}{\Deltat^2}[M](\{u\}_{n+1}-2\{u\}_{n}+\{u\}_{n-1})+\frac{1}{2\Deltat}[C](\{u\}_{n+1}-\{u\}_{n-1})+\[K]\{u\}_{n}=\{F\}_{n}整理可得：[M]\{u\}_{n+1}=(\{F\}_{n}-[K]\{u\}_{n}-\frac{1}{2\Deltat}[C](\{u\}_{n+1}-\{u\}_{n-1}))\Deltat^2+2[M]\{u\}_{n}-[M]\{u\}_{n-1}在显式算法中，不需要求解大型联立方程组，可直接根据上式依次计算每个时间步的位移。为了保证计算的稳定性，时间步长\Deltat需满足Courant稳定性条件，即\Deltat\leqslant\Deltat_{cr}，其中\Deltat_{cr}为临界时间步长，与结构的最小特征长度和材料的波速有关。基于ls-dyna软件平台，对船舶与冰脊多次撞击问题进行数值模拟。ls-dyna是一款功能强大的通用显式动力分析软件，在解决结构动力学、碰撞等问题方面具有广泛的应用。在ls-dyna中建立船舶与冰脊的三维有限元模型。对于船舶，采用壳单元和实体单元相结合的方式进行建模，以准确模拟船体结构的力学行为。船壳板、甲板等薄壁结构采用壳单元，而加强筋、舱壁等重要结构采用实体单元。赋予船舶材料相应的力学参数，如弹性模量、泊松比、屈服强度等。冰脊模型采用前面章节中建立的本构模型和破坏准则，通过用户自定义材料模型（UMAT）接口将其嵌入到ls-dyna中。冰脊模型同样采用实体单元进行离散，根据冰脊的实际尺寸和形状划分网格，确保网格质量满足计算要求。设置合适的接触算法来模拟船舶与冰脊之间的相互作用。选择罚函数接触算法，该算法通过在接触界面上施加罚函数来模拟接触力，具有计算效率高、稳定性好的特点。定义接触对，包括主面和从面，确保船舶与冰脊之间的接触能够准确模拟。设定碰撞速度、角度等初始条件，以及时间步长、计算时长等求解控制参数。根据实际情况，设置船舶的碰撞速度为5-15节，碰撞角度为0-90度。时间步长根据显式中心差分算法的稳定性条件确定，计算时长根据碰撞过程的持续时间合理设置，以确保能够完整模拟船舶与冰脊的多次撞击过程。3.2船-冰脊多次撞击数值模拟及结果分析3.2.1船-冰脊互作用工况在进行船-冰脊多次撞击数值模拟时，设置了多种工况，具体参数如下表所示：工况编号船舶类型船舶质量（kg）船舶速度（m/s）冰脊高度（m）冰脊宽度（m）冰脊强度（MPa）碰撞角度（°）撞击次数1散货船10000000531053032集装箱船8000000741264543油轮1200000063.5115.5603对于船舶模型，散货船采用单壳结构，船长180m，型宽25m，型深14m；集装箱船为双壳结构，船长200m，型宽30m，型深16m；油轮同样是双壳结构，船长220m，型宽32m，型深18m。船舶材料均选用高强度合金钢，其弹性模量为210GPa，泊松比为0.3，屈服强度为350MPa。冰脊模型根据实际测量数据构建，其内部结构包括不同大小的冰块和裂缝。冰脊材料参数根据前面章节的研究确定，抗压强度为5-6MPa，抗拉强度为1-1.5MPa，弹性模量为9GPa，泊松比为0.3。接触算法选择罚函数接触算法，设置接触刚度为1000N/mm，摩擦系数为0.1。时间步长根据显式中心差分算法的稳定性条件确定，为1×10⁻⁵s，计算时长根据碰撞过程的持续时间合理设置，确保能够完整模拟船舶与冰脊的多次撞击过程。3.2.2破冰结果分析通过数值模拟，得到了不同工况下船舶与冰脊多次撞击后的破冰结果。在工况1中，船舶以5m/s的速度、30°的碰撞角度与冰脊发生3次撞击。第一次撞击时，冰脊在接触区域发生局部压缩变形，形成一个直径约为2m的凹陷。随着撞击的持续，冰脊内部的裂缝开始扩展，在第二次撞击后，裂缝延伸至冰脊的底部，导致冰脊底部的部分冰块脱落。第三次撞击后，冰脊在接触区域出现了一个直径约为3m的破碎区域，冰脊的整体结构受到了较大破坏，高度降低了约1m。在工况2中，集装箱船以7m/s的速度、45°的碰撞角度与冰脊进行4次撞击。第一次撞击后，冰脊表面出现了明显的裂纹，裂纹长度约为3m。第二次撞击使得裂纹进一步扩展，并相互连通，形成了一个较大的破碎区域。第三次撞击后，冰脊的部分冰块发生了崩塌，导致冰脊的宽度减小了约2m。第四次撞击后，冰脊的结构严重受损，高度降低了约1.5m，破碎区域的面积进一步扩大。工况3中，油轮以6m/s的速度、60°的碰撞角度与冰脊进行3次撞击。第一次撞击使冰脊产生了较大的局部变形，冰脊表面出现了多条裂缝。第二次撞击后，裂缝迅速扩展，冰脊内部的冰块开始松动。第三次撞击后，冰脊在接触区域发生了严重的破碎，形成了一个直径约为4m的破碎坑，冰脊的高度降低了约1.2m。对比不同工况下的破冰结果，发现船舶速度和碰撞角度对破冰效果有显著影响。船舶速度越大，冰脊受到的冲击力越大，破碎程度越严重。碰撞角度越大，冰脊在碰撞方向上的受力面积越小，单位面积上的冲击力越大，冰脊更容易发生破碎。冰脊的高度和宽度也会影响破冰结果，冰脊高度越高、宽度越大，其抵抗船舶撞击的能力相对越强，破碎程度相对较小。3.2.3冰载荷结果分析分析不同工况下的冰载荷时程曲线，能够清晰地了解冰载荷的变化规律。在工况1中，第一次撞击时，冰载荷迅速上升，在0.05s时达到峰值12000kN，随后逐渐下降。第二次撞击时，冰载荷在0.15s时达到峰值10000kN，峰值较第一次撞击有所降低。第三次撞击时，冰载荷在0.25s时达到峰值8000kN，峰值进一步下降。这是因为随着撞击次数的增加，冰脊的结构逐渐被破坏，其抵抗船舶撞击的能力减弱，导致冰载荷峰值逐渐降低。工况2中，第一次撞击时，冰载荷在0.04s时达到峰值15000kN。第二次撞击时，冰载荷在0.13s时达到峰值13000kN。第三次撞击时，冰载荷在0.22s时达到峰值11000kN。第四次撞击时，冰载荷在0.31s时达到峰值9000kN。同样，随着撞击次数的增加，冰载荷峰值呈下降趋势。工况3中，第一次撞击时，冰载荷在0.06s时达到峰值13000kN。第二次撞击时，冰载荷在0.16s时达到峰值11000kN。第三次撞击时，冰载荷在0.26s时达到峰值9000kN。对比不同工况下的冰载荷峰值，发现船舶速度和碰撞角度对冰载荷峰值有显著影响。船舶速度越大，碰撞时的动能越大，冰载荷峰值越高。碰撞角度越大，冰脊对船舶的反作用力在船舶运动方向上的分量越大，冰载荷峰值也越高。冰脊的强度和尺寸也会影响冰载荷，冰脊强度越高、尺寸越大，冰载荷峰值相对越高。3.2.4吸能结果分析计算不同工况下船舶和冰脊的吸能情况，对于评估碰撞过程中的能量传递和结构损伤具有重要意义。在工况1中，船舶在三次撞击过程中的吸能分别为2000kJ、1500kJ和1000kJ，冰脊的吸能分别为1000kJ、800kJ和600kJ。随着撞击次数的增加，船舶和冰脊的吸能都逐渐减少，这是因为冰脊的结构在撞击过程中逐渐被破坏，其吸收能量的能力降低，同时船舶在碰撞后的剩余动能也逐渐减小。工况2中，船舶在四次撞击过程中的吸能分别为2500kJ、2000kJ、1500kJ和1000kJ，冰脊的吸能分别为1200kJ、1000kJ、800kJ和600kJ。同样，船舶和冰脊的吸能随着撞击次数的增加而逐渐减少。工况3中，船舶在三次撞击过程中的吸能分别为2200kJ、1800kJ和1400kJ，冰脊的吸能分别为1100kJ、900kJ和700kJ。对比不同工况下的吸能结果，发现船舶速度和碰撞角度对吸能有显著影响。船舶速度越大，碰撞时的动能越大，船舶和冰脊的吸能都越高。碰撞角度越大，船舶在碰撞过程中的能量转化越复杂，吸能情况也会发生变化。冰脊的强度和尺寸也会影响吸能，冰脊强度越高、尺寸越大，其吸收能量的能力相对越强，船舶的吸能相对减少。在船舶设计和极地航行安全评估中，应充分考虑这些因素，以提高船舶的抗冰能力和航行安全性。四、船-冰脊碰撞结构强度评估4.1破冰船与冰脊碰撞艏部结构动力响应分析以某型号破冰船为研究对象，该破冰船专为极地航行设计，船长130m，型宽25m，型深12m，排水量10000吨。船艏采用特殊的球鼻艏结构，这种结构能够有效降低船舶在冰区航行时的阻力，提高破冰能力。球鼻艏的形状经过优化设计，其曲率半径和角度等参数都是根据大量的实验和模拟结果确定的，在与冰脊碰撞时，球鼻艏能够将碰撞力分散，减少对船体结构的集中冲击。船艏结构由船壳板、肋骨、纵骨和舱壁等构件组成，船壳板厚度为30mm，采用高强度合金钢材料，这种材料具有优异的强度和韧性，能够在低温环境下保持良好的力学性能。肋骨间距为600mm，纵骨间距为300mm，舱壁采用加厚设计，厚度达到40mm，以增强船艏结构的整体强度和稳定性。破冰船材料选用高强度合金钢，其主要材料属性如下：弹性模量E=210GPa，泊松比\nu=0.3，屈服强度\sigma_y=400MPa，密度\rho=7850kg/m^3。这些材料属性是通过严格的材料测试和实验确定的，确保在船冰碰撞的复杂工况下，材料能够满足结构强度和稳定性的要求。在低温的极地环境中，材料的弹性模量和屈服强度会发生一定变化，在数值模拟中，考虑了温度对材料属性的影响，通过实验数据拟合得到材料属性随温度变化的函数关系，在模拟过程中根据实际温度条件实时调整材料属性，以保证模拟结果的准确性。在船冰模拟中，设定如下边界条件：在破冰船的艉部，约束其三个方向的平动自由度和三个方向的转动自由度，模拟船舶在与冰脊碰撞时艉部相对固定的状态。这是因为在实际碰撞过程中，艉部受到船舶整体惯性和推进系统的作用，其运动相对较小，通过这种边界条件的设定，能够更准确地模拟船艏与冰脊碰撞时的受力和变形情况。在冰脊的底部，约束其三个方向的平动自由度，模拟冰脊与海底或冰层的连接状态。冰脊在海洋中虽然处于漂浮状态，但底部与周围的冰层或海底存在一定的相互作用，通过约束底部平动自由度，能够考虑这种相互作用对冰脊力学行为的影响。在船冰接触面上，采用罚函数法来模拟船冰之间的接触力，设置接触刚度为10^8N/m，摩擦系数为0.1。罚函数法是一种常用的接触模拟方法，通过在接触面上施加罚函数来模拟接触力，接触刚度和摩擦系数的取值是根据相关研究和实验结果确定的，能够较好地模拟船冰之间的接触和摩擦行为。4.2船-冰脊碰撞艏部结构动态响应分析在船舶与冰脊碰撞的过程中，船艏结构会经历复杂的动态响应，包括结构损伤、冰载荷变化以及能量吸收等方面，这些响应对于评估船舶的安全性和结构完整性至关重要。在碰撞过程中，船艏结构会出现不同程度的损伤。船壳板作为船艏的最外层结构，直接承受冰脊的撞击力，容易发生局部凹陷和撕裂。当船舶以较高速度与冰脊碰撞时，船壳板在接触区域会产生明显的凹陷，凹陷深度可达数十厘米，甚至会出现撕裂裂缝，导致船壳板的水密性受到破坏。在一次模拟船舶以10m/s的速度与冰脊碰撞的数值计算中，船壳板在碰撞后0.1s时，接触区域的凹陷深度达到了30cm，且在凹陷边缘出现了长度约为1m的撕裂裂缝。肋骨作为支撑船壳板的重要构件，在碰撞时会发生弯曲和断裂。由于肋骨需要承受船壳板传递的冰脊撞击力，当撞击力超过肋骨的承载能力时，肋骨就会发生弯曲变形，严重时会断裂。在船舶与冰脊多次碰撞的情况下，肋骨的损伤会逐渐累积，进一步削弱船艏结构的强度。在某一工况下，经过三次碰撞后，部分肋骨的弯曲角度达到了15°，且有两根肋骨发生了断裂。纵骨主要承受纵向的载荷，在碰撞过程中，纵骨可能会发生屈曲和断裂。冰脊的撞击力会使纵骨受到轴向压力和弯曲力矩的作用，当这些力超过纵骨的临界屈曲载荷时，纵骨就会发生屈曲变形，进而导致断裂。在一些数值模拟中发现，当船舶与冰脊的碰撞角度较小时，纵骨更容易发生屈曲和断裂，因为此时纵骨受到的轴向压力更大。冰载荷在船-冰脊碰撞过程中呈现出复杂的变化规律。在碰撞初期，冰载荷会迅速上升，这是因为船舶与冰脊刚接触时，碰撞速度较大，能量瞬间释放，导致冰载荷急剧增加。在0.01s内，冰载荷就可能从0上升到峰值的50%以上。随着碰撞的持续，冰脊开始发生破碎和变形，吸收部分能量，冰载荷会逐渐下降。当冰脊破碎到一定程度后，冰载荷会趋于稳定，但仍会在一定范围内波动，这是由于船艏与冰脊破碎后的碎冰继续相互作用，以及船舶自身的运动等因素导致的。在碰撞后期，冰载荷的波动范围可能在峰值的20%-30%之间。船舶与冰脊碰撞速度对冰载荷峰值有显著影响，碰撞速度越大，冰载荷峰值越高。当船舶碰撞速度从5m/s增加到10m/s时，冰载荷峰值可能会增加1-2倍。这是因为碰撞速度增加，船舶的动能增大，在碰撞瞬间传递给冰脊的能量也增多，从而导致冰载荷峰值升高。碰撞角度也会影响冰载荷，当碰撞角度增大时，冰脊与船艏的接触面积减小，单位面积上的冰载荷增大，冰载荷峰值也会相应提高。当碰撞角度从30°增加到60°时，冰载荷峰值可能会提高30%-50%。在船-冰脊碰撞过程中，船舶和冰脊都会吸收能量。船舶吸收的能量主要用于船艏结构的变形和损伤，包括船壳板的凹陷、肋骨的弯曲、纵骨的屈曲等。冰脊吸收的能量则主要用于自身的破碎和变形。在碰撞初期，船舶的动能迅速转化为船艏结构的变形能和冰脊的破碎能，随着碰撞的进行，能量逐渐耗散。通过对不同工况下的数值模拟结果进行分析，发现船舶吸收的能量与碰撞速度和碰撞角度密切相关。碰撞速度越大，船舶吸收的能量越多；碰撞角度越大，船舶吸收的能量也会相应增加。在某一工况下，当碰撞速度为8m/s，碰撞角度为45°时，船舶吸收的能量为5000kJ；当碰撞速度提高到12m/s，碰撞角度不变时，船舶吸收的能量增加到8000kJ；若碰撞速度保持在8m/s，碰撞角度增大到60°，船舶吸收的能量则增加到6000kJ。冰脊吸收的能量也会随着碰撞速度和碰撞角度的变化而改变，碰撞速度越大、碰撞角度越大，冰脊吸收的能量越多。4.3破冰船-冰脊碰撞结构响应影响因素分析研究船速和冰脊厚度对碰撞结构响应的影响，有助于深入了解船冰碰撞的力学过程，为极地船舶的安全航行和结构设计提供科学依据。为探究船速对碰撞结构响应的影响，设置船舶速度分别为5节、10节、15节，其他参数保持不变，进行多组数值模拟。在模拟过程中，详细记录碰撞力、船艏结构变形和能量吸收等关键数据。碰撞力方面，随着船速的增加，碰撞力峰值显著增大。当船速为5节时，碰撞力峰值为3000kN；船速提升至10节，碰撞力峰值跃升至6500kN；船速达到15节时，碰撞力峰值更是高达10000kN。这是因为船速增加，船舶的动能增大，在与冰脊碰撞时，能够释放出更多的能量，从而导致碰撞力峰值大幅上升。船艏结构变形也随船速的增加而加剧。船速较低时，船艏结构的变形相对较小，主要表现为局部的轻微凹陷；随着船速的提高，船艏结构的变形范围逐渐扩大，凹陷深度增加，甚至出现船壳板的撕裂和肋骨的断裂。在15节船速下，船艏的最大凹陷深度达到0.5m，部分肋骨的弯曲角度超过30°，船壳板出现多处撕裂裂缝，长度最长可达1.5m。能量吸收情况同样受到船速的显著影响。船速越快，船舶和冰脊吸收的能量越多。当船速为5节时，船舶吸收的能量为1000kJ，冰脊吸收的能量为800kJ；船速提升至10节，船舶吸收的能量增加到3000kJ，冰脊吸收的能量达到2000kJ；船速达