ansys整體剛度矩陣
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時間：2023-03-08
ansys整體剛度矩陣的實例教程
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元剛度和質(zhì)量矩陣的提取.rar
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩陣
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 !設(shè)置文件名稱和剛度矩陣類型(剛度，質(zhì)量，阻尼等)
nsel,all !選擇所有節(jié)點
m,all,all !定義所有節(jié)點自由度為主自由度
solve !求解
selist,matname,3 !列出整體剛度矩陣
展開 、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個用戶開發(fā)程序例子編
了一個程序（附件中）。
開發(fā)環(huán)境：compaq fortran 6.5
運行環(huán)境：win2000。
一個主文件：self.for，
另一個文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫文件，將其
引入項目中（也可直接扔進(jìn)debug目錄）
.full文件由子空間迭代模態(tài)分析獲得。
運行編譯后的可執(zhí)行文件.exe
2、單元剛度和質(zhì)量矩陣的提取。
/DEBUG命令。詳細(xì)說明可由以下轉(zhuǎn)載文章
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,, !
展開 在上面討論中已經(jīng)指出，有限元求解方程的系數(shù)矩陣具有奇異性，必須引入適當(dāng)?shù)奈灰萍s束條件，以消除這種奇異性，亦即消除彈性體的剛體位移。消除了整體剛度矩陣的奇異性后，才能從方程組（32）求解結(jié)點位移。在一般情況下，所考慮問題的邊界往往已有一定的約束條件，排除了剛體運動的可能性。否則，可適當(dāng)指定某些結(jié)點的位移值，以避免計算機(jī)存儲作大的更動。下面就介紹兩種比較簡單的引入已知結(jié)點位移的方法。
1、對角元素改l法
這種方法是把結(jié)點的指定值置入方程給(32)，保持方程仍是2nX2n階，而將K和P修正。例如，若指定結(jié)點i在y方向位移vi的值，則令K中的元素Ki,i為1，而第i行和i列的其余元素都為零。P中的第i個元素則用位移v的已知值代入，P中的其他各行元素都減去結(jié)點位移的指定值和原來K中這行的相應(yīng)行元素的乘積。
為了說明這一引進(jìn)結(jié)點已知位移的過程，我們來考察下面只有四個方程的簡單例子。方程(32)展開成如下的形式
設(shè)這個系統(tǒng)中結(jié)點位移u1和u2被指定為
當(dāng)引用上述方法后，方程(50)就變成
然后，就用這組維數(shù)不變的方程來求解所有的結(jié)點位移。顯然，其解答為u1=β1、u2=β3；
v1、v2仍為原方程的解答。
這種方法最適用于給定零位移，此時除將給定的零值位移修改對應(yīng)的載荷陣元(如例中令Px1=0，Px2=0)外，其他載荷陣中的元素不必作任何修正。
2、對角元素乘大數(shù)法
此法是將K中與指定結(jié)點位移有關(guān)的主對角元素乘上一個大數(shù)，例2xl015，同時將p的
對應(yīng)元素?fù)Q上結(jié)點位移指定值與同一個大數(shù)的乘積。實際上，這種方法就是使得K中相應(yīng)行的修正項遠(yuǎn)大于非修正項。
展開 1.引論
經(jīng)常使用Ansys、Abaqus等一系列有限元分析軟件進(jìn)行計算、學(xué)習(xí)的學(xué)生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質(zhì)量矩陣的重要性。但是由于軟件的黑盒性質(zhì)，大家往往在實際使用十分成熟的商業(yè)化軟件的過程中慢慢忽視了有限元及其衍生出的商業(yè)軟件背后的原理與方法。
這時，不管是在學(xué)習(xí)中還是在工程應(yīng)用中往往都會遇到一個同樣的問題，那么就是如何將Ansys APDL運行中的產(chǎn)生的各種數(shù)據(jù)（例如：剛度矩陣、質(zhì)量矩陣）導(dǎo)出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實際工作或?qū)W習(xí)中需要用到此類技能的同學(xué)、同事們，讓大家更了解Ansys APDL背后的工作原理與數(shù)據(jù)導(dǎo)出方式。
當(dāng)然，在社區(qū)中早就有大佬回答過了這個問題，并給大家制作了相應(yīng)的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質(zhì)量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進(jìn)一步學(xué)習(xí)了解軟件背后機(jī)理的群體，并在此基礎(chǔ)上保留教學(xué)的簡潔性，提供導(dǎo)出矩陣與轉(zhuǎn)換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因為優(yōu)化不完全導(dǎo)致的運行bug。
2.有限元軟件導(dǎo)出剛度矩陣與質(zhì)量矩陣的方法
在使用APDL進(jìn)行求解時，每次在求解完成后都會在工作路徑下生成一個.full文件，而這個文件十分關(guān)鍵，其正是剛度矩陣與質(zhì)量矩陣的所在之處。
展開 
ansys整體剛度矩陣的相關(guān)專題、標(biāo)簽、搜索
ansys整體剛度矩陣的最新內(nèi)容
在有限元分析中，ANSYS 可以導(dǎo)出大規(guī)模稀疏矩陣（如剛度矩陣、質(zhì)量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對后續(xù)二次開發(fā)、動力學(xué)分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個 MATLAB 函數(shù)，可直接從 ANSYS 導(dǎo)出的 HB 矩陣文件中讀取并重構(gòu)成 MATLAB 稀疏矩陣：
1.引論
經(jīng)常使用Ansys、Abaqus等一系列有限元分析軟件進(jìn)行計算、學(xué)習(xí)的學(xué)生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質(zhì)量矩陣的重要性。但是由于軟件的黑盒性質(zhì)，大家往往在實際使用十分成熟的商業(yè)化軟件的過程中慢慢忽視了有限元及其衍生出的商業(yè)軟件背后的原理與方法。
這時，不管是在學(xué)習(xí)中還是在工程應(yīng)用中往往都會遇到一個同樣的問題，那么就是如何將Ansys
最近在考慮自己編寫的程序和商用軟件的驗證問題，有限元結(jié)構(gòu)分析中最關(guān)鍵的一環(huán)就是剛度矩陣的獲得，如果涉及到模態(tài)分析，還有質(zhì)量矩陣。考慮到商業(yè)軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來看自己編寫的程序是否正確，因此如何提取ANSYS中結(jié)構(gòu)的剛度矩陣，并進(jìn)行隨后的驗證或者二次開發(fā)是一個問題。
https://www.yqgqt.org.cn/content/post/1796144
1. 背景
從事結(jié)構(gòu)振動控制、車橋耦合振動、結(jié)構(gòu)健康監(jiān)測傳感器優(yōu)化布置、結(jié)構(gòu)動力性能分析等等一系列研究的同仁們應(yīng)該都面臨過一個同樣的問題—“怎么把結(jié)構(gòu)的剛度和質(zhì)量矩陣建立出來？”。這對于那些數(shù)值分析高手和專家可能不是什么問題；但是對于科研剛?cè)腴T的新手來說，這個難度還是相當(dāng)大的。如果都靠自己寫程序來建立有限元模型，則對理論基礎(chǔ)、編程水平都有很高的要求，甚至程序做出來也未必能保證其正確性，是一個很讓人頭疼的問題
就ansys如何提取剛度矩陣、如何解讀提取的文檔以及利用Python進(jìn)行解析。
在workbench中實現(xiàn)整個過程的參數(shù)化過程除了前幾次文章介紹的模型與網(wǎng)格，還應(yīng)該包括材料參數(shù)的參數(shù)化定義。利用Python進(jìn)行二次開發(fā)能夠?qū)崿F(xiàn)材料參數(shù)的自由定義，比如來源于excel表格或者文檔的數(shù)據(jù)，通過Python代碼的自動讀取，參與到實際的有限元分析進(jìn)程中。
結(jié)構(gòu)有限元最后的求解過程總是歸結(jié)到求解一個大型矩陣方程
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個用戶開發(fā)程序例子編
了一個程序（附件中）。
開發(fā)環(huán)境：compaq fortran 6.5
運行環(huán)境：win2000。
一個主文件：self.for，
另一個文件matrixout.f90用于矩陣輸出
binlib.lib
1、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個用戶開發(fā)程序例子編
了一個程序（附件中）。
開發(fā)環(huán)境：compaq fortran 6.5
運行環(huán)境：win2000。
一個主文件：self.for，
另一個文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫文件
在上面討論中已經(jīng)指出，有限元求解方程的系數(shù)矩陣具有奇異性，必須引入適當(dāng)?shù)奈灰萍s束條件，以消除這種奇異性，亦即消除彈性體的剛體位移。消除了整體剛度矩陣的奇異性后，才能從方程組（32）求解結(jié)點位移。在一般情況下，所考慮問題的邊界往往已有一定的約束條件，排除了剛體運動的可能性。否則，可適當(dāng)指定某些結(jié)點的位移值，以避免計算機(jī)存儲作大的更動。下面就介紹兩種比較簡單的引入已知結(jié)點位移的方法。
1、對角元素改
1.
Which matrix you would like? element stiffness matrix or full stiffness matrix?
element stiffness is within file.emat. full stiffness matrix is within file.full
A simple way to dump the matrix
