塑性加工过程CAE-02FEM基础.ppt

第三章 板料成形的有限元基础 §3.1 序言 §3.1 序言 圆面积的近似求解： 工程中有限元的应用 机械、航空、土木、汽车工程 结构分析（静态/动态、线性/非线性） 热/液体流动 电磁 地质力学 生物力学 … 有限元的简要发展历史 1943 —— Courant 1956 —— Tuner, Clough, Martin, Topp (Stiffness) 1960 —— Clough (“Finite Element” 平面问题) 1970s —— 在大型机上得到应用 1980s —— 在微机上得到应用，前后处理软件 1990s —— 大型机构系统分析 §3.1.1 基本概念 结构分析中的有限元（过程） 将结构体划分成小片（单元，节点） 形成描绘物理特性的单元刚度 将单元组装成一个整体结构的近似方程组 求解已经引入位置物理量（位移）的方程组 计算单元中用户所关心的物理量（应变，应力） 计算机 前处理（建立有限元模型，载荷和约束） 有限元求解器（组装和求解系统方程） 后处理（显示计算结果） 商业化有限元软件包 ANSYS，NASTRAN，ALGOR(通用目的) ABAQUS(非线性动力分析) PATRAN，HyperMesh(前后处理) LS-DYNA(碰撞动力分析) DynaForm(前后处理) … 本章有限元基础课的目的 理解有限元分析的基本原理和思路 掌握本章中所涉及到的单元模型的推导和适用范围 对一个给定的问题能够建立适当的有限元模型 能够解释并评价有限元分析结果的优劣（知道问题的物理意思） 明白有限元的局限性（不要错误地应用有限元—数值工具） 线性代数方程组系统： 矩阵形式： A称为n×n矩阵，x和b分别是n维的列向量 矩阵的加法减法： 矩阵的乘法： 矩阵的转置： 对称矩阵： 单位矩阵： 矩阵行列式的值： 奇异矩阵： 如果 det A =0，那么系统存在问题（非唯一解，发散等） 矩阵的逆： 线性方程组系统的求解技术： 高斯消去法 迭代法 正定矩阵： 对于所有非零向量X，有XTAX0，A为正定矩阵 正定矩阵为非奇异矩阵 矩阵的导数和积分： 弹簧元： 考虑弹簧元的力的平衡条件： 节点i： 节点j： 矩阵形式： 注意：K为对称矩阵，K是非奇异的还是奇异的？ 弹簧元系统： 单元1： 单元2： 对整个系统进行单元刚度矩阵的组装： 考虑节点力的平衡条件： 节点1： 节点2： 节点3： 矩阵形式： K为该弹簧系统的刚度矩阵（结构矩阵） 单元刚度矩阵另一种组装方法： 分别扩大单元1和2的刚度矩阵 这与根据节点力平衡得出的矩阵是一样的。 引入边界条件和力的条件： 假设u1=0，F2=F3=P，那么 未知量为： u2，u3和F1 求解方程可得： 检查计算结果： 结构变形后的形状 外力平衡 有关弹簧元的注意事项： 适于刚度分析的计算 不适合用于弹簧本身的应力分析计算 在弹簧元的横向是否具有刚度，弹簧元是否具有扭转刚度 例子1： 已知： 求： (a) 整体刚度矩阵 (b) 节点2和3的位移 (c) 节点1和4的支反力 (d) 弹簧2的力 （a）问：分别求出单元刚度矩阵： 单元1： 单元2： 单元3： 组装后： 组装后： 注意到整体刚度矩阵是对称并带状分布的。 该系统的平衡方程为： （b）问： 将边界条件（u1=u4=0）应用到平衡方程中，去掉1行1列，4行4列后： 求解得： （c）问： 从平衡方程组中的1和4可得： （d）问： 弹簧元2的平衡方程为： 式中i=2，j=3，可以计算出弹簧力为： 例子2： 问题描述：对上述一个具有任意弹簧元节点和单元的系统，求其整体刚度矩阵。 单元拓扑关系： 上表中为每个弹簧元的局部节点号和整体节点号的对应关系。 然后依次求出每个弹簧元的刚度矩阵： 组装后： 整体刚度矩阵是对称并带状分布的。 §3.2 杆和梁单元 大部分结构分析问题都可以看作是线性静力分析问题，它们都基于以下假设： 1. 小变形（加载方式不会因为变形而改变） 2. 线弹性材料（不存在塑性和破裂） 3. 静力载荷（在结构上的载荷是慢速平稳地施加上去的） 线性静力分析可以解决结构分析问题中大部分的问题，对于大多数的结构分析问题，静力分析可以得到一个近似的结果。 线性静力分析是非线性分析的基础。 杆单元： 刚度矩阵——直接法： 假设位移u沿着杆的轴向线性分布： K即为杆单元的刚度系数，杆单元和弹簧元类似。 单元刚度矩阵为： 单元的平衡方程组为： 节点自由度： 对一维的杆单元，每个节点就只有一个自由度 刚度矩阵K中系数的物