对于梁单元ppt课件

1、三角单元的单元节点力的问题。这些 力是怎么来的？在那三个节点上不是只应当有应力吗？ 如果单元节点上是力，那么那三条边上的应力又应当怎样理解？ 相邻单元间的协调只是节点上的协调吗？相邻边的协调呢？两个相邻单元在相邻边上的应力连续吗,第五章 杆系结构单元,本章主要内容是： 结构离散为单元的有关问题 单元（局部）坐标系和结构（整体）坐标系 单元坐标系中各类杆件单元的特性：单元刚度矩阵、等价节点力矩阵等。 结构坐标系中的单元特性及坐标变换矩阵,杆系结构是由一些杆件单元组成。主要结构类型有：梁、拱、框架、桁架等，如图（5-1）所示,图（5-1,2、编号 （1）节点编号 节点编号应按正整数不间断逐点编号。

2、编号时应力求单元两端点号差最小，以便使结构刚度矩阵元素集中在主对角线附近，后面结构刚度矩阵组集中有详细说明,5.1 结构离散,1、离散方法 取杆件与杆件交点、集中力作用点、杆件与支承的交点为节点。 相邻两节点间的杆件段是单元; 杆件结构的单元一般只有2个节点,2）单元编号 单元也要逐个依次编号。谁前谁后按实际情况而定,3、记录基本信息 应建立一个数据文件（DATA.*）来记录基本信息，以便计算时调用。基本信息包括： （1）单元总数（NE）、节点总数（NJ）、节点自由度数(NDF)。 （2）弹性模量（E）、波桑系数（AMU）。 （3）单元I端节点号IO（NE）、 J端节点号JO（NE） （4）有

3、约束的节点数（ NRJ ）、有约束的节点号（KRJ(NRJ)）、受约束的自由度（KRL(NDF,NRJ),5）单元截面面积（A）、截面惯性矩（ZI） （6）节点坐标：X(NJ)、Y(NJ) （7）分布力荷载集度qx(NE)、qyi(NE)、qyj(NE) （8） 受集中力作用的节点数（MJL）、受集中力作用的节点号（NJL（MJL）、集中力数值（VJL（NDF，MJL,DATA.FRA (1) NE、NJ、NDF 25, 18, 3 (2）E、 AMU 3.25e7, 0.15 (3）IO（NE）、JO（NE） 1, 4; 4, 7; 7, 10; 10,13; 13,16; 2, 5; 5,

4、 8; 8, 11; 11,14; 14,17; 3, 6; 6, 9; 9, 12; 12,15; 15,18; 4, 5; 5, 6; 7, 8; 8, 9; 10,11; 11, 12,; 13,14; 14,15; 16,17; 17,18 (4）NRJ 、KRJ(NRJ)、 （KRL(NDF,NRJ)） 3， 1,2,3; 9*1 A(NE)、ZI(NE) (6) X(NJ)、Y(NJ) (7)qx(NE)、qyi(NE)、qyj(NE) (8) MJL, NJL(MJL), VJL(NDF,MJL) 数据填写顺序应和程序对应,4、示例,练习 对平面铰接桁架进行结构离散，并写出数据文

5、件,已知：E=2.1109kN/m2 , P1=P3=10kN, P2=50kN, 15cm2 (斜杆） A = 65cm2 (上下弦杆） 40cm2 (竖杆,4.2 单元（局部）坐标系,杆系结构单元在结构中的位置是复杂的。如图（5-2）桁架所示,图（5-2,如果每一个单元都在统一的整体坐标系XY中写单元刚度矩阵, 可能导致结构中处于不同位置的同一类型单元，其单元刚度矩阵不相同。这不利于计算机编程运算,容易理解，采用适合于单元具体方位的坐标系将会改善上述状况，得出规格化的结果。这种属于每个单元的坐标系称为单元坐标系，也称局部坐标系,杆系结构单元主要有铰接杆单元和梁单元两种类型, 它们都只有2个

6、节点i、j,约定：单元坐标系的原点置于节点i；节点i到j的杆轴（形心轴）方向为单元坐标系中x轴的正向。 y轴、z轴都与x轴垂直，并符合右手螺旋法则（图5-3,图5-3,对于梁单元， y轴和z轴分别为横截面上的两个惯性主轴,为了便于对单元坐标系中的单元特性进行识别，引入以下符号,e单元坐标单元位移 Fe单元坐标单元力 ke单元坐标单元刚度矩阵,下面，开始讨论几种杆系结构单元在单元坐标中的一些特性,5.3 铰结杆单元,图5-4示出了一维铰结杆单元，横截面积为A，长度为l，弹性模量为E，轴向分布载荷为qx。单元有2个结点i，j，单元坐标为一维坐标轴x,1、一维铰接杆单元,5-1,单元力向量为,5-2

7、,1）位移模式和形函数,位移模式,单元结点位移向量为,因为只有2个结点，每个结点位移只有1个自由度，因此单元的位移模式可设为,5-3,式中a1、a2为待定常数，可由结点位移条件确定 x=xi时， u=ui ； x=xj时， u=uj,由此可确定a1、a2 。再将其代入式（5-3），得,5-4,形函数,将式（5-4）改写为下列形式,5-5,式中e由式（5-1）确定，形函数N为,5-6,2）应变矩阵,一维铰结杆单元仅有轴向应变,将式（5-5）、（5-6）代入上式，得,上式也可写为,5-7,式中B为应变矩阵,5-8,由应力应变关系,3）应力矩阵,将式（5-7）代入上式，得,5-9,式中S为应力矩阵,

8、5-10,4) 等价节点力,单元上作用分布力qx，则等价节点力计算公式仍为以下形式,当分布力集度qx为常数时，有,5-11,式（5-11）概念是将分布力引起的合力按静力等效原则分配到单元节点上。由于位移模式是线性函数，因此按公式（5-11）计算结果与静力等效分配是一致的,5) 单元坐标单元刚度矩阵,单元坐标单元刚度矩阵仍式（2-33）推出,2-33,对于等截面铰接杆单元，dv=Adx，A 单元截面面积。有,将式（5-8）代入上式，得,5-12,5-8,2、平面铰结杆单元,1）单元坐标单元位移向量,2）位移模式和形函数,由于平面铰接杆单元只有轴向力。位移模式同式（5-3）、（5-4,形函数,5-

9、13,3）应变矩阵,位移模式,5-7,应变矩阵 B为,5-14,4）应力矩阵,5-9,应力矩阵 S为,5-15,5) 等价节点力,5-16,6) 单元坐标单元刚度矩阵,对于等截面铰结杆单元,3、空间铰接杆单元,1）单元坐标单元位移向量,图5-6,5-18,2）形函数,5-19,3）应变矩阵,5-20,4）应力矩阵,5-21,5) 等价节点力,5-22,6) 单元坐标单元刚度矩阵,对于等截面铰接杆单元,5.4 梁单元,1、两端承受剪力、弯矩的平面梁单元,图5-7,l,1）单元坐标下的单元位移和单元力,单元位移,5-24,其中,vy方向位移，即挠度。 角位移。 对于梁： =dv/dx （5-25,

