个，文献一篇

[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
=
K_s*grad(h))

这里的K_s别直接取文献值，实测某红黏土的参数得先做三轴试验标定。

有个坑是当表面开始积水时，渗透系数会骤降两个数量级，这时候自适应网格必须跟上。

/>

收敛***方面，三点血泪经验：

1. 初始场别用零压力，拿稳态渗流场当初始条件
2. 时间步长用指数递增：0.1s→1s→10s→100s
3. 牛顿迭代次数上限调到50次，别信默认的20次

后处理要玩出花，试试在截面云图上叠加流线动画。

导出数据用Python处理更灵活：

import
matplotlib.pyplot
标记降雨1小时

当看到40mm/h工况下入渗率曲线在1小时后突然"躺平"，说明边坡表面已经开始哗哗地流汤了。

这时候的湿润锋推进速度比4mm/h工况快了近7倍，孔隙水压力监测点数值跟过山车似的上下颠簸。

/>

最后说个反直觉的现象：特大暴雨时边坡底部反而会出现负压区。

这是因为表层土壤渗透不及，形成临时隔水层，导致下部产生抽吸效应。

这种现象在传统达西定律模型中完全无法捕捉，必须上混合边界才行。