悬索大挠度风振响应的有限体积法分析.pdfVIP

第23卷第4期 计算力学学报 V01．23，No．4 of Mechanics 2006 ChineseJournal August 2006年8月 Computational 悬索大挠度风振响应的有限体积法分析 陈水福。， 李常虹 (浙江大学土木工程系，浙江杭州310027) 摘 要：基于Lagrange原理。建立了一套新的悬索大挠度动力特性和动力响应分析的有限体积法列式，推导了结 点力向量、质量矩阵和单元刚度矩阵的显式表达式。该列式的一个显著特点是直接利用工程应变定义结构变形， 其物理意义明确，列式简单，适用于各种垂庋和荷载情况的悬索大挠度动力分析。实例动力特性和随机风振响应 分析表明，该有限体积列式不仅计算效率高，而且具有良好的计算精度。 关键词：悬索结构；有限体积法；几何非线性；大挠度；风振响应 中图分类号：TU357；TU311．3文献标识码：A 索大挠度(几何非线性)动力特性和风振响应的分 引 言 析。该方法沿索长方向将悬索划分为若干索段，每 一索段称为一个有限体积，然后在变形后的构形上 悬索具有质量轻、柔性大、阻尼低等特点，在正 常工作情况下多处于大位移、小应变的工作状态， 方法定义应变和应力，并进一步求得应变能和动能 对低频荷载如风荷载等十分敏感。悬索在风荷载 的计算式，再根据Lagrange原理导出了悬索大挠 作用下的响应问题是一个较为复杂的几何非线性 度振动的有限体积离散方程。利用该方法进行了 动力问题。 悬索动力特性和风振响应的实例分析。有限体积 目前分析悬索大挠度动力特性和风振响应的 法原是计算流体力学中类似于差分法的一种数值 主要方法有有限单元法(FEM)和其他半解析方 离散方法[5]，本文借鉴了原方法的区域离散格式， 法[1-4]。有限单元法多采用多结点模型(通常为3 但列式推导和离散方程的建立则与有限单元法的 ～5结点)‘2’33或悬链线模型[43模拟索元的变形特 思路一致。本文第一作者及其合作者曾采用改进 性。为考虑大挠度的影响，一般采用Lagrange应 后的有限体积法对纺织物的复杂大挠度变形进行 变和Kirchhoff应力等定义索在某一状态的变形 了静力分析[6’7]，取得了很好的计算效果。 和受力情况。这种应力、应变定义物理意义不够明 确，由此推出的有限元列式较为复杂，计算效率不 2基本假定和索元划分 很高，不利于针对大跨结构进行较长时间段的动力 首先将悬索沿索长方向划分为若干个小索段 响应时程分析。文献[1]从索振动的微分方程出 如图1所示。每一索段为一有限体积(两相邻短线 发，忽略索的纵向振动并假定索的振动可以由其面 围成)；每一有限体积包含一个结点，位于该体积的 内和面外的第一阶振型线性表示，利用Ga