«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2018, Vol. 39 Issue (10): 1589-1597  DOI: 10.11990/jheu.201704095
0

引用本文  

王庆凯, 雷瑞波, 李志军. 融冰期北极海冰单轴压缩强度的试验研究[J]. 哈尔滨工程大学学报, 2018, 39(10), 1589-1597. DOI: 10.11990/jheu.201704095.
WANG Qingkai, LEI Ruibo, LI Zhijun. Experimental study on the uniaxial compressive strength of the Arctic sea ice during melt season[J]. Journal of Harbin Engineering University, 2018, 39(10), 1589-1597. DOI: 10.11990/jheu.201704095.

基金项目

全球变化研究国家重大科学研究计划项目(2015CB953901);国家自然科学基金项目(51639003,51579028);南北极环境综合考察与评估专项(CHINARE-0301/0402)

通信作者

李志军, E-mail:lizhijun@dlut.edu.cn

作者简介

王庆凯(1991-), 男, 博士研究生;
李志军(1960-), 男, 教授, 博士生导师

文章历史

收稿日期:2017-04-26
网络出版日期:2018-06-06
融冰期北极海冰单轴压缩强度的试验研究
王庆凯1, 雷瑞波2, 李志军1    
1. 大连理工大学 海岸和近海工程国家重点实验室, 辽宁 大连 116024;
2. 中国极地研究中心 国家海洋局极地科学重点实验室, 上海 200136
摘要:为了评估融冰期北极海冰的单轴压缩强度,其中在中国第七次北极科学考察期间,现场钻取冰芯,测量冰层温度、盐度、密度和晶体结构;在低温实验室内将冰芯加工成标准试样并测量单轴压缩强度,其中试验应变速率为10-7~10-2 s-1,温度分别为-3℃、-6℃和-9℃,加载方向为垂直冰面加载。试验结果表明:融冰期北极海冰表现出高温、低盐、低密度的物理性质;随应变速率增加,单轴压缩强度先增加再降低最后稳定,试样破坏形态由主裂缝贯穿转变为产生较多微裂缝;随孔隙率增加,单轴压缩强度降低。试验建立了以应变速率和孔隙率评估海冰单轴压缩强度的计算模型,通过冰层物理参数计算孔隙率,给出了融冰期北极冰层垂直方向的极限单轴压缩强度参考值。
关键词海冰    单轴压缩强度    应变速率    孔隙率    北极    融冰期    试验研究    
Experimental study on the uniaxial compressive strength of the Arctic sea ice during melt season
WANG Qingkai1, LEI Ruibo2, LI Zhijun1    
1. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China;
2. Key Laboratory for Polar Science of the State Oceanic Administration, Polar Research Institute of China, Shanghai 200136, China
Abstract: To estimate the uniaxial compressive strength of the Arctic sea ice during melt season, the ice temperature, salinity, density, and crystal structure were measured during the seventh Chinese National Arctic Research Expedition. The ice cores were processed into standard specimens to test the uniaxial compressive strength with a strain rate of 10-7 s-1 to 10-2 s-1. The test temperatures were set at -3℃, -6℃, and -9℃, and the loading direction was vertical to the ice surface. The results showed that the Arctic sea ice exhibits high temperature, low salt, and low density during the melt season. With increasing strain rate, the uniaxial compressive strength first increased, and then declined, and finally stabilized, and the failure modes were transformed from the main breakthrough crack to micro cracks. With an increase in the porosity, the uniaxial compressive strength decreased. An equation was established to calculate the uniaxial compressive strength using strain rate and porosity. Based on the porosities of the ice covers determined using measured physical properties, the uniaxial compressive strength of Arctic ice cover in vertical direction during the melt season is inferred.
Keywords: sea ice    uniaxial compressive strength    strain rate    porosity    Arctic    melt season    experiment study    

