賀遠(yuǎn)松 伍友軍 朱 凌 蔡 偉

(中國船舶及海洋工程設(shè)計研究院1) 上海 200011) (武漢理工大學(xué)高性能船舶技術(shù)教育部重點實驗室2) 武漢 430063) (武漢理工大學(xué)交通學(xué)院3) 武漢 430063)

0 引 言

現(xiàn)代艦船甲板受到輪印載荷作用的情況越來越多，如艦船甲板遭受直升機輪印載荷作用，車輛甲板遭受重載車輛輪印載荷作用等.在一些極端輪印載荷作用下甲板結(jié)構(gòu)會發(fā)生彈塑性變形，這將會影響甲板結(jié)構(gòu)的工作性能和安全性能.

Jackson等[1]開展了對輪載作用下甲板板架結(jié)構(gòu)的塑形響應(yīng)進(jìn)行了系列實驗研究，得到不同板格尺寸下加載力-變形歷程曲線，給出了不同永久變形準(zhǔn)則下板厚計算設(shè)計曲線.Hughes[2]基于剛塑性理論，提出了板格在均布載荷和集中位置均布載荷作用下的半經(jīng)驗公式，并基于不同永久變形準(zhǔn)則給出板厚計算設(shè)計圖譜.Konieczny等[3]研究了局部均布載荷作用下甲板板格的非線動力響應(yīng)，基于載荷與永久變形之間的關(guān)系，提出了板厚設(shè)計公式.彭興宇等[4]基于Hughes公式，根據(jù)許可永久變形準(zhǔn)則提出了輪載作用下甲板板格的板厚設(shè)計圖譜.曾鳴等[5-6]通過參考現(xiàn)有船級社規(guī)范和有限元直接計算結(jié)果，制定出車輛甲板的結(jié)構(gòu)設(shè)計和強度計算標(biāo)準(zhǔn)大綱，此外還針對氣墊船著陸墊與車輛甲板之間構(gòu)成非線性接觸問題，利用有限元計算方法分析不同工況下由接觸導(dǎo)致的甲板各構(gòu)件的最大應(yīng)力和最大位移值.張文強等[7]根據(jù)了相關(guān)船級社規(guī)范的規(guī)定，提出了多輪載荷等效為單輪印的方法.王維舟等[8]依據(jù)各船級社規(guī)范，分析了設(shè)計板厚下板格在輪印載荷作用下的應(yīng)力水平及離散程度，并分析了車輛甲板考慮塑性承載能力的程度.王智慧等[9]開展了加筋板輪印載荷分配的實驗和數(shù)值研究，分析了幾種典型輪印位置、輪印大小和輪印橡膠厚度下的輪印載荷分布特性，并且研究了波紋夾層板的輪印載荷分配及極限承載能力.劉聰?shù)萚10]通過在橡膠塊上加壓的方式模擬輪印載荷，測量橡膠塊與加筋板結(jié)構(gòu)之間的接觸壓力，揭示輪印載荷分布特性.張玉龍等[11]針對均布輪印載荷下的波紋型夾層板格結(jié)構(gòu)，建立了應(yīng)力分析計算方法，分析了輪印載荷下夾層板格總的彎曲特性.康杰豪等[12]針對多輪印載荷工況作用于船舶多跨梁結(jié)構(gòu)時，提出一種將遺傳算法與有限元方法相結(jié)合求解任意多跨梁上有多種輪印載荷作用時最危險工況的方法.朱廣紀(jì)等[13]提出一種輪印載荷下多跨梁裝載方案分級方法.Shi等[14]開展了一系列剛性壓頭準(zhǔn)靜態(tài)加卸載下甲板板格彈塑性響應(yīng)的試驗和數(shù)值研究，分析了不同板格尺寸下加載力-變形曲線和板格塑性變形模態(tài).基于前期的甲板板格彈塑性響應(yīng)研究工作，Zhu等[15]總結(jié)了輪印載荷作用下甲板結(jié)構(gòu)塑性設(shè)計準(zhǔn)則和原理，并提出了甲板結(jié)構(gòu)塑性設(shè)計方法.

由于輪印作用在甲板結(jié)構(gòu)位置不同，甲板板格產(chǎn)生彈塑性變形也會有所不同，文中針對典型輪印載荷作用工況，開展了輪印載荷作用下甲板板架結(jié)構(gòu)彈塑性響應(yīng)數(shù)值仿真研究，討論了不同工況下相同輪印載荷大小下甲板板架結(jié)構(gòu)彈塑性響應(yīng)，并基于許可永久變形準(zhǔn)則，提出了特定輪印尺寸和板架尺寸下甲板板厚設(shè)計曲線.

1 計算模型及工況

1.1 模型尺寸

根據(jù)實船的甲板結(jié)構(gòu)，選取3×4局部甲板板架模型為計算對象，見圖1.圖2為直升機的起降輪組尺寸.將中心處的板格標(biāo)記為Plate-1和Plate-2，板格Plate-1和Plate-2處四周縱向骨材角鋼依次標(biāo)記為L1、L2和L3.

