第四篇:精準建模:HJC 材料模型與 EFP 彈丸成型技術之深度解析
第四篇:精準建模:HJC 材料模型與 EFP 彈丸成型技術之深度解析
在現代防護工程與毀傷彈藥設計領域,針對非均質脆性材料(如混凝土)與金屬大變形物理過程(如爆炸成型彈丸)的數值模擬,已成為科研與工程應用的核心手段。要準確預測高速衝擊載荷下材料的動態響應,不僅需要對連續介質力學有深刻理解,更需要掌握特定組構模型(Constitutive Models)的參數標定與數值算法的邊界條件。
本篇將深度探討專為混凝土設計的 Holmquist-Johnson-Cook (HJC) 材料模型,並結合 爆炸成型彈丸 (Explosively Formed Projectile, EFP) 的物理機制,詳解 3D 數值建模過程、參數敏感性分析及其在技術諮詢中的應用價值。
1. HJC 混凝土材料模型:描述脆性材料的數學語言
針對混凝土這類具備高度壓力依賴性、應變率敏感性以及體積壓縮特性的材料,Holmquist、Johnson 和 Cook 於 1993 年提出的 HJC 模型提供了一個兼顧計算效率與物理精度的框架。該模型主要由三個核心部分組成:
(1) 等效強度模型 (Yield Surface)
HJC 強度模型的物理核心在於將標準化等效應力 $\sigma^$ 定義為壓力、應變率與損傷的耦合函數。其表達式精確捕捉了混凝土在高壓環境下剪切強度提升的特點: $$ \sigma^ = [A(1 - D) + BP^{N}] [1 + C \ln \dot{\epsilon}^] $$
- $A$ (標準化內聚強度):決定了材料在無壓力狀態下的基礎強度。
- $B$ (壓力硬化係數) 與 $N$ (硬化指數):描述了隨著靜水壓力增加,混凝土變得越「硬」的物理現象(即圍壓效應)。
- $C$ (應變率係數):控制動態載荷下的強度增益。這解釋了為何在高速撞擊下,混凝土表現出的強度遠高於靜態測試。
(2) 損傷累積機制 (Damage Evolution)
混凝土的毀傷並非瞬時發生,而是微裂紋演化與孔隙塌陷的累積過程。HJC 模型透過累積等效塑性應變與塑性體積應變來計算損傷度 $D$(從 0 到 1)。 $$ D = \sum \frac{\Delta \epsilon_p + \Delta \mu_p}{\epsilon_p^f + \mu_p^f} $$ 這種雙變量累積機制使得模型能夠區分「剪切破壞」與「壓實破壞」,這對於模擬穿甲彈藥在混凝土靶板內形成的漏斗狀坑洞至關重要。
(3) 三階段狀態方程 (Equation of State, EOS)
為了模擬炸藥近場產生的衝擊波傳播,HJC 模型採用了分段式的狀態方程來描述壓力 $P$ 與體積應變 $\mu$ 的關係:
- 彈性階段 ($P \le P_{crush}$):材料表現為線性彈性,孔隙尚未塌陷。
- 過渡壓實階段 ($P_{crush} < P < P_{lock}$):模擬混凝土內部的空氣孔隙被逐步壓實排出。此階段伴隨著顯著的體積能量耗散,是混凝土抗炸吸能的關鍵。
- 完全壓實階段 ($P \ge P_{lock}$):孔隙完全消失,材料表現為緻密固體的非線性壓縮行為。
2. 爆炸成型彈丸 (EFP):藥型罩的極限變形藝術
爆炸成型彈丸 (EFP) 是一種特殊的聚能裝藥結構,其核心目標是透過藥型罩(Liner)的受控變形,形成一個具備良好空氣動力學形狀、高度集中且具備極高動能的單體彈丸。
成型原理:從「碟」到「彈」的瞬間蛻變
與形成細長金屬射流(Jet)的錐形罩不同,EFP 通常採用球缺形或淺盤狀藥型罩。當炸藥起爆,爆轟波以超過 8000 m/s 的速度撞擊藥型罩時,會發生以下物理過程:
- 衝擊加速 (Impact Phase):爆轟波壓力極高(>20 GPa),迫使金屬材料在微秒級時間內發生塑性流動。
- 閉合與翻轉 (Closing Phase):由於藥型罩邊緣較薄且受到爆轟產物的側向擠壓,邊緣速度往往高於中心速度。這導致藥型罩向軸線方向潰縮並向後翻轉(或向前摺疊),將原本的盤狀結構「捏」成一個實心的金屬彈丸。
- 拉伸成型 (Stretching Phase):由於頭部與尾部的速度梯度,彈丸在飛行初期會被拉長,形成類似「鳥頭」或「長桿」的形狀。
3. 3D 數值模型的建立與參數設定
在 LS-DYNA 或 AUTODYN 等顯式動力學軟體中,建立 EFP 侵徹混凝土的高精度模型需要精細的策略。
(1) 網格策略:Lagrangian, Eulerian 與 SPH 的共舞
- EFP 成型階段:推薦使用 ALE (Arbitrary Lagrangian Eulerian) 算法。它允許網格隨邊界調整,能有效解決藥型罩在大變形過程中的單元畸變問題,精確預測彈丸初速與形狀。
- 侵徹混凝土階段:推薦使用 SPH (Smoothed Particle Hydrodynamics) 算法。混凝土在高壓下會發生劇烈的破碎與飛濺,傳統網格容易崩潰。SPH 無網格法能自然地模擬出混凝土的飛散碎片雲與深孔侵徹過程。
(2) 關鍵材料參數設定
參數的準確性直接決定了模擬的可信度。以下是基於文獻的典型參數設定參考:
-
混凝土 (HJC 模型):
- 密度 ($\rho$):約 2.34 - 2.50 g/cm³。
- 剪切模數 ($G$):對於普通強度混凝土約 14-15 GPa;對於高強混凝土(>140 MPa)可達 22 GPa。
- 抗壓強度 ($f_c'$):依據實際材料設定(如 48 MPa 或更高)。
- 狀態方程參數:$P_{crush}$(壓碎壓力)通常設定為 $f_c'/3$,$\mu_{lock}$(鎖定體積應變)約為 0.1,反映孔隙率。
-
藥型罩 (銅, Johnson-Cook 模型):
- 需定義 $A$ (屈服強度), $B$ (硬化模數), $C$ (應變率常數), $m$ (熱軟化指數) 等參數。
- 例如 OFHC 銅的典型參數:$A=90$ MPa, $B=292$ MPa, $n=0.31$。
-
炸藥 (JWL 狀態方程):
- 使用 JWL 方程描述爆轟產物的壓力-體積關係。例如常用炸藥 Octol 的爆速約 8480 m/s,密度約 1.821 g/cm³。
結語
透過 HJC 模型與 EFP 成型技術的深度整合,工程師可以在計算機中「預演」極端衝擊下的物理場景。這不僅能大幅減少實彈測試的高昂成本,更能透過參數化分析(如調整藥型罩厚度或混凝土骨材強度),設計出更高效的毀傷元或更堅固的防護結構。
參考文獻:
- Strain-rate effects associated with the HJC concrete model - EPJ Web of Conferences
- Evaluation of deflagration fracturing effect based on the HJC constitutive model - PMC
- Johnson-Holmquist Concrete Model - LS-DYNA Theory Manual
- 碩士論文:彈體以高速撞擊鋼纖維混凝土之物理行為分析 (Tab 1, Tab 6)
- Numerical simulations of the formation behavior of explosively formed projectiles
- 精準建模:HJC 材料模型與 EFP 彈丸成型技術之深度解析
- Formation Behaviors of Coated Reactive Explosively Formed Projectile - PMC
留言
張貼留言