近年来,随着全球气候变暖,北极海冰持续减少[1-2],使得北极通航的可能性大大增加[3-4],北极航行的相关领域也愈发受到国际瞩目。北极航道是连接亚欧之间最短的海上通道,其航运经济价值毋庸置疑[5]。而且北冰洋蕴藏丰富的石油和天然气资源,对世界能源格局的重分布具有重要影响[6]。尽管北极航行集中于夏季融冰期,但是海面仍漂浮大量浮冰,严重威胁过往船舶的安全,因此研发安全可靠的高新技术破冰船迫在眉睫。冰对船的作用力是破冰船设计需要的主要荷载。破冰船冰荷载的研究是一个极其复杂的问题,一方面需要了解冰-船之间的响应机制,另一方面需要详实的冰的原位物理力学参数[7]。我国对于破冰船冰荷载的研究尚处于起步阶段。郭春雨等[8-10]依托船舶拖曳水池,进行了破冰船在碎冰中航行阻力方面的模型试验和数值模拟研究,王超等[11]提出了破冰船在平整冰中航行性能预报的数值方法;对于原位海冰物理力学性质研究,雷瑞波等[12-13]对北极海冰物理性质进行了现场调查,韩红卫[14]对北极海冰力学行为进行了初步测试。国外对北极海冰物理力学性质的研究起步较早,Perovich等[15-16]分别对北极海冰的温度和盐度变化进行了长期观测,Timco[17]总结并给出了北极不同类型海冰密度的一般分布范围;Chen等[18-21]对北极不同地区海冰强度进行原位和实验室试验,并给出了海冰发生韧性破坏时的强度规律。虽然对北极海冰的物理及力学性能已有一定探究,但仍缺少对北极海冰快速融化新形式下其物理及力学性质相应变化的认识。

破冰船的连续破冰过程可以认为是“接触-挤压-弯曲”的反复循环[22]。其中,海冰单轴压缩强度是冰-船挤压过程中的重要力学参数。海冰的单轴压缩强度受外部加载条件和海冰内部物理性质(温度、盐度、密度和晶体结构)的影响。应变速率是影响海冰单轴压缩强度的主要外部因素,国内外学者做了大量单轴压缩强度的应变速率效应相关方面研究[23-25]。对于控制海冰单轴压缩强度的物理参数,卤水体积曾被视为关键因素[26]。然而,对于高温海冰,由于卤水排泄严重,卤水泡被空气占据,冰内空气含量增加,仅使用卤水体积已经不能准确反映内部物理状态对其力学行为的影响[27]。海冰孔隙率是卤水体积比和空气体积比之和,李志军等[28]

提出用孔隙率来评价海冰的单轴压缩强度,并将其成功地应用到渤海海冰强度的评估。此外,北极海冰快速消融,其物理性质随之变化,孔隙率由海冰温度、盐度和密度决定[29],可以综合反映北极环境快速变化下的海冰力学行为变化。对于破冰船设计,需要的是整个冰层的强度作为设计参数。受现场条件限制,冰层单轴压缩试验难以实施,Timco等[21]曾提出模型将评估试样单轴压缩强度的计算方法应用到评估整个冰层的强度,并得到现场试验的验证[18]

本文根据中国第七次北极科学考察期间实测冰层物理参数以及海冰试样单轴压缩强度,建立了以应变速率和孔隙率评估试样单轴压缩强度的计算模型;并参考Timco和Frederking模型[21],根据实测融冰期海冰物理参数,将试样强度计算方程应用到计算冰层强度,以期为北极航行破冰船设计提供参考。

1 现场观测及试验

北极海冰融冰期每年自4月开始至9月结束,并在9月达到当年海冰面积的最小值[30]。中国第七次北极科学考察利用雪龙船为平台,于2016年7月25日进入北极冰区(72.3°N,165.1°W)至8月24日退出冰区(73.8°N,168.8°W),冰区作业为期31天。走航期间进行6个短期冰站考察(SIC01~SIC06),最北达(82.7°N,166.9°W)并建立长期冰站(LIC01)。雪龙船冰区航线及短、长期冰站位置如图 1所示。短期冰站在表面平坦的平整冰区域钻取冰芯,现场测量冰芯温度和盐度,并运回船载低温实验室测量冰芯密度和晶体结构;长期冰站除测量冰芯物理性质外,钻取冰芯并在低温实验室内加工成标准单轴压缩试样测量其单轴压缩强度。

Download:
图 1 中国第七次北极考察冰区航迹及冰站位置 Fig. 1 Trajectory during the 7th Chinese Arctic expedition in the ice zone and locations of ice camps
1.1 冰层物理性质观测

现场用直径为9.0 cm的冰芯钻钻取冰芯后,立即测量冰芯温度,方法为用手电钻以10.0 cm间隔从冰芯表面钻孔至冰芯内部,将温度探针伸入孔内并用冰屑埋实,示数稳定后读取冰芯温度;盐度测量使用另一根完整冰芯,冰芯取出后以10.0 cm间隔切成段并装入塑料容器运回船上,冰芯段融化后用盐度计测量融冰水盐度。冰芯密度测量和晶体结构观测在低温实验室内完成。冰芯密度测量采用质量-体积法,将冰芯以10.0 cm间隔切成段并将冰芯段切出2块边长为4.0 cm的立方块,测量立方块质量并计算密度,取立方块密度平均值作为该冰芯段的密度;冰芯晶体结构观测先将冰芯以10.0 cm切成段,再将冰芯段加工成薄冰片并冻结在玻璃片上,薄冰片厚度切削至小于0.1 cm后放置在费氏台上,通过正交偏振光观测海冰晶体结构。

