COMSOL模拟流固耦合井筒周围应力分布。 此案列介绍在井筒壁周围施加径向荷载(孔压和地应力），分析其径向应力、环向应力以及孔压变化，附有详细的建模说明书。
打开COMSOL新建模型，先选"多物理场耦合"里的达西流-固体力学接口。这里有个坑要注意：别忘了勾选"孔隙弹性"选项，不然流体和固体就各玩各的了。模型树里先定义砂岩的弹性模量（建议25GPa）和泊松比（0.25），渗透率参数用9e-13 m²比较符合常规储层。
材料参数设置代码片段： model.param.set('E_rock', '25e9[Pa]'); //岩石弹性模量 model.param.set('nu', 0.25); //泊松比 model.param.set('k0', '9e-13[m^2]'); //初始渗透率几何建模阶段，先画个半径0.1m的井眼（圆），外围扩展10米形成地层域。网格划分要讲究——井壁附近用边界层网格加密，远处用常规三角形网格。这里有个小技巧：在边界层设置里把厚度比例调到0.2，层数给8层，这样既能捕捉应力集中又不会卡死计算。
边界条件设置是重头戏。地应力用各向同性压力15MPa，井筒内压设置8MPa的梯度下降（生产时的压降过程）。流固耦合边界要特别注意：在井壁勾选"流体通量"和"位移约束"，选择压力边界类型为瞬态变化。
边界条件核心代码： //井壁压力边界 model.physics('dl').feature('pr1').set('p0', '8e6[Pa]*(1-exp(-t/3600))'); //地应力加载 model.physics('solid').feature('bnd2').set('loadType', 'Pressure'); model.physics('solid').feature('bnd2').set('p0', '15e6[Pa]');求解器设置建议用瞬态分析，时间步长用对数刻度：前10秒用0.1秒步长，之后逐步放大到100秒步长。计算完成后别急着关软件，先检查位移云图是否出现"蝴蝶形"应力分布——这是判断耦合是否成功的标志。
后处理阶段要玩转截面绘图：在径向方向切出0°、45°、90°三个剖面线，对比不同方向的环向应力。有个隐藏功能可能你不知道——在结果里添加"最大剪切应力"变量，用公式sqrt((sr-sθ)^2/4 + tau_rθ^2)，能更直观看到潜在破裂区。
后处理变量定义： //最大剪切应力自定义变量 model.result().dataset("cpt1").feature("expr1").set("expr", "sqrt(((solid.sr-solid.st)^2)/4 + solid.trt^2)");最后导出数据时，建议把井壁周围0.5米范围内的数据单独存为文本文件。用MATLAB或Python做个极坐标图，比软件自带的图表更专业。注意看环向应力的"马鞍形"分布，这个特征能说明地应力是否达到岩石抗张强度。
整个模型跑下来大概需要20分钟（配置i7+32G内存），如果发现应力云图出现不连续跳变，八成是网格在井壁附近不够密。下次可以试试在物理场控制里勾选"自适应网格"，让软件自动加密关键区域，效果立竿见影。
