基于LS—DYNA的流冰撞击坝体仿真计算.pdfVIP

基于 LS—DYNA的流冰撞击坝体仿真计算 张凤德 李广一 (黑龙江省水利工程建设局 哈尔滨 150081) (中水东北勘测设计研究有限责任公司 长春 130021) 【摘 要】 寒冷地区河流的流冰问题一直为工程界和学术界所关注，尤其是 “武开江”河流春季流冰的冲击力巨 大，常造成水工建筑物的冰损事故。本文运用非线性动力有限元分析软件Ls—DYNA模拟流冰与丁坝的撞击过 程，得到 了与实际情况比较符合的挤压破坏冰力时程曲线。用Ls—DYNA进行该类仿真分析直观、简便，可用于模 拟复杂的流冰撞击问题。 【关键词】 流冰 有限元 Ls～DYNA 仿真分析 1 前 言 数种程序合成一个软件包，称为Ls—DYNA，大大加强 我国北方寒冷地区河流春季凌汛期普遍存在流冰 了处理碰撞问题的功能。LS—DYNA程序是功能齐全 现象，流冰对水工建筑物的撞击破坏作用很大，撞击作 的几何非线性 (大位移 、大转动和大应变)、材料非线 用受冰强度 、冰破坏形式、风速 、水流速、流冰方向、水 性 (140多种材料动态模型)和接触非线性 (50多种) 工建筑物的形式及破冰体的形式等诸多因素的影响， 程序。它以Lagrange算法为主，兼有 ALE和 Euler算 其计算十分复杂。本文运用 Ls—DYNA分析软件 ，对 法；以显式求解为主，兼有隐式求解功能；以结构分析 流冰撞击丁坝过程进行仿真分析 ，探讨流冰对固床撞 为主，兼有热分析、流体一结构耦合功能；以非线性动 击力有限元仿真模型的建模方法 ，为工程技术人员解 力分析为主，兼有静力分析功能。它是一款军用和民 决类似问题提供参考。 用相结合的通用结构分析非线性有限元程序。 2 流冰撞击坝体建模 2．2 模型基本参数 2．1 Ls—DYNA简介 由于冰排漂浮于水面之上，因此可以只考虑风力 、 Ls—DYNA是世界上著名的通用显式动力分析程 水流拖曳力等水平方向荷载，忽略重力、浮力等垂向荷 序 ，能够模拟真实世界的各种复杂问题，特别适合求解 载。水平方向荷载由冰排的初始速度体现，不同大小、 各种二维 、三维非线性结构的高速碰撞、爆炸和金属成 厚度的冰排以不同的速度撞击丁坝时，撞击力是不同 型等非线性动力冲击问题，同时可以求解传热、流体及 的，建模时以冰排为模拟对象计算流冰对丁坝的撞击 流固耦合问题。在工程应用领域被广泛认可为最佳的 作用，模型基本参数如下： 分析软件包。与实验的无数次对比证实了其计算的可 a．模型尺寸：冰排尺寸为80m×40m，厚度0．7m。 靠性。1976年，美 国的 LLNL实验室 (LawrenceLiver— b．模型边界条件：坝体底面及丁坝一端为全约 moreNationalLaboratory)在 J．O．Hallquist主持下研制 束，其他均为 自由表面。 开发了第一代显式非线性动力分析通用有限元程序 C．模型荷载：冰排水平方向初速度 0．8m／s，作用 DYNA3D和DYNA2D。这是碰撞动力有限元软件发展 在冰排上的风力、水流推动力等均由初始速度体现。 的起点，并且由于这一程序的源代码对外开放，现有的 d．模型单元：丁坝与冰排均采用SOLID164—8节 全部显示动力分析有限元软件均起源于这两个程序。 点六面体单元。 1988年，J．O．Hallquist创建了LSTC公司，对 DYNA系 e．材料模型参数：坝体的参数选用散抛块石的弹 列软件进行商业化工作，1997年，LSTC公司将旗下的 模、密度和泊松比；冰在不同温度下其力学参数差别较