1.2 海冰试样单轴压缩试验

现场采集的冰芯装入冰芯袋运回船载低温实验室,在-15 ℃室温环境中保存。加工标准单轴压缩试样时,先将冰芯沿长度方向以20.0 cm间隔切成冰芯段,之后用冰样专用车床将冰芯段直径加工至7.0 cm,用锯骨机将冰芯段两端切除并加工成长度为17.5 cm的标准圆柱体试样。根据试样在冰层中的位置以及冰层晶体观测结果判断试样晶体类型并分别保存于恒温箱内。

单轴压缩试验在低温实验室内进行。试验设备如图 2所示。试验压力机使用液压油泵加载系统,并配有力传感器和激光位移传感器同步记录加载荷载和冰样位移。将调速电机与油泵节流阀连接,通过控制调速电机转速调节试验加载速率。试验前预先将试样置于恒温箱中恒温至设计试验温度24 h以上,试验时调节低温实验室室温至试验温度。为计算压缩试样的孔隙率,每个试样加载前用游标卡尺和电子天平测量直径、高度和质量;加载后将破碎块装入玻璃容器,融化后用盐度计测量融冰水盐度。

Download:
图 2 试样单轴压缩试验设备 Fig. 2 Equipment of the uniaxial compression test of samples

如果以融冰期冰层实际温度范围作为压缩试验设计温度,则试样处在相变状态不易于试验进行,因此试验仍选择在相对较低的温度条件下进行测试。同时,为获得较宽的孔隙率范围,设计试验温度为-3 ℃、-6 ℃和-9 ℃。柱状冰为各向异性材料,沿其晶体长、短轴方向加载表现不同强度。北极冰层较厚,水平方向上大面积的完整取样比较困难,无法进行平行冰面方向的压缩试验。因此,对于采用冰芯钻竖直方向取样的方式,试样加载方向为垂直冰面方向。海冰在低、高应变速率下表现不同的力学性质,因此设计试验应变速率区间为10-7~10-2 s-1,使冰试样表现出不同的力学行为。

力传感器和激光位移传感器同步记录加载荷载-时间曲线和试样位移-时间曲线。如图 3所示,当荷载达到极大值Ff时,试样发生破坏,失去承载力,对应破坏时间tf和破坏位移Df。定义试样单轴压缩强度σ=Ff/S,应变速率$\dot \varepsilon $=Df/(L·tf),其中S为试样横截面积,L为试样高度。

Download:
图 3 试样单轴压缩试验荷载-时间、位移-时间曲线 Fig. 3 Temporal variations of the load and displacement during a uniaxial compression test
2 试验结果及分析 2.1 冰层物理性质

北极融冰期冰层表现高温、低盐、低密度的特点。图 4给出了各短、长期冰站冰层温度、盐度和密度的观测结果。融冰期冰层温度较高,在0~-1.0 ℃变化。冰层表层冰温受气温影响而波动,SIC01~04冰站表层冰温为-0.2 ℃左右,SIC05、06和LIC01冰站表层冰温由于气温较低而小于-0.5 ℃,冰层中下部冰温随深度增加而降低,底层冰温为-1.0 ℃左右;融冰期冰层盐度较低,由于表层融化剧烈,卤水下渗,导致盐度随深度增加,表层盐度接近于0,底层盐度小于4.0 PSU。冰层密度随深度增加而增加,平均密度为829.9 kg/m3

Download:
图 4 冰层温度、盐度和密度剖面 Fig. 4 Profiles of temperature, salinity and density of ice floes at ice camps

海冰的晶体结构反映了海冰的生长状态,也决定海冰是否为各向异性。图 5给出了各冰站冰芯晶体结构竖直剖面的观测结果。冰芯晶体以柱状冰为主。SIC04冰站表层有10.0 cm左右粒状冰层,其余冰站表层粒状冰层可能融化殆尽。SIC02底层冰芯和LIC01中层冰芯由于生长期温度变化而形成粒状冰。

Download:
图 5 冰层晶体结构剖面 Fig. 5 Profiles of crystal texture of ice floes at ice camps
2.2 海冰试样单轴压缩试验 2.2.1 试样破坏特征

