弹性力学有限元第七章 FEM部分之一 有限单元法入门.ppt

第七章 有限单元法入门 §7-3 建立有限元计算模型的常用方法 点在杆上的位置 精确解 配置法 子域法 伽辽金法 最小二乘 0 0 0 0 0 0 2.5 0.001027 0.001076 0.001029 0.001041 0.001027 5 0.002213 0.002259 0.002209 0.002220 0.002208 7.5 0.003615 0.003660 0.003618 0.003624 0.003618 10 0.005333 0.005384 0.005330 0.005342 0.005331 第七章 有限单元法入门 §7-4 有限元分析的错误来源 类似于ANSYS的有限元综合软件包已经成为工程师应该掌握的工具。 由于缺乏应有的训练或对基本概念缺乏足够的理解，经常得出错误的结果，其来源主要有： 河海大学 机电工程学院 力学教研室 第七章 有限单元法入门 第六章空间问题的基本理论和解答 第七章 有限单元法入门 §7-1 引言 实际工程问题，我们一般得不到它们的精确解，这主要是因为微分方程组的复杂性，以及难以确定的边界条件和初值条件； 有限元FEM(Finite Element Method),FEA(Finite Element Analysis) 是一种用于求解各类工程问题的数值计算方法。 有限元方法（基于变分原理的差分方法），是一种用于求解微分方程组或积分方程组数值解的数值技术. 这一解法基于完全消除微分方程, 即将微分方程转化为代数方程组(稳定情形); 或将偏微分方程(组)改写为常微分方程(组)的逼近, 这样可以用标准的数值技术(例如欧拉法,龙格－库塔法等)求解. 第七章 有限单元法入门 §7-1 引言 有限单元法的起源 有限元法最初起源于土木工程和航空工程中的弹性和结构分析问题的研究. Alexander Hrennikoff(1941) Richard Courant(1942) 利用网格离散化将一个连续区域转化为一族离散的子区域, 通常叫做元.Hrennikoff 的工作离散用类似于格子的网格离散区域; Courant 的方法将区域分解为有限个三角形的子区域, 用于求解来源于圆柱体转矩问题的二阶椭圆偏微分方程. Courant 的贡献推动了有限元的发展, 绘制了早期偏微分方程的研究结果. 有限元方法的发展开始于五十年代中后期使用在机身框架和结构分析上，并于六十年代通过斯图加特大学的en:John Argyris和柏克莱加州大学的en:Ray W. Clough在土木工程中的应用工作中积累经验。 第七章 有限单元法入门 §7-1 引言 有限单元法分析的基本过程 建立求解域，离散化成单元(Element)和节点(Node) 假定描述单元属性的形函数(Shape Function)，即用一个近似的连续函数描述每个单元的解 建立单元刚度矩阵 组装单元，构造总体刚度矩阵 应用边界条件和初值条件，并施加荷载 求解线性或非线性微分方程组得到节点值 根据所求的节点值(位移)，计算其他关心的数据(应变，应力)。 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 假有一端承受荷载P的变截面杆，一端固定，一端承受荷载。求杆件受载后沿着杆件长度方向上不同点变形的大小。(体力不计) 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 5 Nodes, 4 Element 每个单元的横截面积，由构成单元节点处的横截面的平均面积表示 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 平均应力 平均应变 物理关系 由三者可得 该式与 线性弹簧的方程 F = kx 相似 等截面杆轴向受拉可以看做为一个弹簧，其等价刚度为： 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 作为近似，可以将杆看成是一系列轴向受拉且具有不同横截面积的构件组合。 即由4个弹簧串联而成的模型。 单元的弹性行为(第i个单元，第i,i+1个节点)： 得出弹簧单元的等效刚度： 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 利用上述模型，假定力施加在节点上，按照静力平衡条件有： 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 重组方程组，表示为矩阵形式 [反作用力矩阵]=[刚度矩阵][位移矩阵]－[荷载矩阵] 第七章 有限单元法入门 §7-2 有限单元法的‘Hello, World’ 引入边界条件 u1=0, 系统方程可以变