1.1、打開ANSYS工作臺，創建一個“顯式動力學”分析，檢查各個單元。我們將使用默認的結構鋼作為鈑金，并添加一種雙線性各向同性硬化，屈服強度為470MPa，切線模量為1000MPa。 1.2、導入幾何體(見圖1)。 圖 1 鈑金成型模型的幾何形狀 1.3、網格化模型。金屬板材初始厚度為3毫米。將機器部件改為剛體，僅保留鈑金作為柔性體。

功能梯度材料分析 在功能梯度材料（FGM）梁的彎曲分析中，單元通過材料屬性的厚度方向插值，可模擬彈性模量的連續變化，準確計算沿厚度方向的應力梯度。單元計算的正應力（σ_x）和橫向切應力（τ_xz）與超細網格高階實體單元（C3D20）結果的偏差均小于 3%，但同時后者會增加巨大的計算量！

屈服強度360 MPa，硬化模量0（理想塑性） 如果材料有硬化（如 H' = 1000 MPa）： TBDATA,,360,1000 ! 屈服后模量1000 MPa （此時，塑性階段的切線模量 = E_t = H'） 總結： TB,BKIN,2 → 定義材料2為雙線性隨動強化模型。 TBDATA,,360,0 → 設置屈服強度360 MPa，無硬化（理想塑性）。

時間積分方案策略： 過程為： 通過線性方程組單次迭代求解滿足收斂判據的pk2應力以及位錯密度 收斂判據 為了避免位錯密度波動對數值穩定性的影響，修正收斂判據為： 收斂之后更新柯西應力： 一致性切線模量： 數值模型驗證： 材料參數： 不同取向單晶鋁合金拉伸模擬與實驗結果對比

而Ansys、Abaqus暫時沒有虛擬質量方法，都是直接采用基于聲學有限元的流固耦合來求濕模態。 本章只介紹基于虛擬質量的濕模態計算。

https://www.yqgqt.org.cn/content/post/1195061 第十七篇：幾何非線性的物理含義。 https://www.yqgqt.org.cn/content/post/1198459 第十八篇：幾何非線性的應變。

對于求 解循環對稱模態，ANSYS也提供了專用的求解宏指令（不可直接用solve命令），該 指令格式為：CYCSOL,NDMIN,NDMAX,NSECTOR, LOW 各參數含義如下： NDMIN、NDMAX:計算的上下節徑范圍,NDMIN最小為0,NDMAX 對偶數 最大可取n/2,對奇數最大可取（n-1）/2。 NSECTOR:循環對稱的扇區數，這里為4。

溫度網格單元的節點為標量，自由度為1，形函數物理含義明確，單剛和總剛計算都比較簡單，不像電熱耦合雙向傳遞，從熱場到結構是單向傳遞，也符合由易到難的順序。

插入靜態結構 材料 編輯結構鋼 在塑性下，添加雙線性各向同性硬化 屈服強度，2.07e8 Pa 切線模量，6.92e8 Pa 創建新材料，稱為“Balloon” 在超彈性下，添加 Mooney-Rivlin

材料楊氏模量=1.5e6，泊松比=0.3。