每一单轴压试样均可得到一条应力-应变曲线,不同应变速率对应不同形式的曲线。与渤海海冰类似,高应变速率加载对应脆性破坏特征,低应变速率加载对应韧性破坏特征。如图 6(a)所示,当应变速率较快时,加载后应力随应变线性上升,至峰值后立刻回落,试样破坏时没有明显变形产生,多为一条主裂缝贯穿试样;当应变速率较慢时,如图 6(b),试样加载后先表现出弹性特征,应力随应变线性增加,在弹性变形达到极限后应力进入屈服阶段,此时试样内部产生许多微小裂纹,屈服阶段结束之后应力下降,试样由于裂纹过多而崩坏,产生的变形也较大。值得注意的时,由于融冰期北极海冰内部孔隙多,裂缝得以充分发展,因此与渤海海冰相比,韧性破坏时融冰期北极海冰有更长的屈服阶段。

Download:
图 6 单轴压缩试验应力-应变曲线及其破坏形态 Fig. 6 Stress-strain curves and failure modes of the uniaxial compression test
2.2.2 试样单轴压缩强度与应变速率的关系

试样单轴压缩强度随应变速率变化而发生变化。如图 7所示,根据试样压缩破坏形式可将应变速率分为3个分区:试样均为韧性破坏的韧性区、韧性破坏和脆性破坏同时存在的韧脆转变区和试样均为脆性破坏的脆性区。图 7将试验数据汇总于强度和应变速率双对数坐标系,并给出强度上下限。韧性区试样单轴压缩强度随应变速率的增加呈增加趋势,韧脆转变区试样强度随应变速率的增加呈减少趋势,脆性区强度随应变速率增加呈稳定趋势。试样强度受试样内部物理条件的影响在上下限值之间波动,极限强度值出现在由韧性区向韧脆转变区过渡的临界应变速率。

Download:
图 7 试样单轴压缩强度与应变速率的关系 Fig. 7 The relationship between the uniaxial compressive strength of samples and strain rate
2.2.3 试样单轴压缩强度与孔隙率的关系

孔隙率是影响海冰单轴压缩强度的主要内部物理条件。Cox等[29]根据海冰相成份,利用海冰温度、盐度和密度推导出-2.0 ℃~-22.9 ℃的孔隙率拟合公式;Leppäranta在此基础上,将拟合式的适用范围扩展至高温(0~-2.0 ℃)海冰[31]

$ v = {v_a} + {v_b} $ (1)
$ {v_a} = 1 - \rho /{\rho _i} + \rho {S_i}{F_2}\left( T \right)/{F_1}\left( T \right) $ (2)
$ {v_b} = \rho {S_i}/{F_1}\left( T \right) $ (3)
$ {\rho _i} = 0.917 - 1.403 \times {10^{ - 4}}T $ (4)

式中:v为孔隙率,va为气泡体积比,vb为卤水体积比,ρ为海冰密度,ρi为纯冰密度,Si为海冰盐度,T为海冰温度,F1(T)和F2(T)为关于T的三次多项式,其系数由表 1给出。

表 1 F1(T)和F2(T)的多项式系数 Tab.1 Coefficients for the polynomials of F1(T) and F2(T)

根据试样在压缩试验前后测量的密度和盐度及对应的试验温度可以得到试样的孔隙率,其中柱状冰试样盐度和密度分别为0~2.8 PSU和756.2~989.8 kg/m3,粒状冰试样盐度和密度分别为0.7~2.9 PSU和759.7~934.2 kg/m3图 8汇总了试样单轴压缩强度和孔隙率的关系,并给出强度上限。试样强度随孔隙率的增加呈减少趋势。当试样孔隙率为20‰时,柱状冰和粒状冰试样强度极大值均达到10 MPa左右;当孔隙率增加至100‰左右时,压缩试样强度锐减,柱状冰极大强度仅在4 MPa左右,粒状冰极大强度仅在5 MPa左右。

Download:
图 8 试样单轴压缩强度与孔隙率的关系 Fig. 8 The relationship between the uniaxial compressive strength of samples and porosity
2.2.4 海冰试样单轴压缩强度计算模型

海冰单轴压缩强度受内部物理条件和外部加载条件的影响,学者们一直试图选用合适的内外影响因子来评估海冰的单轴压缩强度。应变速率作为外部加载条件的代表因素已经被普遍认可;海冰内部物理条件最初选用温度单因子,后来选用由温度和盐度表示的卤水体积,现在引入由温度、盐度和密度确定的孔隙率,逐渐完善海冰内部物理条件与力学行为之间的联系。

