弹性力学用差分法和变分解平面问题.pptVIP

（6）不规则边界内节点、虚节点? 值 h h h 1 0 9 5 6 B ? h 边界外 y x q 3qh 3qh 19 21 20 16 17 18 22 23 24 25 26 1 2 3 4 7 10 13 5 8 11 14 6 9 12 15 A D C B I H G F E J K L M §5-3 应力函数的差分解的实例 1. 问题 设一正方形的混凝土深梁（边长 6h），上边界受有均布压力 q ，并下角点处的两反力维持平衡。试由应力函数的差分解法，求各节点的应力分量。 2. 求解 由于对称性，如图建立坐标系，并取其一半分析。求解过程： （1） 适当划分差分网格、编节点号； （2） 选取基点A，并计算边界节点的 ? 及其导数值； 计算公式： y x q 3qh 3qh 19 21 20 16 17 18 22 23 24 25 26 1 2 3 4 7 10 13 5 8 11 14 6 9 12 15 A D C B I H G F E J K L M 4.5qh2 4.0qh2 2.5qh2 0 0 0 0 0 ? 0 0 0 – – – 0 0 – – – – 3qh – – 0 M L K J E、F、G、H、I D B、C A 结点 （3） 计算边界外一行虚节点的 ? 值； §5-1 差分公式的推导 主 要 内 容 §5-2 应力函数的差分解 §5-3 应力函数的差分解的实例 §5-4 位移变分方程 第五章 用差分法和变分法解平面问题 §5-5 位移变分法 §5-6 位移变分法的例题 §5-0 弹性力学的数值计算方法简介 工程问题 （力学、物理等） 建立一组 基本方程 控制微分方程 定值条件 常微分方程 偏微分方程 位移边界条件 力的边界条件 初始条件 求解 精确解 近似解 （数值解） （均质、边界条件简单） （1）有限差分法 （2）等效积分法（包括变分法） （3）有限单元法 （4）边界单元法 …… （1）有限差分法（FDM） 要点： 差分 微分； 代替 差分方程 代替 微分方程。 （代数方程） y x 0 1 3 h h 优点： y x 0 1 3 h h 收敛性好、程序设计简单、 非线性适应好。 代表性软件：FLAC 缺点： 当边界几何形状复杂时，解的精度受到限制。 （2）等效积分法 控制微分方程 边值条件 建立等效的 积分方程 假设未知函数 整个区域内 近似求解 （a）加权余量法（加权残值法） （配点法、子域法、最小二乘法、力矩法、Galerkin法、等） （b）变分法 当原问题存在某个泛函时，则原问题等价于求该泛函的驻值。如：Ritz 法等。 特点： 在整个区域内，假设未知函数。 适用于边界几何形状简单的情形。 x y （3）有限单元法（FEM） —— 加权余量法、变分法的推广。 基本思想： 整个区域 分成若干个单元 区域离散 假设未知函数 在单元上 由变分原理等求出单元结点上值 （近似解） 主要有限元软件： SAP， ADINA NASTRAN、 ANSYS 、 ABAQUS、ASKA、SAFE、MARC等 —— 早期的软件 要 点: 将微分方程转变成差分方程。 基本思想: 将基本方程和边界条件（一般为微分方程）近似地用改用差分方程（线性代数方程）表示，把解微分方程的问题变成求代数方程的问题。 §5-1 差分公式的推导 1. 中心差分公式 设函数： 为弹性体 内的某个函数（应力分量、位移分量、应力函数 ? 、温度等）。 y x 在弹性体上用相隔等间距 h 且平行于坐标轴的两组平行线组成网格，称为差分网格。网格线的交点称为节点（结点）。 0 1 3 h 则函数 f = f (x，y) 在平行于 x 轴的网格线上，如节点：3-0-1 上，它只随 x 而变化。 考察结点 0 处，函数 f = f (x，y) 的变化，可展开成 Taylor 级数： （a） h 若略去三次幂以上各项，式（a）变为： （b） 节点3及1的 x 坐标： 将其代入式（b），有： （c） （d） 联立求解，得： （5-1） （5-2） y x 0 1 3 h h 2 4 y x 0 1 3 h h （5-1） （5-2） 同理，在网格线4-0-2上取 （e） 类似于，x方向的讨论，有 （5-3） （5-4） 式（7-1）?（7-4）称为基本差分公式。 5 6 7 8 混合二阶导数的差分公式： （5-5） 进一步可导出四阶偏导数的差分公式： y x 0 1 3 2 4 h h 进一步可导出四阶偏导数的差分公式： 9 10 11 12 6 y x 0 1 3 2 4