abaqus剛度矩陣縮減
關注創建者：王靖雯 創建時間：2023-02-27
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 剛度矩陣 ￥1000
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致。
=======概述=======
EMM(Export Matlab Matrix)是集成在ABAQUS/CAE中的一個插件，能夠一鍵輸出Abaqus模型的單元及全局剛度、質量、載荷矩陣，并自動轉換為MATLAB矩陣。
Abaqus可以輸出剛度、質量等單元或者全局矩陣。但需要手動添加關鍵詞，較為麻煩，且輸出的.mtx矩陣格式比較亂，不易閱讀。
本工具能一鍵實現以下功能：
1. 輸出Abaqus模型的剛度、質量、載荷矩陣到.mtx，包括單元和全局的；
2. 將上述矩陣轉換為MATLAB的.mat文件；
3. 用Matlab打開上述.mat文件。
如果用戶有任何問題或者需要合作，歡迎聯系我們。email: SnowWave02@qq.com
=======使用過程=======
打開Abaqus模型，切換到Job模塊，在Abaqus菜單欄的Plug-ins里看到iSolver插件的菜單。
點擊iSolver->Export Matlab Matrix…，彈出EMM界面。
勾選需要輸出的矩陣和Set，點擊EMM界面右側的Submit后將提交Abaqus計算，計算完畢得到如下類似矩陣。
=====演示視頻======
http://www.yqgqt.org.cn/college/video/c11343
=====下載======
用到的模型
Job-Cabin-T.rar
用戶手冊
Export Matlab Matrix工具說明(SnowWave02 20191217).pdf
EMM插件
展開 本期給大家推薦一款由木木自研的小工具：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 工具搜索出這兩個的安裝位置，復制進去即可。
展開 劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 
abaqus剛度矩陣縮減的相關專題、標簽、搜索
abaqus剛度矩陣縮減的最新內容
本期給大家推薦一款由木木自研的小工具：TransferMatrix，主要用于導出 Abaqus 中各種剛度矩陣，方便和自研的程序進行實時對比，基于 Pyside6 搭建的軟件界面，用戶僅需導入 inp 文件，就可以自動調用電腦內的 Abaqus 進行計算分析，不受限于 Abaqus 版本，可以導出：
單元剛度矩陣
單元質量矩陣
單元分布節點荷載列陣
整體剛度矩陣
ABAQUS 剛度矩陣11個月前
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致
參考
連接器振動脫落_abaqus重啟動分析_顯式隱式切換_插拔力預應力模態_TeeSim天深科技
Abaqus預應力模態分析 附Abaqus 分析用戶手冊材料卷下載- 技術鄰
Natural frequency extraction - SIMULIA 用戶幫助 2020 User
（原創，轉載請注明出處）
==概述==
本系列文章研究成熟的有限元理論基礎及在商用有限元軟件的實現方式，通過
（1） 基礎理論
（2） 商軟操作
（3） 自編程序
三者結合的方式將復雜繁瑣的結構有限元理論通過簡單直觀的方式展現出來，同時深層次的學習有限元理論和商業軟件的內部實現原理。
有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論
（原創，轉載請注明出處）
==概述==
本系列文章研究成熟的有限元理論基礎及在商用有限元軟件的實現方式。有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論，只是在實際應用過程中，商用CAE軟件在傳統的理論基礎上會做相應的修正以解決工程中遇到的不同問題，且各家軟件的修正方法都不一樣，每個主流商用軟件手冊中都會注明各個單元的理論采用了哪種理論公式，但都只是提一下用什么方法修正
==概述==
本系列文章研究成熟的有限元理論基礎及在商用有限元軟件的實現方式。有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論，只是在實際應用過程中，商用CAE軟件在傳統的理論基礎上會做相應的修正以解決工程中遇到的不同問題，且各家軟件的修正方法都不一樣，每個主流商用軟件手冊中都會注明各個單元的理論采用了哪種理論公式，但都只是提一下用什么方法修正，很多沒有具體的實現公式
==概述==
在CAE領域，從學校、實驗室的自研算法到實現真正的商業化軟件是一條無比漫長的道路。我們不研究有限元的新方法、新理論，只是研究商用有限元軟件的實現方式。有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論，只是在實際應用過程中，商用軟件在這些傳統的理論基礎上會做相應的修正以解決工程中遇到的不同問題，且各家軟件的修正方法都不一樣，每個主流商用軟件手冊中都會注明各個單元的理論采用了哪種理論公式
=======概述=======
EMM(Export Matlab Matrix)是集成在ABAQUS/CAE中的一個插件，能夠一鍵輸出Abaqus模型的單元及全局剛度、質量、載荷矩陣，并自動轉換為MATLAB矩陣。
Abaqus可以輸出剛度、質量等單元或者全局矩陣。但需要手動添加關鍵詞，較為麻煩，且輸出的.mtx矩陣格式比較亂，不易閱讀。
本工具能一鍵實現以下功能：
1. 輸出
==概述==
本系列文章研究成熟的有限元理論基礎及在商用有限元軟件的實現方式。有限元的理論發展了幾十年已經相當成熟，商用有限元軟件同樣也是采用這些成熟的有限元理論，只是在實際應用過程中，商用CAE軟件在傳統的理論基礎上會做相應的修正以解決工程中遇到的不同問題，且各家軟件的修正方法都不一樣，每個主流商用軟件手冊中都會注明各個單元的理論采用了哪種理論公式，但都只是提一下用什么方法修正，很多沒有具體的實現公式
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
