钢结构稳定理论-3..ppt

由此可得： 将上式代入（1）中可得单元刚度矩阵为： 其中： 无轴向力时梁单元刚度矩阵 则单元刚度矩阵可以表达为： 利用边界条件，形成结构整体刚度方程后： 可以对上式求解，得到节点变形，及单元内部各点变形（内力等）。 考虑轴力作用的几何刚度矩阵 有关讨论： 几何刚度矩阵[kg]与P、l有关；而[ke]与EI、l有关，与P无关。 方程{q}=[k]{δ}这时还不能求解，应为单元刚度矩阵是奇异的，只有引入边界条件后才能求解。 当结构由若干单元构成时，还应根据单元间变形和内力协调条件，集合成总刚[k总]。 求结构屈曲荷载时，利用能量法，令|[k总]|=0，可得到一个关于P的方程，从而求解，这种方法又称为特征值法；如果外荷载已知，可以求出此时的结构变形。 本节所讲内容只是针对平面梁柱单元，其它单元的求解方法和原理类似，只是刚度矩阵有所不同。 §3-5 有限积分法 将在下一章结合压弯构件弹塑性屈曲进行介绍。 * 则总势能为： 利用势能驻值原理（即总势能一阶变分为0） 利用分部积分和边界条件 可知上式第一项和第二项分别为： 将上二式代入 得 由于δy(l)、δy’(0)、δy’(l)、δy均为边界上不为0的任意值，所以上式等于0的条件为其系数衡等于0： 前三项为边界上的弯矩和横向剪力，即自然边界条件； 而第四项为平衡方程； 可见势能驻值与平衡方程等价； 结论： 很多复杂结构很难建立平衡方程，这时可先写出总势能表达式，令其一阶变分为0，即可得到平衡方程； 可以假设构件的挠曲线函数，（必须满足几何边界条件），将其代入总势能表达式，通过一阶变分为0，求解屈曲荷载，这就是著名的瑞利－里兹法（Rayleigh,L., Ritz, W.）； 如果假设的挠曲线函数既符合构件的几何边界条件，又符合自然边界条件，也可直接利用(4)式求解屈曲荷载，这就是著名的迦辽金法（Galerkin, B.G.）。 §3-2 瑞利－里兹法 假定满足几何边界条件的挠曲线为： ai－待定系数； fi(x)－满足边界条件的函数(至少满足几何边界条件)； 可见挠曲线y为一个泛函（函数的函数）。 将y的表达式代入总势能表达式中——总势能表达为系数ai的函数。 使用势能驻值原理 ，可写成： δa1、δa2、δa3…不能同时为0，则可得到： 给出几种满足几何边界条件的常用挠曲函数 一般取前一、二项就有良好精度 例题1，如图所示变截面悬臂柱，求其轴心受压临界承载力（在工业厂房中常用此结构形式） 根据前面表格，设挠曲线方程为： 满足y(0)＝0，y’(0)=0 应变能： 代入y，并分步积分得： 应变势能或外力功： 代入y，并分步积分得： 则总势能为： 利用势能驻值原理： 讨论如下： 如果利用平衡微分方程求解，精确解为 误差仅为1.25％。 若变形函数只取一项y=a f (x)时，即只有一个参变量a时，可以证明： 能量法结果偏高。 因为挠曲函数与真是变形不完全相符，相当于对真实杆件人为添加了若干横向水平约束，提高了压杆的抗失稳能力。 例题2，连续变截面杆件，假设忽略腹板的存在。求解其临界力。 下部截面I1、h1已知，其它截面特征为： 解： （1）用静力平衡方程求解： 变系数的微分方程，很难求解。 （2）用能量法求解： 首先根据边界条件选用合适的变形曲线 若选用的变形曲线只有一个参变量，则： 为了方便求解，令m＝0.5，则h2/h1=0.5，即I2=0.25I1。 设挠曲线为： 可见虽然变形曲线中的参数a是未知的，最后也没有求出a的具体数值，但确通过此方法求出了临界力。 §3-3 迦辽金法(Galerkin, B. G.) 迦辽金法直接利用势能驻值原理中的平衡微分方程求解，不再需要写出总势能表达式。 但这样做的前提是挠曲线必须既满足几何边界条件，又同时满足自然边界条件。 由前面可知平衡微分方程为： 令： 则： (1) (x1,x2更适于普遍情况) 设位移函数 其一阶变分为： (2) 将(2)代入(1)式 由于δa1、δa2、…δan都为不等于0的微小的任意值，而(1)式是恒等式