三維有限元素熱機耦合框架下 316L 不銹鋼管冷彎與深冷殘留應力演化及水下與船岸系統全壽命疲勞評估 (Evolution of Residual Stress in Cold Bending and Cryogenic Treatment of 316L Stainless Steel Pipes under a 3D Thermo-Mechanical Coupled Finite Element Framework and Full Life-Cycle Fatigue Assessment for Subsea and Ship-to-Shore Systems)

一、摘要與緒論

1.1 研究背景：全球能源轉型與深冷傳輸技術

在全球能源轉型的宏觀背景下，液化天然氣（LNG）與液態氫（LH2）的跨洋運輸與水下傳輸系統，已成為支撐現代潔淨能源供應鏈的關鍵基礎設施。為了實現高效能的體積壓縮，LNG 必須儲存與輸送於 -162°C（111 K）的深冷狀態，而液態氫更需低達 -253°C（20 K）。在這些極端工程應用中，包含海底管線（Subsea Pipelines）、浮式液化天然氣船（FLNG）、常規 LNG 運輸船（LNGC）以及碼頭裝卸臂（Marine Loading Arms），結構材料的選擇至關重要¹。

316L 沃斯田鐵不銹鋼（Austenitic Stainless Steel）憑藉其面心立方（FCC）晶體結構在極低溫下不發生延脆轉移（Ductile-to-Brittle Transition）的特性，兼具極高的斷裂韌性與抗腐蝕能力，被廣泛指定為深冷壓力容器與傳輸管線的標準結構材料³。國際海事組織（IMO）的《國際散裝液化氣體船舶構造與設備章程》（IGC Code）與相關指導原則，亦嚴格規範了此類船舶與離岸設施在設計與建造時，必須將深冷環境與動態載荷導致的材料疲勞、微觀相變與斷裂風險降至最低⁵。

1.2 工程挑戰：多物理場熱機耦合歷史

然而，LNG 船岸系統與海底管線在整個全壽命週期中，經歷了極為嚴苛且複雜的多物理場熱機耦合（Thermo-mechanical Coupling）歷史。首先，在船體內部狹窄空間、碼頭棧橋或海床佈局中，管線必須經過感應冷彎成形（Cold Induction Bending，如 3D 或 5D 彎曲半徑）以適應複雜的幾何限制。此一強烈塑性變形過程不可避免地會引入顯著的壁厚減薄、截面橢圓化以及極高量級的殘留應力（Residual Stress）⁸。其次，當系統進行 LNG 裝卸作業或深海輸運，並通入低達 -162°C 的流體時，管壁材料將經歷劇烈的熱收縮與頻繁的熱循環。在極低溫與高殘留應力疊加的臨界狀態下，316L 內部微觀組織會發生強烈的應變誘發馬氏體相變（Strain-Induced Martensitic Transformation, SIMT），從而徹底重塑管線的宏觀與微觀應力場¹⁰。最後，船舶航行時的六自由度運動慣性力、碼頭裝卸臂的機械位移與旋轉扭矩，以及水下懸空段承受的波浪與海流交變載荷，疊加在深冷熱應力與冷彎殘留應力之上，極易引發材料的塑性棘輪效應（Ratcheting）與多軸疲勞損傷¹。

1.3 研究目的與方法

為彌補現有研究往往將「製程殘留應力」、「深冷微觀相變動力學」與「船/岸/海動態疲勞」孤立探討的重大缺口，本研究建立了一套嚴謹的三維有限元素分析（FEA）熱機耦合框架。透過引入溫度相依的 Chaboche 混合硬化模型與 Olson-Cohen 相變動力學模型，精確捕捉 316L 鋼管在冷彎及深冷狀態下的應力集中演化機制。最終，結合船舶運動學、碼頭管線受力特徵與 Airy 波浪理論，並嚴格基於 ASME B31.3、ASME BPVC Section VIII Division 2 以及 DNVGL-RP-C203 規範，導入 Smith-Watson-Topper (SWT) 參數進行全壽命週期的疲勞評估。本研究旨在為新一代 LNG 船岸與水下管線的結構完整性設計提供堅實的學術基礎與工程指導。

二、 316L 不銹鋼之深冷力學行為與相變動力學

要精確模擬 316L 不銹鋼在複雜熱機載荷下的力學響應，首要任務是深入探討其在深冷環境下的微觀變形機制。材料的層錯能（Stacking Fault Energy, SFE）是決定其塑性變形機制的關鍵熱力學參數，並且強烈依賴於環境溫度與化學成分¹⁰。

2.1 層錯能演化與微觀變形機制之轉變

