有限元分析的基本术语与步骤

2016/7/1702.有限元分析的基本术语与步骤2.有限元分析基本术语与步骤1单元：有限元模型中每一个小的块体；线、三角形、四边形、四面体、六面体。。。节点：确定单元形状、表述单元特征、连接相邻单元的点；载荷：外在施加的力或力矩；不同的学科有所区别；集中力、分布力、力矩、温度、磁场。。。边界条件：结构在边界上受到的外加约束；初始条件：结构响应前所施加的初始速度、初始温度、预应力。。。有限元常用术语1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤2节点和单元节点:空间中的坐标位置，具有一定自由度和存在相互物理作用。单元:一组节点自由度间相互作用的数值、矩阵描述（称为刚度或系数矩阵)。单元有线、面或实体以及二维或三维的单元等种类。有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。载荷载荷1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤3自由度（DOFs）自由度(DOFs)用于描述一个物理场的响应特性。结构DOFs结构位移热温度电电位流体压力磁磁位方向自由度ROTZUYROTYUXROTXUZ1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤4节点和单元(续)每个单元的特性是通过一些线性方程式来描述的。作为一个整体，单元形成了整体结构的数学模型。尽管梯子的有限元模型低于100个方程（即“自由度”），然而在今天一个小的ANSYS分析就可能有5000个未知量，矩阵可能有25，000，000个刚度系数。历史典故早期ANSYS是随计算机硬件而发展壮大的。ANSYS最早是在1970年发布的，运行在价格为＄1，000，000的CDC、由Univac和IBM生产的计算机上，它们的处理能力远远落后于今天的PC机。一台奔腾PC机在几分钟内可求解5000×5000的矩阵系统，而过去则需要几天时间。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤5节点和单元(续)信息是通过单元之间的公共节点传递的。分离但节点重叠的单元A和B之间没有信息传递（需进行节点合并处理）具有公共节点的单元之间存在信息传递...AB........AB...1node2nodes1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤6节点和单元(续)节点自由度是随连接该节点单元类型变化的。JIIJJKLILKIPOMNKJIL三维杆单元(铰接)UX,UY,UZ三维梁单元二维或轴对称实体单元UX,UY三维四边形壳单元UX,UY,UZ,三维实体热单元TEMPJPOMNKJIL三维实体结构单元ROTX,ROTY,ROTZROTX,ROTY,ROTZUX,UY,UZ,UX,UY,UZ1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤7单元形函数FEA仅仅求解节点处的DOF值。单元形函数是一种数学函数，规定了从节点DOF值到单元内所有点处DOF值的计算方法。因此，单元形函数提供出一种描述单元内部结果的“形状”。单元形函数描述的是给定单元的一种假定的特性。单元形函数与真实工作特性吻合好坏程度直接影响求解精度。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤8真实的二次曲线.节点单元二次曲线的线性近似(不理想结果).2单元形函数(续)节点单元DOF值二次分布..1节点单元线性近似(更理想的结果)真实的二次曲线.....3节点单元二次近似(接近于真实的二次近似拟合)(最理想结果)..41有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤9单元形函数(续)遵循原则:当选择了某种单元类型时，也就十分确定地选择并接受该种单元类型所假定的单元形函数。在选定单元类型并随之确定了形函数的情况下，必须确保分析时有足够数量的单元和节点来精确描述所要求解的问题。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤连续体离散化单元分析整体分析确定约束条件有限元方程求解结果分析与讨论有限元法的分析过程1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤1.连续体离散化连续体：是指所求解的对象（如物体或结构）。离散化（划分网格或网络化）：是将所求解的对象划分为有限个具有规则形状的微小块体，把每个微小块体称为单元，相邻两个单元之间只通过若干点互相连接，每个连接点称为节点。相邻单元只在节点处连接，载荷也只通过节点在各单元之间传递，这些有限个单元的集合体，即原来的连续体。单元划分后，给每个单元及节点进行编号；选定坐标系，计算各个节点坐标；确定各个单元的形态和性态参数以及边界条件等。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤1.连续体离散化所示为将一悬臂梁建立有限元分析模型的例子。该悬臂梁划分为许多三角形单元；三角形单元的三个顶点都是节点。悬臂梁及其有限元模型1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤2.单元分析连续体离散化后，即可对单元体进行特性分析，简称为单元分析。单元分析工作主要有两项：(1)选择单元位移模式(位移函数)用节点位移来表示单元体内任一点的位移、应变和应力，就需搞清各单元中的位移分布。一般是假定单元位移是坐标的某种简单函数，用其模拟内位移的分布规律，这种函数就称为位移模式或位移函数。通常采用的函数形式多为多项式。根据所选定的位移模式，就可以导出用节点位移来表示单元体内任一点位移的关系式。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤进行单元力学特性分析，将作用在单元上的所有力（表面力、体积力、集中力）等效地移置为节点载荷；采用有关的力学原理建立单元的平衡方程，求得单元内节点位移与节点力之间的关系矩阵单元刚度矩阵。(2)分析单元的特性，建立单元刚度矩阵3.整体分析把各个单元的刚度矩阵集成为总体刚度矩阵，以及将各单元的节点力向量集成总的力向量，求得整体平衡方程。集成过程所依据的原理是节点变形协调条件和平衡条件。2.单元分析1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤4.确定约束条件由上述所形成的整体平衡方程是一组线性代数方程，在求解之前，必修根据具体情况分析，确定求解对象问题的边界约束条件，并对这些方程进行适当修正。5.有限元方程求解通过求解整体平衡方程，即可求得各节点的位移，进而根据位移可计算单元的应力及应变。6.结果分析与讨论1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤2.有限元法基本求解方法(1)位移法以节点位移作为基本未知量，通过选择适当的位移函数，进行单元的力学特性分析。在节点处建立单元刚度方程，再组合成整体刚度矩阵，求解出节点位移后，进而由节点位移求解出应力。位移法优点是比较简单，规律性强，易于编写计算机程序。所以得到广泛应用，其缺点是精度稍低。(2)力法该法是以节点力作为基本未知量，在节点处建立位移连续方程，求解出节点力后，再求解节点位移和单元应力。力法的特点是计算精度高。(3)混合法此法是取一部分节点位移和一部分节点力作为基本未知量，建立平衡方程进行求解。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤17矩阵位移法基本思路简单概括为：“先分再合，拆了再搭”根据位移条件和平衡条件将离散的单元组合成原结构，进行整体分析——建立结点力与结点位移之间的关系（结构刚度方程）。将结构离散成有限的单元，进行单元分析——建立杆端力与杆端位移之间的关系（单元刚度方程）。解算刚度方程，完成结构计算。1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤18“先分再合，拆了再搭”21yxo)(P11)(P22)(P331i2i23123211i2i)(F(1)1(1)1)(F(1)2(1)2)(F(2)2(2)2)(F(2)1(2)1)(P11)(P22)(P331矩阵位移法的算例EIil1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤SpaceStructureResearchCenter,HIT,CHINA1911i)(F(1)1(1)1)(F(1)2(1)2(1)21(1)11(1)2(1)21(1)11(1)1δ4iδ2iFδ2iδ4iF22i)(F(2)2(2)2)(F(2)1(2)1(2)22(2)12(2)2(2)22(2)12(2)1δ4iδ2iFδ2iδ4iF单元1单元2(e)(e)(e)δkF写成矩阵形式(1)2(1)11111(1)2(1)1δδ4i2i2i4iFF(2)2(2)12222(2)2(2)1δδ4i2i2i4iFF单元1单元21、单元分析（物理条件）1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤203(2)22(2)1(1)21(1)1ΔδΔδδΔδ0FPM0FFPM0FPM(2)233(2)1(1)222(1)111位移条件平衡条件23211i2i)(F(1)1(1)1)(F(1)2(1)2)(F(2)2(2)2)(F(2)1(2)1)(P11)(P22)(P3312、整体分析1有限元分析常用术语2有限元分析步骤3(2)22(2)1(1)21(1)1ΔδΔδδΔδ位移方程平衡方程(1)21(1)11(1)2(1)21(1)11(1)1δ4iδ2iFδ2iδ4iF(2)22(2)12(2)2(2)22(2)12(2)1δ4iδ2iFδ2iδ4iF物理方程将位移方程代入物理方程后再代入平衡方程，可得：3322223222211112111P)Δ4iΔ(2iP)Δ2iΔ(4i)Δ4iΔ(2iP)Δ2iΔ(4i0FPM0FFPM0FPM(2)233(2)1(1)222(1)1111有限元分析常用术语2有限元分析步骤将上方程组写成矩阵的形式32132122221111PPPΔΔΔ4i02i2i)4i(4i2i02i4i简写为：PΔK称为“整个结构的刚度方程”。3322223222211112111P)Δ4iΔ(2iP)Δ2iΔ(4i)Δ4iΔ(2iP)Δ2iΔ(4i1有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤23结论：将单元刚度矩阵中的元素或子块，按其整体编码的下标，“对号入座、同号相加”组集整体刚度矩阵。按矩阵形式来表达：单元刚度矩阵整体刚度矩阵(1)22(1)21(1)12(1)11(1)kkkkk1221(2)33(2)32(2)23(2)22(2)kkkkk2332(2)33(2)32(2)23(2)22(1)22(1)21(1)12(1)11kk0kkkkkkK0233211222211114i2i02i4i4i2i2i4i02332111有限元分析常用术语2有限元分析步骤2.有限元分析基本术语与步骤1有限元分析常用术语2有限元分析步骤243、支承条件的处理:（1）后处理法概念：（2）支承条件的引入“主1副零”法32132122221111MMMθθθ4i2i02i)4i(4i2i02i4i原刚度方程：引入支承条件后212121111MMθθ)4i(4i2i2i4i为便于编程，保持原矩阵行列不变先不考虑支承条件建立整个结构的刚度方程，而后再引入支承条件修改刚度方程的方法。