劉婧蓉，謝良喜，熊 肖，胡 騰

（1.武漢科技大學(xué)機械自動化學(xué)院，湖北 武漢 430081；2.湖北誠宇信諾科技有限公司，湖北 武漢 430081）

1 引言

目前，護舷是各個港口和船舶所必不可少的承載裝置。當(dāng)船舶靠泊時，護舷通過緩沖、降噪以及隔震的方式，保護碼頭和船舶免受損壞。隨著我國經(jīng)濟和科技的發(fā)展，護舷由木質(zhì)和鋼制結(jié)構(gòu)演化為橡膠和鋼混合形式的結(jié)構(gòu)，且廣泛應(yīng)用于現(xiàn)代船舶制造業(yè)、航運業(yè)以及港口碼頭。但由于傳統(tǒng)橡膠具有易老化和顏色單調(diào)等特點，造成傳統(tǒng)橡膠護舷存在使用壽命縮短、不夠美觀及易在船舶表面形成擦痕等問題。聚氨酯作為一種新興橡膠材料，比傳統(tǒng)的橡膠更加耐磨、耐切割、抗撕裂強度高而且無需增強材料，制造上改變了傳統(tǒng)橡膠可能會出現(xiàn)的斷面、密度分布不均勻，且避免了復(fù)雜的硫化工藝等缺點。因此，用聚氨酯護舷取代橡膠護舷，具有廣闊的應(yīng)用前景。但聚氨酯與傳統(tǒng)橡膠具有不同的力學(xué)屬性，若直接沿用傳統(tǒng)橡膠護舷的截面形狀生產(chǎn)聚氨酯護舷，其反力與吸能性能能否達到相關(guān)標(biāo)準(zhǔn)要求，尚需要進一步研究[1]。

傳統(tǒng)的護舷結(jié)構(gòu)設(shè)計采用經(jīng)驗與試驗方法，有一定的局限性。目前，愈來愈多地采用有限元方法對橡膠護舷進行仿真分析，進而進行性能及壽命分析等，為優(yōu)化護舷產(chǎn)品的結(jié)構(gòu)設(shè)計提供有效的理論依據(jù)。但橡膠為非線性超彈性材料，護舷分析中面臨大變形非線性和接觸非線性的問題，使有限元分析面臨較大技術(shù)難度。ABAQUS 有限元分析軟件具有強大的非線性求解能力，文獻[2]根據(jù)材料拉伸試驗結(jié)果，采用OGDEN 三階模型定義橡膠的超彈性特性，建立鼓型橡膠護舷非線性有限元分析模型，并進行了試驗驗證和性能分析。文獻[3]采用附加質(zhì)量模型，對工程船的靠泊碰撞進行有限元數(shù)值模擬。研究了停泊船及橡膠護舷在靠泊船沖擊力下的響應(yīng)。文獻[4]采用有限元方法模擬了鋼板撞擊橡膠泡沫型護舷的全過程，分析了鋼板在撞擊過程中速度改變量、不同時刻護舷的應(yīng)變量以及護舷等效應(yīng)力。綜上所述，已有不少學(xué)者對傳統(tǒng)橡膠材料護舷進行有限元仿真分析并取得良好效果，但聚氨酯護舷延用傳統(tǒng)護舷截面形狀，是否能夠維持傳統(tǒng)橡膠良好的力學(xué)性能并投入工程使用需要進一步研究。為此，以SA400H-1000L 聚氨酯護舷為研究對象，首先對聚氨酯材料進行拉伸和壓縮試驗，根據(jù)試驗數(shù)據(jù)回歸預(yù)測聚氨酯材料性能的最優(yōu)本構(gòu)模型，其次建立聚氨酯護舷的有限元模型并開展非線性分析，最后對分析結(jié)果進行解讀并提出合理建議。

2 材料測試數(shù)據(jù)及其本構(gòu)模型

2.1 材料測試數(shù)據(jù)

聚氨酯屬于超彈性材料，其力學(xué)性能受材料配方比例、加載歷程、環(huán)境溫度等影響，導(dǎo)致應(yīng)力-應(yīng)變關(guān)系呈非線性。為了準(zhǔn)確定義聚氨酯的材料屬性，必須確立適合該材料的本構(gòu)模型。根據(jù)《GBT 7757-2009 硫化橡膠或熱塑性橡膠壓縮應(yīng)力-應(yīng)變性能的測定》和《GBT 528-2009 硫化橡膠或熱塑性橡膠拉伸應(yīng)力-應(yīng)變性能的測定》等相關(guān)規(guī)范，對聚氨酯材料進行單軸拉伸與單軸壓縮試驗，試驗圖，如圖1、圖2所示。

