有限差分法学案.pptVIP

第十章 有限差分法 弹性力学的经典解法存在一定的局限性，当弹性体的边界条件和受载情况复杂一点，往往无法求得偏微分方程的边值问题的解析解。因此，各种数值解法便具有重要的实际意义。;第一节 差分方程; 工程中常用得数值解法有有限单元法和差分法。; 我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，Δx=Δy=h，如图。; 我们将只考虑离开结点０充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为：;联立(c),(d),解得差分公式： ; 以上(１-１)~(１-４)是基本差分公式,从而可导出其它的差分公式如下：; 差分公式(１-１)及(１-３)是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。;第二节 应力函数的差分解; 一旦求得弹性体全部节点??φ值后，就可按应力分量差分公式（对节点0）算得弹性体各节点的应力。; 可见，用差分法解平面问题，共有两大任务： 一、建立差分方程 将（1-6）代入双调和方程; 二、联立求解这些线性代数方程，就能求得各内结点处的值。;代入上式，即得： ;关于边界上任一点处; 将此式亦从A点到B点沿s进行积分，就得到边界上任一点B处