ansys輸出結(jié)果的含義的案例
APDL實(shí)現(xiàn)ANSYS的結(jié)果輸出
APDL實(shí)現(xiàn)ANSYS的結(jié)果輸出
長(zhǎng)安CAE
1 概述
ANSYS作為通用有限元仿真計(jì)算軟件應(yīng)用非常普遍，其強(qiáng)大的參數(shù)化功能更是使得ANSYS的使用方便很多。在計(jì)算過(guò)程中，通常需要對(duì)計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì)，ANSYS計(jì)算的模型大多時(shí)候節(jié)點(diǎn)單元數(shù)目很多，結(jié)果數(shù)據(jù)也很多，因此在GUI界面的查詢操作不太方便，工作量太大，而且結(jié)果不好記錄。而且很多時(shí)候需要借助第三方軟件比如MATLAB等數(shù)據(jù)處理軟件對(duì)結(jié)果進(jìn)行處理，此情況下需要ANSYS能夠?qū)⒂?jì)算結(jié)果輸出到文件中，再通過(guò)對(duì)結(jié)果文件進(jìn)行處理得到想要的結(jié)果。
APDL不僅在參數(shù)化建模、加載和求解有很大的發(fā)揮，在后處理時(shí)同樣有發(fā)揮。本次借助APDL編寫后處理命令流，將需要的計(jì)算結(jié)果以設(shè)置好的格式輸出為txt文件。
2 應(yīng)用
將編寫好的前處理和求解APDL命令流在ANSYS里面運(yùn)行，命令流文件如圖1所示，該段命令流為兩根梁的求解計(jì)算，如圖2所示。
圖1 前處理與求解APDL命令流
圖2 模型
首先將所有節(jié)點(diǎn)的坐標(biāo)位置輸出，命令流如圖3所示。
圖3 輸出節(jié)點(diǎn)位移APDL命令流
ANSYS還可以將結(jié)果保存為數(shù)組，如圖4所示，數(shù)組可以仿照?qǐng)D3的循環(huán)部分輸出，數(shù)組中的數(shù)據(jù)可以繼續(xù)在ANSYS里面反復(fù)調(diào)用。
圖4 輸出數(shù)組
同樣的，其他的結(jié)果參數(shù)也可以輸出。如圖5所示，為節(jié)點(diǎn)應(yīng)力結(jié)果的輸出。
圖5 APDL獲取其他的結(jié)果
基本思路是通過(guò)循環(huán)語(yǔ)句獲取所有節(jié)點(diǎn)的參量，保存成數(shù)組，也可以不保存，然后新建一個(gè)文件用于輸出節(jié)點(diǎn)的參量，主要的命令是*CREATE和*GET。
展開 ANSYS中的自動(dòng)化參數(shù)研究，自動(dòng)建模/分網(wǎng)/多參數(shù)求解/自動(dòng)輸出云圖/自動(dòng)輸出所需結(jié)果
前言
我們經(jīng)常會(huì)進(jìn)行一些具有
可變參數(shù)的有限元模型
的求解，以觀察某些結(jié)果量對(duì)這些參數(shù)的敏感性。在ANSYS中有很多方法可以實(shí)現(xiàn)這一點(diǎn)。當(dāng)然，最簡(jiǎn)單粗暴的就是一個(gè)參數(shù)建一次模型，求解一次。
本文給出的教程案例是通過(guò)使用數(shù)組將參數(shù)的各種值存儲(chǔ)在第一列中，
然后，使用*do命令循環(huán)遍歷數(shù)組中的3個(gè)值，對(duì)于數(shù)組中的每個(gè)值，分別進(jìn)行一次求解。
本教程案例提取模型最大主應(yīng)力存儲(chǔ)在陣列的第二列中，繪制最大主應(yīng)力等值線圖，參數(shù)值作為標(biāo)簽在圖上標(biāo)出。繪圖存儲(chǔ)為jpeg圖片文件，對(duì)研究的參數(shù)的3個(gè)值中的每一個(gè)重復(fù)此操作。
最后，寫入文本文件，并將其與所產(chǎn)生的最大主應(yīng)力一起列出的參數(shù)回顯到屏幕上。
通過(guò)使用具有更多列的數(shù)組，此技術(shù)可以擴(kuò)展到多個(gè)參數(shù)，這項(xiàng)技術(shù)可以自動(dòng)化參數(shù)分析，并產(chǎn)生高生產(chǎn)率收益。
雖然，workbench也可以進(jìn)行這樣子的參數(shù)分析，但還是沒有命令流方便，
這里也顯示了ANSYS APDL命令流建模分析相對(duì)于經(jīng)典界面操作和workbench的一個(gè) 優(yōu)點(diǎn)。
關(guān)注
芷行說(shuō)
公眾號(hào)，后臺(tái)私信獲取完整命令流。
在本教程案例中，我們研究的是如下圖模型，左邊界固支約束，右邊界施加面載荷。
模型建立
通過(guò)以下命令，定義塊體的大小幾何參數(shù)，塊體中孔的位置參數(shù)，同時(shí)定義了3行兩列數(shù)組，其中第一列儲(chǔ)存要研究的3個(gè)孔直徑參數(shù)。
展開 ANSYS如何批量輸出結(jié)果圖
ANSYS提供了很多圖片格式，但有些格式只適用于特定的操作系統(tǒng)或者Device，且有些圖片格式是不能通過(guò)APDL語(yǔ)言導(dǎo)出來(lái)的(ANSYS導(dǎo)出圖片的命令流有三種，水哥只推薦/image，其他兩種使用起來(lái)太繁瑣，不做推薦)，只能采用GUI操作。本人就比較偏愛的三種格式簡(jiǎn)單如下：
1）JPG
此種格式清晰度較高，且可以根據(jù)自己需要設(shè)置圖片質(zhì)量高低，經(jīng)常不涉及批量出圖，只需一兩張結(jié)果圖時(shí)，我便會(huì)使用這種格式。但這種格式不能通過(guò)/image命令導(dǎo)出來(lái)，且只會(huì)保留圖片到你剛開始指定的工作目錄下，文件名字不能更改。
GUi路徑如下：
2）BMP、PNG
這兩種格式均可通過(guò)GUI和命令流輸出，GUi輸出和上述JPG的輸出方法如出一轍，這里不再介紹。既然可以使用命令流輸出，那么就可以采用循環(huán)的方式批量出圖了，后面會(huì)做例子簡(jiǎn)單演示。
但此兩種格式的圖片空間較大，一般在2M左右，但對(duì)于如今的磁盤空間來(lái)講，這點(diǎn)缺點(diǎn)微不足道了。
3）emf
這個(gè)格式不用說(shuō)了，對(duì)寫論文的童鞋來(lái)講在熟悉不過(guò)了。此格式唯一不好的地方是不能采用命令流輸出，只能GUi操作。
回歸正題，如何批量出圖。
使用命令流：/image，該命令流使用格式如下：
label一般選擇save，F(xiàn)name 為文件的名字，當(dāng)采用循環(huán)存儲(chǔ)的時(shí)候?yàn)榱吮苊鈭D片相互之間因?yàn)橥桓采w，此處需要用到將數(shù)字轉(zhuǎn)為字符的命令%_%，可通過(guò)引號(hào)指定圖片保存的位置，ext為圖片格式，此處可以為bmp，png。
例：
某等厚度圓環(huán)孔板的模態(tài)分析，已經(jīng)求得其前十階的振動(dòng)模態(tài)，現(xiàn)需要批量導(dǎo)出其前十階模態(tài)。
展開 MatlabGUI界面調(diào)用Ansys計(jì)算并輸出計(jì)算結(jié)果
.*'},'File Selector'); strh = [Pnameh,Fnameh];
pathname = Pnameh;
set(handles.text1,'String',strh);
[temp1,temp2] = xlsread(strh);
set(handles.uitable1,'Data',temp1);
% Update handles structure
guidata(hObject, handles);
為了讀取圖示方框中的數(shù)據(jù)，并用到ANSYS的APDL文件中，需要字符串的讀取和合并，首先需要使用str2num函數(shù)把字符串轉(zhuǎn)換成數(shù)值，如果沒有輸入值時(shí)，使用缺省值。
將兩個(gè)txt合并成test3.mac作為APDL語(yǔ)言開始的參數(shù)定義，生成test3.mac之后再使用system函數(shù)調(diào)用ANSYS的求解器，并讀取test3.mac進(jìn)行計(jì)算
在計(jì)算之前，是不能生成圖片的，這時(shí)需要設(shè)置只有點(diǎn)擊“開始重構(gòu)”按鈕之后，其他按鈕才可用。
點(diǎn)擊按鈕開始計(jì)算之后，會(huì)分別輸出兩個(gè)名為residualstress.jpg和deformation.jpg的圖片，對(duì)應(yīng)的語(yǔ)句為
/image,save,'E:\GUIRStest\residualstress',jpg
設(shè)置當(dāng)點(diǎn)擊“生成殘余應(yīng)力云圖”和“生成角變形云圖”時(shí)，會(huì)讀取圖片的路徑并使用imshow生成圖片。
至此，一個(gè)簡(jiǎn)易的MatlabGUI界面調(diào)用ANSYS計(jì)算并輸出圖片就完成了。
展開 
Moldex3D模流分析之結(jié)果輸出至Ansys Workbench
Moldex3D 提供射出成型結(jié)果中纖維配向、初始應(yīng)力 (翹曲應(yīng)力)、纖維濃度以及縫合線的輸出。從 Moldex3D 輸出的檔案可直接由 Ansys Workbench 讀取，并可與 Ansys Material Designer 提供的材料模型進(jìn)行整合，以利于纖維強(qiáng)化復(fù)材件的射出模擬。以下是使用Moldex3D 輸出結(jié)果項(xiàng)至 Ansys Workbench 的操作流程介紹：
?使用纖維強(qiáng)化復(fù)材時(shí)，須注意在前處理的計(jì)算參數(shù)中，有啟用「執(zhí)行纖維配向計(jì)算」功能。若使用者想預(yù)測(cè)纖維濃度 (Fiber Volume Fraction) 在塑件內(nèi)的分布，則須額外勾選進(jìn)階選項(xiàng)中「計(jì)算填料濃度」功能。
?請(qǐng)注意：初始應(yīng)力(*.ist )對(duì)應(yīng)到射出成型分析中，由充填與保壓階段形成在模穴內(nèi)部的應(yīng)力，頂出后所造成的翹曲變形。因此若要輸出初始應(yīng)力，在設(shè)定分析順序時(shí)，應(yīng)選擇有包含充填( F )、保壓 ( P ) 以及翹曲 ( W ) 的分析。
?完成 Moldex3D 模擬分析。
?當(dāng)分析完成后，在 FEA 接口的功能選項(xiàng)中，選取欲輸出的檔案：纖維濃度 (fiber volume fraction) (*.fcd)、初始應(yīng)力 (*.ist)、縫合線 (*.nwd) 或纖維配向 (*.o2d)。除了勾選的功能外，軟件同時(shí)也會(huì)輸出 Ansys 可讀的網(wǎng)格檔(*.ans 與 *.cdb)。
?此處輸出之檔案可直接以 Ansys Workbench Toolbox 中的 Injection Molding Data 讀取。
?在 Ansys Workbench 環(huán)境中，提供一套短纖維復(fù)合材料仿真流程，可預(yù)測(cè)短纖維強(qiáng)化復(fù)材之射出成型件的熱機(jī)械行為。
展開 74基于matlab的PSO-ELM的多輸入，單輸出結(jié)果預(yù)測(cè)，輸出訓(xùn)練集和測(cè)試機(jī)預(yù)測(cè)結(jié)果及誤差 ￥35.9
基于matlab的PSO-ELM的多輸入，單輸出結(jié)果預(yù)測(cè)，輸出訓(xùn)練集和測(cè)試機(jī)預(yù)測(cè)結(jié)果及誤差，適應(yīng)度值。數(shù)據(jù)可更換自己的，程序已調(diào)通，可直接運(yùn)行。
70基于matlab的BP神經(jīng)網(wǎng)絡(luò)多輸入單輸出數(shù)據(jù)結(jié)果預(yù)測(cè)，輸出結(jié)果包括均方根誤差，決定系數(shù)。 ￥25.9
基于matlab的BP神經(jīng)網(wǎng)絡(luò)多輸入單輸出數(shù)據(jù)結(jié)果預(yù)測(cè)，輸出結(jié)果包括均方根誤差，決定系數(shù)。數(shù)據(jù)可更換自己的，程序已調(diào)通，可直接運(yùn)行。
原創(chuàng)：第一方向殘余應(yīng)力結(jié)果中名義厚度0.969含義。
在Moldflow軟件中，1代表產(chǎn)品表面，0是產(chǎn)品中心層，Normalized thickness =0.969，就是非常靠近表皮的斷面位置，只差0.31%的厚度位置。
不了解的話，看下面的圖片說(shuō)明。Additional information, Why 0.969, not 1 or 0.938，0.969 is middle between 1 and 0.938.(也有人問過(guò)這個(gè)問題)
單元應(yīng)力求解方法簡(jiǎn)述、結(jié)果應(yīng)力種類選項(xiàng)值的含義以及計(jì)算舉例 ￥1
后處理求解結(jié)果應(yīng)力種類選項(xiàng)值的含義與計(jì)算舉例：
unveraged
areraged
nodal difference
nodal fraction
elmemntal difference
elemental mean
elemental fraction
【ANSYS文件輸出】ANSYS怎么批量文件輸出？
這期是ANSYS不定期更新，這期主要講一下ANSYS怎么批量文件輸出，還有批量輸出前的一些預(yù)處理操作。</span></p><p><strong>01</strong></p><p><strong>怎么去除右上角的時(shí)間?</strong></p><p>/PLOPTS,DATE,0 !都不顯示</p><p>/PLOPTS,DATE,1 !只顯示日期</p><p>/PLOPTS,DATE,2 !顯示日期和時(shí)間</p><p><strong>02</strong></p><p><strong>怎么去除右上角的LOGO？</strong></p><p>/plopts,logo,on !以logo形式</p><p>/plopts,logo,off !以文本形式</p><p> 用上面兩種形式，logo還是會(huì)存在的，還有一種方法就是在界面，鼠標(biāo)箭頭放在logo上，右鍵點(diǎn)擊data，REPLOT一下，這樣logo就會(huì)不再出現(xiàn)。</p><p><strong>03</strong></p><p><strong>怎么把背景變成白色？</strong></p><p>!
展開 abaqus二次開發(fā)：后處理批量提取場(chǎng)輸出和歷程輸出結(jié)果（源碼帶注釋） ￥198
abaqus二次開發(fā)：后處理批量提取場(chǎng)輸出和歷程輸出結(jié)果（購(gòu)買后需要加V私聊獲取使用教程）
購(gòu)買后，請(qǐng)加V：wzd_1021_
更多插件內(nèi)容，請(qǐng)私聊
同時(shí)提供定制hypermesh/hyperview/abaqus前后處理插件開發(fā)。
ANSYS WORKBENCH中弱彈簧的含義
許多CAE朋友在ANSYS WOKKBENCH中進(jìn)行靜力學(xué)和瞬態(tài)動(dòng)力學(xué)仿真時(shí)，都遇到過(guò)弱彈簧（weak spring）的問題，我們發(fā)現(xiàn)，在求解結(jié)束以后，ANSYS經(jīng)常提到它給我們加了一個(gè)弱彈簧，并建議我們檢查一下模型，這是什么意思呢？弱彈簧是好還是不好，對(duì)于結(jié)果有沒有影響，該不該加，如何加呢？ANSYS加弱彈簧的目的又是什么呢？
我們先考察一個(gè)超級(jí)簡(jiǎn)單的例子，然后通過(guò)該例子來(lái)考察ANSYS所施加的弱彈簧的含義。一根矩形截面梁，長(zhǎng)度為1米，橫截面是100mm*100mm，左邊固定，右邊加10kN的力，現(xiàn)在考慮加力后它的變形和應(yīng)力。
顯然，這是一個(gè)簡(jiǎn)單的拉伸問題，在軸線方向上應(yīng)力都是一樣的，按照拉伸的應(yīng)力公式，可以計(jì)算其大小為1Mpa。我們先對(duì)該問題進(jìn)行建模，然后修改邊界條件，來(lái)考察弱彈簧的含義。
1. 創(chuàng)建靜力學(xué)分析系統(tǒng)。
2. 創(chuàng)建梁的三維模型。
雙擊geometry單元格，進(jìn)入到DM中，然后創(chuàng)建一根三維梁
其尺寸設(shè)置如下
即長(zhǎng)度為1000mm，而截面尺寸是100mm*100mm。
3. 劃分網(wǎng)格得到有限元模型。
雙擊model,進(jìn)入到mechanical中，并自動(dòng)劃分網(wǎng)格，結(jié)果如下。
4. 施加邊界條件。
左端面施加位移邊界條件，三個(gè)方向的位移都為零。
在右端面上施加10KN的拉力。
5.求解并得到結(jié)果。
計(jì)算完畢后，沒有任何警告或者錯(cuò)誤信息，而X方向的位移結(jié)果是
即拉伸了0.00498mm左右。
其應(yīng)力大小是
由于在左邊存在應(yīng)力集中，此處有輕微的變化。而桿件的絕大部分應(yīng)力是1Mpa，這與實(shí)際情況是吻合的。
6.改變位移邊界條件，改變成力的邊界條件。
在上圖中，軸線方向是X方向。該軸也只是在X方向上受力。從理論上看來(lái)，對(duì)于左端面，可以只約束X方向，而Y方向和Z方向可以是自由的。
展開 ABAQUS 結(jié)果文件輸出
結(jié)果文件有output database（.odb）文件和result file（.fil）文件。區(qū)別是
a.odb文件默認(rèn)生成，是二進(jìn)制文件，abaqus才能訪問，python后處理也是通過(guò)abaqus訪問。
b.fil文件不是默認(rèn)生成，需要在inp中添加關(guān)鍵字（*FILE FORMAT, ASCII），添加的位置可以是（*Restart）命令之后。是文本文件，文本編輯器可以訪問。
要輸出的量可以通過(guò)cae界面選擇，也可以通過(guò)inp添加。
輸出到fil文件屬于歷史輸出范圍，輸出的是所有增量步的值。fil文件的大小似乎是固定的，在cae界面中，歷史輸出可以設(shè)置只輸出最后一個(gè)increment的值，對(duì)應(yīng)在inp文件中表示為（*Output, history, frequency=99999），輸出頻率非常大，故只輸出最后一幀。這個(gè)命令針對(duì)fil的輸出無(wú)用，這里姑且認(rèn)為fil只能輸出的所有增量步的值。
*NODE FILE 輸出節(jié)點(diǎn)變量，位移（UCOORD-101），節(jié)點(diǎn)坐標(biāo)（COORD-107）
*EL FILE 輸出單元變量，MISES應(yīng)力（S），應(yīng)力分量（SINV-12）
要使用URDFIL子程序，首先必須先創(chuàng)建.fil文件，子程序讀取.fil文件，F(xiàn)ortran代碼實(shí)現(xiàn)結(jié)果的格式化輸出。
展開 6混淆矩陣，根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣。將自己標(biāo)簽數(shù)據(jù)帶進(jìn)去即可得到結(jié)果，可更改顏色，更改預(yù)測(cè)的個(gè)數(shù) ￥8.9
混淆矩陣，根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣。將自己標(biāo)簽數(shù)據(jù)帶進(jìn)去即可得到結(jié)果，可更改顏色，更改預(yù)測(cè)的個(gè)數(shù)，基于MATLAB平臺(tái)，程序已調(diào)通可直接運(yùn)行。
根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣 ￥9.9
混淆矩陣，根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣。將自己標(biāo)簽數(shù)據(jù)帶進(jìn)去即可得到結(jié)果，可更改顏色，更改預(yù)測(cè)的個(gè)數(shù)，基于MATLAB平臺(tái)，程序已調(diào)通可直接運(yùn)行。標(biāo)價(jià)為程序價(jià)格，不包含售后。程序保證可直接運(yùn)行。
