0708第一学期计算机在材料科学与工程中的应用答案.docVIP

07/08第一学期计算机在材料科学与工程中的应用答案 一、 写出四阶龙格库塔的计算公式，并画出计算机计算流程图 四阶龙格库塔的计算公式为： 式中h为步长。 二、将y’=a’ebx,并写出求解系数的计算公式。 y=lny’ , a=lna’ b=b, x=x; y=a+bx 解得： 相关系数r 式中 ; r0 拟合曲线斜率为正 r0 斜率为负 r=0 则x和y无线性关系 三、Ax=B，写出追赶法求解过程 设存在克劳特分解 A=LU 其中 L和U分别为: 将其代入分解式 = = 对比等式两边的元素可得 存在克劳特分解? 设: Ax=f； 其中A=[ai, bi, Ci]n1 右端向量 f=(f1,f2,…,fn)T Ly=f ； Ux=y； y1=f1/l1; yi=(fi-mi yi-1) I=2,3,…, n Xn=yn, xi=yi-uixi+1 I=n-1,n-2,…,1 四、建模： 由于对称性，取厚度的一半处理；h=0.006/6=0.001; 时间按1秒，初始条件，T[i]=293K, 右边界条件：Tn+11=293+（700-20）/（4G60G60）G?τ 内部传热为： 左边界条件： 离散： 内部节点离散： 将右边界条件代入： 整理成： （1） 内部离散整理为： （2） 右边界条件离散： 写成 代入内部节点： （3） 写成矩阵：