一种垃圾填埋场失稳流滑灾害数值仿真方法技术

本发明专利技术涉及一种移动半隐式粒子法（简称MPS）的垃圾填埋失稳流滑灾害的数值仿真方法。针对原始MPS方法存在严重的压力震荡，从三个方面对源程序进行改进，包括：核函数、Poisson压力方程以及自由表面判别。为刻画流滑灾害体的非牛顿流变特性，创新性的引入宾汉姆非牛顿流体本构关系。建立等效粘度概念，将宾汉姆流体本构转化为变粘度的牛顿流体本构，并编程实现。通过平衡压力震荡和引入宾汉姆非牛顿流体本构关系发明专利技术了改进的MPS方法。该发明专利技术能够精确仿真垃圾填埋场失稳流滑的动力过程。

【技术实现步骤摘要】

本专利技术涉及一种垃圾填埋场失稳流滑灾害数值仿真方法，属于岩土计算力学、流体力学与环境科学

技术介绍
垃圾填埋场在我国处理城市固体废弃物中广泛使用。城市固体废弃物在填埋环境作用下，在降雨、地震、渗滤液回灌等外部因素的影响下，容易液化失稳。液化使原有结构遭到破坏，并形成新的流动结构，最终演化为流滑灾害事故，造成废弃物、渗滤液及发酵气体大范围扩散，严重污染城市地表及地下环境。因此分析研究垃圾填埋场失稳流滑灾害具有重大的经济效益和环境效益。传统基于网格的数值仿真方法处理流滑等大变形灾害时，在不规则区域网格布置、自由表面追踪等方面面临极大挑战与困难。以MPS为代表的粒子法，通过一系列粒子离散问题域，能够克服网格类方法面临的网格扭曲、畸变等问题。因此，采用MPS方法开展垃圾填埋产失稳流滑灾害的数值仿真具有极大的优势。然而传统的MPS方法存在极大的压力震荡，影响数值仿真的准确性。同时失稳流滑的垃圾体流动已超出牛顿流体的范畴，急需发展一种能准确刻画该现象的本构关系，以捕捉其流体动力特性，为垃圾填埋场的规划选址提供有力的科学依据。
技术实现思路
本专利技术的目的在于提供一种可用于准确仿真垃圾填埋场失稳流滑灾害的数值方法。为达到上述目的，本专利技术采取了一下技术方案：一种垃圾填埋场失稳流滑灾害数值仿真方法，具体步骤如下：1)确定垃圾填埋场的研究区域及物理参数，所述物理参数包括填埋场尺寸、>垃圾体密度及强度参数；2)根据步骤1)确定的垃圾填埋场的研究区域计算粒子的初始分布；3)根据步骤2)得到的粒子的初始分布考虑粘性力项和重力项以计算粒子的临时速度；所述粘性力项采用下式(1)宾汉姆流体本构模型计算得到； τ = η γ · + τ m i n τ > τ min γ · = 0 τ ≤ τ min - - - ( 1 ) ]]>其中：τ为剪应力，η是粘度系数，为应变速率，τmin为屈服剪切强度；为考虑垃圾土的特性，将摩尔-库伦准则引入上述宾汉姆流体本构模型以界定垃圾填埋场失稳流滑的初始状态，则得到屈服剪切强度τmin如下式(2)：τmin＝c+Ptanφ(2)c表示垃圾土粘聚力，P表示垃圾土正应力，φ表示垃圾土内摩擦角。同时为简化数值计算，将等效粘度引入η′(如下式(3))引入宾汉姆流体本构模型得到修正后的宾汉姆流体本构模型如下式(4)； η ′ = η + c + P t a n φ | γ · | ( 1 - e - m γ · ) - - - ( 3 ) ]]> τ = η ′ γ · - - - ( 4 ) ]]>其中m为压力增长系数；4)根据步骤3)得到的粒子的临时速度计算粒子的临时位置；5)根据步骤4)得到的粒子的临时位置计算粒子的临时粒子数密度ni，公式如下式(5)： n i = Σ j ≠ i w ( | r j - r i | ) - - - ( 5 ) ]]>其中：rj和ri代表粒子j和i的坐标，w(r)表示核函数，r表示粒子间距离；由于传统MPS采用的核函数，当r＝0时，取值为无穷大，会导致压力出现震荡。为平衡压力震荡，采用下式(6)所属的核函数代替传统MPS采用的核函数； 本文档来自技高网...

