ansys中單元剛度矩陣、總剛矩陣的提取

《ansys中單元剛度矩陣、總剛矩陣的提取》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、Ansys中單元剛度矩陣、總剛矩陣的提取一、單元剛度矩陣的提取/DEBUG命令詳細(xì)說明：finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX,,,2e11UIMP,1,NUXY,,,0.3,UIMP,1,DAMP,,,0.2,UIMP,1,DENS,,,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMES

2、H,allNSEL,S,LOC,X,0.0D,all,,,,,,ALL,,,,,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,,!將輸出信息送到cp.out文件/debug,-1,,,1!指定輸出單元矩陣/SOLUSOLVEfinish/OUTPUT,TERM!將輸出信息送到outputwindows中這時用編輯器打開cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩陣二、整體剛度矩陣的提?。ㄓ腥N方法：用戶程序法、超單元法、HBMAT命令法）1、用戶程序法：需要二次開發(fā)（略）

3、2、超單元法/soluantype,7!substructuring分析類型seopt,matname,1!設(shè)置文件名稱和剛度矩陣類型(剛度，質(zhì)量，阻尼等)nsel,all!選擇所有節(jié)點m,all,all!定義所有節(jié)點自由度為主自由度solve!求解selist,matname,3!列出整體剛度矩陣3、HBMAT命令法提取整體矩陣命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---輸出矩陣的路徑和文件名，缺省為當(dāng)前工作路徑和當(dāng)前工作文件名。ext---輸出矩陣文件的擴展名，缺省為.mat

4、rix。form---定義輸出矩陣文件的格式，其值可取：=ASCII：ASCII碼格式；=BIN：二進制格式。matrix---定義輸出矩陣的類型，其值可?。?STIFF：輸出剛度矩陣。可用于寫入了.FULL文件的任何類型的分析。=MASS：輸出質(zhì)量矩陣?？捎糜谔卣髦登?、子結(jié)構(gòu)分析、模態(tài)分析。=DAMP：輸出阻尼矩陣。僅用于有阻尼的模態(tài)分析。rhs---右邊項輸出控制（右邊項指用矩陣所表示方程的等號右端矢量，這里可為節(jié)點荷載向量），如rhs=YES則輸出，如rhs=NO則不輸出。模態(tài)分析時，因僅LANB和QR法可生成完整的質(zhì)量

5、矩陣，因此也僅采用這兩種方法時才可使用HBMAT命令得到質(zhì)量矩陣文件。⑵Harwell-Boeing文件格式用HBMAT命令可輸出結(jié)構(gòu)剛度矩陣、質(zhì)量矩陣和阻尼矩陣，其文件記錄格式為大型稀疏矩陣的標(biāo)準(zhǔn)交換格式，采用索引存儲方法僅記錄矩陣的非零元素。文件基本格式是前面有4或5行描述數(shù)據(jù)，其后為單列矩陣元素值，說明如下：第1行：格式（A72），為文件頭的字符型解釋，如剛度矩陣或質(zhì)量矩陣等標(biāo)題。第2行：格式（5I14），分別表示該文件的總行數(shù)（不包括文件頭）、矩陣列指針的總行數(shù)、矩陣行索引的總行數(shù)、矩陣元素數(shù)值的總行數(shù)、右邊項總行數(shù)。第

6、3行：格式（A3,11X,4I14），分別為矩陣類型、矩陣行數(shù)、矩陣列數(shù)、矩陣行索引數(shù)（對組裝后的矩陣，該值等于矩陣行索引數(shù)）、單元元素數(shù)（對組裝后的矩陣此值為0）。第4行：格式（2A16,2A20），分別表示列指針格式、行索引格式、系數(shù)矩陣數(shù)值格式、右邊項數(shù)值格式。第5行：格式（A3,11X,2I14），A3各列分別表示右邊項格式、應(yīng)用高斯起始矢量、應(yīng)用eXact求解矢量；兩個整數(shù)分別表示右邊項列數(shù)、行索引數(shù)。三個字符中的第1個字符可?。篎---全部存貯（如節(jié)點荷載向量的全部元素）、M---與系數(shù)矩陣相同方法。第6行后：矩陣元

7、素值（單列）。矩陣類型用3個字符表示，第1個字符可?。篟---實數(shù)矩陣、C---復(fù)數(shù)矩陣、P---僅矩陣結(jié)構(gòu)（無元素數(shù)值）；第2個字符可?。篠---對稱矩陣、U---不對稱矩陣、H---Hermitian矩陣、Z---病態(tài)對稱矩陣；R---帶狀矩陣；第3個字符可?。篈---組裝的矩陣、E---單元矩陣（未組裝）。對稱矩陣只存儲下三角元素，如結(jié)構(gòu)剛度矩陣為對稱矩陣，Harwell-Boeing格式則僅記錄下三角元素。根據(jù)Harwell-Boeing文件格式，可讀取矩陣的任意行列元素的數(shù)值，也可編程還原為滿矩陣存儲，以便它用，很顯然

8、這種提取方式比較方便。如當(dāng)生成.FULL文件后，可采用命令/AUX2$FILE,mywork,full$HBMAT,mystiff,txt,ASCII,STIFF,YES$FINISH將二進制mywork.full文件輸出為ASCII碼文件mystiff.txt