abaqus殼單元劃分的案例
【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
此外，工程中的復雜結構常涉及不規則網格劃分，擬協調元的網格適應性使其成為理想選擇。
非線性擬協調固體殼單元的應用
非線性擬協調固體殼單元憑借其高精度、高效率及良好的適應性，在多個工程領域和學術研究中展現出廣泛的應用前景，主要包括以下幾個方面：
（一）幾何非線性問題分析
大變形薄板殼結構
在薄板的大撓度彎曲、薄殼的失穩分析中，非線性擬協調固體殼單元能準確捕捉結構的幾何非線性響應。例如，對固支方板在均布載荷作用下的大變形分析（后期推文介紹，敬請期待！），單元通過共旋坐標法分離剛體運動與彈性變形，結合 von Karman 非線性板理論，可精確模擬載荷 - 位移曲線中的 “階躍” 現象。即使在粗網格（4×4×2）下，單元計算結果與解析解的誤差仍小于 5%，顯著優于傳統 C3D8R/Solid45 單元。
將擬協調單元CSS8與 ANSYS 的 Solsh190、ABAQUS 的 SC8R進行對比，從精度、效率、穩定性三方面評估優勢。例如，在 薄膜分析中，CSS8 單元在 2×2×2 網格下的位移誤差為 5.2%，優于 Solsh190 的 17.3%，SC8R的25%。
復雜曲面殼結構
對于含初始曲率的殼結構（如半球殼、圓柱殼），單元能有效避免曲率厚度鎖定，準確描述雙曲率變形。在頂部開孔半球殼的大變形分析中，八節點擬協調固體殼單元（CSS8）在 16×16×2 網格下的位移計算誤差僅為 3.2%，而傳統殼單元（如 Abaqus C3D8）誤差高達 15% 以上。
結構失穩與后屈曲分析
在淺殼結構的失穩分析中，單元結合弧長法可追蹤完整的后屈曲路徑，準確預測臨界載荷和失穩模式。
展開 ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
hypermesh tcl腳本 自動劃分2D殼單元網格 ￥6.66
自己編寫的基于batchmesh自動劃分2D單元tcl腳本，需提前準備好criteria和param文件；
通過tcl腳本選擇component抽取中面并創建殼單元，然后輸入材料、自動識別厚度并重命名。
【JY】Abaqus殼單元概述與應用（一）
殼法線還定義了施加在單元上正壓力載荷的方向，并可以在 Abaqus/Post 中顯示。在創建殼單元時，應注意法線方向是否正確，否則可能導致接觸定義錯誤或載荷方向相反。
殼截面定義：殼單元的橫截面特性可以由沿厚度方向的數值積分確定（*SHELL SECTION），或在分析開始時應用計算的橫截面剛度（*SHELL GENERAL SECTION）。*SHELL GENERAL SECTION 非常有效，但僅用于線性材料，而SHELL SECTION 可用于線性和非線性材料。數值積分在沿殼厚度方向的一系列積分點上進行，這些積分點就是單元變量可以被輸出的位置。最外層的積分點位于殼單元的表面。
復合材料鋪層定義：對于復合材料層合殼，殼截面定義變得更加復雜，需要定義各鋪層的材料、厚度、方向和順序。在 Abaqus/CAE 中，可以通過 "Composite Layup" 工具方便地定義復合材料鋪層。常規殼方法將層合板建模為一系列薄層，同時忽略厚度方向的復雜剪切變形；連續殼方法則使用三維網格離散化，能夠更準確地描述厚度方向的剪切變形。
4.2 殼單元網格劃分最佳實踐
網格密度控制：殼單元的網格密度應根據結構的幾何特征、載荷條件和分析目標確定。在應力集中區域、幾何突變處或載荷復雜的區域，應使用更精細的網格。對于以彎曲為主的結構，沿厚度方向應至少劃分 2-3 個單元，以捕捉厚度方向的應力梯度。對于高精度分析，可能需要更多的單元。
單元形狀優化：殼單元應盡量保持規則的形狀，如接近正方形的四邊形單元。應避免使用高度扭曲或長寬比過大的單元，這些單元可能導致計算精度下降或收斂困難。對于三角形單元，應盡量使其接近等邊三角形，避免出現尖銳角度。
展開 ABAQUS殼單元輪胎模型仿真案例
圖1子午線輪胎結構分布圖
目前不少工作對輪胎的建模通常采用軸對稱單元，在充氣后通過修改INP文件將輪胎置于路面上令其滾動觀察響應，三維實體單元的輪胎建模方法可見ABAQUS三維輪胎充氣滾動案例_輪胎仿真 ABAQUS-技術鄰，本文介紹一種采用殼單元對輪胎進行建模的方法，相比三維實體，殼單元的計算速度更快，建模方式更簡便，但相對的殼單元的計算精度與模擬的準確性上有時會不太理想。
1 建模
輪胎模擬的一個難點是其內部加強層的模擬。通常的軸對稱單元與實體單元采用rebar layer的方式進行建模，并采用內嵌區域的方法將加強層嵌入到輪胎主體中。但殼模型無法作為主體區域，因此本研究采用復合層的截面定義方式對機輪殼模型進行截面賦予，對機輪不同區域定義不同的復合層數及相應的厚度與材料屬性。如鋼線圈區域，共指派了十一層，并按照橡膠-內面層-橡膠-鋼線圈-橡膠-鋼線圈-橡膠-鋼線圈-橡膠-內面層-橡膠的排布方式賦予了該區域相應的截面屬性，每一層的厚度與旋轉角均與輪胎本身的定義保持一致，鋼線圈區域的復合層定義與層堆疊繪圖見表1與圖2所示。機輪其余區域的截面定義方式與鋼線圈類似。
展開 ABAQUS實體、殼、梁單元的軸力、剪力、彎矩的提取方式及準確性驗證 ￥8
在ABAQUS中，對結構或者構件進行受力分析除了分析應力云圖之外，通常還需要對部件的軸力、剪力或彎矩的變化趨勢進行分析。本帖基于以下的實體solid、殼shell、梁/beam（truss）模型，分別提取這三類模型的軸力、剪力、彎矩，并與理論計算相結合，驗證提取結果的準確性，并解釋相應有限元的計算原理。
計算模型
梁單元計算結果
實體單元計算結果
殼單元計算結果
帖子內容概況
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喵星人教你學會鋼筋混凝土殼單元的前處理與后處理
ABAQUS中的殼單元大家通常用于模擬鋼板等鋼結構，對于混凝土板殼，新手可能對內部的配筋方式，以及前后處理方法可能存在各種問題。實際上，ABAQUS提供了鋼筋混凝土板配筋的接口，這種“寫入式”而不進行直接建模的方法通常比較冷門且后處理相對不主流。今天喵星人就通過一個教程教你學會鋼筋混凝土殼單元的前處理與后處理。
0.前提
使用板殼單元的有限元模擬必須有兩個前提：
1、板殼力學及殼單元通常應用于一個方向尺寸遠小于另外兩個方向（通常不超過1/5）的結構。
喵星人點評：大家總有一個誤區，總覺得實體單元的精度最高，實則不然。對于板殼結構，由于其采用了Kirchhoff板假定，在此情況下相比實體單元，殼單元形函數更加逼近實際結構，其計算精度與計算代價均優于采用實體單元。
2、由于采用Kirchhoff板假定，即忽略混凝土板中鋼筋的粘結滑移行為，因此在精細化的鋼筋混凝土滯回模型中通常不再適用。
1、前處理
1.1 縱橫方向與局部坐標系
配筋的板殼單元，尤其是兩個平面方向差異配筋的板殼單元，必須指定坐標系，且喵星人建議使用局部坐標系。這是為了避免在裝配件中因旋轉導致整體坐標系的變換。本案例中的坐標系指派如圖所示。需要注意的是，鋼筋縱橫方向與局部坐標系方向直接掛鉤。
1.2 配筋面積/間距/方向
殼單元的配筋方法需在“編輯截面”中完成，不能直接建立線單元鋼筋。采用“寫入式”的建模方法，如下圖所示。
其實這種方法很像設計軟件中的操作，即通過加勁的方式考慮配筋混凝土。
展開 
Abaqus中考慮橫向剪切的復合材料厚殼單元分析
在前面的兩篇文章中，已經對Abaqus復合材料殼單元分析模型的傳統建模方法和快捷建模方法做了詳細的介紹。熟悉Abaqus復合材料分析的人都知道，在采用二維Lamina復合材料模型配合殼單元進行分析時，材料參數中除了輸入兩個方向模量E1,E2，面內泊松比及面內剪切模量G12之外，還要額外輸入兩個面外剪切模量G13和G23，如下圖所示。
這里的面外剪切模量G13和G23就是用來考慮橫向剪切變形的。
一般，針對薄板結構（跨厚比大于20），通常做以下假設（Kirchhoof假設）：
(1)平行于中面 的各層互不擠壓：即垂直于中面法向的正應力很小，可以忽略。
(2)直法線假設：變形前垂直于中面的直線段，在變形后仍保持是直線，且仍垂直于變形后的中面。
(3)撓度沿板厚度方向的變化可以忽略，即統一厚度各點的撓度都 等于中面的撓度
(4)板的中面無伸縮和剪切變形
根據上述假設，板的橫向變形為零，相當于垂直于中面的各個面內剪切模量無窮大。薄板理論的假設在求解薄板問題時，精度足以滿足工程計算要求。
但對于中厚板或者厚板、集中力作用點附近、薄板邊界周圍以及開孔周圍，上述理論將不再適用，誤差大甚至會導致錯誤的結果，因此為了解決此類問題，便有了考慮剪切變形的中厚板理論。
那么在Abaqus分析中怎么考慮橫向剪切變形的影響呢？Abaqus默認的復合材料模型定義及截面屬性定義中是已經考慮了橫向剪切的，軟件會自動計算橫向剪切剛度。
而薄板、中厚板的區分在于單元類型，如下圖所示，S8R5為薄殼單元的一種，S8R為厚殼單元的一種，在設置單元屬性時會有明確的說明：
現在，測試一下薄殼與厚殼計算結果的差異有多大。
展開 [軟件使用]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，如下圖所示。
展開 