有限元与数值方法-讲稿2.pptVIP

1有限元与数值方法第二讲授课教师：刘书田TelEmail：教室：研究楼351时间：2013年3月22日：8：00-9:35;9:45-10：30

2固体力学控制方程的根本形式〔弹性力学为例〕一般方程根本未知数〔15〕：应力〔6〕；应变〔6〕位移〔3〕根本方程〔15〕：平衡方程〔3〕，本构方程〔6〕，几何方程〔6〕边界条件：应力边界条件；位移边界条件位移解法根本未知数：位移分量〔3个〕根本方程〔3〕：用位移表示的平衡方程〔3〕边界条件：应力边界条件；位移边界条件应力解法根本未知数〔6〕：应力〔6〕根本方程：平衡方程〔3〕，相容方程〔3+3〕边界条件：应力边界条件；位移边界条件

3根本控制方程：微分方程的边值问题求响应量平衡方程几何方程物理方程〔热弹性〕〔各项同性〕

4求解方法及相对应的控制方程： 力法：以应力为根本位移量 平衡方程：变形协调方程 位移法：以位移为根本未知数 平衡方程：位移表示的平衡方程对应原理：变分原理微分方程的积分形式，泛函变分与根本方程的对应。建立各种问题所对应的变分原理。总之：固体力学是建立固体变形规律所必需满足的规律以及数学模型。为各种求解策略提供理论根底

5求解策略物理模型的简化杆〔梁〕平面问题平面应力问题，薄板平面应力问题：〔很厚〕柱体板〔壳〕开展数值求解策略直接求解方法：差分方法积分形式的控制方程问题：有限元，加权余量等

6杆件(包括杆和梁)杆:只在轴向受力梁:考虑弯曲变形几何特点:结构两个方向的尺度相近,但远小于第三方向尺度;根本假定:垂直于杆件轴线的平断面变形后保持平断面;三维问题简化为一维;无分布力时，杆的变形用线性函数描述即可

7平面应力问题几何特点:结构两个方向的尺度相近,但远大于第三方向尺度;受力特点:只受到在平面内的力;平面外的应力分量为0;三维问题简化为二维xyz

8平面应变问题几何特点:结构两个方向的尺度相近,但第三方向尺度无穷大;受力特点:受到的力在第三方向是均匀的;平面外的应变分量为0;三维问题简化为二维xyz

9xyzq(x,y)几何特点:结构两个方向的尺度相近,但远大于第三方向尺度;受力特点:只受到在平面外的力;直法线假定;三维问题简化为二维平板弯曲问题

10第三章：数值方法的根本理论力学问题可表示成微分方程的边值问题。因此，力学问题的求解可参照微分方程的求解方法构造。根据问题的不同，可采用微分形式，也可采用积分形式。不同的描述方式，可构造不同的数值方法。计算机最容易处理的是数值计算、求解代数方程。因此，根据弹性力学理论建立的控制微分方程组，需要离散化。所谓离散化是将一个函数用假设干点的函数值近似表示。Y=f(x)xyy0y1y2y3y4y5离散化成为数值方法的关键

11§3.1有限差分法〔FiniteDifferencemethod〕1.差分格式：有限差分方法的根本思想是用差分代替微分，将微分方程的求解转化为代数方程的求解，求得的是待求函数〔定义在定义域的每一点〕在离散的网格点上的值。为了说明差分方法的思想，先回忆一维问题的差分和微分

12§3.1有限差分法〔FiniteDifferencemethod〕差分近似：考虑一个二维问题。假设：是一物理问题的解，定义在一个矩形区域a1034567891011122xyhhb为了使用差分法，首先要把所关心的区域剖分为网格，对网格节点给以编号。对于图示的矩形区域，用等距的均匀网格是最容易的

13在０点附近沿x轴线上点的函数值如果网格很细密，很小在节点1和3上：分别为和联立求解，可获得：1034567891011122xyhhb差分近似

14同理，沿y方向取点，可得：差分公式：同理，可得到混合二阶偏导数：差分近似1034567891011122xyhhb

15同理，可获得更高阶的导数。例如差分近似1034567891011122xyhhb求解以上4式，得到

16同理，可得差分近似1034567891011122xyhhb导数的阶次越高，涉及的周围节点数越多以上的格式是中心差分格式

17向前差分例如：用，，表示：向后差分：不同的差分方式，对应不同的差分方法1034567891011122xyhhab差分近似

18简单问题对以下边值问题建立差分格式物理背景：有热源、包含指定温度和指定热流边界的一维稳态温度场；有分布力作用下的杆的变形

19简单问题的解以上最后一式代入4点处的第2类边界条件，得到一式各式代入1-4点处的微分方程，得到线性方程组