ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖

2017-02-23  by:CAE仿真在線  來(lái)源:互聯(lián)網(wǎng)



編者按

遇到幾次問(wèn)詢關(guān)于ANSYS重新生成云圖的問(wèn)題,今天以一個(gè)例子作個(gè)比較詳細(xì)的說(shuō)明。



正文


在實(shí)際應(yīng)用中經(jīng)常會(huì)需要將ANSYS的計(jì)算結(jié)果進(jìn)行二次運(yùn)算,并重新顯示新的結(jié)果云圖。也即是ANSYS的現(xiàn)有計(jì)算結(jié)果不能滿足實(shí)際需求,需要自己在一次計(jì)算結(jié)果的基礎(chǔ)之上編寫(xiě)計(jì)算方法。

在ANSYS里,這個(gè)過(guò)程是通過(guò)修改節(jié)點(diǎn)或者單元的結(jié)來(lái)實(shí)現(xiàn)的,筆者以前在ABAQUS里面也實(shí)現(xiàn)過(guò)這樣的過(guò)程,不過(guò)在ABAQUS里面不是通過(guò)直接修改節(jié)點(diǎn)/單元解實(shí)現(xiàn)的,而是可以重新定義新的結(jié)果變量。

ANSYS修改節(jié)點(diǎn)解釋通過(guò)DNSOL命令完成的,命令解釋如下:

DNSOL, NODE, Item, Comp, V1, V2, V3, V4, V5, V6

其中Item和Comp這兩個(gè)量是需要修改的變量名稱,例如需要修改位移X,則Item應(yīng)為U,Comp應(yīng)為X,后面的V1-V6就是新的值。

可以看到,利用DNSOL命令每次只能修改一個(gè)節(jié)點(diǎn)的值,因此,很多情況下是需要對(duì)全部的節(jié)點(diǎn)值進(jìn)行修改,故通常需要和遍歷命令一起使用。

下面以一個(gè)實(shí)際的例子說(shuō)明具體的使用方法。

新建一個(gè)簡(jiǎn)單的模型,加載求解得到以下的結(jié)果,分別為x方向的位移和y方向的位移。

X方向的位移如下:

ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖ansys圖片1

Y方向的位移如下:

ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖ansys圖片2

現(xiàn)在需要將X和Y方向的位移進(jìn)行重新計(jì)算,假設(shè):

新的Ux=Ux**2+0.1

新的Uy=Uy**2-0.1

具體實(shí)現(xiàn)過(guò)程是先依次讀取計(jì)算得到的Ux和Uy,保存在自定義的數(shù)組中,然后定義新的數(shù)組,將前面的數(shù)組的數(shù)據(jù)分別處理后保存到新的數(shù)組之中,最后依次讀取新的數(shù)組的數(shù)據(jù)再通過(guò)DNSOL命令進(jìn)行修改,修改完成即可顯示新的結(jié)果。

!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!讀取并保存原始節(jié)點(diǎn)解

alls,all

*get,node_num,node,,count

*dim,dis_data,array,node_num,3

*do,i,1,node_num,1

*get,Ux_temp,node,i,u,x

*get,Uy_temp,node,i,u,y

dis_data(i,1)=i

dis_data(i,2)=Ux_temp

dis_data(i,3)=Uy_temp

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!重新計(jì)算節(jié)點(diǎn)解并保存

*dim,dis_data_new,array,node_num,3

*do,i,1,node_num,1

dis_data_new(i,1)=dis_data(i,1)

dis_data_new(i,2)=dis_data(i,2)**2+0.1

dis_data_new(i,3)=dis_data(i,3)**2-0.1

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!逐一修改節(jié)點(diǎn)解

*do,i,1,node_num,1

dnsol,i,u,x,dis_data_new(i,2)

dnsol,i,u,y,dis_data_new(i,3)

*enddo


這里面需要注意一個(gè)問(wèn)題,利用DNSOL命令是必須在激活FULL graphics的情況下,因此必須關(guān)閉POWRGRPH。

之后重新顯示x和y方向的位移,分別如下:

X方向的位移如下:

ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖ansys分析圖片3

Y方向的位移如下:

ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖ansys分析圖片4

可以看到,云圖形狀發(fā)生了變化,數(shù)值也發(fā)生了變化。

修改之前:

X方向位移:最大0.506,最小0

Y方向位移:最大0.079,最小-0.079

修改之后:

X方向位移:最大0.356,最小0.1

Y方向位移:最大-0.0793,最小-0.1


APDL命令流


/prep7

et,1,182

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,EX,1,,196e3

MPDATA,PRXY,1,,0.3

MPDATA,DENS,1,,7.85E-9

RECTNG,-50,50,-50,50,

wprot,0,0,90

asbw,all

wprot,0,90

asbw,all

alls,all

lsel,r,length,,50

lesize,all,,,10

alls,all

amesh,all

nummrg,node,,,,low

/solu

alls,all

lsel,r,loc,x,-50

aplot

dl,all,,all

alls,all

lsel,r,loc,x,50

sfl,all,pres,-1000

alls,all

solve

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/post1

alls,all

*get,node_num,node,,count

*dim,dis_data,array,node_num,3

*do,i,1,node_num,1

*get,Ux_temp,node,i,u,x

*get,Uy_temp,node,i,u,y

dis_data(i,1)=i

dis_data(i,2)=Ux_temp

dis_data(i,3)=Uy_temp

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*dim,dis_data_new,array,node_num,3

*do,i,1,node_num,1

dis_data_new(i,1)=dis_data(i,1)

dis_data_new(i,2)=dis_data(i,2)**2+0.1

dis_data_new(i,3)=dis_data(i,3)**2-0.1

*enddo

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*do,i,1,node_num,1

dnsol,i,u,x,dis_data_new(i,2)

dnsol,i,u,y,dis_data_new(i,3)

*enddo



開(kāi)放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:ANSYS的結(jié)果進(jìn)行二次運(yùn)算并顯示云圖 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識(shí) ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn) 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過(guò)下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國(guó)服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com