matlab/ansys协同的一个例子

前段时间找了点Matlab与ANSYS协同工作的资料,在我看来所谓协同是ANSYS将结果写到文件然后Matlab读,或者Matlab写数据到文件ANSYS读,协同是通过读写文件实现的。我找了个瞬态分析的例子,写了APDL,将网上找到的ANSYS结果输出部分代码加了进来,实现三个节点位移时间数据写到文件,然后Matlab读取该文件并作三节点位移-时间图。说明一下:ANSYS生成的文件都在Matlab当前工作目录下。

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 !结果文件
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 运行过程概括
执行Matlab的mat2ans.m,由于第一条命令启动ANSYS,控制权由Matlab交给ANSYS,ANSYS以批处理执行transient.mac的APDL代码,并生成想要的结果文件result1.txt然后自动退出,Matlab得到控制权接着执行余下的命令。
在这里插入图片描述
matlab/ansys协同的一个例子

4 结果展示

4.1 有限元模型位移图
在这里插入图片描述
matlab/ansys协同的一个例子

4.2 ANSYS作三节点位移-时间曲线
在这里插入图片描述
matlab/ansys协同的一个例子
4.3 Matlab作三节点位移-时间曲线
在这里插入图片描述
matlab/ansys协同的一个例子http://blog.sina.com.cn/s/blog_c0d2007a0102vhf2.html

                                </div><link href="https://csdnimg.cn/release/phoenix/mdeditor/markdown_views-b6c3c6d139.css" rel="stylesheet"><div class="more-toolbox"><div class="left-toolbox"><ul class="toolbox-list"><li class="tool-item tool-active is-like "><a href="javascript:;"><svg class="icon" aria-hidden="true"><use xlink:href="#csdnc-thumbsup"></use></svg><span class="name">点赞</span><span class="count"></span></a></li><li class="tool-item tool-active is-collection "><a href="javascript:;" data-report-click="{&quot;mod&quot;:&quot;popu_824&quot;}"><svg class="icon" aria-hidden="true"><use xlink:href="#icon-csdnc-Collection-G"></use></svg><span class="name">收藏</span></a></li><li class="tool-item tool-active is-share"><a href="javascript:;"><svg class="icon" aria-hidden="true"><use xlink:href="#icon-csdnc-fenxiang"></use></svg>分享</a></li><!--打赏开始--><!--打赏结束--><li class="tool-item tool-more"><a><svg t="1575545411852" class="icon" viewBox="0 0 1024 1024" version="1.1" xmlns="http://www.w3.org/2000/svg" p-id="5717" xmlns:xlink="http://www.w3.org/1999/xlink" width="200" height="200"><defs><style type="text/css"></style></defs><path d="M179.176 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5718"></path><path d="M509.684 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5719"></path><path d="M846.175 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5720"></path></svg></a><ul class="more-box"><li class="item"><a class="article-report">文章举报</a></li></ul></li></ul></div></div><div class="person-messagebox"><div class="left-message"><a href="https://blog.csdn.net/qq_33588707"><img src="https://profile.csdnimg.cn/F/A/8/3_qq_33588707" class="avatar_pic" username="qq_33588707"><img src="https://g.csdnimg.cn/static/user-reg-year/1x/4.png" class="user-years"></a></div><div class="middle-message"><div class="title"><span class="tit"><a href="https://blog.csdn.net/qq_33588707" data-report-click="{&quot;mod&quot;:&quot;popu_379&quot;}" target="_blank">小黄三爷</a></span></div><div class="text"><span>发布了3 篇原创文章</span> · <span>获赞 2</span> · <span>访问量 2143</span></div></div><div class="right-message"><a href="https://im.csdn.net/im/main.html?userName=qq_33588707" target="_blank" class="btn btn-sm btn-red-hollow bt-button personal-letter">私信</a><a class="btn btn-sm  bt-button personal-watch" data-report-click="{&quot;mod&quot;:&quot;popu_379&quot;}">关注</a></div></div></div>@[TOC](这里写自定义目录标题)

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

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

相关文章

security center拒绝访问_Steam被曝出0day提权漏洞,但厂商拒绝修复

近期&#xff0c;全球流行的Steam游戏客户端被曝出0day提权漏洞&#xff0c;影响全球一亿多Steam用户。该漏洞可让机器上的低权限用户以SYSTEM权限运行程序。这意味着恶意软件很可能利用这个漏洞对受害者的机器进行深度破坏。考虑到Steam是一款总用户1.2亿&#xff0c;日常在线…

赵老师_01

1、赵4 老师 课程&#xff0c;第2课 (VS2010) (Console的 空项目) 1.1、项目-->属性-->配置属性-->C/C -->预编译头&#xff1a;  不使用预编译头 1.2、项目-->属性-->配置属性-->C/C -->代码生成-->基本运行时检查&#xff1a;  默认值 1.3、项…

人工智能第一股“旷视科技”赴港递交IPO材料 半年亏损52亿估值280亿

来源&#xff1a; 人工智能和大数据8月25日&#xff0c;人工智能领域“独角兽”旷视科技正式向港交所递交IPO材料。招股书显示&#xff0c;旷视拟采用“同股不同权”方式上市&#xff0c;“同股不同权”股东为创始人兼CEO印奇以及两位联合创始人唐文斌、杨沐&#xff0c;三人占…

C++中set按降序排序

