悬索结构大位移分析改进的两节点索单元模型.pdfVIP

第 27 卷第 4 期 河 海 大 学 学 报 Vol 27 No. 4 1999 年 7 月 JOURNAL OF HOHAI UNIVERSITY Jul. 1999 悬索结构大位移分析改进的两节点索单元模型 唐建民　卓家寿 (河海大学土木工程学院 　南京 　210098) 摘 　要 　基于非线性弹性理论 ,提出了适合悬索结构大位移分析的改进的两节点索单元有限元模 型. 在位移函数的表达式中计及了索元自重垂度的影响 ,从而克服了目前悬索结构普遍采用的两节 点索元无法直接计及自重垂度影响的缺点. 采用Lagrange 轴向应变的表达式建立了索元轴向应力、 应变的直接关系 ,并基于Lagrange 描述和虚功原理得到了索元的非线性平衡方程与切线刚度矩阵. 采用 NewtonRaphson 法进行了实例计算 ,结果表明本文方法的精度明显优于两节点直线索单元的 精度 ,适合于中、大跨度悬索结构的大位移分析. 关键词 　悬索结构 ;大位移 ;非线性有限元 中图号 　TU351 悬索结构在工程中有着广泛的应用, 如悬索桥、拉索屋面、拉索塔架以及游览索道等等. 悬索结构在力学 上主要表现为高度柔性和极强的几何非线性. 目前对中、小跨度悬索结构, 当各索段的跨度不大时, 主要是利 用两节点直线单元进行有限元分析, 否则采用多节点单元. 常规的两节点杆单元由于忽略了索元柔性和自重 垂度引起的非线性影响, 当结构中索元的预拉力不大或受有沿索长分布的分布荷载作用时, 会带来不精确的 结果. 采用多节点曲线元固然可以提高精度, 但计算自由度会显著增加. 因此, 本文提出了适合悬索结构大位 移分析的改进的两节点索单元有限元模型. 在位移函数的表达式中计及了索元自重垂度的影响, 采用 Lagrange 轴向应变的表达式建立了索元轴向应力、应变的直接关系, 并基于 Lagrange 描述和虚功原理得到了 索元的非线性平衡方程与切线刚度矩阵. 采用 NewtonRaphson 法进行了实例计算, 结果表明本文方法的精度 明显优于两节点直线索单元的精度, 适合于中、大跨度悬索结构的大位移分析. 1 　位 移 函 数 采用Lagrange 坐标, 建立如图 1 所示的单元局部坐标与整体坐标. ξ( ξ) ( ) 设索元的垂度方程为　　　　　f = 4 1 - f max 1 式中f max为索元跨中最大初始垂度, 且 ξ ( ξ ) ( ) = x ′/ L 　　　 0 ≤ ≤1 2 ( ) Δ ξ( ξ) Δ ( ) 由式 1 , 则 f = 4 1 - f max 3 　　根据索元的几何与物理关系可求得[ 1] 3 f max 图1 　整体坐标与局部坐标下的索单元 Δ Δ ( ) f max = L 4 2 L Fig. 1 　Cable elements under entire and ( ) ( )