【技术保护点】
一种垃圾填埋场失稳流滑灾害数值仿真方法，其特征在于具体步骤如下：1)确定垃圾填埋场的研究区域及物理参数，所述物理参数包括填埋场尺寸、垃圾体密度及强度参数；2)根据步骤1)确定的垃圾填埋场的研究区域计算粒子的初始分布；3)根据步骤2)得到的粒子的初始分布考虑粘性力项和重力项以计算粒子的临时速度；所述粘性力项采用下式(1)宾汉姆流体本构模型计算得到；τ=ηγ·+τminτ>τminγ·=0τ≤τmin---(1)]]>其中：τ为剪应力，η是粘度系数，为应变速率，τmin为屈服剪切强度；为考虑垃圾土的特性，将摩尔‑库伦准则引入上述宾汉姆流体本构模型以界定垃圾填埋场失稳流滑的初始状态，则得到屈服剪切强度τmin如下式(2)：τmin＝c+P tanφ  (2)c表示垃圾土粘聚力，P表示垃圾土正应力，φ表示垃圾土内摩擦角；同时为简化数值计算，将等效粘度引入η′(如下式(3))引入宾汉姆流体本构模型得到修正后的宾汉姆流体本构模型如下式(4)；η′=η+c+Ptanφ|γ·|(1-e-mγ·)---(3)]]>τ=η′γ·---(4)]]>其中m为压力增长系数；4)根据步骤3)得到的粒子的临时速度计算粒子的临时位置；5)根据步骤4)得到的粒子的临时位置计算粒子的临时粒子数密度ni，公式如下式(5)：ni=Σj≠iw(|rj-ri|)---(5)]]>其中：rj和ri代表粒子j和i的坐标，w(r)表示核函数，r表示粒子间距离；由于传统MPS采用的核函数，当r＝0时，取值为无穷大，会导致压力出现震荡；为平衡压力震荡，采用下式(6)所述的核函数代替传统MPS采用的核函数；w(r)=(1-rre)3(1+rre)3(0≤r<re)0(r≥re)---(6)]]>re为影响域半径；6)根据步骤5)得到的粒子的临时粒子数密度建立Poisson压力方程；采用包含一个主项和两个修正项的Poisson压力方程(下式(7))以平衡压力震荡：1ρ▿2Pk+1=n*-2nk+nk-1n0Δt2+aΔtnk-nk-1n0Δt+bΔt2nk-n0n0---(7)]]>其中：ρ表示垃圾体密度，▽表示梯度算子，Pk+1表示第(k+1)步的压力，n*表示临时粒子数密度，nk表示第k步的粒子数密度，nk‑1表示第(k‑1)步的粒子数密度，n0表示初始粒子数密度，△t表示时间增量，a和b为误差项修正系数；7)根据步骤6)求解Poisson压力方程得到压力场分布；8)采用自由表面判别公式(8)确定自由表面粒子，并将自由表面粒子的压力置为零，得到修改后的压力场分布；n*<βn0     (8)n*表示临时粒子数密度，n0表示初始粒子数密度，β为系数；9)根据步骤8)得到的修改后的压力场分布，计算压力梯度修正粒子速度及位置；10)根据步骤9)，判断是否达到预设计算时间，若达到即停止计算，否则重复步骤3)‑9)，直至达到预设计算时间。...

【技术特征摘要】
1.一种垃圾填埋场失稳流滑灾害数值仿真方法，其特征在于具体步骤如下：
1)确定垃圾填埋场的研究区域及物理参数，所述物理参数包括填埋场尺寸、垃
圾体密度及强度参数；
2)根据步骤1)确定的垃圾填埋场的研究区域计算粒子的初始分布；
3)根据步骤2)得到的粒子的初始分布考虑粘性力项和重力项以计算粒子的临
时速度；所述粘性力项采用下式(1)宾汉姆流体本构模型计算得到；
τ = η γ · + τ m i n τ > τ min γ · = 0 τ ≤ τ min - - - ( 1 ) ]]>其中：τ为剪应力，η是粘度系数，为应变速率，τmin为屈服剪切强度；为考
虑垃圾土的特性，将摩尔-库伦准则引入上述宾汉姆流体本构模型以界定垃圾填
埋场失稳流滑的初始状态，则得到屈服剪切强度τmin如下式(2)：
τmin＝c+Ptanφ(2)
c表示垃圾土粘聚力，P表示垃圾土正应力，φ表示垃圾土内摩擦角；同时为简
化数值计算，将等效粘度引入η′(如下式(3))引入宾汉姆流体本构模型得到
修正后的宾汉姆流体本构模型如下式(4)；
η ′ = η + c + P t a n φ | γ · | ( 1 - e - m γ · ) - - - ( 3 ) ]]> τ = η ′ γ · - - - ( 4 ) ]]>其中m为压力增长系数；
4)根据步骤3)得到的粒子的临时速度计算粒子的临时位置；
5)根据步骤4)得到的粒子的临时位置计算粒子的临时粒子数密度ni，公式如
下式(5)：
n i = Σ j ≠ i w ( | r j - r i | ) - - - ( 5 ) ]]>其中：rj和ri代表粒子j和i的坐标，w(r)表示核函数，r表示粒子间距离；由于

\t传统MPS采用的核函数，当r＝0时，取值为无穷大，会导致压力出现震荡；为
平衡压力震荡，采用下式(6)所述的核函数代替传统MPS采用的核函数；
w ( r ) = ( 1 - ...

【专利技术属性】
技术研发人员：黄雨，朱崇强，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海;31