安徽建筑工程学院计算结构力学1.ppt

,*,首页,上页,返回,下页,计算结构力学,授课人：干洪,陶明霞,韦璐,第一章 绪 论,1-1,概 述,结构矩阵分析,利用矩阵代数理论来分析结构力学问题，是随着计算机的迅速发展而兴起的结构分析方法。,计算结构力学,利用计算机来进行结构的力学分析。,计算结构力学的开课目的,本课程属于技术基础课，主要是强化计算机在结构分析方面的应用，是现代结构分析重要的不可缺少的手段，是专业技术能适应现代化需要的组成部分。,本课程主要研究杆系结构，可对六种杆系结构进行分析。,主要采用矩阵位移法或称杆系有限元法进行分析。,主要内容为建立结构刚度方程的矩阵形式及求解，程序设计和上机。,要求掌握杆系有限元进行结构分析的过程；,进行程序设计及计算机应用方面的训练；,要求在其它专业课中融汇贯通，借此达到专业技能与全面素质的提高。,本课程的主要内容和任务,1,、矩阵位移法,(,刚度法,),：以结点位移为基本未知量，建立结构的刚度方程。,2,、矩阵力法,(,柔度法,),：以结点力为基本未知量，建立结构的柔度方程。,3,、矩阵混合法,(,杂交法,),：以部分结点位移、部分结点力为未知量，建立结构的混合法方程。,结构矩阵分析的主要方法,由结构力学内容可知,:,刚度法只需满足平衡条件，在荷载形式一定的情况下自然满足，故普遍得到使用。,柔度法要确立多余约束建立基本结构，并满足位移,协调条件,，要具体分析，故很难规范化统一格式编程，不易实现计算自动化。,所以，工程计算一般采用,矩阵位移法,。,但在梁、板、壳等问题中，所假设的位移场在某些情况下不能满足一些单元的协调性,(C,连续性问题,),，故混合法或柔度法仍得到运用，并能进一步发展，现主要在板壳结构中使用。,本课程主要介绍,矩阵位移法,。,在矩阵位移法中,所有的方程组均采用矩阵的形式表示。,所有的推导和运算均借助于矩阵代数，形式紧凑明了，方便程序设计。,采用矩阵结构分析方法，并不改变结构力学的基本原理和基本假设。如平衡原理、叠加原理、变形协调原理、能量原理等。,本,课程基本假设：,小变形假设；,材料线性行为假设,(,结构联接为理想联结,),。,满足以上两个假设的结构称为线性结构。,1-2,有限单元法简介,结构理想化的概念：,结构理想化是一种简化手段，如同材料力学中的计算简图的概念。,在结构力学中，就是假设结构为连续体，理想连接、均匀各向同性的线性结构。经上述理想化以后，即可画出结构的计算简图，其主要特点有：,1,、以杆件轴线代替实际杆线；,2,、结构联结主要有刚结、铰结、链杆联接等；,3,、支座可简化为活动支座、固定铰支座和固定支座等。,结构矩阵分析所采用的主要方法为,有限单元法,，其基本思想是：把整个结构看成是由有限个单元,(,杆件、平面、壳体、块体等,),所组成的集合体，各个单元由结点相互连结，这就是结构的离散化，由各单元的平衡条件建立单元刚度方程，再利用整体平衡条件将各单元集合在一起，恢复为原结构，得到结构整体平衡方程,(,结构刚度方程,),。,结构刚度方程,形式为线性代数方程组，利用矩阵代数和数值计算方法编制成计算机程序，上机求解未知量。由此可知有限单元法的中心思想是,一分一合,。由于单元的个数有限，故称其为,有限单元法,。,单元的类型主要有：,杆单元平面单元及板单元,壳单元块体单元,杆单元,平面单元及板单元,壳单元,块体单元,本课程主要研究杆系结构，称为,杆系有限元,。,由于采用结点位移为未知量，故称为,有限元位移法,。,在实施中，由单元的刚度方程，依各结点的集约条件，可直接形成结构刚度方程，其方法称为,直接刚度法,。,结构的离散化过程,本课程可以对六种杆系结构进行分析,：,梁 刚架 桁架 排架 框排架 刚铰混合结构,(,或梁桁组合结构,),。无论对哪一种结构，总可以假想地将它拆开，视为有限个单杆在其端点联结，可以自然剖分，亦可以细分，这些单杆称为,单元,，联结点就称为,结点,(,节点,),。,1,2,3,1,2,3,4,5,6,1,2,3,4,5,1,2,3,4,5,6,7,1,2,3,4,5,6,7,8,9,10,11,12,13,14,1,2,3,4,5,6,结构的离散化过程：,(a),梁,(b),刚架,(c),桁架,(d),排架,(e),框排架,(f),梁桁组合,11,12,13,14,15,16,17,18,1-3,结点位移和结点力,结点位移包括：线位移和角位移。,单元两端的结点位移又称单元的杆端位移，或称其为单元结点位移。,已知杆端位移及荷载情况便可了解整个单元的变形状态，用,i,表示结点,i,的结点位移列阵；,表示单元的结点位移列阵；,表示结构的结点位移列阵。,桁架：,连续梁：,刚架：,结点力包括：力和力偶矩。,单元结点力：单元杆端力，这是结构内力，对单元而言是作用在单元两端结点上的外力。,结构结点力：由于汇交于每一结点的各单元杆端力的总和即等于该结点所受的力，故结构结点力是外力，为相应的结点荷载或结点支座的支座反力。,结点,i,的结点力列阵用,F,i,表示；,单元的杆端力列阵用,F,表示；,结构结点力列阵用,P,表示；,反力用,R,表示；,结点位移与结点力的各个分量应相互对应，如：,i,与,F,i,与,P,；,结点位移编号,(,或结点力编号,),与结点编号有关。结点编号是人为的，现已可用程序实现结点自动编号；,在进行结构分析时，首先应编好结点号。结点编号的好坏直接影响计算精度及内存，其原理是应尽量使每个单元两端结点号的差值最小。,1-4,基本未知量,如何确定结构的基本未知量,?,根据有限单元法的离散化要求，各个单元仅在结点处联结，因此只有结点处的力学量,(,结点位移或结点力,),可以作为基本未知量。,由于矩阵位移法采用位移为未知量，故在有限元位移法中采用结点位移作为未知量。,因为荷载已知，此时相应的结点力向量应为已知，这对于一般结点均满足。,关于支座情况，需要进行,约束处理。,因为结点力包括了支座结点反力，这在通常情况下为未知，这点看来与上述要求不符。但由于在不考虑弹性支承情况下，有结点反力的这部分支座位移通常已知,(,零或已知沉降量,),，不需求解，可在结构刚度方程中将这一方程划去，这就是,约束处理。,直接刚度法分为,前处理法,及,后处理法,。,前处理法,在形成结构刚度矩阵之前，也就是在建立结构刚度矩阵之前考虑到实际的约束情况，再形成结构刚度矩阵。,后处理法,在形成结构刚度矩阵之前，先不考虑支承情况，而在形成结构刚度矩阵之后，再根据约束情况对结构刚度矩阵进行修改。,很明显，前处理法可减少存贮单元。,本课程采用前处理法，具体的,约束处理,如下：,对于一般结点,(,指无约束的结点,),的未知量编号，在对结构的全部结点编号后即可确定；,对于整个结构的全部未知量编号，还需要加上支座的未知量编号；,这实际上就是要对结构的支座约束进行处理，这是结构分析中非常重要的环节。,例,1,：对图示一般刚架进行结点和未知量编号。,解：先编结点号，后根据每个结点三个未知量进行未知量编号。,对支座进行约束处理，可通过约束特征数来实现。,支座约束特征数：,表示支座结点的某一位移未知量有无刚性约束的人为赋值数字。,有约束：约束特征数为,1,；,无约束：约束特征数为,0,。,显然，有约束则无位移未知量，无约束则有位移未知量。对于采用刚架单元的平面杆系结构，其,支座约束特征数,如下：,例,2,：对图示结构进行结点未知量编号。,解：结点编号如图所示，未知量编号见表,1,。,若不考虑杆件的轴向变形，则未知量编号见表,2,。,例,3,：刚铰混合结构,解：因铰处有两个转角未知量，因此多编一个结点号！（见表,1,）,若不考虑杆件的轴向变形，则未知量编号见表,2,。,2,例,4,：梁桁组合结构,1-5,轴力杆单元刚度方程,(,桁式单元,),首先介绍坐标系的概念。,1,、结构坐标系,(,整体坐标系,),；,2,、自身坐标系,(,局部坐标系,),。,局部坐标系建立在单元上，由始结点至终结点，用“,”,作记号。,由结构力学位移法可知，杆单元的平衡方程可写成,k=F,的形式，其中,k,为,22,阶的单元刚度矩阵，现予详细介绍。,单元刚度系数,k,ij,的定义：仅当：,时在,i,处所需施加的力。,取单元两端轴向位移为未知量：,u,1,、,u,2,如右图所示。,当,u,1,=1,、,u,2,=0,时，由刚度系数的定义可知：,F,11,=k,11,u,1,=k,11,F,21,=k,21,u,1,=k,21,当,u,1,=0,、,u,2,=1,时，,F,22,=k,22,u,2,=k,22,F,12,=k,12,u,2,=k,12,这里,也可由刚度系数的定义直接得到,当局部坐标系和整体坐标系不重合时，如竖杆、斜杆等，其刚度方程（或刚度矩阵）如何推得？一般可通过坐标变换（后面再详细介绍），现再利用静力法推导杆单元在整体坐标系中的刚度方程。,此时，杆单元每端有两个未知量：,而,u,1,、,v,1,所引起的杆端力可由类似方法得到，令,cos,=l,sin=m,，,最后可将单元的刚度方程写成：,可将上面的单元刚度方程写成矩阵形式：,1-6,刚度法的基本概念,采用结构力学中位移法的基本方法，建立结构的刚度方程，现以下面的超静定桁架为例来说明。,例：求图示结构各杆的轴力，,A=2000mm,2,E=200KN/mm,2,。,解：为未知量，刚度法基本步骤：,1,、列出结点的平衡方程（如有多个结点，均应依次列出）,2,、变形条件（用位移来表示各单元变形）,如图，,e,i,在表示,i,的伸长量,3,、物理量,线刚度,(*),最后结构刚度方程（）是由平衡、变形、物理条件形成，若能由计算机直接形成（），即这三步均可不列，过程得以简化，关键在于求,K,、,P;,求,K,的主要过程在于第,3,步，相当于单元分析，如果能求得各类单元的统一公式，求,K,则只是程序设计而已。,上述二点为刚度法的主要点，由此可设想直接刚度法的具体过程。,讨 论,1-7,坐标系与单元定位向量,确定结构上各点,(,特别是结点,),的位置及几何参数，坐标系类型主要有：,1,、整体坐标系,(,结构坐标系、固定坐标系,),。,2,、局部坐标系,(,自身坐标系、活动坐标系,),。,一、建立坐标系的目的,说明：,图示中所示方向均为正向；,结构坐标的原点可任意设置，通常设在左下角；,单元坐标系的原点规定设在单元的始结点，由此可确定,y,轴的正向，及转角的正方向。,二、单元定位向量,定 义：,它,是按单元结点编号顺序由单元各结点的未知量编号所组成的一列数字,(,列向量,),，其分量即为该单元结点的未知量编号。,1,、直接确定各单元刚度矩阵在结构刚度矩阵中的位置，由此可极其方便且准确无误地形成经约束处理后的,K,；,2,、可将单元的等效结点力,F,e,叠加到,P,中；,3,、可,从,中取出,；,4,、解决一维变带宽存贮的寻址问题；,5,、边界或约束处理；,6,、利用主从关系可模拟各种类型的结构等。,单元定位向量的主要用途,(,在直接刚度法中,),：,可以说，程序设计的自始至终，在每一个环节都离不开,单元定位向量,，计算机自动化计算所需的信息主要由它提供，所以在程序设计中是各模块的组织者，起到主线的作用，占有非常重要的地位，值得深入研究。,对每个单元,，,只要知道,i,、,j,及两端结点的未知量编号，即可确定每个单元的定位向量,MW(6),。,e,e,例：求图示刚架各单元的定位向量,MW(6),。,解：结点及单元编号如图，单元定位向量列表如下：,1-8,形成结构单元定位向量的程序设计,采用平面刚架单元进行分析，首先介绍几个概念：,1,、特殊结点,凡未知量不足三个的结点称特殊结点。,介绍一些常用的数组和变量：,NJ,：,结点,(,总,),数,NJT,：,特殊结点数,2,、约束特征数,表示结点的某一位移方向有无约束或该位移是否是无效未知量的人为赋值数字。,1,：表示有约束，沿该方向位移为零。,0,：表示无约束，沿该方向可以有自由位移。,10001,：表示无效的未知量。,3,、主从关系,表示非独立的结点位移未知量与主结点位移之间的从属关系。,主从关系由约束特征数表示,：即从结点相应未知量的约束特征数可直接填主结点号，当主结点为,1,时，以,1001,代替。,4,、结点未知量编号数组,JW(3,，,NJ),这是形成单元定位向量所必须的数组，要形成,JW(3,，,NJ),，,则要了解结点约束情况，杆系类型,(,无效未知量情况,),，主从结点情况，故应先输入特殊结点信息，以,JTX(4,，,NJT),表示。,以上过程简称,约束处理,，需事先在程序中作为输入语句。,例,1,形成图示刚架结,构的,JTX(4,NJT),。,解：本例特殊结点数,NJT=2,，,特殊结点信,息数组,JTX(4,2),为,凡从属结点都应作为特殊结点处理，要填约束信息，,允许一主多从，而不允许一从多主。,JTX(4,，,NJT),的填写是项很认真的工作。编写好这个数组，就可以用一般平面刚架单元的程序设计来分析各种类型的杆系结构。,注 意：,1,2,3,4,例,1,连续梁,解：本例,NJT=4,例2,桁架,解：本例,NJT=4,1,2,3,4,例3,刚铰混合结构分二种情况,EA=const；EA=,EA=，,NJT=6,解：,EA=const，NJT=3,C,PROGRAM OF MWEXAM,DIMENSION JH(2,20),JTX(4,20),JW(3,20),MW(6),WRITE(*,*)FINDING THE MW(6)OF,*,ELEMENTS,OPEN(1,FILE=MWE.DAT),READ(1,*)NE,NJ,NJT,READ(l,*)(JH(I,J),I=1,2),J=1,NE),READ(1,*),(JTX(I,J,),I=1,4),J=1,NJT),CALL QJW(NJ,NJT,JTX,JW,N),形成单元定位向量的程序设计,DO 10 M=1,NE,CALL QMW(M,NE,NJ,JH,JW,MW),10 WRITE(*,100)M,(MW(I),I=1,6),100 FORMAT(1X,ELEMENT No.=,I5,/6X,MW,*,=,615),STOP,END,SUBROUTINE QMW(M,NE,NJ,JH,JW,MW),DIMENSION JH(2,NE),MW(6),JW(3,NJ),JL=JH(1,M),JR=JH(2,M),DO 100 J=l,3,MW(J)=JW(J,JL),100 MW(J+3)=JW(J,JR),RETURN,END,子程序,QJW,见书,14,页,子程序,QMW,为：,数据文件为：,9,8,2,3,1,4,2,5,3,6,4,7,5,8,6,1,2,3,4,5,6,7,1,1,1,8,1,1,1,习题,2,：就本节三个例题，分别写出相应结构 的单元定位向量,上机计算并打印结果,例,4,：编写本节例,1,三层刚架求,MW,的数据文件。,