上机实验三：Matlab的应用.doc

计算机在材料科学中的应用  上机实验三：Matlab的应用 例、如图所示位一个2*2的矩形区域，其一边（x=0）给定温度，该温度T是y坐标的函数（T=500*y），另一边（x=2）为对流换热的传热边，暴露在温度为20℃的大气中。该表面与空气之间的对流传热系数k=20w/（m2·K）λ=5w/（m·K）9个格点的温度分布。 采用Matlab软件计算其温度场的分布情况。 比较分析上述两种方法的计算结果。 7 8 9 4 5 6 1 2 3 解：（1）采用有限差分法求解9个格点的温度分布： 7’ 8’ 9’ 7 8 9 12 4 5 6 11 1 1’ 2 2’ 3 3’ 10 1. ∵ 上下两条绝热边 ∴ T1=T1’ , T2=T2’ , T3=T3’ , T7=T7’ , T8=T8’ , T9=T9’ 左边那条边从下至上温度依次为0 , 250 , 500 , 750 , 100（K） 对1-9九个点列差分方程为： 250+T1+T2+T4==4*T1, T1+T2+T5+T3==4*T2, T2+T10+T6+T3==4*T3, 500+T1+T7+T5==4*T4, T4+T2+T6+T8==4*T5, T3+T9+T5+T11==4*T6, 750+T4+T7+T8==4*T7, T7+T8+T5+T9==4*T8, T6+T8+T9+T12==4*T9 对10、11、12三个点列差分方程： -10*(T3-T10)==20*(T10-293), -10*(T6-T11)==20*(T11-293), -10*(T9-T12)==20*(T12-293) 5. 综上12个方程利用MMA求解9个格点的温度： T1=70244/203=346.03 T2=70488/203=347.23 T3=63732/203=313.95 T4=3086/7=440.86 T5=2672/7=381.71 T6=2258/7=322.57 T7=108744/203=535.68 T8=84488/203=416.20 T9=67232/203=331.20. （K） （2）采用Matlab软件计算其温度场的分布情况： 步骤： 1.调制选项：Options - Application - Heat Transfer； 绘制2*2的矩形；  定义边界条件； x=0这条边为给定函数边：  y=0以及y=2这两条边为绝热边： （3）x=2这条边为对流传热边：  4.选择PDE模式，根据公式求解；  4.调制选项： Solve - Solve PDE； 5.调制选项： Mesh - Mesh Mode； 调制选项：Plot - Plot Solution； 7.调制选项：I. Mesh - Export Mesh II. Solve - Export Solution III. Mesh - Show Node Labels; 由输出的pet表格找出各点对应的数字、坐标以及温度;  （3）比较分析上述两种方法的计算结果： 1. 利用有限差分法以及MMA求得的温度为 T1=70244/203=346.03 T2=70488/203=347.23 T3=63732/203=313.95 T4=3086/7=440.86 T5=2672/7=381.71 T6=2258/7=322.57 T7=108744/203=535.68 T8=84488/203=416.20 T9=67232/203=331.20. 2. 利用Matlab作图计算出的温度为 总结：对比可以看出，9个点中大致有6个点还是比较对应。 利用差分法以及MMA求解时要注意取值得当，避免误差随步数增加而累积造成解得不精确； 利用Matlab软件求解温度场时要注意边界条件和方程的正确选取以及精确赋值，避免一步错，步步错。 绝热边 750K 500K 对流边 250K 绝热边 中心点差分公式 对流传