圖1 聚氨酯拉伸試樣與試驗Fig.1 Polyurethane Tensile Specimen and Test

圖2 聚氨酯壓縮試樣與試驗Fig.2 Polyurethane Compression Test Specimen and Test

通過試驗采集數(shù)據(jù)為20314組，其中19250組為拉伸試驗且各項數(shù)據(jù)精確到小數(shù)點后5位。試驗數(shù)據(jù)擬合應(yīng)力應(yīng)變曲線如圖2（a）所示。從圖2（b）觀察發(fā)現(xiàn)，應(yīng)變?yōu)椋?00～200）%區(qū)間試驗數(shù)據(jù)出現(xiàn)同一個應(yīng)變對應(yīng)多個應(yīng)力值。分析后有發(fā)現(xiàn)：應(yīng)變?yōu)?%以前符合胡克定律，應(yīng)力應(yīng)變成線性關(guān)系。屈服點后，應(yīng)變軟化（應(yīng)變增加較快而應(yīng)力增加較慢），出現(xiàn)分子鏈運動出現(xiàn)相對滑動，導(dǎo)致試樣與夾具之間可能發(fā)生了輕微相對位移。又因試驗機采集精度較高，故出現(xiàn)一個應(yīng)變對應(yīng)多個應(yīng)力值的現(xiàn)象。

因此，需要對拉伸數(shù)據(jù)進行篩選并再次擬合應(yīng)力-應(yīng)變曲線，分別在19250 組單軸拉伸試驗和1064 組單軸壓縮試驗數(shù)據(jù)中各篩選出100 組有效數(shù)據(jù)，進行了數(shù)據(jù)擬合，并與原始試驗數(shù)據(jù)進行對比，對比圖，如圖4、圖5 所示。對比發(fā)現(xiàn)篩選數(shù)據(jù)所擬合應(yīng)力-應(yīng)變曲線與原數(shù)據(jù)曲線近乎完全重合，具有良好的一致性。

圖4 篩選拉伸數(shù)據(jù)和壓縮數(shù)據(jù)所得應(yīng)力應(yīng)變曲線圖Fig.4 Stress-Strain Curves Obtained by Screening Tensile Data and Compression Data

圖5 單軸壓縮和單軸拉伸材料評估圖Fig.5 Material Evaluation Diagrams for Uniaxial Compression and Uniaxial Tension

2.2 本構(gòu)模型選取及參數(shù)設(shè)置

聚氨酯材料屬于非線性材料，而非線性材料力學(xué)性能常用描述方法主要有兩類，一類是基于熱力學(xué)統(tǒng)計的方法，其中包括Arruda-Boyce模型和Van der Waal模型；另一類是基于非線性材料為連續(xù)介質(zhì)的唯象學(xué)描述方法，其中包括N 次多項式模型（Polynomical）和Ogden 模型。對單軸拉伸與單軸壓縮的篩選試驗數(shù)據(jù)進行擬合，找出模擬聚氨酯材料的最佳本構(gòu)模型。四種模型的預(yù)測結(jié)果，如圖5所示。

結(jié)果表明二階多項式模型（POLY_N2）擬合出的應(yīng)力-應(yīng)變曲線變化趨勢與試驗數(shù)據(jù)吻合較好，表明該本構(gòu)模型可以較為精確地描述材料的力學(xué)性能；Mooney-Rivlin模型（POLY_N1）、三階Ogden模型（OGDEN_N3）和YEOH 模型（R_POLY_N3）擬合出的應(yīng)力-應(yīng)變曲線與試驗數(shù)據(jù)相差較大，不適合用于描述聚氨酯彈性體的力學(xué)性能。因此，選用多項式本構(gòu)模型的二階模型，來描述該配方聚氨酯彈性體的本構(gòu)關(guān)系，此模型的具體系數(shù)，如表1所示。根據(jù)《GB/T 533-91硫化橡膠密度測定》測出的聚氨酯材料的密度為1070kg/m3。

