有限单元法基础.ppt

第一页，共十四页，2022年，8月28日 2．单元特性分析 1)选择位移模式 在有限单元法中，选择节点位移作为基本未知量称为位移法；选择节点力作为基本未知量称为力法；取一部分节点力和一部分节点位移作为基本未知量称为混合法。位移法易于实现计算自动化，所以在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元中的一些物理量如位移、应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法中我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数，如 ψi ，其中αi是待定系数，ψi是与坐标有关的某种函数。 第二页，共十四页，2022年，8月28日 2)分析单元的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。 3)计算等效节点力 物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边界传速到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，也就是用等效的节点力来替代所有作用在单元上的力。 第三页，共十四页，2022年，8月28日 3．单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新联接起来，形成整体的有限元方程 kδ=f (2—1) 式中，k是整体结构的刚度矩阵；δ是节点位移列阵；f是载荷列阵。 4．求解未知节点位移 解有限元方程式(2—1)得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。 通过上述分析，可以看出，有限单元法的基本思想是“一分一合”，分是为了进行单元分析，合则是为了对整体结构进行综合分析。 第四页，共十四页，2022年，8月28日 2.2 结构静力分析的有限单元法 2.2.1 单元特性的导出方法 进行有限元分析的基本步骤之一就是要找出所划分的单元的刚度矩阵和质量矩阵，一般来说，建立刚度阵的方法有：(1)直接方法；(2)虚功原理法；(3)能量变分原理方法；(4)迦辽金法。 下面主要叙述直接方法、虚功原理法及能量变分原理法。 1．直接方法 直接方法是直接应用物理概念来建立单元的有限元方程和分析单元特性的一种方法。这种方法仅能用于简单形状的单元，如梁单元。但它可以帮助理解有限元法的物理概念。 图2—2所示是XOY平面中的平面桁架简图，E为杆的弹性模量，A为截面积。现在，以它为例用直接方法建立单元的刚度矩阵。 第五页，共十四页，2022年，8月28日 第六页，共十四页，2022年，8月28日 简例是一平面杆系问题 由二个单元组成；三个节点；每个节点上有二个位移（u1,v1,u2,v2,u3,v3）、二个外力(fx1,fy1,fx2,fy2,fx3,fy3); 先来看单元一： 单元一有二个节点（1和2）；节点上各有二个节点力分量（fx11,fy11,fx21,fy21）、各有二个位移分量(u11,v11,u21,v21). 这里的下标是节点号；上标是单元号。 列出单元一的力与位移的关系式： 式（1－1） fx11＝k11u11 + k12v11 + k13u21 + k14v21 fy11＝k21u11 + k22v11 + k23u21 + k24v21 (1-1) fx21＝k31u11 + k32v11 + k33u21 + k34v21 fy21＝k41u11 + k42v11 + k43u21 + k44v21 第七页，共十四页，2022年，8月28日 式（1－1）可写成矩阵形式： 式（1－2） 还可简写成：{f}1 = [k]1 {δ}1 （1－3） {f}1 单元1节点力列阵; [k]1 单元1刚度矩阵 ; {δ}1 单元1节点位移列阵。 上面的式中，{f}1 是已知的；{δ}1是要求的未知量； [k]1 是单元刚度矩阵； k11、k12 、 k13 、 k14…… 这些刚度系数是有物理意义的: 令式 (1-1) 中 u11 = 1, v11 = u21 = v21 =0 ; 式 (1-1)就变成： fx11＝k11 fy11＝k21 fx21＝k31 f