CN104517006B - 无粘性土各向异性力学特性微观机理数值模拟方法 - Google Patents

无粘性土各向异性力学特性微观机理数值模拟方法 Download PDF

Info

Publication number
CN104517006B
CN104517006B CN201410599055.9A CN201410599055A CN104517006B CN 104517006 B CN104517006 B CN 104517006B CN 201410599055 A CN201410599055 A CN 201410599055A CN 104517006 B CN104517006 B CN 104517006B
Authority
CN
China
Prior art keywords
sample
stress
anisotropy
samples
parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410599055.9A
Other languages
English (en)
Other versions
CN104517006A (zh
Inventor
张坤勇
罗兴军
钟思成
王耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201410599055.9A priority Critical patent/CN104517006B/zh
Publication of CN104517006A publication Critical patent/CN104517006A/zh
Application granted granted Critical
Publication of CN104517006B publication Critical patent/CN104517006B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明采用颗粒离散元单元方法，开发水平加荷程序和微观力学参数监测程序，给出一种从微观角度研究土体各向异性力学特性的思路和方法：（1）根据室内试验获取无粘性土的应力应变曲线和相关力学参数；（2）根据（1）中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原生各向异性进行模拟；（3）根据（1）中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应力各向异性进行模拟；（4）对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区别与联系，并找出土体两种各向异性产生的微观机理。本发明凭借颗粒流离散元法的在散粒材料微观结构方面研究的优点，为无粘性土各向异性力学特性的数值模拟提供了新的研究思路。

Description

无粘性土各向异性力学特性微观机理数值模拟方法
技术领域
本发明属于岩土力学研究领域，是凭借颗粒流离散元方法在微观结构和微观力学机制研究方面
的优势，来探讨无粘性土这种散粒材料宏观力学特性各向异性的微观机理的数值模拟方法。
背景技术
一方面，在工程实践中，土体的应力路径是复杂多变的，由于加荷方式引起的各向异性的工程问题大量存在，如深、大型基坑的开挖和支护导致的应力状态的改变；真空预压方法抽、卸真空所引起的小主应力的改变；高土石坝蓄水导致的内部单元应力状态的变化等，许多工程实例也证明不考虑后期加荷方式引起的各向异性的计算方法是不完善的；另一方面，现有室内试验手段有限，虽然室内真三轴试验也能实现从不同方向加荷，但是考虑到室内试验的成本较高，可重复性较难，更重要的是室内真三轴试验，无法在得到土体宏观各向异性力学特性的同时，给出无粘性土这种散粒材料自身微观结构和微观力学机制的变化规律，无法找出宏观与微观之间的关系。
虽然目前已有不少专家和学者通过颗粒流离散单元法对土体的力学特性进行了模拟和研究，但是很少集中在土体应力各向异性的研究方面上，而且现有的少量关于土体各向异性的研究，大多是按与试样沉积面的不同角度(如0度代表竖直，90度代表水平)切取试样，导致组成切取后试样的颗粒单元发生了变化，即用来从不同方向上加载的试样的组成成份不完全相同，而且这些研究也主要是针对原生各向异性的模拟。
发明内容
鉴于目前对无粘性土各向异性力学特性方面研究成果的缺乏，以及无法从微观角度研究土体各宏观各向异性力学特性的情况，本发明采用颗粒离散元单元方法，开发水平加荷程序和微观力学参数监测程序，给出一种从微观角度研究土体各向异性力学特性的思路和方法，在得出无粘性土各向异性力学特性表现规律的基础上，研究其形成的微观机理。
如图1-图5所示，一种无粘性土应力各向异性的数值模拟方法，流程图如图1，其具体步骤如下：
(1)根据室内试验获取无粘性土的应力应变曲线和相关力学参数；
(1.1)收集室内试验曲线，得到土体的应力应变曲线和泊松比曲线，以及相关力学参数；
(1.2)建立颗粒流离散元数值模型，进行参数标定，根据室内试验获取试样的细观力学参数(如：颗粒大小，摩擦系数，法向刚度，法切向刚度比值等)使数值模型的宏观应力应变力学特性与室内试验的结果在尽可能的吻合，以保证数值模拟的结果能够对实际工程具有一定的指导意义；
(2)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原生各向异性进行模拟；
(2.1)生成四面墙体，然后在墙体所围区域内以一较大孔隙率生成颗粒集合体，给颗粒单元定义微观参数，然后给颗粒单元施加重力，使颗粒在重力作用下发生沉积，沉积到一定步数后停止，撤去重力方向上的墙体，在新的位置生成墙体，使生成后的试样为正方形试样，试样生成过程如图1；
(2.2)为了模拟重力沉积趋势对力学特性的影响，按(2.1)中方法并给定三种不同的重力沉积步数生成三组试样，其中每组试样由两个完全相同的试样组成，双轴试验装置如图1，试样生成过程如图2；
(2.3)按重力沉积法生成试样后，需要对试样施加一个较小的各向同性初始应力，以保证试样与新生成的墙体间接触均匀，然后对以上三组试样进行加载，其中，每组试样中的两个试样，一个平行重力沉积方向加载，一个垂直重力沉积方向加载，加载方式如图3；
(2.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值，初始弹模和泊松比定义为：
E0:初始弹模；Δε:初始应变增量(当总应变增量小于0.01时,施加应力增量产生的应变增量)；Δσ:应力增量；v0:泊松比；Δεv:体积应变增量；Δε:加载方向应变增量.
(2.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(2.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化，颗粒长轴定向如图4；
(2.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(2.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观各向异性力学特性的微观机理；
(3)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应力各向异性进行模拟，技术路线如图5；
(3.1)生成四面墙体，然后在墙体所围区域内以一较大孔隙率生成颗粒集合体，给颗粒单元定义微观参数，让颗粒进行循环，使颗粒集合体达到平衡状态；
(3.2)，按(3.1)中方法生成3组试样，每组试样包含两个，然后分别从水平和竖直方向上对每组试样施加不等向初始应力。为了研究土体所受初始应力状态对其力学特性的影响及变化规律，每组试样所受的初始应力状态不能相同，3组试样在水平和竖直方向上所受初始应力的差异逐渐增大，但是每组试样中的两个试样所受初始应力保证完全相同；
(3.3)对以上三组试样进行加载，其中，每组试样中的两个试样，一个沿大主应力方向加载，一个沿小主应力方向加载；
(3.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值；
(3.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(3.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化；
(3.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(3.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观应力各向异性力学特性的微观机理；
(4)对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区别与联系，并找出土体两种各向异性产生的微观机理。
本发明无粘性土各向异性力学特性微观机理研究方法优点是：
在采用有限元模拟土体各向异性特性时是将土体当作一种连续介质，室内试验也无法在试验过程中记录土体微观结构和微观力学机制的演化，无法提示散粒材料各向异性力学特性的微观机理。而本发明采用颗粒离散元的模拟方法克服了这些缺陷，既可以方便地模拟无粘性土原生各向异性和应力各向异性的表现规律，进行比较和分析，又可以提示土体各向异性力学特性在微观结构和微观力学特性上的演化机理。本发明一方面针对无粘性土的散粒特点，采用长条单元模拟土体颗粒，与实际土体颗粒形状和颗粒间的咬合力更加接近。按相同方法生成平面双轴数值试样，两组试样具有相同初始应力状态，保证两组试样在加载前的颗粒组成，受力状态均完全一样，加载过程中，一个从小主应力方向加载，另一个从大主应力方向加载，使土体的受荷方式与工程实践更加接近，而且很方便的比较不同加载方向的力学曲线和力学参数，得出土体各向异性力学特性的表现规律，以及两种各向异性的联系与区别；另一方面，凭借颗粒流离散单元法的优势，不断记录加载过程中，土体内部颗粒长轴定向，位移，颗粒间接触力，平均配位数等微观参数的变化，找出土体宏观各向异性力学特性的微观机理。因此，采用颗粒流离散单元法模拟无粘性土的各向异性力学特性具有一定的优势。本发明凭借颗粒流离散元法的在在散粒材料微观结构方面研究的优点，为无粘性土各向异性力学特性的数值模拟提供了新的研究思路。
此外，根据数值模拟结果，对工程实践中因不考虑土体各向异性而引发的工程事故进行解释，并从施工或设计方面提出改进措施，从而对实际工程具有一定的指导意义。
附图说明
图1为本发明实施过程流程图
图2为长条颗粒双轴试样重力沉积过程图
图3为原生各向异性试样加载过程
图4为颗粒单元长轴定向
图5应力各向异性加载过程
具体实施方式
如图1-图5所示，一种无粘性土应力各向异性的数值模拟方法，流程图如图1，其具体步骤如下：
(1)根据室内试验获取无粘性土的应力应变曲线和相关力学参数。
(1.1)收集室内试验曲线，得到土体的应力应变曲线和泊松比曲线，以及相关力学参数。
(1.2)建立颗粒流离散元数值模型，进行参数标定，根据室内试验获取试样的细观力学参数(如：颗粒大小，摩擦系数，法向刚度，法切向刚度比值等)使数值模型的宏观应力应变力学特性与室内试验的结果在尽可能的吻合，以保证数值模拟的结果能够对实际工程具有一定的指导意义；
(2)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原生各向异性进行模拟；
(2.1)生成四面墙体，然后在墙体所围区域内以一较大孔隙率生成颗粒集合体，给颗粒单元定义微观参数，然后给颗粒单元施加重力，使颗粒在重力作用下发生沉积，沉积到一定步数后停止，撤去重力方向上的墙体，在新的位置生成墙体，使生成后的试样为正方形试样，试样生成过程如图1。
(2.2)为了模拟重力沉积趋势对力学特性的影响，按(2.1)中方法并给定三种不同的重力沉积步数生成三组试样，其中每组试样由两个完全相同的试样组成，双轴试验装置如图1，试样生成过程如图2。
(2.3)按重力沉积法生成试样后，需要对试样施加一个较小的各向同性初始应力，以保证试样与新生成的墙体间接触均匀，然后对以上三组试样进行加载，其中，每组试样中的两个试样，一个平行重力沉积方向加载，一个垂直重力沉积方向加载，加载方式如图3。
(2.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值，初始弹模和泊松比定义为：
(2.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(2.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化。颗粒长轴定向如图4。
(2.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(2.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观各向异性力学特性的微观机理；
(3)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应力各向异性进行模拟，技术路线如图5；
(3.1)生成四面墙体，然后在墙体所围区域内以一较大孔隙率生成颗粒集合体，给颗粒单元定义微观参数，让颗粒进行循环，使颗粒集合体达到平衡状态。
(3.2)，按(3.1)中方法生成3组试样，每组试样包含两个，然后分别从水平和竖直方向上对每组试样施加不等向初始应力。为了研究土体所受初始应力状态对其力学特性的影响及变化规律，每组试样所受的初始应力状态不能相同，3组试样在水平和竖直方向上所受初始应力的差异逐渐增大，但是每组试样中的两个试样所受初始应力保证完全相同。
(3.3)对以上三组试样进行加载，其中，每组试样中的两个试样，一个沿大主应力方向加载，一个沿小主应力方向加载。
(3.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值。
(3.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(3.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化。
(3.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(3.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观应力各向异性力学特性的微观机理；
(4)对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区别与联系，并找出土体两种各向异性产生的微观机理。

Claims (1)

1.一种无粘性土各向异性力学特性微观机理数值模拟方法，其具体步骤如下：
(1)根据室内试验获取无粘性土的应力应变曲线和相关力学参数；
(1.1)收集室内试验曲线，得到土体的应力应变曲线和泊松比曲线，以及相关力学参数；
(1.2)建立颗粒流离散元数值模型，进行参数标定，根据室内试验获取试样的细观力学参数，使数值模型的宏观应力应变力学特性与室内试验的结果在尽可能的吻合；
(2)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的原生各向异性进行模拟；
(2.1)生成四面墙体，然后在墙体所围区域内以2倍于试样最终孔隙率的孔隙率生成颗粒集合体，给颗粒单元定义微观参数，然后给颗粒单元施加重力，使颗粒在重力作用下发生沉积，直到试样达到最终孔隙率后，撤去重力方向上的墙体，在新的位置生成墙体，使生成后的试样为正方形试样；
(2.2)为了模拟重力沉积趋势对力学特性的影响，按(2.1)中方法并给定三种不同的重力沉积步数生成三组试样，其中每组试样由两个完全相同的试样组成；
(2.3)按重力沉积法生成试样后，需要对试样施加一个较小的各向同性初始应力,这个初始应力既要能保证试样与新生成的墙体间接触均匀，又不能对试样内部结构产生较大影响.然后对以上三组试样进行加载，其中，每组试样中的两个试样，一个平行重力沉积方向加载，一个垂直重力沉积方向加载；
(2.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值，初始弹模和泊松比定义为：
E0:初始弹模；Δε:初始应变增量(当总应变增量小于0.01时,施加应力增量产生的应变增量)；Δσ:应力增量；v0:泊松比；Δεv:体积应变增量；Δε:加载方向应变增量;
(2.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(2.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化；
(2.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(2.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观各向异性力学特性的微观机理；
(3)根据(1)中颗粒流双轴试验确定的细观力学参数，建立平面双轴试样，对土体的应力各向异性进行模拟；
(3.1)生成四面墙体，然后在墙体所围区域内以2倍于试样最终孔隙率的孔隙率生成颗粒集合体生成颗粒集合体，给颗粒单元定义微观参数，让颗粒进行循环，使颗粒集合体达到平衡状态；
(3.2)，按(3.1)中方法生成3组试样，每组试样包含两个，然后分别从水平和竖直方向上对每组试样施加不等向初始应力，为了研究土体所受初始应力状态对其力学特性的影响及变化规律，每组试样所受的初始应力状态不能相同，3组试样在水平和竖直方向上所受初始应力的差异逐渐增大，但是每组试样中的两个试样所受初始应力保证完全相同；
(3.3)对以上三组试样进行加载，其中，每组试样中的两个试样，一个沿大主应力方向加载，一个沿小主应力方向加载；
(3.4)得到每组试样的应力应变曲线，并取当加载方向应变为1％时的应力应变值，计算初始弹模、泊松比和应力应变曲线上的的应力峰值；
(3.5)比较每一组试样中的两个试样的力学特性曲线以及力学参数，得出原生各向异性的规律，进一步比较不同重力沉积趋势下，原生各向异性的发展趋势和规律；
(3.6)试样加载过程中，编写程序记录微观参数的变化如颗粒单元长轴定向的变化，颗粒间接触力大小和分布的变化以及颗粒间平均配位数的变化；
(3.7)分析每一组试样中的两个试样的微观参数的变化规律，然后将三组试样的微观参数变化规律进行对比分析，得出结论；
(3.8)将土体宏观力学特性的表现和变化规律与其微观参数的变化规律进行比较分析，得出土体宏观应力各向异性力学特性的微观机理；
(4)对比分析土体原生各向异性和应力各向异性的表现规律，得出两种各向异性的区别与联系，并找出土体两种各向异性产生的微观机理。
CN201410599055.9A 2014-10-30 2014-10-30 无粘性土各向异性力学特性微观机理数值模拟方法 Active CN104517006B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410599055.9A CN104517006B (zh) 2014-10-30 2014-10-30 无粘性土各向异性力学特性微观机理数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410599055.9A CN104517006B (zh) 2014-10-30 2014-10-30 无粘性土各向异性力学特性微观机理数值模拟方法

Publications (2)

Publication Number Publication Date
CN104517006A CN104517006A (zh) 2015-04-15
CN104517006B true CN104517006B (zh) 2018-05-08

Family

ID=52792300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410599055.9A Active CN104517006B (zh) 2014-10-30 2014-10-30 无粘性土各向异性力学特性微观机理数值模拟方法

Country Status (1)

Country Link
CN (1) CN104517006B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106124316B (zh) * 2016-06-28 2019-03-22 河海大学 散粒体岩土材料各向异性的室内真三轴试验方法
CN108152308A (zh) * 2017-11-29 2018-06-12 中国电力科学研究院有限公司 一种土体颗粒间化学胶结力的检测方法
CN108844818A (zh) * 2018-04-28 2018-11-20 辽宁工程技术大学 一种确定岩石在加载任意时刻孔隙率的方法
CN111062085B (zh) * 2019-12-31 2023-04-11 浙江科技学院 减少附加竖向荷载对挡墙水平压力影响的方法
CN110909414B (zh) * 2019-12-31 2023-04-11 浙江科技学院 减少附加竖向荷载对挡墙水平压力影响的随机填颗粒方法
CN112163328B (zh) * 2020-09-18 2022-03-04 武汉大学 一种基于深度学习和数据驱动的岩土颗粒材料本构建模方法
CN113008675B (zh) * 2021-03-01 2024-04-26 浙江科技学院 原状砂胶结失效的检测方法
CN113627061B (zh) * 2021-08-12 2022-07-05 江苏中设集团股份有限公司 一种双排桩基坑支护颗粒流数值模拟方法
CN114969892B (zh) * 2022-04-20 2024-05-17 西安建筑科技大学 碎石粒料颗粒接触方式的界定量化方法、***及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101024949A (zh) * 2007-03-09 2007-08-29 武亚军 振动增压快速固结软地基处理的方法
CN101838985A (zh) * 2010-05-19 2010-09-22 东南大学 用固化剂改良过湿黏性土壤的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101024949A (zh) * 2007-03-09 2007-08-29 武亚军 振动增压快速固结软地基处理的方法
CN101838985A (zh) * 2010-05-19 2010-09-22 东南大学 用固化剂改良过湿黏性土壤的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
单调与循环加荷条件下砂土力学性质细观模拟;史旦达;《中国博士学位论文全文数据库 工程科技II辑（月刊）》;20091015;第2009年卷(第10期);参见正文第2-4章 *
土体两种各向异性的区别与联系;张坤勇等;《岩体力学与工程学报》;20050531;第24卷(第9期);参见正文第1-5节 *

Also Published As

Publication number Publication date
CN104517006A (zh) 2015-04-15

Similar Documents

Publication Publication Date Title
CN104517006B (zh) 无粘性土各向异性力学特性微观机理数值模拟方法
Kim et al. DEM simulation of collapse behaviours of unsaturated granular materials under general stress states
CN103323575B (zh) 一种沉积岩石的重构方法及***
Valore et al. The bearing capacity of footings on sand with a weak layer
Flores-Johnson et al. Discrete element simulation of dynamic behaviour of partially saturated sand
Zhao et al. Effects of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of pipeline buried in elasto-plastic seabed foundations
Zhao et al. Parametric study of the wave-induced residual liquefaction around an embedded pipeline
CN110705087A (zh) 一种含封闭应力的岩石离散元试样构建方法
Seyyedan et al. Numerical simulation of direct shear test on granular materials composed of breakable angular particles: A DEM-XFEM approach
CN113205858A (zh) 模拟颤振环境铰链机构运动分子动力学仿真模型及方法
CN112485106B (zh) 一种控制土体状态参数的物理模型分层制备与试验方法
CN110987749B (zh) 多相复合材料的等效渗透系数细观尺度研究方法
Wensrich et al. Shock waves in granular materials: Discrete and continuum comparisons
Zhou et al. Uplift mechanism for a shallow-buried structure in liquefiable sand subjected to seismic load: centrifuge model test and DEM modeling
Shimomura et al. Bubble growth processes in magma surrounded by an elastic medium
Pan et al. DEM simulation and fractal analysis of particulate fouling on coal-fired utility boilers' heating surfaces
Tafreshi et al. Cyclic loading on foundation to evaluate the coefficient of elastic uniform compression of sand
Ter-Martirosyan et al. Rheological properties of sandy soils
Kadau et al. Collapsing granular suspensions
CN103645120B (zh) 一种置换型合金中影响间隙原子溶解行为的表征方法
Shrestha et al. Experimental and Numerical Analysis on the Rate-Dependent Cyclic Lateral Load on a Single Pile in Sand
Ecemis et al. Feasible packing of granular materials in discrete-element modelling of cone-penetration testing
Ye et al. Investigation on the reliquefaction behaviors of sand using shaking table tests
Humaish et al. Design, manufacturing and testing of small shaking table
Zhang et al. Mesoscopic Parameter Calibration of Granite Particles Based on Neural Network and Discrete Element Triaxial Numerical Simulation Test

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant