
發(fā)布
注冊(cè)
/
登錄ansys提取變量的案例
Python提取場(chǎng)變量-干貨 ¥1.9
提取各節(jié)點(diǎn)和單元的場(chǎng)變量,利用函數(shù)查看可提取的場(chǎng)變量有哪些。
在Abaqus的二次開(kāi)發(fā)過(guò)程中,通常需要采用Python腳本語(yǔ)言將Abaqus的計(jì)算結(jié)果進(jìn)行輸出,然后再進(jìn)行處理。Python使Abaqus的內(nèi)核語(yǔ)言,使用較為方便,Abaqus運(yùn)行Python語(yǔ)言的方式有多種,可以直接命令窗口,也可以讀入腳本,還可以采用類(lèi)似批處理的方式。
本次以一個(gè)例子細(xì)說(shuō)Python語(yǔ)言在Abaqus后處理中的應(yīng)用,模型的計(jì)算結(jié)果云圖如圖1所示。
2 輸出所有節(jié)點(diǎn)的Mises應(yīng)力
直接上Python代碼:
import os
myodb=openOdb(path='Job-1.odb')
cpFile=open('artlcF1.txt','w')
RF=myodb.steps['Step-1'].frames[1].fieldOutputs['S'].values
for i in range(len(RF)) :
cpFile.write('%10.3Fn' % (RF[i].mises))
展開(kāi) 三步快速提取路徑單元變量 ¥2
<span style="font-family:'宋體';font-size:12.0pt;white-space:pre-wrap;ABAQUSref="https://www.yqgqt.org.cn/major/abaqus" style="text-decoration: none; cursor: text; color: #333;">ABAQUS后處理
三步快速提取路徑單元變量
例如提取沖擊載荷下,路徑上不同時(shí)刻的位移應(yīng)力等力學(xué)參數(shù)
截面路徑快速建立的方法
Visualization 模塊
1.創(chuàng)立Path路徑
Tool > Path > Create
展開(kāi) Abaqus_Python 后處理結(jié)點(diǎn)數(shù)據(jù)提取、場(chǎng)變量云圖輸出 ¥1.99
每個(gè)case需要輸出幾個(gè)路徑上的部分節(jié)點(diǎn)值,而且需要針對(duì)不同的case、不同的幀、不同的視角截圖、保存云圖,如果人工出圖和提取節(jié)點(diǎn)數(shù)據(jù),工作量會(huì)非常大,費(fèi)時(shí)會(huì)遠(yuǎn)遠(yuǎn)超出計(jì)算所用時(shí)間;并且,同事還要求我根據(jù)case名稱、組件名稱(模型有多個(gè)組件組成,且網(wǎng)格是在part上畫(huà)的,所以同一個(gè)節(jié)點(diǎn)編號(hào)可能對(duì)應(yīng)多個(gè)節(jié)點(diǎn))等保存成不同的文件。
結(jié)點(diǎn)溫度數(shù)據(jù)沒(méi)有保存,因此沒(méi)有輸出。下圖是自動(dòng)輸出的png格式云圖。
比較倉(cāng)促,代碼也比價(jià)粗糙。回頭做一個(gè)詳細(xì)的介紹。
主要代碼見(jiàn)下。
ANSYS Beam188提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù) (解決彎矩圖鋸齒狀) ¥20
在ANSYS中有些數(shù)據(jù)無(wú)法直接訪問(wèn),需要通過(guò)定義單元表完成單元的結(jié)果的訪問(wèn)。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù)的詳細(xì)過(guò)程。
1. 首先需要知道在哪里定義單元表:Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數(shù)據(jù),這里以Beam188的彎矩為例
2.1 啟動(dòng)ANSYS幫助菜單, 在索引框輸入Beam188然后搜索, 在單元輸出介紹找到彎矩的名稱(代號(hào))。
2.2 回到ANSYS界面,比如要輸出Mz, 則需要在添加SMISC,3 和SMISC,16 ,如圖
3. 輸出數(shù)據(jù):Main Menu>General Postproc>Element Table> List E T, 選擇前面定義的SMISC,3 和SMISC,16 輸出單元I和J節(jié)點(diǎn)的Mz數(shù)值,如圖
4. 顯示彎矩云圖:Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res, 這里要注意要在LabI 選SMISC,3 LabJ 選SMSCI,16。
輸出彎矩到這就結(jié)束了,小編突然發(fā)現(xiàn),輸出的彎矩值在每個(gè)單元的I和J處是一樣的(Beam188為2節(jié)點(diǎn)單元),彎矩圖也就成了鋸齒形,于是去問(wèn)了度娘一波,各路盆友給出解決方法,然而并沒(méi)有起作用的,于是乎我又想起來(lái)了“幫助文檔大法”,于是認(rèn)認(rèn)真真將Beam188的幫助文檔閱讀了一遍,功夫不負(fù)有心人,最終。。。
展開(kāi) 
ANSYS的get命令常用操作(信息提取和結(jié)果結(jié)果提取)
ANSYS的get命令常用操作(信息提取和結(jié)果結(jié)果提取)
在ANSYS分析過(guò)程中,*get命令作為一個(gè)提取信息的常用命令,作用非常大,不管是在前處理、求解還是后處理過(guò)程中,都能夠有發(fā)揮的空間,尤其是后處理過(guò)程,對(duì)結(jié)果的批量輸出來(lái)說(shuō)不可缺少。
*get能夠提取的信息相當(dāng)多,其命令語(yǔ)句如下:
*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM
Par:定義的變量名稱,用于存儲(chǔ)提取的數(shù)據(jù);
Entity:關(guān)鍵字,是信息提取的對(duì)象,包括NODE, ELEM, KP, LINE, AREA, VOLU, PDS等;
ENTNUM:當(dāng)前對(duì)象的數(shù)字標(biāo)識(shí),比如節(jié)點(diǎn)的節(jié)點(diǎn)號(hào),單元的單元號(hào)等;
Item1:提取的信息,可用的非常多,后面展開(kāi);
IT1NUM:和Item1配合使用。
由于*get的功用實(shí)在太多,不就一一列舉,單就常用的枚舉。
展開(kāi) ANSYS中的變量
ANSYS中的變量總的來(lái)說(shuō)分為兩大類(lèi):一類(lèi)是標(biāo)量參數(shù);一類(lèi)是數(shù)組參數(shù)。其意義與C語(yǔ)言中的參數(shù)和數(shù)組類(lèi)似。標(biāo)量參數(shù)是指單個(gè)的變量,而數(shù)組則是由一系列具有相同意義的數(shù)據(jù)組成。
一:標(biāo)量參數(shù)
標(biāo)量參數(shù)分為兩種:數(shù)值型的和字符型的。
1:標(biāo)量參數(shù)的定義或修改:可以使用*SET命令或直接寫(xiě)出賦值表達(dá)式。修改時(shí),只需改變數(shù)據(jù)即可。
例如:
*SET,LENGTH,10 !定義了一個(gè)名為L(zhǎng)ENGTH的數(shù)值型變量,并且值為10
*SET,COMMENT,'DISP' !定義了一個(gè)名為COMMENT的字符型變量,并且賦值為DISP
以上兩句等價(jià)于:
LENGTH=10
COMMENT='DISP'
2:標(biāo)量參數(shù)的刪除:同樣可以使用*SET命令或直接寫(xiě)出賦值表達(dá)式或使用*DEL命令。將需要?jiǎng)h除的變量賦空值,注意不是數(shù)字0或空格。
例如:我們將上面定義的兩個(gè)參數(shù)刪除。語(yǔ)句如下:
*SET,LENGTH,
*SET,COMMENT,
或:
LENGTH=
COMMENT=
又或:
*DEL,LENGTH
*DEL,COMMENT
3:標(biāo)量參數(shù)的GUI操作菜單。GUI路徑為:Utility Menu->Parameters->Scalar Parameters。標(biāo)量參數(shù)的定義,刪除和修改都可以在這個(gè)窗口中完成。
二:數(shù)組參數(shù)
數(shù)組參數(shù)按維數(shù)可分為:一維數(shù)組,二維數(shù)組和三維數(shù)組。
按存儲(chǔ)的數(shù)據(jù)類(lèi)型可分為:
1)一般數(shù)組參數(shù),也成為ARRAY Parameter。是一種默認(rèn)的數(shù)據(jù)類(lèi)型,也就是說(shuō),如果用戶不聲明數(shù)據(jù)類(lèi)型,系統(tǒng)就會(huì)自動(dòng)使用該數(shù)據(jù)類(lèi)型。該數(shù)組的行、列、面的索引均為從1開(kāi)始的連續(xù)整數(shù),元素值可以為整數(shù)或?qū)崝?shù)。
2)字符數(shù)組參數(shù),又名CHAR Array Parameter。元素值為不多于8個(gè)的字符或數(shù)字組成。其行,列,面的索引值也是從1開(kāi)始的連續(xù)整數(shù)。
展開(kāi) ANSYS APDL如何求變量的最大值
在其他語(yǔ)言中求最大值非常容易,比如有三個(gè)變量分別是MXS6,MXS7和MXS8,要求他們的最大值賦予MaxS,用到的函數(shù)往往只是一個(gè)函數(shù)MaxS=max(MXS6,MXS7,MXS8)。但是在ANSYS Parametric Design Language Guide第三章第8節(jié) Parametric Functions部分卻找不到這樣的函數(shù)。
無(wú)奈自力更生吧。
MaxS=MXS6
*if,MXS7,GT,MaxS,THEN
MaxS=MXS7
*ENDIF
*if,MXS8,GT,MaxS,THEN
MaxS=MXS8
*ENDIF
這是三個(gè)變量,如果有更多變量方法也是一樣的,變量太多就用個(gè)循環(huán),如果需要留言給我,我給大家把代碼寫(xiě)出來(lái)。
展開(kāi) 如何從Ansys APDL中提取剛度矩陣與質(zhì)量矩陣? ¥69
其可以直接導(dǎo)出full文件中的矩陣數(shù)據(jù),而通過(guò)我下面給大家提供的源代碼便可以直接將導(dǎo)出的質(zhì)量矩陣與剛度矩陣直接轉(zhuǎn)化為matlab中可使用變量,實(shí)現(xiàn)我們的矩陣提取操作。
Hbmat法:
! 提取剛度矩陣
/AUX2
FILE,'file',full ! 將’file’改為自己路徑下的.full文件名
HBMAT, 'Stiffness_mat', dat, , ASCII, STIFF, YES, YES !剛度矩陣
HBMAT, 'Mass_mat', dat, ,ASCII, MASS, YES, YES ! 質(zhì)量矩陣
FINISH
成功導(dǎo)出后你會(huì)在你的工作路徑中看到儲(chǔ)存在Ansys中的剛度矩陣與質(zhì)量矩陣。
2.2GUI方法
目前大部分此類(lèi)教程都忽略了最基本的GUI方法,可能是自動(dòng)帶入了經(jīng)驗(yàn)豐富的工程師角色,但對(duì)于初學(xué)者而言,GUI方法十分有助于理解與學(xué)習(xí)有限元軟件,因此我也將介紹提取剛度質(zhì)量矩陣的GUI方法。
從而我們便能導(dǎo)出所需要的剛度矩陣與質(zhì)量矩陣源文件了。
3.源代碼的使用與轉(zhuǎn)換(使用方法)
我們打開(kāi)所導(dǎo)出剛度矩陣會(huì)發(fā)現(xiàn)里面的數(shù)字比較混亂,一時(shí)間分不出哪些是我們需要的數(shù)據(jù),也沒(méi)有辦法直接拿它用來(lái)計(jì)算,這時(shí),我們便需要使用我們自己的矩陣轉(zhuǎn)換m文件,其可以自動(dòng)幫我們提取出剛度矩陣與質(zhì)量矩陣,并生成matlab中的矩陣變量形式,我們可以直接運(yùn)用這段源代碼來(lái)進(jìn)行操作。
展開(kāi) 【ANSYS經(jīng)驗(yàn)技巧】- 學(xué)會(huì)使用變量數(shù)組函數(shù)(轉(zhuǎn)載)
5.直到你為函數(shù)中所有的狀態(tài)的所有變量提供賦值,才能保存為表格式參數(shù)使用。
注意:在ansys分析中用函數(shù)加載必須的兩個(gè)步驟:
1利用函數(shù)編輯器創(chuàng)建任意方程或函數(shù)
2在利用函數(shù)加載器加載函數(shù),并以函數(shù)定義表參數(shù)
轉(zhuǎn)自:三維網(wǎng)。 作者:4kpolo
原帖鏈接:http://www.3dportal.cn/discuz/viewthread.php?tid=746468&extra=page%3D1%26amp%3Bfilter%3Ddigest
有相互依存關(guān)系的離散變量的ansys與workbench聯(lián)合優(yōu)化分析
需要說(shuō)明的是對(duì)應(yīng)基于離散變量的優(yōu)化,采用不同的響應(yīng)面構(gòu)建方法和優(yōu)化算法,效率相差特別大。即使對(duì)于本問(wèn)題節(jié)點(diǎn)數(shù)目5個(gè),單元數(shù)目6個(gè)。選擇的響應(yīng)面構(gòu)建方法和優(yōu)化算法不同,也有可能計(jì)算幾個(gè)小時(shí)。對(duì)于本問(wèn)題采用Latinhypercube sampling (LHS拉丁超立方體抽樣)生成試驗(yàn)設(shè)計(jì),采用神經(jīng)網(wǎng)絡(luò)方法來(lái)構(gòu)建響應(yīng)面,實(shí)際證明效率較高。
另外對(duì)應(yīng)基于離散變量的優(yōu)化分析,目前workbench只支持篩選法和混合整數(shù)序列二次規(guī)劃優(yōu)化算法。
另外,其實(shí)該問(wèn)題也可以完全采用ansys經(jīng)典完成程序優(yōu)化設(shè)計(jì),利用離散編碼陷阱實(shí)現(xiàn)從連續(xù)變量到離散變量的轉(zhuǎn)變。但是該方法也有很多缺點(diǎn):
1.最終得優(yōu)化的變量依然是連續(xù)的,需要人為后處理,實(shí)現(xiàn)規(guī)格表的編碼。
2.最終得到的優(yōu)化結(jié)果,可能陷入局部最小陷阱。采用首次得到的優(yōu)化結(jié)果為初始值,然后縮小優(yōu)化變量的采用空間,可以一定程度上改善結(jié)果的精度。
3.規(guī)格表的離散區(qū)間步長(zhǎng)對(duì)于求解的效率的影響非常大。因此,需要增大優(yōu)化迭代次數(shù)。
4.系統(tǒng)優(yōu)化過(guò)程中,可能多次在等效解處徘徊。影響求解效率。
5.人為將連續(xù)變量離散化后,基于偏導(dǎo)算法的一階優(yōu)化方法將不能處理該類(lèi)問(wèn)題。
6.最終解碼得到的材料規(guī)格往往需要返回到分析中去,才可以得到真實(shí)的狀態(tài)變量數(shù)值。
完全采用ansys優(yōu)化的具體方法這里不在提供。
這里順便說(shuō)下ansys和workbench優(yōu)化分析的優(yōu)缺點(diǎn):
1.采用ansys可以很方面的實(shí)現(xiàn)網(wǎng)絡(luò)結(jié)構(gòu)的編程和變量提取后控制。對(duì)于類(lèi)似問(wèn)題,如果分析的模型更大,在workbench中建模可以說(shuō)是一件極其痛苦的事情。
2.workbench提供了比ansys更多的優(yōu)化算法。自身就擁有離散變量的優(yōu)化功能。這也或許是現(xiàn)在ansys舍棄經(jīng)典優(yōu)化界面的一個(gè)很大原因。
展開(kāi) ansys優(yōu)化,因變量和目標(biāo)函數(shù)都沒(méi)有變化【急】【急】
ansys優(yōu)化,因變量和目標(biāo)函數(shù)都沒(méi)有變化【急】【急】
ansys優(yōu)化之后,為什么只有自變量發(fā)生了變化,而因變量和目標(biāo)函數(shù)都沒(méi)有變化,還是和初始值一樣?也進(jìn)行了四五十次的迭代,也有顯示最優(yōu)解,只是因變量和目標(biāo)函數(shù)都沒(méi)有變化,疑惑中。

ANSYS模型剛度、質(zhì)量矩陣快速提取小軟件—km_from_Ansys ¥88
通過(guò)求取模態(tài)特征值來(lái)驗(yàn)證提取矩陣的有效性:
>>load('KM.mat')
>> [v,d]=eigs(kk,mm,10,'sm');
>> d=sqrt(d);
>> d=d/2/pi;
>> d=sum(d);
>> d=d';
上圖為Ansys中模態(tài)分析的結(jié)果,下圖為提取矩陣的分析結(jié)果,結(jié)果一致,表明提取矩陣有效。
6. 說(shuō)明
提取工具適用于任意單元,任意結(jié)構(gòu)形式,但僅限于彈性模型。
展開(kāi) ANSYS知識(shí)普及1——如何提取模態(tài)質(zhì)量(ANSYS專(zhuān)家編輯,非原創(chuàng),歡迎轉(zhuǎn)摘)
本人準(zhǔn)備出一個(gè)ANSYS知識(shí)普及系列,將有用的網(wǎng)上資料歸攏,由于知識(shí)水平有限,不對(duì)之處請(qǐng)諒解。也歡迎各位網(wǎng)友提供好的資料分享,讓我們共同完成這個(gè)ANSYS知識(shí)普及系列。
編輯人:技術(shù)鄰ANSYS專(zhuān)家
業(yè)務(wù)咨詢網(wǎng)址:http://www.yqgqt.org.cn/content/other/402981
(打個(gè)小廣告)
聲 明:1、ANSYS知識(shí)普及系列中所有資料均來(lái)自網(wǎng)上;
2、如侵犯知識(shí)產(chǎn)權(quán),請(qǐng)聯(lián)系ANSYS專(zhuān)家本人或者技術(shù)鄰,我將第一時(shí)間刪除。
下面是《ANSYS Verification Manual》中VM89.DAT稍加修改后提取模態(tài)質(zhì)量的例子:
/PREP7
/TITLE, VM89, NATURAL FREQUENCIES OF A TWO-MASS-SPRING SYSTEM
C*** VIBRATION THEORY AND APPLICATIONS, THOMSON, 2ND PRINTING, PAGE 163,EX 6.2-2
ET,1,COMBIN14,,,2
ET,2,MASS21,,,4
R,1,200 ! SPRING CONSTANT = 200
R,2,800 ! SPRING CONSTANT = 800
R,3,.5 ! MASS = .5
R,4,1 ! MASS = 1
N,1
N,4,1
FILL
E,1,2 ! SPRING ELEMENT (TYPE,1) AND K = 200 (REAL,1)
TYPE,2
REAL,3
E,2 ! MASS ELEMENT (TYPE,2) AND MASS = .5 (REAL,3)
TYPE,1
REAL,2
E,2,3 !
展開(kāi) ANSYS workbench中質(zhì)心的提取方法 ¥88
在ansys中,我們經(jīng)常需要提取模型的質(zhì)心,查看其質(zhì)心的位置,以及動(dòng)力學(xué)中查看質(zhì)心的移動(dòng)軌跡,并非某一參考點(diǎn)的移動(dòng)軌跡。在ansys軟件中很容易查看模型,或者某一組件的質(zhì)心如圖所示。只需要在幾何模型中選中模型,然后再property中查看坐標(biāo)即可。但是這種為靜態(tài)結(jié)果,也就是初始狀態(tài)的結(jié)果
如果要查看模型在受力作用下發(fā)生了變形,那么對(duì)應(yīng)的質(zhì)心顯然會(huì)發(fā)生移動(dòng)。但是在這種情況下怎么查看質(zhì)心呢。一個(gè)近似的方法可以如下方式:
在后處理中插入自定義的結(jié)果loc-defy,可以得到模型變形后的坐標(biāo),然后查看average的結(jié)果,就是質(zhì)心的結(jié)果,該方法根據(jù)變形量的平均效果來(lái)得到質(zhì)心的,如下所示
但是真正的質(zhì)心提取方法并非上面的方式,需要準(zhǔn)確獲得,需要采用下面的命令
upcoord,1更新模型單元位置
展開(kāi) ANSYS Workbench中如何提取截面內(nèi)力 ¥3.9
在土木及水利設(shè)計(jì)中,截面內(nèi)力是結(jié)構(gòu)設(shè)計(jì)過(guò)程中極為重要的參數(shù),也是結(jié)構(gòu)穩(wěn)定性的重要依據(jù)。本文重點(diǎn)介紹如何在Workbench平臺(tái)自定義截面并獲得相應(yīng)截面的內(nèi)力,并將其結(jié)果輸出。方法簡(jiǎn)單,操作易上手!最終結(jié)果顯示如下:
具體步驟為:1、自定義創(chuàng)建截面,這里建議采用局部坐標(biāo)系的方法建立截面位置;