一种判断混凝土泵送过程中堵管风险性的方法与流程

本发明属于建筑工程中混凝土施工技术领域，特别涉及一种判断新拌混凝土在泵送过程中堵管风险性的方法。

背景技术：

目前我国工程建设的规模和技术难度居于全球领先水平，各种超高、超长、超深、超大型结构不断推陈出新，例如，世界前10高度的超高层建筑有6座在中国，因此对混凝土施工性能、力学性能和耐久性提出了更高的要求，c60及以上等级的高强混凝土在高层建筑、大跨径桥梁、大型公共建筑等工程中广泛应用，并且混凝土的主要施工方式是泵送，在混凝土泵的压力推动下沿输送管道将混凝土运输至施工面。因此研究泵送高强混凝土的施工稳定性具有重要的理论价值和工程实践意义。

然而，工程实践表明泵送高强混凝土存在一系列问题亟需解决，突出表现在：高强混凝土粘度大、流动性差，导致泵送压力大、泵送施工困难，容易出现堵管现象，给施工企业造成重大损失。已有经验评价方法难以全面预测和评价高强混凝土的泵送性能与可泵性。因此，亟待提出能够准确判断混凝土泵送过程中堵管风险性的方法，解决泵送高强混凝土面临着泵送性能控制困难的问题。

技术实现要素：

目前高强混凝土粘度大、流动性差，导致泵送压力大、泵送施工困难，已有经验评价方法难以全面预测和评价高强混凝土的可泵性。本发明提供了一种判断混凝土泵送过程中堵管风险性的方法，通过基于cfd+dem耦合的形式，建立分析新拌混凝土流变机制和泵送性能的数值仿真模拟方法，预测与判断混凝土泵送过程中堵管风险性。

本发明解决所采用的技术方案是：

步骤一、针对不同强度等级与不同稠度指数的新拌混凝土，分别检测它们的稠度指数k1、雷诺数re、屈服应力σ0和幂率指数n；

步骤二、针对不同强度等级与不同稠度指数的新拌混凝土，确定该次进行cfd+dem耦合模拟的新拌混凝土的精度n。其中cfd是采用基于latticeboltzmannmethod(lbm)方法的开源软件平台palabos，dem是采用开源软件平台liggghts。在linux服务器上部署搭建以上两个软件的工作空间，使其可以正常工作,并可以相互调用后台代码从而为cfd+dem的耦合创建编译运行环境；

步骤三、针对不同强度等级与不同稠度指数的新拌混凝土，以herschel-bulkley本构模型σ＝σ0+k1γ^w为基础，建立新拌混凝土应力-应变率之间的关系，导入精度n、屈服应力σ0、雷诺数re、幂率指数n和稠度指数k1；

步骤四、根据泵管尺寸与泵送高度h，采用开源软件palabos建立3d泵送管道模型。应用cfd+dem耦合模拟混凝土在泵管中的泵送状态；用palabos模拟砂浆，掺入由软件liggghts模拟的石子，使其均匀离散化分布在砂浆中。

步骤五、通过流体动力学牛顿迭代算法，模拟分析石子-砂浆固液两相混合流相互作用，实现泵送混凝土的过程。在每一步迭代运算中，通过调用内置的获取ω的函数可计算出局部网格处的松弛因子ω，再将松弛因子作为实际参数,根据内置bgk动力学类创建bgk动力学模型对象。再通过开源软件内置多块3d类创建网格对象，其实际参数为x，y，z三轴方向上的格子数，以及上一步中创建的bgk动力学对象。格子对象调用内置的局部密度分布函数和局部速度场函数，均以x，y，z三轴方向上的格子数为实际参数。局部速度场与宏观设置的平均速度场进行差值运算，当差值逐渐接近零时，说明计算结果收敛并导出内部石子和砂浆的空间位置场和速度场，便可通过这两个场云图直观地形成一种基于cfd+dem耦合后模拟得到的评价混凝土泵送顺畅性(即堵管风险性)的方法。具体判别标准如下：

1)当石子的某一时刻的空间分布密集程度不均匀，出现局部空间位置混凝土密度骤减，即小于整体平均密度的15％时，该新拌混凝土会发生堵管；

2)当石子速度场云图中某位置速度骤减，不再保持有限范围内的波动，即局部速度减小超过50％时，该新拌混凝土会发生堵管；

3)若新拌混凝土的速度场云图收敛后，即局部速度场u与宏观设置的平均速度场u进行差值运算，当|u-u|逐渐趋近于0时，该新拌混凝土适宜泵送。

