Matlab给Ansys助攻

大家好,今天二狗准备介绍一下如何利用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;
end
if filterindex==2 
c_temp=importdata(str);
c_close=c_temp.Sheet1;
else
c_temp=importdata(str);
c_close=c_temp;
end
[hang,lie]=size(c_close);
if lie==3
mingzinn=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:hangif i==1fprintf(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);elsestr2=[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         
end
fprintf(fp,'%s\r\n','!File generation using Matlab.m');
fclose(fp);
sound(sin(2*pi*25*(1:4000)/100));
else
msgbox('数字矩阵必须为x*3形态');    
end

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

[filename,pathname,filterindex]=uigetfile({'*.sca';'*.csv';'*.xlsx';'*.mat'},'选择数据文件');str=[pathname filename]; 
if filename==0 % 取消
return;
end
if 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==3
mingzinn=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:hangif i==1fprintf(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);elsestr2=[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         
end
fprintf(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));
else
msgbox('数字矩阵必须为x*3形态(不算 *.sca)');    
end

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

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

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

   对于Ansys,懂的看它的Log日志文件很重要,但是此文件里有很多的无效命令,就是前面加了!,感叹号是注释的意思,如何把这些非实质性操作清空呢,以下是二狗写的几句代码:
[h,l]=size(log);
xin=[];
nametxt=['log1','-Ansys','.dat'];
fp=fopen(nametxt,'a');
for p=1:1:hjiancha= log(p,1);l8= strlength(jiancha);if l8>1weizhi=strfind(jiancha,'!');[k1,k2]=size(weizhi);if k1>0 && weizhi(1,1)==1elsexin=[xin;jiancha]; fprintf(fp,'%s\r\n',jiancha);  endelseend jdt=waitbar(p/h);
end
fprintf(fp,'%s\r\n','!File idealization using Matlab.m');%写下来源  嘿嘿
delete(jdt);
fclose(fp);
sound(sin(2*pi*25*(1:4000)/100));

在这里插入图片描述
是不是很方便哦⊙∀⊙!

在学习工作中,有很多操作可以用代码去准确高效的实现,也可直接把其打包成Exe,若需要。嗯。为达目的,“不择手段”。

二狗就写这么多,文笔不好,见谅,有错误的地方请告诉二狗。解释一下,提到的Ansys是支持数据导入导出的,二狗觉得,不好用。这些操作用vb、c、Python、office的宏等等都可以实现,但是matlab,除了生孩子,无所不能,哈哈ヾノ≧∀≦)o

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/490522.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

归并排序的时间复杂度为什么为nlogn

归并排序的递归过程如下,该递归树的高度为log2n(计算过程:假设待排序的数组元素个数为n,设高度为x,x意味着n个元素需要连续二分x次才剩下1个元素,即n/2^x1,xlog2n),每一层…

reflection

function $f new ReflectionFunction(funcName);var_dump($f->getFileName());exit; method $func new ReflectionMethod(classname,function); var_dump($func->getFileName()); 转载于:https://www.cnblogs.com/fenle/p/5573841.html

python二分法查找算法_排序算法和二分法查找

请填充代码&#xff0c;使mySort()能使传入的参数按照从小到大的顺序显示出来。 function mySort() { var tags new Array(); for (var i 0; i < arguments.length; i) { tags.push(arguments[i]); } tags.sort(function sortNum(a, b) { return a - b; }); return tags; …

深度解读达芬奇架构:华为AI芯片的“秘密武器”

本文转自华为官网2019年6月&#xff0c;华为发布全新8系列手机SoC芯片麒麟810&#xff0c;首次采用华为自研达芬奇架构NPU&#xff0c;实现业界领先端侧AI算力&#xff0c;在业界公认的苏黎世联邦理工学院推出的AI Benchmark榜单中&#xff0c;搭载麒麟810的手机霸榜TOP3&#…

FGSM(Fast Gradient Sign Method)python实现

有现成的库cleverhans可以调用 https://github.com/cleverhans-lab/cleverhans 这个库是用tensorflow实现的&#xff0c;但也实现了pytorch版本 因为我调用的是pytorch版本&#xff0c;我接下来大概讲一下如何使用 他们给出了一个调用pytorch版本的例子&#xff1a; https…

记忆化搜索 codevs 2241 排序二叉树

codevs 2241 排序二叉树 ★ 输入文件&#xff1a;bstree.in 输出文件&#xff1a;bstree.out 简单对比时间限制&#xff1a;1 s 内存限制&#xff1a;128 MB 【问题描述】 一个边长为n的正三角形可以被划分成若干个小的边长为1的正三角形&#xff0c;称为单位三角形。 …

窥探当今AI芯片中的类脑模型

