剛度矩陣重組
關注創建者：小白菜瓜 創建時間：2018-05-22
剛度矩陣重組的實例教程
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩陣
3.
其原理很簡單，即使用ansys的超單元即可解決問題。定義超單元，然后列出超單元的剛度矩陣即可。
面是一個小例題，自可明白。
/prep7
k,1
k,2,3000
l,1,2
et,1,beam3
mp,ex,1,2e5
mp,prxy,1,0.3
r,1,5000,2e7,200
lesize,all,,,10
lmesh,all
finish
!----以上正常建立模型，不必施加約束和荷載
/solu
antype,7 !substructuring分析類型
seopt,matname,1 !設置文件名稱和剛度矩陣類型(剛度，質量，阻尼等)
nsel,all !選擇所有節點
m,all,all !定義所有節點自由度為主自由度
solve !求解
selist,matname,3 !列出整體剛度矩陣
展開 其可以直接導出full文件中的矩陣數據，而通過我下面給大家提供的源代碼便可以直接將導出的質量矩陣與剛度矩陣直接轉化為matlab中可使用變量，實現我們的矩陣提取操作。
Hbmat法：
! 提取剛度矩陣
/AUX2
FILE,'file',full ! 將’file’改為自己路徑下的.full文件名
HBMAT, 'Stiffness_mat', dat, , ASCII, STIFF, YES, YES !剛度矩陣
HBMAT, 'Mass_mat', dat, ,ASCII, MASS, YES, YES ! 質量矩陣
FINISH
成功導出后你會在你的工作路徑中看到儲存在Ansys中的剛度矩陣與質量矩陣。
2.2GUI方法
目前大部分此類教程都忽略了最基本的GUI方法，可能是自動帶入了經驗豐富的工程師角色，但對于初學者而言，GUI方法十分有助于理解與學習有限元軟件，因此我也將介紹提取剛度質量矩陣的GUI方法。
從而我們便能導出所需要的剛度矩陣與質量矩陣源文件了。
3.源代碼的使用與轉換（使用方法）
我們打開所導出剛度矩陣會發現里面的數字比較混亂，一時間分不出哪些是我們需要的數據，也沒有辦法直接拿它用來計算，這時，我們便需要使用我們自己的矩陣轉換m文件，其可以自動幫我們提取出剛度矩陣與質量矩陣，并生成matlab中的矩陣變量形式，我們可以直接運用這段源代碼來進行操作。
展開 在有限元分析中，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）。
展開 我記得在patran 中輸出的bdf文件中輸入一條命令，運行一下就可以了。但是記不得了。請教高手，先謝謝了。
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及以后版本。
展開 
剛度矩陣重組的相關專題、標簽、搜索
剛度矩陣重組的最新內容
在有限元分析中，ANSYS 可以導出大規模稀疏矩陣（如剛度矩陣、質量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續二次開發、動力學分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數，可直接從 ANSYS 導出的 HB 矩陣文件中讀取并重構成 MATLAB 稀疏矩陣：
本期給大家推薦一款由木木自研的小工具：TransferMatrix，主要用于導出 Abaqus 中各種剛度矩陣，方便和自研的程序進行實時對比，基于 Pyside6 搭建的軟件界面，用戶僅需導入 inp 文件，就可以自動調用電腦內的 Abaqus 進行計算分析，不受限于 Abaqus 版本，可以導出：
單元剛度矩陣
單元質量矩陣
單元分布節點荷載列陣
整體剛度矩陣
ABAQUS 剛度矩陣11個月前
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys
材料的彈性特征可以使用剛度矩陣來進行完全描述。剛度矩陣中任意元素都可以根據由第一性原理計算得出的聲子色散曲線來確定。
在本次案例研究中，我們介紹了一個用立方體結構來分析硅、金剛石和銅單晶的剛度矩陣的方法。
如圖1，硅是具有金剛石結構的立方晶體，其剛度矩陣如下所示。
由于立方晶體具有對稱性，使用SIESTA模擬軟件計算可得該模型中有3個獨立元素：C11C11，C12C12，和C44C44
基于matlab的實現四節點板單元剛度矩陣求解，振動模態分析。可自主輸入材料參數，板單元長寬厚尺寸。輸出振型結果，輸出多少階可自主設置。最后以可視化的形式展示，程序已調通，可直接運行。
基于matlab的平面桁架結構的總體剛度矩陣計算，最后以圖形形式顯示出桁架結構，程序已調通，可直接運行。
Workbench中剛度矩陣提取方法
Abaqus剛度矩陣提取方法
Optistruct中剛度矩陣提取方法
