CN117235984B - 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 - Google Patents
一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 Download PDFInfo
- Publication number
- CN117235984B CN117235984B CN202311078657.5A CN202311078657A CN117235984B CN 117235984 B CN117235984 B CN 117235984B CN 202311078657 A CN202311078657 A CN 202311078657A CN 117235984 B CN117235984 B CN 117235984B
- Authority
- CN
- China
- Prior art keywords
- sample
- particles
- volume
- membrane
- boundary
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012360 testing method Methods 0.000 title claims abstract description 22
- 238000011439 discrete element method Methods 0.000 title claims abstract description 16
- 239000012528 membrane Substances 0.000 claims abstract description 56
- 239000002245 particle Substances 0.000 claims abstract description 53
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 230000001788 irregular Effects 0.000 claims abstract description 14
- 239000000463 material Substances 0.000 claims abstract description 12
- 239000012798 spherical particle Substances 0.000 claims abstract description 12
- 238000010008 shearing Methods 0.000 claims abstract description 8
- 239000002689 soil Substances 0.000 claims abstract description 8
- 238000007596 consolidation process Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 5
- 238000005096 rolling process Methods 0.000 claims description 4
- 238000003825 pressing Methods 0.000 claims description 3
- 230000000452 restraining effect Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 description 4
- 239000013590 bulk material Substances 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，包括如下步骤：生成膜边界；在膜边界内生成互不接触的球形颗粒或不规则形状颗粒，直到空间内的整体颗粒达到试样预设的孔隙率；根据岩土体离散材料属性，针对试样中颗粒与颗粒之间、颗粒与立方体内壁之间的各个接触位置，分别设置接触模型；在试样的固结和剪切阶段，围压作用在膜边界上的力由围压与膜表面的乘积计算；将试样的体积分为三部分，即试样的上表面或下表面与中心点连接而形成的两个圆锥体，以及圆柱体的其余部分；在每个时间步长重复上两个步骤即可时时记录膜边界的不规则体积变形。本发明可以精细计算试样在三轴剪切过程中任意不规则体积形变。
Description
技术领域
本发明涉及岩土数值模拟技术领域，尤其是一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法。
背景技术
随着计算机技术的普及，通过数值模拟方法研究具有离散特征的岩土体材料物理力学特性受到越来越多的科研工作者的关注。相较于实验室试验周期长，耗时久，人员多，开销大等弊端，数值模拟方法可以大大降低人工成本和试验周期。
离散元是一种常用于岩土领域的数值模拟方法，它通过考虑颗粒间的接触来研究岩土散体材料的物理力学特征，目前在岩土力学方面取得了越来越多的应用。特别是离散元法模拟三轴试验方面，取得了重要的进展。常见的通过离散元法模拟三轴试验的方法包括圆柱边界法，REV局部法，以及膜边界法等。这三种方法各有利弊。圆柱边界法将试验变形近似为圆柱形，计算简单，但难以模拟试样的不规则变形；REV局部法通过试样内部特征区域的局部行为表征试样的整体行为，可以大大减少计算量，但是局部和整体毕竟有差距；膜边界法可以显示剪切过程中样品的不规则变形，与其他边界条件相比，膜边界条件最能体现试样的实际变形情况，但仍有一个关键问题没有得到根本解决，即试样的不规则体积变形的精确计算。由于膜边界的优势在于反应试样在模拟过程中的体积变形，因此变形是不规则的，难以精确跟踪过程中的变化，见图1。目前，常用的计算膜边界体积应变的方法仍为近似法，代表性的膜边界体积应变计算方法包括第一种根据试样半径的变化来计算体积变化。通过确定边界粒子和样本中心线之间的水平距离，获得平均半径，再通过获取试样当前的高度，计算圆柱体的样本体积，该方法本质仍在体积计算时把试样不规则体积变形近似为圆柱形。第二种方法通过平均孔隙比和初始体积来评估体积变化，但是平均孔隙比是通过试样局部REV来获得的，无法如实地代表试样的整体情况。因此，这两种方法尤其是在试样在变形较大的情况下的误差不能忽略。
发明内容
本发明所要解决的技术问题在于，提供一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，可以精细计算试样在三轴剪切过程中任意不规则体积形变。
为解决上述技术问题，本发明提供一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，包括如下步骤：
步骤1、生成膜边界；
步骤2、在膜边界内生成互不接触的球形颗粒或不规则形状颗粒，直到空间内的整体颗粒达到试样预设的孔隙率；
步骤3、根据岩土体离散材料属性，针对试样中颗粒与颗粒之间、颗粒与立方体内壁之间的各个接触位置，分别设置接触模型；
步骤4、在试样的固结和剪切阶段，围压作用在膜边界上的力由围压与膜表面的乘积计算；
步骤5、将试样的体积分为三部分，即试样的上表面或下表面与中心点连接而形成的两个圆锥体，以及圆柱体的其余部分；
步骤6、在每个时间步长重复步骤4和步骤5即可时时记录膜边界的不规则体积变形。
优选的，步骤1中，膜边界条件是用规则排形的球颗粒组成，即膜颗粒，膜颗粒采用球颗粒。
优选的，步骤2中，采用颗粒膨胀法生成互不接触的球形颗粒或不规则形状颗粒。
优选的，步骤2中，针对立方体堆积体，采用夯实试样或振动试样的方式处理，制备满足预设孔隙率的试样。
优选的，步骤2中，针对立方体堆积体的顶部，采用上压板针对立方体堆积体的顶部进行击打，实现夯实试样方式处理，制备满足预设孔隙率的试样。
优选的，步骤2中，针对立方体堆积体的底部位置设置下底板，控制下底板以默认100Hz的频率振动，通过振动试样方式处理，制备满足预设孔隙率的试样。
优选的，步骤3中，接触模型包括线性接触模型，Hertz模型，平行粘接模型，抗滚动模型。根据不同的材料属性进行选择，进行简单估算可选用线性接触模型或Hertz模型，若有粘结的散体材料可选用平行粘结模型，若颗粒不规则形状较为明显，则选用抗滚动模型。
优选的，步骤4中，每三个相邻的膜颗粒的中心形成一个三角形，膜的表面积对所有三个相邻膜颗粒的中心形成的小三角形的求和；然后，通过膜颗粒和材料之间的接触，将膜上的约束力传递到试样上，在样品被大幅剪切后，相邻膜颗粒组成的三角形并不总是等边的。
优选的，计算三角形面积具体包括如下步骤：
步骤41、假设三角形三点坐标分别为P1(x1,y1,z1,),P2(x2,y2,z2,),P3(x3,y3,z3,),则三角形形心P0(x0,y0,z0,)计算公式为：
步骤42、三角形面积计算公式为：
F＝p×A。
优选的，步骤5中，第三部分的体积是所有三维单纯形之和，三维单纯形为三个相邻膜颗粒的三角形和样品中心形成的体积，计算第三部分的体积具体包括如下步骤：
步骤51、圆锥体体积
其中A为上下底面积，H为上下底面分别到质心的距离；
步骤52、区域三内单个三维单纯形体积：
假设某三维单纯形四点坐标
P1(x1,y1,z1,),P2(x2,y2,z2,),P3(x3,y3,z3,),M(x4,y4,z4,),则该三维单纯形的体积计算公式为：
其中J为雅可比行列式：
所有三维单纯形体积为
Vrest＝∫Vsimple
步骤53、总体积计算公式：
Vtotat＝Vcone1+Vcone2+Vrest。
本发明的有益效果为：本发明可以时时精确计算任意体积变形的试样体积应变，解决了基于离散元法的三轴试验膜边界试样模拟剪切形为的痛点问题，通过精确计算试样的体积形变有利于完善数值模型的精确程度，实现虚拟实验对室内实验的复制，模拟材料宏观属性，获取模拟材料的膨胀角、特征角等参数。
附图说明
图1为本发明膜边界不规则形变示意图。
图2为本发明的膜边界及膜颗粒分布示意图。
图3为本发明的体积区域划分及三维单纯形示意图。
具体实施方式
一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，包括如下步骤：
步骤A：生成膜边界。如图1所示，膜边界条件通常是用规则排形的球颗粒组成，简称膜颗粒。在初始阶段是圆柱形，与试样一致。散体材料填充在膜边界内。随着剪切的过程，膜边界上每个膜颗粒受力发生有限自由运动，由于膜颗粒不是统一行动，因此从整体看，膜变形即体现了试样的不规则形状。膜边界条件由相同的重叠膜颗粒组成，高度分布189个圆，每个圆201个颗粒，总共37989个膜颗粒，如图2(a)所示。膜颗粒的数量越大，越能精确表征试样的不规则变形细节，但是计算时间也会相应增加。因此，膜粒径的选择既考虑了计算效率，又考虑了模型精度。膜粒子的接触模型是线性接触键模型。然后进入步骤B；
步骤B：在膜边界内生成互不接触的球形颗粒或不规则形状颗粒,直到空间内的整体颗粒达到试样预设的孔隙率，然后进入步骤C；
采用颗粒膨胀法生成互不接触的球形颗粒或不规则形状颗粒。针对立方体堆积体，采用夯实试样或振动试样的方式处理，制备满足预设孔隙率的试样。针对立方体堆积体的顶部，采用上压板针对立方体堆积体的顶部进行击打，实现夯实试样方式处理，制备满足预设孔隙率的试样。针对立方体堆积体的底部位置设置下底板，控制下底板以默认100Hz的频率振动，通过振动试样方式处理，制备满足预设孔隙率的试样。接触模型包括线性接触模型，Hertz模型，平行粘接模型，抗滚动模型。
步骤C：根据岩土体离散材料属性，针对试样中颗粒与颗粒之间、颗粒与立方体内壁之间的各个接触位置，分别设置接触模型，然后进入步骤D；
步骤D：在试样的固结和剪切阶段，围压作用在膜边界上的力由围压与膜表面的乘积计算。每三个相邻的膜颗粒的中心形成一个三角形，如图2(b)所示。膜的表面积对所有三个相邻膜颗粒的中心形成的小三角形的求和。然后，通过膜颗粒和材料之间的接触，可以将膜上的约束力传递到试样上。在样品被大幅剪切后，相邻膜颗粒组成的三角形并不总是等边的，如图2(c)所示。
步骤D-1：
假设三角形三点坐标分别为P1(x1,y1,z1,),P2(x2,y2,z2,),P3(x3,y3,z3,),则三角形形心P0(x0,y0,z0,)计算公式为：
步骤D-2：
三角形面积计算公式为：
步骤E：将试样的体积分为三部分，即试样的上表面或下表面与中心点连接而形成的两个圆锥体，以及圆柱体的其余部分，如图3所示。其中，第三部分的体积是所有三维单纯形(三个相邻膜颗粒的三角形和样品中心形成的体积)之和。
步骤E-1：圆锥体体积
其中A为上下底面积，H为上下底面分别到质心的距离。
步骤E-2：区域三内单个三维单纯形体积：
假设某三维单纯形四点坐标
P1(x1，y1，z1，)，P2(x2，y2，z2，)，P3(x3，y3，z3，)，M(x4，y4，z4)，则该三维单纯形的体积计算公式为：
其中J为雅可比行列式：
所有三维单纯形体积为
Vrest＝∫Vsimple
步骤E-3：
总体积计算公式：
Vtotal＝Vcone1+Vcone2+Vrest
步骤F：在每个时间步长重复步骤D和步骤E即可时时记录膜边界的不规则体积变形。
Claims (9)
1.一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，包括如下步骤：
步骤1、生成膜边界；
步骤2、在膜边界内生成互不接触的球形颗粒或不规则形状颗粒，直到空间内的整体颗粒达到试样预设的孔隙率；
步骤3、根据岩土体离散材料属性，针对试样中颗粒与颗粒之间、颗粒与立方体内壁之间的各个接触位置，分别设置接触模型；
步骤4、在试样的固结和剪切阶段，围压作用在膜边界上的力由围压与膜表面的乘积计算；
步骤5、将试样的体积分为三部分，即试样的上表面或下表面与中心点连接而形成的两个圆锥体，以及圆柱体的其余部分；圆柱体除了两个圆锥体的其余部分的体积是所有三维单纯形之和，三维单纯形为三个相邻膜颗粒的三角形和样品中心形成的体积，计算第三部分的体积具体包括如下步骤：
步骤51、圆锥体体积
其中A为上下底面积，H为上下底面分别到质心的距离；
步骤52、区域三内单个三维单纯形体积：
假设某三维单纯形四点坐标
P1(x1,y1,z1,),P2(x2,y2,z2,),P3(x3,y3,z3,),M(x4,y4,z4,),则该三维单纯形的体积计算公式为：
其中J为雅可比行列式：
所有三维单纯形体积为V/es0＝∫Vsimple
步骤53、总体积计算公式：
V0o01l＝Vcone1+Vcone2+V/es0；
步骤6、在每个时间步长重复步骤4和步骤5即可时时记录膜边界的不规则体积变形。
2.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤1中，膜边界条件是用规则排形的球颗粒组成，即膜颗粒，膜颗粒采用球颗粒。
3.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤2中，采用颗粒膨胀法生成互不接触的球形颗粒或不规则形状颗粒。
4.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤2中，针对立方体堆积体，采用夯实试样或振动试样的方式处理，制备满足预设孔隙率的试样。
5.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤2中，针对立方体堆积体的顶部，采用上压板针对立方体堆积体的顶部进行击打，实现夯实试样方式处理，制备满足预设孔隙率的试样。
6.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤2中，针对立方体堆积体的底部位置设置下底板，控制下底板以默认100Hz的频率振动，通过振动试样方式处理，制备满足预设孔隙率的试样。
7.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤3中，接触模型包括线性接触模型，Hertz模型，平行粘接模型，抗滚动模型。
8.如权利要求1所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，步骤4中，每三个相邻的膜颗粒的中心形成一个三角形，膜的表面积对所有三个相邻膜颗粒的中心形成的小三角形的求和；然后，通过膜颗粒和材料之间的接触，将膜上的约束力传递到试样上，在样品被大幅剪切后，相邻膜颗粒组成的三角形并不总是等边的。
9.如权利要求8所述的用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法，其特征在于，计算三角形面积具体包括如下步骤：
步骤41、假设三角形三点坐标分别为P1(x1,t1,z1,),P2(x2,y2,z2,),P3(x3,y3,z3,),则三角形形心P-(x-,y-,z-,)计算公式为：
步骤42、三角形面积计算公式为：
F＝p×∑A。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311078657.5A CN117235984B (zh) | 2023-08-25 | 2023-08-25 | 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 |
| LU507010A LU507010B1 (en) | 2023-08-25 | 2023-12-05 | Method for accurately calculating volumetric strain of flexible membrane in triaxial test simulated by discrete element method |
| PCT/CN2023/136474 WO2024120391A1 (zh) | 2023-08-25 | 2023-12-05 | 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311078657.5A CN117235984B (zh) | 2023-08-25 | 2023-08-25 | 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117235984A CN117235984A (zh) | 2023-12-15 |
| CN117235984B true CN117235984B (zh) | 2024-06-04 |
Family
ID=89087042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311078657.5A Active CN117235984B (zh) | 2023-08-25 | 2023-08-25 | 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN117235984B (zh) |
| LU (1) | LU507010B1 (zh) |
| WO (1) | WO2024120391A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117521475B (zh) * | 2024-01-08 | 2024-06-07 | 佛山市交通科技有限公司 | 沥青混合料三轴试验柔性边界的离散元模拟方法及设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107330227A (zh) * | 2017-07-31 | 2017-11-07 | 南京大学 | 考虑膜效应的岩土体三轴试验的离散元建模及数值模拟方法 |
| CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
| WO2019239176A1 (en) * | 2018-06-14 | 2019-12-19 | Total Sa | Method for determining the geometry of an area of a reservoir |
| CN112630040A (zh) * | 2021-01-06 | 2021-04-09 | 合肥工业大学 | 三向刚性加载条件下的单颗粒强度测试装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8560286B2 (en) * | 2011-03-31 | 2013-10-15 | Dem Solutions Limited | Method and apparatus for discrete element modeling involving a bulk material |
| CN109960856B (zh) * | 2019-02-26 | 2023-04-18 | 东南大学 | 一种基于圆柱体伺服条件的假三轴剪切试验离散元数值模拟方法 |
-
2023
- 2023-08-25 CN CN202311078657.5A patent/CN117235984B/zh active Active
- 2023-12-05 WO PCT/CN2023/136474 patent/WO2024120391A1/zh unknown
- 2023-12-05 LU LU507010A patent/LU507010B1/en active IP Right Grant
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107330227A (zh) * | 2017-07-31 | 2017-11-07 | 南京大学 | 考虑膜效应的岩土体三轴试验的离散元建模及数值模拟方法 |
| WO2019239176A1 (en) * | 2018-06-14 | 2019-12-19 | Total Sa | Method for determining the geometry of an area of a reservoir |
| CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
| CN112630040A (zh) * | 2021-01-06 | 2021-04-09 | 合肥工业大学 | 三向刚性加载条件下的单颗粒强度测试装置 |
Non-Patent Citations (3)
| Title |
|---|
| Discrete element modelling of flexible membrane boundaries for triaxial tests;Tongming Qu 等;Computers and Geotechnics;20191231;全文 * |
| 军用车辆松软路面多工况牵引特性数值仿真与试验验证研究;李桥;中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑;20220215;全文 * |
| 土石混合体宏细观力学特性和变形破坏机制的三维离散元精细模拟;金磊 等;岩石力学与工程学报;20180315(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| LU507010B1 (en) | 2024-07-03 |
| WO2024120391A1 (zh) | 2024-06-13 |
| CN117235984A (zh) | 2023-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117235984B (zh) | 一种用于精确计算离散元法模拟三轴试验柔性膜体积应变的方法 | |
| Li et al. | A 3D grain-based model for simulating the micromechanical behavior of salt rock | |
| Stahl et al. | Discrete element simulation of ballast and gravel under special consideration of grain-shape, grain-size and relative density | |
| Nezami et al. | Simulation of front end loader bucket–soil interaction using discrete element method | |
| Yu et al. | Granular discrete element simulation of the evolution characteristics of the shear band in soil–rock mixture based on particle rotation analysis | |
| O’Sullivan et al. | Examination of the response of regularly packed specimens of spherical particles using physical tests and discrete element simulations | |
| Manso et al. | Crushing and oedometer compression of rockfill using DEM | |
| CN107169236B (zh) | 一种基于有限元与离散元耦合的虚拟三轴试验仿真方法 | |
| CN110761342A (zh) | 爆破荷载下基坑施工过程稳定性研究模型试验系统及方法 | |
| CN115659449B (zh) | 基于轻质水泥土桩复合地基的计算方法、系统及应用 | |
| Clément et al. | Sinking during earthquakes: Critical acceleration criteria control drained soil liquefaction | |
| CN115050431B (zh) | 一种水泥稳定再生集料的三维细观结构的建模分析方法 | |
| CN109541186A (zh) | 一种基于形状参数的粗骨料密实度计算方法 | |
| CN114613458A (zh) | 一种基于常规态型近场动力学的颗粒破碎模拟方法及系统 | |
| CN116579195A (zh) | 一种三维结构面岩体多尺度非连续模型建立和仿真方法 | |
| CN111289346A (zh) | 含断层破碎带隧道围岩变形破坏的三维模型试验方法 | |
| CN114218651A (zh) | 砂性土地基平板载荷试验的颗粒离散元仿真方法和设备 | |
| CN114626215A (zh) | 一种基于离散元与合成岩体技术的岩质边坡建模方法 | |
| Orosz et al. | Coupling Finite And Discrete Element Methods Using An Open Source And A Commercial Software. | |
| CN108629089A (zh) | 一种基于三维连续-离散单元耦合的强夯加固地基模拟方法 | |
| CN218765897U (zh) | 一种模拟不同产状断裂带错动桩基承载特性的试验装置 | |
| Ng | Behavior of ellipsoids of two sizes | |
| Ullidtz | Distinct element method for study of failure in cohesive particulate media | |
| Qin et al. | Influence of fabric anisotropy on seismic response of circular foundation | |
| Yohannes et al. | Mechanistic modeling of unbound granular materials |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |