有限单元法基础介绍.ppt

《有限单元法基础介绍.ppt》由会员分享，可在线阅读，更多相关《有限单元法基础介绍.ppt（41页珍藏版）》请在三一办公上搜索。

1、有限单元法基础介绍,南京农业大学工学院机械工程系,2023/8/17,南京农业大学工学院机械工程系,问题的引出,问题一：上图所示为一个悬臂梁，其中长度L=2000mm，高度H100mm，宽度B50mm，在力P1000N作用下，求该梁A点的挠度。已知:材料弹性模量E=2.01011Pa，泊松比为0.3。,解答：根据材料力学公式，其中，为矩形截面的惯性矩，代入已知条件，计算结果为：,问题：,2023/8/17,南京农业大学工学院机械工程系,有限元单元法简介,有限单元法(FEM)是20世纪50年代以来随着计算机的广泛应用而发展起来的一种现代数值解法。该方法首先应用在连续力学领域飞机结构静、动态特性分

2、析中。随后很快就广泛应用于求解传导、电磁场、流体力学等连续性问题。有限元分析是利用数学近似分析方法对真实物理系统（几何、载荷工况）进行模拟，利用简单而又相互作用的元素，即单元，用有限数量的未知量去逼近无限未知量的真实系统。,Definition,历史典故结构分析的有限元方法是在二十世纪五十年代到二十世纪六十年代由一批学术界和工业界的研究者创立的。有限元分析理论现在已成为悬索桥和蒸汽锅炉进行手算评核的基础。,2023/8/17,南京农业大学工学院机械工程系,有限单元法力学基础,各力学学科分支的关系,非变形体（刚体）,变形体,2023/8/17,南京农业大学工学院机械工程系,有限单元法几种典型的分

3、析对象,双向拉索悬索桥,发动机有限元模型,齿轮接触有限元分析,汽车碰撞实验模拟（福特）,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基本思想,有限单元法的思想是将物体（连续体）离散成有限个按一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析法。物体被离散后，通过对其中各个单元进行单元分析，然后再进行整体分析，最终得到整个物体的分析结果。,Definition,2023/8/17,南京农业大学工学院机械工程系,有限单元法的思路及技术路线,以力学为基础，分析任意变形体的受力情况，运用数值近似分析方法

4、，以计算机为工具进行分析计算，以获得该结构所有的力学信息，并使得该方法能够普及、简单、高效、方便，一般工程技术人员熟练掌握后就可以使用。,思路,标准化（研究力学理论：任意复杂问题标准化分解，单元建模有限种标准单元类型）规范化（前处理：CAD几何、力学建模、求解、后处理显示）计算机程序化（标准程序、模块）应用的规模化、普及性（可求解大型计算问题108-1010DOF),技术路线,2023/8/17,南京农业大学工学院机械工程系,物理系统举例,物理系统 几何体 载荷,2023/8/17,南京农业大学工学院机械工程系,有限元模型,有限元模型 是真实系统理想化的数学抽象。,Definition,真实系

5、统,有限元模型,2023/8/17,南京农业大学工学院机械工程系,有限元单元模型中几个重要概念,单元网格划分中每一个小的块体节点确定单元形状、单元之间相互联结的点节点力单元上节点处的结构内力载荷作用在单元节点上的外力（集中力、分布力）约束限制某些节点的某些自由度弹性模量（杨式模量）E泊松比（横向变形系数）密度,单元,单元,载荷,节点,节点力,约束,2023/8/17,南京农业大学工学院机械工程系,单元的自由度（DOFs）,自由度(DOFs)用于描述一个物理场的响应特性。,结构 DOFs,2023/8/17,南京农业大学工学院机械工程系,节点和单元,载荷、约束和力等信息是通过单元之间的公共节点传

6、递的。,分离但节点重叠的单元A和B之间无法进行信息传递（需进行节点合并处理）,具有公共节点的单元之间存在信息传递,.,.,.,A,B,.,.,.,.,.,.,.,.,A,B,.,.,.,2 nodes,2023/8/17,南京农业大学工学院机械工程系,节点和单元(续),节点自由度是随连接该节点 单元类型 变化的。,J,I,I,J,J,K,L,I,L,K,I,P,O,M,N,K,J,I,L,三维杆单元(铰接),UX,UY,UZ,三维梁单元,二维或轴对称实体单元,UX,UY,三维四边形壳单元,三维实体热单元,TEMP,J,P,O,M,N,K,J,I,L,三维实体结构单元,UX,UY,UZ,ROTX

7、,ROTY,ROTZ,UX,UY,UZ,UX,UY,UZ,ROTX,ROTY,ROTZ,2023/8/17,南京农业大学工学院机械工程系,物体变形及受力情况的描述,基本变量（位移）（应变）（应力）基本方程力的平衡方程几何方程物理方程求解方法经典解析半解析传统数值解法现代数值解法（计算机硬件、规范化、标准化、规模化）,三大方面,三大方程,即：,=E E 弹性模量,2023/8/17,南京农业大学工学院机械工程系,三大基本方程,根据静力学、几何学和物理学三方面条件，建立三套方程。平面问题中，根据微分体的平衡条件，建立平衡微分方程：(6)根据微分线段上形变与位移之间的几何关系，建立几何方程：(7)根

8、据应力与形变之间的物理关系，建立物理方程：(8-1)(8-2),2023/8/17,南京农业大学工学院机械工程系,有限单元法解题的一般步骤,物体结构的离散化将连续的结构离散成有限个单元，并在每一单元中设定有限个节点，将连续体看作只在节点处相连接的一组单元的集合体。单元特性分析，选择位移模式选定场函数的节点值作为基本未知量，并在每一单元中假设一近似插值函数已表示单元中场函数的分布规律。单元组集，建立平衡方程求解节点的位移利用力学中的某种变分原理去建立用以求未知量节点的有限单元法方程。求解单元中的应力和应变,LessonObjectives,2023/8/17,南京农业大学工学院机械工程系,有限元

9、单元法分析步骤（一）,结构离散化 将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点连接。结构的离散化是有限元法分析地第一步，关系到计算精度和效率，包括以下三个方面：单元类型的选择。选定单元类型，确定单元形状、单元节点数、节点自由度数等。单元划分。网格划分越细，节点越多，计算结果越精确，但计算量越大。网格加密到一定程度后计算精度提高就不明显，对应应力变化平缓区域不必要细分网格。节点编码。注意：有限元分析的结构已不是原有的物体或结构物，而是由同样材料、众多单元以一定方式连接成的离散物体。所以，用有限元分析计算所获得的结果是近似的。,2023/8/17,南京农业大学工学院机械工程系,有

10、限元单元法分析步骤（二）,单元特性分析 选择未知量模式选择节点位移作为基本未知量时，称为位移法，在有限元计算中位移法应用较多；选节点力作为基本未知量时，称为力法；取一部分节点位移和一部分节点力作为未知量，称为混合法。分析单元力学性质根据单元材料性质、形状、尺寸、节点数目、位置等，找出单元节点力和节点位移关系式，应用几何方程和物理方程建立力和位移的方程式，从而导出单元刚度矩阵。计算等效节点力作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，即用等效力来替代所有作用在单元上的力。,2023/8/17,南京农业大学工学院机械工程系,选择位移模式,在有限元法中，虽然整个连续物体已经变化

11、为离散化结构，但是每个单元仍然作为一个连续的、均匀的，各向同性的弹性体。对于每个单元，要计算内部的应变和应力，要求得出该单元中的位移函数。因此可以对单元假定一个位移插值函数，或称之为位移模式，得到用节点位移表示单元体内任一点的唯一的关系式。有了位移模式，就可利用几何关系和应力-应变关系表出用单元节点位移表示单元中应变和应力的表达式。,Definition,2023/8/17,南京农业大学工学院机械工程系,有限元单元法分析步骤（三）,整体分析集成整体节点载荷矢量。结构离散化后，单元之间通过节点传递力，作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，形成等效节点载荷。将所有节点载

12、荷按照整体节点编码顺序组集成整体节点载荷矢量。组成整体刚度矩阵，得到总体平衡方程：将线性代数方程组，引进边界约束条件，解总体平衡方程可求得所有未知的结点位移。通过上述分析可以看出有限单元法的基本思想是“一分一合”，分是为了进行单元分析，合是为了对整体的结构进行综合分析。,2023/8/17,南京农业大学工学院机械工程系,计算单元中的应变和应力,依据求得的结点位移，由可求得单元中任一点的应变和应力。,LessonObjectives,2023/8/17,南京农业大学工学院机械工程系,用有限元法解决工程问题的方式,利用有限元基本理论及相关数学理论，自行编程实现。特点：仅能解决某类特定问题，且编程繁

13、琐，适用于求解规模不大，或者某些专业性较强的问题。利用现有成熟的专用或通用有限元商业软件进行分析计算。特点：能解决绝大多数工程问题，建模和计算效率较高，但需要对软件进行深入的学习，了解各个功能参数的设置。,2023/8/17,南京农业大学工学院机械工程系,平面有限元解法编制计算机程序界面,2023/8/17,南京农业大学工学院机械工程系,平面有限元解法计算机程序计算结果,2023/8/17,南京农业大学工学院机械工程系,通用有限元计算程序ANSYS计算结果,2023/8/17,南京农业大学工学院机械工程系,通用有限元计算程序ANSYS计算结果,2023/8/17,南京农业大学工学院机械工程系,

14、有限元单元法基础理论（结构静力学问题）,矩阵分析法对杆系结构举例分析水平杆单元刚度矩阵,珩架,水平杆单元,杆单元两端各有一个水平节点位移 和，两端节点力分别为 和。杆的受力情况分为两种情况考虑：，此时j点被固定；，此时i点被固定。,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,，j 点被固定时；单元应变：单元应力：（材料力学中以拉应力为正，而有限元中，以向右的节点力为正）单元左端节点力：单元右端节点力：，i点被固定时，与上状态相反，因此单元左端节点力：单元右端节点力：,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学

