有限元基础理论教程lecture.ppt

1.3 有限单元分析的基本步骤 有限单元分析归纳为以下三个基本步骤： 网格划分（离散化） 单元分析 整体分析 网格划分 平面桁架结构、平面刚架结构的的组成部分自然地就是单元，铰接点和刚梁直段的端点就是结点。 网格划分（续） 对于连续体，先进行必要的简化，再划分为有限个单元组成的离散体。单元之间通过单元结点相连接。 由单元、结点、结点连线构成的集合称为网格。 通常把三维实体划分成四面体（Tetrahedron）或六面体（Hexahedron）单元的网格，把平面问题划分成三角形（triangular)或四边形(quadrilateral)单元的网格。 网格划分（续） 网格划分（续） 带方孔的长方形钢制等厚度薄板，受到均布轴向拉力q的作用。薄板的厚度为d，拉力沿厚度方向均匀分布。 网格划分（续） 单元分析 单元分析就是在离散化得到的单元上寻找待解问题的近似解，建立由单元结点变量之间的关系式。 平面桁架结构中的直杆只承受轴向力作用，单元分析就是建立结点力和结点位移之间的关系式。 弹性力学问题的单元分析，就是建立各个单元的结点位移和结点力之间的关系式。 整体分析 由单元分析得到了局部的近似解，把全部单元组成的整体，由局部近似解得到待解问题的全局近似解，这个过程就整体分析。 对于平面桁架结构来说，在铰接点上外载荷与相关单元的结点力之间存在平衡关系，列出全部平衡方程可以得到整个结构所承受外力与结点位移之间的线性方程组。根据约束条件修改这个方程组，由修改后的方程组可以解出结点位移分量。 例1-2、平面桁架123如图1-15所示。已知各杆件的弹性模量，各杆的截面积均为，杆13的长度为100cm，载荷P=2KN。求平面桁架的位移与内力。 平面桁架单元分析 1）在局部坐标系中建立结点力分量与结点位移分量的关系矩阵。 平面桁架单元分析（续） 2）根据局部坐标与整体坐标之间的关系，建立转换矩阵[T]。 平面桁架单元分析（续） 3）在整体坐标系中建立单元结点力分量与位移分量的关系矩阵，称为单元刚度矩阵。 平面桁架单元分析（续） 平面桁架单元分析（续） 平面桁架整体分析 根据单元结点的局部编号与结点的总体自由度编号的对应关系，把单元刚度矩阵叠加得到总体刚度矩阵。 平面桁架整体分析（续） 平面桁架整体分析（续） 位移边界条件， 在总体刚度矩阵中去掉与零位移对应的行与列，得到方程组。 平面桁架整体分析（续） 单元的内力按以下公式计算， 弹性力学问题单元分析 将单元的结点位移作为基本变量，单元分析首先要为单元内部的位移确定一个近似表达式，然后计算单元的应变、应力，再建立单元中结点力与结点位移的关系式。 弹性力学问题单元分析（续） 弹性力学问题单元分析（续） 整体分析（续） 1.4有限元软件与应用实例 1.4.1 有限元软件 从二十世纪60年代中期以来，进行了大量的理论研究，不但拓展了有限单元法的应用领域，还开发了许多通用或专用的有限元分析软件。 主要的研究领域有： 大型线性方程组的解法，非线性问题的解法，动力学问题计算方法，并行算法； 高精度单元； 多物理场耦合，如热力耦合、电磁耦合； 复杂材料模型，如塑性、多孔、晶粒度等。 有限元软件（续） MSC.Nastran，结构分析 MSC.Dytran，动力学、流固耦合、强非线性分析 MSC.Marc，非线性分析，并行处理 MSC公司（/）发布，最新版本为Nastran2004 ANSYS，结构分析，多场耦合分析 ANSYS/LS-DYNA ，动力学分析 ANSYS/FLOTRAN ，2D/3D流场分析 ANSYS公司（ http：///）发布，最新版本为ANSYS 9.0 有限元软件（续） ABAQUS/Standard，线性、非线性、动力学分析 ABAQUS/Explicit，非线性、显式动力学分析 ABAQUS公司（ /）发布，2005年1月发布Version6.5。 ADINA，高度非线性分析 ADINA-F，流场分析 ADINA-T，传热分析 ADINA RD公司（ / ）发布，最新版本为ADINA8.2 有限元软件（续） DEFORM-2D，二维体积成形问题 DEFORM-3D，三维体积成形问题 DEFORM-HT，热处理 SFTC (/)发布。 AutoForm-Incremental，试模与整个冲压过程 AutoForm Hydroforming，液压胀形 AutoForm Engineering（/ ） 发布，最新版本为AutoForm4.0。 有限元软件（续） DYNAFORM/DFE ，冲压模具设计 DYNAFORM/Formability，成形分析 ETA Software 公司（/）发布，基于LS