在室溫（約 293 K）下，316L 的層錯能相對較高，材料的塑性變形主要由完美的位錯滑移（Dislocation Glide）與交滑移（Cross-slip）主導。此時材料呈現穩定的沃斯田鐵（Austenite）特徵，具備優異的延展性。然而，當環境溫度驟降至深冷區域（如 LNG 的 -162°C 或液氫的 -253°C），316L 的層錯能會顯著降低至臨界極限值以下¹²。層錯能的急劇下降使得全位錯極易分解為 Shockley 不全位錯，並形成寬闊的層錯帶。這導致材料的變形機制由單一的位錯滑移，轉變為孿晶誘發塑性（Twinning-Induced Plasticity, TWIP）與相變誘發塑性（Transformation-Induced Plasticity, TRIP）的複雜混合機制¹⁴。

在深冷變形過程中，316L 會經歷 γ (FCC) → ε (HCP) →α’ (BCC/BCT) 的多階段相變路徑¹⁵。其中，具有六方最密堆積（HCP）結構的 ε-馬氏體通常為過渡相，優先在交叉的剪切帶成核；隨後轉變為體心立方/體心四方（BCC/BCT）結構的α’-馬氏體。α’-馬氏體的硬度極高，其大量生成雖顯著提升了材料的宏觀降伏強度與極限抗拉強度（例如在深冷下其極限抗拉強度可飆升至 1200 MPa 以上），但也伴隨著延展性的局部喪失³。更為關鍵的是，BCC/BCT 結構的比容大於 FCC 結構，此相變過程會產生約 2% 至 3% 的局部體積膨脹（Volumetric Dilatation）¹⁴。這種體積膨脹在受幾何拘束的船體或水下管線結構中，會產生顯著的內部微觀壓應力，從而深刻影響宏觀殘留應力的演化。

2.2 Olson-Cohen 相變動力學模型之導入

為了在 FEA 框架中精確量化 α’-馬氏體隨應變與溫度的動態演化過程，本研究引入了具備明確物理意義的 Olson-Cohen 模型。該模型假設α’-馬氏體優先成核於兩個或多個相交的剪切帶交匯處。馬氏體的體積分數 f^α’ 可表示為等效塑性應變 ε ̅^p 的非線性 S 型（Sigmoidal）函數¹⁴：

f^α’=1-exp{-β[1-exp(-αε ̅^p ) ]ⁿ }

在此方程式中：

α 參數控制著剪切帶隨塑性應變的形成速率。由於層錯能隨溫度下降而降低，剪切帶在深冷下極易形成，因此低溫會顯著放大 α 值¹⁸。
β 參數代表剪切帶交匯處成功轉變為 α’-馬氏體胚核的機率。低溫提供了極大的化學驅動力（Chemical Driving Force），使得 β 值在深冷環境下急劇上升¹⁸。
n指數與微觀組織的幾何特徵相關，對於 316L 通常取值為 2 至 3 ¹⁴。透過自訂的材料子程序（UMAT），FEA 模型能夠在每個積分點動態計算當前溫度與應變下的相變體積分數，並將相伴隨的體積膨脹轉化為相變應變張量，實現真實的熱-機-相變多場耦合分析。

三、三維循環塑性與 Chaboche 混合硬化本構模型

LNG 船岸管線在冷彎成形過程（極大單調塑性應變）以及頻繁的裝卸貨熱循環與波浪載荷（中小幅值交變塑性應變）中，經歷了高度非比例的三維加卸載。傳統的各向同性硬化（Isotropic Hardening）或線性隨動硬化（Kinematic Hardening）模型無法精確描述 316L 不銹鋼在此複雜加載歷史下表現出的包辛格效應（Bauschinger Effect）、瞬態循環硬化/軟化以及棘輪效應²⁰。因此，本研究採用了先進的 Chaboche 混合硬化模型，並將其拓展為溫度與塑性應變相依的形式。

3.1 Chaboche 模型的理論架構

Chaboche 模型的降伏準則基於 von Mises 屈服面理論，引入了背應力張量（Backstress Tensor） X 以表徵屈服面在應力空間中的平移，以及各向同性硬化變量 R 以表徵屈服面的均勻擴張或收縮²¹：

F=√(3/2 (S-X):(S-X) )-σ_Y-R=0

其中，S 代表偏應力張量，σ_Y 為材料的初始降伏應力。為精確擬合 316L 從微小應變到極大塑性應變的寬廣響應範圍，總背應力 X 被分解為多個非線性隨動硬化分量Xi 的疊加（本研究採用 M=3）²¹：

X=∑_i=1^MXi

每一個分量的演化法則遵循 Armstrong-Frederick 的非線性微分方程式：

