- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
用PDE Toolbox求解三类典型方程求解抛物型方程
实验二 用PDE Toolbox求解三类典型方程
----求解抛物型方程
一、实验目的
本次实验在掌握PDE Toolbox的功能和PDE图形用户界面的基础上,求解抛物型方程。
二、实验内容
1、解抛物型方程的一个例子
热传导方程:金属板的导热问题
考虑一个带有矩形孔的金属板上的热传导问题。板的左边保持在100℃,板的右边热量从板向环境空气定常流动,其他边及内孔边界保持绝缘。初始时板的温度为0℃,于是概括为如下定解问题:
的外边界定点坐标为(-0.5,-0.8),(0.5,-0.8),(0.5,0.8),(-0.5,0.8),内边界顶点坐标为(-00.5,-0.4),(00.5,-0.4),(00.5,0.4),(-00.5,0.4)。
使用GUI求解这一问题。在PDE Toolbox窗口的工具栏中选择Generic Scalar.
第一步:区域设置 单击工具,在窗口拖拉出一个矩形,双击矩形区域,在对话框中输入Left为-0.5,Bottom为-0.8,Width为1,Height为1.6,单击OK按钮,显示矩形区域R1。用同样方法做内孔R2,只要设置Left为-00.5,Bottom为-0.4,Width为0.1,Height为0.8即可。然后在栏中键入R1-R2。
第二步:设置边界条件 单击工具,图形边界变红,逐段双击边界,打开Boundary Condition对话框,输入边界条件。对于同一类型的边界,可按shift键,将多个边界同时选择,统一设置边界条件。本题在左边界上选择Dirichlet条件,输入h 为1,r为100;右边界上,选择Neumann条件,输入g为-1,q为0;其他边界上,选择Neumann条件,输入g为0,q为0 。
第三步:设置方程类型 单击PDE菜单中PDE Specification…选项,打开PDE Specification对话框,选择方程类型。本题单击Parabolic(抛物型),输入d为1,c为1,a为0,f为0,然后单击OK按钮。
第四步:网格剖分 单击工具,或者单击Mesh菜单中Initialize Mesh选项,可进行初始网格剖分,如果需要网格加密,再单击,或者单击Mesh菜单中Refine Mesh选项。
第五步:初值和误差设置 单击Solve菜单中Parameters…选项,打开Solve Parameters…对话框,输入Time为0:5,u(t0)为0 ,Relative tolerance为0.01,Absolute tolerance为0.001,然后单击OK按钮。
第六步:解的图形 单击Plot菜单中Parameters…选项,出现Plot Selection对话框,从中可以选择Color, Contour, Arrows, Deformed mesh, Height(3-D plot),还可以设置等值线的数目等。本例中选择Color, Contour和Arrows三项,然后单击Plot按钮,窗口显示出时的解的彩色图形。
第七步:数值解的输出 单击Solve菜单中Export Solution选项,在打开的Export对话框中输入u,单击OK按钮确定。再在MATLAB命令窗口键入u,按回车键确定,即显示按节点编号排列的数值解
u
u =
0 79.6291 95.3774 98.0380 98.5706 98.6395
0 79.6292 95.3774 98.0380 98.5705 98.6395
100.0000 100.0000 100.0000 100.0000 100.0000 100.0000
0 83.3577 96.3858 98.5869 99.0275 99.0845
0 83.3653 96.3874 98.5875 99.0279 99.0849
0 89.8188 97.8111 99.1613 99.4315 99.4664
0 89.8118 97.8096 99.1607 99.4310 99.4660
100.0000 100.0000 100.0000 100.0000 100.0000 100.0000
0 79.3324 95.3139 98.0139 98.5544 98.6243
0 78.5592 95.1462 97.9486 9
您可能关注的文档
- 生物图表曲线类题目的解题技巧.ppt
- 生物安全柜1洁净过程验证.doc
- 生物工程设备复习资料.doc
- 生物圈是最大的生态系统(新教材).ppt
- 生物工程设备考试题目.docx
- 生物技术 发酵工程.ppt
- 生物材料检验各重点.doc
- 生物技术 酶工程.ppt
- 生物学导论biology12.ppt
- 生物污泥的减量化研究.ppt
- 分析let s单元56ago2卷纸zheng unit56.pdf
- 塑胶材料其它分类原料pa9t 12.pdf
- md16x16数字媒体切换器设备.pdf
- 者参考项目发起人学科类型单位序列承包商修订页代码顺序典型.pdf
- 届世界天然气大会阿姆斯特丹2006add10288.pdf
- 期测试记录表每周weekly g1g6 journeys tests level 6 lesson26.pdf
- modernize-whitepaper现代化您应用程序白皮书.pdf
- anybackup产品典型案例分析.pdf
- 约克金融工程课程tfeslide32.pdf
- 广州市妇女儿童医疗中心历份教学药历01tjy.pdf
文档评论(0)