基于直接剪切試驗的海冰粘接力研究

于 馳，王 剛

（1．大連大學 環境與化工學院，遼寧 大連 116622；2. 大連大學 建筑工程學院，遼寧 大連 116622）

基于直接剪切試驗的海冰粘接力研究

于 馳1，王 剛2

（1．大連大學 環境與化工學院，遼寧 大連 116622；2. 大連大學 建筑工程學院，遼寧 大連 116622）

在渤海海冰動力學本構模型研究中，基于Mohr-Coulomb屈服準則的黏彈塑性本構模型已經得到了廣泛的應用。在此本構模型中，表征冰間黏結強度的粘接力是影響本構模型精度的一個重要參數，因次如何確定粘接力是一個重要的問題。本文首先通過理論分析，確定了影響粘接力的兩個重要參數：溫度和鹽度，并建立了粘接力與溫度和鹽度的理論公式，并在實驗室內，通過海冰試件的直接剪切試驗，得到了粘接力與海冰溫度和鹽度的試驗數據，通過數據擬合，分別得到了粘接力與溫度和鹽度的函數關系。由擬合函數可知，海冰的粘接力與海冰溫度呈線性關系，粘接力與海冰的鹽度呈負冪函數關系。

海冰；Mohr-Coulomb屈服準則；粘接力；溫度；鹽度；直接剪切試驗

0 引言

海冰是我國渤海區域冬季常見的一種自然現象。海冰的存在對海上航運、石油氣開發、水產養殖都帶來了較大的影響，由于海冰存在而導致的海冰災害也越來越受到人們的關注，國家海洋局每年發布的《中國海洋災害公報》中，也把海冰災害列為重要的災害預警內容之一。如2013年度、2014年度，由于海冰災害造成的直接經濟損失分別為3.22億和0.24億[1-2]。為了研究渤海海冰的動力運移規律，眾多學者在海冰動力學方面做了許多研究工作，其中，針對渤海海冰運動特點，季順迎等人[3]基于Mohr-Coulomb屈服準則，建立了中小尺度海冰動力學的黏彈-塑性本構模型，并將其應用到渤海海冰的數值模擬中，模擬結果表明該本構模型具有較好的可靠性。但是在此動力學本構模型中，表征冰間凍結作用的海冰粘接力由于影響因素較多，因此如何確定其表達式一直是一個難點問題?；诖?，本文首先通過理論分析，分別建立粘接力與海冰溫度和鹽度的函數關系，并通過實驗室內海冰的直接剪切試驗，得到海冰的剪切強度與溫度和鹽度的試驗數據，并將試驗數據與粘接力的理論函數相擬合，進而得到了粘接力的具體表達形式。

1 Mohr-Coulomb屈服準則

在 σ-τ平面上(如圖 1 所示)，Mohr-Coulomb 準則可寫作：

圖1 Mohr-Coulomb屈服準則

2 粘接力的理論確定

2.1 粘接力與溫度

溫度是影響粘接力最重要的因素之一。對于粘接力與溫度的關系，Fish等人采用式(2)或式(3)來表示粘接力與溫度之間的關系[4]：

這里，Tm是海冰的融解溫度，對于純冰而言等于0℃；c0是在融解溫度冰開始形成時所產生的粘接力分量；是一個特征參數。Vyalov[5]等還采用下式來描述溫度與粘接力之間的關系：

式中，b是一個經驗常數。

由公式(1)可知，粘接力與海冰的剪切強度成線性關系。而粘接力的本質也是海冰強度的一種表現，基于此，結合公式(2)，設定粘接力c與溫度T成線性關系，即：

a、d為待定常數，具體通過海冰直接剪切試驗確定。

2.2 粘接力與鹽度

在渤海海冰研究中，李志軍等[6]根據渤海海冰的現場試驗，統計得到海冰的剪切強度與孔隙率之間的關系如下：

