无限自由度压杆稳定数值解析法程序设计研究.pdfVIP

第二篇 工程建设软件开发技术与应用 无限自由度压杆稳定数值解析法程序设计 于宝林霍学军 李强 (天津大学建工学院) 【摘要】弹性理想压杆为无限自由度体系，本文采用数值解析法直接求解分支点失稳临界状态的 偏微分方程，用FORTRAN90编制了无限自由度理想压杆临界荷载和失稳模态的求解程序，精确 的求解了Euler压杆的临界荷载和相应的失稳模态。通过与精确解、现有软件求解结果的对比分析 发现，本程序所的结果准确、可靠、精度高。 【关键词】理想压杆；临界荷载；失稳模态；常微分方程：COLSYS 1引言… 随着高强轻质材料的发展与建筑高度的不断增加，。结构的长细比愈来愈大，截面的尺 寸日趋减小。在高层、特别是超高层建筑结构中，由于结构的侧向刚度相对减弱，加之高 层建筑巨大的竖向荷载，结构的整体稳定和局部稳定已引起了学术界的广泛关注。压杆作 为结构中常见的承载构件之一，其稳定分析是有重要的工程和理论意义的。 2 常微分方程特征值问题及其变换b1 竖向压杆根据两端约束情况的不同，所建立的微分方程形式也不同。图示Euler压杆 在临界状态下，体系出现了新的平衡形式，如图中曲线所示，由结构力学可知， 2一 (1)(1) 边界条件为：x=0，y(0)=0；X一--1，Y(1)=0 以上所建立的实际上是一个Euler杆临界荷载和失稳模态 广滁 圈 Differential Y 对应上述ODE问题(Ordinary Equation)的特 征值问题。可以用数值解析法，本文中采用常微分方程编制问 题通用程序COLSYS进行求解。在调用COLSYS时，由于其只能 求解标准的线性和非线性ODE问题，并不具备求解特征值问题 一 的功能，难以直接应用。所以我们还需对问题本身进行适当的 变换，将其转换为标准的ODE形式，使其能够为COLSYS所接受． 本问题的解答被重新描述为： ·对于Euler压杆，杆件单元的微分方程式(1)为： 心r}I棼■蝌字酊 第十三届全国工程建设计算机应用学术会议论文集 广东佛山·2006 EO萨2y+毋=。一x∈【o，1】 (2) ·” 式中， A棚。这里及以下，y-y@)，(·)’=d(·)／dx。 ·基于COLSYS的要求，对方程中的未知常数引入平凡ODE． ‘一 名=0 (3) ·为了唯一地确定相应于入的特征函数向量，取以下归一化条件： y2M=贴)：，2啦 。(斗) 定义 +尺(工)=Y。 (5) 则积分(4)式等价为如下的ODE问题： t． R=、，2(abx)j ’(6 其中：’ R(以)=OR(b)=-I+ ·边界条件：对于支座结点(1，，。r)一结点边界约束信息实型数组v=l，表示单元 横向有约束；v--O，表示单元横向自由；