下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、温度场模拟matlab代码:clear,clc,clfL1=8;L2=8;N=9;M=9;% 边长为8cm的正方形划分为8*8的格子T0=500;Tw=100; % 初始和稳态温度a=0.05; % 导温系数tmax=600;dt=0.2; % 时间限10min和时间步长0.2sdx=L1/(M-1);dy=L2/(N-1);M1=a*dt/(dx2);M2=a*dt/(dy2);T=T0*ones(M,N);T1=T0*ones(M,N);t=0;l=0;k=0;Tc=zeros(1,600);% 中心点温度,每一秒采集一个点for i=1:9 for j=1:9 if(i=1|i=9|j=
2、1|j=9) T(i,j)=Tw;% 边界点温度为100 else T(i,j)=T0; end endendif(2*M1+2*M2=1) % 判断是否满足稳定性条件 while(ttmax+dt) t=t+dt; k=k+1; for i=2:8 for j=2:8 T1(i,j)=M1*(T(i-1,j)+T(i+1,j)+M2*(T(i,j-1)+T(i,j+1)+(1-2*M1-2*M2)*T(i,j); end end for i=2:8 for j=2:8 T(i,j)=T1(i,j); end end if(k=5) l=l+1; Tc(l)=T(5,5); k=0; end
3、end i=1:9;j=1:9; x,y=meshgrid(i); figure(1); subplot(1,2,1); mesh(x,y,T(i,j)% 画出10min后的温度场 axis tight; xlabel(x,FontSize,14);ylabel(y,FontSize,14);zlabel(T/,FontSize,14) title(1min后二维温度场模拟图,FontSize,18) subplot(1,2,2); C,H=contour(x,y,T(i,j); clabel(C,H);axis square; xlabel(x,FontSize,14);ylabel(y,FontSize,14); title(1min后模拟等温线图,FontSize,18) figure(2); xx=1:600; plot(xx,Tc,k-,linewidth,2) xlabel(时间/s,FontSize,14);ylabel(温度/,FontSize,14);title(中心点的冷却曲线,FontSi
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年新闻采辑项目立项申请报告模板
- 2024年工控电源项目申请报告
- 2024年减肥项目申请报告范本
- 2024年医疗健康电子产品项目规划申请报告范本
- 2024年光电子器件及激光器件项目规划申请报告
- 2024年工业VOC治理系统项目申请报告范样
- 产业聚集形成的原因和影响的研究
- 健康预期寿命提高如何促进经济增长基于跨国宏观数据的实证研究
- 【可行性报告】2023年止痛药项目可行性研究分析报告
- 【可行性报告】2023年潮玩行业项目可行性分析报告
- 第13课+学会依法维权(课时1)【中职专用】中职思想政治《职业道德与法治》高效课堂(高教版2023·基础模块)
- 印刷排版设计与操作
- 电商运营提成方案
- 2024年宁夏宁东融资担保有限公司招聘笔试参考题库含答案解析
- 深圳的餐饮行业现状分析报告
- 人工智能赋能教育行业培训
- 学生心理健康小品剧本
- 泌尿编码课件
- 服装陈列课件
- 2023气管插管意外拔管的不良事件分析及改进措施
- 《高中数学余弦定理》课件
评论
0/150
提交评论