专利名称：一种岩土材料流滑灾变的计算模拟方法
技术领域：
本发明涉及一种岩土材料流滑灾变的计算模拟方法，属于岩土计算力学、地质灾害防治和地质环境保护技术领域。
背景技术：
近年来，在地震、降雨等自然因素和人类活动的综合作用下，地质灾害的发生频率越来越高，对人类生命财产造成的危害也越来越严重。这些地质灾害，如地震诱发的高速远程滑坡、土体液化后的侧向流动、泥石流、填埋体失稳流滑等，从本质上来说都属于岩土材料的大变形流动破坏问题。要揭示这些地质灾害发生及发展的动力学特征并提出相应的防治措施，就必须对岩土材料的复杂流动行为有充分的认识。现有的计算方法都属于基于网格划分的固体力学计算方法，无法处理岩土材料失稳流滑所涉及的极大变形、高度非线性以及多相耦合等科学难题，无法实现岩土材料流滑灾变的全过程模拟。因此，迫切需求一种切实有效的计算方法来模拟岩土材料的流滑灾变过程，捕捉其流体动力学特征，为流滑灾害高风险地区的工程设计、工程建设和防灾减灾等提供有力的科学依据。
发明内容
本发明的目的是针对当前高速远程滑坡、液化土体流动、泥石流、垃圾堆填体失稳流滑等岩土材料流滑灾害频发的现状，以及现有计算方法无法模拟大变形流动破坏全过程的局限性，提出一种岩土材料流滑灾变的计算模拟方法，基于不可压缩流体本构和下负荷面剑桥本构模型，并结合固液完全耦合理论，建立考虑固液双相耦合的流滑灾变计算模型，实现对岩土材料大变形流动破坏全过程的模拟再现。为了达到上述目的，本发明提出的岩土材料流滑灾变的计算模拟方法，具体步骤如下
(I)，粒子信息输入及初始化
首先确定问题域，并将问题域离散成间距相等的粒子，将所有粒子的坐标、速度和加速度等信息输入程序并初始化，然后输入计算问题域的标识、模型参数以及所有粒子的初始密度、外力等信息。(2)，邻近粒子搜索
确定粒子的影响半径，采用Verlet neighbor list方法搜索所有影响范围内的所有邻近粒子。首先以一个粒子作为计算粒子，然后根据粒子编号按顺序确定周围粒子，并分别计算周围粒子与所述计算粒子之间的间距，并进行判断，若粒子间距小于光滑长度，则确定为邻近粒子，记录计算列表，参与密度、应力应变等物理量的计算。如此循环往复，直至确定所有的邻近粒子。粒子搜索流程如图2所示。(3)，计算核函数及其微分形式
选用三次B样条函数作为核函数，以2倍粒子间距为光滑长度，根据粒子编号，按顺序计算所有粒子对应的核函数以及核函数的微分形式。核函数表达式如下，其中R为粒子间距，A为光滑长度，a =15/7nh2a
权利要求
1.一种岩土材料流滑灾变的计算模拟方法，其特征在于具体步骤如下 (1)，粒子信息输入及初始化 首先确定问题域，并将问题域离散成间距相等的粒子，将所有粒子的坐标、速度和加速度信息输入程序并初始化，然后输入计算问题域的标识、模型参数以及所有粒子的初始密度、外力信息； (2)，邻近粒子搜索 确定粒子的影响半径，采用Verlet neighbor list方法搜索所有影响范围内的所有邻近粒子；首先以一个粒子作为计算粒子，然后根据粒子编号按顺序确定周围粒子，并分别计算周围粒子与所述计算粒子之间的间距，并进行判断，若粒子间距小于光滑长度，则确定为邻近粒子，记录计算列表，参与密度、应力应变物理量的计算；如此循环往复，直至确定所有的邻近粒子； (3)，计算核函数及其微分形式 选用三次B样条函数作为核函数，以2倍粒子间距为光滑长度，根据粒子编号，按顺序计算所有粒子对应的核函数以及核函数的微分形式；核函数表达式如下，其中W为粒子间距，A为光滑长度，a =15/7 Jih2 ；
全文摘要
本发明涉及一种岩土材料流滑灾变的计算模拟方法，属于计算岩土力学、地质灾害防治和地质环境保护技术领域。针对现有基于网格的计算方法都只适用于小变形分析的局限性，本发明公开了一种基于光滑粒子流体动力学(简称SPH)的岩土材料流滑灾变全过程的计算模拟方法。该方法以Navier-stokes方程为控制方程，采用下负荷面剑桥模型和不可压缩流体本构模型分别描述岩土材料的固液两相，并引入多孔介质理论和达西定律计算固液耦合作用力，建立了考虑水-土完全耦合的岩土材料流滑灾变的计算模型。本发明能够有效再现岩土材料大变形流动破坏的全过程，捕捉其流态化特征，从而为流滑灾害高风险地区的工程设计、工程建设和防灾减灾等提供有力的科学依据，同时有力地推动了计算岩土力学在实际工程建设和地质灾害防治等技术领域的应用。
文档编号G06F17/50GK102819650SQ201210292028
公开日2012年12月12日 申请日期2012年8月16日 优先权日2012年8月16日
发明者黄雨, 戴自立, 张卫杰, 成花林 申请人:同济大学