来源&#xff1a;脑人言撰文丨邓 磊&#xff08;加州大学-圣塔芭芭拉分校UCSB 博士后&#xff09;责编丨高茂森 Soma排版丨夏獭科学研究的动机&#xff0c;并非仅仅是出于能看到多酷炫的应用&#xff0c;而是因为我们尚有诸多未知世界需要探寻&#xff0c;这是一切应用的前提。…

mysql调优_MYSQL企业常用架构与调优经验分享

一、选择Percona Server、MariaDB还是MYSQL1、Mysql三种存储引擎MySQL提供了两种存储引擎&#xff1a;MyISAM和 InnoDB&#xff0c;MySQL4和5使用默认的MyISAM存储引擎。从MYSQL5.5开始&#xff0c;MySQL已将默认存储引擎从MyISAM更改为InnoDB。MyISAM没有提供事务支持&#xf…

shutdown-t-s3600定时关机

winr定时关机shutdown-t-s3600 单位 秒

pycharm无缘无故打开Nosetest

一个正常的代码&#xff0c;本来前几天可以正常run&#xff0c;突然有一天点击run之后不是运行程序&#xff0c;而是报错一堆&#xff1a; Launching Nosetest with arguments。。。 usage: _jb_nosetest_runner.py 。。。 _jb_nosetest_runner.py: error: unrecognized arg…

python按比例生成数据组_基于python中的一个值生成“正态分布”数据

通过施加总和temp100你介绍的依赖至少两个数据点之间&#xff0c;因此无法建立一套独立的采样随机数据点。 一个简单的例子&#xff1a; 想象一下投币。系统中的随机性正好是一个二进制结果&#xff0c;或1位。 想象一下两个硬币翻转。系统中的随机性恰好是两个二进制结果或2位…

SpringMVC学习--参数绑定

spring参数绑定过程 从客户端请求key/value数据&#xff0c;经过参数绑定&#xff0c;将key/value数据绑定到controller方法的形参上。springmvc中&#xff0c;接收页面提交的数据是通过方法形参来接收。而不是在controller类定义成员变更接收。 默认支持的类型直接在controll…

matlab与ie交互

clc clear ieactxserver(internetexplorer.application); ie.Navigate(网址); while ~strcmp(ie.readystate,READYSTATE_COMPLETE) pause(.01) end loginnameie.document.getElementById(输入用户名的LABEL的id); loginname.value用户名; passwordie.document.getElementById(输…

MySQL windows安装

https://www.cnblogs.com/winton-nfs/p/11524007.html

现代物理学7大经典问题,你能理解几个?或许一个都不能理解

来源&#xff1a;搜狐新闻算法与数学之美相对论相对论是物理学中两大著名理论之一&#xff0c;两者都是阿尔伯特爱因斯坦提出的。1905年爱因斯坦出版了狭义相对论&#xff0c;后者确定最终宇宙速度极限&#xff1a;光速。并称时间因某物体移动的速度而实现加速或者减慢。1916年…

unity asset store下载不了_Unity手游实战:从0开始SLG——资源管理系统-基础篇(三)AssetBundle原理...

先用一句话介绍一下AssetBundle吧。AssetBundle系统提供了一种压缩文件的格式&#xff0c;可以把1到多个文件进行索引和序列化。Unity项目在交付安装之后&#xff0c;会通过AssetBundle对不包含代码的资源进行更新。这就允许开发人员先提交一个小的应用程序包&#xff0c;将运行…

《人月神话》阅读笔记(三)

我又一次在相关书籍看到了文档的重要性 数据的表现形式是编程的根本。创造出自精湛的技艺&#xff0c;精炼、充分和快速的程序也是如此。 技艺改进的结果往往是战略上的突破&#xff0c;而不仅仅是技巧上的提高。 这种战略上突破有时是一种新的算法&#xff0c;如快速傅里叶变…

matlab中GUIDE的UItable居中方法

需要借助JAVA&#xff0c;代码示意如下&#xff1a; import javax.swing.*; %导入javax.swing类 table findjobj(handles.uitable1); %需要findjobj.m文件&#xff0c;matlab中是没有的 table1 get(table,Viewport); jtable get(table1,View); renderer jtable.getCellRe…

MySQL Workbench运行脚本

首先在SCHEMAS区域右击&#xff0c;选择创建“create schema” 填入自定义名称&#xff0c;点击“apply” 在弹出来的框中点击“apply”&#xff0c;再在弹出来的框中点击finish 点击菜单栏“file”——“run SQL script”&#xff0c;选择sql文件&#xff0c;在下面选择刚刚创…

自动驾驶事故权威揭秘:算法和软件错误太多

本文来源&#xff1a;AutoR智驾自动驾驶发展就好比是登月计划。从传感器到人工智能&#xff0c;经典的电子供应链已经形成了一个协作矩阵&#xff0c;致力于实现自动驾驶车辆的安全性。为此&#xff0c;还需进行大量硬件和软件开发工作&#xff0c;以确保驾驶员、乘客和行人受到…