宋 陽， 馬 駿

(大連理工大學 船舶工程學院， 遼寧 大連 116024)

LPG船大型分段運輸動力學分析

宋 陽， 馬 駿

(大連理工大學 船舶工程學院， 遼寧 大連 116024)

以83 000 m3超大型全冷式液化氣船雙層底分段作為研究對象，根據(jù)實際模型圖紙數(shù)據(jù)建立分段運輸模型，主要分析其在不同運輸工況下動力學特性。在設置運輸環(huán)境時，主要分析路面不平整和緊急剎車兩種運輸工況對分段結構可靠性的影響。通過分析出的大型分段運輸最佳工況，為提升船廠建造LPG船等高端船舶的效率提供參考。

LPG船；大型分段；運輸動力學；非線性有限元法；接觸分析

0 引言

船體分段運輸動力學是研究船體結構可靠性的學科，指在船舶建造過程中，由平板車運送船體分段過程中，對船體分段的強度、變形狀況及穩(wěn)定性的研究。運輸動力學研究起始于20世紀50年代，現(xiàn)階段主要集中在汽車業(yè)，汽車運輸動力學方面有代表性的著作是德國MITSCHKE和WALCENTOWITZ[1]編著的汽車動力學。美國GILLESPIE[2]教授也全面介紹了車輛動力學的基本工程原理及汽車性能的分析方法。國內有代表性的著作是張洪欣[3]教授的汽車系統(tǒng)動力學和喻凡、林逸[4]教授所著的汽車系統(tǒng)動力學。合肥工業(yè)大學的史沛瑤[5]基于Adams/CAR軟件，根據(jù)路面激勵的分類建立了整車運輸包裝系統(tǒng)動力學模型。在對LPG船分段的研究中，大多采用隱式有限元分析法，通過建立LPG 船艙段的三維有限元模型，對其進行動力學分析[6]。

雖然大型貨物運輸動力學的相關研究已經(jīng)在國內外取得了一些成果，但是在這些研究成果中多數(shù)是偏宏觀和定性的陳述。針對大型貨物運輸可靠性而開發(fā)的許多系統(tǒng)，也多是面向大型貨物運輸前期工作中的路線規(guī)劃。對于實際運輸過程中影響大型貨物運輸可靠性的相關因素卻鮮有涉及，而為數(shù)不多的涉及安全穩(wěn)定性的研究也只是處于理論推導階段。實際上，大型貨物運輸過程中存在著多種安全問題，例如：貨車上坡時其牽引力和產(chǎn)生的摩擦力能否滿足不發(fā)生相對移動的條件；下坡時制動力多大才能保障可靠性；路面不平整度對運輸過程的影響；等。本文選取LPG船分段建立計算模型，針對運輸過程中路面不平整激勵、不同車速緊急剎車激勵等特殊運輸工況以及如何施加載荷做詳細的計算與分析。

1 數(shù)值計算模型

1.1 數(shù)值方法

采用非線性有限元法求解接觸載荷，在求解非線性問題的有限元方法中，Lagrange描述增量法是較常用的方法之一。在使用Lagrange方法中，節(jié)點作為形成網(wǎng)格的基礎成員分布在研究對象上，各個節(jié)點之間通過有序連接形成特定單元，網(wǎng)格則由這些特定單元組合而成。在變形分析中，節(jié)點、單元都隨著結構材料一起移動，單元質量恒定。

1.2 模型的建立

LPG船大型分段運輸模型主要包含：LPG船分段、墩木和平板車3個部分。其中，LPG船分段選取的是83 000 m3全冷式LPG船的左舷典型雙層底大型分段。該分段模型重176.8 t，縱向尺寸為16.9 m，橫向尺寸為18.3 m，垂向尺寸為6.75 m。墩木邊長1.2 m，寬0.8 m，高0.5 m，共15個。墩木作為分段與平板車之間的過渡支撐，位于旁桁材與框架肋板交界處下方，如圖1所示。

