寫在前文

  嗨！老朋友們~~~又再一次與大家分享！隔了這么久沒冒泡，大家還好嗎？筆者近期在整理相關研究資料時，系統梳理了 Abaqus 中實體單元的分類邏輯、理論基礎及不同場景下的選擇策略，發現現有實踐中有粉絲仍存在單元類型誤用、特性理解不充分等問題。鑒于此，本文將從單元分類、選擇原則、特定場景應用及最佳實踐等方面展開論述，旨在為從事 Abaqus 仿真分析的研究者與工程技術人員提供系統性參考，助力提升數值模擬的可靠性與科學性。

  在現代工程仿真領域，Abaqus 作為一款功能強大的有限元分析軟件，提供了豐富的單元庫來滿足各種結構力學分析需求。特別是在三維應力分析場景中，正確選擇和使用單元類型對于獲得準確、高效的計算結果至關重要。Abaqus 單元庫按照單元族分類，主要包括實體單元 (C)、殼單元 (S)、梁單元 (B)、桁架單元 (T)、剛體單元 (R) 等。每個單元族又包含多種具體單元類型，適用于不同的幾何特征、載荷條件和分析目標。

   本篇是三維應力單元，即實體單元篇！

【相關閱讀】

【JY】有限元分析的單元類型分享一波~

【JY】有限單元分析的常見問題及單元選擇

1 實體單元分類與理論基礎

  實體單元是 Abaqus 中最基礎也最常用的單元類型，可在其任何表面與其他單元連接，能夠精確地反映幾何形狀的復雜性，適用于模擬具有復雜形狀的結構。在 Abaqus 中，實體單元根據節點位移插值階數、積分方式和特殊功能可分為多種類型。

1.1 按節點位移插值階數分類

根據節點位移插值階數，實體單元主要分為三類：

1. 線性單元（一階單元）：節點僅布置在單元角點，各方向采用線性插值，計算相對簡單，適用于變形較簡單的情況。線性單元的主要優勢在于計算效率高，但精度相對較低，適用于對精度要求不高或初步分析階段。
2. 二次單元（二階單元）：每條邊上有中間節點，采用二次插值，能更好地擬合復雜變形，精度更高，用于對精度要求高、變形復雜的結構分析。二次單元能夠更準確地模擬應力集中區域的應力變化，特別適合高精度分析場景。

1. 修正的二次單元：僅三角形 (Tri) 和四面體 (Tet) 單元有此類型，邊上有中間節點，采用修正的二次插值，針對三角形和四面體單元特性優化。修正的二次單元在保持高精度的同時，改善了常規二次單元的一些不足，如對接觸分析和大變形的適應性。

1.2 按積分方式分類

根據積分方式，實體單元可分為：

1. 完全積分單元（C3D8、C3D20等）：使用足夠數量的高斯積分點精確積分單元剛度矩陣中的多項式。當單元形狀規則時，完全積分單元能夠提供精確的應力結果，但承受彎曲載荷時會出現剪切自鎖現象，導致單元過剛，即使網格細密，計算精度也不佳，不適用于主要承受彎曲載荷的結構。
2. 縮減積分單元（XXX-R，如C3D8R）：相比普通完全積分單元，在每個方向少用一個積分點，減少計算量，提升計算效率。線性縮減積分單元僅在單元中心有一個積分點，存在沙漏數值問題，使單元表現過于柔軟。在模擬承受彎曲載荷結構時，沿厚度方向至少劃分四個單元，可一定程度減輕沙漏效應影響。
3. 非協調模式單元（XXX-I，如C3D8I）：僅在 Abaqus/Standard 中存在，可解決線性完全積分單元的剪切自鎖問題。與其他軟件中的 4 節點四邊形單元或 8 節點六面體單元相似，在 Abaqus 中選用該單元，計算結果與其他軟件有較好的一致性。非協調單元成功克服了剪切自鎖問題，在單元扭曲程度較小時，位移和應力計算結果精確。
4. 雜交單元（XXX-H，如C3D8RH）：僅在 Abaqus/Standard 中存在，每種實體單元都有對應的雜交單元，Abaqus/Explicit 中沒有雜交單元。雜交單元用于模擬不可壓縮材料（如泊松比為 0.5 的橡膠）或近似不可壓縮材料（泊松比大于 0.475）。普通單元在模擬此類材料響應時（除平面應力問題外），單元中應力不確定，而雜交單元可有效解決該問題。

