上机实验三：matlab应用.docVIP

计算机在材料科学中的应用  上机实验三：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个格点的温度分布： 绝热边 绝热边 750K 7’ 750K 8’ 9’ 500K 7 500K 8 9 对流边 12 对流边 250K 4 250K 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.调制选项： = 1 \* ROMAN \* MERGEFORMAT I. Mesh - Export Mesh = 2 \* ROMAN \* MERGEFORMAT II. Solve - Export Solution = 3 \* ROMAN \* MERGEFORMAT 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作图计算出的温度