与现有技术相比，本发明具有以下优点：

(1)本发明提供的一种判断混凝土泵送过程中堵管风险性的方法，不需要进行实际泵送施工检测，只需要通过计算机模拟导出泵管中混凝土内部石子与砂浆的速度场和空间场云图，就能够判断泵送过程中堵管的风险性，该方法科学有效、方便快捷且经济性好、实用性好，能有效满足工程实际需要；

(2)本发明可以多次模拟得到各种新拌混凝土的综合成本、稠度指数、雷诺数等参数，进而拟合出各种混凝土综合成本与稠度指数、雷诺数之间的关系曲线，确定适合泵送的混凝土稠度指数、雷诺数等关键指标上限值，避免混凝土实际泵送施工时泵机选型错误、设备能耗大等问题，从而大幅提高混凝土泵送施工安全性。

附图说明

图1是实施例1的空间分布密度云图：

图2是实施例2的空间速度场云图：

图3是实施例3的空间速度场云图：

具体实施方式一：

一种判断混凝土泵送过程中堵管风险性的方法，下面结合具体实施方式，对本方法作进一步详细说明：

步骤一，一种新拌混凝土，检测其稠度指数k1＝2.0、雷诺数re＝2000、屈服应力σ0＝180pa和幂率指数n＝1.8；

步骤二，确定该次进行cfd+dem耦合模拟的新拌混凝土的精度n＝100。其中cfd是采用基于latticeboltzmannmethod(lbm)方法的开源软件平台palabos，dem是采用开源软件平台liggghts。在linux服务器上部署搭建以上两个软件的工作空间，使其可以正常工作,并可以相互调用后台代码从而为cfd+dem的耦合创建编译运行环境；

步骤三，以herschel-bulkley本构模型σ＝σ0+k1γⁿ为基础，建立新拌混凝土应力-应变率之间的关系，导入精度n＝100、屈服应力σ0＝180pa、雷诺数re＝2000、幂率指数n＝1.8和稠度指数k1＝1.8；

步骤四，根据泵管尺寸与泵送高度h＝500m，采用开源软件palabos建立3d泵送管道模型。应用cfd+dem耦合模拟混凝土在泵管中的泵送状态；用palabos模拟砂浆，掺入由软件liggghts模拟的石子，使其均匀离散化分布在砂浆中。

步骤五，通过流体动力学牛顿迭代算法，模拟分析石子-砂浆固液两相混合流相互作用，实现泵送混凝土的过程。在每一步迭代运算中，通过调用内置的获取ω的函数可计算出局部网格处的松弛因子ω，再将松弛因子作为实际参数,根据内置bgk动力学类创建bgk动力学模型对象。再通过开源软件内置多块3d类创建网格对象，其实际参数为x，y，z三轴方向上的格子数，以及上一步中创建的bgk动力学对象。格子对象调用内置的局部密度分布函数和局部速度场函数，均以x，y，z三轴方向上的格子数为实际参数。局部速度场与宏观设置的平均速度场进行差值运算，当差值逐渐接近零时，说明计算结果收敛并导出内部石子和砂浆的空间位置场和速度场，便可通过这两个场云图直观地形成一种基于cfd+dem耦合后模拟得到的评价混凝土泵送顺畅性(即堵管风险性)的方法。

由图1可发现上部空间位置混凝土密度较小，比整体平均密度低16.54％，因此判断该新拌混凝土会发生堵管。

具体实施方式二：

一种判断混凝土泵送过程中堵管风险性的方法，下面结合具体实施方式，对本方法作进一步详细说明：

步骤一，一种新拌混凝土，检测其稠度指数k1＝1.8、雷诺数re＝2300、屈服应力σ0＝190pa和幂率指数n＝1.9；

步骤二，确定该次进行cfd+dem耦合模拟的新拌混凝土的精度n＝100。其中cfd是采用基于latticeboltzmannmethod(lbm)方法的开源软件平台palabos，dem是采用开源软件平台liggghts。在linux服务器上部署搭建以上两个软件的工作空间，使其可以正常工作,并可以相互调用后台代码从而为cfd+dem的耦合创建编译运行环境；

步骤三，以herschel-bulkley本构模型σ＝σ0+k1γⁿ为基础，建立新拌混凝土应力-应变率之间的关系，导入精度n＝100、屈服应力σ0＝190pa、雷诺数re＝2300、幂率指数n＝1.9和稠度指数k1＝1.8；