dX_i=2/3 C_idε ̅^p-γ_iX_idε ̅^p

其中， C_i為初始動態學硬化模數，γ_i 為控制該分量達到飽和速率的非線性動態回復參數²⁰。各向同性硬化 R 的演化則定義為：

R=Q(1-exp(-bε ̅^p ) )

此處 Q 代表各向同性硬化的最大飽和值（若為循環軟化則為負值），b 為飽和速率參數。

3.2 316L 溫變循環塑性參數之標定

316L 不銹鋼的循環響應具有強烈的溫度相依性。在室溫（20°C）下，材料在低週疲勞測試中表現出初始短暫硬化後進入穩定的「循環軟化」特徵（Q 值為負）；然而，當環境溫度降至 -162°C 下，由於 α’-馬氏體相變誘發塑性效應的強烈介入，軟化現象完全消失，轉而表現出持續且劇烈的「循環硬化」（Q 值為極大的正值）²⁰。本研究基於文獻實驗數據，標定了室溫與 -162°C 下的 Chaboche 多物理場參數，詳見表 1。

參數類別	物理意義	室溫 (20°C) 參數值	深冷環境 (-162°C) 參數值
基礎熱力學	楊氏模數 E (GPa)	193	205
	蒲松比 v	0.30	0.28
	熱膨脹係數 α_th(/K)	15.9×10^-6	14.5×10^-6
初始降伏	初始降伏應力σ_Y (MPa)	205	540
隨動硬化	分量 1: C₁(MPa),γ₁	189500, 2950	235000, 3200
	分量 2: C₂ (MPa), γ₂	33500, 350	52000, 480
	分量 3: C₃(MPa), γ₃	4500, 15	8200, 25
等向硬化	最大飽和值 Q (MPa)	-60 (循環軟化)	+380 (強烈循環硬化)
	飽和速率參數 b	5.5	8.2

表 1：316L 不銹鋼在室溫與深冷環境下之多物理場與 Chaboche 循環塑性參數設定²⁰。

四、 3D 與 5D 半徑冷彎成形之殘留應力場與幾何解析

為適應 LNG 船舶狹窄的機艙空間、接收碼頭的複雜棧橋或海床起伏，管線常需透過感應彎曲（Induction Bending）等冷成形工法來改變走向。彎曲半徑 R 通常以管線公稱外徑 D₀ 的倍數表示⁸。本研究以公稱管徑 NPS 8（外徑219.1 mm），厚度等級 Schedule 80（壁厚12.7mm）的 316L 鋼管建立三維高精度六面體網格 FEA 模型。

4.1 幾何變異：壁厚減薄與截面橢圓度分析

彎管成形過程本質上是在管線兩端施加極大彎矩，這會導致外側（Extrados）承受巨大張應力而拉伸減薄，內側（Intrados）承受壓應力而壓縮增厚⁸。依據 ASME B31.3 與 B16.49 標準，成品彎管外側的減薄厚度絕對不得低於直管段的最小設計耐壓厚度，且截面橢圓度（Ovality）必須嚴格控制在 3% 以內⁸。

透過解析公式t_thin=t×(R/(R+D_o/2)) 與 FEA 模擬的交叉驗證，本研究量化了厚度變化率²⁹：

3D 冷彎 (R=657.3 mm)：外側壁厚減薄率高達29%（某些實務極端情況可達 18%），內側增厚率則達 20.00%²⁸。
5D 冷彎 (R=1095.5 mm)：外側壁厚減薄率顯著降低至09%（約在 9% 至 12.5% 之間），幾何畸變大幅減緩²⁸。在 3D 彎管極端塑性變形中，316L 強烈的應變硬化能力會迫使截面中性軸顯著向內側偏移。此偏移加劇了外側的塑性流動區域深度與廣度，使其更逼近法規的安全厚度極限。

4.2 三維自平衡殘留應力場的拓撲特徵

成形結束並卸除外部負載後，管線內部的彈性回彈（Springback）無法完全釋放巨大的塑性應變能，從而在管內形成極高梯度的自平衡殘留應力場⁸。對於 3D 彎管，外側區域由於受到周圍材料的彈性收縮拘束，外表面發生應力反轉，最終穩定於極高的軸向拉伸殘留應力，峰值可高達 350 至 400 MPa³⁰。此數值已遠超 316L 在室溫下的初始降伏強度（205 MPa），意味著局部材料的延展性已被嚴重「耗竭」。相較之下，5D 彎管由於彎曲曲率較緩，等效塑性應變峰值大幅降低約 40%，外側拉伸殘留應力峰值有效降至約 220 MPa。在工程實務上，採用 5D 彎曲能大幅降低管線在後續承受深冷衝擊與海洋波浪振動時的應力腐蝕開裂與疲勞萌生敏感度。

