matlab fprintf_Matlab给Ansys助攻

最新推荐文章于 2023-02-16 16:24:45 发布

weixin_39963465

最新推荐文章于 2023-02-16 16:24:45 发布

阅读量116

点赞数

文章标签： matlab fprintf matlab 写入txt matlab写入txt

大家好，今天二狗准备介绍一下如何利用matlab写一些简单的脚本去解决一些复杂的事情，比如为大型有限元软件Ansys服务(类比)包括数值模拟的前处理和后处理，都大有作为。

用过Ansys(没有Ctrl+z)的童鞋知道，它有GUI和命令流俩个操作模式，优缺点就不多说了。前处理就是建模，点线面体几何拓扑。如果有几千个坐标点(XYZ)，成点很麻烦，线面体更麻烦。简单，可回退，可传播分享的一种方式就是利用脚本批量写入命令流到文本文件(Txt-dat-sca-csv...)

下面是导入点的源码，欢迎参考:

%三维数据到ANSYS[filename,pathname,filterindex]=uigetfile({'*.csv';'*.xlsx';'*.mat'},'选择数据文件');    str=[pathname filename]; if filename==0 % 取消return;endif filterindex==2 c_temp=importdata(str);c_close=c_temp.Sheet1;elsec_temp=importdata(str);c_close=c_temp;end[hang,lie]=size(c_close);if lie==3mingzinn=char(filename);ll=strfind(mingzinn,'.');ll0=mingzinn(1:ll-1);nametxt=[ll0,'2ANSYS','.dat'];%扩展名可改fp=fopen(nametxt,'a');str111='k , ';str222=' , ';str333=' ';str444='LP1 = _return';%str555='L,LP1,LP2';%str666='LP1 = LP2';%str777='LP2 = _return';for i=1:1:hang   if i==1    fprintf(fp,'%s\r\n','/prep7');      fprintf(fp,'%s\r\n','CSYS,0');    str1=[str111,num2str(i),str222,num2str(c_close(1,1)),str222,num2str(c_close(1,2)),str222,num2str(c_close(1,3))];    fprintf(fp,'%s\r\n',str1);   else    str2=[str111,num2str(i),str222,num2str(c_close(i,1)),str222,num2str(c_close(i,2)),str222,num2str(c_close(i,3))];       fprintf(fp,'%s\r\n',str2);      end         endfprintf(fp,'%s\r\n','!File generation using Matlab.m');fclose(fp);sound(sin(2*pi*25*(1:4000)/100));elsemsgbox('数字矩阵必须为x*3形态');    end

这个是导入点自动成线的源码:

[filename,pathname,filterindex]=uigetfile({'*.sca';'*.csv';'*.xlsx';'*.mat'},'选择数据文件');    str=[pathname filename]; if filename==0 % 取消return;endif filterindex==3 c_temp=importdata(str);c_close=c_temp.Sheet1;elseif filterindex==1 c_temp=importdata(str);c_close=c_temp.data;  c_close=c_close(:,2:4);elseif filterindex==2 c_temp=importdata(str);c_close=c_temp;elseif filterindex==4 c_temp=importdata(str);c_close=c_temp;end[hang,lie]=size(c_close);if lie==3mingzinn=char(filename);ll=strfind(mingzinn,'.');ll0=mingzinn(1:ll-1);nametxt=[ll0,'-Ansys','.dat'];%扩展名可改fp=fopen(nametxt,'a');str111='K,, ';str222=' , ';str333=' ';str444='LP1 = _return';str555='L,LP1,LP2';str666='LP1 = LP2';str777='LP2 = _return';for i=1:1:hang   if i==1    fprintf(fp,'%s\r\n','/prep7');      str1=[str111,num2str(c_close(1,1)),str222,num2str(c_close(1,2)),str222,num2str(c_close(1,3)),str333];       fprintf(fp,'%s\r\n',str1);    fprintf(fp,'%s\r\n',str444);   else    str2=[str111,num2str(c_close(i,1)),str222,num2str(c_close(i,2)),str222,num2str(c_close(i,3)),str333];       fprintf(fp,'%s\r\n',str2);       fprintf(fp,'%s\r\n',str777);       fprintf(fp,'%s\r\n',str555);    fprintf(fp,'%s\r\n',str666);   end         endfprintf(fp,'%s\r\n','!File generation using Matlab.m');fprintf(fp,'%s\r\n','!By 阿辉');fclose(fp);sound(sin(2*pi*25*(1:4000)/100));elsemsgbox('数字矩阵必须为x*3形态(不算 *.sca)');    end

成面与体，看自己实际需求，也可直接在ansys里写命令流。Ansys求解后的结果该如何导出呢，举个栗子，比如各节点的x方向位移(需要把Ansys里对应的List的文本复制出来，比如到csv)，源码如下:

% load a88[hang,lie]=size(a88);ping=[];for i=1:1:hangjiancha=a88(i,1);kankan=str2num(jiancha);if isempty(kankan)==1    else    ping=[ping;kankan]; endh=waitbar(i/hang);enddelete(h);msgbox('OK');

这样得到的就是数据文件而不是字符串，且没有空白，没有注释。

对于Ansys，懂的看它的Log日志文件很重要，但是此文件里有很多的无效命令，就是前面加了!，感叹号是注释的意思，如何把这些非实质性操作清空呢，以下是二狗写的几句代码：

[h,l]=size(log);xin=[];nametxt=['log1','-Ansys','.dat'];fp=fopen(nametxt,'a');for p=1:1:h   jiancha= log(p,1);   l8= strlength(jiancha);   if l8>1      weizhi=strfind(jiancha,'!');      [k1,k2]=size(weizhi);     if k1>0 && weizhi(1,1)==1              else       xin=[xin;jiancha];        fprintf(fp,'%s\r\n',jiancha);       end   else          end    jdt=waitbar(p/h);endfprintf(fp,'%s\r\n','!File idealization using Matlab.m');%写下来源  嘿嘿delete(jdt);fclose(fp);sound(sin(2*pi*25*(1:4000)/100));

之前是：