摘要: 下面的有限差分程序计算2D稳态热传导方程。其模型和文章有限元计算2D稳态热传导方程的模型一样。从后面也可以看到,其结果也和该文章的结果一样。即有限差分和有限元法,解决同一问题,得到了相同的解。程序见下面:-----------------------------------------------------------------------------------% Solution of 2D steady heat conduction equation(Poisson equation): d2FI/dx2+d2FI/dy2=f % with finite differences o 阅读全文
posted @ 2012-11-23 15:22 heaventian 阅读(734) 评论(0) 推荐(0) 编辑
摘要: 本文程序,使用有限元计算2D稳态热传导方程。本程序只使用2D三角形单元求解,没有四边形单元求解。不过,也非常容易就可做到。主程序为:Heat_Conduction_steady.m需要用到的子程序有:coord3.m 用来生成节点坐标,节点坐标与单元的链接。边界条件等等。生成文件coordinates.dat 包含两列,第一列节点x坐标,第二列,节点y坐标elements3.dat 包含三列:分别是各个单元的i,j,k节点的总体编号。dirichlet.dat 包含两列:分别是第一类边界条件的两端点节点的总体编号。另外,coord4.m与coord3相同,但是,生成的是四边形单元。即:... 阅读全文
posted @ 2012-11-23 14:59 heaventian 阅读(3189) 评论(0) 推荐(0) 编辑