斜坡式防波堤前海冰堆积数值模拟.pdfVIP

水 利 学 报 2003 06 SHUILI XUEBAO 6 年 月 第 期 文章编号：0559-9350（2003）06-0105-06 斜坡式防波堤前海冰堆积数值模拟 1 1,2 1 1 王永学 ，李春花 ，孙鹤泉 ，李志军 (1.大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116024；2.国家海洋环境预报中心，北京 100081) 摘要：本文将海冰在斜坡式防波堤前的堆积过程简化为二维的动力学问题，假定海冰为粘弹性体，在多边块体运动 模拟系统基础上，采用离散元数值方法对冰排爬升堆积过程进行了实时仿真模拟。通过与二维理论断裂模型及物理 模型试验结果对比，证明了此数值模型的合理性与可行性。 关键词：海冰堆积；防波堤；离散元法 中图分类号：P731.15 文献标识码：A 海洋动力因素造成的海冰变形断裂、重叠与堆积过程的机理以及海冰与结构物的相互作用问题长期以 来一直是海岸工程领域的困难课题。国内外当前海冰堆积的理论成果基本上是假定冰排为一均质的半无限 长弹性板，在风和流的驱动力下与斜面结构相互作用，冰排在爬升堆积过程中只发生弯曲断裂。由于斜面 倾角较大时，屈曲或挤压破坏的成分不可避免，所以早期国际上发展的斜面冰堆积和冰作用力模式都是针 对小倾角斜面。对于海冰堆积过程的数值模拟，曾有学者试图利用有限元和边界元法来处理，但发现此类 方法只能模拟冰排的初始断裂[1,2]。目前，离散元法被认为是模拟海冰堆积过程比较有效的一种数值方法， 它能弥补有限元法和边界元法的不足，此方面的代表性工作有：Hocking等人假定海冰为弹塑性体，模拟 了冰排在坡度较小的人造岸坡上的爬升与堆积过程[3,4]；Hopkins还对小角度斜坡上海冰堆积过程进行了数 [5] 值模拟，并与实验结果进行了对比，但模拟现象与实验现象有所出入 。 本文模型在多边块体运动模拟系统基础上[6]，假定海冰为粘弹性体，设立了海底及斜坡式防波堤坡面 边界条件，利用离散元法来实时模拟斜坡式防波堤前海冰的堆积过程。此模型可以模拟海冰在不同倾角坡 面上的堆积过程，其中包括大倾角坡面，另外，所设立的海底边界条件还能够模拟冰堆触底后的堆积过程。 为了检验此数值模型，本实验室利用DUT1非冻结可破碎合成模型冰材料对斜坡式防波堤前海冰的堆积过程 进行了对比试验[7]。根据相同参数下数值模拟结果与物模试验结果对比，二者模拟的爬升堆积过程现象非 常相似，对建筑物的力荷载过程线趋势基本一致，只是力的幅值与频率有所差异；此外，数值模拟中冰排 初始断裂对应的水平荷载和垂向荷载比值与理论计算值几乎完全一致。 1 控制方程 图1所示为冰排与斜坡式防波堤接触前的状态，此时冰排以恒定的速度u 沿x方向前进，当冰排与防波 0 堤表面接触后，冰排将沿着坡面上爬，并产生变形断裂破坏，断裂的块体或是沿坡面继续上爬或滑入水中 在堤前形成堆积。 如图1所示，斜坡堤及海床为两个固定边界，假定冰排由若干冰块粘结组成，根据动力学原理，第i个 冰块在运动及与结构物相互作用或块体间相互作用时，满足以下方程， 收稿日期：2002-02-08 基金项目：国家杰出青年基金；国家自然科学基金。 作者简介：王永学(1955-)，男，辽宁大连人，教授，博士生导师，从事波浪、海冰方面研究。 105 水 利 学 报