第十一章 桥梁结构几何非线性.ppt

第十一章 桥梁结构几何非线性 计算理论 11.1 概 述 11.2大跨度桥梁几何非线性分析的有限元方法 11.3桥梁结构分析常用单元的切线刚度矩阵 11.4桥梁结构几何非线性分析特殊问题的讨论 11.5非线性方程组的求解 11.6小结 11.1 概 述 Oden说过“我们生活在一个非线性世界里”。早在十九世纪未，科学家就发现，固体力学的经典线性理论在许多情况下并不适用，于是开始了对非线性力学问题的研究。二十世纪中，科学家奠定了非线性力学的理论基础。但由于计算繁复，许多非线性微分方程的边值问题无法求解，用解析法解决非线性工程问题仍显得无能为力。直到二十世纪六十年代末，有限元法与计算机相结合，才使工程中的非线性问题逐步得以解决。 11.2 大跨度桥梁几何非线性 分析的有限元方法 或： (11－60) 其中, 是轴力N的函数，称为稳定函数，其值随?的变化而变化。以稳定函数表达的刚度系数包含了轴力对弯曲刚度的影响，相当于前面切线刚度阵中弹性刚度系数与几何刚度系数之和。图11.11给出了稳定函数与单元切线刚度系数随?变化的情况，说明 时，二者基本一样，当 时，二者区别也不大。 时，二者基本一样，当 时， 二者区别也不大。 图11.11 稳定函数与单元切线刚度系数随?的变化曲线 事实上，将C展开成级数形式有： (11－61) ) 而 (11－62) 对比表明，几何刚度阵系数就是稳定函数忽略高阶项的轴力影响系数。 从以上讨论可以看到：当 时，随着?的增大， 几何 刚度矩阵的误差也增大。但由于?与 成正比，有限元分析中，只要减小单元长度，就可避免使用几何刚度阵产生的这种误差。 11.4.2 弯矩对轴向刚度的影响 杆单元的弯曲将引起杆件计算长度(杆件两端节点的距离)的改变，从而影响杆件的轴向刚度。在杆件微段上，弯曲引起的杆件轴线计算长度的改变量 为： (11－63) (11－64) (11－65) 在外力作用下，杆件总的缩短量?为： (11－66) 因此 (11－67) 式中：? ?? 弯矩引起的轴向刚度修正系数。 由于: (11－68) (11－69) 式中：?由式(11－52)定义，式(11－69)代入式(11－65)积分，由式(11－66)、(11－67)可得 (11－70) 其中： 当 P0 (压杆)时 (11－71) 当P0 (拉杆) 时 (11－72) 11.4.3活载的几何非线性分析 为了讨论方便,将非线性状态下荷载最不利加载区域称为影响区。活载几何非线性分析，会遇到如下问题： ? 线性叠加原理失效，无法再用传统的影响线加载法进行活载分析。 ? 确定影响区本身是一个非线性问题。仅用恒载初始状态计算活载，会带来影响区范围改变和不正确载位引起的误差。 　 ? 单位强迫变位产生的等效力很大，用机动法求解影响区将破坏指定状态结构影响区的真实形状。 考虑到结构在确定初始状态下，影响线函数的大小，仍代表单位荷载作用对关心截面计算参数的影响。因此，可以按如下方法计算非线性活载最不利响应：以结构恒载状态为初态，计算影响区函数，用相应的最不利活载作为试探，求出第一次近似，将前一次试探活载与恒载共同作用时的状态代替初态，重新计算影响区和最不利荷载。当本次活载效应与上次活载效应的误差落在某一允许范围时，计算收敛。这样,每叠代一次，都加强了活载对计算状态的影响。这一计算可归结为如下步骤： 1)将结构恒载受力状态作为求解影响区的初始状态，计算出初始影响函数。 2)用动态规划加载法，找出最不利加载位置，并作好记录。 3)以恒载受力状态为计算初态，将活载按最不利载位一次性作用于结构，分析恒、活载共同作用下的结构受力状态和关心截面力学量。 4)将恒、活载共同作用下的结构状态作为求解下一步影响区函数新的初态，重复1)～3)的计算。经过数次迭代计算就可以得到活载作用下关心力学量的最值。 求解活载影响区可用机动法,但单位强迫变位必须取用一个很小的数(如 ），这样，可以保证确定的影响区不失真，使动态规划法找到的载位即为相应内力状态下的最不利载位 。 11.4.4 几何非线性调值计算 桥