桥梁遭受船舶撞击的有限元仿真分析.pdf

桥梁遭受船舶撞击的有限元仿真分析 阙水杰 张 南 樊文才 (南京工业大学土木工程学院 南京 210009) 摘 要 桥梁在船舶碰撞时受到的动力荷载是复杂 的动力非线性 问题 ，近代非线性有限元技术为该 问题提供 了有效 的工具 。文中简述 了该技术 的基本原理 ，运用 ANsYs／IS—DYNA对一艘 10000t 船和桥梁 的碰撞进行 了模拟 ，并建立 了3个不同的桥墩有限元模型，从碰撞力、桥墩应力、变形等不同 角度分析了不同的建模方式对碰撞结果的影响。结果表明，3种模型的计算结果具有较好 的一致性， 可 以满足一般的工程精度要求 。 关键词 船桥碰撞 ；结构非线性 ；建模方式 ；丁程精度 中图分类号 ：U66．14 文 献标志码 ：A DOI：10．3963／j．ISSN 1674—4861．2010．01．030 近几年国内发生的船桥碰撞事故给人类的生 分法求解碰撞 问题要注意时问步长的选取 ，因为 命 、财产 以及社会和环境造成 了极大的捐审 ■ 中心差分法是条件稳定的，其 时间步长不能超过 影响到水运 、公路和铁路运输的安全性。因此，对 临界时间步长。实用中常以最小有限单元网格的 于船桥碰撞 问题的关 注 由来 已久 ，目前船桥碰撞 特征长度除以应力波速来近似临界时间步长，即 的计算主要有 以下几种方法 ：Minorsky理论、汉 r △f≤At。minf ＼ 斯一德鲁彻理论、近似数值解法 、简化解析法 、试 L／， 验方法、有限元法。 1．2 碰撞 中的接触算法 目前 ，已经有一些优秀的软件可 以用于结构 Is—DYNA3D程序 中缺省的接触算法为对 碰撞分析 ，于是如何建立 1个准确而又经济高效 称罚函数法。其基本 的分析原理为 ：每一时步检 的船桥碰撞计算模型成为数值仿真的关键 。本文 查从节点是否穿透主面，没有穿透则对从节点不 利用非线性有限元技术 ，实现船舶与桥梁碰撞过 做任何处理 ；如果穿透则在该从节点与被穿透的 程的数值仿真 ，并对不同的桥墩建模方式进行 比 主面之 问引入 1个较大 的界面接触力 ，其大小与 较，探讨碰撞 中的碰撞力 ，桥墩及与桩基础 内部应 穿透深度 、主面刚度成正比，这个接触力亦称为罚 力分布与变化情况。 函数值。接触力 由下面公式计算 ： fs一 志△ 1 碰撞 问题的非线性计算理论 式中：尼 为接触 刚度 因子 (由单元 尺寸和材料特 1．1 非线性有 限元控制方程 性等确定)；A 为穿透量。 碰撞问题的运动方程可 以一般地表示为： 定义了接触关系后 ，就能得到最大摩擦力值 Mn+ Cll，+Kd= F 为 F 一 l-厂I，式 中： 为摩擦因数，厂 为从节点 式中：M 为船桥质量矩阵；：C为阻尼矩阵；K 为 的法 向接触力 。根据作用力与反作用力原理 ，计 刚度矩阵；a为加速度 向量；为速度 向量 ；d为位 算分配到对应主段上各个主节点的摩擦力。 移 向量 ；F 为包括碰撞力在 内的外力 向量 。碰撞 2 有