當海冰的溫度在-0.5～-22.9 ℃之間時，孔隙率可利用鹵水體積νb來代替，其表達式為[7]：

式中，Si為海冰的鹽度。

季順迎等人[8-9]對環渤海海冰的剪切強度試驗結果也得到了類似的規律，并給出海冰的剪切強度與鹵水體積的計算公式為：

這里需要說明的是，式(6)～式(8)所描述的海冰強度指的是實驗室范圍內的微小尺度下的海冰力學性質。Zhang[10]等人指出，在不同的空間尺度下，運動的海冰具有不同的材料性質。對于漂浮在海面上的破碎海冰而言，不同的碎冰其力學性質是不同的。因此在海冰堆積研究中，海冰的力學參數指的是中小尺度下海冰的平均狀態下的參數。這就需要建立不同尺度下海冰力學性質的聯系。Lepparanta[11]給出了應力狀態與尺度之間一個簡單的比例關系 σ∝L-0.5。因此，假定在海冰堆積過程中，中小尺度的海冰力學參數與式(6)～式(8)具有相同的函數形式。因此，將粘接力表示冪函數形式為：

式中，f和g為經驗常數，需要通過海冰的直接剪切試驗來獲得。

3 直接剪切試驗

在實驗室內，采用帶側限的單面剪切試驗方法來獲取海冰試件的試驗數據。為了獲得冰溫與剪切強度的對應關系，采取必要的保溫與控溫措施。

試驗中采用兩塊厚鋼板來限制住冰塊在受剪切荷載作用下產生側向膨脹，防止冰樣底部產生拉應力，克服彎曲破壞，保證實驗的精確性。施壓方向采集剪切方向力和位移的大小。剪切實驗的加載方式和室內實驗裝置如圖 2所示。海冰試件剪切破壞后如圖 3所示。

圖2 海冰室內剪切實驗裝置

圖3 海冰試件直接剪切破壞圖

試驗進行時，由計算機采集壓頭壓力數據 T及壓頭的位移，并同步測量和記錄海冰試樣的溫度和鹽度。根據材料力學剪切強度公式計算獲得海冰的剪切強度：

式中，Tmax為海冰剪切破壞時的最大荷載，t、h分別為試樣受力面的長與寬。

試驗所用海冰試樣分別取自渤海遼東灣三個海冰監測區域，海冰試件尺寸為90 cm×90 cm×50 cm和90 cm×88 cm×50 cm，剪切面面積分別為90 cm×50 cm和88 cm×50 cm。冰溫按照試驗要求控制在-1～-10 ℃，并形成冰溫梯度，同時獲取每個冰樣的鹽度，試驗共測試49組海冰試件。

4 數據擬合

通過海冰粘接力的理論研究，并結合海冰室內直接剪切強度試驗，通過數據擬合，分別確定了海冰粘接力與海冰溫度和鹽度的函數形式。

4.1 粘接力與溫度關系確定

在純剪切狀態下，海冰的剪切強度即代表了粘接力。通過海冰的直接剪切試驗，可以得到海冰的粘接力與冰溫的關系如圖4所示。

圖4 海冰溫度與壓縮強度

結合公式(5)，通過試驗數據線性擬合，得到粘接力與溫度的擬合曲線為：

由上式可以看出，隨著海冰溫度的降低，冰間的黏結強度增大，從而使冰間粘接力增大。

4.2 粘接力與鹽度關系確定

通過海冰的剪切強度試驗，可以得到海冰的粘接力與鹽度的關系如圖5所示。

圖5 粘接力與海冰鹽度

并通過數據擬合，得到海冰的粘接力與鹽度之間的函數關系為：

由上式可以看出，隨著海冰鹽度的增大，冰塊之間的粘接力減小，海冰不易形成較大的塊體，以離散的形式存在。

5 結論