10、单元力,5-26,其中， Q剪力 M弯矩 对于梁,2）位移函数和形函数,设单元坐标位移模式为,5-28,位移模式,形函数,由单元两端点的节点位移条件，解出式（5-28）中的a1、a2、a3、a4。再代入该式，可将位移模式写为以下形式,5-27,5-29,式中,5-30,5-31,3）应变矩阵,单元弯曲应变b与节点位移e的关系。 由材料力学知，梁单元上任一点的应变和该点挠度之间关系为,5-32,将式（5-29）代入（5-32），得单元弯曲应变和单元位移之间关系,5-34,5-33,4）应力矩阵,5-35,5) 等价节点力,对于梁上作用的集中力或集中力矩，在划分单元时可将其作用点取为结点，按结构的

11、节点载荷处理。 这里考虑的是把单元上的横向分布载荷转化为等价节点力问题。当梁单元上作用有横向分布荷载qy(x)时（图5-8,图5-9,qy (x)dx,横向分布荷载qy(x)的势能Vq为,5-36,形函数矩阵由式（5-30）和（5-31）给出。对于具体问题，只要将qy(x)代入上式进行积分即可。表1给出了几种特殊情况的等价节点力,q,几种横向分布荷载等价节点力 表 1,6) 单元坐标单元刚度矩阵,梁单元刚度矩阵公式为,将式（5-34）代入上式进行积分，并注意到,Iz梁截面对Z轴（主轴）的惯性矩,得单元坐标单元刚度矩阵ke,5-37,单元刚度矩阵式(5-38)适合于连续梁分析,2、两端承受轴力、

12、剪力、弯矩的平面梁单元,图5-10,1）单元坐标单元位移和单元力,单元位移,5-39,其中,ux方向（轴向）位移。 vy方向位移，即挠度。 角位移,单元力,5-40,其中： N轴向力 Q剪力 M弯矩,对于小变形问题，可以认为轴向变形和弯曲变形互不影响，因此，位移模式和形函数可以分别按5.3节一维铰结杆单元和5.4节两端承受剪力、弯矩的平面梁单元的结果（式5-3和式5-28）简单集合而成,2）位移函数和形函数,位移模式,5-41,以下形函数和一些基本矩阵都可按此思路推演,形函数,式中形函数N为,5-42,5-43,其中,3）应变矩阵,单元弯曲应变与节点位移e的关系 承受轴向力、剪力、弯矩的梁单元

