第二部分 平面杆系结构的线弹性有限元法.docVIP

§3 平面杆系结构的线弹性有限元法 §3.1　概论 在有限元法中，可以采用位移法，也可以采用力法或混合法。其中提出最早并且应用最广的是位移法。对于平面杆系结构来说，位移法实际上就是结构力学中的矩阵位移法（也称刚度法），在计算时以结点位移作为基本未知量。 杆系结构的矩阵分析实际上就是有限元法。其基本思路是：先把结构离散成有限个数目的单元，然后再考虑某些条件，将这些离散的单元重新组合在一起进行分析计算。这样使一个复杂的计算问题转化为简单的单元分析和集合问题。根据这个思路，杆系结构的有限元法可分为两大步骤： （1）单元分析。研究单元的受力与变形之间的关系； （2）整体分析。研究如何将这些离散的单元重新组合得到与实际问题相符合的（如边界条件、外界荷载等等）的计算模型—整体刚度方程。 在有限元中，一般采用矩阵形式进行分析求解，因为矩阵运算不仅使公式非常紧骤，而且形式统一，易于编程，适合在电子计算机上进行自动求解。因此，在有限元法的一般格式中，应尽量采用矩阵形式进行运算。 　§3.2 局部坐标系下的单元刚度矩阵 1　单元的划分。在杆系结构的有限元法中，一般将由相同材料、具有相同横截面的一根杆件（即等截面直杆）当成一个单元，整个结构就是由有限个杆件单元组成的集合体。杆件单元具有2个结点，即首结点和末结点，但一般是先确定结点的位置，结点一旦确定，则结点之间的单元也就确定了。在进行杆系结构的单元划分时，应注意如下事项： 结点位置的确定。结点一般选在杆件的如下位置：杆件的转折点、杆件汇交点、支承点、截面或材料的突变点，这些点都是结构的构造点，有时为了使结构只承受结点荷载，在集中荷载的作用处也设置一个结点。 结点的编号。为了使集合以后的总刚的带宽最小，一般应遵循尽量使相关结点（有单元相连的结点）编号差值的最大值最小的原则进行。 2　单元刚度矩阵 考虑一等截面的平面梁单元，单元首末结点分别为，单元长为，单元抗弯刚度为，为材料的弹性模量，是截面的抗弯惯矩，取轴为沿梁单元中心轴，轴与轴成90o，如图1所示。 位移模式和形函数 如果不考虑杆件的轴向变形与横向弯曲变形的互相影响，且设轴向的位移（即单元轴向位移）取为的线性函数，而对于轴向的位移（即单元横向位移，亦即梁单元挠度）取为的三次多项式表示。于是有： 　　 (1) 不考虑剪切变形的影响，即假定，可推导出位移模式又可表示为 　　（2），其中； 式中、、、、、反映了单元的位移形态，称为形态函数，简称形函数，称为形函数矩阵。 用结点位移表示的应变和应力 如单元产生拉压变形和弯曲变形，则其纵向纤维的线应变可分成两个部分： （1）－称为拉压应变，也称为轴向应变，整个梁截面都相同； （2）－弯曲应变，沿梁高因存在曲率（称为广义弯曲应变）而不同。 假设不考虑剪切变形，由梁单元的几何方程可知 　　（3） 在杆件截面高度处，有曲率引起的应变为 　　 （4） 由（3）和（4）式可写成 （5） 式中表示对求一阶导，表示对求二阶导，其余说明类似。 上式可写成 　（6） 由虎克定律，就可以得到利用结点位移表示的应力表达式 （7） 式中称为弹性矩阵，且 由式（6）可得 　　（8） 式中，称为应力矩阵。 由虚位移原理导出梁单元的刚度方程 应用虚功原理于梁单元上，可得到梁单元保持平衡的刚度方程为 　（9） 令 　 （10） 　　 （11） 式中称为单元刚度矩阵，为由于分布荷载而移置的等效结点荷载，为直接作用在结点上的荷载。 于是 　　　 （12） 令 　　　　　（13） 则 　　　　　　　（14） 上式即为局部坐标系中的单元刚度方程。 将的具体的表达式代入（12）式，经积分和矩阵运算可得到平面梁单元的单元刚度矩阵为,具体表达式如下: 　　（15） 式(13)中元素为平面杆系单元刚度矩阵元素，即： 显然单元刚度矩阵为对称矩阵。§3.3 整体坐标系下的单元刚度矩阵 在前面的分析之中，单元刚度矩阵是在单元的局部坐标系中形成的，由于各个单元的局部坐标系不同，因此必须将每个单元的单刚转换到同一个公共的坐标系下，这个公共坐标系就是整体坐标系。为了区别起见，在局部坐标系下的杆端分量符号顶上加“－”。下面首先介绍转换矩阵的概念，然后据之建立整体坐标系下的单元刚度矩阵。 一　转换矩阵 　如图3－2所示，任一单元的首端结点力在两种坐标系中的分量。其中图(a)表示在局部坐标系中的三个分量、和，图(b)表示在整体坐标系中的三个分量、和。为了导出、、与、、之间的关系式，在图中将两个力、分别投影在和轴上，可得出下式中的前两式： 其中第三式表明，在两个坐标系中的力偶仅是仍彼此相等。表示由轴转到轴的角度，以逆时针方向为正。同理对单元的另一端力也可得出类似的关系： 把以上两个方程组组合成一个矩阵方程得 或简写成 式中是