abaqus實現屈曲的案例
有限元理論基礎及Abaqus內部實現方式研究系列48：屈曲分析（1）-理論
（原創，轉載請注明出處）
1 概述
本系列文章研究成熟的有限元理論基礎及在商用有限元軟件的實現方式，通過
（1） 基礎理論
（2） 商軟操作
（3） 自編程序
三者結合的方式將復雜繁瑣的結構有限元理論通過簡單直觀的方式展現出來，同時深層次的學習有限元理論和商業軟件的內部實現原理。
有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論，只是在實際應用過程中，商用CAE軟件在傳統的理論基礎上會做相應的修正以解決工程中遇到的不同問題，且各家軟件的修正方法都不一樣，每個主流商用軟件手冊中都會注明各個單元的理論采用了哪種理論公式，但都只是提一下用什么方法修正，很多沒有具體的實現公式。商用軟件對外就是一個黑盒子，除了開發人員，使用人員只能在黑盒子外猜測內部實現方式。
一方面我們查閱各個主流商用軟件的理論手冊并通過進行大量的資料查閱猜測內部修正方法，另一方面我們自己編程實現結構有限元軟件iSolver，通過自研CAE軟件和商軟的結果比較來驗證我們的猜測，如同管中窺豹一般來研究的修正方法，從而猜測商用有限元軟件的內部計算方法。我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。iSolver包括完整的前后處理和有限元求解器，功能如下，有興趣可直接在下面網址下載：
百度網盤鏈接: https://pan.baidu.com/s/10d6jHdZ01SBY2JxiS6bffw 提取碼: 6fdf
2 屈曲分析
結構失效的方式有兩種：1.
展開 ANSYS經典案例在Workbench中實現之薄壁結構的屈曲與后屈曲分析
在靜力學分析中，尤其是在結構局部屈曲分析和時間相關非線性材料分析中，非線性穩定性控制（nonlinear stabilization）是一個非常好的選擇。
此外，在進行重啟動計算時，盡量避免在計算發散的時間點重啟動計算。比如本例中，計算在第十四個子步發散，在重啟動計算時，選擇在第十個子步設置nonlinear stabilization并啟動計算。如果計算依舊不收斂，則應選擇更靠前的時間點。
在后屈曲分析中，需要設置非線性穩定性控制（nonlinear stabilization）以保證計算收斂。通常有兩種方式實現：定義阻尼系數或定義能量比。用戶需要選擇通過哪種方式進行控制，然后輸入適當的值。理論上，用戶應該根據結構失穩類型、網格尺寸與類型，以及時間步長的綜合考慮進行設置。用戶可能需要經過多次嘗試才能尋找到相對合理的輸入。總之，nonlinear stabilization的目的是以最小的附加阻尼力獲得整個計算的收斂性。
（1）阻尼方法：如果用戶使用阻尼方法，意味著阻尼將定義在所有單元上。一旦阻尼值過大，附加在結構上的控制也會很大，因此結構剛度相應偏大，即結構很難產生變形；如果阻尼值過小，那么該設置不會起到控制計算收斂的作用。綜上所述，阻尼控制方法不適合用于局部非線性屈曲失效的結構。
（2）能量方法：由于本案例中的薄壁結構的失效是由于局部屈曲失效引起的，所以能量方法更加適合本案例。能量方法的原理是在不同的單元處施加不同的阻尼，所以程序會給容易失穩的單元施加較大的阻尼，而給不易失穩的單元施加較小的阻尼，這樣既可以控制計算的穩定性，同時也可以將附加控制力降低到最小。本案例中取能量比為0.000143進行計算。
展開 SOL101線性靜力求解，SOL103模態求解，SOL105屈曲求解，sol106非線性靜力求解設置在hypermesh軟件界面中操作實現 ￥18.8
常見的求解類型包括SOL101線性靜力求解，SOL103模態求解，SOL105屈曲求解，sol106非線性靜力求解，sol145顫振分析求解，sol129非線性動力求解，sol107轉子復特征值分析（轉子臨界轉速）求解。
其中SOL145、SOL129、SOL107求解設置無法全部通過hypermesh軟件進行設置，建議在MSC PATRAN中設置后存為對應的求解bdf模板，供后續參考，其他建議通過hypermesh軟件設置后存為求解模板。
本文主要介紹SOL101線性靜力求解，SOL103模態求解，SOL105屈曲求解，sol106非線性靜力求解設置在hypermesh軟件界面中如何操作實現。
展開 ABAQUS與MidasCivil 在屈曲分析的對比
<a href="/major/ABAQUS與Midas Civil 在屈曲分析的對比
[摘耍]本文是基于Abaqus和Midas Civil采用梁單元，對實腹矩形截面構件在軸心受壓情況下發生彎曲失穩的線性屈曲分析。通過考慮材料線性得出構件發生彎曲失穩的特征值。通過保持構件的截面、長度和荷載不變，只改變邊界條件，進而實現得到不同邊界條件彎曲失穩的特征值，用兩種仿真軟件進行比較，供計算屈曲的參考。
[關鍵詞] 特征值 屈曲
1、 計算機配置情況
2、 計算時間
第一種工況
第二種工況
3、 模型參數：
構件尺寸（單位：mm）：
1500×1000×250
材料屬性：
彈性模量:1.0×10-2tonf/mm
荷載：
軸向載荷集中荷載1tonf，不考慮材料自重。
展開 
Abaqus 非線性屈曲分析方法 附ABAQUS分析手冊分析卷下載
當然，對于方筒這類實際上是通過顯示方法實現的，更準確的講是動力屈曲分析，所以我們還得判斷動能、塑形耗散等能量參數，才能使結果更加準確。
下載地址：ABAQUS分析手冊分析卷
abaqus非線性屈曲
線性屈曲分析
*buckle
用于估計最大臨界載荷和屈曲模態，無法查看屈曲后狀態。可用作引入缺陷的之前的計算分析步，需要加載荷；屈曲特征值與載荷相乘就是屈曲載荷。主要用于缺陷不敏感結構。
非線性屈曲分析
*static, riks
用于計算最大臨界載荷和屈曲以后的后屈曲響應，可以查看后屈曲狀態，用弧長量代替時間量。載荷比例因子與載荷相乘就是屈曲載荷。可以用于缺陷敏感結構，如果結構存在接觸，容易出現收斂問題。
通用靜力分析
*static
用于計算結構剛度不變或結構剛度增大的結構，如果結構出現屈曲或者垮塌，很容易出現不收斂問題，無法計算后屈曲狀態。
通用靜力分析+阻尼穩定
*static, stabilize
在靜力分析步中加阻尼，有助于收斂，計算的結束點可以比通用靜力分析要后一些，但要注意阻尼不能加得過大。
隱式動力分析
*Dynamic
將屈曲問題作為隱式動力問題來處理，適合接觸脫開的問題，但是假如結構接觸對較多，很容易出現收斂問題。這種分析類型使用的是隱式積分方法。
顯式動力分析
*dynamic, explicit
將屈曲問題作為顯式動力問題來處理，適合接觸脫開的問題，能夠適應復雜的模型，復雜的接觸對， 收斂效果較好。但是計算量較大，計算時間較長，計算完以后需要評估計算結果是否可靠。這種分析類型使用的是顯式積分方法。
展開 abaqus易拉罐受壓屈曲
abaqus易拉罐受壓屈曲，需要原文件的關注抖音abaquser，私信易拉罐即可
Abaqus復合材料非線性屈曲分析
對一個復合材料結構進行非線性屈曲仿真分析求解他的極限承載能力，是不是需要對該材料的塑性屬性參數設置？這個屬性一般用什么方法設置啊？各項同性材料就是根據應力應變曲線獲得他的應力以及對應的塑性應變，各向異性材料也是一樣嘛？
ABAQUS非線性屈曲分析
屈曲分析主要用于研究結構在特定載荷下的穩定性以及確定結構失穩的臨界載荷，屈曲分析包括： 線性屈曲和非線性屈曲分析。線彈性失穩分析又稱特征值屈曲分析；線性屈曲分析可以考慮固定的預載荷，也可使用慣性釋放；非線性屈曲分析包括幾何非線性失穩分析， 彈塑性失穩分析（材料非線性失穩分析）， 非線性后屈曲分析（包含幾何非線性和材料非線性）。
ABAQUS屈曲分析有三種方法：
1、直接施加極值載荷，拉出力-位移曲線，查看區區狀態。這種方式不適合對稱結構，如一塊板、或圓筒，軸向加載時分析不出屈曲效果；
2、特征值屈曲分析方法，可以評估結構的屈曲臨界值，但是只能是線性分析；
3、Riks法，這種方法可以計算最大臨界載荷和屈曲后的后屈曲響應，可查看后屈曲狀態，可以考慮材料非線性、幾何非線性及初始缺陷的影響，其中初始缺陷通過特征值屈曲模態、振型及一般節點位移來表述。
我們此次課程中采用屈曲分析方式，先計算屈曲模態，也就是先做特征值屈曲分析，此分析為線性屈曲分析，在小變形的情況下進行，得出臨界載荷（一般取一階模態的eigenvalue乘以加載的單位載荷1），且需要在inp文件中輸入如下圖字符，輸入次字符的目的是將初始缺陷的節點輸出為.fil文件；然后將1階屈曲模態做為初始缺陷引入極限載荷后屈曲分析，后屈曲分析可以定義非線性材料及幾何非線性，所以risk屈曲分析也成為非線性屈曲分析.
展開 abaqus屈曲模態分析教程詳解 ￥10
abaqus屈曲模態分析教程詳解
視頻下方附帶工程文件inp，大家可以自行下載學習參考
Abaqus 非線性屈曲分析方法
在有限元分析中，我們主要通過屈曲分析（Buckling Analysis）去判斷發生屈曲的臨界載荷大小。而這其中根據實際結構和要求的不同又分為線性屈曲分析（通常直接簡稱為屈曲分析）和后屈曲分析。當然，如何涉及非線性問題，后屈曲分析是必要的，不過對于后屈曲分析的實現方式也會更加麻煩一些，因為需要局部調整inp關鍵字達到目的，但只要掌握了關鍵點，依葫蘆畫瓢還是非常湊效的。
在Abaqus中對于屈曲的計算考慮則依據結構的復雜性而定，簡單的可以只考慮線性屈曲分析預估臨界載荷大小；對于較復雜的模型，則可以考慮Riks法進行后屈曲計算，從而可獲取屈曲以后的結構響應情況；但對于涉及接觸脫開等特別復雜的問題可能得借助Explicit來實現；而對于局部褶皺問題需要借助Static,Stabilize來實現。
1
線性屈曲分析
線性屈曲分析用于預估臨界失穩載荷和失穩模態；所求得的屈曲特征值與所加載的載荷大小相乘就是臨界失穩載荷；當然，對完善結構的屈曲問題，線性屈曲分析也是為后屈曲分析引入缺陷（擾動）做好準備，這是非常關鍵的。
在Abaqus中進行線性屈曲分析的方法是通過Buckle進行的。
一般線性屈曲分析只需要關注第一階屈曲模態，并根據計算所得的第一階屈曲載荷因子預估使結構發生屈曲所需要的臨界載荷是多大。但通常而言線性屈曲分析得到的臨界失穩載荷大小是保守的，偏大的。為了獲取更加準確的結果，特別是復雜模型，就需要進行非線性屈曲分析（或稱為后屈曲分析）。
因此通常會在線性屈曲分析中考慮添加關鍵字作為后屈曲分析的擾動引入參數。
展開 
求助abaqus 管道非線性屈曲分析
要求一段管道，內部有內壓，側面有側壓，分析受力變形。 有哪位大佬有教程或者模型，可有償。
abaqus屈曲分析實例操作視頻教程
分享abaqus屈曲分析實例操作視頻資料，讓高手見笑了，資料很不錯很詳細，值得一看
abaqus屈曲分析實例操作視頻教程 2.rar
abaqus屈曲分析實例操作視頻教程 1.rar
---------------------------------------------------------------------------------------------------------------------------------------------
技術鄰推薦：
Abaqus 中創建零厚度cohesive單元的幾種方法
ABAQUS焊接模擬-移動熱源（DFLUX）-平板對接不帶生死單元圖文介紹
ABAQUS規格體網格陣列使用技巧
展開 abaqus 2017屈曲分析后處理odb轉vtu python文件 ￥100
（1）abaqus 2017屈曲分析后處理odb轉vtu python文件
（2）單元介紹
##############################################后處理函數
# CAX3: 三節點三角形單元，用于二維和三維分析。
# CAX4R: 四節點四邊形單元，用于二維和三維分析。
# C3D8: 八節點六面體單元，用于三維分析。
# C3D8R: 八節點六面體單元，具有簡化的積分方案，用于三維分析。
# C3D6: 六節點楔形單元，用于三維分析。
# C3D4: 四節點四面體單元，用于三維分析。
# S4: 四節點殼單元，用于二維和三維分析。
# S3: 三節點殼單元，用于二維和三維分析。
展開 利用ABAQUS進行屈曲梁負剛度分析
定義材料參數
設置分析步
第一分析步中，限制梁兩端節點只能沿梁軸線方向滑動，同時給梁中部節點一個向上的微小位移，這里定為0.1mm；第二步里，給梁兩端的節點一個相對的位移，均為0.75mm，以模擬梁的壓縮狀態；第三步里，分析步選用非線性屈曲分析，在第二步的基礎上，給中部節點施加一個向下的力，定為1.
輸出設置
定義載荷
載荷的加載與分析步相對應，BC-1和BC-2為固定梁兩端五個自由度，僅保留水平方向自由度，BC-3為給垂向一個位移0.1mm；BC-3和BC-4為給梁兩端一個0.75mm的位移；在第二分析步中已經令梁兩端有了一個相對位移，默認到達移動位置后固定，因此BC-6和BC-7為冗余約束，此處作了無效處理。
加載
單元選擇
選擇二次單元
提交計算
計算之前修改分析步第三步的步長為固定值0.0005，共500步，如圖所示。
計算結果
兩端均壓縮0.75mm后梁的變形為：
中部加載前梁的應力分布云圖
加載后梁的應力分布云圖：
將數據導出，在MATLAB中繪圖，得到力位移曲線：
由圖中數據計算可知梁的負剛度為：
按柔性設計手冊里的公式：
可見兩者相差極小，可認為仿真結果有效。取上圖局部放大，得到壓力為0時梁的位移情況如下圖所示。
由于仿真計算采樣點的限制且誤差不可避免，通過上圖可知，位移為0時，受力也為接近0。
展開 