2 三維實體單元選擇策略

2.1 基本選擇原則

在選擇三維實體單元時，應遵循以下基本原則：

1. 幾何形狀優先原則：對于三維區域，盡可能采用結構化網格劃分技術或掃掠網格劃分技術，從而得到六面體 (Hex) 單元網格，減小計算代價，提高計算精度。當幾何形狀復雜時，也可以在不重要的區域使用少量楔形 (Wedge) 單元。
2. 單元類型優選：對于大多數工程問題，推薦使用的是六面體單元(Hex)、非協調單元，以及線性縮減積分單元。這些單元在精度和計算效率之間取得了較好的平衡。特別注意：Abaqus的所有單元均可用于動態分析，選取單元的一般原則與靜力分析相同。但在使用Abaqus/Explicit模擬沖擊或爆炸載荷時，應選用線性單元，因為它們具有集中質量公式，模擬應力波的效果優于二次單元所采用的一致質量公式。
3. 
4. 網格質量控制：網格質量直接影響計算結果的準確性，應盡量使用高質量的網格，避免畸形單元。單元形狀越接近規則形狀（如正方體），網格質量越高。長寬比過大的單元可能導致計算結果不準確，扭曲的單元會導致計算結果不穩定。
5. 單元密度優化：網格密度應根據應力梯度變化情況進行調整，應力梯度大的區域，網格密度應相應提高。在關鍵區域（如應力集中區域）進行局部加密，可以提高精度，同時節省計算資源。

2.2 不同求解器下的單元選擇差異

Abaqus/Standard 和 Abaqus/Explicit 在單元選擇上存在一定差異：

1. Abaqus/Standard：實體單元多樣，有二維和三維的線性、二次單元，積分方式可選完全積分或縮減積分。還包含修正的二次 Tri、Tet 單元，以及非協調模式單元和雜交單元。
2. Abaqus/Explicit：實體單元包括二維和三維的線性縮減積分單元，以及修正的二次 Tri 和 Tet 單元，但不存在二次完全積分實體單元。

2.3 特定場景下的單元選擇

針對不同的工程場景，三維實體單元的選擇策略也有所不同：

1. 應力集中問題：對于應力集中問題，盡量不要使用線性縮減積分單元，可使用二次單元來提高精度。如果在應力集中部位進行了網格細化，使用二次縮減積分單元與二次完全積分單元得到的應力結果相差不大，而二次縮減積分單元的計算時間相對較短。
2. 
3. 彈塑性分析：對于彈塑性分析，如果材料是不可壓縮性的（例如金屬材料），則不能使用二次完全積分單元，否則會出現體積自鎖問題，也不要使用二次 Tri 單元或 Tet 單元。推薦使用的是修正的二次 Tri 單元或 Tet 單元、非協調單元，以及線性縮減積分單元。如果使用二次縮減積分單元，當應變超過 20%~40% 時要劃分足夠密的網格。
4. 接觸分析或大扭曲變形：如果模型中存在接觸或大的扭曲變形，則應使用線性四邊形 (Quad) 或六面體 (Hex) 單元，以及修正的二次 Tri 單元或 Tet 單元，而不能使用其他的二次單元。
5. 彎曲為主的問題：對于以彎曲為主的問題，如果能夠保證在所關心部位的單元扭曲較小，使用非協調單元（例如 C3D8I 單元）可以得到非常精確的結果。
6. 不可壓縮材料：除了平面應力問題之外，如果材料是完全不可壓縮的（例如橡膠材料），則應使用雜交單元；在某些情況下，對于近似不可壓縮材料也應使用雜交單元。

3 各類三維實體單元詳解

3.1 線性完全積分單元

