ansys單元質量提取的案例
ANSYS中整體、單元剛度和質量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
『分享』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 APDL中提取剛度矩陣與質量矩陣？ ￥69
2.1命令流方法
（1）用戶程序法
（2）超單元法
（3）Hbmat法
這三種在網絡上都能很輕易的查到其使用的命令流與方法，但前兩種并不推薦大家使用，不過感興趣的朋友可以自行搜索盡可能多的了解一些。在這里我主要給大家介紹的方法是Hbmat法，也是我推薦大家使用的方法。其可以直接導出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方法。
從而我們便能導出所需要的剛度矩陣與質量矩陣源文件了。
展開 ?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模型剛度、質量矩陣快速提取小軟件—km_from_Ansys ￥88
提取方法
Ansys中其實提供了幾種單元剛度矩陣和整體剛度矩陣的提取接口，先做一下簡單介紹：
（1）單元剛度矩陣提取，集成整體剛度矩陣
首先選擇要提取的單元，通過/debug命令輸出選擇集中的單元剛度矩陣，再根據結構的單元連接信息，組合成整體剛度矩陣，命令如下：
/OUTPUT,cp,out,, ! 將輸出信息送到cp.out文件
/debug,-1,,,1 ! 指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
感受：這個方法很直接，但稍顯麻煩，需要提取到結構的單元組成信息、節點信息以及單元剛度矩陣等，對于稍大一些結構來說，可能需要1個小時甚至更長時間才能完成。
（2）超單元方法
/solu
antype,7 !substructuring分析類型
seopt,matname,1 !設置文件名稱和剛度矩陣類型(剛度，質量，阻尼等)
nsel,all !選擇所有節點
m,all,all !定義所有節點自由度為主自由度
solve !求解
selist,matname,3 !列出整體剛度矩陣
這種方法是可以提取到結構整體剛度矩陣的，但是問題在于需要指定主自由度，對于一般結構而言，這個方法還是適用的，問題是如果結構中存在耦合關系，結構自由度存在依賴關系，該方法就不太適用了，可能還有其他解決辦法，但是我后來還是放棄這種方法。
（3）HBMAT命令方法提取整體矩陣
這種方法的適用性比較廣，也是個人比較推崇的一種方法。
展開 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知識普及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中的剛度矩陣和質量矩陣到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命令得到質量矩陣文件。
展開 hypermesh-ansys聯合仿真之質量單元
質量單元屬于0維單元，ANSYS提供了質量單元mass21，該單元有6個自由度，3個平動自由度和3個繞軸的轉動自由度，可以分別設置不同方向上的不同質量和轉動慣量，但是一般3個平動方向上的質量是相同的，而3個轉動方向上的轉動慣量可能分別不同。轉動慣量可能對某些非轉動模態影響較小甚至可以忽略，但是對某些模態影響比較明顯，所以在較容易獲得部件轉動慣量的情況下盡量將部件簡化為質量單元時輸入每個方向上的轉動慣量參數。
質量單元的另一個功能是作為輔助單元使用，在利用hypermesh為ANSYS求解器建模前處理時，涉及到不同零部件單元之間的連接裝配，此時在一些連接單元的節點上需要安放一個單元才能在導入ANSYS計算時正常進行，下面舉實例說明。
上圖是一個板通過4個紫色的柱焊接在板的4個孔上，建模是通過CERIG單元將板與柱在焊接位置剛性連接，然后在4個柱的頂端安裝在其他部件上，這里將柱的頂端連接到同一個節點上（節點號為4417），然后在該節點上施加固定約束邊界條件。建好模型后導出CBD文件并讀入ANSYS進行模態求解，開始求解時報出如下圖錯誤。
報錯信息顯示為，約束方程1有未使用的節點4417。主要原因是hypermesh中的CERIG單元轉化到ANSYS是約束方程。在建立節點耦合時，比如將若干單元的節點自由度耦合到一個新建的節點時，這個新建節點比如依附于某個單元，否則求解時就會報出上述錯誤信息，這里的解決方案就是在節點4417處建立一個mass21單元，為了消除mass21單元對求解結果的影響需將mass21的質量屬性設置到非常低，特別是在模態求解時，質量會嚴重影響模態求解結果，效果如下圖。
展開 ANSYS提取具體三維單元的體積，面元的面積和線的長度
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號為100的線的長度
以上對應的GUI操作： Utility Menu>Parameters>Get Scalar Data 如果要一次性提取多個元素的相關參數，可以用命令 *VGET, ParR, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM, KLOOP 對應GUI操作：Utility Menu>Parameters>Get Array Data
輸入命令 alist,p 出選項框，選你要看的那個面積，提取選中的單元面積。
*cfopen,'area','txt',
*GET,MaxEleNum,ELEM,,NUM,MAX
*GET,MinEleNum,ELEM,,NUM,MIN
*do,i,MinEleNum,MaxEleNum,1
*if,esel(i),eq,1,then
*get,volu,elem,i,volu
*vwrite,i,volu
(f5.0,f15.12)
*end if
*enddo
*cfcols
展開 
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸 ￥20
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸
問題：
在FKM關于結構疲勞評估計算方法中指出：零部件特征尺寸，影響疲勞結果評估。原因是材料的應力壽命曲線是由標準試樣進行試驗測試獲得的。當零部件的特征尺寸與測試樣件不一致時，需要考慮零部件的特征尺寸這一因素。（一般而言，當零部件的尺寸大于材料標準測試樣件時，零部件的表面或內部缺陷發生的概率會增加，從而導致零部件尺寸越大，疲勞壽命越低）
對與規則幾何形狀的零部件，有相應的經典公式提供特征尺寸的計算；例如圓形細長桿的特征尺寸是直徑；薄板零部件的特征尺寸是板厚等；但是實際工作中的零部件幾何形狀千差萬別，沒有統一的經典公式可以提供特征尺寸的計算；在FKM手冊中給出了一個通用公式，用于估計零部件疲勞危險區域的局部特征尺寸；
FKM關于循環載荷的疲勞評估中，提及可以使用循環載荷下的有限元應力結果進行疲勞損傷估計。此時，除了需要由應力結果估計危險疲勞區域，提取危險點的應力結果外，還需要給出危險疲勞區域的特征尺寸。在Ansys Workbench中，用戶可以方便的查看應力結果云圖，從而大體評估出危險疲勞區域。并且用戶可以通過選取高應力區域的單元體，再通過特征尺寸一般計算公式，來估計高應力區域的特征尺寸，進行進行合理的FKM疲勞評估。
但是，Ansys Workbench中，當用戶選中了某個/某些體單元后，在選擇信息欄中并不能直接給出單元體積和表面的有效信息輸出。并且通過查詢資料，即使在APDL經典界面中對與體單元也是僅僅只能輸出體積（沒有體單元表面的輸出）；并且對與FKM特征尺寸的一般計算公式中，關于表面積A，也并不是指每個體單元所有面的表面積的總和。
展開 