五、 -162°C 深冷縮熱機耦合與裝卸熱循環機制

當 LNG 船舶靠港進行裝卸作業，或水下 LNG 管線啟動通入流體時，整個系統（包含碼頭旋轉接頭、裝卸臂與懸空管段）將經歷從環境溫度（約 20°C）驟降至 -162°C 的劇烈熱衝擊，作業完成後又會逐漸回溫至常溫²。此頻繁且極端的熱循環帶來了嚴峻的熱機耦合挑戰。

5.1 宏觀熱應力與殘留應力的非線性疊加效應

由於管線系統通常受到海底固定塊、法蘭或裝卸臂剛性支架的強烈軸向幾何拘束，阻礙了管線的自由熱收縮。若暫不考慮材料的塑性鬆弛，依據純彈性理論，由ΔT=-182K 溫差所引發的完全拘束熱拉應力可表示為σ_th=E_cryo⋅α_th⋅|ΔT|。將 316L 的深冷參數（E_cryo = 205GPa, 平均 α_th≒ 14.5*10^-6/K）代入，理論熱拉應力高達 541 MPa²⁹。

在嚴格的三維熱機耦合 FEA 框架中，此高達 541 MPa 的宏觀熱應力疊加在 3D 彎管外側原有的 400 MPa 拉伸殘留應力之上。這導致外側區域的等效應力迅速突破 316L 提升後的深冷屈服點（約 540 MPa），進而引發極為顯著的「二次塑性流動」（Secondary Plastic Flow）。這種局部的塑性屈服雖然發揮了應力鬆弛（Stress Relaxation）的作用，避免應力無限攀升，但卻在極低溫下嚴重耗散了該區域材料所剩無幾的殘餘疲勞延展度。

5.2 應變誘發馬氏體相變（SIMT）的雙面刃機制

在二次塑性流動區，極高的應力與 -162°C 的深冷環境共同強烈激活了γ→ α’ 的馬氏體相變機制¹⁰。根據前述 Olson-Cohen 模型的動態預測，熱應力集中區局部的α’-馬氏體體積分數會迅速攀升至 10% 甚至 15%。此相變對應力演化展現出複雜的「雙面刃」機制：

宏觀：相變誘發塑性（TRIP）之應力鬆弛。FCC 轉 BCC/BCT 伴隨的體積膨脹，在受拘束的外側區域產生了等效的「相變壓應力」。這股微觀膨脹力局部抵銷了宏觀的熱收縮應變，成功抑制了宏觀拉伸應力峰值的進一步發散¹⁵。
微觀：相界面錯配應力集中。新生成的馬氏體硬相與未相變的沃斯田鐵軟基體之間，存在巨大的彈性模數與降伏強度差異。在兩者的相界面（Phase Interfaces）處將產生極高的微觀錯配應力。這些微觀應力集中點，在後續的反覆裝卸熱循環與波浪載荷下，極易成為低週疲勞微裂紋（Micro-crack）萌生的溫床¹¹。

六、隨機波浪與船體/碼頭動態載荷評估

除了上述的深冷熱應力，LNG 系統在全壽命週期中必須承受複雜的動態載荷。海底管線需承受直接的海水流體動力，而 LNG 船舶內部管線與碼頭裝卸臂則需承受巨大的慣性力與機械強制位移¹。

6.1 海底管線之 Airy 波浪流體動力學與 Morison 模型

針對敷設於海床或具懸空段（Free Spans）的水下管線，由於其外徑（D₀=0.219m）遠小於典型海洋波浪的波長，流體載荷可直接適用經典的 Morison 方程式進行計算³³：

F(t)=⏟(1/2 ρ_wC_DD_ou(t)|u(t)| )┬曳力 (Drag Force) +⏟(ρ_wC_M(πD_o²)/4 u ̇(t) )┬慣性力 (Inertia Force)

本研究選取台灣高雄外海（假設水深 h=30m）的典型海況數據進行 Airy 線性波浪理論分析²⁹：

日常疲勞海況（顯著波高H_S=1.5m，週期T=6.0s）：計算得出波長L≒1 m，海床最大水質點速度僅u_max=0.0545 m/s。此時與速度平方成正比的曳力極小（0.33 N/m），流體力完全由慣性力（4.41 N/m ）主導，總波浪力為4.41 N/m ²⁹。
極端颱風海況（顯著波高H_S=6.0m，週期T=12.0s）：計算得出波長L≒0 m，海床最大水質點速度急劇增加至u_max=1.23 m/s。非線性的曳力項發生爆發性增長（169.75 N/m），使得總波浪動力載荷激增近 40 倍，達到173.40 N/m ²⁹。這些劇烈的橫向交變載荷疊加在原本就處於臨界狀態的高平均拉應力場之上，將對管線系統造成毀滅性的彎矩衝擊。

