ansys 剛度矩陣輸出的案例
nastran 中如何輸出質量矩陣和剛度矩陣
我記得在patran 中輸出的bdf文件中輸入一條命令，運行一下就可以了。但是記不得了。請教高手，先謝謝了。
MSC Nastran模型剛度矩陣和質量矩陣的輸出方法
1 概述
MSC Nastran模型的剛度矩陣和質量矩陣，可以輸出為文本文件。工程實際中，工程師可以校核、集成矩陣，進行第二次開發，完成商用軟件和自研程序的完美集成。例如：工程師有一個計算線性動力學方程組的瞬態python程序，可以集成MSC Nastran的剛度矩陣和質量矩陣。
2 剛度矩陣和質量矩陣的輸出方法
1) 剛度矩陣和質量矩陣輸出至punch（.pch）文件
如果需要在其他MSC Nastran計算中，重用MSC Nastran模型的矩陣，可以將MSC Nastran矩陣輸出至Punch文件，方法為：
l 在MSC Nastran卡片中，添加參數：PARAM,EXTOUT,DMIGPCH
注：Punch文件中的矩陣，Patran不支持
2) 剛度矩陣和質量矩陣輸出至f06（.f06）文件
如果想直接在f06中查看輸出的矩陣，可以使用如下方法：
l 在執行控制部分（CEND前），添加如下卡片：
COMPILE EXTOUT $
ALTER 'RETURN'(,-1) $
MATPRN KAA,,,,// $
MATPRN MAA,,,,// $
l 添加如下參數（BEGIN BULK），例如：PARAM,EXTOUT,DMIGPCH
實例：
輸入文件：
剛度矩陣
質量矩陣
3 參考信息
適用版本：MSC Nastran 2005及以后版本。
展開 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插件
展開 x.y MSC Nastran模型剛度矩陣和質量矩陣的輸出方法
1 概述
MSC Nastran模型的剛度矩陣和質量矩陣，可以輸出為文本文件。工程實際中，工程師可以校核、集成矩陣，進行第二次開發，完成商用軟件和自研程序的完美集成。例如：工程師有一個計算線性動力學方程組的瞬態python程序，可以集成MSC Nastran的剛度矩陣和質量矩陣。
2 剛度矩陣和質量矩陣的輸出方法
1) 剛度矩陣和質量矩陣輸出至punch（.pch）文件
如果需要在其他MSC Nastran計算中，重用MSC Nastran模型的矩陣，可以將MSC Nastran矩陣輸出至Punch文件，方法為：
l 在MSC Nastran卡片中，添加參數：PARAM,EXTOUT,DMIGPCH
注：Punch文件中的矩陣，Patran不支持
2) 剛度矩陣和質量矩陣輸出至f06（.f06）文件
如果想直接在f06中查看輸出的矩陣，可以使用如下方法：
l 在執行控制部分（CEND前），添加如下卡片：
COMPILE EXTOUT $
ALTER 'RETURN'(,-1) $
MATPRN KAA,,,,// $
MATPRN MAA,,,,// $
l 添加如下參數（BEGIN BULK），例如：PARAM,EXTOUT,DMIGPCH
實例：
輸入文件：
剛度矩陣
質量矩陣
3 參考信息
適用版本：MSC Nastran 2005及以后版本。
展開 
NASTRAN輸出剛度矩陣的方法
在特殊情況下有可能需要結構的剛度矩陣，簡單的可以通過編程或者手算得到，但對于較為復雜的結構，這樣做顯然不劃算。NASTRAN提供了輸出單元剛度矩陣的方法，下面就輸出方法做簡單介紹。
1.輸出剛度矩陣的設置
首先你需要建立有限元模型，并設置自由-自由模態分析工況。你可以通過卡片PARAM\UNSUPPORTED_PARAMS設置輸出剛度矩陣，也可以通過在輸出的bdf文件中添加指定代碼PARAM,EXTOUT,DMIGPCH,將剛度矩陣輸出到pch文件中，pch文件中的輸出效果見圖1.
圖1 輸出的剛度矩陣
2.剛度矩陣解釋
圖1出自某個不具名字的大神，解釋的很清楚。小翼在這里簡單描述下我的理解。首先第一行是剛度矩陣的介紹，詳細的含義可以參見nastran的幫助文檔。對小翼有用的僅僅是矩陣總列數，132這個數字說明該結構共有132/6=22個節點，剛度矩陣是132行X132列。每個節點6個自由度，把總剛度矩陣想象成22*22個6X6的子矩陣，一個子矩陣代表一個節點，這樣容易理解數值怎么往總剛度矩陣里面填寫。
從第二行開始是剛度矩陣的具體數值（黑色下劃線數值-剛度系數）和位置信息，紅色下劃線數字代表列節點編號Ci，綠色代表行節點編號Ri(i=1,2,...,22)，由Ci和Ri即可確定剛度系數屬于哪一個子矩陣。藍色數字代表剛度系數在子矩陣中的列數cm，粉紅色下劃線代表黑色下劃線數值在子矩陣中的行數rm，（m=1,2,...,6）行數列數配合可以確定剛度系數在子矩陣中的位置。例如圖一中帶下劃線的剛度系數在22X22矩陣中的K11子矩陣中的k11位置。
展開 如何從Ansys APDL中提取剛度矩陣與質量矩陣？ ￥69
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys APDL運行中的產生的各種數據（例如：剛度矩陣、質量矩陣）導出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實際工作或學習中需要用到此類技能的同學、同事們，讓大家更了解Ansys APDL背后的工作原理與數據導出方式。
當然，在社區中早就有大佬回答過了這個問題，并給大家制作了相應的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進一步學習了解軟件背后機理的群體，并在此基礎上保留教學的簡潔性，提供導出矩陣與轉換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因為優化不完全導致的運行bug。
2.有限元軟件導出剛度矩陣與質量矩陣的方法
在使用APDL進行求解時，每次在求解完成后都會在工作路徑下生成一個.full文件，而這個文件十分關鍵，其正是剛度矩陣與質量矩陣的所在之處。
展開 基于 MATLAB 的 ANSYS Harwell-Boeing 格式稀疏矩陣提取工具 —— 剛度矩陣與質量矩陣 ￥30
在有限元分析中，ANSYS 可以導出大規模稀疏矩陣（如剛度矩陣、質量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續二次開發、動力學分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數，可直接從 ANSYS 導出的 HB 矩陣文件中讀取并重構成 MATLAB 稀疏矩陣：
1.剛度矩陣提取函數
輸入：ANSYS 導出的剛度矩陣 HB 文件（stiff.txt）
輸出：MATLAB 稀疏矩陣 K，可直接用于動力學計算或驗證
支持自動對稱化，保證數值正確
2.質量矩陣提取函數
輸入：ANSYS 導出的質量矩陣 HB 文件（mass.txt）
輸出：MATLAB 稀疏矩陣 M
使用與剛度矩陣同樣的解析邏輯，無需額外修改
案例說明：
本文以高速鐵路接觸網結構為例，展示了如何將 ANSYS 中導出的稀疏剛度矩陣和質量矩陣，在 MATLAB 中完整展開，并進行后續動力學分析準備。
通過該方法，可以將大規模有限元矩陣快速轉化為 MATLAB 可操作形式，為自定義振動分析、模態分析及其他科研或工程應用提供基礎。
優勢與應用：
支持大規模稀疏矩陣解析
自動對稱化，保證數值精度
適用于剛度矩陣、質量矩陣、其他 HB 格式矩陣
可作為動力學求解器或后處理工具的基礎模塊
使用方法：
1.使用以下代碼對ansys中生成的質量及剛度矩陣進行提取，file,5,full（5為工作目錄下full文件的文件名，例如：filename.full）。
展開 【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 提取ANSYS中的剛度矩陣和質量矩陣到MATLAB中進行二次開發 ￥88
最近在考慮自己編寫的程序和商用軟件的驗證問題，有限元結構分析中最關鍵的一環就是剛度矩陣的獲得，如果涉及到模態分析，還有質量矩陣。考慮到商業軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來看自己編寫的程序是否正確，因此如何提取ANSYS中結構的剛度矩陣，并進行隨后的驗證或者二次開發是一個問題。
https://www.yqgqt.org.cn/content/post/1796144
受上述帖子的啟發，使用MATLAB提取ANSYS中的剛度和質量矩陣，并進行模態分析驗證提取的矩陣的正確性。
首先，在ANSYS中使用HBMAT命令方法提取整體矩陣。
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---輸出矩陣的路徑和文件名，缺省為當前工作路徑和當前工作文件名。
ext---輸出矩陣文件的擴展名，缺省為.matrix。
form---定義輸出矩陣文件的格式，其值可取：
=ASCII：ASCII碼格式；
=BIN：二進制格式。
matrix---定義輸出矩陣的類型，其值可取：
=STIFF：輸出剛度矩陣。可用于寫入了.FULL文件的任何類型的分析。
=MASS：輸出質量矩陣。可用于特征值屈曲、子結構分析、模態分析。
=DAMP：輸出阻尼矩陣。僅用于有阻尼的模態分析。
rhs---右邊項輸出控制（右邊項指用矩陣所表示方程的等號右端矢量，這里可為節點荷載向量），如rhs=YES則輸出，如rhs=NO則不輸出。
模態分析時，因僅LANB和QR法可生成完整的質量矩陣，因此也僅采用這兩種方法時才可使用HBMAT命令得到質量矩陣文件。
展開 ANSYS中整體、單元剛度和質量矩陣的提取
將輸出信息送到cp.out文件
/debug,-1,,,1 ! 指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
ANSYS模型剛度、質量矩陣快速提取小軟件—km_from_Ansys ￥88
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---輸出矩陣的路徑和文件名，缺省為當前工作路徑和當前工作文件名。
ext---輸出矩陣文件的擴展名，缺省為.matrix。
form---定義輸出矩陣文件的格式，其值可取：
=ASCII：ASCII碼格式；
=BIN：二進制格式。
matrix---定義輸出矩陣的類型，其值可取：
=STIFF：輸出剛度矩陣。可用于寫入了.FULL文件的任何類型的分析。
=MASS：輸出質量矩陣。可用于特征值屈曲、子結構分析、模態分析。
=DAMP：輸出阻尼矩陣。僅用于有阻尼的模態分析。
rhs---右邊項輸出控制（右邊項指用矩陣所表示方程的等號右端矢量，這里可為節點荷載向量），如rhs=YES則輸出，如rhs=NO則不輸出。
模態分析時，因僅LANB和QR法可生成完整的質量矩陣，因此也僅采用這兩種方法時才可使用HBMAT命令得到質量矩陣文件。
這個是Ansys幫助文件中對HBMAT命令的解釋，在Ansys分析中，會在根目錄下形成.full的二進制文件，里面存儲的就是已經組合好的剛度和質量矩陣，只是是以二進制文件的形式存在，通過HBMAT命令可以將二進制文件轉換為可以使用的十進制文件，十進制文件是以txt的形式存放在根目錄下。
打開生成的十進制文件，發現并沒有得到我們所想象的剛度、質量矩陣的形式，而是一長串的數字。
展開 
ANSYS剛度矩陣的提取與解析(python解析）
就ansys如何提取剛度矩陣、如何解讀提取的文檔以及利用Python進行解析。
在workbench中實現整個過程的參數化過程除了前幾次文章介紹的模型與網格，還應該包括材料參數的參數化定義。利用Python進行二次開發能夠實現材料參數的自由定義，比如來源于excel表格或者文檔的數據，通過Python代碼的自動讀取，參與到實際的有限元分析進程中。
結構有限元最后的求解過程總是歸結到求解一個大型矩陣方程Ax=b，對于一些情況還需要考慮質量矩陣M和阻尼矩陣C。有限元程序在組裝完所有單元的剛度矩陣后，考慮模型所施加的約束和載荷，最終將剛度矩陣進行一些處理，例如乘大數法，變成Ax=b的形式，其中A是剛度矩陣，b是節點載荷，x為待求的節點位移，A和b全為已知量。
基本上各類有限元軟件均能夠提取模型的剛度矩陣，此次針對剛度矩陣的提取與解析做一個例子，采用的軟件是ANSYS經典。
在ANSYS中建立一個簡單的模型，劃分網格后共12個節點，定義材料參數，施加約束和載荷后求解。有限元模型如下所示。
待求解結束后，會在工作目錄下生成一個后綴為full的文件，之后即可進行剛度矩陣的提取。
通過主菜單，如下所示。
選擇Matrix后，彈出如下所示的界面。
其中，File to be read需要指定工作目錄下生成的full文件，Name of file to write為所導出剛度矩陣的文件名稱；Output matrix file format表示文件格式，還有Binary，生成的是文檔文件，選擇Ascii即可；Matrix to write表示輸出的是剛度矩陣/質量矩陣還是阻尼矩陣；RHS選項表示是否同時輸出右端項，也即是Ax=b中的b。
打開生成的剛度矩陣文檔，如下所示。
展開 『分享』ANSYS中整體、單元剛度和質量矩陣的提取
將輸出信息送到cp.out文件
/debug,-1,,,1 ! 指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩陣
