河冰单轴压缩破坏过程细观数值仿真-水利学报.PDF

水 利 学 报 2018年11月 SHUILI XUEBAO 第49卷 第11期 文章编号：0559-9350（2018）11-1339-07 河冰单轴压缩破坏过程细观数值仿真 1，2 2 3 邓 宇 ，王 娟 ，李志军 （1．黄河水利委员会 黄河水利科学研究院，河南 郑州 450003；2．郑州大学 水利与环境学院，河南 郑州 450001； 3．大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116024） 摘要：为了解并认知河冰材料的内在属性，对黄河封冻期的河冰进行了现场取样观测并开展了单轴抗压破坏过程 的物理试验，依据物理试验结果，在细观层面将河冰看成由晶粒、晶界、初始缺陷构成的非均质材料，给出了 随机缺陷界面弹簧元模型模拟河冰单轴压缩强度的数值方法，模拟了河冰在单轴受压荷载作用下的开裂破坏过 程，分析了初始缺陷分布和含量、晶粒尺寸等各相组分对河冰力学性能的影响。结果表明：随着初始缺陷含量的 增加，河冰单轴压缩强度呈下降趋势，但下降幅度逐渐变缓；随着晶粒尺寸的增加，河冰单轴压缩强度呈降低趋 势，且与晶粒尺寸 d -1/2 呈线性相关关系。与物理试验对比发现，该计算模型能较好地模拟黄河冰的单轴压缩强 度和开裂过程。 关键词：河冰；单轴压缩强度；数值仿真；黄河 中图分类号：TV131.4，P332.8 文献标识码：A doi：10.13243/ki.slxb 1 研究背景 河冰内部细观微裂纹的变化是导致其复杂宏观表象的重要因素，对河冰内部细观微裂纹开裂与 ［1］ 扩展机理的研究将有助于从本质上解释河冰开河期的宏观力学行为 。随着细观力学理论的不断完善 和数值模拟方法的成熟，利用数值仿真技术，研究冰的细观结构与其宏观力学特性的关系已成为冰 ［2］ ［3］ ［4］ 研究的一个热点。相关学者围绕冰的晶粒位错堆积 、晶界滑移 和晶体弹性各向异性的机制 等细 ［5］ ［6］ ［7］ ［8］ 观结构理论已开展了大量的研究工作。近年来，关于高原冰 、冰川冰 、极地冰 、海冰 、大气 ［9］ 冰 等有关冰材料细观结构观测研究成果也不断涌现，与上述冰类型相比较，河冰由于自身材料的特 殊性、形成及演化过程的多变性及有效数据获取的季节性和时效性，其细观结构、物理力学性能及 开裂演变规律的研究仍处于起步阶段［10-11］。相关学者在冰的数值计算方面也取得了长足的进步，王军 ［12］ ［13］ ［14］ 等 对冰塞形成及其演变机理问题、杨开林等 对渠道冰问题、季顺迎等 对海冰与海洋平台结构 ［15］ ［16］ 物相互作用问题、李志军等 对海冰断裂问题、韩雷等 对冰板与锥面结构相互作用中的弯曲破坏 ［17］ 问题、茅泽育等 对冰盖下水流流动的掺混特性问题均开展过数值模拟研究，但基于细观冰力学理 论，采用数值模拟方法，从细观尺度出发对河冰结构特性、材料特性及微裂纹的演变规律进行研 究，进而构建河冰宏细观多尺度的响应机制，现有研究成果仍有待于进一步丰富。 本文采用数值模拟与试验实测资料分析相结合的方法，以黄河河道内封冻期的天然河冰为研究 对象，构建了黄河冰强度的细观计算模型，从细观尺度对河冰的单轴压缩强度进行了分析，对不同 初始缺陷含量变化、不同晶粒尺寸变化下河冰压缩破坏过程进行了模拟。该方法不仅能够直观地描 收稿日期：2018-