原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
用有限差分法和Matlab计算二维热加工温度场分析最新版本
用有限差分法和Matlab计算二维热加工温度场分析
Eg1 薄板焊接过程温度场分析。
取焊件的一半为模型进行离散化，起始点为o点，以后以v速度沿x轴运动。
根据题意为二维不稳态导热，二维不稳态导热方程为：
图 二维焊接离散化
题目可化为以下微分方程组：(以y轴正方向为上，x轴正方向为右)
需要的参数（均已cm,cal,g为单位，所以不必换算）：
参数 数值 备注 ρ 7.82 g/cm3 密度 v 0.4 cm/s 焊接速度 h 1cm 板厚度 Qm 4000 cal/cm3 热源分布密度 β 0.0008cal/cm2·s·℃ 换热系数 Te，T0 20℃ 周边介质温度，初始温度 k 0.1 cal/cm·s·℃ 导热系数
用PDE Tool解题步骤如下：
0. pdetool
区域设置
单击工具，在窗口拉出一个矩形，双击矩形区域，在Object Dialog对话框输入Left为0，Bottom为0，Width为2,Height为2。
与默认的坐标相比，图形小的看不见，所以要调整坐标显示比例。方法：选择Options-Axes Limits,把X，Y轴的自动选项打开。
设置Options-Application为Heat Transfer （设置程序应用热传输模型）
设置边界条件
单击 ，使边界变红色，然后分别双击每段边界，打开Boundary Conditions对话框，设置边界条件（根据边界条件）。所有的边界都为Neumann条件。输入值见下表：
边界 g (Heat Flux) q (Heat Transfer Coefficient) 左边界 0 0 右边界 0.0008*20 0.0008 下边界（x轴） -0.0008*20 -0.0008 上边界 0.0008*20 0.0008
设置方程类型
单击，打开PDE Specification对话框，设置方程类型为Parabolic（抛物型），
C（比热）为0.16，a（导热系数）为0.1，d(密度）为7.82，f（热源）为4000*exp(-3*(x.^2+(y-0.4*t).^2)/0.49)。
网格划分
单击，或者加密网格，单击。
初值和误差的设置
单击Solve菜单中Parameters…选项，打开Solve Parameters对话框，输入Time为0:0.5:5，u(t0)为20,其他不变。
解方程
单击，开始解方程。
整理数据
单击Mesh-Export Mesh…输出 p e t 的数值，单击Solve-Export Solution…输出u.
回到Matlab主窗口执行下面两条命令：
u1=[p,u(:,7)] %将节点坐标和其在3s时的温度组成新矩阵
u2=sortrows(u1,3) %将u1按温度值大小升序排列。
u1=[p,u(:,4)] %将节点坐标和其在1.5s时的温度组成新矩阵
u2=sortrows(u1,3) %将u1按温度值大小升序排列。
温度场分布：
点plot parameters 勾选相关项目，右下角的时间选为1.5S,然后点plot 即可。
1.5s时
1.5s时的三维图
3s 时的温度场：
3s时的三维图：
4
文档评论（0）