ansys單元質量提取
關注創建者：王靖雯 創建時間：2023-03-07
ansys單元質量提取的實例教程
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
、整體剛度和質量矩陣的提取。
該功能需要進行二次開發，由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,, !
展開 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提取單元的剛度矩陣和質量矩陣：
*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中其實提供了幾種單元剛度矩陣和整體剛度矩陣的提取接口，先做一下簡單介紹：
（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單元質量提取的相關專題、標簽、搜索
ansys單元質量提取的最新內容
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸
問題：
在FKM關于結構疲勞評估計算方法中指出：零部件特征尺寸，影響疲勞結果評估。原因是材料的應力壽命曲線是由標準試樣進行試驗測試獲得的。當零部件的特征尺寸與測試樣件不一致時，需要考慮零部件的特征尺寸這一因素。（一般而言，當零部件的尺寸大于材料標準測試樣件時，零部件的表面或內部缺陷發生的概率會增加
在有限元分析中，ANSYS 可以導出大規模稀疏矩陣（如剛度矩陣、質量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續二次開發、動力學分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數，可直接從 ANSYS 導出的 HB 矩陣文件中讀取并重構成 MATLAB 稀疏矩陣：
1.引論
經常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學習的學生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質量矩陣的重要性。但是由于軟件的黑盒性質，大家往往在實際使用十分成熟的商業化軟件的過程中慢慢忽視了有限元及其衍生出的商業軟件背后的原理與方法。
這時，不管是在學習中還是在工程應用中往往都會遇到一個同樣的問題，那么就是如何將Ansys
最近在考慮自己編寫的程序和商用軟件的驗證問題，有限元結構分析中最關鍵的一環就是剛度矩陣的獲得，如果涉及到模態分析，還有質量矩陣。考慮到商業軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來看自己編寫的程序是否正確，因此如何提取ANSYS中結構的剛度矩陣，并進行隨后的驗證或者二次開發是一個問題。
https://www.yqgqt.org.cn/content/post/1796144
質量單元屬于0維單元，ANSYS提供了質量單元mass21，該單元有6個自由度，3個平動自由度和3個繞軸的轉動自由度，可以分別設置不同方向上的不同質量和轉動慣量，但是一般3個平動方向上的質量是相同的，而3個轉動方向上的轉動慣量可能分別不同。轉動慣量可能對某些非轉動模態影響較小甚至可以忽略，但是對某些模態影響比較明顯，所以在較容易獲得部件轉動慣量的情況下盡量將部件簡化為質量單元時輸入每個方向上的轉動慣量參數
1. 背景
從事結構振動控制、車橋耦合振動、結構健康監測傳感器優化布置、結構動力性能分析等等一系列研究的同仁們應該都面臨過一個同樣的問題—“怎么把結構的剛度和質量矩陣建立出來？”。這對于那些數值分析高手和專家可能不是什么問題；但是對于科研剛入門的新手來說，這個難度還是相當大的。如果都靠自己寫程序來建立有限元模型，則對理論基礎、編程水平都有很高的要求，甚至程序做出來也未必能保證其正確性，是一個很讓人頭疼的問題
在ANSYS中有些數據無法直接訪問，需要通過定義單元表完成單元的結果的訪問。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數據的詳細過程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數據，這里以Beam188
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
在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
本人準備出一個ANSYS知識普及系列，將有用的網上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術鄰ANSYS專家
業務咨詢網址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網上；
