matlab/ansys协同的一个例子
前段時(shí)間找了點(diǎn)Matlab與ANSYS協(xié)同工作的資料，在我看來所謂協(xié)同是ANSYS將結(jié)果寫到文件然后Matlab讀，或者M(jìn)atlab寫數(shù)據(jù)到文件ANSYS讀，協(xié)同是通過讀寫文件實(shí)現(xiàn)的。我找了個(gè)瞬態(tài)分析的例子，寫了APDL，將網(wǎng)上找到的ANSYS結(jié)果輸出部分代碼加了進(jìn)來，實(shí)現(xiàn)三個(gè)節(jié)點(diǎn)位移時(shí)間數(shù)據(jù)寫到文件，然后Matlab讀取該文件并作三節(jié)點(diǎn)位移-時(shí)間圖。說明一下：ANSYS生成的文件都在Matlab當(dāng)前工作目錄下。
1 Matlab的代碼（mat2ans.m）
!“D:\Program Files (x86)\Ansys Inc\v100\ANSYS\bin\intel\ansys100.exe” -b -i “E:\WP\m2a\transient.mac” -o “E:\WP\m2a\ans.out”
%save(‘E:\WP\m2a\matdat.dat’,‘data’,’-ascii’)
%load(’-ascii’,‘E:\WP\m2a\ans.rst’)
load ‘result1.txt’
t=result1(:,1);
u1=result1(:,2)*1e3;
u2=result1(:,3)*1e3;
u3=result1(:,4)*1e3;
plot(t,u1,‘k’,t,u2,‘c’,t,u3,‘g’)
legend(‘Node 100’,‘Node 146’,‘Node 200’)
xlabel(‘time/s’),ylabel(‘disp/mm’)
2 ANSYS APDL代碼（transient.mac）
finish
/cwd,‘E:\WP\M2A\ans2mat’
/filename,ansys2matlab
/title,ansys result output formatted
/units,SI
/prep7
et,1,shell63
et,2,beam4
r,1,0.02
r,2,2e-4,2e-8,2e-8,0.01,0.02
mp,ex,1,2e11
mp,prxy,1,0.3
mp,dens,1,7.83e3
rectng,0,2,0,1
kgen,2,1,4,1,-1
do,i,1,4,1
l,i,i+4
enddo
lsel,all
lesize,all,0.1
asel,all
amesh,1
mat,1
real,2
lmesh,5,8,1
fini
!
/output,cp,out!Output Info to cp.out
/debug,-1,1 !element mass & stiff matrix
!
/solu
antype,trans
alphad,5
!dkdele,all,all
ksel,s,kp,5,8,1
dk,all,all
allsel
sbctran
outres,all,all
!
time,1
autots,on
deltim,0.2,0.05,0.5,1
kbc,0
sfa,1,1,pres,10000
lswrite,1
!
time,2
lswrite,2
!
time,4
kbc,1
sfa,1,1,pres,5000
lswrite,3
!
time,6
sfa,1,1,pres,0
lswrite,4
allsel
lssolve,1,4,1
!
fini
!*
/output,term !To output window
/post1
!file,ansys2matlab,rst !結(jié)果文件
set,last
*GET,N_SET,ACTIVE,0,SET,NSET
!n_set=5
*DIM,NODES,ARRAY,3
*DIM,TIMES,ARRAY,N_SET
*DIM,RESULT,ARRAY,3,N_SET
NODES(1)=100,146,200
*DO,I,1,N_SET,1
SET,I
*GET,TI,ACTIVE,0,SET,TIME
TIMES(I)=TI
*DO,J,1,3,1
*GET,RESULT(J,I),NODE,NODES(J),U,Z
*ENDDO
*ENDDO
*cfopen,result1,txt
!*vwrite
!(’ ')
!*vwrite,nodes(1),nodes(2),nodes(3)
!(‘Time Node’,f9.0,2f15.0)
*do,i,1,N_SET
*vwrite,times(i),result(1,i),result(2,i),result(3,i)
(f8.5,8x,3f15.10)
*enddo
*cfclose
fini
！*
!dump a matrix
!/aux2
!fileaux2,ansys2matlab,emat
!form,long
!dump,all
!*
3 運(yùn)行過程概括
執(zhí)行Matlab的mat2ans.m，由于第一條命令啟動(dòng)ANSYS，控制權(quán)由Matlab交給ANSYS，ANSYS以批處理執(zhí)行transient.mac的APDL代碼，并生成想要的結(jié)果文件result1.txt然后自動(dòng)退出，Matlab得到控制權(quán)接著執(zhí)行余下的命令。
matlab/ansys協(xié)同的一個(gè)例子
4 結(jié)果展示
4.1 有限元模型位移圖
matlab/ansys協(xié)同的一個(gè)例子
4.2 ANSYS作三節(jié)點(diǎn)位移-時(shí)間曲線
matlab/ansys協(xié)同的一個(gè)例子
4.3 Matlab作三節(jié)點(diǎn)位移-時(shí)間曲線
matlab/ansys協(xié)同的一個(gè)例子http://blog.sina.com.cn/s/blog_c0d2007a0102vhf2.html
總結(jié)
以上是生活随笔為你收集整理的matlab/ansys协同的一个例子的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： LaTeX的表格标题位置
- 下一篇： Ansys节点数据批量一键导出脚本生成C