圖1 LPG船大型分段、墩木和簡化平板車模型

在大型分段建模過程中，通常將橫截面較小的梁結構使用線單元建立，通過定義截面屬性的方式來模擬梁的效果，但是這樣就會影響模型中與其相近肘板等結構的尺寸。為了讓模型分析結果更加準確，在建立船體分段模型的過程中，全部使用面單元Shell 163建立，平板車的承載面也由Shell 163單元建立，考慮平板車自身特性，將此模型中的平板車視為剛體。墩木模型及平板車承載面下方簡化輪胎模型所用單元為Solid 164單元[7]。

2 多種運輸工況仿真計算

2.1 分段運輸中基礎參數(shù)設置

在設置分段運輸模型的基礎參數(shù)時，重力加速度為9.8 m/s2，墩木與平板車之間的最大靜摩擦系數(shù)為0.4，滑動動摩擦系數(shù)為0.35。圖2為平板車靜止情況下的計算結果，以供比較。

圖2 平板車靜止時分段模型3個方向的時程曲線

由圖2a可以看出，水平方向(X向和Y向)兩個加速度保持為零，而豎直方向(Z向)的加速度由于初始階段受重力作用影響，墩木和分段受到?jīng)_擊而發(fā)生振蕩。

2.2 分段運輸中緊急剎車工況

特大型平板車的剎車時間通常約1 s，緊急剎車時間約為0.4 s，最大行駛速度為5 km/h。由限制滑移的條件可知，緊急剎車時加速度a<3.92 m/s2?？紤]LPG船分段的特殊性，設定4種工況分別對模型進行分析。緊急剎車工況中約束平板車簡化輪胎模型在豎直方向上的位移，平板車運行總時間為10 s。4種工況如表1所示。

表1 緊急剎車激勵載荷

2.3 分段運輸中承受路面不平整激勵工況

假定平板車運輸分段過程中，在各種路面不平整激勵載荷作用時行駛速度均為v=0.8 m/s，平板車運行總時間均為10 s，根據(jù)路面不平整激勵的級別，設定了4種激勵載荷，詳細數(shù)值如表2所示。

表2 路面不平整激勵載荷

3 結果分析

截取分析過程中的前4 s結果，如圖3～圖6所示。

3.1 緊急剎車激勵對分段運輸可靠性的影響

圖3為分段模型在3個方向上加速度的變化情況以及靜止狀況下分段3個方向上的加速度比較。從中可以看出，在緊急剎車激勵工況下分段模型只在Y方向上的加速度有較大變化，可知緊急剎車激勵對Y方向的動力學特性影響較大。因此，本文以剎車激勵對Y方向上動力學特性的影響為例，分析分段運輸過程中的動力學特性。

圖4為不同行駛速度下平板車與分段的時程曲線。從中可以看出，平板車和分段在1.8 s(剎車)～2.2 s，速度曲線保持不變，這表明兩者受到相同的加速度作用，速度減小趨勢相同，平板車與分段之間沒有產(chǎn)生滑移，此過程為平穩(wěn)剎車。

圖3 不同行駛速度分段模型加速度時程曲線

圖4 不同行駛速度時平板車與分段速度時程曲線

圖5 不同車速緊急剎車激勵對分段行駛方向最大應力及應變影響

除了宏觀地分析分段運輸過程中的穩(wěn)定性之外，還需要定性地分析分段模型在運輸過程中的受力及變形情況。表3列出了4種緊急剎車工況下分段運行方向的最大應力與應變值。

表3 各工況結果數(shù)據(jù)統(tǒng)計

圖6 不同路面激勵下分段模型加速度時程曲線

根據(jù)表3的數(shù)據(jù)繪制變化曲線，得到不同車速下緊急剎車激勵對分段行駛方向最大應力及最大應變的影響，如圖5所示。

從圖5中可以看出，在緊急剎車時間同為0.4 s的條件下，隨著平板車行駛速度的增加，分段行駛方向最大應力和最大應變均不斷增加。由最大應力曲線可以看出，在考慮緊急剎車時間為0.4 s時，此分段運輸過程應控制運行最大速度為0.8 m/s，以保證分段緊急剎車中最大應力不超出屈服極限235 MPa而產(chǎn)生塑性變形。

3.2 路面不平整激勵對分段運輸可靠性的影響

在ANSYS/LS-DYNA中計算4種路面不平整激勵工況，結果如下。

由圖6中后3個工況的加速度變化情況與工況一分段3個方向上的加速度變化情況比較可知，分段模型在Z軸方向上的加速度有較大變化，故路面不平整激勵工況對Z軸方向上的動力學特性影響較大。因此，本文以路面不平整激勵對Z軸方向上動力學特性的影響為例，分析分段運輸過程中的動力學特性。由計算結果列出4種工況下豎直方向最大應力值及最大應變值，如表4所示。

表4 4種工況結果數(shù)據(jù)統(tǒng)計

[][]

根據(jù)表4的數(shù)據(jù)，在Excel軟件中繪制變化曲線，得出不同路面不平整激勵對分段豎直方向最大應力及最大應變的影響，如圖7所示。

圖7 不同路面激勵對分段豎直方向最大應力及最大應變影響

從圖7中可以看到，隨著豎直方向路面不平整激勵的增加，分段在豎直方向上的最大應力在不斷增加，在豎直方向上的最大應變也在不斷增加。由于制造分段所用的鋼材屈服極限為235 MPa，假設平板車行駛速度保持為v=0.8 m/s，則當路面豎直方向激勵值大于22.5 mm時，分段結構強度容易遭到破壞，行駛過程應小心避免。

4 結語

以LPG大型分段為研究對象，采用有限元法對其進行動力學分析，得到如下結論：

(1) 為保證緊急剎車時LPG船分段的可靠性，平板車行駛速度應限制在0.8 m/s以內。

(2) 當平板車以0.8 m/s速度行駛時，大于22.5 mm的豎直方向激勵易使分段遭到破壞，應謹慎避免。

[1] (德)MITSCHKE M, WALLENTOWITZ H.汽車動力學[M].4版.北京：清華大學出版社,2009.

[2] (美)GILLESPIE T D.車輛動力學基礎[M].趙六奇,金達鋒，譯.北京：清華大學出版社,2006.

[3] 張洪欣.汽車系統(tǒng)動力學[M].上海：同濟大學出版社,1996.

[4] 喻凡,林逸.汽車系統(tǒng)動力學[M].北京：機械工業(yè)出版社,2005.

[5] 史沛瑤.重型貨車運輸過程動力學建模分析及貨物可靠性影響因素研究[D].合肥：合肥工業(yè)大學,2013.

[6] 梅明.全壓式LPG船結構強度直接計算研究[D].武漢：武漢理工大學,2013.

[7] 張紅松,胡仁喜,康士延,等.ANSYS 14.5/LS-DYNA非線性有限元分析實例指導教程[M].3版.北京：機械工業(yè)出版社,2013.

Dynamics Analysis of Large Block of LPG Ship in Transport

SONG Yang, MA Jun

(Department of Naval Architecture, Dalian University of Technology, Dalian 116024, Liaoning， China)

Taken the block size of 83 000 m3double bottom LPG ship as the research target, the establishment of transportation model block is based on the data of actual model drawings. The dynamic characteristics of the block in different transportation conditions are mainly analyzed. When the transportation environment is set up, the uneven pavement excitation and the emergency braking are taken into account. The best transportation environment of large block size can provide references for the efficiency improvement of advanced shipbuilding such as LPG.

LPG ship; large block; transport dynamics; nonlinear finite element method; contact analysis

宋 陽(1990-)，男，碩士研究生，研究方向為船舶與海洋結構物振動及強度分析研究。

U673

A