| 查看: 6428 | 回复: 22 | ||||||||||
| 【奖励】 本帖被评价6次，作者cpu3012增加金币 4.8 个 | ||||||||||
[资源]
一种简单已用的热加工图画法
|
||||||||||
|
想必很多研究工作者在得到了热压缩模拟实验的真应力应变曲线后，想计算并绘制出热加工图，但是具体的操作步骤又不是那么轻松吧！分享给大家一种直接由温度、ln应变速率、对应应变量的应力值三者直接计算并绘制出的热加工图的方法！如下： Matlab平台的具体计算程序 应变量为0.5时的案例 clc;clear; X=[600,660,730,775,860,980] %X轴,温度 Y=[-3,-2,-1,0,1,1.70] %Y轴,log应变速率 Z=[108.3,78.5,57.1,36.1,24.3,10.9, 152.4,112.9,71.7,48.8,35.4,11.4, 184.9,154.4,112.9,77.8,52.6,15.2, 239.4,215.4,152.0,107.6,67.6,34.7, 269.8,249.1,221.7,178.3,110.1,61.8, 280.5,264.1,247.1,211.8,112.6,89.7] %X列，Y行 Zlog=log10(Z) %log应力 i=50 %设置插值矩阵的行数 j=50 %设置插值矩阵的列数 Xi=(600: (980-600)/(j-1):980) %赋值给X插值点 Ylogi=(-3: (1.70-(-3))/(i-1):1.70) %赋值给Y插值点 Yi=10.^Ylogi %应变速率的插值向量，10的Ylogi元素次方 Zlogi=interp2(X,Y,Zlog,Xi,Ylogi','spline') %用spline方法插值(log应力)矩阵 Zi=interp2(X,Y,Z,Xi,Ylogi','spline') %用spline方法插值应力矩阵 [Fx,M]=gradient(Zlogi,1,(1.70-(-3))/(i-1)) %计算log力/log速率=M。如果F是二维矩阵，返回F的二维数值梯度。[FX,FY]=gradient(F,HX,HY)。HX,HY参数表示各方向相邻两点的距离。 S=cumtrapz(Yi,Zi) %通过递归循环计算积分,计算累积梯形积分,应力在log应变速率上的点积分 P=zeros(i,j) %设置全0矩阵 Es=-3 %赋值log应变速率最小值 Smin=Zi(1,: ).*10^Es./(M(1,: )+1) %计算应变速率最小时的第一项定值,(3,: )代表第三行，所有的列。(:,3)代表所有行,第三列。(:,1,1)代表所有第一维,第二维和第三维的第一个元素。 for n=1:j G(:,n)=S(:,n)+Smin(n) %每项的累积积分值加上最小值-G值 P(:,n)=Zi(:,n).*Yi(: ) %每项应力和应变速率的乘积-P值 end eta=2*(1-G./P)*100 %计算η值矩阵 ksai=2*M./(eta/100)-1 %计算ζ值矩阵 contour(Xi,Ylogi,eta) %η二维图，contour(Xi,Ylogi,eta,10) ,10表示10根等高线，不加10表示自动 surf(Xi,Ylogi,eta) %η三维图，三维图容易造成Matlab死机 附带告诉大家一些Matlab常用的命令，可以用到的 hold on命令保持当前图形，并将此后绘制的所有图形添加到当前的图形窗口中，如果新的曲线所对应的坐标极限值与原图不一致，系统将自动进行调整。 hold off取消hold on的效果，hold命令在hold on与hold off之间切换。 cabel 使用等值矩阵生产标注，并将标注显示在当前图形 contour显示矩阵Z的二维等高线图 contour3显示矩阵Z的三维等高线图 contourf显示矩阵Z的二维等高线图，并在各等高线之间用实体颜色填充 contourc用于计算由其他等高线函数调用的等值矩阵 meshc创建一个与二维等高线图匹配的网线图 surfc创建一个与二维等高线图匹配的曲面图 阴影fill函数 clc; clear all; close all; figure; box on; hold on; contour(peaks(20),10); fill([1:2:20 20:-2:1], [repmat(18, 1, 10) repmat(20, 1, 10)], 'r'); [ Last edited by cpu3012 on 2011-9-25 at 16:07 ] 由于很久脱离科研，很久没来未来小木虫，看到有些童鞋在楼下回复怎么弄出阴影。 本人回忆下仅存的一些记忆，当初是用了hold on这个命令，然后图阴影叠加上去的，也就是说两个图叠加的。 [ Last edited by cpu3012 on 2015-8-31 at 10:49 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 软件与教程 | 科研工具资源Nick | 资源收集 | 金属材料科研工具 |
| 加工图 | 模拟仿真技术贴 | 博士后 | 热加工图 |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助各位，这样的柱形图画法，谢谢
已经有4人回复
- 请高手帮忙，热加工图绘制,并把两图合并起来。
已经有12人回复
- 求助热处理方面的一个问题和一个冷热加工的问题？？？
已经有12人回复
- 材料热加工领域的物理模拟技术（已查重）
已经有103人回复
- 如何描述/表征金属热加工过程中的动态回复现象？
已经有4人回复
- 有关纳米材料的文献中的各种示意图的画法，你知道吗
已经有38人回复
- 只有压缩实验数据，请问作热加工图得花多长时间一般！？
已经有7人回复
- 热加工图的使用
已经有6人回复
- 用matlab做热加工图
已经有12人回复
- 求助一个origin图的画法
已经有8人回复
- [求助]热加工图制作的具体步骤
已经有9人回复
- 【讨论】软质电加热套的正确使用方法
已经有6人回复
- 【求助】热处理与冷加工的顺序问题
已经有16人回复
- 863计划现代农业领域“食品非热加工技术与设备” 课题在京通过验收
已经有4人回复
3楼2012-02-25 19:58:07
4楼2013-06-22 09:05:24
2楼2012-01-10 14:45:12
5楼2013-09-22 17:31:41
6楼2014-06-27 17:07:26
7楼2014-08-25 15:29:35
8楼2015-01-14 16:15:09
9楼2015-03-25 14:29:18
10楼2015-05-17 20:53:46
11楼2015-06-12 14:28:21
12楼2015-08-31 10:52:05
13楼2015-08-31 10:53:19
14楼2015-08-31 10:54:09
|
楼主，您好，请问代码最开始的Z值和XY是怎么对应的？ X=[600,660,730,775,860,980] %X轴,温度 Y=[-3,-2,-1,0,1,1.70] %Y轴,log应变速率 Z=[108.3,78.5,57.1,36.1,24.3,10.9, 152.4,112.9,71.7,48.8,35.4,11.4, 184.9,154.4,112.9,77.8,52.6,15.2, 239.4,215.4,152.0,107.6,67.6,34.7, 269.8,249.1,221.7,178.3,110.1,61.8, 280.5,264.1,247.1,211.8,112.6,89.7] |
15楼2015-12-15 11:13:11
16楼2015-12-15 15:28:20
17楼2015-12-15 15:28:51
★★★★★ 五星级,优秀推荐
|
请教一下大神，在Zlogi=interp2(X,Y,Zlog,Xi,Ylogi','spline') %用spline方法插值(log应力)矩阵之后，出现 Error using griddedInterpolant The grid vectors do not define a grid of points that match the given values. Error in interp2/makegriddedinterp (line 214) F = griddedInterpolant(varargin{:}); Error in interp2 (line 127) F = makegriddedinterp({X, Y}, V, method,extrap); 我该怎么处理？ |
18楼2016-03-28 11:05:04
19楼2016-03-28 14:16:50
20楼2017-09-13 10:32:25
21楼2018-05-24 15:26:49
22楼2018-05-25 10:13:48
简单回复
向宽处行23楼
2018-09-28 10:42
回复
五星好评 
5