set一般插入元素时&#xff0c;默认使用关键字类型的< 运算符来比较两个关键字&#xff0c;故一般插入后为升序&#xff0c;举个例子&#xff1a; int main(int argc, char* argv[]) {set<int> s;s.insert(2);s.insert(1);s.insert(3);set<int>::iterator it …

matlab调用ansys-------1

1.版本&#xff1a;MATLAB R2014a ANSYS14.5 / MATLAB R2017a ANSYS17.0 2.思路&#xff1a;将MATLAB中的数据输出到shuju.txt中&#xff0c;用ANSYS读取shuju.txt后将计算计算结果输出到jieguo.txt中并将其返回到MATLAB中进行运算分析。 3.程序&#xff1a;在MATLAB中输入 …

python可以代替plc吗_PLC可以代替安全控制器吗?电气工程师告诉你

导读 国内的企业越来越重视设备的安全防护&#xff0c;可事故却频频发生&#xff0c;这是为什么呢&#xff1f;主要还是源自于大家对于安全防护知识上存在的诸多误区&#xff0c;比如很多人认为&#xff0c;稳定性高、性能强的PLC即可替代安全控制器的作用&#xff0c;其实这就…

php抽象 与接口

抽象方法和抽象类 PHP5支持抽象类和抽象方法。抽象类不能直接被实例化&#xff0c;你必须先继承该抽象类&#xff0c;然后再实例化子类。抽象类中 至少要包含一个抽象方法。如果类方法被声明为抽象的&#xff0c;那么其中就不能包括具体的功能实现。 继承一个抽象类的时候&…

C++ string转字符串

假设有一个string字符串&#xff0c;定义为string s "123"&#xff0c;则以下两种方法均可将其转换为char* s.c_str()s.data() 区别在于data除了返回字符串内容外&#xff0c;不附加结束符’\0’&#xff0c;而c_str()返回一个以‘\0’结尾的字符数组。 注意&…

Nature展示迄今为止最详细的“人脑零部件清单”

来源&#xff1a;中国生物技术网近日&#xff0c;美国艾伦脑科学研究所的一项新研究为我们展示了迄今为止最详细的“人脑零部件清单”&#xff0c;阐明了人类大脑皮层的细胞结构&#xff0c;它对于理解我们的认知能力和对疾病的易感性至关重要。该研究近日已发表在《Nature》上…

matlab调用ANSYS--------2

最近使用ANSYS计算温度场方面的东西&#xff0c;每次都要手动运行&#xff0c;自己处理计算结果。可怜无助的我深感劳累啊&#xff0c;最主要的是怎么能把我宝贵的时间用在重复的劳动中。so&#xff0c;出于懒的需要&#xff0c;查询了一些关于matlab调用ANSYS的资料&#xff0…

beautifulsoup网页爬虫解析_爬虫第三天——初识Xpath

解析神器Xpath&#xff1a;1. 什么是XpathXPath即为XML路径语言&#xff08;XML Path Language&#xff09;&#xff0c;它是一种用来确定XML文档中某部分位置的语言。XPath基于XML的树状结构&#xff0c;提供在数据结构树中找寻节点的能力。起初XPath的提出的初衷是将其作为一…

全球知名大学课件下载地址汇总

本文集合了美国众多知名大学的网上课程课件&#xff0c;不用花一分钱就能获得知名大学的前沿知识&#xff01;聪明的同学早就收藏转发啦 美国大学课件下载地址汇总 1、加州大学伯克利分校http://webcast.berke…

C++将01数组转换为二进制对应的数值

假设有一个32维的01数组array&#xff0c;假设array的值为 01001010001011110100101000101111 如果需要将其转换为用该数组值代表的二进制所对应的数值&#xff0c;即通过该数组得到01001010001011110100101000101111这个int型数值&#xff08;int型有4字节&#xff0c;一共4…

机器人4.0时代来临!四大核心技术助推大规模商用部署

来源&#xff1a;北京物联网智能技术应用协会当前&#xff0c;全球机器人市场规模持续扩大&#xff0c;工业机器人市场增速稳定&#xff0c;服务机器人增速突出。2018 年&#xff0c;全球机器人市场规模达 298.2 亿美元&#xff0c; 2013-2018 年的平均增长率约为 15.1%。 在装…

Chapter2 MSP430硬件结构

1、总架构 MSP430单片机采用冯诺依曼结构&#xff0c;包含16位RISC CPU、存储器、片内外设、时钟系统、仿真系统以及数据总线和地址总线。 2、MSP430X CPU(CPUX) DB    AB          ALU位数 MSP430F1XX 16位  16位(最大寻址64KB) …

怎么用python编写个apk_新手用tkinter写了一个APK提取工具

[Python] 纯文本查看 复制代码""" author:qh datetime:2019-3-15 mood:<(*&#xffe3;▽&#xffe3;*)/ """ import os import base64 import requests import json import shutil import hashlib from tkinter import * from tkinter impo…

Matlab给Ansys助攻

大家好&#xff0c;今天二狗准备介绍一下如何利用matlab写一些简单的脚本去解决一些复杂的事情&#xff0c;比如为大型有限元软件Ansys服务&#xff08;类比&#xff09;包括数值模拟的前处理和后处理&#xff0c;都大有作为。用过Ansys&#xff08;没有Ctrlz&#xff09;的童鞋…

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

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

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; …