圖1 3×4局部甲板板架模型(單位：mm)

圖2 車輪輪組尺寸(單位：mm)

1.2 計算工況

選取四種危險工況為典型位置工況進(jìn)行分析，見圖3，主要分為兩種情形：①車輪縱向布置；②車輪橫向布置.

圖3 計算工況(單位：mm)

1.3 有限元計算模型

采用商業(yè)軟件ABAQUS對以上幾種工況進(jìn)行有限元計算，在有限元建模時，采用shell單元模擬板材、T型材及等效球扁鋼的腹板，采用beam單元模擬等效角鋼的球頭；板架模型采用雙線性材料模型，屈服強度590 MPa，彈性模量E為206 GPa，泊松比為0.3，切變模量Gt=E/1 000.在進(jìn)行網(wǎng)格劃分時，總體單元劃分為規(guī)則四面體單元，輪印邊界部分采用不規(guī)則的四面體單元及三角形單元，并輪印接觸區(qū)域進(jìn)行網(wǎng)格加密處理，模型整體部分網(wǎng)格尺寸為50 mm，網(wǎng)格加密部分網(wǎng)格尺寸為20 mm，其他網(wǎng)格過度部分為20～50 mm；加載時，采用局部均布載荷來描述輪印載荷，單個輪印設(shè)計載荷p為2.2 MPa；邊界條件為板架四周剛性固定.四種工況下有限元計算模型見圖4.

圖4 有限元計算模型

2 四種工況計算結(jié)果及分析

對上述四種工況下的有限元模型進(jìn)行計算，單個輪印載荷p為2.2 MPa，在計算中設(shè)置加載步和卸載步，對板格進(jìn)行加卸載.下面選取板厚為14 mm為分析對象，對比四種工況下的板架的典型位置處應(yīng)力、應(yīng)變狀態(tài)，根據(jù)等效塑性應(yīng)變云圖(PEEQ)來判斷模型是否進(jìn)入塑性階段，當(dāng)PEEQ>0時，即模型發(fā)生屈服，此外，并以Mises應(yīng)力來描述典型位置處的應(yīng)力狀態(tài).圖5為四種工況下的輪印載荷加載到2.2 MPa時甲板板架模型的Mises應(yīng)力云圖，由圖5可知，輪印載荷加載到2.2 MPa時，工況1和工況2下甲板板架上Mises應(yīng)力主要集中在輪印載荷施加區(qū)域及板格Plate-1、Plate-2邊緣與中間角鋼L2相交處，而工況3和工況4甲板板架上Mises應(yīng)力主要集中在輪印載荷施加區(qū)域及板格Plate-2邊緣與角鋼相交處.表1為四種工況下輪印載荷完全卸載后甲板板架上應(yīng)力、應(yīng)變結(jié)果，由表1可知，輪印載荷卸載到0時，工況1和工況2下甲板板架仍然有較大的殘余應(yīng)力，且已經(jīng)進(jìn)入塑性，然而，工況3和工況4下甲板板架幾乎處于彈性變形階段，沒有出現(xiàn)塑性應(yīng)變情況.因此，相同輪印載荷大小時，工況1和工況2下在卸載后板架上仍然有較大的殘余應(yīng)力和塑性應(yīng)變，比工況3和工況4下更加危險.

圖5 輪印載荷加載到2.2 MPa時板架的Mises應(yīng)力云圖

表1 四種工況下輪印載荷完全卸載后甲板板架上應(yīng)力、應(yīng)變結(jié)果對比

為了進(jìn)一步分析板架模型在加缷載過程中的載荷-變形曲線變化關(guān)系，以甲板板厚為14和16 mm為例，單個輪印載荷p為2.2 MPa，對比四種工況下的板架模型變形隨載荷的變化關(guān)系，見圖6.由圖6可知，工況1和工況2下的加載到2.2 MPa時板架模型的最大變形小于工況3和工況4下的最大變形，卸載到0后工況1和工況2下板架模型的永久變形大于工況3和工況4下的永久變形.圖7為不同板厚下板架模型的最大和永久變形值，由圖7可知，相同條件下，甲板板架板厚越大，加卸載過程中最大變形值和永久變形值越小.此外，不同板厚下工況1和工況2下即輪印縱向布置時的最大變形值比較接近，工況3和工況4下即輪印橫向布置的永久變形值比較接近.在四種工況中，工況1下板架模型的最大變形值最小但是永久變形值最大.由于在進(jìn)行車輛甲板結(jié)構(gòu)設(shè)計中，常常關(guān)注的是輪印載荷對甲板產(chǎn)生的塑性變形值，由表2和圖7b)可知，工況1下板架模型的永久變形值最大，故此種情況是在甲板設(shè)計過程需要考慮的最危險工況.

圖6 四種工況下板架模型載荷-變形曲線

圖7 不同板厚下板架模型的最大和永久變形值

3 工況1下甲板板架板厚設(shè)計圖譜實例

