（原创）利用Matlab对ANSYS数据进行后处理实例

shifang157300

（原创）利用Matlab对ANSYS数据进行后处理实例

问题探讨：作者 - 师访，QQ - 1549221758，Email- pomato157300@126.com

Research gate：https://www.researchgate.net/profile/Fang_Shi8

仿真论坛：http://home.simwe.com/space-uid-334526.html

个人网站：http://phipsi.top/author.html

   本帖子给出一个利用Matlab对ANSYS进行后处理的2D实例。处理的数据为ANSYS的y方向位移场，Matlab绘制位移场矢量图和云图。模型宽10m，高5m，上端面受拉，下端面固定约束。本贴的关键是Matlab云图的绘制，对于已知节点坐标和节点处位移、温度、应力等标量的问题，常用的Matlab云图绘制命令contour和contourf无法使用，此时可以用griddata和pcolor绘制云图。

1导出ANSYS数据至文本文件

finish

/clear

/FILNAME,simple-case-data-matlab

/Titie,SHIFang,CUMT

/PREP7

/DSCALE,ALL,1.0

/eshape,on

!---------------------

et,1,42

!---------------------

KEYOPT,1,3,2   

KEYOPT,1,6,0

keyopt,1,10,0

!---------------------

mp,ex,1,10e9

mp,nuxy,1,0.25

!---------------------

k,1,0,0

k,2,10,0

k,3,10,5

k,4,0,5

!---------------------

a,1,2,3,4

!---------------------

LESIZE,all,1,, , ,, , ,,

!---------------------

MSHKEY,0

AMESH,all

!---------------------

DL,1,,all,      

!---------------------

SFL,3,PRES,-1e6,

!---------------------

FINISH

!---------------------

/SOLU

SOLVE

!---------------------

/post1

PLNSOL,U,Y, 0,1.0  

!------------------------------

!--------提取数据--------------

!------------------------------

*get,NodeNum,NODE,0,COUNT            !得到模型的所有节点数

*dim,nodes,array,NodeNum             !定义一个存储节点的数组

*dim,xyz,array,NodeNum,2             !定义一个存储节点坐标的二维数组

!*dim,xyz,array,NodeNum,3             !定义一个存储节点坐标的三维数组

*dim,disY,array,NodeNum,1            !定义一个存储节点Y坐标的数组                                                                        

*do,i,1,NodeNum,1

    xyz(i,1)=NX(i)                  !在xyz(i,1)中存储节点的X坐标

    xyz(i,2)=NY(i)                  !在xyz(i,1)中存储节点的Y坐标

   !xyz(i,3)=NZ(i)                  !在xyz(i,1)中存储节点的Z坐标

    *Get,disY(i),NODE,i,U,Y

   !*Get,sxy(i),NODE,i,S,XY

*enddo

!------------------------------

!--------导出文件--------------

!------------------------------

*cfopen,x-coord,txt                   !打开文件，写入数据

*vwrite,xyz(1,1)

(e16.8)

*cfclos

!

*cfopen,y-coord,txt                   !打开文件，写入数据

*vwrite,xyz(1,2)

(e16.8)

*cfclos

!

*cfopen,disY,txt                   !打开文件，写入数据

*vwrite,disY(1)

(e16.8)

*cfclos

ANSYS命令流：

2 Matlab后处理

m文件内容如下：

%Written By:: Shi Fang, China University of Mining and Technology (2013)

%Email: betaxfem@gmail.com or pomato157300@126.com

%QQ: 1549221758

clearall; close all; clc; format compact; format long;

%Settings

%------------------------------

SizeFont=20;

set(0,'defaultfigurecolor','w')

%Read data

%------------------------------

disp('Readinginput file....')

x  = load('x-coord.txt');

x=x';

y  = load('y-coord.txt');

y=y';

Uy= load('disY.txt');

Uy=Uy';

disp('Readdone.')

%Plot the vector

%------------------------------

disp('Plottingthe vector....')

num=size(x,2);

fori=1:num

    Ux(i)=0;

end

