桁架结构几何非线性问题的力密度法-应用力学学报.PDF

第28 卷 第5 期 应 用 力 学 学 报 Vol.28 No.5 2011 年 10 月 CHINESE JOURNAL OF APPLIED MECHANICS Oct. 2011 文章编号：1000- 4939(2011) 05-0509-05 桁架结构几何非线性问题的力密度法 王勇 高冀峰 （华南理工大学 510640 广州） 摘要：桁架结构大都具有较强的几何非线性，受荷载后易出现较大的几何变形。提出一种基于力 密度的针对桁架结构几何大变形问题的解法；引入杆件单元的力密度矩阵，推导出相应非线性方 程的Jacobi 矩阵；与有限单元法集成求解的思想相同，采用力密度矩阵建立结构变形后整体的精 确非线性平衡方程。研究结果表明：应用Newton-Raphson 迭代法求解，采用适当的迭代收敛精度 可得到较精确的桁架结构位移解。 关键词：桁架结构；几何非线性；力密度矩阵； Newton-Raphson 迭代法 中图分类号： TU31 文献标识码：A 相对于上述求解方法，本文提出基于力密度的 1 引 言 桁架结构几何非线性问题的求解方法。力密度方法 是由 Shek[6]最早提出的，后由 Grundig[7]等人进行了 空间桁架结构是工程中常用的一种简化结构， 完善和发展，其主要是用于解决膜结构的找形问 如空间网架、空间网壳、塔架结构等。这些结构由 题，在求解过程中力密度值为定值。 于构造或形状的特殊性以及工作荷载的特点，工作 本文基于力密度思想，应用变化的力密度值， 时常常产生相对较大的几何变形，但材料仍然处于 由杆单元的力密度矩阵，建立了杆单元变形后的几 弹性阶段，即结构处于“大位移、小应变”的几何 何非线性平衡方程；采用 Newton-Raphson 法进行 非线性状态。在空间杆系结构分析的有限元法中， 迭代求解，并通过算例验证此方法对桁架结构位移 研究人员已给出了一些具体的实现方法。 解计算的精确度。 文献[1]通过结构变形后的状态建立平衡方程， 推导出空间桁架大位移的全量控制方程，但其刚度 2 力密度方法 矩阵不对称。文献[2]从连续体力学的角度出发，基 于修正的总体 Lagrange 坐标描述法推导了空间桁 架弹性大位移问题的增量平衡控制方程，得到了一 对称的双刚度矩阵，但在每次迭代下需要对其中的 Euler 应力项进行修正。文献[3-4]在推导切线刚度矩 阵时均进行了一些简化。文献[5]基于能量原理推导 出建立在变形后的空间桁架单元平衡方程，并给出 了桁架单元精确的切线刚度矩阵，计算效果较好。 图 1 杆单元模型 来稿日期：2010-11-30 修回日期：2011-06-07 第一作者简介：王勇，男，1956 年生，硕士，华南理工大学土木与交通学院，副教授；研究方向——索膜结构、张弦结构等大跨度空间结构。 E-mail: cvywang@scut.edu.cn 510 应 用 力 学 学 报 第 28 卷 如图 1 所示，设杆件单元节点力分量为 e e e e