15、问题）,把以上两种结果叠加起来，得到左、右两端都有位移情况下单元节点力：,写成矩阵形式：,称为单元刚度矩阵。,单元轴力可写为：,称为单元应力矩阵。,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,实际，节点i和j 除了水平位移外，还可以产生垂直位移(但在小变形下，垂直节点位移对杆内力无影响）。引入垂直节点位移，和垂直节点力，单元刚度矩阵扩展为：,或,称为单元刚度矩阵。,节点力，,节点位移,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,倾斜单元刚度矩阵,为正交矩阵，其中,局部坐标,与整体坐标x,y之间的位移 与

16、 之间存在如下变换关系：,式中，转换矩阵,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,局部坐标系中节点力 与整体坐标系中节点力 之间的关系为：,其中，i点 节点力,得出：,局部坐标系中节点力,其反映单元节点位移与单元节点力关系，称为单元刚度方程。,将上式记为,，j点 节点力,i点 节点位移,，j点 节点位移,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,刚度系数Kij的意义是节点j的单位节点位移在节点i上产生的节点力,从一个珩架中取一节点i，设环绕该点有三个单元，即ij、im、ip。该节点承受水平和垂直载

17、荷分别为Xi和Yi，即,节点平衡方程与整体刚度矩阵,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,根据力的平衡原理，节点受力平衡方程为：,杆单元ij在节点i的节点力为：,其它单元施于结点i的节点力为：,其它节点都有以上平衡方程，对于全部节点i1,2,N的结构,得到2N阶线性方程组,其中，,，为全部节点位移组成的列阵；,，为全部节点载荷组成的列阵；,，为整体刚度矩阵。,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,总体刚度矩阵的合成,采用大域变换矩阵法把单元刚度矩阵合成结构整体刚度矩阵。,结构总体刚度矩阵 与

18、单元刚度矩阵 之间的关系为,其中，为单元大域变换矩阵，对平面珩架结构，单元自由度m=4,节点自由度h2，整个结构有n个节点，则该单元大域变换矩阵为m(hn)维。其中ij单元假定为全局单元编号中第3个，其大域变换矩阵为,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,总体结构的载荷矢量、位移矢量与单元载荷矢量、位移矢量之间的关系为：,边界条件的处理,边界条件指结构边界上所受到的外加约束。边界上的节点通常有两种情况：一种可以自由变形，如图中5、6、7、8；如果节点3作用外载荷Q，可以令该点的载荷为Q；另一种边界上的节点，已经规定了节点的位移数值。如：,20

19、23/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构静力学问题）,把结构平衡方程组重新排列，得到如下方程：,式中，是已知节点位移，是未知节点位移；,相应地，是已知节点载荷，是未知支点反力。,只要已给出的位移 足以阻止结构的刚体移动，则子矩阵 将是非奇异的，可以解出未知的节点位移：,进而求出未知支点反力：,但在有限单元法中，未知量的个数通常有几百个，甚至几十万个，一般利用计算机求解。,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构动力学问题）,结构动力学问题有限元方法运动状态中各节点的动力平衡方程为：,式中，分别为惯性力、阻尼力和动力载荷，为弹性

20、力。,弹性矢量可用节点位移 和刚度矩阵 表示为：,根据达朗贝尔原理，可用质量矩阵 和节点加速度 表示惯性力如下：,设结构具有粘滞阻尼，可用阻尼矩阵 和节点速度 表示阻尼力如下：,得到运动方程：,即,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构动力学问题）,单自由度系统的阻尼单自由度系统的自由振动方程为：,式中，质量；阻尼系数；为刚度系数；为位移。,两边除以m,得到：,设初始条件为：当t=0时，符合这些初始条件的解为：,系统的自振频率为，其振幅随着时间而逐渐衰减。,其中，称为阻尼比，为系统的自振频率（角频率）。,大多数结构的阻尼比都比较小，较多为，阻尼对结构自振频

21、率的影响是很小的，通常取。,2023/8/17,南京农业大学工学院机械工程系,有限元单元法基础理论（结构动力学问题）,结构自振频率与振型结构自由振动方程为：,实际中，阻尼对结构自振频率和振型影响不大，因此可忽略阻尼力，得到无阻尼自由振动方程：,设结构作简谐运动：,对于每个自振频率，可确定一组各节点的振幅值它们之间应保持固定的比值，但绝对值可以任意变化，它们构成一个矢量，称为特征矢量，在工程上通常称为结构的振型。,代入上式得到方程：,在自由振动时，结构中各节点的振幅 不全为零，所以结构自振频率方程为：,上式是关于 的n次代数方程，可求出结构的自振频率：,2023/8/17,南京农业大学工学院机械工程系,THE END,