渤海海冰动力学中的粘弹塑性本构模型-水科学进展.PDF

第 13 卷 第 5 期 水 科 学 进 展 Vol13 ,No5 2002 年 9 月 　　 ADVANCES IN WATER SCIENCE 　　Sep . ,2002 　 渤海海冰动力学中的粘弹塑性本构模型 季顺迎 , 岳前进 , 姚 　征 (大连理工大学工业装备结构分析国家重点实验室 , 辽宁 大连 　116023) 摘要 : 在粘塑性海冰本构模型的基础上 , 将 KelvinVoigt 粘弹性本构理论引入到海冰动力学中 , 进而建立了粘弹塑性海冰本构模型 。该模型可较好地反映渤海海冰在小应变和小应变率条件下 的粘弹性力学行为 , 同时还考虑了大应变率下的海冰粘塑性力学行为 。对渤海辽东湾海冰进行 了48h 数值模拟 。结果表明: 粘弹塑性海冰本构方程较 Hibler 的粘塑性本构模型可更好地处理渤 海的冰间相互作用 , 提高海冰数值模拟的计算精度 。 关 　键 　词 : 海冰动力学 ; 本构方程 ; 粘塑性模型 ; 粘弹塑性模型 ; 海冰数值模拟 ( ) 中图分类号 : P 73115 　　　文献标识码 : A 　　　文章编号 2002 0559906 海冰数值模拟中的海冰动力学是影响海冰漂移和形变的重要因素 , 其中表征冰间相互作用 的海冰本构模型最为复杂 。在早期的海冰动力学研究中 , 人们将海冰视为二维连续体 , 分别建 立了线粘性 、弹塑性和粘塑性海冰本构模型[ 1～4 ] 。在极区及其冰缘区的中、长期海冰数值模拟 中 , Hibler 的粘塑性海冰本构模型得到了最广泛的应用 。为提高海冰长期数值模拟的计算效 率 , 该粘弹塑性海冰本构模型忽略了海冰的弹性变形 , 仅将海冰变形视为海冰应变率的函 数[4～6 ] 。基于粘塑性本构模型 , Flato 和 Hibler ( 1992) 为研究极区海冰对全球气候的影响 , 提 出了 “空化流体”模型以增强海冰数值模拟的计算效率[7 ] 。Hunke 等 ( 1997) 建立了弹粘塑性 海冰本构模型 , 以提高短期海冰数值模拟中海冰内力的计算精度[8 ] 。Lu 等 ( 1998) 考虑冰体 粘弹性变形行为 , 对 Hibler 的粘塑性本构模型进行了改进以研究冰塞的形成机理[9 ] 。 在渤海海冰的热力 - 动力耦合的中、小尺度数值模拟中 , 为进一步提高海冰动力学的计算 精度 , 应充分考虑小应变和小应变率下海冰形变中的粘弹性力学行为 。另外 , 尽管 Hibler 的粘 塑性海冰本构模型给出了海冰屈服函数 , 但在以往海冰数值模拟中又将一特定的应变速率值或 块体粘性系数作为海冰屈服判据 , 这使得海冰在屈服前后的应力状态不连续 。 本文将充分考虑海冰在小应变和小应变率下的粘弹性力学行为 , 并引入 Kelvin - Voigt 粘弹 模型 , 在 Hibler 粘塑性海冰本构模型的基础上建立粘弹塑性海冰本构模型 , 进而对渤海海冰动 力学进行较为精确的数值计算 。 1 　粘塑性海冰本构模型 海冰漂移的动力方程是按牛顿定律来描述的。单位面积海冰的运动平衡方程为 d Vi τ τ ξ ( σ) ( ) M = - Mf K ×V