改进的SPH边界处理方法与土体大变形模拟

来源：计算力学学报杂志2018年第3期北京时间：

作者：骆钊,汪淳

单位：中交第二航务工程局有限公司, 武汉 430040,中交公路长大桥建设国家工程研究中心有限公司, 北京 250073,上海交通大学 船舶海洋与建筑工程学院, 上海 200240

摘要：沙土滑坡往往会造成重大的人身财产损失，研究这类土体大变形问题对防灾工程具有指导意义。光滑粒子流体动力学SPH（Smoothed Particle Hydrodynamics）方法是一种拉格朗日型无网格粒子法，十分适用于模拟大变形问题。在SPH方法中，合适的边界处理方法一直是个难点，传统的边界虚粒子法或排斥力法较难模拟复杂边界。本文引入了一种能处理任意形状边界的方法——统一半解析壁面边界条件处理方法USAW（unified semi-analytical wall boundary conditions），通过在控制方程中引入修正因子并保留边界面积分项来弥补边界缺失。为了更准确模拟问题域边界，提出无质量边界粒子的新概念。利用该方法成功模拟了土体滑坡算例，验证了方法的可靠性，并避免了边界零粒子层问题。通过数值模拟，分析了内摩擦角和黏聚力等土体物性参数对滑坡过程的影响。最后，应用该方法研究了滑坡冲击楔形体时的压力响应。

关键词：SPH,边界处理,土体,Drucker-Prager模型,大变形,

基金资助：国家重点基础研究计划（（973计划，2014CB046801）子课题）资助项目.