Kovacs[20]曾参考Glen公式[32]并认为强度与孔隙率之间呈幂函数关系,建议用如下公式来描述单轴压缩强度、应变速率和孔隙率之间的关系[27]

$ \sigma = A{{\dot \varepsilon }^n}{v^m} $ (5)

式中:Anm为经验系数。考虑到当孔隙率为0时,式(5)无法保证其数学合理性,建议用(1+ν)代替ν来表述孔隙率与单轴压缩强度的关系。试验应变速率覆盖韧性区、韧脆转变区和脆性区,韧性区强度随应变速率增加而增加,且认为二者符合Glen公式[32];韧脆转变区强度随应变速率增加而减小;脆性区强度不随应变速率发生变化;整个应变速率区间内强度随孔隙率的增加而减小,且认为二者为幂函数关系。因此,结合不同应变速率区间内海冰试样单轴压缩强度与应变速率和孔隙率的关系,给出海冰试样单轴压缩强度的计算模型:

$ \begin{array}{c} \frac{\sigma }{{{\sigma _1}}} = \\ \left\{ \begin{array}{l} {A_d}{(\dot \varepsilon /{{\dot \varepsilon }_1})^{{n_d}}}{(1 + v/{v_1})^{{m_d}}},\dot \varepsilon \le {{\dot \varepsilon }_d}\\ \left[ {{A_t} + \frac{{{B_t}}}{{\dot \varepsilon /{{\dot \varepsilon }_1}}} + \frac{{{C_t}}}{{{{(\dot \varepsilon /{{\dot \varepsilon }_1})}^2}}}} \right] \cdot {(1 + v/{v_1})^{{m_t}}},\\ \quad \quad \quad \quad \quad \quad \quad \quad {{\dot \varepsilon }_d} < \dot \varepsilon \le {{\dot \varepsilon }_b}\\ {A_b}(1 + v/{v_1}),\quad \quad \quad \quad \dot \varepsilon > {{\dot \varepsilon }_b} \end{array} \right. \end{array} $ (6)

式中:σ1=1 MPa,${{{\dot \varepsilon }_1}}$=1 s-1v1=1‰,${{\dot \varepsilon }_d}$${{\dot \varepsilon }_b}$分别为由韧性区向韧脆转变区过渡和由韧脆转变区向脆性区过渡的临界应变速率;ABCnm为拟合系数;下标dtb分别代表韧性区、韧脆转变区和脆性区。拟合系数取值及拟合优度由表 2给出。

表 2 式(6)拟合系数及拟合优度R2(p < 0.05) Tab.2 The fitting coefficients and coefficients of determination for Eq. (6) (p < 0.05)

拟合系数nm分别代表应变速率和孔隙率对试样单轴压缩强度的影响。由表 2发现,应变速率对柱状冰试样韧性区强度的影响小于对粒状冰韧性区强度的影响;孔隙率对柱状冰韧性区、韧脆转变区和脆性区强度的影响相同,对粒状冰韧性区、韧脆转变区和脆性区强度的影响程度依次递减,且孔隙率对柱状冰强度的影响比对粒状冰强度的影响大。

由此,建立了用应变速率和孔隙率评估海冰试样单轴压缩强度的计算模型,并给出了试样强度随应变速率和孔隙率变化的三维曲面,如图 9。当应变速率为${{\dot \varepsilon }_d}$时,试样表现出极限强度;当应变速率大于${{\dot \varepsilon }_b}$时,试样表现出脆性强度,${{\dot \varepsilon }_d}$${{\dot \varepsilon }_b}$不随孔隙率发生移动,其取值见表 3

Download:
图 9 试样单轴压缩强度随应变速率和孔隙率的变化曲面 Fig. 9 The surfaces of the uniaxial compressive strength of samples varying with strain rate and porosity
表 3 试样极限强度和脆性强度临界应变速率 Tab.3 The critical strain rates of ultimate strength and brittle strength of samples
3 评估融冰期冰层单轴压缩强度

以上建立了评估海冰试样单轴压缩强度计算方程,根据Timco和Frederking模型[21],由实测冰层温度、盐度和密度,可以将方程应用到评估冰层尺度上的单轴压缩强度。

3.1 Timco和Frederking模型

Timco和Frederking在总结众多现场海冰单轴压缩强度试验结果和海冰试样单轴压缩强度计算方程的基础上,提出以试样单轴压缩强度评估冰层单轴压缩强度的模型。Timco和Frederking首先根据试验结果建立用孔隙率和应变速率计算试样强度的参数化方程,之后将冰层沿厚度等分9层,每层分别计算孔隙率并计算强度,9层强度的平均值作为整个冰层的强度。由于缺少实测冰层物理参数,Timco和Frederking将冰层温度、盐度和密度与气温和冰厚建立关系,通过后者计算孔隙率。

3.2 利用实测数据评估冰层单轴压缩强度

本文依据现场实测海冰物理参数计算各短、长期冰站冰层的孔隙率。将冰层沿厚度等分9层,由图 4所示冰层物理参数测量结果,各分层会包含多个测量段,取分层内所有测量段冰温、盐度和密度的平均值,采用Leppäranta给出的高温海冰孔隙率计算公式[31],计算各分层的平均孔隙率。图 10给出各短、长期冰站冰层的平均孔隙率剖面。

Download:
图 10 冰层平均孔隙率剖面和极限压缩强度剖面 Fig. 10 Profiles of average porosity and ultimate uniaxial compressive strength of ice floes at ice camps

柱状冰和粒状冰试样强度计算公式的拟合系数不同,因此需要根据图 5冰层晶体结构观测结果判断各分层海冰的晶体类型。当分层为柱状冰、粒状冰混合层时,以两种晶体类型的各自占比为权重计算强度加权平均作为该层冰层的强度值。

由此,根据海冰孔隙率剖面以及与晶体类型相应的强度计算公式,再选用合适的应变速率,即可评估冰层垂直冰面方向的单轴压缩强度。当冰层与船舶或海洋结构物发生碰撞时,应变速率是浮冰速度、接触宽度和冰层厚度等多因素共同影响下的参数。尽管已有学者提出应变速率的计算方法[33-34],但仍没有得到一致的结论。在破冰船和抗冰结构物的设计中,需要选取最不利情况为设计保障,即选择冰层的极限强度。根据表 3所列海冰试样极限强度对应的应变速率,代入式(6)计算各分层海冰极限强度,取各分层极限强度的平均值作为整个冰层的极限单轴压缩强度。图 10给出了各短、长期冰站冰层的极限强度剖面,表 4列出了各短、长期冰站冰层垂直方向极限单轴压缩强度值。取各冰站冰层极限强度的平均值3.32 MPa作为垂直方向北极融冰期冰层极限单轴压缩强度值。

表 4 冰层垂直方向极限单轴压缩强度 Tab.4 The vertical ultimate uniaxial compressive strength of ice floes at ice camps
4 结论

1) 在全球气候变暖的条件下,北极海冰在融冰期表现出高温、低盐、低密度的特点,冰温介于-1.0~0.0 ℃,盐度介于0~4.0 PSU,平均密度为829.93 kg/m3