理論基礎：線性完全積分單元在每個方向上使用足夠的高斯積分點，以精確積分單元剛度矩陣中的多項式。當單元形狀規則時，能夠精確計算單元剛度矩陣。

適用場景：線性完全積分單元適用于模擬變形較簡單的結構，如主要承受拉伸或壓縮載荷的結構。由于其計算簡單，在對精度要求不高的初步分析中可以考慮使用。

優缺點分析：

• 優點：計算相對簡單，計算時間少，在規則網格和簡單變形情況下能提供合理結果。
• 缺點：承受彎曲載荷時會出現剪切自鎖現象，導致單元過剛，計算精度下降，不適用于主要承受彎曲載荷的結構。此外，對于復雜幾何形狀的適應性較差。

使用注意事項：

• 避免在主要承受彎曲載荷的結構中使用線性完全積分單元。
• 當必須使用時，應確保單元形狀盡可能規則，避免過度扭曲。
• 不適用于高精度分析場景，特別是存在應力集中或復雜變形的區域。

3.2 二次完全積分單元

理論基礎：二次完全積分單元在每條邊上有中間節點，采用二次插值函數，能夠更好地擬合復雜變形。通過足夠數量的高斯積分點精確積分單元剛度矩陣，能夠準確捕捉應力變化。

適用場景：二次完全積分單元適用于對精度要求高、變形復雜的結構分析，特別是存在應力集中的區域。它們能夠精確計算應力集中區域的應力分布，適用于裂紋擴展、缺口效應等需要高精度應力分析的場景。

優缺點分析：

• 優點：應力計算精準，能細致捕捉應力集中區域的應力變化；一般不會出現剪切自鎖現象，在多種受力工況下能合理反映單元的變形情況；計算結果可靠性高。
• 缺點：計算時間較長，計算成本高；不適用于接觸分析，接觸問題中接觸界面的復雜力學行為難以用該單元準確模擬；對于不可壓縮材料（如金屬）的彈塑性分析，易產生體積自鎖；當單元發生扭曲或彎曲應力存在梯度時，可能出現自鎖問題。

使用注意事項：

• 不適用于接觸分析或存在大變形的場景。
• 在彈塑性分析中，避免用于不可壓縮材料，否則易發生體積自鎖。
• 單元扭曲會嚴重影響計算精度，應確保單元形狀盡可能規則。
• 對于高精度分析，特別是應力集中區域，可以考慮使用二次完全積分單元，但需要權衡計算成本和精度需求。

3.3 線性縮減積分單元（默認設置）

理論基礎：線性縮減積分單元在每個方向上比完全積分單元少用一個積分點，通常在單元中心有一個積分點，屬于常應力單元。這種積分方式減少了計算量，但可能導致沙漏現象。

適用場景：線性縮減積分單元適用于各種結構分析，特別是存在彎曲載荷的結構。它們對位移求解結果較精確，能較好反映結構實際位移情況，適用于以位移為關鍵分析目標的場景。此外，在網格發生扭曲變形（如四邊形單元角度大幅偏離 90°）時，分析精度受影響較小，在處理復雜幾何模型時更具優勢。

優缺點分析：

• 優點：計算效率高，計算時間少；對位移求解結果較精確；在彎曲載荷下不易發生剪切自鎖；網格扭曲對計算精度影響較小。
• 缺點：積分點應力較精確，但經外插值和平均得到的節點應力精度差；存在沙漏數值問題，使單元表現過于柔軟；不適用于以應力集中部位節點應力為分析目標的情況。

使用注意事項：

• 在模擬承受彎曲載荷結構時，沿厚度方向至少劃分四個單元，以減輕沙漏效應影響。
• 對于以應力集中部位節點應力為分析目標的情況，不建議使用線性縮減積分單元。
• 可通過沙漏控制選項（如增強沙漏控制、剛度沙漏控制等）減輕沙漏問題，但應謹慎使用，過度控制可能導致單元過剛。
• 在 Abaqus/Explicit 中，線性縮減積分單元是動態分析的首選，因為它們具有集中質量公式，模擬應力波的效果優于二次單元所采用的一致質量公式。

