第三节微分方程的数值计算原理机械加工.pptVIP

(3-40) 由式（3-27），对于六点差分格式（即C-N格式）有： (3-41) 将（3-38）式、（3-39）式、(3-40)式和(3-41)式描绘成函数图形，如图3-3所示。可以看出，当f值很小时，各种差分格式很接近于精确解；但当f值逐渐增大时，不同差分格式的情况就大相径庭。在f较大，如f=1/6时，C-N格式与精确解最为接近；当f继续增大时，C-N格式与完全隐式格式比较接近精确解。 需要指出的是，上述分析只能说明由于对 项采取不同的差分近似而导致的不同差分方程的相对准确性，并不能说明一个具体差分格式的绝对精度。事实上，差分格式的准确性除取决于差分方程本身外，还可能受到边界条件差分化的影响。 图3-3 不同差分格式的精度和稳定性 §3.4直接差分法与边界条件的差分化 §3.4.1直接差分法 §3.4.2边界条件的差分化 §3.4.1直接差分法 除了基于导热微分方程用差分代替微分来建立差分方程外，还可将热传导定律直接应用于经网格剖分后的单元体，通过能量平衡列写差分方程。这种方法就是所谓的直接差分法，有时又被称为单元热平衡法。如图3-4所示的二维节点（i，j），其能量平衡为： 当λx=λy，△x=△y时，整理得： 这是二维显式差分方程，与用差分代替微分由微分方程得到的结果相同，同样用这种方法也可得到各种