一种泥浆消纳植桩承载力分析方法及系统

1.一种泥浆消纳植桩承载力分析方法，其特征在于，包括：获取地质勘测数据，基于所述地质勘测数据构建泥浆消纳植桩区域的三维地质模型，其中，所述三维地质模型包括泥浆、土体和植桩在内的三维复杂模型；将所述三维地质模型离散化，得到多个模型小单元；确定各模型小单元对应的材料特性和初始条件，并基于所述材料特性和初始条件对模型小单元进行多尺度数值模拟计算，以对各模型小单元求解，得到承载力信息；确定各模型小单元对应的材料特性和初始条件，并基于所述材料特性和初始条件对模型小单元进行多尺度数值模拟计算，以对各模型小单元求解，得到承载力信息的步骤，包括：利用有限元分析和流体力学对所述模型小单元进行分析，其中，在有限元分析中，每个模型单元将被分配材料特性和初始条件，然后根据所选择的有限元方法进行离散求解，通过求解局部微分方程和施加边界条件来计算出消纳植桩内部的应力、应变和变形情况，其中，利用有限元分析和流体力学对所述模型小单元进行分析的步骤，包括：利用流体力学分析所述模型小单元时，利用CFD模拟泥浆在土体孔隙中的流动和扩散过程，包括根据纳维-斯托克斯方程和质量守恒方程确定泥浆在孔隙中的流动情况，其中，泥浆在孔隙中的流动情况表达式如下：利用流体力学分析所述模型小单元时，利用CFD模拟泥浆在土体孔隙中的流动和扩散过程，包括根据纳维-斯托克斯方程和质量守恒方程确定泥浆在孔隙中的流动情况，其中，泥浆在孔隙中的流动情况表达式如下： ；其中， 为泥浆密度， 为速度场， 为压力场，H为泥浆粘度， 为重力加速度；引入古代模型修改纳维-斯托克斯方程，以描述泥浆在土体孔隙中的流动和扩散过程，其中，古代模型将流体的粘度与剪切速率关联，其表达式为： ；其中， 是流体的动力粘度， 是流体的流动指数， 是古代指数， 是速度场的剪切率；然后，将固体颗粒视为另一种相，并在流体相和颗粒相之间建立耦合关系，得到修改后的纳维-斯托克斯方程如下所示： ；其中， 是泥浆密度， 是速度场， 是压力场， 是流体的动力粘度， 是速度梯度张量， 是重力加速度；另外，颗粒相的运动通过牛顿第二定律来描述的表达式如下： ；其中， 是颗粒的质量，i是颗粒的速度， 是由流体对颗粒的拖曳力， 是重力对颗粒的作用力；其次，利用流体力学分析所述模型小单元时，离散化的模型单元将被用来表示流体域的空间网格，速度、压力和其他流体参数在这些模型单元之间进行插值，并在每个模型单元内求解流体流动和传输方程，以模拟泥浆在土体孔隙中的流动和扩散过程。 2.根据权利要求1所述的泥浆消纳植桩承载力分析方法，其特征在于，利用有限元分析和流体力学对所述模型小单元进行分析的步骤，包括：利用有限元法模拟植桩周围土体的应力和变形情况，土体中任意点的位移u可以通过求解以下弹性力学方程得到： ；其中， 为应力张量，则 、 、 分别表示在土地中浆消纳植桩的某一点在x、y、z方向上的应力分量， 为应变张量，则 、 、 分别表示在土地中浆消纳植桩的某一点在x、y、z方向上的应变分量， 为弹性模量， 为泊松比， ， ， 分别为x方向、y方向、z方向上的分量。 3.根据权利要求2所述的泥浆消纳植桩承载力分析方法，其特征在于，利用有限元法模拟植桩周围土体的应力和变形情况的步骤之后，还包括：引入非线性效应和材料非均质性，其中，Mohr-Coulomb本构模型的应力-应变关系可以表示为： ；其中， 为应力张量， 为剪切模量， 1 为体积模量， 为等效应变， 为体积应变；在这个模型中，等效应变 通常根据von Mises屈服准则计算： ；其中， 为剪应力张量的偏差张量， 表示剪切模量，表示了土体抗剪切变形的能力， 1 为体积模量，表示了土体的抗压性能， 表示等效应变，代表土体的塑性变形程度。 4.根据权利要求3所述的泥浆消纳植桩承载力分析方法，其特征在于，还包括：用各向异性本构模型，并引入动态响应方程,以考虑到土体的各向异性和动态响应；对于各向异性的本构模型，使用Transversely Isotropic Mohr-Coulomb（TIMC）模型，以表示土体在径向和轴向上的不同力学性质，该模型的应力-应变关系可以表示为： ；其中， 为应力张量， 为应变张量， 为弹性刚度张量，其表达式为： ；各向异性的材料参数表示如下： ， ， ：表示在径向的弹性模量， ， ， ：表示在径向和轴向之间的剪切模量， ， ， ：表示径向和轴向的Poisson比；对于动态响应方程，本发明实施例可以考虑采用地震工程中常用的弹性波方程： ；其中，R为位移场， 为波速， 表示拉普拉斯算子。 5.一种泥浆消纳植桩承载力分析系统，其特征在于，包括：获取模块，用于获取地质勘测数据，基于所述地质勘测数据构建泥浆消纳植桩区域的三维地质模型，其中，所述三维地质模型包括泥浆、土体和植桩在内的三维复杂模型；离散模块，用于将所述三维地质模型离散化，得到多个模型小单元；确定模块，用于确定各模型小单元对应的材料特性和初始条件，并基于所述材料特性和初始条件对模型小单元进行多尺度数值模拟计算，以对各模型小单元求解，得到承载力信息；所述确定模块用于执行以下步骤：利用有限元分析和流体力学对所述模型小单元进行分析，其中，在有限元分析中，每个模型单元将被分配材料特性和初始条件，然后根据所选择的有限元方法进行离散求解，通过求解局部微分方程和施加边界条件来计算出消纳植桩内部的应力、应变和变形情况，其中，利用有限元分析和流体力学对所述模型小单元进行分析的步骤，包括：利用流体力学分析所述模型小单元时，利用CFD模拟泥浆在土体孔隙中的流动和扩散过程，包括根据纳维-斯托克斯方程和质量守恒方程确定泥浆在孔隙中的流动情况，其中，泥浆在孔隙中的流动情况表达式如下：利用流体力学分析所述模型小单元时，利用CFD模拟泥浆在土体孔隙中的流动和扩散过程，包括根据纳维-斯托克斯方程和质量守恒方程确定泥浆在孔隙中的流动情况，其中，泥浆在孔隙中的流动情况表达式如下： ；其中， 为泥浆密度， 为速度场， 为压力场，H为泥浆粘度， 为重力加速度；引入古代模型修改纳维-斯托克斯方程，以描述泥浆在土体孔隙中的流动和扩散过程，其中，古代模型将流体的粘度与剪切速率关联，其表达式为： ；其中， 是流体的动力粘度， 是流体的流动指数， 是古代指数， 是速度场的剪切率；然后，将固体颗粒视为另一种相，并在流体相和颗粒相之间建立耦合关系，得到修改后的纳维-斯托克斯方程如下所示： ；其中， 是泥浆密度， 是速度场， 是压力场， 是流体的动力粘度， 是速度梯度张量， 是重力加速度；另外，颗粒相的运动通过牛顿第二定律来描述的表达式如下： ；其中， 是颗粒的质量，i是颗粒的速度， 是由流体对颗粒的拖曳力， 是重力对颗粒的作用力；其次，利用流体力学分析所述模型小单元时，离散化的模型单元将被用来表示流体域的空间网格，速度、压力和其他流体参数在这些模型单元之间进行插值，并在每个模型单元内求解流体流动和传输方程，以模拟泥浆在土体孔隙中的流动和扩散过程。 6.一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行权利要求1到4任一所述的泥浆消纳植桩承载力分析方法。 7.一种可读存储介质，其特征在于，所述可读存储介质中存储有计算机程序，所述计算机程序包括用于控制过程以执行过程的程序代码，所述过程包括根据权利要求1到4任一所述的泥浆消纳植桩承载力分析方法。