弹性力学简化 平面问题的差分解.pptxVIP

弹性力学简化 平面问题的差分解第1页/共171页 §7-1 差分公式的推导主 要 内 容 §7-2 稳定温度场的差分解§7-3 不稳定温度场的差分解§7-4 应力函数的差分解§7-5 应力函数的差分解的实例§7-6 温度应力问题的应力函数差分解§7-7 位移差分解§7-8 位移差分解实例§7-9 多连体问题的位移差分解§7-10 温度应力问题的位移差分解第2页/共171页 §7-0 弹性力学的数值计算方法简介工程问题（力学、物理等）建立一组基本方程控制微分方程定值条件常微分方程偏微分方程位移边界条件力的边界条件初始条件求解精确解近似解（数值解）（均质、边界条件简单）（1）有限差分法（2）等效积分法（包括变分法）（3）有限单元法（4）边界单元法……（1）有限差分法（FDM）要点：差分微分；代替差分方程代替微分方程。（代数方程）yx013hh第3页/共171页 优点：yx013hh收敛性好、程序设计简单、非线性适应好。代表性软件：FLAC缺点：当边界几何形状复杂时，解的精度受到限制。（2）等效积分法控制微分方程边值条件建立等效的积分方程假设未知函数整个区域内近似求解（a）加权余量法（加权残值法）（配点法、子域法、最小二乘法、力矩法、Galerkin法、等）（b）变分法 当原问题存在某个泛函时，则原问题等价于求该泛函的驻值。如：Ritz 法等。特点： 在整个区域内，假设未知函数。适用于边界几何形状简单的情形。xy第4页/共171页 （3）有限单元法（FEM）—— 加权余量法、变分法的推广。 基本思想：整个区域分成若干个单元区域离散假设未知函数在单元上由变分原理等求出单元结点上值（近似解）主要有限元软件：SAP， ADINANASTRAN、 ANSYS 、 ABAQUS、ASKA、SAFE、MARC等—— 早期的软件第5页/共171页 1. 中心差分公式§7-1 差分公式的推导设函数：为弹性体内的某个函数（应力分量、位移分量、应力函数 ? 、温度等）。yx 在弹性体上用相隔等间距 h 且平行于坐标轴的两组平行线组成网格，称为差分网格。网格线的交点称为节点（结点）。013h 则函数 f = f (x，y) 在平行于 x 轴的网格线上，如节点：3-0-1 上，它只随 x 而变化。 考察结点 0 处，函数 f = f (x，y) 的变化，可展开成 Taylor 级数：（a）h第6页/共171页 若略去三次幂以上各项，式（a）变为：（b）节点3及1的 x 坐标：将其代入式（b），有：（c）（d）联立求解，得：（7-1）（7-2）yx013hh第7页/共171页 24yx013hh（7-1）（7-2）同理，在网格线4-0-2上取（e）类似于，x方向的讨论，有（7-3）（7-4）式（7-1）?（7-4）称为基本差分公式。第8页/共171页 5678混合二阶导数的差分公式：（7-5）进一步可导出四阶偏导数的差分公式：yx01324hh第9页/共171页 进一步可导出四阶偏导数的差分公式（7-6）以两侧节点处的函数值表示中间节点处的一阶导数值，称为中点导数值，这种差分公式称为中心差分公式。说明：第10页/共171页 2. 端点差分公式向前差分公式yx0132456789101112h把导数：用函数值：f0 f1 f9表示；把导数：用函数值：f0 f2 f10表示。由：（b）得到：第11页/共171页 yx0132456789101112h（7-7）（e）同理，对 y 方向，有：由此解得：（7-9）式（7-7）、（7-9）称为前差公式。第12页/共171页 向后差分公式yx0132456789101112h把导数：用函数值：f0 f3 f11表示；用函数值：f0 f4 f12表示。由：把导数：（b）得到：由此解得：（7-8）第13页/共171页 yx0132456789101112h同理，对 y 方向，有：（e）可解得：（7-10）（7-8）式 (7-8) 、(7-10) 称为后差分公式第14页/共171页 与中心差分公式类似，由式（7-7）~（7-10）可推出高阶导数的差分公式。（1）中心差分（导数）公式与 端点差分（导数）公式比较，前者的精度较高。所以尽可能应用中心差分公式。说明：（2）在前面差分公式的推导中，应用了近似式：（b）略去了三次幂以上的各项，其实质：在（x – x0）的区间上，将f（x，y）沿 x 方向用抛物线函数代替。所以，式（7-7）~（7-10）称为抛物线差分公式。（3）若在差分公式的推导中，应用线性近似关