波浪作用下斜坡上护面块体断裂破坏的数值模拟.pdfVIP

第 36卷第 1期 水 道 港 口 V01．36 No．1 2015 年 2 月 JournalofWaterwayandHarbor Feb．2015 波浪作用下斜坡上护面块体断裂破坏的数值模拟 崔 焱 ，姜 峰 (大连理工大学建设工程学部，大连116023) 摘 要：基于FEM／DEM方法建立原型尺度数值模型，模拟波浪荷载作用下斜坡上护面块体内部的应力 分布和断裂破坏过程。其中结构物所受的波浪力采用微幅波的理论解计算 ，块体之间的运动、接触以及 块体变形采用FEM／DEM方法模拟 。采用基于势函数的罚函数法计算块体之间的接触力 ，采用Single／ Smeared破坏模型模拟混凝土块体开裂，采用中心差分法的显式积分算法进行数值求解计算。分别采用 文中数值模型与ANSYS模型对 自重作用下块体的应力进行比较分析 ，给出了块体碰撞破坏的两种形 式，验证 了数值模型应力和变形计算的精度。通过数值模拟给出了不同波浪荷载下斜坡上护面块体的 运动、内部应力分布、裂缝形成和断裂过程，给出了不同点的应力历时曲线，探讨了块体的应力变化和断 裂破坏的特性。 关键词：斜坡堤 ；断裂；FEMDEM；护面块体 中图分类号：TV143；0242．1 文献标识码：A 文章编号：1005—8443(2015)01—0026—08 斜坡式防波堤的护面层是一个由离散护面块体相互勾连组成的非连续结构体系 ，块体护面层既具有 很好的消波性能又能在波浪的作用下保持稳定。然而在极端海况下，易发生块体移动或者滚动，导致护面 层水力失稳。同时相邻护面块体之间相互挤压和碰撞，导致个别块体局部可能达到极限状态而断裂破坏口， 从而发生结构失稳。 目前关于斜坡堤护面块体稳定性的研究工作大多基于物理模型试验，如Hudson公式和VanDerMeet~ 公式等经验公式，或是根据实际工程资料在实验室的波浪水池内进行的稳定性实验。 由于物理模型试验中比尺效应的影响，现有模型尺度的物模试验不能正确反映真实结构的强度，同时 原型物模试验又受制于实验条件的限制。因此有必要开展以原型尺度模拟护面块体强度和断裂破坏的数 学模型。近年来 ，随着非连续体计算力学的快速发展 ，已有研究把模拟结构应力和变形的有限单元法和能 够跟踪块体运动的离散单元法结合起来称为FEM／DEM方法 ，该方法保留了有限元和离散元各 自的优势， 可以用于解决多体动力学问题[6和破坏力学问题 。如Latham等人 应用基于FEM／DEM方法的二维数值模 型模拟了矩形断面的规则形状 (方形和十字形)块体在突然失稳状态下的瞬时应力分布。 由于波浪与斜坡堤护面块体的相互作用是一个涉及到多种物理效应的复杂过程 ，关于斜坡堤稳定性的 数值模拟是一项很有挑战性的工作 ，相关的研究工作刚刚起步 。本文基于FEM／DEM方法建立了模拟护 面块体应力分布、裂缝形成及断裂过程的二维数学模型。其中结构物的水动力荷载采用微幅波的理论模 拟，考虑到耦合数学模型的复杂性和计算效率，本文暂且忽略结构物对波浪的作用。结构物的运动和变形 采用FEM／DEM方法模拟。混凝土块体开裂采用Single／Smeared破坏模型。应用所建立的数值模型对线性波 作用下斜坡上规则摆放的护面块体内部应力变化特性和断裂破坏过程进行了模拟分析。 收稿 日期：2014—03—05；修回日期：2014—04—21 基金项 目：国家 自然科学基金项 目 作者简介 ：崔焱(1988一)，男，河北省人 ，硕士研究生，主要从事斜坡堤护面块体断裂破坏的数值模拟。 Biography：CUIYan(1988一)，male，masterstudent． 28 水 道 港 口 第36卷第 1期 阿FFT一] (8)