3.4 二次縮減積分單元

理論基礎：二次縮減積分單元在二次完全積分單元的基礎上減少積分點數量，通常在每個方向上使用較少的高斯積分點，從而減少計算量。這種積分方式繼承了線性縮減積分單元的優點，同時對沙漏和剪切自鎖問題不敏感。

適用場景：二次縮減積分單元適用于對精度要求高但又需要控制計算成本的分析場景。它們對剪切自鎖和沙漏都不敏感，適用于各種復雜受力工況。在不需要模擬非常大的應變或進行復雜接觸條件變化的問題時，應優先采用二次縮減積分單元。

優缺點分析：

• 優點：位移結果精確；節點應力精度雖低于完全積分單元，但仍能滿足大多數工程需求；計算時間相對二次完全積分單元較少；對剪切自鎖和沙漏都不敏感；在復雜應力狀態下表現良好。
• 缺點：不能用于接觸分析；不適用于大應變問題，在材料發生大變形時難以保證計算精度；節點應力精度低于二次完全積分單元。

使用注意事項：

• 不適用于接觸分析場景，接觸問題中應選擇其他單元類型。
• 對于大應變分析，特別是材料發生大變形的情況，二次縮減積分單元可能無法提供足夠的精度，應考慮使用線性縮減積分單元或其他適合大變形的單元。
• 在應力集中區域，如果使用二次縮減積分單元，應確保網格足夠精細，以捕捉應力梯度變化。
• 對于需要精確節點應力結果的情況，尤其是應力集中區域的節點應力，二次縮減積分單元可能不如二次完全積分單元精確，需謹慎使用。

3.5 非協調模式單元

理論基礎：非協調模式單元在標準線性單元的基礎上引入了非協調位移模式，這些附加模式允許單元在不連續的情況下仍能保持協調性，從而克服了剪切自鎖問題。非協調模式單元通過增強變形梯度的方式，使單元交界處不會出現重疊或開洞問題，便于拓展到非線性、有限應變位移分析中。

適用場景：非協調模式單元適用于主要承受彎曲載荷的結構，特別是那些單元扭曲較小的區域。在彎曲問題中，厚度方向只需很少單元就能達到與二次單元相當的結果，但計算成本明顯降低。它們特別適合于模擬薄板、薄殼等以彎曲為主的結構。

優缺點分析：

• 優點：成功克服了剪切自鎖問題；在單元扭曲程度較小時，位移和應力結果精確；彎曲問題中厚度方向只需少量單元就能達到與二次單元相當的效果，大幅降低計算成本；能夠適應一定程度的單元扭曲。
• 缺點：若關注部位的單元扭曲程度大，尤其是出現交錯扭曲時，分析精度會下降；不適用于嚴重扭曲的網格；在某些情況下，特別是大變形分析中，可能需要更精細的網格劃分。

使用注意事項：

• 若在模型中采用非協調單元，應使網格扭曲減至最小，避免在扭曲嚴重的區域使用。
• 非協調單元在 Abaqus/Standard 中可用，但在 Abaqus/Explicit 中不可用，應根據求解器選擇合適的單元類型。
• 在彎曲問題中，非協調單元沿厚度方向只需 1-2 個單元就能獲得較好結果，但在關鍵區域仍建議進行網格細化以提高精度。
• 對于高精度分析，特別是需要精確應力結果的區域，即使使用非協調單元，也應確保網格質量和密度足夠。

3.6 四面體單元

理論基礎：三角形 (Tri) 和四面體 (Tet) 單元是基于線性或二次插值函數的單元類型，它們能夠適應任意幾何形狀，特別適合于復雜幾何模型的網格劃分。修正的二次 Tri 和 Tet 單元采用修正的二次插值函數，針對三角形和四面體單元特性進行了優化。

適用場景：