2) 应变速率和孔隙率是影响海冰单轴压缩强度的主要因子。强度随应变速率增加先增加再减小最后不变,随孔隙率的增加而减小。柱状冰和粒状冰分别在应变速率为1.2×10-5和1.9×10-5 s-1时达到极大值。应变速率对柱状冰韧性区强度的影响小于对粒状冰韧性区强度的影响;孔隙率对柱状冰强度的影响大于粒状冰强度的影响。

3) 建立用应变速率和孔隙率来评估海冰单轴压缩强度的计算公式,基于Timco和Frederking模型将该公式应用到评估整个冰层的强度,并根据实测冰层物理参数,得到融冰期冰层垂直方向极限单轴压缩强度为3.32 MPa。

4) 本次试验对融冰期北极海冰竖直方向的单轴压缩强度进行了测试,但对于抗冰结构而言,水平方向的压缩强度更为重要。未来的工作将着力于开发北极冰层水平方向取样的新技术以及探究海冰水平方向和竖直方向单轴压缩强度之间的关系。

参考文献
[1]
COMISO J C. Large decadal decline of the arctic multiyear ice cover[J]. Journal of climate, 2012, 25(4): 1176-1193. DOI:10.1175/JCLI-D-11-00113.1 (0)
[2]
DAY J J, HARGREAVES J C, ANNAN J D, et al. Sources of multi-decadal variability in Arctic sea ice extent[J]. Environmental research letters, 2012, 7(3): 034011. DOI:10.1088/1748-9326/7/3/034011 (0)
[3]
KHON V C, MOKHOV I I. Arctic climate changes and possible conditions of Arctic navigation in the 21st century[J]. Izvestiya, atmospheric and oceanic physics, 2010, 46(1): 14-20. DOI:10.1134/S0001433810010032 (0)
[4]
SMITH L C, STEPHENSON S R. New Trans-Arctic shipping routes navigable by midcentury[J]. Proceedings of the national academy of sciences of the United States of America, 2013, 110(13): E1191-E1195. DOI:10.1073/pnas.1214212110 (0)
[5]
LIU Miaojia, KRONBAK J. The potential economic viability of using the Northern Sea Route (NSR) as an alternative route between Asia and Europe[J]. Journal of transport geography, 2010, 18(3): 434-444. DOI:10.1016/j.jtrangeo.2009.08.004 (0)
[6]
朱明亚, 平瑛, 贺书锋. 北极油气资源开发对世界能源格局和中国的潜在影响[J]. 海洋开发与管理, 2015, 32(4): 1-7.
ZHU Mingya, PING Ying, HE Shufeng. The potential impact of Arctic oil and gas resources development on the world energy pattern and China[J]. Ocean development and management, 2015, 32(4): 1-7. DOI:10.3969/j.issn.1005-9857.2015.04.001 (0)
[7]
薛彦卓, 倪宝玉. 极地船舶与浮体结构物力学问题研究综述[J]. 哈尔滨工程大学学报, 2016, 37(1): 36-40.
XUE Yanzhuo, NI Baoyu. Review of mechanical issues for polar region ships and floating structures[J]. Journal of Harbin Engineering University, 2016, 37(1): 36-40. (0)
[8]
郭春雨, 谢畅, 王帅, 等. 碎冰条件下冰区船阻力性能试验研究[J]. 哈尔滨工程大学学报, 2016, 37(4): 481-486.
GUO Chunyu, XIE Chang, WANG Shuai, et al. Resistance of ships in pack ice conditions[J]. Journal of Harbin Engineering University, 2016, 37(4): 481-486. (0)
[9]
郭春雨, 李夏炎, 王帅, 等. 冰区航行船舶碎冰阻力预报数值模拟方法[J]. 哈尔滨工程大学学报, 2016, 37(2): 145-150, 156.
GUO Chunyu, LI Xiayan, WANG Shuai, et al. A numerical simulation method for resistance prediction of ship in pack ice[J]. Journal of Harbin Engineering University, 2016, 37(2): 145-150, 156. (0)
[10]
HUANG Yan, LI Wei, WANG Yinghui, et al. Experiments on the resistance of a large transport vessel navigating in the Arctic region in pack ice conditions[J]. Journal of marine science and application, 2016, 15(3): 269-274. DOI:10.1007/s11804-016-1362-x (0)
[11]
王超, 康瑞, 孙文林, 等. 平整冰中破冰船操纵性能初步预报方法[J]. 哈尔滨工程大学学报, 2016, 37(6): 747-753.
WANG Chao, KANG Rui, SUN Wenlin, et al. Preliminary forecasting method for icebreaker maneuverability on level and smooth ice[J]. Journal of Harbin Engineering University, 2016, 37(6): 747-753. (0)
[12]
LEI Ruibo, LI Zhijun, LU Peng, et al. Crucial physical characteristics of sea ice in the Arctic section of 143°~180°W during August and early September 2008[J]. Acta oceanologica sinica, 2012, 31(4): 65-75. DOI:10.1007/s13131-012-0221-0 (0)
[13]
HUANG Wenfeng, LEI Ruibo, ILKKA M, et al. The physical structures of snow and sea ice in the Arctic section of 150°-180°W during the summer of 2010[J]. Acta oceanologica sinica, 2013, 32(5): 57-67. DOI:10.1007/s13131-013-0314-4 (0)
[14]
HAN Hongwei, LI Zhijun, HUANG Wenfeng, et al. The uniaxial compressive strength of the Arctic summer sea ice[J]. Acta oceanologica sinica, 2015, 34(1): 129-136. DOI:10.1007/s13131-015-0598-7 (0)
[15]
PEROVICH D K, ELDER B C. Temporal evolution of Arctic sea-ice temperature[J]. Annals of glaciology, 2001, 33: 207-211. DOI:10.3189/172756401781818158 (0)
[16]
NAKAWO M, SINHA N K. Growth rate and salinity profile of first-year sea ice in the high Arctic[J]. Journal of glaciology, 1981, 27(96): 315-330. DOI:10.1017/S0022143000015409 (0)
[17]
TIMCO G W, FREDERKING R M W. A review of sea ice density[J]. Cold regions science and technology, 1996, 24(1): 1-6. DOI:10.1016/0165-232X(95)00007-X (0)
[18]
CHEN A C T, LEE J. Large-scale ice strength tests at slow strain rates[J]. Journal of offshore mechanics and arctic engineering, 1988, 110(3): 302-306. DOI:10.1115/1.3257066 (0)
[19]
SINHA N K. Field test 1 of compressive strength of first-year sea ice[J]. Annals of glaciology, 1983, 4: 253-259. DOI:10.1017/S0260305500005565 (0)
[20]
KOVACS A. Estimating the full-scale flexural and compressive strength of first-year sea ice[J]. Journal of geophysical research, 1997, 102(C4): 8681-8689. DOI:10.1029/96JC02738 (0)
[21]
TIMCO G W, FREDERKING R M W. Compressive strength of sea ice sheets[J]. Cold regions science and technology, 1990, 17(3): 227-240. DOI:10.1016/S0165-232X(05)80003-5 (0)
[22]
王钰涵, 李辉, 任慧龙, 等. 连续破冰模式下破冰船的冰力研究[J]. 海洋工程, 2013, 31(4): 68-73.
WANG Yuhan, LI Hui, REN Huilong, et al. Study of ice force about icebreaker based on continuous breaking pattern[J]. The ocean engineering, 2013, 31(4): 68-73. (0)
[23]
SINHA N K. Constant strain-and stress-rate compressive strength of columnar-grained ice[J]. Journal of materials science, 1982, 17(3): 785-802. DOI:10.1007/BF00540376 (0)
[24]
王庆凯, 卢燕, 郭可彩, 等. 辽河口附近海冰单轴压缩强度的试验与统计特征[J]. 数学的实践与认识, 2014, 44(24): 157-162.
WANG Qingkai, LU Yan, GUO Kecai, et al. Experiment and statistical characteristics on uniaxial compressive strength of sea-ice in Liaohe River estuary region[J]. Mathematics in practice and theory, 2014, 44(24): 157-162. (0)
[25]
韩红卫, 雷瑞波, 李丙瑞, 等. 夏季南极固定冰单轴压缩强度及试样破碎规律研究[J]. 极地研究, 2016, 28(4): 505-513.
HAN Hongwei, LEI Ruibo, LI Bingrui, et al. Study on the uniaxial compressive strength and crushing regulation of the landfast Antarctic sea ice in Austral summer[J]. Chinese journal of polar research, 2016, 28(4): 505-513. (0)
[26]
LAINEY L, TINAWI R. The mechanical properties of sea ice-A compilation of available data[J]. Canadian journal of civil engineering, 1984, 11(4): 884-923. DOI:10.1139/l84-106 (0)
[27]
FRANSSON L, STEHN L. Porosity effect on measured strength of warm ice[C]//Proceedings of the 12th International Conference on Port and Ocean Engineering under Arctic Conditions. Hamburg, 1993: 23-26. (0)
[28]
LI Zhijun, ZHANG Limin, LU Peng. Experimental study on the effect of porosity on the uniaxial compressive strength of sea ice in Bohai Sea[J]. Science China technological sciences, 2011, 54(9): 2429-2436. DOI:10.1007/s11431-011-4482-1 (0)
[29]
COX G F N, WEEKS W F. Equations for determining the gas and brine volumes in sea-ice samples[J]. Journal of glaciology, 1983, 29(102): 306-316. DOI:10.1017/S0022143000008364 (0)
[30]
王维波, 赵进平. 累积海冰密集度及其在认识北极海冰快速变化的作用[J]. 地球科学进展, 2014, 29(6): 712-722.
WANG Weibo, ZHAO Jinping. Accumulation sea ice concentration and its action on understanding arctic sea ice dramatic change[J]. Advances in earth science, 2014, 29(6): 712-722. (0)
[31]
LEPPÄRANTA M, MANNINEN T. The brine and gas content of sea ice with attention to low salinities and high temperatures[R]. Helsinki: Finnish Institute of Marine Research, 1988. (0)
[32]
GLEN J W. The creep of polycrystalline ice[J]. Proceedings of the royal society a:mathematical, physical and engineering sciences, 1955, 228(1175): 519-538. DOI:10.1098/rspa.1955.0066 (0)
[33]
PONTER A R S, PALMER A C, GOODMAN D J, et al. The force exerted by a moving ice sheet on an offshore structure:Part 1. The creep mode[J]. Cold regions science and technology, 1983, 8(2): 109-118. DOI:10.1016/0165-232X(83)90002-2 (0)
[34]
SODHI D S, HAEHNEL R B. Crushing ice forces on structures[J]. Journal of cold regions engineering, 2003, 17(4): 153-170. DOI:10.1061/(ASCE)0887-381X(2003)17:4(153) (0)