无粘性土各向异性力学特性微观机理数值模拟方法
【专利摘要】本发明采用颗粒离散元单元方法，开发水平加荷程序和微观力学参数监测程序，给出一种从微观角度研究土体各向异性力学特性的思路和方法：（1）根据室内试验获取无粘性土的应力应变曲线和相关力学参数；（2）根据（1）中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原生各向异性进行模拟；（3）根据（1）中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应力各向异性进行模拟；（4）对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区别与联系，并找出土体两种各向异性产生的微观机理。本发明凭借颗粒流离散元法的在在散粒材料微观结构方面研究的优点，为无粘性土各向异性力学特性的数值模拟提供了新的研究思路。
【专利说明】无粘性土各向异性力学特性微观机理数值模拟方法
【技术领域】
[0001] 本发明属于岩土力学研宄领域，是凭借颗粒流离散元方法在微观结构和微观力学 机制研宄方面
[0002] 的优势，来探讨无粘性土这种散粒材料宏观力学特性各向异性的微观机理的数值 模拟方法。
【背景技术】
[0003] -方面，在工程实践中，土体的应力路径是复杂多变的，由于加荷方式引起的各向 异性的工程问题大量存在，如深、大型基坑的开挖和支护导致的应力状态的改变；真空预压 方法抽、卸真空所引起的小主应力的改变；高土石坝蓄水导致的内部单元应力状态的变化 等，许多工程实例也证明不考虑后期加荷方式引起的各向异性的计算方法是不完善的；另 一方面，现有室内试验手段有限，虽然室内真三轴试验也能实现从不同方向加荷，但是考虑 到室内试验的成本较高，可重复性较难，更重要的是室内真三轴试验，无法在得到土体宏观 各向异性力学特性的同时，给出无粘性土这种散粒材料自身微观结构和微观力学机制的变 化规律，无法找出宏观与微观之间的关系。
[0004] 虽然目前已有不少专家和学者通过颗粒流离散单元法对土体的力学特性进行了 模拟和研宄，但是很少集中在土体应力各向异性的研宄方面上，而且现有的少量关于土体 各向异性的研宄，大多是按与试样沉积面的不同角度（如0度代表竖直，90度代表水平）切 取试样，导致组成切取后试样的颗粒单元发生了变化，即用来从不同方向上加载的试样的 组成成份不完全相同，而且这些研宄也主要是针对原生各向异性的模拟。
【发明内容】
[0005] 鉴于目前对无粘性土各向异性力学特性方面研宄成果的缺乏，以及无法从微观角 度研宄土体各宏观各向异性力学特性的情况，本发明采用颗粒离散元单元方法，开发水平 加荷程序和微观力学参数监测程序，给出一种从微观角度研宄土体各向异性力学特性的思 路和方法，在得出无粘性土各向异性力学特性表现规律的基础上，研宄其形成的微观机理。
[0006] 如图1-图5所示，一种无粘性土应力各向异性的数值模拟方法，流程图如图1，其 具体步骤如下： (1)根据室内试验获取无粘性土的应力应变曲线和相关力学参数；
[0007] (1. 1)收集室内试验曲线，得到土体的应力应变曲线和泊松比曲线，以及相关力学 参数；
[0008] (1. 2)建立颗粒流离散元数值模型，进行参数标定，根据室内试验获取试样的细观 力学参数（如：颗粒大小，摩擦系数，法向刚度，法切向刚度比值等）使数值模型的宏观应力 应变力学特性与室内试验的结果在尽可能的吻合，以保证数值模拟的结果能够对实际工程 具有一定的指导意义；
[0009] (2)根据⑴中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体 的原生各向异性进行模拟；
[0010] (2. 1)生成四面墙体，然后在墙体所围区域内以一较大孔隙率生成颗粒集合体，给 颗粒单元定义微观参数，然后给颗粒单元施加重力，使颗粒在重力作用下发生沉积，沉积到 一定步数后停止，撤去重力方向上的墙体，在新的位置生成墙体，使生成后的试样为正方形 试样，试样生成过程如图1 ;
[0011] (2. 2)为了模拟重力沉积趋势对力学特性的影响，按（2. 1)中方法并给定三种不 同的重力沉积步数生成三组试样，其中每组试样由两个完全相同的试样组成，双轴试验装 置如图1，试样生成过程如图2;
[0012] (2. 3)按重力沉积法生成试样后，需要对试样施加一个较小的各向同性初始应力， 以保证试样与新生成的墙体间接触均匀，然后对以上三组试样进行加载，其中，每组试样中 的两个试样，一个平行重力沉积方向加载，一个垂直重力沉积方向加载，加载方式如图3 ;
[0013] (2.4)得到每组试样的应力应变曲线，并取当加载方向应变为1%时的应力应变 值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值，初始弹模和泊松比定义为：
【权利要求】
1. 一种无粘性土各向异性力学特性微观机理数值模拟方法，其具体步骤如下： (1) 根据室内试验获取无粘性土的应力应变曲线和相关力学参数； (1. 1)收集室内试验曲线，得到土体的应力应变曲线和泊松比曲线，以及相关力学参 数； (1. 2)建立颗粒流离散元数值模型，进行参数标定，根据室内试验获取试样的细观力学 参数，使数值模型的宏观应力应变力学特性与室内试验的结果在尽可能的吻合； (2) 根据（1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原 生各向异性进行模拟； (2. 1)生成四面墙体，然后在墙体所围区域内以2倍于试样最终孔隙率的孔隙率生成 颗粒集合体，给颗粒单元定义微观参数，然后给颗粒单元施加重力，使颗粒在重力作用下发 生沉积，直到试样达到最终孔隙率后，撤去重力方向上的墙体，在新的位置生成墙体，使生 成后的试样为正方形试样； (2.2) 为了模拟重力沉积趋势对力学特性的影响，按（2.1)中方法并给定三种不同的 重力沉积步数生成三组试样，其中每组试样由两个完全相同的试样组成； (2.3) 按重力沉积法生成试样后，需要对试样施加一个较小的各向同性初始应力，这 个初始应力既要能保证试样与新生成的墙体间接触均匀，又不能对试样内部结构产生较大 影响.然后对以上三组试样进行加载，其中，每组试样中的两个试样，一个平行重力沉积方 向加载，一个垂直重力沉积方向加载； (2.4) 得到每组试样的应力应变曲线，并取当加载方向应变为1 %时的应力应变值，计 算初始弹模、泊松比和应力应变曲线上的的应力峰值，初始弹模和泊松比定义为：
Etl:初始弹模；△ ε :初始应变增量（当总应变增量小于〇. 〇1时，施加应力增量产生 的应变增量）；△ σ :应力增量；V(l:泊松比；△ ε ν:体积应变增量；△ ε :加载方向应变增 量。 (2. 5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异 性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律； (2. 6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗 粒间接触力大小和分布的变化以及颗粒间平均配位数的变化； (2. 7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观 参数变化规律进行对比分析，得出结论； (2. 8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分 析，得出土体宏观各向异性力学特性的微观机理； ⑶根据⑴中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应 力各向异性进行模拟； (3. 1)生成四面墙体，然后在墙体所围区域内以2倍于试样最终孔隙率的孔隙率生成 颗粒集合体生成颗粒集合体，给颗粒单元定义微观参数，让颗粒进行循环，使颗粒集合体达 到平衡状态； (3. 2)，按（3. 1)中方法生成3组试样，每组试样包含两个，然后分别从水平和竖直方向 上对每组试样施加不等向初始应力，为了研宄土体所受初始应力状态对其力学特性的影响 及变化规律，每组试样所受的初始应力状态不能相同，3组试样在水平和竖直方向上所受初 始应力的差异逐渐增大，但是每组试样中的两个试样所受初始应力保证完全相同； (3. 3)对以上三组试样进行加载，其中，每组试样中的两个试样，一个沿大主应力方向 加载，一个沿小主应力方向加载； (3.4)得到每组试样的应力应变曲线，并取当加载方向应变为1 %时的应力应变值，计 算初始弹模、泊松比和应力应变曲线上的的应力峰值； (3. 5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异 性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律； (3. 6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗 粒间接触力大小和分布的变化以及颗粒间平均配位数的变化； (3. 7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观 参数变化规律进行对比分析，得出结论； (3. 8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分 析，得出土体宏观应力各向异性力学特性的微观机理； (4)对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区 别与联系，并找出土体两种各向异性产生的微观机理。
【文档编号】G06F17/50GK104517006SQ201410599055
【公开日】2015年4月15日 申请日期:2014年10月30日 优先权日:2014年10月30日
【发明者】张坤勇, 罗兴军, 钟思成, 王耀 申请人:河海大学