[2018年最新整理]Chapter3有限元分析的数学求解原理.pptVIP

一般说来，求解方程的途径有两大类： 1）直接针对原始方程进行求解（解析法、半逆解法以及有限差分法） 2）间接针对原始方程进行求解（加权残值法、虚功原理、最小势能原理以及变分方法等） 直接解法——解析法： 解析法从力平衡关系、几何关系以及物理关系出发，推导出一个或一组关于应力或者关于应变、有时是同时含有应力、应变的微分方程或偏微分方程，通过求解微分方程，解出应力、应变和变形量。工程中，常采用的解析方法有材料力学中对杆件的分析，弹性力学中平面问题的求解，板壳理论等。 解析法的很多基本理论是建立在一些简化的假设基础之上的，经过大量的工程实践，被证明能很好的符合构件实际工作情况，已成为成熟的理论。解析法得到的结果是未知量（应力、应变等）的函数解，可直接得到结构中任意点的精确解。解析法在分析理论问题以及一些工程问题时起着重要作用。但是解析法在应用到一些形状复杂或应力分布复杂的结构时，往往由于数学上的问题而显得无能为力，因而使解析法在应力分析中的应用受到限制。 直接解法——有限差分法： 有限差分法：微分方程和积分微分方程数值解的方法。其基本思想是： 有限差分方法(finite difference method )是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 有限差分格式 格式精度：一阶格式、二阶格式和高阶格式。 差分的空间形式：中心格式 时间因子：显格式、隐格式、显隐交替格式等。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 间接解法——加权残值法： 是一种应用广泛的求解微分方程的方法.该方法先假定一族带有待定参数的定义在全域上的近似函数,该近似解不能精确满足微分方程和边界条件,即存在残差.在加权平均的意义下消除残差,就得到加权残值法的方程.由于试函数定义在全域上,所得方程的系数矩阵一般为满阵.选取不同的权函数,可得到不同的加权参量法. 间接解法——最小势能原理： 间接解法——变分原理： 把一个物理学问题用变分法化为求泛函极值（或驻值）的问题，后者就称为该物理问题 的变分原理。如果建立了一个新的变分原理，它解除了原有的某问题变分原理的某些约束条件，就称为该问题的广义变分原理；如果解除了所有的约束条件，就称为无条件广义变分原理，或称为完全的广义变分原理。 在当代，变分原理已成为有限元法的理论基础，而广义变分原理已成为混合和杂交有限元的理论基础。在实际应用中，通常很少能求出精确的解析解，因此大多采用近似计算方法。 工程问题无论是几何形状、受力方式还是材料特性都是前变万化的，因此一种求解方法是否有优势，其判断标准应该是 具有良好的规范性（不需要太多的经验和个人技巧） 具有良好的适应性（可以处理任何复杂的工程实际） 具有良好的可靠性（计算结果收敛稳定，精度高） 具有良好的求解可行性（计算工作量） 本章主要内容 3.1简明问题的解析求解 3.2弹性力学问题近似求解的加权残值法 3.3弹性力学近似求解的虚功原理、最小势能原理及其变分原理 3.4各种求解方法的特点及比较 3.1简明问题的解析求解——一维拉压杆问题 基本变量：ux(x),σx(x),εx(x) 基本方程： 对三大方程直接进行求解得 讨论1 若用材料力学的经验方法求解，则需先作平面假设，即假设应力为均匀分布 σx=P/A 由广义胡克定律得 εx=P/EA 右端的伸长量为 Δu=εxl=Pl/EA 讨论2 应变能 动能 势能 3.1简明问题的解析求解——平面梁问题 简支梁的特征： 梁为细长梁，因此可以只用x坐标来刻画； 主要变形是垂直于x的挠度，可以只用挠度来描述位移场。 针对这两个特征，可以做出以下两个假定： 直法线假定 小变形和平面假设 直法线假定：一垂直中面的直线(称为法线)，变形时不伸缩，并且仍为弹性曲面的法线。 平截面假定：平截面假定是材料力学中最基本的假定之一。这个假定认为所有与杆轴线垂直的截面在杆件变形后仍保持为