abaqus 厚殼單元的案例
Abaqus中考慮橫向剪切的復合材料厚殼單元分析
在前面的兩篇文章中，已經對Abaqus復合材料殼單元分析模型的傳統建模方法和快捷建模方法做了詳細的介紹。熟悉Abaqus復合材料分析的人都知道，在采用二維Lamina復合材料模型配合殼單元進行分析時，材料參數中除了輸入兩個方向模量E1,E2，面內泊松比及面內剪切模量G12之外，還要額外輸入兩個面外剪切模量G13和G23，如下圖所示。
這里的面外剪切模量G13和G23就是用來考慮橫向剪切變形的。
一般，針對薄板結構（跨厚比大于20），通常做以下假設（Kirchhoof假設）：
(1)平行于中面 的各層互不擠壓：即垂直于中面法向的正應力很小，可以忽略。
(2)直法線假設：變形前垂直于中面的直線段，在變形后仍保持是直線，且仍垂直于變形后的中面。
(3)撓度沿板厚度方向的變化可以忽略，即統一厚度各點的撓度都 等于中面的撓度
(4)板的中面無伸縮和剪切變形
根據上述假設，板的橫向變形為零，相當于垂直于中面的各個面內剪切模量無窮大。薄板理論的假設在求解薄板問題時，精度足以滿足工程計算要求。
但對于中厚板或者厚板、集中力作用點附近、薄板邊界周圍以及開孔周圍，上述理論將不再適用，誤差大甚至會導致錯誤的結果，因此為了解決此類問題，便有了考慮剪切變形的中厚板理論。
那么在Abaqus分析中怎么考慮橫向剪切變形的影響呢？Abaqus默認的復合材料模型定義及截面屬性定義中是已經考慮了橫向剪切的，軟件會自動計算橫向剪切剛度。
而薄板、中厚板的區分在于單元類型，如下圖所示，S8R5為薄殼單元的一種，S8R為厚殼單元的一種，在設置單元屬性時會有明確的說明：
現在，測試一下薄殼與厚殼計算結果的差異有多大。
展開 ANSYS中薄殼厚殼分類及單元特性
板殼分類
按板面內特征尺寸與厚度之比劃分：
當 L/h < (5~8) 時為厚板，應采用實體單元。
當 (5~8) < L/h < (80~100) 時為薄板，可選 2D 實體或殼單元
當 L/h > (80~100) 時為薄膜，可采用薄膜單元。
殼類結構按曲率半徑與殼厚度之比劃分：
當 R/h >= 20 時為薄殼結構，可選擇薄殼單元。
當 6 < R/h < 20 時為中厚殼結構，選擇中厚殼單元。
當 R/h <= 6 時為厚殼結構。
上述各式中 h 為板殼厚度， L 為平板面內特征尺度，R 為殼體中面的曲率半徑。
2. 薄板理論的基本假定
薄板所受外力有如下三種情況：
① 外力為作用于中面內的面內荷載。彈性力學平面應力問題。
② 外力為垂直于中面的側向荷載。薄板彎曲問題。
③ 面內荷載與側向荷載共同作用。
所謂薄板理論即板的厚度遠小于中面的最小尺寸，而撓度又遠小于板厚的情況，也稱為古典薄板理論。
薄板通常采用 Kirchhoff-Love 基本假定：
① 平行于板中面的各層互不擠壓，即 σz = 0。
② 直法線假定：該假定忽略了剪應力和所引起的剪切變形，且認為板彎曲時沿板厚方向各點的撓度相等。
③ 中面內各點都無平行于中面的位移。
薄板小撓度理論在板的邊界附近、開孔板、復合材料板等情況中，其結果不夠精確。
3. 中厚板理論的基本假定
考慮橫向剪切變形的板理論，一般稱為中厚板理論或 Reissner（瑞斯納）理論。
展開 【JY】Abaqus“殼”單元概述與應用（二）——固體殼單元
寫在前文
在有限元分析中，單元類型的選擇對計算結果的精度和效率有著決定性影響，尤其對于復合材料結構和薄壁結構的分析更是如此。
Abaqus 作為主流的有限元分析軟件，提供了多種固體殼單元類型以滿足不同工程需求。連續實體殼單元 (CSS8)、非協調元 (C3D8I) 和連續殼單元 (SC8R) 是 Abaqus 中常用于復合材料和薄壁結構分析的三種單元類型，各自具有獨特的理論基礎和適用場景。
相關閱讀：
【JY】Abaqus殼單元概述與應用（一）
除了上述采用類實體單元的“殼”單元外，還有完全的殼單元，如S4R 單元，是 Abaqus 中最常用的常規殼單元之一，為 4 節點減縮積分殼單元，基于經典殼理論，適用于各類薄壁結構的線性與非線性分析，尤其在大變形和接觸問題中表現穩定，將該單元作為對比基準，對上述實體類“殼”單元進行對比分析。
本文旨在對這三種單元類型進行深入比較研究，從理論基礎、自由度、材料本構、積分方案、閉鎖敏感性、計算成本等多個維度展開分析，為工程實踐中的單元選擇提供參考。特別是針對復合材料分析、金屬薄壁結構模擬以及混合建模等應用場景，探討這三種單元的適用性差異，并分析它們在幾何非線性情況下的計算成本和精度表現。
單元類型基本原理與特點
2.1 連續實體殼單元 (CSS8)
連續實體殼單元 (CSS8) 是一種介于 C3D8I (非協調元) 和 SC8R (連續殼單元) 之間的特殊一階單元，由 Vu-Quoc 和 Tan 于 2003 年提出，后集成于 SIMULIA 2017 及以后的版本。它是一種三維單元，具有以下基本特點：
幾何與自由度：CSS8 為 8 節點六面體單元，僅有位移自由度 (無轉動自由度，與實體單元一致)，與實體單元混合建模時易于處理連接過渡。
展開 【JY】Abaqus“殼”單元概述與應用（三）——非線性擬協調固體連續殼單元CSS8
【相關閱讀】
【JY】Abaqus殼單元概述與應用（一）
【JY】Abaqus 三維應力單元解析、選擇與應用指南
【JY】Abaqus“殼”單元概述與應用（二）——固體殼單元
傳統固體殼單元在處理幾何非線性、材料非線性及復雜邊界條件時，存在諸多難以克服的缺陷，這促使研究者探索新的單元構造方法。非線性擬協調固體殼單元的提出，正是為了突破這些局限，其研究動因主要源于以下幾方面：
（一）傳統固體單元的固有缺陷
自鎖現象普遍存在
傳統固體單元（如C3D8R）在模擬薄板殼結構時，易出現剪切自鎖、薄膜自鎖、體積自鎖等問題。剪切自鎖源于單元位移插值無法準確表征純彎曲狀態下的零剪切應變，導致計算結果剛度偏高；薄膜自鎖則因低階形函數無法捕捉不可伸縮彎曲模式下的面內應變分布，使位移被低估；體積自鎖多見于近不可壓縮材料分析，由于單元無法準確描述等體積運動，導致體積變化被過度約束。這些自鎖現象嚴重影響計算精度，尤其是在粗網格或大長高比結構中表現更為突出。
計算效率與精度的矛盾
為克服自鎖問題，需要采用增強假設應變法（EAS）、假設自然應變法（ANS）或雜交應力法等，這些方法往往需要引入額外的內部參數或復雜的數值積分，使得單元列式復雜、相對殼單元計算成本增加。
幾何非線性處理的局限性
現有非線性固體殼單元多基于連續體變形梯度的極分解處理幾何非線性，該方法不僅計算量大，且在 Cartesian 坐標系下難以保證旋轉描述的準確性。在大變形、大轉動問題中，極分解可能導致切線剛度矩陣奇異，影響迭代收斂性。此外，傳統單元在處理不規則網格或畸變網格（如C3D8I）時，精度衰減明顯，難以滿足工程對復雜結構分析的需求。
展開 
【JY】Abaqus殼單元概述與應用（一）
連續殼單元不適用于需要考慮厚度方向法向應力的情況，此時應使用連續實體方法。
3.4 特殊殼單元類型
Abaqus 還提供了一些特殊類型的殼單元，適用于特定的應用場景：
軸對稱殼單元：如 SAX1、SAX2、SAX2T 等，用于模擬軸對稱殼結構，如圓柱殼、圓錐殼等。這些單元能夠顯著減少模型規模，提高計算效率，特別適用于具有軸對稱幾何和載荷的問題。
薄殼單元：如 S4R5、S8R5 等，強化了基爾霍夫條件，適用于非常薄的殼結構。這些單元假設垂直于殼中面的平面在變形后保持垂直，忽略了橫向剪切變形，適用于厚度與其他尺寸之比小于 1/30 的薄殼。
厚殼單元：如 S8R、S8RT 等，考慮了橫向剪切變形，適用于厚度與其他尺寸之比在 1/20 至 1/10 之間的厚殼結構。這些單元能夠更準確地模擬厚殼的剪切變形和應力分布。
可擴展殼單元：如 S6、S9 等，允許用戶定義自定義的殼截面形狀和材料分布，適用于特殊截面或非均質材料的殼結構。
使用注意事項：
軸對稱殼單元僅適用于幾何形狀和載荷均軸對稱的問題，否則會導致嚴重的計算誤差。
薄殼單元和厚殼單元的選擇應基于結構的厚跨比，通常厚跨比小于 1/30 時使用薄殼單元，介于 1/30 和 1/10 之間時使用一般殼單元，大于 1/10 時考慮使用厚殼單元或實體單元。
對于復合材料層合板，特殊殼單元（如連續殼單元）通常比普通殼單元更合適，因為它們能夠更好地模擬層間剪切變形和各向異性材料特性。
在使用任何特殊殼單元之前，應充分了解其理論基礎和適用范圍，確保其適合特定的工程問題。
4 殼單元使用注意事項
4.1 殼單元方向與截面定義
殼單元法線方向：殼單元法線方向決定了單元的正和負表面，為了正確地定義接觸和解釋輸出數據，必須知道其對應的是哪個面。
展開 ABAQUS殼單元輪胎模型仿真案例
圖1子午線輪胎結構分布圖
目前不少工作對輪胎的建模通常采用軸對稱單元，在充氣后通過修改INP文件將輪胎置于路面上令其滾動觀察響應，三維實體單元的輪胎建模方法可見ABAQUS三維輪胎充氣滾動案例_輪胎仿真 ABAQUS-技術鄰，本文介紹一種采用殼單元對輪胎進行建模的方法，相比三維實體，殼單元的計算速度更快，建模方式更簡便，但相對的殼單元的計算精度與模擬的準確性上有時會不太理想。
1 建模
輪胎模擬的一個難點是其內部加強層的模擬。通常的軸對稱單元與實體單元采用rebar layer的方式進行建模，并采用內嵌區域的方法將加強層嵌入到輪胎主體中。但殼模型無法作為主體區域，因此本研究采用復合層的截面定義方式對機輪殼模型進行截面賦予，對機輪不同區域定義不同的復合層數及相應的厚度與材料屬性。如鋼線圈區域，共指派了十一層，并按照橡膠-內面層-橡膠-鋼線圈-橡膠-鋼線圈-橡膠-鋼線圈-橡膠-內面層-橡膠的排布方式賦予了該區域相應的截面屬性，每一層的厚度與旋轉角均與輪胎本身的定義保持一致，鋼線圈區域的復合層定義與層堆疊繪圖見表1與圖2所示。機輪其余區域的截面定義方式與鋼線圈類似。
展開 Abaqus如何使用殼單元建模分析
Abaqus如何使用殼單元建模分析
前幾天突然需要用到Abaqus的殼單元，本以為會和ANSYS似的，直接修改單元類型即可，自己試了試發現完全不是這回事兒。沒辦法網上查了查，居然沒有Abaqus殼單元方面的實際操作，大多都是說殼單元的結果輸出之類的，看來筆者真是知道的太少，無奈之下還是只能自己試。
Abaqus的殼單元做分析在單元類型里面無法直接定義，而是通過材料屬性進行賦予的，但是材料屬性賦予的時候還得和模型的類型有關。下面大致說一下Abaqus用殼單元做分析的過程。
如圖1所示，建立Part時需要指定part類型，筆者想建一個平面，有厚度，用殼單元賦予厚度。那么Modeling Space必須是3D，如果選了2D那么就無法賦予殼單元屬性，雖然建模的時候確實只是建一個平面，但是還是3D，這個理解起來就只能是考慮有厚度，殼單元模型代表的還是3維實體模型。這個和ANSYS的概念還真不一樣，ANSYS沒這么繞。
圖1
之后建立了一個平面矩形，進入材料模塊。添加一個材料屬性后，需要創建一個Section，如圖2所示。
圖2
Section的Category指定為Shell，點擊Continue后，如圖3.
圖3
圖3中的Value指定殼單元的厚度，之后給模型賦予建立的Section，如圖4所示。
圖4
其中的Shell Offset下面有五個選項，這個意義很好理解，殼單元厚度的定義方式，中面底面頂面等。
再到Mesh模塊下面，即可發現有殼單元選項Shell，如圖5所示。
圖5
Abaqus的殼單元類型S4R(縮減積分單元)，還可以通過Quadratic指定為二次單元S8R。
再往后的過程就和其他一致，不作贅述。
展開 ABAQUS實體、殼、梁單元的軸力、剪力、彎矩的提取方式及準確性驗證 ￥8
在ABAQUS中，對結構或者構件進行受力分析除了分析應力云圖之外，通常還需要對部件的軸力、剪力或彎矩的變化趨勢進行分析。本帖基于以下的實體solid、殼shell、梁/beam（truss）模型，分別提取這三類模型的軸力、剪力、彎矩，并與理論計算相結合，驗證提取結果的準確性，并解釋相應有限元的計算原理。
計算模型
梁單元計算結果
實體單元計算結果
殼單元計算結果
帖子內容概況
ABAQUS喵星人教你學會鋼筋混凝土殼單元的前處理與后處理
ABAQUS中的殼單元大家通常用于模擬鋼板等鋼結構，對于混凝土板殼，新手可能對內部的配筋方式，以及前后處理方法可能存在各種問題。實際上，ABAQUS提供了鋼筋混凝土板配筋的接口，這種“寫入式”而不進行直接建模的方法通常比較冷門且后處理相對不主流。今天喵星人就通過一個教程教你學會鋼筋混凝土殼單元的前處理與后處理。
0.前提
使用板殼單元的有限元模擬必須有兩個前提：
1、板殼力學及殼單元通常應用于一個方向尺寸遠小于另外兩個方向（通常不超過1/5）的結構。
喵星人點評：大家總有一個誤區，總覺得實體單元的精度最高，實則不然。對于板殼結構，由于其采用了Kirchhoff板假定，在此情況下相比實體單元，殼單元形函數更加逼近實際結構，其計算精度與計算代價均優于采用實體單元。
2、由于采用Kirchhoff板假定，即忽略混凝土板中鋼筋的粘結滑移行為，因此在精細化的鋼筋混凝土滯回模型中通常不再適用。
1、前處理
1.1 縱橫方向與局部坐標系
配筋的板殼單元，尤其是兩個平面方向差異配筋的板殼單元，必須指定坐標系，且喵星人建議使用局部坐標系。這是為了避免在裝配件中因旋轉導致整體坐標系的變換。本案例中的坐標系指派如圖所示。需要注意的是，鋼筋縱橫方向與局部坐標系方向直接掛鉤。
1.2 配筋面積/間距/方向
殼單元的配筋方法需在“編輯截面”中完成，不能直接建立線單元鋼筋。采用“寫入式”的建模方法，如下圖所示。
其實這種方法很像設計軟件中的操作，即通過加勁的方式考慮配筋混凝土。
展開 [軟件使用]abaqus殼單元局部坐標系，你學會了嗎？
殼單元是一種結構單元，該結構一個方向的尺度(厚度)遠小于其它方向的尺度，并忽略沿厚度方向的應力。例如,壓力容器結構的壁厚小于典型整體結構尺寸的1/10,一般就可以用殼單元進行模擬。
在使用abaqus進行有限元分析的工作中，確定殼單元局部坐標系是一項重要的工作，其原因之一在于在abaqus中，殼單元的位移輸出基于整體坐標系，應力應變輸出基于局部坐標系，因此如果不能準確地確定殼單元的局部坐標系，在后處理查看計算結果時可能會無法準確理解計算結果。
通常情況下，殼單元的局部坐標系如下圖所示，其包含平面內的1，2軸和平面法線的n軸（3軸）。顯然，n軸由殼單元所在平面確定，但是其有兩種選擇，即由“殼內指向殼外”和由“殼外指向殼內”。
那么在abaqus中，殼單元的局部坐標系依據以下規則定義：
（1）對于一個3節點/4節點殼單元，按照右手定則，拇指指向即為n軸方向。
殼單元節點順序為1-2-4-3時的n軸方向。
（2）確定好n軸之后，接下來的1軸和2軸按照以下規則確定：
將整體坐標系的X軸投影到殼單元上，投影方向即為1軸。再按照右手定則，1-2-n軸形成右手坐標系，即右手拇指指向n軸時，其余4指的旋轉方向從1軸轉向2軸，具體圖解如下：右側為整體坐標系，左手為局部坐標系。
按照上述規則必然會存在一種特殊情況，即整體1軸與殼單元垂直，則此時整體1軸投影到殼單元上會是一個點，無法確定局部1軸方向，在這種情況下，abaqus采用整體3軸投影到殼單元上作為局部1軸方向。
以上就是殼單元局部坐標系的確定過程，下面以一個例子，來表明殼單元局部坐標系確定的具體作用。
展開 abaqus做用殼單元節段鋼箱拱時的幾何錯誤問題
部件1.png
路徑1.png
部件3.png
路徑3.png
截面.png
幾何診斷.png
因為拱是鋼箱做的，按照工程給的CAD圖逐個描點繪制路徑，然后采用同一截面掃掠，結果出現有的掃掠能出現正常部件，有的掃掠出現某些單元丟失的現象，經查找幾何診斷，發現可能是兩段的幾何沖突，請問這種問題該如何解決，有快速解決的辦法嗎？
我想過每一段單獨做出來然后合并，但是這樣工作量太大，因為該橋有四個拱，而且合并會導致尺寸和原來不一致，后面的工作也不再準確。
Abaqus復合材料殼單元建模—姊妹篇1：常規建模step-by-step
采用商業有限元軟件Abaqus進行復合材料結構建模時，一般有兩種建模方法：常規建模方法和Composite layup快速建模方法，主要差異在創建屬性、賦屬性和指定鋪層坐標系方面，常規建模方法和一般商業軟件類似，將創建材料、創建屬性、賦屬性和指定鋪層坐標系四個步驟分離，通用性較強，尤其是對于包含UMAT/VUMAT子程序開發的復合材料分析模型或者是三維實體單元顯式動力學分析模型，僅支持該類建模方法；Composite layup快速建模方法將創建屬性、賦屬性和指定鋪層坐標系三部分內容集成在一起，可一次性完成設置，效率較高。本文先從最基本的常規建模方法講起。
一般對于大尺寸復合材料結構，跨厚度比例大，滿足板殼理論的假設，采用殼單元就能獲得高的求解精度。殼單元計算效率高，結合二維損傷起始判據判據（Hashin, Tsai-W, Maxe, Maxs等）可以預測結構的危險區域和危險程度，另外，Abaqus自身還內嵌了二維Hashin的漸進損傷分析模型，采用Hashin失效判據去判斷損傷起始，損傷起始以后采用基于能量演化的連續退化準則對材料剛度進行退化。
Abaqus中常用的殼單元類型有S4、S4R、S8R等。以下介紹復合材料開孔板殼單元模型的建模步驟。
第1步：繪制幾何
在Part模塊下繪制幾何，幾何類型為3D-Deformable- Shell，草圖如下：
繪制完草圖后，退出草圖，得到開孔板的幾何模型，如下：
第2步：創建材料
與復合材料殼單元對應的是2D材料模型Lamina，將視圖切換至Property模塊，點擊創建材料按鈕，在跳出窗口中選擇Mechanical→Elasticity→Elastic選項，在材料類型下拉框中選擇Lamina，如下圖所示。
展開 Abaqus纖維復合材料蜂窩板三點彎曲仿真模型!采用了連續殼單元 ￥20
image_process=/format,webp" data-initial-src="https://img.jishulink.com/202505/attachment/e153f348ff064d30a136d2bf132cbf47.png" data-extentions-extra-ocr-id="f1e17ed722694393082635677f52cd9b">
</figure>
</figure><div contenteditable="false" width="100%">
Abaqus纖維復合材料蜂窩板三點彎曲仿真模型!采用了連續殼單元
</div><div contenteditable="false" width="100%">
有錄制整個建模操作視頻，可贈送復合材料層合板快速建模插件及蜂窩建模插件!
</div><div contenteditable="false" width="100%">
內附cae，inp文件及ODB文件，操作視頻
</div><p><br></p>
展開 Abaqus復合材料殼單元建模—姊妹篇2：layup快捷建模step-by-step
采用商業有限元軟件Abaqus進行復合材料結構建模時，一般有兩種建模方法：常規建模方法和Composite layup快速建模方法，主要差異在創建屬性、賦屬性和指定鋪層坐標系方面，常規建模方法和一般商業軟件類似，將創建材料、創建屬性、賦屬性和指定鋪層坐標系四個步驟分離，通用性較強，尤其是對于包含UMAT/VUMAT子程序開發的復合材料分析模型或者是三維實體單元顯式動力學分析模型，僅支持該類建模方法；Composite layup快速建模方法將創建屬性、賦屬性和指定鋪層坐標系三部分內容集成在一起，可一次性完成設置，效率較高。本文先從最基本的常規建模方法講起。
上一篇已經講解了最基本的常規建模方法，本篇將繼續介紹采用Composite layup實現快速建模，兩篇有明顯差異的地方用紅色字體進行了標注，以利于區分。同樣先介紹復合材料殼單元模型快速建模方式。
第1步：繪制幾何
在Part模塊下繪制幾何，幾何類型為3D-Deformable- Shell，草圖如下：
繪制完草圖后，退出草圖，得到開孔板的幾何模型，如下：
第2步：創建材料
與復合材料殼單元對應的是2D材料模型Lamina，將視圖切換至Property模塊，點擊創建材料按鈕，在跳出窗口中選擇Mechanical→Elasticity→Elastic選項，在材料類型下拉框中選擇Lamina，如下圖所示。
表格中的6個數據分別為縱向（沿纖維方向）彈性模量、橫向（垂直于纖維方向）彈性模量、面內泊松比以及三個方向的剪切模量。與其他商業有限元軟件不同的是，即使是對于二維材料模型，仍然需要輸入面外的剪切模量G13和G23，這兩項數據是用于定義殼的橫向剪切行為。
一個復合材料分析模型中可以包含多種材料模型，例如不同的鋪層采用不同的材料。
展開 abaqus算得到底準不準？殼單元幾何非線性驗證算例一
板殼單元是有限元中應用廣泛而又具有難度的單元類型，其在相當一段時間內是計算力學研究者廣泛研究的對象。常見的有K J bathe的MITC殼單元，龍馭球院士的廣義協調殼單元，Belytschko的belytschko-tsay殼單元等。殼單元在幾何非線性下的響應評估準確度和計算效率，是有限元軟件幾何非線性計算能力的重要體現。
本文演示一個殼單元幾何非線性驗證算例的abaqus操作。后續還可能采用其他軟件進行該模型及其他幾個殼單元幾何非線性的驗證對比。
展開 