6.2 LNG 船舶與碼頭裝卸臂之動態力學特徵

對於船岸系統，動態載荷來源更為複雜：

船舶管線 (LNGC / FLNG)：在極端海況下，船體會經歷劇烈的六自由度運動（橫搖、縱搖、垂盪等）。管線系統不僅需承受內部超低溫流體隨船舶搖晃產生的巨大流體慣性衝擊，還需吸收船體結構大尺度變形所帶來的邊界強制位移¹。
碼頭裝卸臂 (Marine Loading Arms)：根據 OCIMF 規範，裝卸臂（包含多個 3D/5D 彎管與旋轉接頭 Swivel Joints）必須在懸臂狀態下，吸收 LNG 船舶因風浪造成的漂移與吃水變化。這意味著彎管處必須同時承受內部高壓、-162°C 深冷熱應力、接頭旋轉摩擦扭矩，以及來自風載與船舶運動的三維外部動態彎矩²。

七、全壽命週期多軸疲勞壽命評估

LNG 船岸與水下系統所承受的疲勞載荷呈現顯著的雙模態（Bi-modal）特徵：一是頻率極低但應變幅極大的裝卸貨熱循環與極端風暴衝擊（低週疲勞 LCF）；二是頻率高達千萬次以上的日常海浪持續激振與船體引擎振動（高週疲勞 HCF）。

7.1 低週疲勞（LCF）：SWT 模型與 ASME 規範懲罰機制

傳統的 Coffin-Manson 低週疲勞模型僅依賴塑性應變幅進行預測，忽略了冷彎與深冷產生的高「平均拉應力」（Mean Stress）對裂紋張開的加速效應，從而會嚴重高估系統壽命⁴⁰。為此，本研究導入國際公認對平均應力最為敏感的 Smith-Watson-Topper (SWT) 臨界平面參數模型⁴³：

P_SWT=σ_max Δε/2=(σ_‘f )²/E (2N_f)^2b+σ’_f ε’_f(2N_f )^b+c

其中，σ_max包含了平均應力效應，能準確捕捉 3D 彎管外側高拉應力下的疲勞損傷。更為嚴峻的是，針對壓力容器與管線設計指標，本研究引入 ASME BPVC Section VIII Division 2 Part 5 的彈塑性疲勞評估標準。根據 Bree Diagram（布雷圖）評估法⁴⁷，由於 3D 彎管外側區域在深冷熱循環與動態波浪載荷強烈耦合下，明確發生了交變循環塑性，系統被強制劃入 P1 或 P2 區域（Plastic Shakedown / Ratcheting）。在這種狀態下，ASME 規範強制要求在傳統的彈塑性疲勞計算中，人為乘上一個極其嚴苛的疲勞懲罰因子（Fatigue Penalty Factor）K_e,k ⁴⁷。這個K_e,k 因子的引入，將局部交變應力幅成倍放大，極大程度地削弱了 3D 彎管的理論抗疲勞能力，使其幾乎無法滿足任何海上營運的安全性標準。

7.2 高週疲勞（HCF）：DNVGL-RP-C203 規範評估

針對日常海況下高達百萬次以上的連續波浪激振或船體振動，本研究嚴格依據海事工程權威標準 DNVGL-RP-C203 進行 S-N 曲線評估⁵¹。在具腐蝕性的海洋環境中，即便水下管線配備有犧牲陽極等陰極防護（Cathodic Protection）系統，其疲勞曲線仍必須相較於空氣環境進行降級。對於管線的對接銲道（Girth Welds）或高應力集中的彎管區，應適用 DNV 規範中降級的 C1、C2 或 D 等級 S-N 曲線⁵⁴。其基本設計方程為：

logN=loga ̅ -mlog(Δσ(t/t_ref )^k )

分析表明，厚度指數 k（Thickness Effect）的引入會進一步放大厚壁彎管的等效應力幅 Δσ。綜合高低週疲勞結果，未經處理的 3D 冷彎管在海洋深冷環境下的疲勞壽命將呈現斷崖式縮減。

八、結論與工程建議

本研究建立的高保真度三維有限元素熱機耦合框架，全面且深入地解析了 316L 不銹鋼管在「冷彎成形—深冷裝卸—海洋服役」全壽命週期中的多物理場行為，涵蓋了水下管線、LNG 船舶與碼頭裝卸臂系統。得出以下核心結論與工程建議：