步骤四，根据泵管尺寸与泵送高度h＝600m，采用开源软件palabos建立3d泵送管道模型。应用cfd+dem耦合模拟混凝土在泵管中的泵送状态；用palabos模拟砂浆，掺入由软件liggghts模拟的石子，使其均匀离散化分布在砂浆中。

步骤五，通过流体动力学牛顿迭代算法，模拟分析石子-砂浆固液两相混合流相互作用，实现泵送混凝土的过程。在每一步迭代运算中，通过调用内置的获取ω的函数可计算出局部网格处的松弛因子ω，再将松弛因子作为实际参数,根据内置bgk动力学类创建bgk动力学模型对象。再通过开源软件内置多块3d类创建网格对象，其实际参数为x，y，z三轴方向上的格子数，以及上一步中创建的bgk动力学对象。格子对象调用内置的局部密度分布函数和局部速度场函数，均以x，y，z三轴方向上的格子数为实际参数。局部速度场与宏观设置的平均速度场进行差值运算，当差值逐渐接近零时，说明计算结果收敛并导出内部石子和砂浆的空间位置场和速度场，便可通过这两个场云图直观地形成一种基于cfd+dem耦合后模拟得到的评价混凝土泵送顺畅性(即堵管风险性)的方法。

由图2可发现上部速度场云图中速度骤减，局部速度减小了94.11％，判断该新拌混凝土会发生堵管。

具体实施方式三：

一种判断混凝土泵送过程中堵管风险性的方法，下面结合具体实施方式，对本方法作进一步详细说明：

步骤一，一种新拌混凝土，检测其稠度指数k1＝2.1、雷诺数re＝2100、屈服应力σ0＝170pa和幂率指数n＝1.7；

步骤二，确定该次进行cfd+dem耦合模拟的新拌混凝土的精度n＝100。其中cfd是采用基于latticeboltzmannmethod(lbm)方法的开源软件平台palabos，dem是采用开源软件平台liggghts。在linux服务器上部署搭建以上两个软件的工作空间，使其可以正常工作,并可以相互调用后台代码从而为cfd+dem的耦合创建编译运行环境；

步骤三，以herschel-bulkley本构模型σ＝σ0+k1γⁿ为基础，建立新拌混凝土应力-应变率之间的关系，导入精度n＝100、屈服应力σ0＝170pa、雷诺数re＝2100、幂率指数n＝1.7和稠度指数k1＝2.1；

步骤四，根据泵管尺寸与泵送高度h＝450m，采用开源软件palabos建立3d泵送管道模型。应用cfd+dem耦合模拟混凝土在泵管中的泵送状态；用palabos模拟砂浆，掺入由软件liggghts模拟的石子，使其均匀离散化分布在砂浆中。

步骤五，通过流体动力学牛顿迭代算法，模拟分析石子-砂浆固液两相混合流相互作用，实现泵送混凝土的过程。在每一步迭代运算中，通过调用内置的获取ω的函数可计算出局部网格处的松弛因子ω，再将松弛因子作为实际参数,根据内置bgk动力学类创建bgk动力学模型对象。再通过开源软件内置多块3d类创建网格对象，其实际参数为x，y，z三轴方向上的格子数，以及上一步中创建的bgk动力学对象。格子对象调用内置的局部密度分布函数和局部速度场函数，均以x，y，z三轴方向上的格子数为实际参数。局部速度场与宏观设置的平均速度场进行差值运算，当差值逐渐接近零时，说明计算结果收敛并导出内部石子和砂浆的空间位置场和速度场，便可通过这两个场云图直观地形成一种基于cfd+dem耦合后模拟得到的评价混凝土泵送顺畅性(即堵管风险性)的方法。

由图3可发现新拌混凝土的速度场云图收敛，局部速度场u与宏观设置的平均速度场u进行差值运算，|u-u|＝0.0000154，判断该新拌混凝土适宜泵送。

综上所述，本发明提供了一种科学、准确、有效、经济的评价新拌混凝土泵送顺畅性的预测方法，该方法基于cfd+dem耦合的原理，参考稠度指数、雷诺数、屈服应力、幂率指数n、模拟精度n等控制指标，开发了一种有效且经济性好的判断混凝土泵送过程中堵管风险性的方法。最大程度地避免了泵送施工风险，降低了成本，能够满足现代混凝土泵送施工需求，显著降低堵管、爆管等事故的发生。