ansys質量提取的案例
如何從Ansys APDL中提取剛度矩陣與質量矩陣？ ￥69
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys APDL運行中的產生的各種數據（例如：剛度矩陣、質量矩陣）導出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實際工作或學習中需要用到此類技能的同學、同事們，讓大家更了解Ansys APDL背后的工作原理與數據導出方式。
當然，在社區中早就有大佬回答過了這個問題，并給大家制作了相應的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進一步學習了解軟件背后機理的群體，并在此基礎上保留教學的簡潔性，提供導出矩陣與轉換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因為優化不完全導致的運行bug。
2.有限元軟件導出剛度矩陣與質量矩陣的方法
在使用APDL進行求解時，每次在求解完成后都會在工作路徑下生成一個.full文件，而這個文件十分關鍵，其正是剛度矩陣與質量矩陣的所在之處。
展開 ANSYS模型剛度、質量矩陣快速提取小軟件—km_from_Ansys ￥88
背景
從事結構振動控制、車橋耦合振動、結構健康監測傳感器優化布置、結構動力性能分析等等一系列研究的同仁們應該都面臨過一個同樣的問題—“怎么把結構的剛度和質量矩陣建立出來？”。這對于那些數值分析高手和專家可能不是什么問題；但是對于科研剛入門的新手來說，這個難度還是相當大的。如果都靠自己寫程序來建立有限元模型，則對理論基礎、編程水平都有很高的要求，甚至程序做出來也未必能保證其正確性，是一個很讓人頭疼的問題。
對于一些簡單的被動控制裝置或簡單的動力學分析，當然也可以在有限元分析軟件中構造出裝置組成直接分析（剛度+阻尼類型），但是對于稍復雜一些的控制裝置和耦合分析等問題，會受到平臺功能上的限值，尤其是對于主動和半主動等涉及控制算法的研究來說，基本很難在有限元軟件平臺上實現分析。再加上如果需要對裝置進行參數優化，需要進行多次重復計算，難度就更大。
Ansys、ABAQUS等軟件平臺給我們提供了比較穩定有效的有限元模型建立平臺，通過借助商業軟件來建立模型，再將其中的剛度、質量矩陣導出，是非常可取的一種方法。如果能夠提取出模型的矩陣，明晰計算原理，就能夠很容易的通過自己的程序設計對計算過程進行補充、調整，來達到自己定制的計算分析目的。其實，不僅對于振動控制，比如結構靜動力分析、車橋耦合分析、結構傳感器優化配置方案設計等，都有應用需求。因此，一個能夠便捷的提取結構矩陣的方法就顯得至關重要。
技術鄰平臺已經有大佬提供了ABAQUS軟件剛度和質量矩陣的導出方法。這里補充一下在ANSYS中導出質量和剛度矩陣的方法和小軟件。
2.
展開 ANSYS中整體、單元剛度和質量矩陣的提取
這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
ANSYS知識普及1——如何提取模態質量（ANSYS專家編輯，非原創，歡迎轉摘）
本人準備出一個ANSYS知識普及系列，將有用的網上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術鄰ANSYS專家
業務咨詢網址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網上；
2、如侵犯知識產權，請聯系ANSYS專家本人或者技術鄰，我將第一時間刪除。
下面是《ANSYS Verification Manual》中VM89.DAT稍加修改后提取模態質量的例子：
/PREP7
/TITLE, VM89, NATURAL FREQUENCIES OF A TWO-MASS-SPRING SYSTEM
C*** VIBRATION THEORY AND APPLICATIONS, THOMSON, 2ND PRINTING, PAGE 163,EX 6.2-2
ET,1,COMBIN14,,,2
ET,2,MASS21,,,4
R,1,200 ! SPRING CONSTANT = 200
R,2,800 ! SPRING CONSTANT = 800
R,3,.5 ! MASS = .5
R,4,1 ! MASS = 1
N,1
N,4,1
FILL
E,1,2 ! SPRING ELEMENT (TYPE,1) AND K = 200 (REAL,1)
TYPE,2
REAL,3
E,2 ! MASS ELEMENT (TYPE,2) AND MASS = .5 (REAL,3)
TYPE,1
REAL,2
E,2,3 !
展開 
基于 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中整體、單元剛度和質量矩陣的提取
、整體剛度和質量矩陣的提取。
該功能需要進行二次開發，由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中的剛度矩陣和質量矩陣到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命令得到質量矩陣文件。
展開 提取整體剛度矩陣、質量矩陣及阻尼矩陣的三種方法
這時用編輯器打開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 !列出整體剛度矩陣
展開 ?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
ANSYS Beam188提取彎矩為例介紹ANSYS定義單元表提取數據 (解決彎矩圖鋸齒狀） ￥20
在ANSYS中有些數據無法直接訪問，需要通過定義單元表完成單元的結果的訪問。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數據的詳細過程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數據，這里以Beam188的彎矩為例
2.1 啟動ANSYS幫助菜單, 在索引框輸入Beam188然后搜索, 在單元輸出介紹找到彎矩的名稱（代號）。
2.2 回到ANSYS界面，比如要輸出Mz, 則需要在添加SMISC,3 和SMISC,16 ,如圖
3. 輸出數據：Main Menu>General Postproc>Element Table> List E T, 選擇前面定義的SMISC,3 和SMISC,16 輸出單元I和J節點的Mz數值，如圖
4. 顯示彎矩云圖：Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res, 這里要注意要在LabI 選SMISC,3 LabJ 選SMSCI,16。
輸出彎矩到這就結束了，小編突然發現，輸出的彎矩值在每個單元的I和J處是一樣的（Beam188為2節點單元），彎矩圖也就成了鋸齒形，于是去問了度娘一波，各路盆友給出解決方法，然而并沒有起作用的，于是乎我又想起來了“幫助文檔大法”，于是認認真真將Beam188的幫助文檔閱讀了一遍，功夫不負有心人，最終。。。
展開 ANSYS的get命令常用操作（信息提取和結果結果提取）
ANSYS的get命令常用操作（信息提取和結果結果提取）
在ANSYS分析過程中，*get命令作為一個提取信息的常用命令，作用非常大，不管是在前處理、求解還是后處理過程中，都能夠有發揮的空間，尤其是后處理過程，對結果的批量輸出來說不可缺少。
*get能夠提取的信息相當多，其命令語句如下：
*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM
Par：定義的變量名稱，用于存儲提取的數據；
Entity：關鍵字，是信息提取的對象，包括NODE, ELEM, KP, LINE, AREA, VOLU, PDS等；
ENTNUM：當前對象的數字標識，比如節點的節點號，單元的單元號等；
Item1：提取的信息，可用的非常多，后面展開；
IT1NUM：和Item1配合使用。
由于*get的功用實在太多，不就一一列舉，單就常用的枚舉。
展開 
ANSYS workbench中質心的提取方法 ￥88
在ansys中，我們經常需要提取模型的質心，查看其質心的位置，以及動力學中查看質心的移動軌跡，并非某一參考點的移動軌跡。在ansys軟件中很容易查看模型，或者某一組件的質心如圖所示。只需要在幾何模型中選中模型，然后再property中查看坐標即可。但是這種為靜態結果，也就是初始狀態的結果
如果要查看模型在受力作用下發生了變形，那么對應的質心顯然會發生移動。但是在這種情況下怎么查看質心呢。一個近似的方法可以如下方式：
在后處理中插入自定義的結果loc-defy，可以得到模型變形后的坐標，然后查看average的結果，就是質心的結果，該方法根據變形量的平均效果來得到質心的，如下所示
但是真正的質心提取方法并非上面的方式，需要準確獲得，需要采用下面的命令
upcoord,1更新模型單元位置
展開 ANSYS Workbench中如何提取截面內力 ￥3.9
在土木及水利設計中，截面內力是結構設計過程中極為重要的參數，也是結構穩定性的重要依據。本文重點介紹如何在Workbench平臺自定義截面并獲得相應截面的內力，并將其結果輸出。方法簡單，操作易上手！最終結果顯示如下：
具體步驟為：1、自定義創建截面，這里建議采用局部坐標系的方法建立截面位置；
ANSYS如何提取某一節點的應力時程 ￥100
在對結構進行時程分析后，我們經常提取的是全時程最大von Mises stress。
那么如何提取某一個節點的von Mises stress呢？
首先明確ANSYS的節點附加在單元上，可以通過選擇單元上節點的方法提取節點應力。
1 確定節點所在單元，顯示節點編號。
例單元號8560，節點號8678。
2 進入TimeHist Postpro, 定義變量。
3變量顯示。
付費內容為相關命令流。
ansys Workbench螺栓載荷提取時，如何計算載荷偏心距離（VDI2230） ￥10
使用上述截面力矩提取方法，計算VDI2230中的初始條件，載荷偏心距a的插件，使用介紹。
插件主要實現，
1. 自動循環移動截面位置，提取X軸彎矩數據；
2. 繪制彎矩曲線圖；
3. 插值計算彎矩0點位置；
4. 在零點附近增補提取截面，精確插值結果;