13、上任一点的应变，应为该点挠度（v）引起的应变和轴向位移（u）引起的应变之和,考虑到式（5-8）和（5-34），单元应变矩阵为,5-44,5) 等价节点力,i,j,图5-11,4）应力矩阵,5-46,应力矩阵形式同式（5-35,将式（5-36）、（5-11）膨胀成61矩阵后相加，并注意到式（5-43），有,5-11,最后得等价节点力矩阵,5-47,表2给出了几种特殊情况的等价节点力,几种横向分布荷载等价节点力 表 2,6) 单元坐标单元刚度矩阵,梁单元刚度矩阵公式为,式(5-48)用于分析平面框架,3、两端承受扭矩和面外剪力、弯矩的平面梁单元,图5-12,此类单元适用于格栅以及受面外荷载的平面框

14、架, 之所以仍称为平面梁单元，是由于结构本身是平面结构，而节点也是3个自由度,x、Mx截面绕扭心轴的扭转角和相应扭矩。 y、My截面绕y轴的弯曲转角和相应弯矩。 w、Q 截面形心的横向位移和相应横向剪力,如果截面形心和扭心不重合，则弯曲和扭转之间是相互不独立的。例如，横向剪力不通过扭心，它会引起对扭心轴的扭矩,这里只讨论截面形心与扭心重合或可以近似认为重合的情形。则弯曲和扭转之间是相互独立的,此外，这里的扭转限于纯扭转或称均匀扭转。其特点是扭矩和扭率（单位长度上的相对扭转角）成正比。即,5-49,扭矩平衡条件,5-50,由此得,5-51,式中GJ为截面扭转刚度,只需要将式（5-48）中的Iz换

15、成Iy，并注意编号次序。同时考虑到式（5-51），即得,4、空间梁单元,空间梁单元，每个节点有6个自由度，单元自由度为12。图5-13给出了空间梁单元节点位移分量的正方向及其编号。单元力的正向及其编号与单元位移相同,图5-13a,图5-13b,综合前述结果，得空间梁单元单元坐标单元刚度矩阵(式5-53）（式5-57,5-54,5-55,5-56,5-57,5.5 单元特性在两类坐标系中的转换,在5.3、5.4节中，单元位移和单元力都是按单元坐标系的坐标轴分量定义的，由此建立的单元刚度矩阵属于单元坐标单元刚度矩阵,进行系统分析时，需要把单元力按统一的结构坐标轴的分量表示出来，以便建立节点平衡方程

16、。因此，在进行系统分析之前，必须把单元坐标系中的单元力以及单元刚度矩阵都转换到结构坐标系中去。此外，还需要把结构坐标系中的节点位移转换到单元坐标系中去，以计算结构内力。因此了解坐标变换是必要的,设XYZ为结构坐标（整体）系，xyz为单元（局部）坐标系,约定下列符号,结构坐标单元位移 F结构坐标单元力 k结构坐标单元刚度矩阵,1、坐标变换矩阵定义,把单元位移从结构坐标系转换到单元坐标系的变换矩阵定义为坐标变换矩阵，用符号R表示。有,继续使用前单元坐标中的符号,e单元坐标单元位移 Fe单元坐标单元力 ke单元坐标单元刚度矩阵,式（5-58）给出了结构坐标单元位移转换为单元坐标单元位移的转换式，同时

17、是坐标变换矩阵R的定义式。 本节只了解R的存在和概念，有关坐标变换矩阵R 的具体形式、内容留在5.6节中专门讨论,2、结构坐标单元力,单元力在单元位移上作的功，不因其坐标系的改变而变。则有,5-58,将式（5-58）代入,对上式两端进行转置，注意到,消去，得,即得,5-59,式（5-59）表明：结构坐标单元力等于单元坐标单元力前乘坐标变换矩阵的转置,必须指出：式（5-58）是从整体（结构）坐标系到局部（单元）坐标系的变换式；式（5-59）是从局部（单元）坐标系到整体（结构）坐标系的变换式,在单元坐标系中，有,3、结构坐标单元刚度矩阵,上式两端左乘RT,注意到式（5-58）、（5-59），有,k

18、结构坐标单元刚度矩阵。并有,5-60,式（5-60）给出了把单元坐标单元刚度矩阵转换为结构坐标单元刚度矩阵的转换式,5.6 坐标变换矩阵,坐标变换矩阵因单元类型不同而异,1、平面铰接杆单元,设OXY为结构坐标，oxy为单元坐标。为任意单元 i 端的任一矢量。它在结构坐标系中的分量为X、Y；在单元坐标系中的分量为x、y。 X、Y 在单元坐标x轴上投影的代数和给出x 。同理， X、Y 在单元坐标y轴上投影的代数和给出y,由图5-14得,图5-14,5-61,写成矩阵形式,取,得,5-58,上式可写成,坐标变换矩阵R的具体内容为,用节点坐标描述方向余弦,5-62a,式中，（Xi，Yi）和（Xj，Yj）分别为节点i和节点j在结构坐标系中的坐标值,2、空间梁单元,Displacement components at node 1,Displace