工业产品材料力学设计》-第2章-拉压有限元法.pptVIP

* 工业产品材料力学设计 机械设计制造及其自动化专业 二年级第2学期 主讲教师：王一军 2011.3 TEL:687648 EMAIL:WANGYJ67@SINA.COM 2.4 拉压问题的有限单元法 第2章 轴向拉伸和压缩 力学是数学的乐园， 因为我们在这里获 得了数学的果实。 -Leonardo de Vinci 第2章 轴向拉伸和压缩 2.4.1 有限元法概述 2.4 拉压问题的有限单元法 Finite Element Method 有限元法就是将复杂的几何或受力对象划分为数量有限、形状简单、具有共同特征的子域，这些子域叫单元。 其核心思想是离散化，即将实际可变形体假想地离散为单元，然后采用能量原理进行单元分析，最后进行单元组装，完成复杂几何或受力对象的位移、应变和应力分析。 element(单元) node (节点) 第2章 轴向拉伸和压缩 2.4.1 有限元法概述 2.4 拉压问题的有限单元法 常用的一些典型单元 一维杆单元。可采用材料力学理论进行分析 二维板壳或三维体单元。以结构力学和弹性力学为理论基础。 第2章 轴向拉伸和压缩 2.4.1 有限元法概述 2.4 拉压问题的有限单元法 因此，在研究拉压杆问题时引入有限元的概念和方法，是学习有限元最好和最容易的方式，并将为以后对复杂的几何或受力对象开展有限元分析奠定坚实的理论基础。 由于： （1）杆单元是最简单的有限元分析单元； （2）拉压问题是最简单的变形体分析问题； 2.4.2 拉压杆的有限元分析 (1)单元划分与节点编号 ●单元划分：根据“三同原则” →标准杆单元 ①②③ ●节点编号：每单元两端有1个节点→节点号1,2,3,4 ●单元隔离：假设节点力和节点位移→ 2.4.2 拉压杆的有限元分析 (2)杆单元分析 (2-1)单元几何与节点描述 基本变量：节点位移 (2-2)单元位移场分析 在有限元分析中，要求必须事先给出或设定位移函数。位移函数的设定方法有2种：直接法和试函数法(Trial function) 。 位移场：即位移函数u(x)，它描述了单元各截面位移与其位置之间的关系。 直接法：就是解析法。 自主推导 试函数法(Trial function) 就是在分析前假设一个位移函数，该函数称为“试函数”。 对于如图所示的线弹性杆单元，根据解析法，我们知道其位移函数是1个线性函数，因此，我们假设的“试函数”为: 其中，a1,b1为待定系数，可根据边界条件求出。 N(x)称为形状函数矩阵 显然这是一个真实的位移函数 (2-3)单元应变场分析 (2-4)单元应力场分析 (2-5)单元变形能计算 (2-5)单元刚度方程 变形能法：卡氏第一定理 (2-5)单元刚度方程 虚功法： 虚功原理 两边分别对u1和u2求偏导数 (3)计算各单元的单元刚度方程 (4)组装各单元刚度方程 即 若已知 (5)处理边界条件 (6)求支座反力 (7)求应变与应力 经过试验：你得到如图所示的位移曲线。 你发现：无法用1个函数表达式来描述该位移曲线。 （1）关于试函数的讨论 假设你是一名企业的工程师，采用一种新型材料设计了如图所示的杆件。 你的问题是： 如果将该杆件作为一个杆单元处理，如何构造1个位移函数？ 如何分析其应变、应力？ 如何进一步构造单元刚度方程？ 2.4.3 有限元法的相关讨论 思路：无法得到真实的位移函数→构造1个来近似 数学“逼近”思想：一个复杂的函数，可以通过一系列基底函数的组合来“近似”，称为函数逼近。 近似方法一：全域逼近 基底函数 逼近函数 瑞利-里兹（Rayleigh- Ritz）法 近似方法二：分段逼近 有限元法 ☆将杆件分成n段（即划分为n个单元） ☆将每段简化为线性单元，用线性函数 定义。 ※显然，各单元也可以用其他试函数来定义，如 显然，无论全域还是分段，得到的逼近函数就是1个试函数。其中，分段逼近的优势非常明显，它可“化繁为简”： 可以把复杂形状简化为简单形状； 可以把复杂函数问题简化为简单函数问题； 可以把求解微分方程问题简化为简单的线性方程求解问题； 如果单元划分得足够小时，其分析结果将相当精确。 这就为我们分析复杂形状的物体和复杂结构提供了一个“化繁为简”的思路和方法。 在工程技术领域内，经常会遇到两类典型的问题。 第一类问题，可以归结为有限个已知单元体的组合。例如，材