ansys形變和剛度的案例
ANSYS中整體、單元剛度和質(zhì)量矩陣的提取
這時(shí)用編輯器打開(kāi)cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元剛度和質(zhì)量矩陣的提取.rar
『分享』ANSYS中整體、單元剛度和質(zhì)量矩陣的提取
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫(kù)文件，將其
引入項(xiàng)目中（也可直接扔進(jìn)debug目錄）
.full文件由子空間迭代模態(tài)分析獲得。
運(yùn)行編譯后的可執(zhí)行文件.exe
2、單元剛度和質(zhì)量矩陣的提取。
/DEBUG命令。詳細(xì)說(shuō)明可由以下轉(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,, !
展開(kāi) 提取ANSYS中的剛度矩陣和質(zhì)量矩陣到MATLAB中進(jìn)行二次開(kāi)發(fā) ￥88
最近在考慮自己編寫的程序和商用軟件的驗(yàn)證問(wèn)題，有限元結(jié)構(gòu)分析中最關(guān)鍵的一環(huán)就是剛度矩陣的獲得，如果涉及到模態(tài)分析，還有質(zhì)量矩陣。考慮到商業(yè)軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來(lái)看自己編寫的程序是否正確，因此如何提取ANSYS中結(jié)構(gòu)的剛度矩陣，并進(jìn)行隨后的驗(yàn)證或者二次開(kāi)發(fā)是一個(gè)問(wèn)題。
https://www.yqgqt.org.cn/content/post/1796144
受上述帖子的啟發(fā)，使用MATLAB提取ANSYS中的剛度和質(zhì)量矩陣，并進(jìn)行模態(tài)分析驗(yàn)證提取的矩陣的正確性。
首先，在ANSYS中使用HBMAT命令方法提取整體矩陣。
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---輸出矩陣的路徑和文件名，缺省為當(dāng)前工作路徑和當(dāng)前工作文件名。
ext---輸出矩陣文件的擴(kuò)展名，缺省為.matrix。
form---定義輸出矩陣文件的格式，其值可取：
=ASCII：ASCII碼格式；
=BIN：二進(jìn)制格式。
matrix---定義輸出矩陣的類型，其值可取：
=STIFF：輸出剛度矩陣。可用于寫入了.FULL文件的任何類型的分析。
=MASS：輸出質(zhì)量矩陣。可用于特征值屈曲、子結(jié)構(gòu)分析、模態(tài)分析。
=DAMP：輸出阻尼矩陣。僅用于有阻尼的模態(tài)分析。
rhs---右邊項(xiàng)輸出控制（右邊項(xiàng)指用矩陣所表示方程的等號(hào)右端矢量，這里可為節(jié)點(diǎn)荷載向量），如rhs=YES則輸出，如rhs=NO則不輸出。
模態(tài)分析時(shí)，因僅LANB和QR法可生成完整的質(zhì)量矩陣，因此也僅采用這兩種方法時(shí)才可使用HBMAT命令得到質(zhì)量矩陣文件。
展開(kāi) 