结构数值分析有限单元法基础.ppt

结构数值分析 有限单元法基础,第一章 有限单元法的总体思路、方法概说,1.1 引言,0、为什么要学有限单元法,1m,2m,P/2,P/2,t,A,求点A的位移?,结构力学方法：,P,弹性力学方法：,当弹性板计算,结构不能简化,（1） 是一种以电子计算机为计算工具的结构分析数值方法。 （2） 是以剖分插值和能量原理为基础的一种结构计算方法。,1、什么是有限单元法, 剖分插值把结构剖分（离散）为有限个单元（小局部），利用“插值函数”研究单元的平衡和协调；再把这有限个离散单元集合（还原）成结构，保证被还原的结构满足平衡和变形协调条件。,图1-1 弹性悬臂板剖分与集合,单元、节点需编号, 能量原理研究单元平衡和结构整体平衡时，都用能量原理描述。, 广泛使用 “矩阵”这一数学工具。,可以说有限元法的 基本思路是：结构剖分和单元集合。 有关平衡的表述采用：能量原理。 主要优点：计算规格化、数值化、 适合计算机编程。,有限元法避免了经典弹性力学连续解的困难，使大型、复杂结构的计算容易地在计算机上完成。,2、 有限元法的分析思路 （1）用点、线或面把结构剖分为一系列离散单元。进行单元分析，使每个单元都满足平衡条件和变形连续条件：,1、F1,2、F2,3、F3,4、F4,1、F1,2、F2,3、F3,4、F4,图1-2 两端固定梁剖分,单元的节点上有位移和力F,（2）再把所有被离散的单元集合起来。进行系统（结构整体）分析，保证系统在单元与单元间连接点处的平衡条件及变形协调条件得到满足。,最终实现对结构的完整分析计算。,图1-3,3、节点和节点自由度数概念 单元与单元间、单元和支座间的连接点称节点。 节点自由度数指：计算时对一个节点应考虑的独立位移数。,节点位移 节点力F,图1-4 平面应力板,一个单元的每个节点上都有用来描述其变形的广义位移和相应的广义力。广义一语意味着所论位移既可以是线位移，也可以是角位移。力除了代表力以外，还可以代表力矩。,4、单元节点位移和单元节点力,上述单元的节点位移和节点力是对给定的坐标系来说的，这里，暂假定各单元具有统一的坐标系XY。,图1-5,把单元上所有节点的位移（或力）依次集合起来排列成一个列向量（或F ），称（或F ）为单元节点位移（或单元节点力），可简称为单元位移（或单元力）。,单元上节点位移总数称单元的自由度数，等于单元节点数乘节点自由度数。 研究单元时F是外力，但就整个结构而论，它是内力。,研究单元时F是外力， 就整个结构而论它是内力。,5、简要归纳,（1）剖分结构时应对单元、节点分别用连续正整数编号。,（2）节点独立位移数是节点自由度数。,平面桁架节点：2（ux、uy）,空间桁架节点：3 （ux、uy 、uz）,平面梁节点：2（ uy 、z ）,空间梁节点：6（ ux、uy 、uz 、 x、 y、z）,平面应力板节点：2（ ux 、 uy ）,节点的独立位移数与结构的变形性质有关。,（3）从结构中取出的任何单元，应有单元号，单元上的节点应有节点号。这些编号应和结构剖分时的编号对应。（十分重要，定义有关数据地址）,杆件单元,板单元,图1-7,（4）每个单元都有单元位移 、单元力F 。它们是把单元上所有节点的位移（或力）依次集合起来排成的一个列向量 （或F ）。,图1-9 平面应变板单元,（5）确定单元位移、单元力应有坐标概念。如果是在特定的单元坐标系中得出的单元位移、单元力 称单元坐标单元位移、单元坐标单元力。 （6）对单元：F是外力； 对结构：F是内力。,1.2 单元分析单元刚度矩阵概念,单元分析的主要工作是：通过研究单元力和单元位移之间关系，建立单元刚度矩阵。 对任意单元而言，描述单元力和单元位移之间关系的一个方阵，称单元刚度矩阵。以图1-10示出的平面梁单元为例。坐标系XY如图所示。,图1-10,（1-1）,描述单元力和单元位移之间关系的矩阵式一般可写为：,其中，kij（i=1、4，j =1、4）称刚度系数， 矩阵:,称为第（e）号单元的单元坐标单元刚度矩阵，可简称为单元刚度矩阵，简写为k。即,（1-2）,引入单元位移、单元力符号： 、,公式（1-1）被缩写为：,（1-3）,从式（1-3）和（1-1）看出：,（1）单元刚度矩阵表明了单元力和单元位移之间关系。由于后面推演这套关系中的刚度系数kij时，保证了单元内部的平衡和协调。所以，引用单元刚度矩阵就意味着单元的平衡和协调条件已经得到满足。,（2）单元刚度矩阵中的任意元素kij是单位位移j=1、其它位移为零时的Fi。,单元分析的主要任务就是确定单元刚度矩阵k。,1.3 系统分析 1、系统分析力学概念 系统分析就是将离散的单元集合起来还原成结构，对结构进行整体分析。 这一过程中，要保证结构在节点处是平衡的；同时相交于同一节点的所有单元在节点处的位移是协调的。,以图1.11所示连续梁为例:,在3号节点处,图1-12,F4 + F2 =0 F3 + F1 + P2 =0 平衡条件,3 = 1 =结构上3点挠度 4 = 2 =结构上3点转角 协调条件,2、系统分析内容 （1）确定受约束自由度的位置、数量； （2）确定结构自由度总数； （3）组集结构刚度矩阵； （4）求解系统节点综合方程获得节点位移解答。,3、确认受约束自由度的位置、数量 约束指结构中那些使节点位移为零的刚性约束，如刚性支座约束等。约束一定发生在节点上（不考虑非节点约束），但不一定约束节点的全部位移。,图1-13 简支梁,简支梁节点自由度为2，1、4点的竖向位移受约束，但角位移不受约束。,悬臂平面应力板剖分为4个单元后，只在节点1、2有约束。,受约束位移的位置：节点1、4的竖向位移受约束的自由度数（ nr ）： nr =2,4、计算结构自由度总数,（1）不计入约束时 结构自由度总数(nf)等于节点总数(nj)乘节点自由度数( ndf ) 。 nf =nj×ndf （1-4）,（2）计入约束后 结构自由度总数（ N ）等于不计约束时的结构自由度总数（ nf ）减受约束的自由度数（nr）。 N = nf - nr （1-5）,例如: 图示平面框架，有4个节点。 每个节点有3个自由度(线位移2、角位移1)。 节点总数 nj=4 节点自由度数 ndf=3 受约束自由度数 nr=6 结构自由度总数 不计约束： nf =4×3=12 计入约束： N =12-6=6,图1-14,5、组集结构刚度矩阵K 建立起系统节点力P与节点位移之间的关系满足结构整体平衡和协调条件的综合方程。,图1-15示出的受集中力作用的两端固定梁 系统的节点位移如图1-16所示。,图1-15,1,2,1,2,3,4,x,y,图1-16,4,3,系统节点力与节点位移关系为：,（1-6）,或缩写为：,P=K （1-7）,其中 K结构刚度矩阵,实际运算中，结构刚度矩阵是由单元刚度矩阵集合而成。集合结构刚度矩阵的过程就是使系统的平衡条件和变形协调条件得到满足的过程。第3章将对此进行深入讨论。,6、关于等价节点荷载,从图1-15、16和式1-7看出，综合方程中的荷载都是节点荷载。然而，工程中的荷载不只是节点荷载，还有线荷载、面荷载。有限元计算中，必须把线荷载、面荷载化为与之等价的节点荷载。,7、求解结构整体平衡方程,从方程（1-7）可解出节点位移。 由此可见，用有限单元法计算结构的力学状态是以节点位移为未知量的。,P=K （1-7）,1.4 计算结构内力和应力 返回每个单元逐一分析 1、根据变形协调条件，从中找出相应的单元位移，并把它换算成单元坐标的单元位移。 2、应用公式 F =k 计算F 。 3、结构内力=F+等价分布荷载过程的附加内力。,总 结,1.1 引言,1.2 单元分析单元刚度矩阵概念,1、什么是有限单元法,2、 有限元法的分析思路,3、节点和节点自由度数概念,4、单元节点位移和单元节点力,1.3 系统分析 1、系统分析力学概念,2、系统分析内容,3、确认受约束自由度的位置、数量,4、计算结构自由度总数,5、组集结构刚度矩阵K,6、关于等价节点荷载,7、求解结构整体平衡方程,1.4 计算结构内力和应力,1.5 有限元思路框图,下接第2章,