• 線性 Tri 和 Tet 單元：精度較差，不建議在關心部位及附近區域使用，但計算時間少，可用于模型中不重要的區域。
• 二次 Tri 和 Tet 單元：精度好，用于模擬任意幾何形狀，但計算時間較多，不適用于接觸分析及大變形大應變場景。
• 修正的二次 Tri 和 Tet 單元：適應于接觸分析和大變形大應變場景，在 Abaqus/Explicit 中應選擇修正的 Tet 單元 C3D10M，在 Abaqus/Standard 中可以選擇 C3D10，但如果有大的塑性變形或模型中存在接觸，也應選擇修正的 Tet 單元 C3D10M。

優缺點分析：

• 優點：能夠適應任意復雜幾何形狀；網格生成相對容易；修正的二次單元能夠處理接觸和大變形問題。
• 缺點：線性單元精度較差；二次單元計算成本高；非修正的二次單元不適用于接觸分析和大變形場景；同等數量單元下，四面體單元的應力結果通常不如六面體單元精確。

使用注意事項：

• 優先選擇四邊形 (Quad) 或六面體 (Hex) 單元，盡量避免使用 Tri 或 Tet 單元，尤其是在線性形式下。
• 如果必須使用 Tet 單元，應選擇二次或修正的二次單元，避免使用線性 Tet 單元。
• 在 Abaqus/Explicit 中，應選擇修正的 Tet 單元 C3D10M，以適應接觸分析和大變形場景。
• 當使用自由網格劃分技術時，Tet 單元的類型應選擇二次單元，在 Abaqus/Explicit 中選擇修正的 Tet 單元 C3D10M，在 Abaqus/Standard 中可以選擇 C3D10，但如果有大的塑性變形或接觸，也應選擇修正的 Tet 單元 C3D10M。
• 對于重要區域或需要高精度結果的區域，應盡量避免使用 Tet 單元，或在使用后通過網格細化和結果驗證確保計算精度。

3.7 雜交單元

理論基礎：雜交單元是一種特殊的單元類型，其特點是在單元內部使用不同的插值函數來描述位移和壓力場，從而能夠更準確地模擬不可壓縮材料的行為。雜交單元通過引入額外的自由度來表示壓力或體積變形，從而避免了常規單元在模擬不可壓縮材料時出現的體積自鎖問題。

適用場景：雜交單元主要用于模擬不可壓縮材料（如泊松比為 0.5 的橡膠）或近似不可壓縮材料（泊松比大于 0.475）。除了平面應力問題外，普通單元在模擬此類材料響應時，單元中應力不確定，而雜交單元可有效解決該問題。

優缺點分析：

• 優點：能夠準確模擬不可壓縮材料的行為；避免了體積自鎖問題；在處理幾乎不可壓縮材料時表現出色；提供更準確的應力和應變結果。
• 缺點：計算成本高于常規單元；需要更多的計算資源；不適用于可壓縮材料；在某些情況下，特別是非均勻材料或復雜幾何形狀中，可能需要更精細的網格劃分。

使用注意事項：

• 僅在 Abaqus/Standard 中可用，Abaqus/Explicit 中不支持雜交單元。
• 僅在材料為完全不可壓縮或近似不可壓縮時使用，對于普通可壓縮材料（如大多數金屬），常規單元通常更有效。
• 雜交單元需要適當的網格劃分，避免過度扭曲或不規則形狀，以確保計算精度。
• 在使用雜交單元時，應確保材料模型正確定義了不可壓縮特性，如橡膠材料通常使用超彈性模型并設置泊松比為 0.5。
• 對于高精度分析，特別是涉及不可壓縮材料的復雜結構，雜交單元是首選，但需要進行充分的網格敏感性分析以確保結果的可靠性。

4 實體單元使用最佳實踐

4.1 網格劃分最佳實踐

