(精)CAD课件6——最新.ppt

6.1 引言 6.2 有限元单元类型及形函数 6.3 热传导有限元分析 6.4 弹性有限元法 6.5 弹塑性有限元法 6.6 刚塑性有限元法 6.1 引言 6.1.1 材料成形CAE概述 6.1.2 材料成形过程数值方法 6.1.1 材料成形CAE概述 材料成形主要是指现代金属材料和非金属材料的加工，即采用铸造、锻压、焊接和注射等方法将原材料加工成所需的形状、尺寸，并达到一定的组织性能要求，这一过程称为材料成形。 计算机模拟即是通过数值计算得到用微分方程和边值问题来描述材料成形问题中工件和模具的速度场(位移场)、应变场、应力场、温度场等，据此预测工件成形过程中组织性能的变化以及可能出现的缺陷。 6.1.2 材料成形过程数值方法 1. 有限元法的概念 2. 有限元法的基本步骤 3. 直接公式法 4. 最小总势能公式法 5. 加权余数法 有限元法的概念 有限元法是用有限个单元将连续体离散化，通过对有限个单元作分片插值求解各种力学、物理问题的数值方法。 有限元法把连续体离散成有限个单元：单元是各种形状(如三角形、四边形、六面体等)的单元体。有限元法已被用于求解线性和非线性问题，并建立了各种有限元模型。 2. 有限元法的基本步骤 (1) 建立求解域并将之离散化成有限元，即将问题分解成节点和单元。 (2) 假设代表单元物理行为的形函数，即假设代表单元解的近似连续函数。 (3) 对单元建立方程。 (4) 将单元组合成总体的问题，构造总体刚度矩阵。 (5) 应用边界条件、初值条件和负荷。 (6) 求解线性或非线性微分方程组，以得到节点的值。 (7) 后处理，得到其他重要的信息。 3. 直接公式法 下面用一实例说明直接公式法解题的步骤和过程。 【例6.1】考虑带有负荷P的变横截面杆，如图6.1所示。 1) 将问题域离散成有限的单元 首先将问题分解成节点和单元。用4个节点和3个单元的模型代表杆，如图6.2所示。 2) 假设近似单元行为的近似解 为了研究典型单元的行为，考虑一个等横截面为A的实体的位移量，单元的长度为l，承受的外力为F，如图6.3所示。 实体的平均应力 (6-1) 实体的平均应变 (6-2) 在变形区域内，应力和应变服从虎克定律 (6-3) 式中，E为物体的弹性模量。 结合(6-1)、(6-2)、(6-3)式，有 (6-4) 横截面相等的实体可以视为一个弹簧，其等价刚度为 (6-5) （6-6) (6-7) 表示成矩阵形式，有 (6-8) 或 （6-9) 即 {反作用力矩阵}＝[刚度矩阵]{位移矩阵}-{负荷矩阵} 3) 对单元建立方程 对每个单元建立两个方程。如图6.5所示。 4) 将单元组合起来表示整个问题 将方程(6-11)描述单元的方法应用到所有单元，并把它们组 合起来即可得到总体刚度矩阵。单元(1)的刚度矩阵表示如下 类似地，对于单元(2) 和(3)，有 5) 应用边界条件和负荷 边界条件u1＝0。在节点4处应用外力P (6-13) 有限元公式可写成如下的一般形式 [刚度矩阵]{位移矩阵}={负荷矩阵} 6) 求解代数方程 7) 获取其他信息 在例6.1中，按图6.6所示截取截面可得到单元的平均应力 (6-15) 应变 (6-16) 4. 最小总势能公式法 当实体拉伸量为时，物体内储存的能量为 (6-17) 将方程写成标准应变和应力形式 (6-18) (6-19) 最小总势能原理表明对于一个稳定的系统，平衡位置会出现位移以使系统的总势能最小。或者说在所有几何可能位移中，真实位移使得总势能取最小值。 对于例6.1，用最小总势能公式和用直接公式法得到的总体 刚度矩阵是完全相同的 进一步，应用边界