quiver(x,y,Ux,Uy);

disp('Plotdone.')

%Plot the contour

%------------------------------

disp('Plottingthe contour....')

pp1=0:0.1:10;              % Range of the x axial

pp2=0:0.01:5;             % Range of the y axial

[X,Y]=meshgrid(pp1,pp2);

T=griddata(x,y,Uy,X,Y,'cubic');

%T(pp>2,pp<3)=nan;         % Blank zone if necessary

figure;hold on;

title('\itDisplacement y','FontName','Times New Roman','FontSize',SizeFont);

pcolor(X,Y,T);

axisequal;  

axis([min(x)max(x) min(y) max(y)]);

set(gca,'XTick',[],'YTick',[],'XColor','w','YColor','w');

colorbar('FontAngle','italic','FontName','TimesNew Roman','FontSize',SizeFont);

shadinginterp;

disp('Plotdone.')

m文件：

Matlab绘制的y方向位移矢量图：

Matlab绘制的y方向位移云图：

我的帖子汇总：

（1）扩展有限元法(XFEM)漫谈（原创_持续更新）：http://forum.simwe.com/thread-1245357-1-1.html

（2）ANSYS UPFs二次开发教程（原创，持续更新）：http://forum.simwe.com/thread-1021888-1-1.html

（3）比特币(Bitcoin)价格波动曲线监控程序_基于Python编写__原创：http://forum.simwe.com/thread-1254734-1-1.html

（4）地震作用下水箱晃动分析（原创，附命令流和K文件）：http://forum.simwe.com/thread-1098856-1-1.html

（5）ANSYS常用单元特性总结及简单实例（原创）：http://forum.simwe.com/thread-1019487-1-1.html

（6）ANSYS计算受压闭合裂纹的应力强度因子实例命令流（原创）：http://forum.simwe.com/thread-1102153-1-1.html

（7）（原创）利用Matlab对ANSYS数据进行后处理实例：http://forum.simwe.com/thread-1094497-1-1.html

（8）（原创）ANSYS跌落分析全命令流模板：http://forum.simwe.com/thread-1093360-1-1.html

（9）钢球入水后反弹CEL流固耦合3D模拟，附CAE文件（原创）:http://forum.simwe.com/thread-1107014-1-1.html

（10）ANSYS16.0新增扩展有限元XFEM裂纹扩展仿真简介（原创）：http://forum.simwe.com/thread-1121373-1-1.html

（11）FEAPpv简介、Linux下的编译及运行方法_原创: http://forum.simwe.com/thread-1136956-1-1.html

（12）基于ANSYS的结构地震分析教程 (静力+模态+反应谱+时程) http://forum.simwe.com/thread-1176187-1-1.html

wegean

这方面的资料挺少的，感谢您的分享！

shijc

很需要，谢谢分享！

陆天琪

厉害啊！多学习。

zcc861007

好帖！:lol:lol

tangwenwen

你好，我用实体单元画了一护栏，结果已分析出来，我想请教一下怎么用matlab得到位移-梁长的曲线，例如梁长为10cm处的位移值，20cm处的位移值，等等等。我是初学者ANSYS用的是菜单法，所以没有命令流。谢谢

zcc861007

非常感谢楼主fenxiAANG

jibinnici

不错的资料，很实用，谢谢楼主无私分享

nimo479

谢谢楼主

独自凭栏

高手啊，什么都会用

shifang157300

独自凭栏 发表于 2014-3-20 17:20
高手啊，什么都会用

谢谢支持

springer_

顶楼主，还出了本ansys二次开发的书，

shifang157300

springer_ 发表于 2014-7-26 11:07
顶楼主，还出了本ansys二次开发的书，

sjzdh

不错的资料

nieyiguo

ansys后处理的时候，如果需要处理数据，太慢了，像是提刚阵的时候，还是matlab快一点

wangfire

多谢楼主

weldinghui

楼主好强大，mark了以后再来有用到再来学习

回头看世界

默默关注中

hangxiaolin

学习了~

beifengabc

支持下，先下啦，慢慢研究。