燕山大学Matlab三级项目.doc_第1页
燕山大学Matlab三级项目.doc_第2页
燕山大学Matlab三级项目.doc_第3页
燕山大学Matlab三级项目.doc_第4页
燕山大学Matlab三级项目.doc_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

基于Matlab的导热问题的数值解法 班级：13级机械装备2班 指导教师：任素波 引言传热学在工程中应用非常广泛，特别是冶炼、轧钢行业，处处有传热学的身影。虽然传热过程是容易做数学处理的一种热量传递方式，但其往往涉及到常微分方程、偏微分方程、线性及非线性方程组的求解，对于数学知识掌握不深的专业人员而言，解题有一定难度。高级计算机语言的出现给数值计算带来了希望，但用这些计算机语言编制程序，需要对所用的语言语法和相关算法有相当的熟练度。Matlab软件是一种功能强、效率高、便于进行科学和工程计算的交互式软件包。它包含了一般数值计算、矩阵计算等应用程序，易学易用，大大降低了对使用者的数学基础和计算机语言知识的要求，而且编程效率和计算效率极高。因此，对于数学模型复杂，涉及较多运算方法传热过程，Matlab是一种非常重要的辅助工具。1. 实际问题 有一个边长为1m的正方形物体如下图所示，其左侧面恒温为100，顶面恒温为500，其余两侧均暴露在对流环境中，环境温度为100。已知物体的导热系数为，求出各节点的温度值。 2. 问题分析 求解导热问题，实际上就是对导热微分方程在给定条件下进行求解，找出具体的函数关系式，从而得出导热过程的数学表达式。本案例是一个很常见的导热问题，我们可以使用温度节点离散化线性方程组，并用矩阵形式来表示。3. 数学模型 对上述9x9阶矩阵进行求解，相当于解九个方程，最终求出九个未知数，其难度以及计算量对于工程人员来说是一个很大的挑战。但如果使用Matlab进行编程计算，则会大大节省计算时间，减低错误的出现几率，对于工程运算来说如虎添翼。4. Matlab编程此方程组属于大型稀疏线性方程组，求解时广泛使用LU分解法，求解速度快，占用计算机内存小，便于储存。下面为求解过程：functionLU a1=-41010;1-4101;02-4.6700;100-41;0101-4; a2=0000;0000;1000;0100;1010; a3=00102;00020;00002;00000; a4=-4.67001;0-4.6710;01-4.671;101-2.67; a=a1,a2;a3,a4; b=-600,-500,-567,-1000,-66.7,-167,-66.7,-66.7; L,U=lu（a）; t=U（Lb） 执行，可以得到九个节点温度依次为： 5. 结果分析 使用Matlab对矩阵进行数值运算，其结果和人工计算相同。但与人工计算相比，其速度和准确率都大大提高。而且Matlab的操作工程并不复杂，仅需要编写几行程序就可以进行运算。因此，对于数学模型较为复杂的工程求解问题，都可以考虑采用Matla
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论