基于pdetool热传导数值计算.pdfVIP

　第 24 卷第 2 期 　　　　　佳 木 斯 大 学 学 报 ( 自 然 科 学 版 ) 　　 Vol. 24 No. 2 　2006 　年 04 月 　　　　　Journal of Jiamusi University (Natural Science Edition) 　　 Apr. 　2006 ( ) 文章编号 :1008 - 1402 2006 02 - 0270 - 03 ① 基于pdetool 的热传导数值计算 徐梓斌 , 　闵剑青 (浙江树人大学 ,浙江 杭州 310015) 摘　要 : 　分析解法无法解决复杂的几何形状或边界条件的导热问题 ,本文运用 pdetool 工具箱 ,对工 程中经常遇见的周期性变化边界、直肋及三维非稳定导热等问题进行了数值求解. 结果表明应用 pdetool 在解决此类问题是高效、直观的 ,它不仅能解二维导热问题 ,通过变化公式和设置合适的边界 ,可以解决一 维和某些三维导热问题 ,为工程实际问题的解决提供了有效的工具. 关键词 : 　导热 ;pdetool 工具箱 ;数值模拟 中图分类号: 　TK124 　　　文献标识码 : 　A 　　导热问题的分析解法包括直接积分法和分离 变量法 ,即使在解简单的导热问题时也是相当繁杂 的[1 ] ,因此应用分析解法有时会采用简化分析的方 法 ,这会带来一定的计算误差 ,对复杂的几何形状 或边界条件复杂的问题根本就无能为力. 而基于有 限元的MATLABpdetool 数值计算则是十分有效的 方法 ,pdetool 是一个用来研究与求解偏微分方程的 工具集合 ,它为二维空间上的偏微分方程的研究与 求解提供了一个强大而灵活的工作环境. pdetool 在传热中主要用于解如下两类偏微分 方程 : 1 　算例 ( ) α 椭圆型 PD E : - ·c u + ·u = f u 1. 1 　周期性变化边界条件下的导热问题 ( ) α 抛物线型 PD E : d t - ·c u + ·u = f 某地的地面年最高温度30 . 5 ℃,最低 - 3. 5 ℃, α 式中 , u 为求解变量 , t 为时间变量 , d , c , ,f 将温度的周期性波动看成简谐波 ,即地面平均温度 为常数或变量. T = 13. 5 ℃,地面年温度振幅 A = 17 ℃,年温度波 pdetool 定义了两类边界条件 : m ( ) 动频率 F = 124 3600 365 Hz ,土壤的热扩散 Dirichlet 边界条件 h ·u = r - 6 2 率 a = 0 . 617 ×10 m s ,设过余温度 u = T - T ; ( Δ)