abaqus降低剛度的案例
ABAQUS/Explicit顯式計算成本的降低
加速模擬過程降低計算成本
為了減少所需的增量步數，n，可以加快模擬（相對于事件的實際時間），即人為地縮短事件的時間周期，但是這可能引發兩個錯誤：
（1） 如果模擬速度過快，增加的慣性力將嚴重影響之后的分析結果。避免該情況的唯一方法是適當的提速，不能過大。
（2） 除了慣性力，其他因素，例如材料特性也可能與速率有關。在這種情況下，所模擬事件的實際時間不能更改。
這種方法很常見，例如在模擬貫入問題時，現場實際的時間可能是幾個小時，貫入速度為幾毫米每秒，而模擬時通常將分析時間設置為幾秒或幾分鐘，來降低計算成本。
Step time（/s）
10
5
2.5
速度（mm/s）
50
100
200
CPU time（/s）
38.4
26.6
19.3
注：此處使用了相同的多線程和質量縮放。
加速模擬過程降低計算成本效果顯著，但是速度必須合適。
3. 質量縮放降低計算成本
人為地增加材料密度ρ，假定放大因子為f2，則增量步數由n降低到n/f，T降低到T/f，這個概念稱為“質量縮放”。它降低了事件時間與波動在單元間傳播的時間之比，同時使事件時間保持不變。這允許在分析中包含與速率相關的特性。但必須謹慎地使用，以確保慣性力不會對分析結果產生顯著影響。質量縮放對慣性力的影響與加速模擬時間的方式相同。質量縮放可以通過改變密度來完成，但是ABAQUS中提供了更多的方法來縮放整個模型或模型中的特定單元集。
展開 ABAQUS 剛度矩陣 ￥1000
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致。
【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 TransferMatrix：導出Abaqus剛度矩陣
本期給大家推薦一款由木木自研的小工具：TransferMatrix，主要用于導出 Abaqus 中各種剛度矩陣，方便和自研的程序進行實時對比，基于 Pyside6 搭建的軟件界面，用戶僅需導入 inp 文件，就可以自動調用電腦內的 Abaqus 進行計算分析，不受限于 Abaqus 版本，可以導出：
單元剛度矩陣
單元質量矩陣
單元分布節點荷載列陣
整體剛度矩陣
整體質量矩陣
可選是否導出后打開 Matlab ，并保存為 mat 文件
默認保存的格式為 csv，用戶可以直接在TransferMatrix 中查看
可選是否保存整體剛度/質量矩陣為稀疏矩陣/全矩陣形式
軟件界面見下圖：
若勾選導出后打開matlab，可直接查看：
軟件實現原理
使用注意事項
軟件啟動速度較慢，需要有點耐心，鄙人對于 PySide 的打包技術認識有限
若 inp 文件內沒有密度參數，是不能導出質量矩陣的，建議在動力類型的分析步中進行導出質量矩陣
非協調單元和雜交單元只能輸出質量矩陣
在輸出載荷列陣時，Abaqus 產生的是單元節點的分布載荷，如果沒有這中載荷類型，也是無法導出的
建議再導入 inp 文件之前，確保這個文件能在 Abaqus 上跑通，本軟件只是負責轉換剛度矩陣
不支持在TransferMatrix 中查看整體剛度/質量矩陣，因為一般情況下全局矩陣尺寸很大，可選擇導出至 Matlab 中進行查看
如果出現風險提示，請選擇信任：
軟件同級目錄放置的config.json文件，里面配置的是本地的 abaqus.bat 和 matlab.exe 目錄，建議用 everything 工具搜索出這兩個的安裝位置，復制進去即可。
展開 
基于ABAQUS的橡膠懸置膠合件剛度仿真計算
本文將以一個懸置膠合件仿真的實例講解一下如何利用ABAQUS來獲取其三個方向的靜態特性。所用膠合件的數模圖如圖1所示。其設計圖紙上標注的三向剛度如表1所示，膠料硬度是邵氏50±5度。
圖1 膠合件結構
表1設計要求
1、 網格劃分
采用HYPERMESH對圖一懸置進行網格劃分到的有限元模型如圖2所示。
2、材料設置
把劃分好的網格導入ABAQUS中，設置其材料參數，由于不同本構模型對橡膠懸置膠合件剛度計算結果有一定的影響。結合何小靜，上官文斌發表的《橡膠隔振器靜態力－ 位移關系計算方法》一文的研究結果表明，Mooney-Rivlin 模型的計算精度最高，其相對誤差均小于10%，所以本文采用M-R模型進行計算。50度膠料的M-R材料常數C10=0.2969，C01=0.0584。
3、剛度求解
3.1求解X方向剛度
按表 1要求，做如下設置：在Z方向先預載8mm，再在X向加載500N。取值0~5.6mm,對X向靜剛度進行求解。
求得的力和位移關系見表2所示，用表中數據進行畫圖差值可得到圖3所示的X向靜剛度為38N/mm,與設計值非常接近，其變形云圖見圖4所示
表2 X向力和位移關系表
圖3 X向剛度差值結果
圖4 X向云變形圖
3.2求解Y方向剛度
按表 1要求，做如下設置：在Z方向先預載8mm，再在X向加載1000N。取值2~4mm,對Y向靜剛度進行求解。
求得的力和位移關系見表3所示，用表中數據進行畫圖差值可得如圖5到Y向靜剛度為98N/mm,與設計值80N/mm有一定差異，見圖4。其變形云圖見圖6.
展開 ABAQUS疑難雜癥之沙漏剛度
源文件可在公眾號:‘易木木響叮當’ 內回復“沙漏剛度”自動獲取。
今晚直播 | 基于ABAQUS的橡膠減震件剛度分析
“
為了提高仿真工程師、高校師生實際工程能力，技術鄰特開展2021年ABAQUS系列直播課，我們甄選了四個熱門方向（巖土、二次開發、橡膠分析、混凝土）的基礎入門課，助力小伙伴們夯實有限元基礎。
第三期直播《基于ABAQUS的橡膠減震件剛度分析》將于今晚開啟，歡迎大家關注學習！
”
橡膠材料由于其獨特的物理和化學的特性（如超彈性，粘彈性且柔軟性、耐磨性、絕緣性和阻隔性等），使得其在工程上得到了非常廣泛應用，這一點在汽車行業尤為明顯。縱觀過去近200年的歷史，硫化橡膠的誕生直接推動了汽車革命。
如今在我們的汽車中，橡膠制品早已是“汽車的半條命”。就拿我們常見的桑塔納轎車來說，其就擁有270多個橡膠密封制品，而這些橡膠組件的性能直接決定了汽車的性能和安全。
展開 基于ABAQUS的鋼筋混凝土梁的剛度分析
ABAQUS 是一套功能強大的工程模擬的有限元軟件，其解決問題的范圍從相對簡單的線性分析到許多復雜的非線性問題。 ABAQUS 包括一個豐富的、可模擬任意幾何形狀的單元庫。并擁有各種類型的材料模型庫，可以模擬典型工程材料的性能，其中包括金屬、橡膠、高分子材料、復合材料、鋼筋混凝土、可壓縮超彈性泡沫材料以及土壤和巖石等地質材料，作為通用的模擬工具， ABAQUS 除了能解決大量結構（應力 / 位移）問題，還可以模擬其他工程領域的許多問題，例如熱傳導、質量擴散、熱電耦合分析、聲學分析、巖土力學分析（流體滲透 / 應力耦合分析）及壓電介質分析。
本文采用ABAQUS簡單搭建鋼筋混凝土結構有限元模型，解決工程實際問題。本文對某鋼筋混凝土模型進行模擬，采用實體單元（混凝土）與Beam單元（鋼筋）相結合的方式進行模擬。首先導入模型，劃分實體網格和一維梁單元。建立材料和截面，分為實體和一維梁單元。然后將材料賦予個模型，進行網格劃分。然后加載，簡支梁模型，約束兩端下部，在上部施加均勻壓力。最后求解計算。
由結果處理可知，鋼筋混凝土梁的最大位移為2.55cm，可以達到模擬效果。ABQUS可以有效模擬鋼筋混凝土模型，對工程實際有很大的參考指導意義。
設備基本情況：i5 4核 ；耗時：10分鐘內。
圖1 鋼筋混凝土梁CAD模型
圖2 線框模型，中間為鋼筋模型
圖3-4 網格模型
圖5 施加載荷示意圖
圖6 位移云圖
展開 利用ABAQUS進行屈曲梁負剛度分析
計算結果
兩端均壓縮0.75mm后梁的變形為：
中部加載前梁的應力分布云圖
加載后梁的應力分布云圖：
將數據導出，在MATLAB中繪圖，得到力位移曲線：
由圖中數據計算可知梁的負剛度為：
按柔性設計手冊里的公式：
可見兩者相差極小，可認為仿真結果有效。取上圖局部放大，得到壓力為0時梁的位移情況如下圖所示。
由于仿真計算采樣點的限制且誤差不可避免，通過上圖可知，位移為0時，受力也為接近0。
基于ABAQUS超彈性材料橡膠襯套的剛度計算 附基于Abaqus的橡膠和粘彈性建模下載
對于一些結構簡單的橡膠制品，我們可以基于一些理論推導或工程經驗算法在設計初期來獲取其靜剛度特性。但由于橡膠具有非線性粘彈性與超彈性，這種理論計算結果往往與試驗存在一定誤差，并且這種誤差在一般情況下是不可以忽略不計的，其具有一定的工業應用價值。
為減小誤差或實現零誤差的前期預測，我們引入了有限元仿真分析技術，其可以通過控制模型參數與網格質量實現較小誤差的預測計算。其價值也在各個行業實際的生產中得到了很好的驗證。本文基于減振襯套簡單講訴一下基于ABAQUS軟件的橡膠制品靜剛度仿真分析過程。
仿真分析過程可分為三個大過程：前處理、求解計算和后處理。本文基于ABAQUS軟件設定的分析步驟，不再重點區分分析的三個過程，將操作過程拆分為：部件、屬性、裝配、分析步與輸出設置、相互作用、網格、加載、作業提交與監管以及計算結果的可視化處理九個模塊，下面講訴橡膠襯套靜剛度仿真分析過程。
一、部件
由于本文主旨是為介紹橡膠剛度仿真的過程，所以選用了結構較為簡單的橡膠襯套為例，直接借助ABAQUS軟件的部件模塊常見如圖1所示的幾何模型。
圖1、幾何模型結構圖
二、屬性
為了使仿真結果更接近與實驗值或真實值，除了需要一個適合的仿真求解器和一個高質量的網格文件，更需要選擇一個合適的橡膠本構模型，在ABAQUS軟件中內置了許多相對成熟的橡膠本構模型（如圖2所示），我們可以通過指定相關的系數來實現本構模型的定義,當然我們還可以直接提交我們的試驗數據，交由ABAQUS軟件進行擬合，得出相對精準的參數。
展開 (一維彈簧單元的直接剛度方法)Python編程和ABAQUS結果對比
這個在我去年的時候就已經發過帖子了，但是沒有講解代碼，沒有和ABAQUS有限元軟件對比。這次發帖子就是想講解代碼以及和ABAQUS結果對比，發現我又更深層次理解了有限元方法。
這次寫得是最簡單的模型：一維線彈簧單元。采用的是直接剛度法求解。
直接剛度法的求解思路如下所示，其中粉紅色的是輸入，淡藍色是輸出。主要是要區分齊次邊界條件和非齊次邊界條件，非齊次邊界條件的話就要修改【F】。
下面將貼出我用python寫得一維彈簧單元的直接剛度法：
例子計算：
如圖是一個彈簧系統，單元節點信息如下，5節點受到一個強制位移20mm，明顯這是一個非齊次邊界條件問題。
Python編程輸入信息如下：
結果如下：
可以看到，輸出結果和書上的答案一致。
下面進行ABAQUS模擬：
添加彈簧單元
添加邊界條件
顯示單元編號、節點編號如下所示，紅色表示單元編號，黃色是節點編號
ABAQUS計算結果如下：
首先是變形圖前后對比
反力云圖如下所示，基本和直接剛度法計算的結果一致
位移云圖如下所示，基本和直接剛度法計算的結果一致
整體剛度矩陣如下所示，因為ABAQUS彈簧單元是三維的，每一個節點有3個自由度，15x15，原味的剛度矩陣如下
我們把剛度矩陣轉化為一維的，方便和編程的結果對比
從結果可以看出，ABAQUS的整體剛度矩陣和直接剛度法計算出來的整體剛度矩陣有些差異，如圖標紅的所示。
那么在整體剛度矩陣上為什么ABAQUS會和直接剛度法的整體剛度矩陣有差異呢，到底ABAQUS的整體剛度矩陣對不對呢？答案將在下一期揭曉。歡迎大家積極討論。
展開 
續集(一維彈簧單元的直接剛度方法)Python編程和ABAQUS結果對比
我們編程的結果如下所示
ABAQUS結果和編程結果對比
所以ABAQUS提取的整體剛度矩陣實際是經過補償后的剛度矩陣，嚴格來說是不正確的，但是并不影響力和位移的結果。
正確的剛度矩陣如下
如果疑問，歡迎交流和指正。
補償法的部分代碼如下
?abaqus提取單元的剛度矩陣和質量矩陣
abaqus提取單元的剛度矩陣和質量矩陣：
*Output, history, variable=PRESELECT,
*File Format,ASCII
*Element Matrix Output,Elset=Beam-1.Set-2,
File Name=shuchu,Frequency=1,Output File=User Defined,Stiffness=Yes
*End Step
*Step
*Matrix generate，stiffness
*End Step
有限元理論基礎及Abaqus內部實現方式研究系列6：General梁單元剛度矩陣 ￥1
本文首先簡單介紹梁單元的基本理論，分析了每一部分剛度的來源，并研究了Abaqus中General梁的B31單元的剛度矩陣的修正方式，采用這些修正方式可以得到和Abaqus完全一致的剛度矩陣，自研軟件的正確性證明是非常困難的，無論你測試多少標準算例或者和試驗結果對比，都很難讓用戶踏實的相信你的結果和商業軟件是一致的，但如果一個模型僅有General梁組成，那么無論這個模型多么的復雜，我們一般都可以得到和Abaqus完全一致的分析結果。最后我們使用了一個簡單的算例綜合驗證，針對該算例我們對比采用了同樣修正的自編程序iSolver和Abaqus的結果，可以發現結果如預期的一樣，沒有任何誤差。
==演示視頻==
https://www.yqgqt.org.cn/college/video/c12884?chapter=1
梁的靜力分析結果校核視頻演示任意分段梁的靜力分析，證明iSolver的分析結果和Abaqus沒有誤差，有兩個算例：Generalized梁和L型梁。
==總結==
General梁的B31單元的剛度矩陣在Timoshenko梁理論基礎上的修正如下表：
項次
剛度
修正
不修正
說明
1
軸向拉伸剛度
√
2
橫向彎曲剛度
√
采用減縮積分
3
軸向拉伸和橫向彎曲耦合剛度
√
根據形心偏置進行修正
4
軸向扭轉剛度
√
剪切中心的偏置會影響抗扭剛度系數J
5
橫向剪切剛度
√
增加了一個幾何因子，使得細長梁的時候該項趨于0。
6
軸向扭轉和橫向剪切耦合剛度
√
根據剪切中心的偏置進行修正。
展開 基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 