1、有限元法基础与ANSYS简介,通过介绍有限元法的基本概念、理论、方法与软件，使学生能够掌握使用其求解力学问题的特点、解题过程，熟悉一种有限元软件，初步具备使用有限元方法解决工程设计中实际问题的分析能力。,本课程讲授目的,提 纲,1 绪论 2 有限元法的基本概念与求解方法 3 有限元法常用单元介绍 4 ANSYS软件介绍与基本使用方法 问题与讨论 附录A：弹性力学中的基本概念与方法 附录B：平面问题的基本理论,1.1 有限元法的一般概念 1.2 有限元法的发展简介 1.3 有限元法与其他课程之间的关系,绪 论,有限元法是求解数理方程的一种数值计算方法，是解决工程实际问题的一种有力的数值计算工具，
2、最初这种方法被用来研究复杂的飞机结构中的应力，是将弹性理论，计算数学和计算机软件有机的结合在一起的一种数值分析技术。由于这一方法的灵活，快速和有效性，使其迅速发展成为求解各领域的数理方程的一种通用的近似计算方法，目前已在许多学科领域和工程问题中得到广泛的应用。 常用数值分析方法：差分法，有限元法，边界元法,有限元法的一般概念,有限元法的一般概念,将一个连续的求解域（连续体）离散化即分割成彼此用节点（离散点）互相联系的有限个单元，在单元体内假设近似解的模式，用有限个结点上的未知参数表征单元的特性，然后用适当的方法，将各个单元的关系式组合成包含这些未知参数的代数方程，得出个结点的未知参数，再利用插
3、值函数求出近似解。是一种有限的单元离散某连续体然后进行求解得一种数值计算的近似方法。由于单元可以被分割各种形状和大小不同的尺寸，所以它能很好的适应复杂的几何形状，复杂的材料特性和复杂的边界条件，再加上它有成熟的大型软件系统支持，使它已成为一种非常受欢迎的，应用极广的数值计算方法。,有限元法的一般概念,单元：分割连续体的小区域，有线、面或实体等种类。 节点：连接单元的空间点（由空间坐标确定），具有一定自由度。 自由度：用于描述一个物理场（位移）的响应特性的参量。,有限元法的一般概念,J,有限元法的一般概念,节点编号节点编号分为局部节点编号和总体节点编号两种，如下图中的矩形，分为个节点，个单元，其
4、中，为总体节点编号。而对于任一单元 中，为局部节点编号，在公式推导中使用i，j，m编号。单元编号按从小到大顺序依次排列。,有限元法的发展简介,50年代，发展与萌生，单一功能程序，简单单元； 60年代，数学基础与证明，单一功能程序，多种单元； 70年代，单元库丰富，线性到非线性通用程序，如SAP，NONSAP等； 80年代，多种功能扩大，大型通用程序如ADINA等； 90年代，领域扩大，前后处理功能增强，大型商用软件，如ANSYS、MARC、NASTRAN等； 目前，面向工程，与CAD结合成为CAE（计算机辅助工程）软件。 说明：有限元法的发展与计算机学科的发展紧密相关。,有限元法的发展简介,现
5、有有限元软件的计算功能与应用可求解结构位移场、温度场、电磁场、流场、耦合场等多种问题。,有限元法的发展简介,现有有限元软件的计算功能与应用可求解结构位移场、温度场、电磁场、流场、耦合场等多种问题。,有限元法与其他课程的关系,有限元法与其他课程的关系,材料力学：研究杆状构件在拉压、剪切、弯曲、扭转作用下的应力和位移。其中引入了构件形变状态或应力分布的假设 。 结构力学：在材料力学基础上研究杆状构件所组成的结构（杆件系统），例如行架等的应力与位移。 弹性力学：针对非杆状结构（如板和水坝等实体结构）以及对杆状构件作进一步较精确的分析。以微元体为研究对象，通过建立应力、形变与位移间的关系进行求解。 计
6、算力学：是结构力学、弹性力学、计算数学、计算机学的结合，提供近似的数值计算方法解决问题。有限元法是其中的一种方法。上述各种方法最终目标是确立研究对象的应力、形变和位移，用以校核其是否有所需要的强度和刚度。,有限元法的基本概念与求解方法,2.1 结构离散化与刚度矩阵 2.2 位移函数与形函数 2.3 单元刚度方程 2.4 载荷移置与等效节点载荷 2.5 结构刚度方程 2.6 位移边界条件的处理 2.7 应力计算 2.8 有限元法的普遍公式 2.9 有限元方程组的解法,结构离散化与刚度矩阵,结构离散化：1）网格划分将结构划分为有限个单元；2）载荷移置将作用在结构上的非节点载荷等效地移置为节点载荷；
7、3）简化约束把结构边界上的约束，用适当的节点约束代替。,结构离散化与刚度矩阵,有限元网格划分原则,有限元中单元的网格剖分原则）各节点必须相连。如图所示中（a）是正确的，而（b）是错误的。,结构离散化与刚度矩阵,有限元网格划分原则,结构离散化与刚度矩阵,有限元网格划分原则,）单元不能奇异，也就是单元中的边长不能相差太大，或者有过大的钝角或过小的锐角，如图示：,结构离散化与刚度矩阵,有限元网格划分原则,）单元的大小、数目取决于计算精度要求和计算容量限制分网时首先满足计算精度的要求，同时可利用结构的对称性、循环对称性的特点，从厚结构中取出一部分进行分析，或者对有应力集中的构件，采用疏密不同的网格剖分
8、。也可以采用子结构法。 ）同一单元内的结构，几何特性与材料特性相同，也就是不要把厚度不同或材料不同的区域划分在同一个单元里。,结构离散化与刚度矩阵,网格划分示例,结构离散化与刚度矩阵,结构离散化与刚度矩阵,刚度矩阵描述单元特性的矩阵，表示了单元抵抗变形的能力。它由刚度系数组成，由单元节点的个数和自由数决定规模。如图平面三角形三节点单元中，有3个节点，每个节点有2个自由度，故刚阵中的元素个数为36个。 刚度系数Kij相当于一维弹簧的刚度K的含义。即产生单位位移时需要的作用力的大小。,结构离散化与刚度矩阵,位移函数,结构离散化后，要对单元进行力学特性分析，也就是确定单元 节点力与节点位移之间的关系
9、，这时就需要把单元内的任一点的位 移分量表示成坐标的某种函数。这种函数就叫位移函数。,位移函数与形函数,位移函数的一般介绍,.定义：把单元中任一点的位移分量与坐标的函数关系叫位移函数或叫位移模式。 .选择位移函数的原因 （）决定了单元的力学特性。（意义） （）反映了单元的位移形态。（物理意义） （）它是利用位移法求解问题的开始。（基础） .位移函数必须具备的条件 （）在节点上的值应等于节点的位移 （）所采用的函数必须保证有限元的解收敛于真实解,位移函数与形函数,位移函数的一般形式,位移函数一般为多项式形式，这样处理是从两方面出发的（）进行数学运算（如微分，积分）较简单（）任意阶次的多项式可以近
10、似地表示精确解，其一般形式为：u=u(x,y)=1+ 2x+ 3y + 4x2+ 5xy + 6y2 + + mynv= v(x,y)=m+1+ m+2x+ + 2myn （-） 式中: ,其中1 m为待定系数。 式中的也称为广义坐标，这种描述方式又称为广义坐标形式。（一维形式多项式u(x)=1+ 2x+ x2+ + nxn）,位移函数与形函数,位移函数的一般形式,（-）式也可以参照帕斯卡三角形来确定,位移函数与形函数,三节点三角形单元的位移函数,.位移函数形式就是最简单的情况而言，可以选取位移为坐标的线性函数形式，也就是：u(x,y)=1+ 2x+ 3yv(x,y)=4+ 5x+ 6y （-
11、）对于图中的三角形单元，为了确定（-）式中的待定系数16，可以将节点i，j，m的位移值及坐标值代入上式，得到方程组：ui=1+ 2xi+ 3yivi=4+ 5xi+ 6yi （ i=i，j，m） （-） 式中 ui ， vi节点位移 xi ， yi节点坐标,位移函数与形函数,三节点三角形单元的位移函数,这是一个一阶线性方程组，可使用克来姆法则求解。,位移函数与形函数,三节点三角形单元的位移函数,.克来姆法则设有一线性方程组：a11x1+ a12x2 + a1nxn =b1a21x1+ a22x2 + a2nxn =b2an1x1+ an2x2 + annxn =bn （a11 ann系数）当其
12、系数行列式 不等于零时上述的方程组有唯一解： （j，n）其中 是将中第j列元素替换为右端项而得到的行列式,位移函数与形函数,三节点三角形单元的位移函数,.待定系数1 的求解如果用节点位移（ui, vi）,（uj, vj）,（um, vm）及节点坐标 （xi, yi）,（xi, yi）,（xi, yi）代入（2-3）式可以得到：ui=1+ 2xi+ 3yiuj=1+ 2xj+ 3yjum=1+ 2xm+ 3ymvi=4+ 5xi+ 6yivj=4+ 5xj+ 6yjvm=4+ 5xm+ 6ym,位移函数与形函数,三节点三角形单元的位移函数,由克来姆法则可知：当2 0，上述方程有唯一解：,位移函数
13、与形函数,三节点三角形单元的位移函数,为了描述方便，引入系数ai= xj ym - xmyj bi= yj - ym ci= -xj + xm aj= xmyi - xi ym bj= ym - yi cj= -xm + xiam= xi yj - xj yj bm= yi - yj cm= -xi + xj,位移函数与形函数,三节点三角形单元的位移函数,代入上式后可以得到,位移函数与形函数,三节点三角形单元的位移函数,.位移函数的插值函数形式 假设这样一个函数：(i=i , j , m) 代入（-）式后可得u=iui +juj+mum v=ivi +jvj+mvm 式中：i，j ，m被称为单
14、元的形状函数，简称形函数或插值函 数。,位移函数与形函数,三节点三角形单元的位移函数,把（-）式写成矩阵形式：简写为：f= Ne （-）,位移函数与形函数,三节点三角形单元的位移函数,式中的矩阵反映了单元的位移形态，又是坐标的函数，我们 称之为形函数矩阵，这种描述方式称为位移函数的插值函数形式。通过上面的推导，我们得到了两种形式的位移函数，（-）式与（-）式后一种描述更简单，更直观，通常采用。这样我们就建立了单元中任一点的位移和单元节点位移之间的关系。,位移函数与形函数,位移函数及其性质,当节点位移一定时，单元形态完全决定于i，j ，m这时形 函数就具有如下的性质：.形函数i在节点i处的值为，
15、而在其他两个节点（j，m） 处的值为零。即： i （xi ，yi ）=1 而i （xj ，yj ）=i （xm ，ym ）=0 同样的j （xi ，yi ）=0 j （xj，yj ）=1 i （xm ，ym ）=0 m （xi ，yi ）=0 m （xj，yj ）=0 m （xm ，ym ）=1,位移函数与形函数,位移函数及其性质,.在单元任一节点处，三个形函数之和等于。 证明如下： i（x，y）+j（x，y）+m（x，y）= （ai+ bix+ ciy+ aj+ bjx+ cjy+am+ bmx+ cmy）（ ）=（ ai + aj + am ）+ （ bi + bj + bm ）x+ （
16、ci + cj + cm ）y（ ）=（ +0+0）/ （ ）=1 此外，形函数与位移函数是同样类型的函数。 如：位移函数 u=1+ 2x+ 3y 形函数 i=（ ai + bix + ciy）（ ）,位移函数与形函数,位移函数与解的收敛性,选择位移函数时，为保证有限元法的收敛性，必须满足以下 个条件： .位移函数必须包含单元的常量应变 .位移函数必须包含单元的刚体位移 .位移函数在单元内部必须是连续函数（连续性要求） .位移函数应使得相邻单元间的位移协调（保续性要求）上述四个条件中，若全部满足，这样的位移函数构成的单元称 为协调单元，若只满足前三条，则称为非协调单元,位移函数与形函数,位移函
17、数与解的收敛性,下面我们用以下四个条件来考察三角形常应变单元的位移函数 （）由=x , y , xyT =2 , 6 , 5 + 3 T 因2 , 6 , 5 + 3都是常数，与某坐标无关，因此含有常应变项 （）将位移函数可改写成,位移函数与形函数,位移函数与解的收敛性,当发生刚体位移时： x = x = xy =0 也就是2 = 6 = 5 + 3 = 0这时：其中u0 , v0为平动位移分量。0为单元绕垂直于x，y平面的轴线作刚体转动时的角位移，它 表示了刚体位移。,位移函数与形函数,位移函数与解的收敛性,（）位移函数（-）或是x，y的单值连续函数，故满足连续性 要求。 （）位移函数（-）
18、式是线性函数，由于相邻单元在公共节点 处的位移值相等，而通过两个节点可以连成一直线，其连线上的位 移相同，因此边界上各点的位移是连续的，不会出现：综上所述，三角形常应变单元属于协调元,位移函数与形函数,面积坐标,面积坐标是利用三角形的面积关系表示三角形单元任一点位置 的一种方法。优点：简明，方便。,位移函数与形函数,面积坐标,对于图中三角形单元任一点P（x , y）可用下三个比值来确定：Li , Lj , Lm 称为P点的面积坐标，显然面积坐标具有以下性质： 性质. Li + Lj + Lm =1 （ i+ j + m = ） 性质. 平行于三角形jm边的直线上所有点其Li相同AB变化时， h
19、i不变, 故i不变， Li 不变,位移函数与形函数,面积坐标,性质. Li =1 Lj =0 Lm =0 （i）Li =0 Lj =1 Lm =0 （j）Li =0 Lj =0 Lm =1 （m） 性质. Li = Lj = Lm=1/3 在三角形形心处 面积坐标与形函数的关系：,位移函数与形函数,面积坐标,同理： Lj = Nj ， Lm = Nm 所以，面积坐标与形函数相同（量值）但意义不同,位移函数与形函数,单元刚度方程,对单元进行力学特性分析目的在于确定单元节点力与节点位移的 关系，并称之为单元刚度方程：e e =e式中：Fe ， e 单元节点力及节点位移列阵e 单元刚度矩阵,基本方法
20、,单元刚度方程,基本方法,建立上述方程时可采用的方法 （）直接刚度法 （）虚位移原理或最小势能原理位移型有限元 （）余虚功原理或最小余能原理力型有限元 （）变分法（非结构问题）,单元刚度方程,基本方法,单元特性分析的步骤 （）假设位移函数 （）建立应力，应变与节点位移间的关系 （）由能量原理，建立单元节点力与节点位移间的关系 （）得到单元刚阵,单元刚度方程,三角形平面单元的单元刚度矩阵,（）上节的知识可以知道 位移函数为： u=iui +juj+mum v=ivi +jvj+mvm 式中 i=（ ai + bix + ciy）（ ） （i=i,j,m） （）应力应变与节点位移的关系 对三节点三
21、角形单元，节点位移e =ui , vi , uj , vj , um , vmTFe =Fix , Fiy , Fjx , Fjy , Fmx , FmyT,单元刚度方程,三角形平面单元的单元刚度矩阵,由弹力知识可知，几何方程为：,单元刚度方程,三角形平面单元的单元刚度矩阵,令：=Bi , Bj , Bm且 （i=i,j,m）方程可简写为： = e,单元刚度方程,三角形平面单元的单元刚度矩阵,我们称单元的几何矩阵，其物理意义反映了单元任一点 的应变与单元位移之间的关系。对于一个给定的单元，节点坐标一定，系数bi，ci也随之确定，也为常数，所以几何矩阵为常量矩阵，这也证明节点三角形单 元是一种常
22、应变单元。由弹性理论中关于平面问题的物理方程可知，当不考虑变温影 响时，单元中任一点的应力为: =D 式中为弹性矩阵，反映了单元材料方面的特性。,单元刚度方程,三角形平面单元的单元刚度矩阵,由上面应变与节点位移之间的关系代入后可得 = D = DBe 若令S= DB 则 = Se 式中，S称为单元的应力矩阵物理意义：反映了单元中任一点的应力与节点位移之间的关系， 对于节点三角形单元D，B为常量矩阵， S也为常量矩阵，这种常 应变单元，也是一种常应力单元，回顾一下，平面应力问题：,单元刚度方程,三角形平面单元的单元刚度矩阵,而对于平面应变问题如果采用： 代入,单元刚度方程,三角形平面单元的单元刚
23、度矩阵,两种问题具有相同的描述形式，只是对材料的弹性模量与泊松 比进行相应的代换，则在计算中可以采用同样形式的弹性矩阵。 （）单元节点力与节点位移之间的关系在位移型有限元法中，对单元的力学特性分析，最终是需要建 立节点位移和节点力之间的关系，也就是确定单元的刚度矩阵。应 用虚位移原理来建立这种关系式。设某单元发生一虚位移，则该单元各节点上的虚位移为e ， 相应地单元内任一点处的虚应变为： 。根据与间的关系有： =B e,单元刚度方程,三角形平面单元的单元刚度矩阵,这时单元体在节点力作用下处于平衡状态，根据虚位移原理， 当虚位移发生时节点力在虚位移上所做的功等于单元的虚应变能， 即：式中：e为单
24、元的体积，上式称为单元的虚功方程。 把 = DBe和 =B e代入上式得由于节点位移e及节点虚位移e均为常量，提出积分外，有：,单元刚度方程,三角形平面单元的单元刚度矩阵,进一步可得：令： 则上式可写为求得了我们所要的形式的方程，称之为单元刚度方程，式中的 e称为单元的刚度矩阵，反映了节点力与节点位移之间的关系。同样，可采用最小势能原理来建立单元节点力与节点位移的关 系式。我们得到的单元刚度矩阵e是普遍公式，适用于各种类型的单 元，对于三角形常应变单元的具体表达式见下。,单元刚度方程,三角形平面单元的单元刚度矩阵,（）三角形常应变单元刚度矩阵的显式：由于普遍公式中，均为常量矩阵，可以提出积分符
25、号， 而d是单元的微元体体积且d=t dx dy式中t为单元的厚度，同一单元，厚度t为常数，故单元体积（ 为单元的面积）普遍公式就可写为：为了便于计算利用B=Bi Bj Bm将上式展开,单元刚度方程,三角形平面单元的单元刚度矩阵,单元刚度方程,式中子刚阵为：Krs =tBrTDBs （r,s= i,j,m）,三角形平面单元的单元刚度矩阵,Krs是一个阶矩阵，因此三角形常应变单元的刚度方程为 的方程，也就是单刚阶数单元的自由度数。对与平面应力问题：将： B=Bi Bj Bm及 代入,单元刚度方程,三角形平面单元的单元刚度矩阵,（r,s= i,j,m）,单元刚度方程,三角形平面单元的单元刚度矩阵,
26、简写为：相应的：,单元刚度方程,单元刚度矩阵的性质,（）单元刚度矩阵是对称矩阵 （）单元刚度矩阵的主对角元素恒为正值 （）单刚为奇异阵 （）单元刚度仅与单元的几何特性（）及材料特性有关（） 而与外力无关。上述四条性质，与杆系的单刚性质相同,单元刚度方程,.由于在进行有限元分析中，单元和单元之间仅通过节点相互联系 当外载不是直接作用在节点上，那么需要将非节点载荷向节点移置， 也就是真实外载 （理想化） 节点上的集中载荷 移置后的载荷称之为等效节点载荷。,非节点载荷移置,载荷移置与等效节点载荷,非节点载荷移置,.结构的非节点载荷移置将各单元所受的非节点外载荷分别移置到各单元的相应节点上， 在公共节
27、点处应用载荷叠加原理，就可以得出 .载荷移置的原则能量等效的原则单元的实际载荷与移置后的等效节点载荷在相应的虚位移上所 做的虚功相等。 .单元载荷移置的方法 （）直接法：利用能量等效原则，直接进行单元载荷移置只适用于线性位移函数的单元,载荷移置与等效节点载荷,非节点载荷移置,（）普遍公式法：根据能量等效原则，推导出普遍公式适用于各种类型的单元 说明：由圣维南原理可知，载荷移置后，只会在结构的局部产生误 差。对整个结构的变形或应力状态的影响不大，由于有限元分析中， 单元一般都很小，移置的结果不会带来很大的误差。,载荷移置与等效节点载荷,载荷移置的普遍公式,.集中力的移置公式： 设（）单元i,j,
28、m中任意一点（x，y）作用集中载荷P=Px，Px（）各节点上的等效节点载荷向量为：Re=Rix, Riy, Rjx, Rjy, Rmx, Rmy（）发生微小位移时，集中力作用点相应的虚位移为：f*=u, v（）各节点相应的虚位移为：*e =ui* , vi* , uj* , vj* , um* , vm* ,载荷移置与等效节点载荷,载荷移置的普遍公式,推导：（）根据单元内位移与节点位移关系 f=u, v=N e f*= N *e P（）根据能量等效原则： *e Re =f*P *e Re =(N *e ) P = *e NP Re = NP 这就是集中力的移置公式，式中为单元的形函数矩阵。,载
29、荷移置与等效节点载荷,载荷移置的普遍公式,.体力g的移置公式 设：单元ijm上作用有体力 g=gx, gy 推导（）将单元体tdxdy上的体积力gdxdy当作集中力，应用集中 力的移置公式 Re = NP 有微元体上dRe = N g tdxdy（）积分在整个单元上有Re = N g tdxdy = t N g dxdy,载荷移置与等效节点载荷,载荷移置的普遍公式,3.表面力q的移置公式 设单元ijm的jm边上作用有表面力q=qx，qy 可将微元面积上tds上的面力qtds当作集中力 则 Re = sjmdRe = sjm qt ds = tsjm q ds 上述公式适用于任何单元及任意坐标方
30、向。,载荷移置与等效节点载荷,载荷移置举例,以单元自重（或作用在单元形心处的集中力）为例。设一个均质等厚的三角形单元ijm，其厚度为t，面积为，材料 比重为，则单元的自重为：W=t,且其作用在单元形心c处,载荷移置与等效节点载荷,载荷移置举例,思路：（）欲求哪个节点在哪个方向上的载荷分量，就在该方向加 一单位虚位移，其他自由度为。（）利用线性位移函数的特点导出几何关系。（）根据能量等效原则列出虚功相等，解出节点载荷分量,载荷移置与等效节点载荷,载荷移置举例,.直接法求解先求，i在y方向的等效节点载荷设 vi* =1 而ui* =uj* = vj* = um* = vm* =0相当于上图 由于，
31、单元具有线性位移函数，当vi* =1时变形情况见上图jm边不动，点b亦不动，由几何关系可知：,载荷移置与等效节点载荷,载荷移置举例,又根据能量等效原则：式中“”号表示与y轴方向相反 同理可得：,载荷移置与等效节点载荷,载荷移置举例,类似可以得出：各节点沿x方向的等效节点载荷 Rix= Rjx= Rmx= 0 Re=Rix , Riy , Rjx , Rjy , Rmx , RmyT= - t 0 1 0 1 0 1 T/3 上式也表明对三角形单元（均厚，等厚）所受重力，只需将自重平 均的移置到节点上，方向与重力方向相同。,载荷移置与等效节点载荷,载荷移置举例,.普遍公式法求解对于此处为集中力P
32、=0 - W T作用在形心c处由Re=N T P 可知,载荷移置与等效节点载荷,载荷移置举例,可以证明：在三角形形心c处，有Ni= Nj= Nm= 1/3代入上式可得Re =-W0 1 0 1 0 1 T/3= - t 0 1 0 1 0 1 T/3 可见采用上述两种方法移置的结果相同 说明：单元具有线性位移函数时，采用直接法移置较简单单元具有非线性位移函数时，只能采用普遍公式法进行。,载荷移置与等效节点载荷,结构刚度方程,通过单元特性分析，可建立单元刚度矩阵Ke同时得到单元刚度方程 Kee = Fe通过单元载荷移置，可建立节点载荷列阵Re集合成结构刚度方程的三个方面的内容是：（）单元的节点位
33、移e 结构的节点位移列阵（）单元的节点载荷列阵Re 结构的节点载荷列阵R（）单元的单刚Ke 结构的总刚K得到 K = R （2-45）,结构刚度方程,结构刚度方程,上为结构刚度方程，表示了节点载荷与节点位移间的关系，是一个以节点位移为未知量的线形代数方程组，可求得 ，进一步求出应变，应力。,结构刚度方程,集合的基本原则,（）在相互连接的公共节点处，各单元的节点位移必须 相等，即必须满足变形协调条件。i =i =i =i 所以，节点位移不须按单元来区分。,结构刚度方程,集合的基本原则,（）公共节点处，各单元对节点的作用力，与作用在该节点上的外载荷Ri之间，必须满足静力平衡条件。Ri =Fi +F
34、i +Fi +Fi 所以，若Ri = 0，则有Fi +Fi +Fi +Fi = 0,结构刚度方程,结构刚度方程的建立,例：1.结构的节点位移列阵根据公共节点处的变形协调条件，不同单元在公共节点处的位 移相等，则有节点位移列阵 =1 2 3 4 T=u1 v1 u2 v2 u3 v3 u4 v4 T （按总体节点编号顺序写出）,结构刚度方程,结构刚度方程的建立,2.结构的节点载荷列阵：（1）若存在非节点载荷，须进行单元载荷移置，并按移置后 的等效节点载荷进行叠加即： Ri=Rix Riy T =Ri + + + （2）不考虑约束反力的作用（3）与节点位移相对应，结构的节点载荷列阵R亦按总体节点
35、编号顺序排列 那么，对于上例R =R1 R2 R3 R4 T=0 0 0.5P 0 0.5P 0 0 0 T,结构刚度方程,结构刚度方程的建立,式中约束反力 R1 = R4 = 0 0 0 0 T 3.结构刚度方程3节点三角形单元的自由度数为6，单刚Ke为66阶矩阵。 总体K由单刚Ke组合而成 总刚K的阶数=结构的自由度数对于图示的例子，4个节点，共8个自由度，结构刚度矩阵为 8 8的方阵。 把图中的两个单元离散开为：,结构刚度方程,结构刚度方程的建立,图中，数字为总体节点编号1，2，3，4 字母i , j, m为局部节 点编号。,结构刚度方程,结构刚度方程的建立,对应关系，单元 i , j,
36、 m 1，2，3单元 i , j, m 1，3，4 则单元 的节点力列阵为： F =F1x F1y F2x F2y F3x F3y T = F1 F2 F3 T 单元 节点力列阵为： F = F1 F3 F4 T 由节点i处的静力平衡条件可知 R1 = F1 + F1 R2 = F2 R3 = F3 + F3 R4 = F4 （2-49）,结构刚度方程,结构刚度方程的建立,上式中： Ri=Rix Riy T ， Fi=Fix Fiy T （i=1 , 2 , 3 , 4） 又 Fe = Kee 对于单元 有： F1 = K11 1 + K12 2 + K13 3 F2 = K21 1 + K2
37、2 2 + K23 3 F3 = K31 1 + K32 2 + K33 3对于单元 有： F1 = K11 1 + K13 3 + K14 4 F2 = K31 1 + K33 3 + K34 4 F3 = K41 1 + K43 3 + K44 4,结构刚度方程,结构刚度方程的建立,代入2-49式R1 = F1 + F1 =（K11 + K11）1 +K12 2 +（K13 + K13 ）3 +K14 4 R2 = F2 = K21 1 +K22 2 + K23 3 R3 = F3 + F3 =（K31 + K31）1 +K32 2 +（K33 + K33 ）3 +K34 4 R4 = F
38、4 = K411 + K43 3 +K44 4 可把上式写成矩阵形式，并进一步简写为：K =R 称为结构刚 度方程。表示了结构的节点载荷R与节点位移列阵之间的关系。K为结构刚度矩阵或总体刚度矩阵，简称总刚。,结构刚度方程,形成总刚的常用方法,上面是通过节点的平衡关系导出结构刚度方程的，这种做法优点 在于力学概念明确。缺点： 繁琐 不便于程序实现。所以通常 采用下面的两种方法：1.按单元形成总刚 做法：A .先将总刚充0，阶数为44，按节点 节点,结构刚度方程,形成总刚的常用方法,B .从单元 开始，计算单刚 Ke ，送入总刚相应位置，然后进行下一 个单元。,结构刚度方程,形成总刚的常用方法,2
39、 .按节点形成总刚 A .方法同前，总刚充零。 B .从节点1开始，检查该节点与哪几个节点相邻确定总刚的元素 Krs。并与哪几个单元相联系确定元素Krs由几个单元相加 , Krs + C .重复上述工作，直到最后一个节点。 即同样就得到了总刚，即结构刚度方程,结构刚度方程,总刚的性质及其应用,1 .总刚为对称方阵单刚对称阵 叠加后总刚也必然对称应用：在程序设计中只需存储上三角或下三角的元素 2 .总刚是奇异矩阵物理：没有约束，存在刚体位移数学：不存在逆矩阵只有引入位移边界条件后，消去奇异性成为正定矩阵才能求解。,结构刚度方程,总刚的性质及其应用,3 .总刚是稀疏矩阵每个节点只与少数几个节点相关
40、存在大量的元素是有大量0元素的稀疏矩阵,结构刚度方程,总刚的性质及其应用,两种结构的总体节点编号方式，使总刚元素的排列方式不同， 计算表明用带状稀疏矩阵，可节省计算存储量，提高计算效率。在编号时，应使同一单元节点号比较接近，最大节点号差尽 可能小,结构刚度方程,总刚的性质及其应用,4 .总刚仅与结构的尺寸，几何形状及材料性能有关，与外载无关。 利用这一性质，可先计算总刚，再考虑外载的作用,结构刚度方程,位移边界条件的处理,对于结构刚度方程： K =R K为奇异阵 的解不唯一因此，就必须引入位移边界条件，以消除K的奇异性数学：唯一解的必要条件物理：限制刚体位移下面介绍常用的几种方法,位移边界条件
41、的处理,总刚的奇异性,对于结构 而言，其结构刚度方程为：,位移边界条件的处理,总刚的奇异性,位移边界条件的处理,总刚的奇异性,由于结构处于平衡状态，所以：将（2-51）式中代入上式：并将两式相加：,位移边界条件的处理,总刚的奇异性,对任意的vi ，ui上式恒等于0，因此其系数分别等于0即：由对称性可知：上式表示了总刚中各行元素之和均为零，即K对应的行列式的各 行线性相关，根据行列式的性质：“若行列式的各行线性相关，或某 一行是其余各行的线性组合，则行列式=0” 有：总刚是奇异矩阵,位移边界条件的处理,处理位移边界条件的常用方法,由于K和R中的各元素均已按照一定的顺序分别存储在相应的数 组中，在
42、对K及R处理时，应尽量不打乱原有的存储顺序，并希望处 理的元素越少越好。常用的方法有3种。 1 .降阶法：降低结构刚度方程阶次的方法若结构刚度方程为：K =R在节点位移中，令A未知位移， B已知位移，利用矩阵分块：,位移边界条件的处理,处理位移边界条件的常用方法,式中RB为未知载荷，并按第一行展开：KAA A + KAB B =RA令则： （阶次比K =RA低）若B为零位移，则上式变为：上式相当于在原结构刚度方程中，将与零位移约束对应的行与列 划去得到，由该式可以解出未知位移 A 采用计算机解题时，降阶会打乱原K及R的存储顺序，且需重新 安排KAA ，RA及A 在程序设计中一般不采用,位移边界
43、条件的处理,处理位移边界条件的常用方法,2 .对角置一法： 对结构刚度方程,位移边界条件的处理,处理位移边界条件的常用方法,已知位移边界条件： ，将其引入刚度方程。 为了不改变列数，处理第i列诸元素。,位移边界条件的处理,处理位移边界条件的常用方法,为了改变行数，对第i行处理，使其体现,位移边界条件的处理,处理位移边界条件的常用方法,对于平面问题，只要进行三次即可使k成为非奇异阵，故可求出 待解位移。,位移边界条件的处理,处理位移边界条件的常用方法,3 .对角元乘大数法 若： 对第i行的主对角元Kii乘以一个大数，如1020，并将对应的Ri改为：其它各行元素均保持不变。这样将第i行展开得到：同
44、除1020得：,位移边界条件的处理,处理位移边界条件的常用方法,若 则在Kii处乘以大数，Ri处置零即可 这种方法应用的最为普遍。,位移边界条件的处理,基本公式,由结构刚度方程解出后可得到用 e单元内应变、应力与节点 位移间的关系。就可得出单元中任一点处的应变与应力：当不考虑温度影响时：或者 =Be （B几何矩阵）,应力计算,由弹性力学可知：代入上式 根据上式就可得出单元中任一点的应力。,基本公式,应力计算,1 .变温等效节点载荷：设弹性体温度由T1升至T2，则弹性体的变温为T= T2- T1变温T为x,y的函数。 若弹性体内不受任何约束，存在变温T时正应力为T。为线 膨胀系数。,变温应力的计
45、算,应力计算,因此在平面应力状态下：若三角形带应变单元，3个节点处的变温为Ti、 Tj 、 Tm则变温T 可由：T=NiTi+ NjTj + NmTm 为简单起见 T=（Ti+ Tj + Tm）/3 那么当考虑温度变化时的平面问题物理方程为：式中：为单元任一点的总应变。o为该点的初应变（或自由热应 变）。应用单元的虚功方程：得：,变温应力的计算,应力计算,进一步得若令 Rte 为节点的变温等效节点载荷上式成为：即： 2、结构刚度方程：Rt变温等效节点载荷列阵,变温应力的计算,应力计算,由式求出后由：对于平面应力问题：对于平面应变问题：代入后，同平面应力公式。在上面的计算中，应力计算均应按单元进
46、行，同样应变亦应按单 元进行。,变温应力的计算,应力计算,计算得到的是单元的应力与应变，由其表示节点处的应力值的两种常用方法：1 .绕节点平均法：将环绕某一节点的各单元加以平均以平均值表示该节点处的应力。若i节点周围有n个单元，则i点应力为：,应力的表示方法,应力计算,2 .按单元面积的加权平均法：以交集于节点i的各单元的面积作为加权系数来计算i接点处的应 力对于一般情况。这种方法要精确一些。,应力的表示方法,应力计算,根据材料力学公式，可以求出平面问题中节点的主应力、主方向 和当量应力，如下式：当s逆时针转至1时，为正值。,主应力和主方向,应力计算,位移型有限元法求解线弹性问题的普遍公式将三角形常应变单元求解平面问题的公式推广为下列矩阵形式的 普遍公式。 1 .单元中任一点处的位移、应变与应力关系：,有限元法的普遍公式,有限元法的普遍公式,2 .单元刚度方程，单刚及其子刚阵,有限元法的普遍公式,有限元法的普遍公式,3 .单元载荷移置公式集中力体积力表面力变温等效节点载荷 4 .结构刚度方程上述公式适用于各种类型的单元及各种类型的问题，称之为位移 型有限元法求解线弹性静力问题的普遍公式。,
