線性單元
關注創建者：cadenas 創建時間：2019-10-28
線性單元的視頻教程
ABAQUS鄰近點匹配算法批量建立連接器單元/非線性彈簧模擬鋼筋混凝土粘結滑移
ABAQUS快速建立鋼筋與混凝土粘結滑移（非線性彈簧單元/連接器單元） 在ABAQUS中，通常采用非線性彈簧單元或連接器單元考慮鋼筋（線單元）與混凝土（實體單元）的粘結滑移，但已有的方法存在以下問題： 1.手動點選結點建模復雜； 2.鋼筋與混凝土結點需嚴格對應，網格劃分困難，對于曲線鋼筋或斜鋼筋幾乎無法完成； 3.插件跨平臺安裝復雜等問題。
￥400 57分鐘 2289播放
查看
線性單元的實例教程
共旋非線性單元的好處就是線性單元通過扣除剛體轉動的映射可以直接變成非線性單元，特別適用于處理大轉動小變形的幾何非線性問題。
參考文獻：
1、Felippa C A, Haugen B. A unified formulation of small-strain corotational finite elements: I. Theory[J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194(21-24): 2285-2335.
2、Sze K Y, Liu X H, Lo S H. Popular benchmark problems for geometric nonlinear analysis of shells[J]. Finite elements in analysis and design, 2004, 40(11): 1551-1569.
展開 為什么完全積分線性單元在彎曲載荷下會剪切自鎖？
以一個平面應力問題的四節點矩形單元為例。
單元的坐標系建立在中心。對于這樣一種線性單元，在構造剛度矩陣的時候，需要進行下式所示的積分。
（四節點矩形單元應該是8×8)
其中B矩陣是單元形函數對空間坐標的相關偏導，D矩陣是本構矩陣。該積分中的被積矩陣（8×8）的每一個元素都是一個三元函數，其針對單元域的積分值成為一個剛度系數。如上單元在高斯積分方案下的減縮積分就是取被積函數在積分域中心點的函數值乘以2（曾攀04P178），實際上就是梯形積分公式。
在純彎曲變形加載模式下，該剛度矩陣得出的節點位移向量解具有一定的特征，莊茁P65的圖示（本文圖1）也表示了這種特征：四個節點在2方向的位移相等，1、3節點在1方向上的位移相等，2、4節點在1方向上的位移相等，且它們互為相反數，也即我們可以得到如下形式的一個節點位移向量：
但是需注意，只有在純彎曲加載模式下，才會得到這樣形式的位移向量。
針對上面的線性矩形單元，其應變矩陣如下圖所示：
在減縮積分模式下，例如積分點(0,0)，并將得到的節點位移代入，可以得到該積分點下的應變值為：
可以看出，在該積分點處，應變的三個分量都為0。在非線性分析中，當前增量步得到積分點上的應力應變值需要代入本構曲線中，更新本構數據，進而構造下一個增量步迭代所需要的初始切線剛度矩陣。如果使用了減縮積分的線性單元，即使不是在純彎曲加載模式下，其得到的應力應變值相比理論預示值應該要小（我推測的^_^，沒空詳細證實），所以用這樣的數據構造的切線剛度矩陣相比其他單元構造的切線剛度矩陣要小，這也許就是通常所說的出現沙漏問題的單元“太軟”的緣故。
結語：本文算不得什么，只是從公式上加深了商業軟件使用者對沙漏這一現象的了解，稍微知其所以然罷了。
展開 針對上面的線性矩形單元，其應變矩陣如下圖所示：
在完全積分模式下，例如針對第四個積分點(a/√3,b/√3)，并將得到的節點位移代入，可以得到該積分點下的應變值為：
如圖中所見，該點的剪切應變不為0，這顯然不是純彎曲加載模式所要求的結果。然而需要注意，該現象是在純彎曲加載得到的節點位移和完全積分所對應的B矩陣的共同作用下得到的，如果不是純彎曲加載，那么節點位移不會有相關特征，完全積分線性單元得到的結果和相關加載模式也是符合的（莊茁P64倒數第二段）；如果純彎曲加載下的線性單元實行減縮積分，也不會出現剪切自鎖問題，但是會帶來沙漏現象，我們將在下一篇筆記中對該現象一探究竟。
結語：本文算不得什么，只是從公式上加深了商業軟件使用者對剪切自鎖這一現象的了解，稍微知其所以然罷了。如果要進一步探究如何防止剪切自鎖，要構造怎樣的位移模式，需要更多功夫，可見如下博文：
易木木響叮當，公眾號：易木木響叮當
有限元編程中如何避免剪切自鎖？（非協調單元詳解）
參考資料：
《有限元分析基礎教程》曾攀，清華大學出版社，2008.
《有限元分析及應用》曾攀，清華大學出版社，2004.
《基于ABAQUS的有限元分析和應用》莊茁等，清華大學出版社2008.
《數值分析》歐陽潔等，高教社2009.
展開 為了提供更好的客戶服務，RK Rose_Krieger與 CADENAS 專家團隊一起擴展了針對線性單元的新型產品配置器。工程師現在可以以簡單的方式組裝所需的線性單元，并根據需要在后臺進行合理性測試。在訪問RK Rose+Krieger的擴展電子產品目錄的同時，還可以使用Move-Tec“ E”新型交互式產品配置器訪問帶有主軸驅動的直線單元、升降柱和電動缸等產品組件。
概述
在使用ABAQUS中的非線性彈簧單元研究鋼筋混凝土粘結滑移、土體和樁的非線性剛度等問題時，需要在樁基和土體間建立彈簧單元。手動操作不太現實，因此本文使用python開發了腳本，可用于快速生成彈簧單元。
2.效果演示
3.核心代碼
給出核心代碼如下供大家參考，如想快速獲取需完整代碼可聯系小編（扣q1871858827）。
4.非線性彈簧單元
ABAQUS/CAE中暫時僅支持定剛度彈簧單元，如需創建非線性彈簧單元，需要在inp文件中修改關鍵字。
inp文件中修改示意（僅供參考）。
5.參考文獻
Abaqus Example Problems Guide (6.14)
Abaqus Analysis User's Guide (6.14)
展開 
線性單元的相關專題、標簽、搜索
線性單元的最新內容
采用 5mm 全局網格尺寸及線性單元完成模型網格劃分。
圖 2 模型所定義旋轉關節示意圖
5、定義分析設置并施加邊界條件。相機實際工作載荷的頻率大概率處于低頻區間，因此將分析頻率范圍設定為 0~30Hz。設置 30 個求解間隔，采用完全求解法，并設定恒定結構阻尼系數為 0.02。
綁定、無摩擦與摩擦接觸的對比分析1個月前
采用線性單元，使總節點數低于學術版軟件許可的限制。設置全局網格尺寸為 25 mm，對螺栓和節點區域采用局部網格尺寸 10 mm，對孔洞采用5 mm 的網格尺寸。網格劃分后的模型示意圖如圖 2 所示。
圖 2 網格模型的示意圖
3、定義各部件之間的接觸關系。軟件會自動在相互鄰近的部件之間設置綁定接觸。
瞬態熱傳導有限元求解器開發3個月前
單元方程
使用三角形線性單元對應的插值函數：
有些教材中，會把面積項提取出來，寫成以下這種形式，所以有的教材上剛度矩陣結果用a、b、c表達的時候，會存在差異，但是本質都是一樣的。
最終單元方程如下，其中M是熱容矩陣，K是傳導矩陣，F是熱載荷。
熱容矩陣乘的是溫度的導數。
使用高階單元
改用三角形二次單元作為速度單元，三角形線性單元作為壓力的單元，并且不使用罰函數，速度和壓力耦合求解，結果就對了。
改進后速度、壓力結果
文獻結果
對于實際應用中承受非線性彈簧單元Combin39的實際應用。
在ANSYS Workbench里提供了兩種方法，一種是WB的雙向彈簧，輸入數據表格，其本質上采用是LINK8單元進行模擬，而不是非線性彈簧combin39。
因科研需要，一直在研究一些單元算法，看著網上相關資料很多，但是和商軟對標的非線性單元技術相對較少。非線性這方面ABAQUS比較受人認可，所以打算用空余時間研究一下ABAQUS的單元技術，推導編寫一下相關程序供大家討論。本人水平十分有限，主要是學習ABAQUS的文檔，力學理論和代碼方面的問題請大家不吝賜教。
接下來，課程進入SAP2000的實際操作建模，你將學習如何正確定義非線性連接單元、分配阻尼器屬性，并將粘滯流體阻尼器集成到結構系統中。你將執行線性和非線性時程分析，解釋關鍵結果，并通過位移、加速度和力響應來評估阻尼器的有效性。
在高級部分，課程涵蓋Perform3D建模，重點關注基于性能的抗震評估、非線性組件、阻尼器校準和驗收標準。
這里不再贅述，感興趣同學可以在我的B站主頁查詢教程“ABAQUS快速建立鋼筋與混凝土粘結滑移（非線性彈簧單元/連接器單元）”。</p><p>對于輸出，用戶手冊嚴謹的給出了力與位移的輸出模式：</p><p>S11：Force in the spring.
本案例來源于喵星人課程“ABAQUS快速建立鋼筋與混凝土粘結滑移（非線性彈簧單元/連接器單元）”，模型如下圖。
幾何非線性分析：殼單元（SC8R、S4R）表現優異，收斂性和精度均優于實體類單元；CSS8 表現良好，適合復合材料非線性分析；C3D8I 表現中等，對網格質量敏感。
屈曲分析：殼單元（特別是 S4R）在屈曲載荷預測方面表現最佳，誤差 < 5%；CSS8 次之，誤差 < 7%；C3D8I 誤差較高（約 10%）。