粘接力是Mohr-Coulomb型海冰動力學本構模型中的一個重要影響參數，由于其影響因素較多，因此一直未有明確的函數形式。本文針對海冰溫度和鹽度兩個對粘接力有影響的參數，通過理論分析給出了粘接力與溫度和鹽度的函數關系，并結合實驗室內獲得的海冰直接剪切試驗數據，擬合分析得到了粘接力與溫度和鹽度相應的函數關系式，擬合結果表明，海冰的粘接力與海冰溫度呈線性關系，海冰溫度越低則冰間的粘接力越強。粘接力與海冰的鹽度呈負冪函數關系，鹽度越高，冰間的粘結力越小。本文的研究結果為進一步完善粘彈塑性海冰動力學本構模型并提高海冰數值模擬的精度提供了一定的幫助。

[1]國家海洋局. 2014年中國海洋災害公報[R]. 北京∶ 2015.

[2]國家海洋局. 2013年中國海洋災害公報[R]. 北京∶ 2014.

[3]季順迎, 沈洪道, 王志聯, 等. 基于Mohr-Coulomb準則的黏彈—塑性海冰動力學本構模型[J]. 海洋學報∶ 中文版,2005(04)∶ 19-30.

[4]Fish A M, Zaretsky Y K. Temperature effect on strength of ice under triaxial compression[C]//Proc. of the 7th IOPE conference, USA∶ Honolulu, 1997∶ 415-421.

[5]Zaretskii Y K，Fish A M. Effect of temperature on the strength and viscosity of ice [J]. Soil Mechanics and Foundation Engineering, 1996, 33(2)∶ 46-52.

[6]李志軍, Devinder S S, 盧鵬. 渤海海冰工程設計參數分布[J]. 工程力學, 2006, 23(6)∶ 167-172.

[7]Frankenstein G, Garner R. Equations for determining the brine volume sea ice from -0.5℃ to -22.0℃[J]. Journal of Glaciology, 1967, 6(48)∶ 943-944.

[8]王安良, 許寧, 季順迎. 渤海沿岸海冰單軸壓縮強度的基本特性分析[J]. 海洋工程, 2014(04)∶ 82-88.

[9]季順迎, 李鵬飛, 劉宏亮, 等. 渤海萊州灣和黃河口附近海冰剪切強度的試驗研究[J]. 海洋通報, 2013(03)∶ 241-245.

[10]ZHANG Z H. On modeling ice dynamics of semi-enclosed seasonally ice-covered seas [R]. Report Series in Geophysics 43, Helsinki, 2000.

[11]Lepparanta M. Sea ice dynamics theory and downscaling [R].Report Series in Geophysics 40, Helsinki, 1998.

Study on Sea Ice Cohesion based on Direct Shear Test

YU Chi1, WANG Gang2
（1. College of Environmental and Chemical Engineering, Dalian University, Dalian 116622, China; 2. College of Civil and Architecture Engineering, Dalian University, Dalian 116622, China)

The viscoelastic-plastic constitutive model based on the Mohr-Coulomb yield criterion has been widely used in the study of dynamics constitutive model about Bohai sea ice. In this constitutive model, the cohesion which presents cohesion strength between ices is one of the important parameters affecting the accuracy of constitutive model. However, it is difficult to determine the cohesion. Firstly, the theoretical analysis indicates that the temperature and salinity are two main parameters for cohesion, and an expression of cohesion against temperature and salinity is presented. Moreover, through direct shear test of sea ice specimens in laboratory, the experiment data of shear strength against temperature and salinity are obtained. Through data fitting functions of cohesion against temperature and salinity are established respectively, which shows that cohesion increases linearly with ice temperature, but gives a negative power functions with the salinity of sea ice.

sea ice; Mohr-Coulomb yield criterion; cohesion, temperature; salinity; direct shear test

P731.15

A

1008-2395(2015)06-0017-04

2015-07-14

國家自然科學基金項目(40906011)；遼寧省教育廳資助項目（L2013473）。

于馳（1978-)，女，博士，講師，研究方向：流體力學數值計算研究。