表1 評估參數(shù)系數(shù)表Tab.1 Evaluation Parameter Coefficient Table

3 超級護舷有限元建模及分析

3.1 幾何模型及邊界條件

對聚氨酯護舷進行非線性分析時，需要考慮三類因素。其一護舷在與船體碰撞過程中會產(chǎn)生較大變形，最大變形量可達到52.5%；其二聚氨酯材料為超彈性材料具有不可壓縮性和明顯的材料非線性；其三典型的邊界非線性問題，即超級護舷內(nèi)部自接觸不連續(xù)。只有充分考慮以上三類因素才可以真實模擬超級拱形護舷的大變形。超級拱形護舷模型是根據(jù)《HG/T 2866-2016》結(jié)構(gòu)要求所建立的，護舷總高度為400mm，底部厚度為40mm，護舷頂部寬度為300mm，具體模型，如圖6所示。

圖6 模型和網(wǎng)格劃分示意圖Fig.6 Schematic Diagram of Model and Meshing

上端采用解析剛體，將護舷底部固定，在剛體與護舷頂部設(shè)置接觸并對護舷內(nèi)部設(shè)置自接觸，剛體垂直向下壓縮護舷。采用Abaqus內(nèi)部模塊進行網(wǎng)格劃分并求解分析。其中不可壓縮聚氨酯材料采用CPS4R單元模擬，解析剛體采用CPS3單元模擬，整體模型的單元數(shù)為2774。

3.2 計算結(jié)果及分析

分析不同壓縮變形量下聚氨酯超級拱形護舷的反力，得到反力-變形量的關(guān)系,如圖7 所示。超級拱形護舷在變形量為52.5%時的應(yīng)力分布云圖，如圖8所示。

圖7 52.5%護舷壓縮變形應(yīng)力云圖Fig.7 Stress Cloud Diagram of 52.5% Fender Compression Deformation

圖8 聚氨酯超級拱形護舷性能曲線圖Fig.8 Performance Graph of Polyurethane Super Fender

由圖7可看出MISES應(yīng)力主要集中在階梯下方60mm處，最大應(yīng)力值約為80MPa，而該聚氨酯材料的抗拉強度小于47MPa，表明該護舷產(chǎn)品在壓縮過程中局部最大應(yīng)力過大。

由圖8可知，在護舷壓縮變形量為（0～17.64）%時，吸能量緩慢增加。在此之后吸能快速增加，至變形量為52.5%時，吸能達到192.42741kN ?m，符合標(biāo)準(zhǔn)《HG/T 2866-2016》的要求（52.5%變形量時，吸能應(yīng)大于46kN ?m）。

護舷的反力變化為三個階段：第一階段，反力穩(wěn)定增加，在17.64%變形時達到第一次峰值1703.21kN。進入第二階段，即變形量為（17.64～19.95）%之間，反力變化較為穩(wěn)定，反力值在（1703.21～1689.91）kN之間波動。在變形量為19.95%以后，由于聚氨酯超級拱形護舷壓縮變形，導(dǎo)致反力迅速減少，當(dāng)變形量為52.5%時，反力最小為426.98392kN，不符合標(biāo)準(zhǔn)《HG/T 2866-2016》要求（52.5%變形量時，反力應(yīng)小于275.0kN）。具體數(shù)據(jù)，如表2所示。

表2 護舷變形反力表Tab.2 Fender Deformation Reaction Table

4 結(jié)論

（1）采用二階多項式模型（POLY_N2）作為聚氨酯超級拱形護舷的本構(gòu)模型，是模擬聚氨酯超級拱形護舷的最佳本構(gòu)模型。

（2）該設(shè)計方案的聚氨酯超級拱形護舷最大應(yīng)力主要集中在階梯下方60mm處，最大應(yīng)力80MPa，超出材料的強度極限，可能導(dǎo)致產(chǎn)品在此處發(fā)生破裂；每米長度的吸能量達192.43kN?m，為標(biāo)準(zhǔn)值的3.5倍，吸能效果非常顯著；但每米長度的最大反力達1708.41kN，為標(biāo)準(zhǔn)值6倍，反力過大易造成船舶和碼頭損壞。

（3）該護舷產(chǎn)品的設(shè)計存在一定缺陷，有必要按照標(biāo)準(zhǔn)要求，對其截面結(jié)構(gòu)進一步優(yōu)化。