ansys 提取載荷矩陣的案例
如何從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）。
展開 提取ANSYS中的剛度矩陣和質量矩陣到MATLAB中進行二次開發 ￥88
模態分析完成后，在ansys中運行HBMAT命令得到十進制文件：
! 獲取矩陣十進制文件
/aux2
file,Frame,full !與文件名一致
hbmat,Stiff,txt,,ascii,stiffness,yes,yes !提取剛度矩陣 存儲為Stiff.txt
hbmat,Mass,txt,,ascii,mass,yes,yes !提取質量矩陣 存儲為Mass.txt
Finish
提取其前10階模態如下：
通過MATLAB對文件提取并處理后，可以使用spy命令查看稀疏矩陣的分布，剛度矩陣分布和各變量維度如下：
稀疏剛度矩陣的具體分布如下：
稀疏質量矩陣的具體分布如下：
計算矩陣的特征值，利用提取的剛度矩陣和質量矩陣做模態分析。
與ANSYS中的結果一致，因此驗證了正確性。
在Stiff.mapping文件中，格式是這樣的：
在MATLAB中使用相同的方法提取并整理，如下：
展開 ANSYS模型剛度、質量矩陣快速提取小軟件—km_from_Ansys ￥88
驗證
Matlab打開生成的矩陣文件，一共有三個矩陣，分別是kk,mm,address,其中kk文件和mm文件分別是結構剛度和質量矩陣文件，以稀疏矩陣形式存儲，address文件是矩陣的自由度說明，其中第二列是對應節點號，第三列是對應自由度編碼，1,2,3,4,5,6,分別代表X,Y,Z,ROTX,ROTY和ROTZ。便于對矩陣進行操作。
通過求取模態特征值來驗證提取矩陣的有效性：
>>load('KM.mat')
>> [v,d]=eigs(kk,mm,10,'sm');
>> d=sqrt(d);
>> d=d/2/pi;
>> d=sum(d);
>> d=d';
上圖為Ansys中模態分析的結果，下圖為提取矩陣的分析結果，結果一致，表明提取矩陣有效。
6. 說明
提取工具適用于任意單元，任意結構形式，但僅限于彈性模型。
展開 
ansys Workbench螺栓載荷提取時，如何計算載荷偏心距離（VDI2230） ￥10
問題：
VDI2230關于螺栓的計算中對于螺栓載荷的提取沒有過多的涉及，本文針對偏心載荷的提取問題進行簡單說明。
VDI2230中，對于載荷偏心距a的定義如下，虛擬軸線到截面彎矩為0的點之間的距離。
對于實際螺栓連接問題，幾何結構和載荷狀態復雜多變，使用經驗公式估計并不理想。本文介紹使用有限元仿真的方法確定載荷偏心距離。
示例：
以VDI2230中的案例5為例進行對比計算，依據案例5的幾何信息創建仿真模型。
約束筒體底面，在內表面施加20Mpa壓力載荷，同時給螺栓施加約150KN的預緊力（加不加結果變化不大），連接面設定為摩擦面。
將兩個側面設定為，frictionless Support，等效對稱邊界。（這里沒有使用圓周循環對稱邊界，是因為圓周對稱邊界不能支持截面彎矩提取）
注意，在輸出控制中 打開“Nodal Forces”，用于端蓋截面的彎矩提取。
計算完成后，在結果提取中，插入Probe——Moment Reaction——使用surface類型進行端蓋截面彎矩載荷的提取，這里只需要關注X軸彎矩。
依次變更截面位置，就可以獲得一條彎矩隨位置變化的曲線，讀取彎矩為0位置的距離值，再進一步處理加上螺栓偏心距Ssym，就可以換算到載荷偏心距a。
個人認為仿真結果17.535，除了在循環對稱設置上與案例給出條件不同外，其余均能反應案例邊界。
補充案例：
以機械設計手冊兩端固支梁，在均布載荷下的反彎點計算模型為例進行驗證。
仿真結果
公式計算值42.2mm，仿真結果42.23mm。
展開 ANSYS中整體、單元剛度和質量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
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中整體、單元剛度和質量矩陣的提取
、整體剛度和質量矩陣的提取。
該功能需要進行二次開發，由ansys形成
的二進制文件.full提取整體剛度和質量
矩陣。
基于ansys的一個用戶開發程序例子編
了一個程序（附件中）。
開發環境：compaq fortran 6.5
運行環境：win2000。
一個主文件：self.for，
另一個文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫文件，將其
引入項目中（也可直接扔進debug目錄）
.full文件由子空間迭代模態分析獲得。
運行編譯后的可執行文件.exe
2、單元剛度和質量矩陣的提取。
/DEBUG命令。詳細說明可由以下轉載文章
finish
/clear
PI=3.1415926
w1=3
w2=10
w3=6
w4=1.2
r=.8
t=0.08
/PREP7
!*
ET,1,SHELL63
R,1,t
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all
NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, !
展開 Ansys Workbench提取螺栓連接面載荷方法記錄 ￥10
問題：
在使用理論方法對螺栓強度進行評估時，需要輸入螺栓所受的載荷作為計算輸入。螺栓載荷在復雜工況下，通常使用有限元仿真的方式進行模擬。此時需要準確提取螺栓位置的載荷大小用后續理論校核。
示例：
如下圖所示，兩個零件一端鉸接一端使用螺栓連接。在螺栓側端面施加2000N載荷（無螺栓預緊力）。需要提取螺栓在連接面處所受到的載荷包括：力和力矩。
載荷提取結果：
1.螺栓連接面位置作用力
2.螺栓連接面位置因載荷分布不均產生的彎矩
詳細步驟：
1.螺栓連接面位置的載荷提取，需要在結果輸出中打開節點力輸出項“Nodal Forces-Yes”
2.需要在螺栓連接面位置創建局部坐標系和虛擬結構面
展開 ANSYS經典提取螺栓軸向載荷的方法 ￥10
Beam188軸向力的提取方法