1. 優先使用六面體單元：三維情況應盡可能采用塊狀單元（六面體）。當幾何形狀復雜時，完全采用塊體單元構造網格會很困難，因此可能有必要采用楔形和四面體單元，但應盡量少用，并遠離需要精確求解的區域。
2. 合理控制單元尺寸：網格密度應根據模型的幾何特征和應力集中區域相匹配。在關鍵區域（如應力集中處、幾何突變處）應進行網格細化，而在次要區域可以使用較粗的網格，以平衡計算精度和效率。
3. 避免過度扭曲：單元扭曲會嚴重影響計算精度，應盡量避免。對于線性單元，單元角度應盡量接近 90°，避免出現尖銳角度或過度拉伸的形狀。對于二次單元，雖然對扭曲的容忍度較高，但仍應保持合理的單元形狀。
4. 使用結構化網格：對于規則幾何區域，應采用結構化網格劃分技術，生成高質量的六面體網格。結構化網格不僅質量高，而且計算效率更高，能夠提供更準確的結果。
5. 混合單元類型的使用：當三維實體幾何形狀復雜，無法用結構化或掃略網格劃分技術得到六面體單元網格時，可以對模型中不重要部分采用自由網格劃分技術生成四面體單元網格，而對關注部分用結構化或掃略網格劃分技術生成六面體單元網格。生成此類網格時，Abaqus 會提示將生成非協調網格，并在不同單元類型交界處自動創建 Tie 綁定約束，以保證模型的連續性。

4.2 單元選擇驗證方法

1. 網格敏感性分析：在確定最終單元類型和網格密度之前，應進行網格敏感性分析，比較不同網格密度和單元類型下的計算結果差異。如果結果對網格密度敏感，說明需要更精細的網格或更高階的單元。
2. 單元類型對比分析：對于關鍵區域或重要分析結果，建議進行不同單元類型的對比分析，觀察分析結果對單元類型的敏感性，并最終選擇合適的單元類型。例如，在應力集中區域，可以比較二次完全積分單元和二次縮減積分單元的結果差異。
3. 參考解析解或實驗數據：如果可能，將數值模擬結果與解析解或實驗數據進行比較，驗證單元選擇和網格劃分的合理性。這對于高精度分析尤為重要。
4. 使用子模型技術：對于復雜模型，可以使用 Abaqus 的子模型技術，先進行整體模型的粗網格分析，然后在關注區域創建子模型進行精細網格分析，從而在保證精度的同時控制計算成本。

4.3 應力結果評估與解釋

1. 應力不連續問題：在不同類型單元網格的交界處，即便單元角部節點重合，仍可能出現應力不連續的情況。而且，交界處應力有可能大幅增大。這是因為不同類型單元的插值方式、精度等存在差異，導致在連接區域應力傳遞不順暢。
2. 線性與二次單元混合使用問題：當在同一實體中混合使用線性和二次單元時，同樣會面臨類似應力不連續和應力增大的問題。因為線性單元和二次單元對結構變形和應力的描述能力不同，二者銜接處易產生計算差異。
3. 節點應力與積分點應力：不同單元類型提供的應力結果位置和精度不同。一般來說，積分點應力（高斯點應力）比節點應力更準確，但節點應力更便于結果展示和比較。對于線性縮減積分單元，節點應力是通過積分點應力外插和平均得到的，精度可能較低。
4. 結果驗證策略：在混合使用不同類型單元時，應使交界處遠離模型中重點關注的區域，減少應力不連續等問題對關鍵部位分析結果的影響。同時，要仔細檢查分析結果是否正確，通過對比、后處理等方式，評估交界處對整體結果的影響程度，確保分析的可靠性。

結論

  綜上所述，Abaqus 實體單元的選擇需以工程問題特性為核心，綜合考量幾何形狀、載荷類型、材料屬性、求解器特性及精度需求等多維度因素。線性單元與二次單元的取舍需平衡計算效率與精度，積分方式的選擇則需規避自鎖、沙漏等數值問題，而雜交單元、非協調單元等特殊類型的應用需嚴格匹配不可壓縮材料、彎曲主導等場景。

  實踐中，應優先采用六面體單元并控制網格質量，針對應力集中、接觸分析、大變形等特定場景靈活調整單元類型，同時結合求解器差異優化選擇策略。本文所梳理的分類體系與選擇框架，可為工程仿真中實體單元的合理應用提供理論支撐與實踐指導，后續研究可進一步結合具體工程案例開展單元性能的量化對比分析，以深化對復雜場景下單元選擇機制的理解。

完

更多精彩，關注建源學堂！