abaqus梁單元類型的案例
ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個多層混凝土框架結構，一般除計算整體指標外，我們在計算具體荷載作用時（如風荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時可用殼單元模擬，主梁、次梁采用梁單元模擬，此時變為梁單元包含在殼面內的情況，當然此類情況是否需要考慮截面偏置，可根據具體工程而定。
對這中梁單元包含在殼單元面內的情況，只需要將梁單元與殼單元共用節點即可，而無須格外建立約束方程。
三、梁單元在殼單元內但不包含
此種情況為梁與殼位于同一面內，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區域兩種方式連接。剛性梁采用MPC184單元，剛性區域采用Cerig命令，具體使用方法下期文章討論。
展開 ANSYS各類型單元連接專題講解（四）之2D梁與2D實體單元剛接
前面文章主要講解了梁單元與其他類型單元鉸接的情況，從本篇文章開始，將主要講解梁單元與各類單元剛接的情況，而這也是我們日常工程中比較常見的一種連接方式。
首先從2D平面單元單元開始說起。
盡管現在的ANSYS版本已經摒棄了很古老的2D梁單元，改用Beam18x系列單元代替，但為究其連接方法，這類方面仍具有一定的講解價值，例如我們計算一榀框架的時候多數時候是采用2D平面單元的。
2D梁單元包括：beam3、beam23、beam54
2D實體單元：plane單元
一般來講，2D梁單元與2D實體單元剛接一般分為三種方法：
1）約束方程法；2）偽梁法；3）MPC法。
三種方法的連接原理無非是建立自由度之間的關系方程，但值得注意的是由于采用了局部區域的節點，因而在建立關系的局部區域內可能會有應力集中的情況，后處理當中應格外注意。
約束方程法后續講解3D梁單元連接時會詳細說明，此處簡單說下偽梁法與MPC法。
其實偽梁法與MPC法原理基本一樣，構造一個虛擬梁單元，使虛擬梁單元與外部梁單元剛接，虛擬梁單元與內部實體單元強制剛接，從而間接實現外部梁單元與實體單元的剛接效果。
使用偽梁法需注意的是，在建立虛擬梁單元時，虛擬梁單元應至少與實體單元的兩個節點相連，剛度可取為無窮大或者實際梁單元的10^5倍。
下面以一個小案例來演示。
如上圖所示，兩塊小鋼板中間靠一小鋼梁連接，小鋼梁上有均布荷載，尺寸如上所示，均以mm計，中間鋼梁所受均布荷載為10KN/m,采用ANSYS模擬該情況。
展開 ANSYS各類型單元連接專題講解（二）之桿與梁殼體單元的連接
前一篇文章主要介紹了單元之間連接的主要原則，今天開始主要從具體方面講解連接方法。
按照桿、梁、殼、實體的順序，先說說桿單元與各單元的連接方法。
那么什么時候需要用到桿單元與各種單元的連接呢？水哥稍微列舉下實際工程中需要考慮此類連接的例子。
案例一：工業廠房
此類結構一般橫向跨度較大，屋頂采用鋼結構形式，在具體模擬屋架時，此時各個桿件可看成鉸接，采用桿單元模擬。而下方框架柱則采用梁單元進行模擬，在相交部位則需要用到桿單元與梁單元的連接。
案例二：門廳鋼結構雨棚
在具體模擬該結構時，雨棚上方拉桿采用桿單元模擬，而下方的鋼梁采用梁單元模擬，混凝土框架柱可采用實體單元模擬。
一直以來，桿單元一般用于模擬桁架結構的時候比較多，其特點是桿件兩端不考慮承受彎矩作用，節點只有平動自由度，是所有單元中最為簡單的一種。
桿單元分為2D桿單元和3D桿單元，2D桿單元節點只有Ux和Uy兩個平動自由度，而3D桿單元除了這兩個，還有Uz。其他單元，梁單元、殼單元、體單元都包含了這三個自由度，且具有相同的物理意義，按照前面一篇文章所介紹的連接總則，桿單元與其他單元連接時只需要共用節點即可，無需建立約束方程。
下面是一個簡單的類似雨棚案例，注意本案例各構件尺寸僅為演示操作需要所擬，未經仔細推敲，各工程大佬可忽略。
某屋外雨棚平面簡化模型如上，長度為4m，折算荷載為10 KN/m，雨棚梁采用工字型鋼I40，系桿截面面積為238.64mm^2，材料均為Q235，采用ANSYS模擬該結構。
下面為建模過程
!
展開 abaqus實體-梁單元，實體-實體單元，梁-梁單元鉸接設置
使用多點約束MPC，實現實體-梁單元，實體-實體單元，梁-梁單元鉸接如何設置，實體單元梁彎矩曲線怎么提取？可下載附件，也可觀看視頻。
https://www.yqgqt.org.cn/college/video/c15810?nagivator=course
abaqus實體-梁單元，實體-實體單元，梁-梁單元鉸接設置.rar
ansys中梁單元截面類型
ansys中梁單元截面類型總共給了12種，如下圖
最后一種“ASEC”，即其他亞類，不需要形狀，只需輸入一些截面的數據即可。
ASEC類型有如下圖幾個參數：
如圖共有11種關于截面屬性的參數：A,Iyy, Iyz, Izz, Iw, J, CGy, CGz, SHy, SHz, TKz,
TKy
各個屬性所代表的參數的意義
A = Area of section 截面面積
Iyy = Moment of inertia about the y axis 對y軸的慣性矩
Iyz = Product of inertia 慣性積
Izz = Moment of inertia about the z axis z軸的轉動慣量
Iw = Warping constant 翹曲慣性矩
J = Torsional constant 扭轉常數
CGy = y coordinate of centroid y坐標的重心
CGz = z coordinate of centroid z坐標的重心
SHy = y coordinate of shear center y坐標的剪切中心
SHz = z coordinate of shear center z坐標的剪切中心
TKz = Thickness along Z axis (maximum height)沿Z軸厚度
TKy = Thickness along Y axis (maximum width)沿Y軸厚度
展開 ANSYS各類型單元連接專題講解（三）之梁與殼體鉸接
前面一篇文章主要講解了桿單元與各類單元連接的基本情況，在很多時候，我們使用梁單元的頻率要遠遠大于桿單元，因而如何處理好梁單元與各類單元的連接是做好仿真模擬的關鍵。
梁單元與桿單元不同之處在于節點除了有平動自由度之外，還附加有轉動自由度。針對2D梁單元，節點具有Ux、Uy以及Rotz三個自由度；針對3D梁單元，節點具有Ux、Uy、Uz以及Rotx、Roty、Rotz以及WaRp(僅Beam18x系列單元)。
板殼單元實際上具有五個自由度，分別為Ux、Uy、Uz以及Rotx、Roty，但很多時候引入了第六個面內轉動Rotz，但值得注意的是該自由度的含義與梁單元的Rotz含義并不相同。
2D實體單元節點自由度僅有Ux、Uy，3D實體單元節點自由度包含Ux、Uy、Uz。
從上面可見，不同單元類型其節點自由度的數目以及含義不一樣，因而在處理單元連接時，需根據實際情況分不同種類來確定其連接方法。但就梁單元而言，與各單元類型的連接可分為如下情況：
1）梁單元與殼、實體單元鉸接；
2）2D梁單元與2D實體單元剛接；
3）3D梁單元與殼單元剛接；
4）3D梁單元與3D實體單元剛接；
本篇介紹梁單元與殼、體單元的鉸接問題。
從上面介紹的三種單元節點自由度類型可見，梁單元與體單元節點的平動自由度物理意義相同，因此如果需實現梁單元與實體單元的鉸接，兩者共用節點即可；也可兩者無共用節點，但具有重合節點時，直接耦合節點的平動自由度。
然殼單元與梁單元的節點自由度除了Rotz有所不同外，其余5個自由度皆具有相同的物理意義，因而當梁單元與殼單元具有公共節點時，可認為是除了Rotz外的一種剛性連接，例如最常見的建筑結構梁板體系的模擬。
展開 Abaqus技巧之變截面梁單元 附使用ABAQUS 生成纖維梁截面下載
變截面梁單元在工程設計中經常使用，例如建筑結構中的懸挑梁就經常采用根部截面大而端部截面小的梁，在一些高聳結構如煙囪，旗桿等，變截面梁也極為常見。
在通用有限元abaqus中，實際上是存在變截面梁單元的，只是其定義方式較為隱蔽而不易被發現，本文給出在abaqus中定義采用變截面梁單元的定義方法。
（1）分別定義變截面梁兩端的profile
（2）建立梁section，選擇截面積分為before analysi，然后選擇截面沿長度變化為Tapered，接著指定start 端和 end 端的profile，并輸入相應的材料屬性。（如果是B31和B32單元需要定義橫向剪切剛度，一般在1e10左右數量級，也可參考幫助文檔的公式進行具體計算，如果需要輸出梁截面的應力，則還需要定義output points坐標作為應力輸出的位置）
其他按照普通梁單元的方式進行定義即可，以上就是定義變截面梁單元的具體步驟，使用變截面梁單元需要注意以下幾點：
（a）即使是變截面梁單元首端和末端截面不能相差太大，如果兩端面積或者慣性矩之比大于10.0，則軟件會報錯表明截面相差太大。
（b）變截面梁單元截面剛度積分只能基于變形前積分。
（c）對于一個幾何梁被劃分為多個梁單元的情況下，需要對每個梁單元分別指定不同的section，如果只定義整個幾何梁的首端和末端，可能會使得實際的梁截面是“鋸齒形”，如下圖所示：
下載地址：使用ABAQUS 生成纖維梁截面
展開 abaqus單元類型
1.單元表征 單元族：單元名字里開始的字母標志著這種單元屬于哪一個單元族。 殼單元（S）； 實體單元（C）； 梁單元（B）； 桁架單元（T）； 剛體單元（R）。 2.殼單元(S) 殼單元：可以模擬有一維尺寸（厚度）遠小于另外兩維尺寸，且垂直于厚度方向的應力可以忽略的結構。 一般殼單元：S4R,S3R,SAX1,SAX2,SAX2T。對于薄殼和厚殼問題的應用均有效，且考慮了有限薄膜應變； 薄殼單元：STRI3,STRI35,STRI65,S4R5,S8R5,S9R5,SAXA。強化了基爾霍夫條件，即：垂直于殼中截面的平面保持垂直于中截面； 厚殼單元：S8R,S8RT。二階四邊形單元，在小應變和載荷使計算結果沿殼的跨度方向上平緩變化的情況下，比普通單元產生的結果更精確； 對于給定的應用，判斷是屬于薄殼還是厚殼問題，一般：如果單一材料制造的各向同性殼體的厚度和跨度之比在1/20－1/10之間，認為是厚殼問題；如果比值小于1/30，則認為是薄殼問題；若介于1/30－1/20之間，則不能明確劃分。由于橫向剪切柔度在復合材料層合殼結構中作用顯著，故比值（厚跨比）將遠小于“薄”殼理論中采用的比值。具有高柔韌中間層的復合材料（“三明治”復合材料）有很低的橫向剪切剛度并且幾乎總是被用來模擬“厚”殼； 橫向剪切力和剪切應變存在于普通殼單元和厚殼單元中。對于三維單元，提供了可估計的橫向剪切應力。計算這些應力時忽略了彎曲和扭轉變形的耦合作用，并假定材料性質和彎曲力矩的空間梯度很小； 殼單元可以使用每個單元的局部材料方向，各項異型材料的數據，如纖維增強復合材料，以及單元輸出變量，如應力和應變，都按局部材料方向而定義。
展開 ABAQUS 單元類型和選型規則
網格單元的出現源于離散化求解，離散化把連續求解域離散為若干有限的子區域，分別求解各個子區域的物理變量，各個子區域相鄰連續與協調，從而達到整個變量場的協調與連續。每個子域內通過數學物理公式描述，單個這樣的子域就稱為單元。
有限元方法不僅應用于力場分析，還可以應用到溫度場、磁場、滲流場等分析領域，對于不同類型場的基本物理定律也是不一樣的，因此就需要用到不同的單元類型，如果需要考慮多場分析，就需要單元同時考慮所需場、以及場之間的耦合關系，也就出現了耦合單元；耦合分析雖然能獲得更準確的計算結果，但計算成本會增大，因此我們需要根據分析的問題來簡化問題，選擇合適的單元類型。
在選單元之前，我們先了解一下ABAQUS 單元編號法則，而了解單元編號法則就不得不提ABAQUS中單元具備的五個基本要素，分別是：
1）單元族群，如下圖所示為力學分析中常用的單元族群，這些族群的主要區別在于幾何特征的差異，適合于研究不同的結構類型，選擇合適的族群可以在不降低計算精度條件下，減少計算量，比如：一座高樓大廈如果全用實體單元建模，可能需要千萬甚至上億個實體單元，但如果將大廈的梁柱簡化為梁單元，墻和樓板簡化為殼單元模擬，單元數量將急劇減少。
單元編號法則1：它們的首字母或前幾位字符通常會作為單元編號的起始字符。比如：‘C3D8’中首
字母‘C’為Continuum elements 的首字母。
2）自由度，是分析過程中計算的基本變量，比如力學分析中的自由度是節點的平移和旋轉自由度；傳熱分析中需要考慮的自由度是節點溫度；滲流分析則是孔隙壓力自由度……
單元編號法則2：單元自由度通常由單元族群和尾部字符確定，比如尾部字符包含T，則表示包含溫度自由度，包含P，則表示包含孔壓自由度。
展開 Abaqus非協調模式單元類型簡介
非協調模式單元（Incompatible modes），以字符I結尾，僅適用于線性四邊形和六面體單元。它把增強單元位移梯度的附加自由度引入線性單元，能克服線性完全積分中的剪切自鎖問題，具有較高的計算精度。
Abaqus中的非協調模式單元和MSC.NASTRAN中的4節點四面體和8節點六面體單元很相似，所以計算結果頁很一致。
非協調模式單元具有如下優點：
（1）克服了剪切自鎖問題，在單元扭曲比較小的情況下，得到的位移和應力結果很精確。
（2）在彎曲問題中，在厚度方向上只需很少的單元，就可以得到與二次單元相當的結果，而計算成本明顯降低。
（3）單元交界不會重疊或開洞，因此很容易擴展到非線性、有限應變的位移。
但是使用這種單元的時候需要注意，如果所關心的部位單元扭曲比較大，尤其出現交錯扭曲時分析精度會降低。
請注意非協調模式和減縮積分單元，兩個只能選擇其一，不能同時選擇。但是同時選擇雜交單元（hybrid）。
轉自：http://blog.sina.com.cn/s/blog_b377d7f70102vew6.html
展開 BCC點陣結構梁單元Abaqus壓縮仿真模擬-顯示動力學質量縮放 ￥19.89
本文通過abaqus顯示動力學的方法對BCC結構進行壓縮仿真模擬，同時為減小計算量，采用梁單元模擬點陣結構，壓頭設置為剛性面，添加質量縮放，加快運算速度，為點陣結構壓縮模擬提供一種便捷方法。
1. 建立BCC點陣模型,以單胞尺寸5X5X5為例。
a.首先建立立方體實體，然后對實體進行處理，得到點陣單胞點陣結構。
b.建立單胞BCC梁單元點陣模型,然后進行刪除面的操作，得到單胞BCC點陣結構，接下來進行陣列操作，得到最大外形尺寸為25x25x25的點陣壓縮模擬試件。
C.建立剛性壓板，設置參考點，模擬萬能試驗機壓頭，剛性單元不參與計算，不影響計算結果，加快運算速度。
2. 裝配，按壓縮試驗進行裝配，從上到下依次為壓板-點陣-壓板。
3.設置材料屬性，本文為鈦合金TC4，密度4.43e-9Tone/mm3,彈性模量為118000MPa，泊松比0.3，應力應變值見下表所示。
設置截面屬性Beam，定義截面半徑0.5mm
指派截面，定義方向[1,2,3]，完成材料屬性設置。
4.設置分析步Dynamic,Explicit，時間設置為5s，以每秒1mm的速度進行壓縮模擬，開啟質量縮放為1e-5，歷程輸出勾選位移和力，以便輸出力-位移曲線，然后計算相應的應力-應變曲線。
5.設置相互作用-切向行為和法向行為，摩擦系數為0.3，設置通用接觸。
以下部分為付費部分
展開 
Abaqus中接觸問題中單元類型的選擇
1.關于單元階次
在接觸分析模擬中一般最好在那些將會構成從面的模型部分使用一階單元，使用二階單元可能會出現問題，這是由接觸算法決定的。
2.單元選擇
較簡單接觸問題：線性減縮積分單元（C3D8R）和非協調單元（C3D8I）。
較復雜接觸問題：修正的二階四面體單元（C3D10M ）是為了應用于復雜的接觸模擬問題而設計的，在模型復雜的接觸分析中推薦使用，但是計算時間也大大增加。
備注：具體內容請參閱莊茁的《基于ABAQUS的有限元分析和應用》，第12章--接觸
展開 abaqus顯示分析中怎么選擇橡膠單元類型
橡膠單元類型只能用雜交嗎，那顯示分析中沒有雜交單元的話用C3D8r可以嘛
ABAQUS喵星人教你看懂不同類型單元的應力方向
<p><span style="color: rgba(0, 0, 0, 0.9);">應力為典型的張量，具有明顯的坐標相關性，大家常用查看單元應力方向的方法為直接通過整體坐標系判斷XYZ方向，但這種方法僅適用于實體單元，對于其他類型單元（例如殼單元、Beam單元、Truss單元、Cohesive單元等）或特殊坐標系下的實體單元則不再適用，若仍然采用整體坐標系判定方向則會限制對后處理結果的解讀。今天喵星人就通過一個教程帶大家學習不同類型單元的應力方向應該如何看。</span></p><p class="ql-align-center"><strong>1.實體單元</strong></p><p><br></p><p>默認的實體單元應力方向服從整體坐標系，若想查看其他坐標系下的應力情況則需定義其他坐標系，建立的方式既可在前處理內定義，也可在后處理內完成，前后處理中坐標系的定義位置如下圖所示。
展開 有限元理論基礎及Abaqus內部實現方式研究系列38： 梁單元差異（2）-梁截面方向
如果是線性問題，那么Nastran和Abaqus的精度誤差主要體現在單元算法、邊界處理、MPC約束關系等，在2017年第二篇：S4殼單元質量矩陣研究文章中我們就曾經分析過Abaqus的S4殼單元和Nastran的Quad4殼單元質量矩陣的內部實現方式和差異，在這里主要研究Abaqus、iSolver與Nastran梁單元差異，由于這三款軟件的梁單元的差異較多，我們分幾篇文章來說明，本篇是Abaqus、iSolver和Nastran梁差異（2）-梁截面方向。
2.1 梁截面方向
有限元是求受力情況下的位移等變形情況，也就是位移等未知量和外力存在一定關系。對于一根三維實體梁，梁實際受到的外力是三維全局空間的，如果直接用全局坐標系下三維的力來求梁的受力分析，那么就需要對梁劃分為三維的體單元求解，網格數目和計算效率比較差，一種簡單方法是對那些細長的梁（Abaqus認為是細長比>8），此時可以用簡單的等效為線單元的形式來表達位移和外力的關系，這樣只要用一個線單元就可以表示這個三維實體梁了，大大簡化了求解矩陣。
實際的加載是多個力的組合，譬如下方采用手輪加載的力、彎矩和扭矩外載荷
但梁的有限元中可以把這個線單元受力關系分為：
（1） 軸向拉伸力
（2） 軸向扭轉
（3） 橫向彎曲力，可以加力載荷或者彎矩
三部分，此時每部分都有簡單的位移和外力的公式，也就是存在一個局部坐標系，簡化梁理論總是先求出梁單元局部坐標系的剛度和質量陣，然后再用三維變換直接轉到全局坐標系下。
對（1）（2）軸向的受力，沿梁的軸向方向，而對（3）彎曲力，沿截面方向。
展開 