abaqus導出矩陣的案例
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提取單元的剛度矩陣和質量矩陣：
*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一鍵輸出剛度矩陣和質量矩陣的插件EMM ￥20
=======概述=======
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插件
展開 【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 
ABAQUS 剛度矩陣 ￥1000
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致。
Abaqus 輸出矩陣的方法
在比較深入的CAE分析過程中，我們可能會關注模型中某些部分的剛度分布情況，因此需要提取剛度矩陣參數，因為使用相對較少，目前的Abaqus前處理模塊CAE中還無法直接實現該設置，因此需要對其inp文件進行局部修改來定義。
Abaqusn輸出剛度矩陣主要分為兩大類，一類是輸出某一部分的單元的剛度矩陣，一類是輸出整體模型的剛度矩陣（可含單元、彈簧、質量、阻尼、力等）。
下面我們分別對這兩中情況分別討論：
1、輸出單元矩陣
輸出單元剛度矩陣的方式相對容易理解，當然你也可以從自帶的英文幫助中需找答案，你只需要在inp文件中添加下述關鍵字即可：
*File Format,ASCII
*Element Matrix Output,Elset=xxx,
File Name=xxx,Frequency=1,Output File=User Defined,Stiffness=Yes
其中第一行的目的是將軟件默認生成的二進制文件轉化為ASCII格式文件，以方便直接查看和讀取避免亂碼。Elset后面的xxx對應的是之前已經定義好的單元集的名稱，也就是你想要獲取的那部分單元。File Name后面的xxx是abaqus計算后所提取出的單元矩陣所保存的文件名稱，后綴為.mtx ，也就是說最終在工作目錄中會生成一個xxx.mtx的文件，可以通過記事打開的剛度矩陣文件。
展開 ABAQUS導出高清圖片的方法
在撰寫論文的時候，我們一般要用到后處理結果的高清圖片，但是ABAQUS默認的輸出是清晰度較低的。導出高清圖片需要手動設置一下即可。
在CAE界面選擇File>print，進入如下界面，選擇好保存路徑以及輸入文件名，當我們選擇了format之后（比如選擇了png格式），右側會默認勾選256色。
需要高清的圖片只需要，點掉256色，并進一步選擇，PNG Format Options即可進入分辨率設定界面。
不難看出，這些操作還是比較繁瑣的，進行這些操作之后，會產生如下的代碼：
#設定分辨率
session.pngOptions.setValues(imageSize=(4096,1792))
#設定圖片格式、名稱和保存路徑
session.printOptions.setValues(reduceColors=False)
session.printToFile(fileName='F:/PublicDemo/damage/damge3',format=PNG,
canvasObjects=(session.viewports['Viewport:1'], ))
如果需要大批量保存的時候，可以直接修改文件名，如上面個的damage3，并將上述代碼復制在界面下方的命令行，這樣新的圖片就可以自動保存在需要的目錄了，可大大提高速度。
更多精彩技術文章請掃碼關注：ABAQUS二次開發
展開 abaqus粒子生成器顆粒導出
導出部件為孤立網格，幾何編輯不能修改顆粒的孤立網格。如何才能把顆粒變成幾何體，用于下一步模擬計算，有償（能解決為前提）
用Python從Abaqus導出txt、dat數據
__biz=MzIzNzcyMDE3OA==&mid=2247484080&idx=1&sn=1b34a15a8bf79e8f9cdf901f43aac9b1&chksm=e8c505addfb28cbb54225ef35bf7c6750ee21673d696b63fe13df9067320f1215cae6861c6a5&token=667787181&lang=zh_CN#rd</strong></a></p><p><strong>以下內容均來自于本人微信公眾號（圖片版可能缺失內容，完整版請關注查看公眾號內容），有興趣請關注一波，共同學習ANSYS和ABAQUS有限元技術，我理解為兩把神器，屠龍刀和倚天劍，用好可行走有限元江湖，做個大俠。</strong></p><p><strong>同濟大學優秀畢業生創立，研究Abaqus、ANSYS在土木領域的科研編程和數值模擬技術，涉及到多尺度、多物理場耦合等計算力學問題，以及基于Python、子程序、INP的Abaqus建模，和ANSYS APDL、Workbench二次開發。
展開 HyperMesh_To_Abaqus接口——模型導入導出問題
方法三：inp導入的都是model，不過可以將不同的part export 不同的inp，導入到abaqus后，在abaqus/CAE中用copy objects將part集合成一個model。
總結：以上就是在Hypermesh和Abaqus之間導入導出可能碰到的一些問題和疑惑。解決了這些問題，至少你可以順利的實現兩個軟件之間的導入導出，如果你對一個軟件比較熟悉，那上手在另一個軟件中進行相應的設置就會更加快一些。
當然，本文只是講解了一些導入導出的問題。要想掌握 hypermesh+abaqus的方法，還是需要更深入的學習?？梢酝ㄟ^看一些高手的做的模型例子或視頻來不斷提高。
本文的例子對應的視頻連接如下：視頻連接
當然，如何在Hypermesh中設置可以直接求解的inp的方法可以學習后續的視頻，目錄文檔：
HypeMesh_For_abaqus.pdf
展開 有限元理論基礎及Abaqus內部實現方式研究系列1：S4殼單元剛度矩陣研究
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
===第一篇：S4殼單元剛度矩陣研究。
Abaqus的殼單元剛度矩陣的理論基礎都是Kirchihoff（薄殼）和Mindlin（厚殼）理論，本章重點研究S4殼單元，該單元基于Mindlin理論，在自編程序中根據Mindlin理論編寫后和Abaqus結果對比，可以發現Mindlin和Abaqus差異很大，然后結合幫助文檔猜測Abaqus的S4單元的內部修正方法。
===S4單元修正方法總結
Abaqus的S4單元的薄膜效應剛度和面外彎曲剛度矩陣是完全積分，面外橫向剪切剛度是減縮積分。
展開 
有限元理論基礎及Abaqus內部實現方式研究系列6：General梁單元剛度矩陣 ￥1
7
其它元素
√
都為0
詳細研究方法，見附件：
有限元理論基礎及Abaqus內部實現方式研究系列6：General梁單元剛度矩陣(SnowWave02 20180830).pdf
==以往的系列文章==
第一篇：S4殼單元剛度矩陣研究。介紹Abaqus的S4剛度矩陣在普通厚殼理論上的修正。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質量矩陣研究。介紹Abaqus的S4和Nastran的Quad4單元的質量矩陣。
http://www.yqgqt.org.cn/content/post/343905
第三篇：S4殼單元的剪切自鎖和沙漏控制。介紹Abaqus的S4單元如何來消除剪切自鎖以及S4R如何來抑制沙漏的。
http://www.yqgqt.org.cn/content/post/350865
第四篇：非線性問題的求解。介紹Abaqus在非線性分析中采用的數值計算的求解方法。
http://www.yqgqt.org.cn/content/post/360565
第五篇：單元正確性驗證。介紹有限元單元正確性的驗證方法，通過多個實例比較自研結構求解器程序iSolver與Abaqus的分析結果，從而說明整個正確性驗證的過程和iSolver結果的正確性。
https://www.yqgqt.org.cn/content/post/373743
展開 怎樣將Abaqus中的材料設置導出成文本形式？
怎樣將Abaqus中的材料設置導出成文本形式？
有限元理論基礎及Abaqus內部實現方式研究系列22： 幾何非線性的剛度矩陣求解
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
自主結構有限元求解器iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第22篇：幾何非線性的剛度矩陣求解==
幾何非線性在界面上是很容易設置的，但商軟內部的處理相當復雜，我們從最基本的剛度矩陣的求解出發，看看在幾何非線性設置后，剛度矩陣具體是怎么實現的。本文首先介紹幾何非線性下的剛度矩陣的理論推導和計算機求解方法，說明理想的求解方式的困難點和猜測Abaqus內部的解決方法。最后利用一個簡單的算例通過對比iSolver和Abaqus的結果，部分驗證我們對Abaqus幾何非線性的剛度矩陣的實現方式的猜測。
1.1 幾何非線性的剛度矩陣推導理論
在前面17章：幾何非線性的物理含義中，我們提到如果是非線性系統，應變能W隨t的變化就是個非線性過程。每個時刻點可以求出一個斜率，這個斜率最終會形成當前時刻點的剛度矩陣。
求導后得到的剛度K：
也就是剛度矩陣將分為兩塊：
（1） 上式的前面一部分稱為材料剛度陣，依然是以前的BDB形式，只不過B換成了當前時刻的應變位移矩陣
（2） 后面新增項一般稱為幾何剛度陣，在Abaqus中稱為初始應力矩陣（initial stress stiffness）。
展開 有限元理論基礎及Abaqus內部實現方式研究系列12： 幾何梁單元的剛度矩陣 ￥1
如果有任何其它疑問或者項目合作意向，也歡迎聯系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
詳細研究方法，見附件：
有限元理論基礎及Abaqus內部實現方式研究系列12：幾何梁單元的剛度矩陣.pdf
以往的系列文章：
第一篇：S4殼單元剛度矩陣研究。介紹Abaqus的S4剛度矩陣在普通厚殼理論上的修正。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質量矩陣研究。介紹Abaqus的S4和Nastran的Quad4單元的質量矩陣。
http://www.yqgqt.org.cn/content/post/343905
第三篇：S4殼單元的剪切自鎖和沙漏控制。介紹Abaqus的S4單元如何來消除剪切自鎖以及S4R如何來抑制沙漏的。
http://www.yqgqt.org.cn/content/post/350865
第四篇：非線性問題的求解。介紹Abaqus在非線性分析中采用的數值計算的求解方法。
http://www.yqgqt.org.cn/content/post/360565
第五篇：單元正確性驗證。介紹有限元單元正確性的驗證方法，通過多個實例比較自研結構求解器程序iSolver與Abaqus的分析結果，從而說明整個正確性驗證的過程和iSolver結果的正確性。
https://www.yqgqt.org.cn/content/post/373743
第六篇：General梁單元的剛度矩陣。介紹梁單元的基礎理論和Abaqus中General梁單元的剛度矩陣的修正方式，采用這些修正方式可以得到和Abaqus梁單元完全一致的剛度矩陣。
展開 