任意几何构型浮体之静水剖面负荷计算方法.pdfVIP

任意幾何構型浮體之靜水剖面負荷計算方法 林宗岳* 陳冠臻* 施囿丞** *財團法人中國驗船中心 研究處研發組 **台灣大學 工程科學及海洋工程學系 摘要 在船舶與海洋工程領域的應用裡，各種幾何構型因任務執行與安全保障需要計算在靜水 中的穩性與強度。傳統用辛普森積分方法計算，沿著某一固定方向切片積分，得到其縱向剖 面負荷，經常因為數值方法誤差而無法保持艏和艉為自由端；對於構件方向任意安排之結構 物，此縱向積分缺少通用性。本研究針對組成浮體系統的所有單元構件的外殼曲面利用眾多 的三角形平面小板加以離散，利用解析的方式精確地計算沒入水中的所有三角形平面小板上 承受的靜水壓力；搭配給定之重量分佈，可算得剖面所受之負荷，亦即剪切、壓縮、彎曲、 扭轉等四個方向負荷，作為後續結構分析使用或縱向強度校核。本文使用箱形船、單體船、 以及雙胴體船對所建立方法進行展示，於給定設計工況下之計算結果與與商用軟體之結果進 行比較，驗證本理論方法的正確性。   關鍵詞：船體彎矩及剪力，剖面負荷，靜浮力，三角化網格     引言 船舶安全的兩大支柱為穩度與強度，於國際防止船舶污染公約 (MARPOL)和船級協會規章 裡皆有嚴格的標準。穩度校核之標的指標為扶正力臂，船舶必需在各裝載條件以及完整/破 損條件下達到要求 [1] ，以保證船體不致翻覆。強度校核之標的指標為剖面容許應力，其計 算方法為剖面負荷除以剖面模數 [2] 。以縱向強度來說，舯部0.4 倍船長內之剖面模數由船 體縱向結構材構成，包含：船殼外板、縱向連續之外板加強材、甲板、內底板、縱隔艙壁等； 當舯剖面設計完成時即決定剖面模數。容許應力則由材料及其安全係數所決定，以鋼材來說 為 175 N/mm2 ，安全係數視鋼材強度而定。是故決定船體縱向強度是否足夠之變因僅剩剖面 負荷，橫剖面彎矩和剪切通常拆解為靜水彎矩/剪切，以及波浪彎矩/剪切，後者由經驗或規 範公式估算，前者依照船殼外形與重量分佈直接計算。 涉及前述穩度與強度直接計算者均為靜水條件，利用船舶外型瘦長形特徵來進行分析， 可將船體分割為多段橫剖面，計算每個二維剖面之靜水特性，繪製成龐琴曲線 [3] ，再組合 回三維船型。此類 2.5 維的計算方式可靠且迅速，已廣泛被採納；然而對於形狀特殊者，或 剖面方向與瘦長型特徵不平行者，此方法則無法適用。再者，由二維特徵組合成三維特徵， 建基於等剖面特性或線性變化，在積分時使用梯形法或辛普森法 [4]造成數值誤差，尤以剖 面特徵變化劇烈處，如艏艉部或雙體船之縱剖面誤差甚巨。本研究嘗試直接以三維空間進行 靜水性能的解析解推導，繞過數值積分法，計算船殼面上之靜水壓力分佈，如下圖 1 ，一來 可用於評估穩度，二來可計算剖面負荷。兩者的理論基礎一致，是為一體的計算方法。 數值方法 定義船體座標系統：+X 為艉向艏，+Y 為左舷，+Z 為基線至甲板。座標原點位於艉垂標 基線上。依照圖 1所描述之三個輸入：船殼曲面離散為三角形非構化面網格，要求為簡單封 閉之拓撲。第二，剖面定義為橫向 (YZ平面)或縱向 (XZ平面) ，每個剖面給定其面上之中性 軸座標位置。第三，重量分佈包括質量、橫向/縱向之跨距和重心位置、重心垂向高度等。 為了達到靜力平衡，首先需算出給定重量構成之總重量與重心位置，並與浮力和浮心平衡。 完成後得到平衡姿勢下的靜壓力分佈；搭配前述重量分佈，經剖面切割後得到剖面負荷。以 下簡述靜浮力模型、剖面定義、以及重量分佈模型。 1 ，靜浮力模型 一浮體於靜水中穩定，依據阿基米德原理，其排開水的重量為浮力，排開水的體積中心 為浮心；重量與浮力大小一樣，方向相反，且作用於同一鉛垂線上。此時水平面與封閉浮體 幾何相交於水線面，而本方法將浮體幾何離散為三角形網格，如圖 2(左)所示。故三角形與 平面之交集狀況如圖 2(中)共五種：全沒水、二頂點沒水、一頂點沒水、全出水、以及與平 面重合。部分沒水狀況定義變數 h ，表示(一頂點)沒水或 (二頂點)出水深度。二頂點沒水之 四邊形可由對角頂