备案期间怎么关闭网站,企业网站建设代理公司,电脑清理优化大师,企业微信app下载安装官网二维钻孔封孔效果模拟案例钻孔封孔这事儿听着简单#xff0c;实际在地下工程里可是个技术活。今天咱们拿MATLAB的PDE工具箱做个二维模拟#xff0c;看看封孔材料怎么影响密封效果。先别急着关页面#xff0c;代码部分我尽量说得像唠嗑#xff0c;保证不催眠。先整点基础设置…二维钻孔封孔效果模拟案例钻孔封孔这事儿听着简单实际在地下工程里可是个技术活。今天咱们拿MATLAB的PDE工具箱做个二维模拟看看封孔材料怎么影响密封效果。先别急着关页面代码部分我尽量说得像唠嗑保证不催眠。先整点基础设置model createpde(); geometryFromEdges(model,circleg); pdegplot(model,EdgeLabels,on) axis equal这段就是创建个圆形钻孔模型circleg是自带的圆形几何模板。注意看生成的边界标签待会儿设置边界条件全靠这些Edge编号。我上次跑这个的时候手滑标错边界结果压力云图画得跟抽象画似的。材料参数得按实际情况来young 3e9; % 弹性模量3GPa poisson 0.25; % 泊松比 applyYoungModulus(model,young); applyPoissonsRatio(model,poisson);这里有个坑——弹性模量单位容易搞错。有次我把GPa输成MPa模拟出来的形变大的像橡皮泥封孔材料直接变橡皮糖效果。泊松比要是没概念的话普通混凝土在0.2左右橡胶材料能到0.49。边界条件设置是关键% 钻孔内壁施加注浆压力 applyPressure(model,2,1e6); % 2号边是内壁 % 外围固定约束 applyDisplacement(model,1,0,0); % 1号边是外围这里1e6帕相当于10个大气压实际操作中注浆压力得看具体设备。注意位移约束别乱加有次我手贱给内壁也加固定约束结果算出来的应力场直接突破天际线。二维钻孔封孔效果模拟案例生成网格时得看计算资源generateMesh(model,Hmax,0.1); mesh model.Mesh; pdeplot(model)Hmax控制网格密度数值越小网格越密。我笔记本跑0.05的网格得等杯咖啡时间工作站能跑到0.02。不过网格太密容易出锯齿状应力集中跟吃薯片掉渣似的。求解器跑起来之后result solve(model); pdeplot(model,XYData,result.VonMisesStress) colormap(jet) colorbar冯米塞斯应力云图一出立马能看出哪里的材料要扛不住了。上次模拟发现最大应力出现在封孔段端部跟现场检测的裂缝位置完全吻合甲方看了直呼内行。重点看压力传导% 提取钻孔径向位移 [ux,uy] evaluateStructuralResult(result); radialDisplacement ux.*cos(theta) uy.*sin(theta); plot(theta,radialDisplacement)径向位移曲线要是出现突变说明该处可能发生界面剥离。有次模拟曲线突然下凹到现场一看果然是注浆管位置留了个气泡空腔。优化参数时可以玩点花的paramSweep linspace(1e9,5e9,10); sealingEfficiency zeros(1,10); for i 1:10 model.MaterialProperties.YoungsModulus paramSweep(i); result solve(model); sealingEfficiency(i) max(result.VonMisesStress); end扫一遍弹性模量参数找应力极值点。有意思的是模量太大反而容易应力集中太小又hold不住压力跟找对象似的得取平衡点。最后说个实战经验模拟显示封孔材料外侧出现压力梯度反转云图颜色从蓝变黄再变蓝这可不是程序出bug。实际是材料塑性变形导致的应力重分布这时候得检查本构模型是不是该换成弹塑性了。下次碰到这种云图别急着重算先去现场取个岩芯看看说不定能发现意外的地质构造。