ansys輸入矩陣的案例
如何在快速在word文檔中輸入多個(gè)矩陣
如何在快速在word文檔中輸入多個(gè)矩陣
放在文前：假如需要我們快速輸入下面這些公式、矩陣數(shù)據(jù)到word文檔中，大家有什么辦法呢？一個(gè)一個(gè)輸入手動(dòng)進(jìn)去？這個(gè)辦法不是不可以，但是有點(diǎn)低效。使用圖像識(shí)別軟件，直接生成mathtype的代碼格式？那如果沒(méi)有現(xiàn)成的公式可以識(shí)別又該怎么辦呢？請(qǐng)大家接著往下看，筆者為你慢慢道來(lái)。
不知道大家有沒(méi)有這么一種令人痛徹心扉的感受：手頭上有一個(gè)報(bào)告急著要交，但是報(bào)告中含有大量的矩陣數(shù)據(jù)。如果一個(gè)一個(gè)手動(dòng)輸入進(jìn)去，相信大家再做的時(shí)候想死的心都有。那么有沒(méi)有一種辦法可以快速在word中輸入多個(gè)矩陣數(shù)據(jù)呢？
上述數(shù)據(jù)是筆者在編寫(xiě)有限元代碼過(guò)程中生成的數(shù)據(jù)，編程語(yǔ)言為MATLAB。數(shù)據(jù)很容易獲取，但是將數(shù)據(jù)轉(zhuǎn)換為mathtpye的格式，還是要費(fèi)點(diǎn)功夫。筆者曾經(jīng)學(xué)過(guò)latex，發(fā)現(xiàn)latex代碼可以轉(zhuǎn)換為mathtype的代碼格式。latex的矩陣代碼如下：
\begin{bmatrix}
1 & 2 & 3\\
4 & 5 & 6 \\
7 & 8 & 9
\end{bmatrix}
那么，接下來(lái)要做的事情就簡(jiǎn)單了，我們僅需要將數(shù)據(jù)寫(xiě)出來(lái)，數(shù)據(jù)分隔符設(shè)置為“&”，再加上頭尾的標(biāo)識(shí)符、換行符就可以了。
MATLAB的輸入寫(xiě)出可以用下列語(yǔ)句：
dlmwrite(filename,data,'delimiter','&','precision',4)
其中，filename為文件名及其所在地址，僅填寫(xiě)文件名會(huì)保存在工作目錄，如：’C:\user\filename.txt’；data為保存的數(shù)據(jù)；4為保存數(shù)據(jù)的位數(shù)，精度高可以將該數(shù)字增大。保存的數(shù)據(jù)為科學(xué)計(jì)數(shù)法。
展開(kāi) 如何從Ansys APDL中提取剛度矩陣與質(zhì)量矩陣？ ￥69
4.實(shí)戰(zhàn)應(yīng)用與范例講解
接上一個(gè)矩陣的例子，其實(shí)際為Ansys中的一個(gè)應(yīng)力集中問(wèn)題模型所導(dǎo)出的剛度矩陣，那么我們?nèi)绾蝸?lái)驗(yàn)證其結(jié)果的準(zhǔn)確性呢，這時(shí)我們就要用到結(jié)點(diǎn)力矩陣來(lái)進(jìn)行驗(yàn)證了，只要所解出來(lái)的位移與Ansys中可展示的結(jié)點(diǎn)位移相同，那么就證明我們的結(jié)果是準(zhǔn)確無(wú)誤的。
以下我們來(lái)進(jìn)行求解，導(dǎo)入結(jié)點(diǎn)力向量矩陣，使用任意求解器來(lái)進(jìn)行求解：
此時(shí)我們便得到了X向量矩陣，也就是結(jié)點(diǎn)的位移矩陣。由于這里用的是高斯消元法，因此計(jì)算運(yùn)行可能會(huì)有些慢。
然后我們打開(kāi)查看結(jié)果：
與Ansys中List導(dǎo)出的結(jié)點(diǎn)位移結(jié)果一致。結(jié)果正確。
5.代碼購(gòu)買(mǎi)說(shuō)明
本源代碼理論上適用于Ansys APDL中導(dǎo)出的各種hb格式矩陣，無(wú)任何限制，購(gòu)買(mǎi)后如有任何問(wèn)題都可以私信本人進(jìn)行答疑，不僅是此代碼方面，任何有關(guān)有限元軟件學(xué)習(xí)的問(wèn)題也都可以向我請(qǐng)教，我也會(huì)盡我所能去幫助大家。
現(xiàn)在購(gòu)買(mǎi)代碼限時(shí)附贈(zèng)求解器程序源代碼以及結(jié)點(diǎn)力向量導(dǎo)出與使用方法、有限元理論教程，歡迎大家一起學(xué)習(xí)與討論。
展開(kāi) DMIG（Direct Matrix Input at Grid Points）直接矩陣輸入技術(shù)
DMIG（Direct Matrix Input at Grid Points）直接矩陣輸入技術(shù)屬于超單元的范疇。在整車(chē)級(jí)的優(yōu)化和多方案計(jì)算中，借助RADIOSS 的直接矩陣輸入（DMIG）技術(shù)，CAE 工程師可以求解超過(guò)計(jì)算機(jī)資源范圍的大問(wèn)題，同時(shí)能夠縮短求解時(shí)間，提高計(jì)算效率。
這項(xiàng)技術(shù)目前應(yīng)用的還不是很廣泛。在此分享一篇基于Radioss的DMIG仿真分析的論文供大家參考。
RADIOSS直接矩陣輸入技術(shù)在整車(chē)多方案計(jì)算中的應(yīng)用.pdf
基于 MATLAB 的 ANSYS Harwell-Boeing 格式稀疏矩陣提取工具 —— 剛度矩陣與質(zhì)量矩陣 ￥30
在有限元分析中，ANSYS 可以導(dǎo)出大規(guī)模稀疏矩陣（如剛度矩陣、質(zhì)量矩陣），通常使用 Harwell-Boeing (HB) CCS 格式。這些矩陣對(duì)后續(xù)二次開(kāi)發(fā)、動(dòng)力學(xué)分析或自定義求解器非常重要，但由于其稀疏和壓縮存儲(chǔ)形式，直接在 MATLAB 中讀取和使用并不方便。
本文提供了 兩個(gè) MATLAB 函數(shù)，可直接從 ANSYS 導(dǎo)出的 HB 矩陣文件中讀取并重構(gòu)成 MATLAB 稀疏矩陣：
1.剛度矩陣提取函數(shù)
輸入：ANSYS 導(dǎo)出的剛度矩陣 HB 文件（stiff.txt）
輸出：MATLAB 稀疏矩陣 K，可直接用于動(dòng)力學(xué)計(jì)算或驗(yàn)證
支持自動(dòng)對(duì)稱化，保證數(shù)值正確
2.質(zhì)量矩陣提取函數(shù)
輸入：ANSYS 導(dǎo)出的質(zhì)量矩陣 HB 文件（mass.txt）
輸出：MATLAB 稀疏矩陣 M
使用與剛度矩陣同樣的解析邏輯，無(wú)需額外修改
案例說(shuō)明：
本文以高速鐵路接觸網(wǎng)結(jié)構(gòu)為例，展示了如何將 ANSYS 中導(dǎo)出的稀疏剛度矩陣和質(zhì)量矩陣，在 MATLAB 中完整展開(kāi)，并進(jìn)行后續(xù)動(dòng)力學(xué)分析準(zhǔn)備。
通過(guò)該方法，可以將大規(guī)模有限元矩陣快速轉(zhuǎn)化為 MATLAB 可操作形式，為自定義振動(dòng)分析、模態(tài)分析及其他科研或工程應(yīng)用提供基礎(chǔ)。
優(yōu)勢(shì)與應(yīng)用：
支持大規(guī)模稀疏矩陣解析
自動(dòng)對(duì)稱化，保證數(shù)值精度
適用于剛度矩陣、質(zhì)量矩陣、其他 HB 格式矩陣
可作為動(dòng)力學(xué)求解器或后處理工具的基礎(chǔ)模塊
使用方法：
1.使用以下代碼對(duì)ansys中生成的質(zhì)量及剛度矩陣進(jìn)行提取，file,5,full（5為工作目錄下full文件的文件名，例如：filename.full）。
展開(kāi) 
提取ANSYS中的剛度矩陣和質(zhì)量矩陣到MATLAB中進(jìn)行二次開(kāi)發(fā) ￥88
最近在考慮自己編寫(xiě)的程序和商用軟件的驗(yàn)證問(wèn)題，有限元結(jié)構(gòu)分析中最關(guān)鍵的一環(huán)就是剛度矩陣的獲得，如果涉及到模態(tài)分析，還有質(zhì)量矩陣。考慮到商業(yè)軟件的成熟性，可以用ANSYS生成的剛度矩陣做參照來(lái)看自己編寫(xiě)的程序是否正確，因此如何提取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：輸出剛度矩陣。可用于寫(xiě)入了.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) Ansys Zemax | 如何使用瓊斯矩陣表面
在序列模式下，該模型表示為“瓊斯矩陣”表面；在非序列模式下該模型表示為“瓊斯矩陣”物體。“瓊斯矩陣”根據(jù)下式描述瓊斯向量（表示電場(chǎng)）：
其中A, B, C, D均為復(fù)數(shù)。您可以在透鏡數(shù)據(jù)編輯器或非序列元件編輯器中分別輸入這些復(fù)參數(shù)的實(shí)部和虛部。
需要特別注意的是，瓊斯矩陣沒(méi)有定義Ez分量。這意味著使用瓊斯矩陣表面或物體的前提假設(shè)是入射光線需垂直于瓊斯矩陣表面，例如將瓊斯矩陣表面放置在平行光束中。該假設(shè)也與大部分實(shí)際應(yīng)用環(huán)境相符：多數(shù)起偏器或波片都是在平行光或發(fā)散角較小的光束中使用的。
如果一束平行光垂直入射至瓊斯矩陣表面，則由于k·E = 0 并且向量k可表示為{0, 0, 1} 因此Ez必須為零，這樣我們就可以只用Ex和Ey分量來(lái)描述偏振。如果入射光的方向向量為其他任意值 {l, m, n}，則OpticStudio將自動(dòng)調(diào)整Ez或{Ex, Ey}以使k·E = 0且E的大小不會(huì)增加。這個(gè)調(diào)整有可能導(dǎo)致E的大小降低，進(jìn)而導(dǎo)致透過(guò)能量的降低。
下表為一些典型偏振器件的瓊斯矩陣參數(shù)，該表格取自用戶手冊(cè)“The Setup Tab”一章:
使用實(shí)例
接下來(lái)是使用瓊斯矩陣模擬四分之一波片的實(shí)例。請(qǐng)聯(lián)系工作人員獲取附件。
需要注意的是：瓊斯矩陣表面不使用曲率半徑這一參數(shù)，該表面類型總是一個(gè)平面。這是因?yàn)樵擃愋捅砻嫱ǔ６际窃诖怪比肷涞钠叫泄庵惺褂谩?em>矩陣的每個(gè)參數(shù)可以在透鏡數(shù)據(jù)編輯器中的參數(shù)欄中輸入。
展開(kāi) Ansys Zemax | 如何使用瓊斯矩陣表面
打開(kāi)該工具，設(shè)置輸入光的偏振態(tài)為左旋圓偏振光：
經(jīng)過(guò)四分之一波片后可以看到左旋偏振光變?yōu)榫€偏光，且透過(guò)率為100%：
如果我們改變瓊斯矩陣的參數(shù)使其改變?yōu)閄方向的半波片（A real=-1，D real=+1，其余均為0），則可以看到輸出光的偏振態(tài)變?yōu)橛倚窆猓ㄗ⒁馄駲E圓的箭頭標(biāo)注）：
如果我們?cè)O(shè)置瓊斯矩陣為X方向的檢偏器（A real=1，其余為0），則只有X方向的偏振光透過(guò)，且透過(guò)率自然而然地降為50%：
需要注意的是：所有在分析 (Analysis) 選項(xiàng)卡 > 偏振 (Polarization) 標(biāo)簽下的分析功能均支持在當(dāng)前工具中直接定義輸入光的偏振態(tài)。如果您想要在其它分析功能中考慮偏振的影響，例如惠更斯PSF，您需要該分析的設(shè)置中勾選“使用偏振 (Use Polarization)”。在這類分析中您無(wú)法直接定義輸入光的偏振態(tài)，您需要在系統(tǒng)全局設(shè)置中進(jìn)行定義。該設(shè)置位于系統(tǒng)選項(xiàng) (System Explorer) > 偏振 (Polarization) 一欄內(nèi)。
小結(jié)
OpticStudio可以詳細(xì)模擬由雙折射或偏振薄膜所引入的偏振影響。然而有時(shí)候，我們也需要一個(gè)快捷的方法直接輸入偏振數(shù)據(jù)。瓊斯矩陣表面（序列模式下）和瓊斯矩陣物體（非序列模式下）可以幫助我們簡(jiǎn)單快速的定義偏振元件。
由于瓊斯矩陣在定義偏振元件時(shí)只使用Ex和Ey分量，該物體需在垂直入射的平行光中使用。OpticStudio也可以計(jì)算非垂直入射或非準(zhǔn)直的光線，但OpticStudio將自動(dòng)調(diào)整Ez或Ex和Ey以保持向量E垂直于向量k，并且在使用該計(jì)算結(jié)果時(shí)您需格外注意。
展開(kāi) ANSYS模型剛度、質(zhì)量矩陣快速提取小軟件—km_from_Ansys ￥88
模態(tài)分析完成后，在ansys中運(yùn)行HBMAT命令得到十進(jìn)制文件：
! 獲取矩陣十進(jìn)制文件
/aux2
file,Frame,full !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!與文件名一致
hbmat,Stiff,txt,,ascii,stiffness,yes,yes !!!!!!!!!!!!!!!!!!!!!!! 提取剛度矩陣 存儲(chǔ)為Stiff.txt
hbmat,Mass,txt,,ascii,mass,yes,yes !!!!!!!!!!!!!!!!!!!!!!! 提取質(zhì)量矩陣 存儲(chǔ)為Mass.txt
Finish
運(yùn)行完成后，在根目錄下生成了四個(gè)文件，這便是轉(zhuǎn)換完成的十進(jìn)制文件，接著需要做的就是將十進(jìn)制文件轉(zhuǎn)換為矩陣形式，就是通過(guò)工具來(lái)實(shí)現(xiàn)。
(3). 通過(guò)工具實(shí)現(xiàn)提取
A. 打開(kāi)運(yùn)行工具，在km_from_Ansys程序目錄下，在Matlab命令行輸入km_from_Ansys，是初始界面
B.點(diǎn)擊“剛度文件路徑”后的“獲取路徑按鈕”，找到剛剛生成的十進(jìn)制文件的存儲(chǔ)路徑，選擇剛度十進(jìn)制文件Stiff.txt：
同樣的方法選擇好質(zhì)量矩陣的文件路徑，和Mapping文件路徑（Mapping文件有兩個(gè)，選擇任意一個(gè)都可以），以及想要把提取完成后的矩陣文件存儲(chǔ)的路徑，全部選擇完成后，點(diǎn)擊“生成文件”。
C. 接著是大約幾秒的工具處理過(guò)程，根據(jù)模型大小不同，所需時(shí)間不等。處理完成后，根據(jù)提示輸入生成的矩陣文件的存儲(chǔ)文件名，比如“KM”。
，點(diǎn)擊“OK”，完成操作，提示信息為文件的存儲(chǔ)路徑，文件格式為Matlab數(shù)據(jù)文件格式。
D.
展開(kāi) 我自己編的ansys輸入文件轉(zhuǎn)到marc輸入數(shù)據(jù)文件的APDL程序 ***
最好在ansys前處理器環(huán)境下進(jìn)行轉(zhuǎn)換操作,還要注意下面的命令要采用批處理方式輸入,不要復(fù)制粘貼.
/COM, ==========================================================
/COM,
/COM, Beijing University of Technology
/COM,
/COM, Beijing, 100022, P.R. China
/COM,
/COM, WITH HONEYCOMBS MODEL EXAMPLE
/COM, Apr. 2006
/COM, ANSYS 10.0/MARC.2005R2
/COM,
/COM, ==========================================================
/COM,
/COM, MA LIANHUA
/COM, School of Mechanical Engineering,
/COM, Beijing University of Technology
/COM, QQ: 29128203
/COM, Email:mark@emails.bjut.edu.cn
/COM,
/COM, ==========================================================
!因?yàn)镸ARC的輸入文件格式是用行與格來(lái)明確定義它所表示的內(nèi)容,所以它不能有任意的空白行或空白鍵出現(xiàn),否則會(huì)發(fā)生讀取資料有誤的情況.
!
展開(kāi) ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時(shí)用編輯器打開(kāi)cp.out文件，可以看到按單元寫(xiě)出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取.rar
<手搓有限元 干翻Ansys> 【1】C++實(shí)現(xiàn)矩陣Matrix類 實(shí)現(xiàn)基本運(yùn)算
本系列文章致力于實(shí)現(xiàn)“手搓有限元，干翻Ansys的目標(biāo)”，基本框架為前端顯示使用QT實(shí)現(xiàn)交互，后端計(jì)算采用Visual Studio C++。
Matrix類
矩陣基本類，用于有限元矩陣計(jì)算。
1、public function
1.1、構(gòu)造函數(shù)與析構(gòu)函數(shù)
構(gòu)造函數(shù)用來(lái)初始化矩陣，析構(gòu)函數(shù)用來(lái)釋放內(nèi)存。
ANSYS剛度矩陣的提取與解析(python解析）
就ansys如何提取剛度矩陣、如何解讀提取的文檔以及利用Python進(jìn)行解析。
在workbench中實(shí)現(xiàn)整個(gè)過(guò)程的參數(shù)化過(guò)程除了前幾次文章介紹的模型與網(wǎng)格，還應(yīng)該包括材料參數(shù)的參數(shù)化定義。利用Python進(jìn)行二次開(kāi)發(fā)能夠?qū)崿F(xiàn)材料參數(shù)的自由定義，比如來(lái)源于excel表格或者文檔的數(shù)據(jù)，通過(guò)Python代碼的自動(dòng)讀取，參與到實(shí)際的有限元分析進(jìn)程中。
結(jié)構(gòu)有限元最后的求解過(guò)程總是歸結(jié)到求解一個(gè)大型矩陣方程Ax=b，對(duì)于一些情況還需要考慮質(zhì)量矩陣M和阻尼矩陣C。有限元程序在組裝完所有單元的剛度矩陣后，考慮模型所施加的約束和載荷，最終將剛度矩陣進(jìn)行一些處理，例如乘大數(shù)法，變成Ax=b的形式，其中A是剛度矩陣，b是節(jié)點(diǎn)載荷，x為待求的節(jié)點(diǎn)位移，A和b全為已知量。
基本上各類有限元軟件均能夠提取模型的剛度矩陣，此次針對(duì)剛度矩陣的提取與解析做一個(gè)例子，采用的軟件是ANSYS經(jīng)典。
在ANSYS中建立一個(gè)簡(jiǎn)單的模型，劃分網(wǎng)格后共12個(gè)節(jié)點(diǎn)，定義材料參數(shù)，施加約束和載荷后求解。有限元模型如下所示。
待求解結(jié)束后，會(huì)在工作目錄下生成一個(gè)后綴為full的文件，之后即可進(jìn)行剛度矩陣的提取。
通過(guò)主菜單，如下所示。
選擇Matrix后，彈出如下所示的界面。
其中，F(xiàn)ile to be read需要指定工作目錄下生成的full文件，Name of file to write為所導(dǎo)出剛度矩陣的文件名稱；Output matrix file format表示文件格式，還有Binary，生成的是文檔文件，選擇Ascii即可；Matrix to write表示輸出的是剛度矩陣/質(zhì)量矩陣還是阻尼矩陣；RHS選項(xiàng)表示是否同時(shí)輸出右端項(xiàng)，也即是Ax=b中的b。
打開(kāi)生成的剛度矩陣文檔，如下所示。
展開(kāi) ansys之——地震波的輸入和求解
對(duì)于地震波的輸入，可以把荷載記錄做成文件，利用apdl的讀取功能讀入倒數(shù)據(jù)庫(kù)中。下面的例子是自己編的一個(gè)小文件。修改一下可以更簡(jiǎn)潔。有用到的朋友自己作一下把。
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波時(shí)程記錄分成了3個(gè)文件，每個(gè)文件是一列。分別記錄x,y,z方向的加速度。acett是時(shí)間記錄。
這樣就可以把加速度記錄讀取倒ansys數(shù)據(jù)庫(kù)中作為數(shù)組。
也可以把加速度記錄做成一個(gè)文件，這樣程序就簡(jiǎn)單多了。大家可以試看看修改一下。
下面是計(jì)算部分語(yǔ)句：
/SOLU
ANTYPE,trans
！
展開(kāi) ANSYS非線性分析MISO模型數(shù)據(jù)輸入的問(wèn)題
在ANSYS10.0及以前版本中，即便有下降段也可以繼續(xù)計(jì)算，但ANSYS12.0以后版本遇到下降段就無(wú)法計(jì)算了。這是因?yàn)槔习姹拒浖皇前堰@個(gè)錯(cuò)誤忽略掉，實(shí)際上并未解決，新版本軟件則老老實(shí)實(shí)地通知了用戶而已。
如何解決這個(gè)問(wèn)題呢？
用上面的實(shí)例來(lái)說(shuō)，就將最后的*0.85去掉即可，即把曲線的下降段換做水平段。
以上材料定義的案例，來(lái)自王新敏老師著《ANSYS工程結(jié)構(gòu)數(shù)值分析》，因?yàn)橐部吹接腥嗽谡搲锇l(fā)帖說(shuō)書(shū)中命令流材料定義有問(wèn)題，試過(guò)之后確認(rèn)書(shū)中內(nèi)容準(zhǔn)確可用。
———————-補(bǔ)充 —————
可能是上面沒(méi)有圖，不形象，所以有的同學(xué)沒(méi)能完全理解。
所以這里還是針對(duì)上面的命令流，用圖形來(lái)表達(dá)。
設(shè)置好MISO屬性后，可以利用TBPLOT命令把這條曲線繪制出來(lái)：
tbplot,miso,1 ; 繪制材料1的miso曲線
在修改前，即最后一行為“tbpt,,0.0033,fc*0.85”的時(shí)候，繪制出來(lái)的曲線如下：
因?yàn)橛邢陆刀危栽谶M(jìn)行分析的時(shí)候悲催的遇到了下面的錯(cuò)誤提示：
于是，將最后一行的0.85改成1以后，不要下降段，材料曲線變成了這樣的：
調(diào)整后，就可以計(jì)算了。
那么，“第1點(diǎn)的斜率”呢？
看到圖上的編號(hào)了吧，第一點(diǎn)的斜率，就是fc*0.19/0.0002，讓這個(gè)數(shù)等于彈性模量就OK了。
更多優(yōu)質(zhì)精彩內(nèi)容，請(qǐng)關(guān)注公眾號(hào)：十千牛
展開(kāi) 『分享』ANSYS中整體、單元?jiǎng)偠群唾|(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、單元?jiǎng)偠群唾|(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) 