COMSOL冻土降雨水热力耦合：内有源文件与参考文献
冻土区工程遇上降雨就像开盲盒——你永远不知道第二天会不会出现融沉滑坡。去年在青藏高原做
路基监测的时候，亲眼见过一场暴雨让原本坚硬的冻土层直接变成"豆腐渣"，监测点的位移传感器直接报
警。这种水-热-力三场联动的鬼畜现象，用COMSOL建模可比现场抢险刺激多了。
先看核心控制方程，这玩意儿就像冻土界的"三体问题"。温度场带着达西流和应力场跳三人转：
```matlab
// 水分迁移方程
theta_w = theta_s*(1 + alpha*(T-T0))^(-1/m);
darcy_flux = -k*grad(h + z);
// 热传导方程带相变潜热
rho*C_p*dTdt = div(k_t*gradT) + L_f*rho_i*dphi_i/dt;
// 弹塑性本构方程
sigma = C : (epsilon - epsilon^p - epsilon^th);
```
特别是这个水分迁移方程里的theta_w，参数m取值能把人逼疯。做过冻融循环试验的都懂，当土体
含冰量超过30%时，渗透系数k会断崖式下跌两个数量级。有次算边坡稳定性，漏设了这个突变条件，结果算
出来的孔隙水压力比现实高出10倍——这要真按这个结果施工，估计现在还在吃牢饭。
耦合边界设置是重灾区。降雨入渗边界建议用Neumann条件配合事件接口：
```matlab
// 降雨边界条件
if t < 7200 // 2小时强降雨
q_water = 5e-5*sin(pi*t/3600); // 正弦波动入渗
else
q_water = 0; // 停止降雨后的蒸发阶段
end
// 热通量耦合
heat_flux = h_conv*(T_ext - T) + L_vapor*mass_flux; // 显热+潜热交换
```
有个坑爹案例：某冻土隧道建模时没考虑空气对流换热，算出来的衬砌温度比实测高8℃。后来发现
是h_conv参数用了默认值，实际高原风速下这个值应该放大3-5倍。改完参数后，温度场分布突然出现蝶形
分叉——原来排水盲管的位置产生了局部热岛效应。