由于工況1(輪印關(guān)于縱向骨材對稱縱向布置情形)下甲板板架處于最危險狀態(tài)，在進(jìn)行甲板結(jié)構(gòu)設(shè)計時常?？紤]最危險載荷工況，以使甲板結(jié)構(gòu)偏向安全設(shè)計.因此，主要針對工況1下最危險情形，運用彈塑性設(shè)計方法，給出甲板板厚設(shè)計曲線實例，以便為甲板設(shè)計提供參考.下面板架模型尺寸仍然為工況1中的計算模型，單個輪印尺寸仍然為460 mm×270 mm，甲板板厚選取為12，13，14，15，16，17和18 mm，輪印載荷大小選取范圍為1～10 MPa.為便于設(shè)計及進(jìn)行多方案比較，減少設(shè)計變量數(shù)目，定義以下量綱一的量參數(shù)：

(1)

(2)

(3)

式中：s為甲板縱骨間距，m；t為甲板板厚，mm；σs為甲板模型材料的屈服應(yīng)力，MPa；E為楊氏模量，MPa;Wp為甲板模型中心的最大永久變形，mm；F=2pab為甲板板架中承受的載荷,103kN,其中:F為雙輪印載荷之和；p為局部均布壓強，MPa；a和b分別為輪印載荷的長度和寬度，m.

針對上述7種板厚下不同載荷下的有限元計算結(jié)果值，并對載荷進(jìn)行無量綱化.圖8為不同板厚下的板架模型經(jīng)過不同大小的輪印載荷P加卸載后的永久變形值WP.圖9為不同柔度系數(shù)Cb下的在無量綱載荷參數(shù)下QP板架模型的WP/t值.由圖9可知，當(dāng)板格的柔度系數(shù)Cb值越大，相同載荷參數(shù)QP下的永久變形值越大，或者相同永久變形下所對應(yīng)的載荷值越小.

圖8 不同板厚下載荷-永久變形曲線

圖9 不同Cb下QP-WP/t曲線

為了得到板厚設(shè)計圖譜，需要確定甲板板架的許可永久變形參數(shù)，通常來說，許可的永久變形值WP為s/150～s/50，通過式(3)變換計算可得，許可永久變形參數(shù)CW為0.125～0.37.為此，在本報告中，選取CW=0.1和CW=0.2為許可永久變形設(shè)計衡準(zhǔn)，由式(1)和式(3)可知，WP/t=CW·Cb，由此可以得到不同柔度系數(shù)下的許可永久變形值WP/t，見表2，然后在圖9中進(jìn)行插值得到許可永久變形值WP/t對應(yīng)的無量綱載荷值QP，從而可以得到兩種設(shè)計衡準(zhǔn)下(CW=0.1和CW=0.2)的無量綱載荷QP～Cb板厚設(shè)計圖譜，見圖10.

表2 不同柔度系數(shù)下的許可永久變形值WP/t

圖10 兩種永久變形衡準(zhǔn)下QP-Cb板厚設(shè)計曲線

假設(shè)有一段由高強度船用鋼組成的甲板分段，其中材料屈服強度為590 MPa，彈性模量為206 000 MPa，甲板縱骨間距s=700 mm，甲板板格長寬比l/s=3.428，單個輪印設(shè)計載荷p=2.4 MPa，單個輪印尺寸仍然為460 mm×270 mm，根據(jù)式(2)可以計算得無量綱載荷參數(shù)QP=0.72，若以CW=0.1為所選取的設(shè)計衡準(zhǔn)，可得Cb=3.5，然后根據(jù)式(1)，可得甲板板厚初步設(shè)計值td=10.7 mm，故甲板板厚設(shè)計值不得低于10.7 mm.

4 結(jié) 論

1) 在相同輪印載荷、輪印尺寸和板架尺寸下，工況1和工況2下在卸載后板架上仍然有較大的殘余應(yīng)力和塑性應(yīng)變，板架模型已經(jīng)進(jìn)入了塑性狀態(tài)，其中，卸載后殘余應(yīng)力和塑性應(yīng)變主要集中在板格Plate-1和Plate-2與中間角鋼L2相交的板格邊緣處，而工況3和工況4下板格上有較小的殘余應(yīng)力值且塑性應(yīng)變值幾乎為0，板架模型可近似認(rèn)為還處于彈性范圍內(nèi).

2) 工況1和工況2下的加載到2.2 MPa時板架模型的最大變形小于工況3和工況4下的最大變形，卸載到0后工況1和工況2下板架模型的永久變形大于工況3和工況4下的永久變形.在進(jìn)行甲板結(jié)構(gòu)設(shè)計中，常常關(guān)注的是甲板產(chǎn)生的塑性變形值，而工況1下板架模型的永久變形值最大，故此種情況是在甲板設(shè)計過程需要考慮的最危險工況.

3) 針對工況1下最危險情形，運用彈塑性設(shè)計方法，基于永久變形準(zhǔn)則給出特定輪印尺寸和板架尺寸下甲板板厚設(shè)計曲線，并給出甲板板厚彈塑性設(shè)計方法實例.