钢结构稳定理论与设计-31..ppt

第三章 稳定问题的近似解法 这是一组齐次方程式，有解的条件是其系数行列式为0 解得(ql)cr=7.838EI/l2 ，与精确解几乎一致。 可见挠曲线由一项改为两项时效果就十分显著。 §3-4 有限元法 Finite Element Method (FEM) 通用程序有很多：Algor, Sap, Ansys, Adina, Abaqus, Nastran, CFD等等 是目前较为常用的数值解法，适用于用计算机对复杂结构系统进行求解 可以分为平面模型、空间模型等 可以使用不同单元求解不同问题，如梁柱单元、壳体单元、实体单元、索单元等等 分析时可以考虑线性、几何非线性、材料非线性等 有限元法的基本思路和过程如下： 选择适合的单元形式 把结构或构件划分成若干单元（网格） 建立单元节点位移和内力的关系,通常以矩阵形式表示(使用能量法或平衡方程) 将单元组合成原系统(单刚→总刚) 求解线性或非线性方程组 得到结构荷载位移曲线、屈曲荷载等 以单元节点位移为未知量 利用边界条件 使用直接求解法或迭代法 例：受轴向力作用的平面梁柱单元 (单元在局部坐标系中变形前后如图) 杆端内力以矩阵形式(列向量)表示为： 1、3剪力，2、4弯矩 与杆端内力对应的位移为： 杆端内力与位移之间存在如下关系： 如何建立单元刚度矩阵[K]是问题的关键 下面根据能量法建立单元刚度矩阵 单元内部应变能： 内外力对单元做功： 根据U＝W得： （1） 用一个三次抛物线代替单元挠曲线（形状函数） 内力q引起 轴力P引起 单元两端几何边界条件为： 利用上述条件求出变形曲线中的a，b，c，d后，代回到变形曲线中去，得到变形曲线为： 哈 尔 滨 工 业 大 学 钢 结 构 稳 定 理 论 为什么要近似求解？ 非等截面构件 压力沿轴线变化的构件 具有变系数的平衡微分方程 压杆的弹塑性屈曲问题 近似求解方法： 能量法、瑞利－里兹法、迦辽金法、有限差分法、有限积分法、有限元法等。 §3-1 能量法 1）几个基本概念 保守系统：体系由平衡位置1变化到平衡位置2时，力系（包括内力和外力）做的功仅与始末位置有关，而与中间过程无关的系统。 能量守恒：如果贮存在结构体系中的应变能等于外力所做的功，则该保守系统处于平衡状态，此谓之能量守恒。 能量准则：当一保守系统处于平衡状态时，其总势能的一阶变分为0。 总能量： 外力势能 内力势能、应变能 一阶变分： 势能驻值原理 外力势能增量 内力势能增量 二阶变分： 为稳定的平衡状态，此时总势能最小 为不稳定的平衡状态 为随遇平衡状态 2）弹性直杆的总势能表达式 以直杆状态为参考状态（总势能为0状态），求微弯后总势能表达式 应变能 外力势能 总势能 E、P和I可能不是常数，P(x)、I(x)不可拿到积分号之外。 3）势能驻值原理 当作用着外力的结构体系，其位移有微小变化而总的能量不变，即总势能有驻值时，则该结构体系处于平衡状态，此即为势能驻值原理。 表达为： 其中： 外力势能 内力势能 外力作功 例：现考查如下结构 上端B点具有弹簧铰rB、抗侧移弹簧kB； 下端A点具有弹簧铰rA。 当由(a)直杆状态过渡到(b)微弯状态时，体系的应变能为： 外力势能为： 则总势能为： 利用势能驻值原理（即总势能一阶变分为0） 利用分部积分和边界条件 可知上式第一项和第二项分别为： 将上二式代入 得 由于δy(l)、δy’(0)、δy’(l)、δy均为边界上不为0的任意值，所以上式等于0的条件为其系数衡等于0： 前三项为边界上的弯矩和横向剪力，即自然边界条件； 而第四项为平衡方程； 可见势能驻值与平衡方程等价； 结论： 很多复杂结构很难建立平衡方程，这时可先写出总势能表达式，令其一阶变分为0，即可得到平衡方程； 可以假设构件的挠曲线函数，（必须满足几何边界条件），将其代入总势能表达式，通过一阶变分为0，求解屈曲荷载，这就是著名的瑞利－里兹法（Rayleigh,L., Ritz, W.）； 如果假设的挠曲线函数既符合构件的几何边界条件，又符合自然边界条件，也可直接利用(4)式求解屈曲荷载，这就是著名的迦辽金法（Galerkin, B.G.）。 §3-2 瑞利－里兹法 假定满足几何边界条件的挠曲线为： ai－