边坡降雨入渗问题中两种边界条件的处理及应用 [1]模型简介：使用数值模拟软件COMSOL，复现...
边坡降雨入渗问题中两种边界条件的处理及应用 [1]模型简介：使用数值模拟软件COMSOL，复现论文（年庚乾,陈忠辉,张凌凡等.边坡降雨入渗问题中两种边界条件的处理及应用[J].岩土力学,2020,41(12):4105-4115.），建立二维边坡模型，应用流量—压力混合入渗边界控制方程，分析了不同降雨强度（4mm/h、40mm/h）下边坡降雨入渗及渗流规律 [2]案例内容：边坡降雨入渗完整数值模型一个（包括边界条件、云图、后处理结果），DXF二维模型一个，文献一篇 [3]模型特色：掌握降雨流量—压力混合入渗边界及渗流边界的处理，掌握模型计算收敛性技巧，锻炼后处理及入渗率、入渗量曲线作图
边坡降雨入渗的数值模拟总让人头疼，尤其是边界条件的处理。最近用COMSOL复现了年庚乾团队提出的流量-压力混合边界模型，这里分享些实战经验。先丢个结论：降雨强度从4mm/h飙到40mm/h时，边坡的渗流场会从温柔渗透直接变成狂暴模式。
建模先从几何开始。建议直接用AutoCAD画好边坡的DXF文件导入，比在COMSOL里手动画省事得多。注意边坡角控制在30-45度之间，坡高15米左右比较符合常规工程场景。导入时别忘了检查单位设置——我上次手滑把毫米当米导入，结果算出来的孔隙水压力直接突破天际。
边坡降雨入渗问题中两种边界条件的处理及应用 [1]模型简介：使用数值模拟软件COMSOL，复现论文（年庚乾,陈忠辉,张凌凡等.边坡降雨入渗问题中两种边界条件的处理及应用[J].岩土力学,2020,41(12):4105-4115.），建立二维边坡模型，应用流量—压力混合入渗边界控制方程，分析了不同降雨强度（4mm/h、40mm/h）下边坡降雨入渗及渗流规律 [2]案例内容：边坡降雨入渗完整数值模型一个（包括边界条件、云图、后处理结果），DXF二维模型一个，文献一篇 [3]模型特色：掌握降雨流量—压力混合入渗边界及渗流边界的处理，掌握模型计算收敛性技巧，锻炼后处理及入渗率、入渗量曲线作图
核心代码藏在"多物理场"模块里：
if 降雨强度 < 饱和渗透系数: 应用Dirichlet边界（压力固定） else: 切换为Neumann边界（流量恒定）这个逻辑在COMSOL里需要写成弱形式PDE。重点是用阶跃函数平滑过渡，否则计算必崩。记得在方程里加个松弛因子，类似：
theta_new = theta_old + dt*(q_rain - K_s*grad(h))这里的K_s别直接取文献值，实测某红黏土的参数得先做三轴试验标定。有个坑是当表面开始积水时，渗透系数会骤降两个数量级，这时候自适应网格必须跟上。
收敛性技巧方面，三点血泪经验：
- 初始场别用零压力，拿稳态渗流场当初始条件
- 时间步长用指数递增：0.1s→1s→10s→100s
- 牛顿迭代次数上限调到50次，别信默认的20次
后处理要玩出花，试试在截面云图上叠加流线动画。导出数据用Python处理更灵活：
import matplotlib.pyplot as plt plt.plot(t, cumulative_infiltration) plt.axvline(x=3600, color='r', linestyle='--') # 标记降雨1小时当看到40mm/h工况下入渗率曲线在1小时后突然"躺平"，说明边坡表面已经开始哗哗地流汤了。这时候的湿润锋推进速度比4mm/h工况快了近7倍，孔隙水压力监测点数值跟过山车似的上下颠簸。
最后说个反直觉的现象：特大暴雨时边坡底部反而会出现负压区。这是因为表层土壤渗透不及，形成临时隔水层，导致下部产生抽吸效应。这种现象在传统达西定律模型中完全无法捕捉，必须上混合边界才行。