大彎曲半徑與熱處理的絕對必要性：無論是船舶機艙內部管線、碼頭裝卸臂還是海底管線，3D 冷彎引入的高達 400 MPa 的殘留應力與嚴重的壁厚減薄，會將 316L 材料推向塑性耗竭極限。強烈建議在工程實務上全面採用 5D 或更大半徑的成形技術。更重要的是，必須嚴格實施彎曲後熱處理（Post-Bend Heat Treatment, PBHT）以徹底釋放成形初始殘留應力。唯有如此，才能避免系統落入 ASME 規範的塑性棘輪區而遭受K_e,k 懲罰因子的致命打擊。
熱循環與微觀相變的雙面刃效應：-162°C 的 LNG 裝卸熱循環會強烈觸發 316L 的γ→ α’ 馬氏體相變。相變伴隨的 TRIP 效應雖能透過體積膨脹來鬆弛宏觀熱應力，但新舊相界面的微觀彈性錯配應力卻大幅增加了低週疲勞條件下的微裂紋萌生風險。
遵循 IGC Code 的進階疲勞設計範式：LNG 船岸系統承受著極端溫度梯度與複雜的三維機械位移。為徹底遵循 IMO IGC Code 與 ASME 規範對安全性的嚴苛要求，設計與分析單位必須摒棄傳統簡化的純彈性分析，改採包含微觀相變動力學（如 Olson-Cohen 模型）與高階循環塑性（如 Chaboche 模型）的進階有限元素框架，以確保跨洋深冷能源傳輸網絡在全壽命週期內的絕對結構完整性與安全性。

參考文獻

