求曲线箱梁桥自振频率的精细传递矩阵法.pdf

求曲线箱梁桥 自振频率的精细传递矩阵法 孙建鹏 ，李青宁 ，张鹏科。 (1．西安建筑科技大学土木工程学院，陕西 西安 710055；2．西部建筑科技 国家重点实验室 (筹)，陕西 西安 710055； 3．陕西市政建筑设计研究院有限公司，陕西 西安 710068) 摘 要：将高精度的精细积分法和力学概念清晰的传递矩阵法结合起来，以微分方程和矩阵分析理论为基础， 提出了一种新 的精细传递矩阵形式，推导了该方法的计算公式，并用于求解 曲线箱形梁桥的 自振频率．与传统 的传递矩阵法相 比，无需对微分方程组进行求解 ，只需按照迭代公式进行计算 ，就可以得到所需要 的传递矩 阵．根据边界条件 ，运用频率搜索方法便可得到结构的 自振频率．提出的方法虽然是条件稳定的，但是稳定性 条件极易满足．算例表明本文方法 的有效性． 关键词 ：精细传递矩阵法 ；指数矩阵运算 ；曲线箱形梁桥 ；自振频率 中图分类号：TU213．2；TU451．4 文献标识码 ：A 文章编号 ：1006—7930(2009)06-0811一O6 固有振动特性是结构的重要性能，它决定 了结构的地震反应行为和抗震性能．由于曲线梁平面内振 动和平面外振动都会引起弯扭耦合反应 ，其结构动力特性非常复杂性 ，一般而言，平面内振动对曲线梁 桥的动力响应影响相对较小[1]．本文主要研究曲线梁平面外的固有振动问题．曲线梁的固有振动方程为 高阶偏微分方程．用传统的数值微分法求解非常繁琐 ，复杂 ，而且随着单元数 目的增加 ，精度有所下降， 甚至出现发散现象．鉴于以上缺陷，本文将精细积分法 与传递矩阵法 相结合 ，提出了一种新的精细 传递矩阵法 (PTMM)，对曲线梁的 自振频率进行求解． 1 精细传递矩阵法的基本原理与公式 曲线梁的计算坐标系如图 1所示 ，轴表示竖 向，轴表示径 向，z轴表示切向．U是竖向位移 ，a是 竖向弯曲转角 ，是扭转角 ，R是曲线半径 ，是弧长对应 圆心角．梁横截面特性及材料特性沿弧长方向 不变 ． 1．1 曲线梁的应变和位移关系 取曲线梁微段进行分析，考虑竖 向剪切变形后， 得 ： a一+塞一 (1) 式中：是竖向剪切转角；相距 ds的两截面绕 轴 的 竖向挠曲率K 为 ： K 一 一 0 (2) 图1 曲线梁计算坐标系 dZ Fig．1 Coordinatesystem ofthecurvedbeam 相距的两截面绕轴的扭曲K 为 ： K 一 dO 十 a (3) d du 翘曲角为 ： )， ~-Rd (4) z 1．2曲线梁 的力和位 移关 系 Q —kGAv (5) My—E1(塞一景) (6) = 一E1dz7+GI dO_T a , ， (7) 一 d2 dZ