從工程應力應變曲線到仿真材料卡片：一位CAE工程師的實戰筆記

國高材高分子材料產業創新中心

2026年5月15日 16:45

材料卡片是仿真分析的"基因"，決定了有限元計算結果的精度上限。

在碰撞仿真、NVH分析、產品可靠性評估等場景中，材料參數設置的準確性直接影響仿真的可信度。然而，實驗室提供的原始材料曲線與仿真軟件所需的有效應力應變曲線之間，存在一道需要跨越的轉化鴻溝。本文基于實戰經驗，系統梳理從材料曲線獲取到仿真材料卡片生成的完整流程，供從事CAE工作的工程師參考。

一、工程應力應變曲線

1.1 材料的關鍵參數

開展有限元分析前，必須明確材料的幾項基礎參數，這些參數構成了材料卡片的骨架。

彈性模量（楊氏模量）是工程應力應變曲線屈服段的斜率，即應力與應變的比值。金屬材料通常為210000 MPa或20600 MPa，塑料材料約為2350 MPa。這一參數直接決定了結構在彈性階段的剛度表現。

圖1 帶引伸計拉伸測試

泊松比是材料在單向受拉或受壓時，橫向正應變與軸向正應變的比值，用于反映材料的橫向變形特性。金屬材料泊松比通常取0.34，塑料材料約為0.39。密度是質量與體積的比值，在碰撞仿真和NVH分析中尤為重要——不同單位制模型中，密度參數容易出現數量級錯誤，導致分析結果嚴重失真。

屈服強度是材料從彈性變形進入塑性變形的臨界點。拉伸過程中，材料在屈服點之前僅產生彈性變形；過了屈服點則進入塑性階段，產生永久不可恢復的變形。塑料材料由于韌性較差，拉伸試驗中基本沒有明顯的屈服階段，工程設計中常以產生0.2%殘余應變時的應力作為條件屈服極限。

抗拉強度是材料應力值的極限點，超過此值材料即被判定破壞失效。斷裂延伸率則是抗拉強度所對應的應變值，塑性應變值超過斷裂延伸率時，材料同樣被視為失效。

圖2 應力應變曲線

1.2 獲取途徑

工程應力應變曲線的獲取主要有三種途徑，各有優劣。

第一種方式是向材料供應商直接索取，這是最理想的信息來源，尤其對于成熟牌號的商業材料，供應商通常能提供完整的測試報告。

第二種方式是委托第三方實驗室進行拉伸試驗，這種方法獲得的數據最為準確可靠，但成本較高，適用于對仿真精度要求極高的關鍵零部件。

第三種方式是通過專業工具從已發表的技術文獻或網絡資源中"白嫖"曲線數據，再利用數字化工具提取坐標點，這種方式成本最低但數據質量參差不齊，僅推薦用于項目前期的快速可行性分析階段。

二、工程曲線→真實曲線→有效曲線

2.1 為什么必須轉換？

實驗室萬能試驗機直接輸出的拉伸曲線稱為工程應力應變曲線，其定義方式為：工程應力 = 力 / 原始截面積；工程應變 = 伸長量 / 原始標距長度。這種表達方式假設樣條在整個拉伸過程中截面積不變，與實際情況存在偏差。

CAE仿真軟件（以LS-Dyna為例）使用的則是有效應力應變曲線，這條曲線需要滿足兩個條件：一是真實反映材料在大變形階段的應力-應變關系；二是曲線形態必須單調遞增，以便于數值計算。因此，從工程曲線到有效曲線需要經過兩次數學轉換。

2.2 第一次轉換：工程曲線→真實曲線

真實應力與工程應力的轉換公式為：

真實應變與工程應變的轉換公式為：

這一轉換的本質是引入瞬時截面積的概念。當材料被拉伸時，樣條的截面積隨著變形而減小，因此真實的應力值實際上高于按原始截面積計算的工程應力值。轉換后的真實應力應變曲線已經呈現出單調遞增的形態。

2.3 第二次轉換：真實曲線→有效曲線

在塑性大變形分析中，有效應力應變曲線采用等效應力的概念進行計算。對于單軸拉伸情況，有效應力與真實應力之間存在以下關系：

經過這兩次轉換得到的有效應力應變曲線，才能真正作為LS-Dyna等仿真軟件的輸入數據使用。

三、聚合物的高非線性與粘塑性力學特征

與具有明確晶格滑移機制的金屬不同，非晶態與半結晶態聚合物（如PC, ABS, PP）的變形源于高分子鏈段的滑移、解纏結與取向。這種微觀機制使得聚合物在宏觀上表現出極其復雜的力學特征：強烈的靜水壓敏感性（拉壓屈服不對稱，壓縮屈服強度往往遠高于拉伸）、顯著的粘彈性/粘塑性耦合響應、極低應變率下的頸縮后冷拉（Cold Drawing）現象，以及伴隨微裂紋（Crazing）與剪切帶（Shear Banding）競爭的損傷演化。

在構建聚合物材料卡片時，傳統的金屬本構模型完全失效。工程界目前傾向于采用兩類策略：

第一類是基于Drucker-Prager或Mohr-Coulomb這類原本用于巖土材料的屈服準則，通過引入靜水壓力項來修正拉壓不對稱性；

第二類則是采用專為聚合物開發的半解析模型，如SAMP-1（Semi-Analytical Model for Polymers）。SAMP-1模型允許用戶直接輸入單軸拉伸、單軸壓縮、雙軸拉伸及純剪切四條不同應力狀態下的屈服曲線，并根據加載路徑自動插值構建動態的三維屈服面。

同時，該模型引入了非關聯流動法則（Non-associated Flow Rule）以準確控制聚合物在塑性變形過程中的體積膨脹（即泊松比隨塑性應變的變化），這對于精確預測塑料扣位的插拔失效與手機外殼的跌落開裂至關重要。

四、材料卡片應用邏輯

各大主流商業求解器（如LS-DYNA, Abaqus, PAM-CRASH）在底層動力學積分算法上殊途同歸，但在材料卡片的關鍵字定義、輸入邏輯與容錯處理機制上存在顯著差異。深入理解這些差異，是避免“垃圾輸入，垃圾輸出”（Garbage In, Garbage Out）的關鍵。

LS-DYNA：顯式碰撞仿真的材料庫開發經驗

LS-DYNA以其在強非線性、大變形和瞬態沖擊領域的絕對統治力而聞名。其內置材料庫極其龐雜，包含逾三百種材料模型，每一種模型都針對特定的物理場景進行了底層代碼級的極速優化。

*MAT_024 (MAT_PIECEWISE_LINEAR_PLASTICITY)： 這是整個汽車工業與消費電子防撞設計中使用頻率最高的通用彈塑性卡片。它的魅力在于卓越的魯棒性。工程師可以通過LCSS表格輸入離散的真實應力-真實塑性應變曲線族來定義應變率效應。實踐經驗警告： LCSS表格的插值機制是基于自然對數插值的。如果輸入的應變率曲線出現交叉（即高應變率下的應力低于低應變率下的應力），或者硬化曲線呈現負斜率（未激活損傷模塊時），求解器的材料剛度矩陣將出現非正定，導致不可控的網格畸變。此外，必須通過外推確保表格覆蓋到極高應變率（如10000 /s），以防求解器在局部高變形區發生錯誤的常數外推。
*MAT_081 (MAT_PLASTICITY_WITH_DAMAGE) 與 MAT_SAMP-1 (MAT_187)： 在涉及破壞仿真的塑料件（如保險杠、手機外殼）中，MAT_187是業界公認的最佳實踐。SAMP-1模型內部的屈服面隨拉壓狀態動態改變，這要求輸入的數據必須在原點保持嚴格的相切連續性。如果單軸拉伸與純剪切曲線在微小塑性應變處的模量不匹配，求解器會在極短時間內由于屈服面不封閉而崩潰。
*復合材料損傷 MAT_054 (MAT_ENHANCED_COMPOSITE_DAMAGE)： 針對碳纖維層合板（CFRP），MAT_054利用Chang-Chang失效準則分別判斷基體與纖維的拉壓破壞。由于復合材料的極度脆性，單元失效極易引發應力波的虛假反射。工程實踐中，必須精細調節DFAIL（失效應變控制）與SOFT（軟化系數控制）參數，同時強制約束單元的最小破壞時間步，以防止仿真因為局部高頻振蕩而中止。

Abaqus：從隱式非線性到用戶子程序的深度定制

Abaqus采用極其模塊化的*MATERIAL關鍵字樹狀結構，使得多物理場耦合特性的定義更加符合人類直覺。

金屬可塑性與損傷演化： 在Abaqus中定義延性金屬，往往通過組合*ELASTIC, *PLASTIC, *RATE DEPENDENCE, 以及延性損傷起始*DAMAGE INITIATION, CRITERION=DUCTILE。Abaqus的強大之處在于其引入了斷裂能正則化機制*DAMAGE EVOLUTION, TYPE=ENERGY。這一機制徹底改變了傳統材料卡片隨網格尺寸變小而急劇變“脆”的網格敏感性缺陷，使得能量耗散成為一個相對客觀的物理不變量。
UMAT / VUMAT 的二次開發： 當標準材料庫無法覆蓋新興材料（如具有形狀記憶效應的鎳鈦合金、相變誘發塑性的TRIP鋼、或者超高周疲勞退化材料）時，最高階的仿真工程師必須依賴Fortran或C++編寫用戶自定義材料子程序（UMAT用于Abaqus/Standard隱式求解，VUMAT用于Abaqus/Explicit顯式求解）。編寫UMAT的核心挑戰并不在于屈服函數的編碼，而在于推導高度復雜的“一致切線剛度矩陣”（Consistent Tangent Modulus）。隱式非線性求解嚴重依賴該矩陣進行牛頓迭代，如果切線剛度推導存在微小誤差，將導致模型在屈服點附近徹底喪失二次收斂性（Quadratic Convergence），陷入無盡的迭代發散死循環。

材料卡片定制