Analysis of tensile mechanical properties of LNG cryogenic flexible pipe – ResearchGate, https://www.researchgate.net/publication/403352179_Analysis_of_tensile_mechanical_properties_of_LNG_cryogenic_flexible_pipe
What are Swivel Joints? Working, Types & Engineering Guide 2026 – EPCLand, https://epcland.com/engineering-guide-swivel-joints/
Cryogenic hydrogen embrittlement of 316plus (EN 1.4420) stainless steel at 77 K and 20 K, https://arxiv.org/html/2604.08751v2
Cryogenic Applications: Best Stainless Steel Grades for Liquid Nitrogen Storage, https://seathertechnology.com/cryogenic-applications-best-stainless-steel-liquid-nitrogen/
RESOLUTION MSC.370(93) (adopted on 22 May 2014) AMENDMENTS TO THE INTERNATIONAL CODE FOR THE CONSTRUCTION AND EQUIPMENT OF SHIPS, https://wwwcdn.imo.org/localresources/en/KnowledgeCentre/IndexofIMOResolutions/MSCResolutions/MSC.370(93).pdf
GUIDELINES ON CARRIAGE OF LIQUEFIED HYDROGEN IN BULK – Irclass.org, https://www.irclass.org/media/8034/guidelines-on-carriage-of-liquefied-hydrogen-in-bulk_sept-2025_new.pdf
International Code for the Construction and Equipment of Ships Carrying Liquefied Gases in Bulk (IGC Code), https://treaties.fcdo.gov.uk/data/Library2/pdf/1998-TS0015.pdf
Induction Bending Pipe for Piping Systems – abter steel pipe manufacturer, natural gas casing and tubing,seamless steel pipe,OCTG,, https://www.abtersteel.com/pipe-fittings/induction-bending-pipe-for-piping-systems/
Hot Induction Bends and long radius 3D – reli engineering, https://www.relieng.com/hot-induction-bends.html
(PDF) Deformation microstructures and martensitic transformation pathways in cryogenically deformed 316L stainless steel – ResearchGate, https://www.researchgate.net/publication/377727076_Deformation_microstructures_and_martensitic_transformation_pathways_in_cryogenically_deformed_316L_stainless_steel
Phase Transformation in 316L Austenitic Steel Induced by Fracture at Cryogenic Temperatures: Experiment and Modelling – PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC7795462/
The behaviour and deformation mechanisms for 316L stainless steel deformed at cryogenic temperatures – Oak Ridge National Laboratory, https://impact.ornl.gov/en/publications/the-behaviour-and-deformation-mechanisms-for-316l-stainless-steel/
Deformation-induced martensitic transformation in fused filament fabrication austenitic stainless steels during tension at wide – IPPT PAN, https://www.ippt.pan.pl/repository/open/o10726.pdf
Martensitic Transformation and Strengthening Mechanism in a 304 Stainless Steel Subjected to Wire Drawing – MDPI, https://www.mdpi.com/1996-1944/19/11/2412
(PDF) Phase Transformation in 316L Austenitic Steel Induced by Fracture at Cryogenic Temperatures: Experiment and Modelling – ResearchGate, https://www.researchgate.net/publication/348089367_Phase_Transformation_in_316L_Austenitic_Steel_Induced_by_Fracture_at_Cryogenic_Temperatures_Experiment_and_Modelling
Microstructural Characterization and Mechanical Properties of L-PBF Processed 316 L at Cryogenic Temperature – PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC8510410/
Damage development during the strain induced phase transformation of austenitic stainless steels at low temperatures – OUCI, https://ouci.dntb.gov.ua/en/works/4Nq0Zo67/
An Alternative Model for Investigating the Kinetics of Deformation-Induced Martensitic Transformation – Iranian Journal of Materials Forming, https://ijmf.shirazu.ac.ir/article_8603_7123dba2ac980cd655924910ba4d0c01.pdf
Deformation-induced martensitic transformation kinetics in X5CrMnN17-8 austenitic stainless steel, https://journal.ut.ac.ir/article_105032_3f8795aca55ce63e71aa08ae95dffdbf.pdf
Cyclic Plasticity and Low Cycle Fatigue of an AISI 316L Stainless Steel: Experimental Evaluation of Material Parameters for Durability Design – PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC8269699/
Cyclic Complex Loading of 316 Stainless Steel: Experiments and Calculations – NC State Repository, https://repository.lib.ncsu.edu/bitstreams/6d61aa04-4ed6-4daf-8976-e632440df6c8/download
The tensile test of 316 L at room temperature, experimental data from Chaboche and Rousselier (1983) – ResearchGate, https://www.researchgate.net/figure/The-tensile-test-of-316-L-at-room-temperature-experimental-data-from-Chaboche-and_fig1_225745297
Cyclic Ramberg–Osgood Parameters of 316SS Weld Metal Under In-Air and PWR-Water Environments – ASME Digital Collection, https://asmedigitalcollection.asme.org/pressurevesseltech/article/143/5/051507/1098844/Cyclic-Ramberg-Osgood-Parameters-of-316SS-Weld
Low-Cycle Fatigue Behavior and Microstructural Damage Mechanisms of 316L Austenitic Stainless Steel in Cryogenic Environments – MDPI, https://www.mdpi.com/1996-1944/19/12/2494
Simulation of Cyclic Behavior of Cold Worked 316L and Solution Annealed 316L, https://www.kns.org/files/pre_paper/16/229%EA%B9%80%EC%A2%85%EB%B2%94.pdf
316 316L 316ti Stainless Steel Grade Complete Guide, https://www.ambicasteels.com/blog/316-316l-316ti-stainless-steel-grade/
Stainless Steel – Austenitic – 1.4404 (316L) Bar and Section – Aalco, https://www.aalco.co.uk/datasheets/Stainless-Steel-14404-316L-Bar-and-Section_39.ashx
Pipe Bending Wall Thinning Calculator | B31.3 Minimum Radius – Industrial Monitor Direct, https://industrialmonitordirect.com/blogs/knowledgebase/calculating-pipe-wall-thinning-and-minimum-bend-radius-per-b313
unknown_url
Characterization of Residual Stress and Plastic Strain in Austenitic Stainless Steel 316L(N) Weldments – Open Research Online, https://oro.open.ac.uk/61483/1/13835622.pdf
LNG Unloading Arms Technical Specs | PDF – Scribd, https://www.scribd.com/document/480279876/L1-R-1132-1003-Att-1-lng-unloading-arms-technical-specification-0
Specification Sheet: Alloy 316/316L – Sandmeyer Steel Company, https://www.sandmeyersteel.com/wp-content/uploads/316-316l-317l-spec-sheet.pdf
Assessment of the Morison Equation. – DTIC, https://apps.dtic.mil/sti/tr/pdf/ADA088185.pdf
Comparison of Potential Theory and Morison Equation for Deformable Horizontal Cylinders, https://journals.nasspublishing.com/index.php/sms/article/download/492/533/3446
Tutorial on Wave Loading Analysis by CAEPIPE, https://www.sstusa.com/tutorials/tutorial_on_wave_loading_analysis_by_caepipe.htm
Analysis of Wind–Wave Relationship in Taiwan Waters – MDPI, https://www.mdpi.com/2077-1312/13/6/1047
Reconstructing hourly coastal total sea levels and assessing current and future extreme sea levels threats to the Coast of China – PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC12586587/
Influences of Nononshore Winds on Significant Wave Height Estimations Using Coastal X-Band Radar Images – Ncku.edu.tw, https://www.comc.ncku.edu.tw/wordpress/wp-content/uploads/2023/02/111%E5%B9%B4%E5%BA%A6-Influences-of-Nononshore-Winds-on-Significant-Wave-Height-Estimations-Using-Coastal-X-Band-Radar-Images.pdf
Analysis of Long-Period Hazardous Waves in the Taiwan Marine Environment Monitoring Service – Frontiers, https://www.frontiersin.org/journals/marine-science/articles/10.3389/fmars.2021.657569/full
(PDF) Cyclic Plasticity and Low Cycle Fatigue of an AISI 316L Stainless Steel: Experimental Evaluation of Material Parameters for Durability Design – ResearchGate, https://www.researchgate.net/publication/352793148_Cyclic_Plasticity_and_Low_Cycle_Fatigue_of_an_AISI_316L_Stainless_Steel_Experimental_Evaluation_of_Material_Parameters_for_Durability_Design
Microstructural mechanisms of cyclic deformation, fatigue crack initiation and early crack growth | Philosophical Transactions of the Royal Society A, https://royalsocietypublishing.org/rsta/article/373/2038/20140132/114789/Microstructural-mechanisms-of-cyclic-deformation
(PDF) An Energy-Based Unified Approach to Predict the Low-Cycle Fatigue Life of Type 316L Stainless Steel under Various Temperatures and Strain-Rates – ResearchGate, https://www.researchgate.net/publication/332178446_An_Energy-Based_Unified_Approach_to_Predict_the_Low-Cycle_Fatigue_Life_of_Type_316L_Stainless_Steel_under_Various_Temperatures_and_Strain-Rates
Multiaxial fatigue behavior of AISI 316L subjected to strain-controlled and ratcheting paths – Re.Public@polimi, https://re.public.polimi.it/bitstream/11311/964505/4/Multiaxial%20fatigue%20behavior%20of%20AISI%20316L_11311-964505_Foletti.pdf
Multiaxial Fatigue in 316L Steel Models | PDF – Scribd, https://www.scribd.com/document/891406266/Critical-Plain
Experimental evaluation and modeling of fatigue of a 316L austenitic stainless steel in high-temperature water and air environments, https://infoscience.epfl.ch/entities/publication/2e7a1821-54b6-4a5a-a927-86761703e678/statistics
Cyclic Deformation and Low-Cycle Fatigue for 316LN Stainless Steel under Non-proportional Loading – MATEC Web of Conferences, https://www.matec-conferences.org/articles/matecconf/pdf/2019/49/matecconf_icmff1218_08002.pdf
ASME BPVC 2025 Bree Diagram Override of Fatigue Exemption – Industrial Monitor Direct, https://industrialmonitordirect.com/blogs/knowledgebase/asme-bpvc-2025-bree-diagram-override-of-fatigue-exemption
Fundamentals of Design by Analysis – The Equity Engineering Group, Inc., https://e2g.com/industry-insights-ar/fundamentals-of-design-by-analysis/
Basics of Design By Analysis in ASME Section VIII, Division 2 – Becht, https://becht.com/becht-blog/entry/basics-of-design-by-analysis-in-asme-section-viii-division-2/
Investigations on Fatigue Life of Tube Connections Based on International Codes of Pressure Vessel – PMC, https://pmc.ncbi.nlm.nih.gov/articles/PMC9822275/
Technical Bulletin: Fatigue Code Update – Impact to Monopile Designs – Wood Thilsted, https://woodthilsted.com/news/fatigue-code-update-technical-bulletin/
DNVGL-RP-0005 RP-C203: Fatigue design of offshore steel structures – UNECE Wiki, https://wiki.unece.org/download/attachments/26902968/R55_11_03-C203_DNVGL-RP-0005_2014-06.pdf?api=v2
How to Simplify DNV RP-C203 Compliance in Offshore Structures – SDC Verifier, https://sdcverifier.com/articles/dnv-rp-c203-compliance-offshore-structures/
SN Curves – Cambridge University Press & Assessment, https://resolve.cambridge.org/core/services/aop-cambridge-core/content/view/7C3E22739F9E3EAFAFE23D4C6092D70B/9781316343982c4_p123-173_CBO.pdf/sn_curves.pdf
How to set classification for DNV RP-C203 – SDC Verifier, https://sdcverifier.com/articles/how-to-set-classification-for-dnv-rp-c203/
Comparison of Fatigue Provisions in Various Codes and Standards -Part 1: Basic Design SN Curves of Non-Tubular Steel Members – Journal of Ocean Engineering and Technology, https://www.joet.org/journal/view.php?number=3010

一、 摘要與緒論