shifang157300 发表于 2013-12-1 13:19:12

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

本帖最后由 shifang157300 于 2017-4-16 17:30 编辑

（原创）利用Matlab对ANSYS数据进行后处理实例问题探讨：作者 - 师访，QQ - 1549221758，Email- pomato157300@126.comResearch 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,0keyopt,1,10,0!---------------------mp,ex,1,10e9mp,nuxy,1,0.25!---------------------k,1,0,0k,2,10,0k,3,10,5k,4,0,5!---------------------a,1,2,3,4!---------------------LESIZE,all,1,, , ,, , ,, !---------------------MSHKEY,0AMESH,all!---------------------DL,1,,all,      !---------------------SFL,3,PRES,-1e6,!---------------------FINISH!---------------------/SOLUSOLVE!---------------------/post1PLNSOL,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)
*cfclosANSYS命令流：

2 Matlab后处理
m文件内容如下：%Written By:: Shi Fang, China University of Mining and Technology (2013)%Email: betaxfem@gmail.com or pomato157300@126.com%QQ: 1549221758clearall; 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;endquiver(x,y,Ux,Uy);disp('Plotdone.') %Plot the contour%------------------------------disp('Plottingthe contour....') pp1=0:0.1:10;            % Range of the x axialpp2=0:0.01:5;             % Range of the y axial=meshgrid(pp1,pp2);T=griddata(x,y,Uy,X,Y,'cubic');%T(pp>2,pp<3)=nan;         % Blank zone if necessaryfigure;hold on;title('\itDisplacement y','FontName','Times New Roman','FontSize',SizeFont);pcolor(X,Y,T);axisequal;axis();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 发表于 2013-12-1 20:18:04

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

shijc 发表于 2013-12-3 21:36:57

很需要，谢谢分享！

陆天琪 发表于 2013-12-4 15:33:04

厉害啊！多学习。

zcc861007 发表于 2013-12-8 09:04:42

好帖！:lol:lol:lol

tangwenwen 发表于 2013-12-8 15:53:44

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

zcc861007 发表于 2014-1-12 01:57:37

非常感谢楼主fenxiAANG

jibinnici 发表于 2014-3-6 17:53:39

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

nimo479 发表于 2014-3-20 11:35:07

谢谢楼主

独自凭栏 发表于 2014-3-20 17:20:03

高手啊，什么都会用:(:D:lol:kiss:

shifang157300 发表于 2014-3-20 18:24:26

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

谢谢支持

springer_ 发表于 2014-7-26 11:07:49

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

shifang157300 发表于 2014-7-27 08:39:33

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

:lol

sjzdh 发表于 2014-8-26 17:47:34

不错的资料

nieyiguo 发表于 2014-8-27 19:48:52

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

wangfire 发表于 2014-9-23 10:09:16

多谢楼主

weldinghui 发表于 2014-9-23 11:08:39

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

回头看世界 发表于 2014-11-5 15:03:35

默默关注中

hangxiaolin 发表于 2014-11-20 00:12:08

学习了~

beifengabc 发表于 2015-4-1 14:30:25

支持下，先下啦，慢慢研究。
页: [1] 2
查看完整版本: （原创）利用Matlab对ANSYS数据进行后处理实例