CN107644120A - 一种用于岩土体三维离散元快速建模和模拟的通用数值模拟箱 - Google Patents
一种用于岩土体三维离散元快速建模和模拟的通用数值模拟箱 Download PDFInfo
- Publication number
- CN107644120A CN107644120A CN201710700733.XA CN201710700733A CN107644120A CN 107644120 A CN107644120 A CN 107644120A CN 201710700733 A CN201710700733 A CN 201710700733A CN 107644120 A CN107644120 A CN 107644120A
- Authority
- CN
- China
- Prior art keywords
- mrow
- particle
- model
- msub
- displacement
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 50
- 239000002689 soil Substances 0.000 title claims abstract description 24
- 239000011435 rock Substances 0.000 title claims abstract description 23
- 239000002245 particle Substances 0.000 claims abstract description 88
- 238000006073 displacement reaction Methods 0.000 claims abstract description 34
- 238000011068 loading method Methods 0.000 claims abstract description 30
- 239000000463 material Substances 0.000 claims abstract description 27
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 239000012798 spherical particle Substances 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 33
- 230000005484 gravity Effects 0.000 claims description 24
- 230000000694 effects Effects 0.000 claims description 16
- 238000009825 accumulation Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000012856 packing Methods 0.000 claims description 5
- 230000008676 import Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 239000013049 sediment Substances 0.000 claims description 3
- 238000004062 sedimentation Methods 0.000 claims description 3
- 230000003139 buffering effect Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 239000008187 granular material Substances 0.000 claims description 2
- 210000003205 muscle Anatomy 0.000 claims description 2
- 230000035772 mutation Effects 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000005469 granulation Methods 0.000 claims 1
- 230000003179 granulation Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种用于岩土体三维离散元建模和模拟的通用模拟箱,1)建立通用模拟箱,生成箱体6个方向上的指定尺寸边界板,边界板由一系列球体颗粒组成;在6个边界板内侧分别生成尺寸稍小的缓冲边界压力板,建立三维模拟箱;2)箱体内生成颗粒堆积体,在箱体内随机生成无重力颗粒;3)建立不同构筑物的三维模型,导入包括分层、分界坐标文件;对轮廓外的颗粒进行移除,形成构筑物三维模型;4)对三维模型赋岩土体材料性质;导入岩土体材料属性文件,根据转换公式将岩土体参数换算为颗粒参数;对颗粒赋参数,建立具有特定弹性性质和破坏性质的三维离散元模型;5)通过边界压力板进行力、位移、震动的加载通过边界施加位移和震动,通过压力板施加应力作用等荷载,进行迭代运算并输出结果。
Description
技术领域
本发明涉及岩土体三维离散元模拟,尤其是高效的建模方法和多功能三维模拟箱。
背景技术
岩土工程和地质工程领域的很多问题涉及到大变形和破坏,如滑坡、崩塌、地面沉降、隧道开挖、水力压裂等。目前对于岩土体结构物的模拟方法大都基于连续介质力学(如有限元法、有限差分法)。而岩土体在微观上是由一系列的颗粒、孔隙和裂隙组成的结构系统,这些方法难以有效模拟这种离散性,以及大变形和破坏。离散元数值模拟方法基于大量的球体颗粒进行计算,是进行大变形和破坏过程模拟的有效手段。目前,岩土离散元建模和模拟还存在着很大的困难,难以应用于工程实践。以国外商业软件PFC为例,建模方法通常为首先在指定几何体内生成颗粒单元,然后建立连接,赋材料参数,最终建立边界条件和调整、计算。这种建模步骤中每一步都涉及大量的复杂操作,并且难以建立特定结构和特定三轴压力的模型。而传统的通过沉积的建模方法又存在着计算量巨大和耗时长等问题。因此,如何快速地建立起特定结构的颗粒离散元几何模型,并施加有效的三轴应力、位移和震动等边界条件,是实现岩土体离散元建模和模拟的关键问题。
发明内容
为了克服目前离散元建模方法效率低、计算量大、操作复杂等问题,本方法目的是,解决离散元高效建模问题,提出一种用于高效建立离散元三维模型的数值模拟方法,建立通用的三维模拟箱,满足不同用途、不同几何尺寸、不同轮廓的模型要求(如:边坡、基坑、隧道以及三轴试验等等)。该方法避免了复杂的离散元编程建模方法和读取图元建模方法,采用重力沉积和夯实完成紧密堆积体的建立,几何文件和材料文件先后从外部导入,实现几何条件与材料赋值程序的隔离,通过边界压力板加载施加荷载,压力板之间交叉但无作用力。由此实现三维离散元模型的高效建立和精确模拟。
本发明的技术方案是:一种用于岩土体三维离散元建模和模拟的通用模拟箱,其步骤包括:(1)建立通用模拟箱。生成6个方向上的指定尺寸边界板,边界板由一系列球体颗粒组成。在6个边界板内侧分别生成尺寸稍小的缓冲压力板,建立三维模拟箱;(2)箱体内生成颗粒堆积体。在箱体内随机生成无重力颗粒;对颗粒施加随机初速度,通过离散元迭代运算,颗粒运动和相互碰撞至随机位置;对颗粒和上部压力板施加重力,使颗粒自然堆积,为使颗粒快速堆积,重力加速度值由下向上增加;在堆积体上方通过柔性板施加周期荷载,压实堆积模型至指定密度。(3)建立不同构筑物的三维模型。导入几何(包括分层、分界坐标)坐标文件;对轮廓外的颗粒进行移除,形成构筑物三维模型;(4)对三维模型赋岩土体材料性质。导入岩土体材料属性文件,根据转换公式将岩土体参数换算为颗粒参数;对颗粒赋参数,建立具有特定弹性性质和破坏性质的三维离散元模型;对岩土体模型进行重新平衡。(5)通过边界压力板进行力、位移、震动等的加载,压力板之间可交叉但不发生作用力,记录加载过程中的各项指标,进行迭代运算并输出结果。
本发明的有益效果是:实现离散元三维模型的高效建模方法,建立通用的三维模拟箱,在箱体内生成大量的无连接无重力的球颗粒,然后重力堆积和循环夯实形成紧密堆积体,此过程计算量大大减小并且模拟了自然沉积过程;在紧密堆积体形成后导入坐标文件进行轮廓切割,在大量迭代运算前设置边界条件;直接导入材料属性文件,将几何设置和材料设置隔离,提高计算效率;通过压力板有效完成各种荷载的施加,压力板之间可以交叉但不发生作用力。此过程实现了离散元三维模型的通用、高效的建模和模拟方法。
附图说明
图1三维离散元高效模拟箱模拟流程示意图;
图2三维模拟箱中紧密堆积体的建立流程图;
图3赋予模型工况条件——以边坡模型为例;
图4对模型施加不同荷载的流程图。
具体实施方式
本发明方法工作流程如图1-图4所示。
图1为三维离散元模型通用模拟箱建模和模拟方法的实现流程图。
步骤10输入三维模拟箱的几何参数,建立模拟箱的外轮廓,包含外边界板和内层的边界压力板,两种板都是由与样品相近粒径的球颗粒组成,边界板为刚性的,以施加边界位移,压力板为柔性,用于缓冲边界位移和产生压力作用。
步骤11建立紧密堆积体。通过重力沉积和循环夯实建立紧密堆积的堆积体。
步骤12赋几何条件。根据构筑物的几何形状和分层分界条件,读取外部坐标文件切割出几何轮廓,移除轮廓外部颗粒。
步骤13赋予材料性质。将岩土体的材料性质、力学参数通过文本文件的形式导入到模型中,通过转换公式将材料参数转换为颗粒参数。
步骤14平衡模型。在对颗粒赋予实际材料参数后,各层单元的受力和单元刚度发生变化,使得单元受力不平衡。因此,模型需要再进行受力平衡。为防止模型破坏,平衡迭代运算时将单元间连接设为不可断裂,使模型单元在重力和弹力作用下运动至稳定状态,将单元间连接恢复至可断裂的胶结状态。
步骤15通过边界对内部堆积模型施加位移、震动等荷载,通过改变压力板所受体力,实现应力的加载;通过调整各边界的位移进行位移的加载,加载后迭代运算完成模拟,输出结果。
图2为三维模拟箱及紧密堆积体的建立流程图。
步骤20输入模拟箱的外观尺寸参数。建立模拟箱的外轮廓,包含6个外边界板和6个内层的边界压力板,两种板都是由与样品同样粒径的球颗粒组成,边界板为刚性的,上压力板为柔性的,其余方向的压力板为刚性的。板的颗粒直径与内部单元的平均直径相同。
步骤21在模型内部生成自由轻颗粒。设定颗粒半径R1=r·f(1),f(1)为均值1正态分布数组,颗粒无重力为轻颗粒。
步骤22建立随机模型。对每一个颗粒施加随机的初速度,颗粒在初速度的作用下运动并发生颗粒之间的碰撞,最终颗粒运动至随机的位置。
步骤23模拟自然沉积建立堆积体。考虑到较高位置的颗粒其沉积位移较大,作用在颗粒上的重力加速度由箱底的零值线性增加到箱顶的10·g。由于位移与所施加的重力加速度成正比,这种线性的重力加速度设置使所有颗粒几乎同时完成沉积。
步骤24夯实堆积模型。将模型颗粒的重力设为零,而对上压力板颗粒施加周期性的大于10·g的重力加速度,对堆积体进行夯实;设置夯实的次数,循环夯实得到不同密度的紧密堆积体。
图3为赋予模型工况条件流程图。
步骤30建立模型图形坐标。在excel软件中进行图形建模工作;将建立好坐标的excel文件导入。
步骤31切割紧密堆积体。导入坐标文件后,将轮廓以外的颗粒移除;在分层层面上将所有连接键断开。
步骤32建立材料参数。在txt文件中设置不同分层岩土体的材料参数,然后导入txt文件。
步骤33赋颗粒材料参数。将宏观的材料参数通过转换公式转换为颗粒间微观参数,模型具有了特定岩土体的弹性性质和破坏性质。
本发明提出的转换公式如下:模型颗粒单元间考虑正向和切向作用力,并由五个参数来确定:正向劲度系数(Kn)、切向劲度系数(Ks)、颗粒间破坏位移(Xb)、颗粒间抗剪强度(Fs0)和颗粒间摩擦系数(μp)。根据平均力学参数建立三维离散元模型,包括五个宏观力学参数,杨氏模量(E)、泊松比(v)、抗拉强度(Tu)、抗压强度(Cu)和泊松比(μi),颗粒间的力学参数通过以下转换公式来确定:
步骤34平衡模型。由于切割模型和赋材料参数后,颗粒的受力和刚度发生变化,模型需要再进行受力平衡;将模型胶结并设置颗粒间连接不可断裂;模型颗粒在重力和弹力作用下运动至稳定状态。
图4表示对模型施加不同荷载的流程图。
步骤40设置荷载的加载条件。主要分为控制应力加载、控制位移加载、震动荷载的加载等。
步骤41控制应力的加载方法。应力加载通过设置施加在压力板上的力来实现。令压力板仅能在板的法向上移动,并锁定其它方向的自由度。将设定的压力板应力乘以颗粒在压力板上的面积得到每个压力板颗粒的体力,例如在垂直方向上通过控制上部压力板的重力实现力的加载。模拟计算时,忽略压力板之间的力的作用,使三个主方向的压力板不相互影响,从而实现模拟箱中三轴应力的施加。
步骤42控制位移的加载。控制某一方向的边界位移实现位移加载,某一方向的边界发生位移时,挤压或拉动相应的压力板,而压力板再与模型中颗粒产生力的作用,压力板起到力的缓冲作用,从而避免模型中出现应力的突变。
步骤43加载震动荷载。对边界(通常为下边界)施加周期性反复的正向和切向位移,在相应压力板的缓冲下,周期性的位移会转化为压力板的正弦式运动,从而模拟出地震作用。
步骤44记录模型的各项状态指标,输出结果。
Claims (5)
1.一种用于岩土体三维离散元建模和模拟的通用模拟箱,其特征是步骤包括:(1)建立通用模拟箱,生成箱体6个方向上的指定尺寸边界板,边界板由一系列球体颗粒组成;在6个边界板内侧分别生成尺寸稍小的缓冲边界压力板,建立三维模拟箱;(2)箱体内生成颗粒堆积体,在箱体内随机生成无重力颗粒;对颗粒施加随机初速度,通过离散元迭代运算,颗粒运动和相互碰撞至随机位置;对颗粒和上部压力板施加重力,使颗粒自然堆积,为使颗粒快速堆积,重力加速度值由下向上增加;在堆积体上方通过柔性板施加周期荷载,压实堆积模型至指定密度;(3)建立不同构筑物的三维模型,导入包括分层、分界坐标文件;对轮廓外的颗粒进行移除,形成构筑物三维模型;(4)对三维模型赋岩土体材料性质;导入岩土体材料属性文件,根据转换公式将岩土体参数换算为颗粒参数;对颗粒赋参数,建立具有特定弹性性质和破坏性质的三维离散元模型;(5)通过边界压力板进行力、位移、震动的加载,压力板之间可交叉但不发生作用力,记录加载过程中的各项指标,进行迭代运算并输出结果。
2.根据权利要求1所述的通用模拟箱,其特征是三维离散元模型通用模拟箱建模和模拟方法的流程步骤为,步骤10输入三维模拟箱的几何参数,建立模拟箱的外轮廓,包含外边界板和内层的边界压力板,两种板都是由与样品相近粒径的球颗粒组成,边界板为刚性的,以施加边界位移,压力板为柔性,用于缓冲边界位移和产生压力作用;
步骤11建立紧密堆积体;通过重力沉积和循环夯实建立紧密堆积的堆积体;
步骤12赋几何条件;根据颗粒堆积的几何形状和分层分界条件,读取外部坐标文件切割出几何轮廓,移除轮廓外部颗粒;
步骤13赋予材料性质;将岩土体的材料性质、力学参数通过文本文件的形式导入到模型中,通过转换公式将材料参数转换为颗粒参数;
步骤14平衡模型;在对颗粒赋予实际材料参数后,各层单元的受力和单元刚度发生变化,使得单元受力不平衡;因此,模型需要再进行受力平衡;为防止模型破坏,平衡迭代运算时将单元间连接设为不可断裂,使模型单元在重力和弹力作用下运动至稳定状态,将单元间连接恢复至可断裂的胶结状态;
步骤15通过边界对内部堆积模型施加位移、震动等荷载,通过改变压力板所受体力,实现应力的加载;通过调整各边界的位移进行位移的加载,加载后迭代运算完成模拟,输出结果。
3.根据权利要求1所述的通用模拟箱,其特征是三维模拟箱及紧密堆积体的建立流程步骤,步骤20输入模拟箱的外观尺寸参数;建立模拟箱的外轮廓,包含6个外边界板和6个内层的边界压力板,两种板都是由与样品同样粒径的球颗粒组成,边界板为刚性的,上压力板为柔性的,其余方向的压力板为刚性的;板的颗粒直径与内部单元的平均直径相同;
步骤21在模型内部生成自由轻颗粒;设定颗粒半径R1=r·f(1),f(1)为均值1正态分布数组,颗粒无重力为轻颗粒;
步骤22建立随机模型;对每一个颗粒施加随机的初速度,颗粒在初速度的作用下运动并发生颗粒之间的碰撞,最终颗粒运动至随机的位置;
步骤23模拟自然沉积建立堆积体;考虑到较高位置的颗粒其沉积位移较大,作用在颗粒上的重力加速度由箱底的零值线性增加到箱顶的10·g;由于位移与所施加的重力加速度成正比,这种线性的重力加速度设置使所有颗粒几乎同时完成沉积;
步骤24夯实堆积模型;将模型颗粒的重力设为零,而对上压力板颗粒施加周期性的大于10·g的重力加速度,对堆积体进行夯实;设置夯实的次数,循环夯实得到不同密度的紧密堆积体。
4.根据权利要求1所述的通用模拟箱,其特征是赋予模型工况条件流程图;步骤30建立模型图形坐标,在excel软件中进行图形建模工作;将建立好坐标的excel文件导入;
步骤31切割紧密堆积体;导入坐标文件后,将轮廓以外的颗粒移除;在分层层面上将所有连接键断开;
步骤32建立材料参数;在txt文件中设置不同分层岩土体的材料参数,然后导入txt文件;
步骤33赋颗粒材料参数;将宏观的材料参数通过转换公式转换为颗粒间微观参数,模型具有了特定岩土体的弹性性质和破坏性质;
转换公式如下:模型颗粒单元间考虑正向和切向作用力,并由五个参数来确定:正向劲度系数(Kn)、切向劲度系数(Ks)、颗粒间破坏位移(Xb)、颗粒间抗剪强度(Fs0)和颗粒间摩擦系数(μp);根据平均力学参数建立三维离散元模型,包括五个宏观力学参数,杨氏模量(E)、泊松比(v)、抗拉强度(Tu)、抗压强度(Cu)和泊松比(μi),颗粒间的力学参数通过以下转换公式来确定:
<mrow>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msqrt>
<mn>2</mn>
</msqrt>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>5</mn>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mi>E</mi>
<mi>d</mi>
</mrow>
<mrow>
<mn>4</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msqrt>
<mn>2</mn>
</msqrt>
<mi>E</mi>
<mi>d</mi>
</mrow>
<mrow>
<mn>4</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>X</mi>
<mi>b</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
</mrow>
<mrow>
<mn>6</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mi>u</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Fs</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mrow>
<mn>6</mn>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>C</mi>
<mi>u</mi>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mn>2</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<mo>+</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<mi>I</mi>
</mrow>
<mrow>
<mn>2</mn>
<mo>+</mo>
<mn>2</mn>
<mi>I</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>I</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤34平衡模型;由于切割模型和赋材料参数后,颗粒的受力和刚度发生变化,模型需要再进行受力平衡;将模型胶结并设置颗粒间连接不可断裂;模型颗粒在重力和弹力作用下运动至稳定状态。
5.根据权利要求1所述的通用模拟箱,其特征是对模型施加不同荷载的流程,步骤40设置荷载的加载条件;主要分为控制应力加载、控制位移加载、震动荷载的加载等;
步骤41控制应力的加载方法;应力加载通过设置施加在压力板上的力来实现;令压力板仅能在板的法向上移动,并锁定其它方向的自由度;将设定的压力板应力乘以颗粒在压力板上的面积得到每个压力板颗粒的体力,例如在垂直方向上通过控制上部压力板的重力实现力的加载。模拟计算时,忽略压力板之间的力的作用,使三个主方向的压力板不相互影响,从而实现模拟箱中三轴应力的施加。
步骤42控制位移的加载。控制某一方向的边界位移实现位移加载,某一方向的边界发生位移时,挤压或拉动相应的压力板,而压力板再与模型中颗粒产生力的作用,压力板起到力的缓冲作用,从而避免模型中出现应力的突变;
步骤43加载震动荷载;对边界(通常为下边界)施加周期性反复的正向和切向位移,在相应压力板的缓冲下,周期性的位移会转化为压力板的正弦式运动,从而模拟出地震作用;
步骤44记录模型的各项状态指标,输出结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710700733.XA CN107644120B (zh) | 2017-08-16 | 2017-08-16 | 一种用于岩土体三维离散元快速建模的通用数值模拟箱的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710700733.XA CN107644120B (zh) | 2017-08-16 | 2017-08-16 | 一种用于岩土体三维离散元快速建模的通用数值模拟箱的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107644120A true CN107644120A (zh) | 2018-01-30 |
CN107644120B CN107644120B (zh) | 2020-10-09 |
Family
ID=61110713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710700733.XA Active CN107644120B (zh) | 2017-08-16 | 2017-08-16 | 一种用于岩土体三维离散元快速建模的通用数值模拟箱的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107644120B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
CN110502825A (zh) * | 2019-08-19 | 2019-11-26 | 青岛理工大学 | 一种提取三维破裂面的方法 |
CN110705067A (zh) * | 2019-09-23 | 2020-01-17 | 湘潭大学 | 一种考虑复杂随机形状的可破碎颗粒三轴试样生成方法 |
CN110738725A (zh) * | 2019-09-24 | 2020-01-31 | 湖北工业大学 | 一种建筑垃圾离散元三维模型构建方法 |
CN111678845A (zh) * | 2020-04-24 | 2020-09-18 | 绍兴文理学院 | 卵砾石材料特征粒径与抗剪强度参数关系模型的构建方法 |
CN112115608A (zh) * | 2020-09-16 | 2020-12-22 | 中国地质大学(北京) | 复合颗粒粒径计算和级配调整配置方法 |
CN112182703A (zh) * | 2020-09-27 | 2021-01-05 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
CN112699587A (zh) * | 2020-12-31 | 2021-04-23 | 三一专用汽车有限责任公司 | 搅拌筒磨损量预测方法、装置和可读存储介质 |
CN112765895A (zh) * | 2021-01-28 | 2021-05-07 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
CN112818611A (zh) * | 2021-01-28 | 2021-05-18 | 南京大学 | 一种单裂隙岩石水力压裂过程流固耦合的数值模拟方法 |
CN113094792A (zh) * | 2021-04-19 | 2021-07-09 | 湘潭大学 | 一种基于离散元的钢纤维混凝土纤维定向分布的方法 |
CN113515878A (zh) * | 2021-07-07 | 2021-10-19 | 重庆交通大学 | 一种基于块石形状和破碎的堆积体离散元三维建模方法 |
CN113887092A (zh) * | 2021-08-31 | 2022-01-04 | 浙江工业大学 | 土工织物包裹泥浆加排水板抽水的离散元模型构建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1808444A (zh) * | 2005-05-28 | 2006-07-26 | 吉林大学 | 基于cad模型的离散元法边界建模方法 |
JP2015118475A (ja) * | 2013-12-17 | 2015-06-25 | 住友ゴム工業株式会社 | フィラー配合ゴムの有限要素モデルの作成方法 |
CN106446402A (zh) * | 2016-09-22 | 2017-02-22 | 南京大学(苏州)高新技术研究院 | 一种土体失水开裂多场耦合离散元快速模拟建模方法 |
-
2017
- 2017-08-16 CN CN201710700733.XA patent/CN107644120B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1808444A (zh) * | 2005-05-28 | 2006-07-26 | 吉林大学 | 基于cad模型的离散元法边界建模方法 |
JP2015118475A (ja) * | 2013-12-17 | 2015-06-25 | 住友ゴム工業株式会社 | フィラー配合ゴムの有限要素モデルの作成方法 |
CN106446402A (zh) * | 2016-09-22 | 2017-02-22 | 南京大学(苏州)高新技术研究院 | 一种土体失水开裂多场耦合离散元快速模拟建模方法 |
Non-Patent Citations (2)
Title |
---|
CHUN LIU等: "Mechanical properties and energy conversion of 3D close-packed lattice model for brittle rocks", 《COMPUTER & GEOSCIENCES》 * |
刘春等: "岩土体大型三维离散元模拟系统的研发与应用", 《工程地质学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
CN110309536B (zh) * | 2019-05-24 | 2023-07-11 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
CN110502825A (zh) * | 2019-08-19 | 2019-11-26 | 青岛理工大学 | 一种提取三维破裂面的方法 |
CN110705067A (zh) * | 2019-09-23 | 2020-01-17 | 湘潭大学 | 一种考虑复杂随机形状的可破碎颗粒三轴试样生成方法 |
CN110738725A (zh) * | 2019-09-24 | 2020-01-31 | 湖北工业大学 | 一种建筑垃圾离散元三维模型构建方法 |
CN110738725B (zh) * | 2019-09-24 | 2023-05-05 | 湖北工业大学 | 一种建筑垃圾离散元三维模型构建方法 |
CN111678845B (zh) * | 2020-04-24 | 2023-04-07 | 绍兴文理学院 | 卵砾石材料特征粒径与抗剪强度参数关系模型的构建方法 |
CN111678845A (zh) * | 2020-04-24 | 2020-09-18 | 绍兴文理学院 | 卵砾石材料特征粒径与抗剪强度参数关系模型的构建方法 |
CN112115608A (zh) * | 2020-09-16 | 2020-12-22 | 中国地质大学(北京) | 复合颗粒粒径计算和级配调整配置方法 |
CN112115608B (zh) * | 2020-09-16 | 2024-03-15 | 中国地质大学(北京) | 复合颗粒粒径计算和级配调整配置方法 |
CN112182703B (zh) * | 2020-09-27 | 2023-08-25 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
CN112182703A (zh) * | 2020-09-27 | 2021-01-05 | 中电建路桥集团有限公司 | 一种基于离散单元法联合荷载结构法的类岩堆体隧道结构分析方法 |
CN112699587A (zh) * | 2020-12-31 | 2021-04-23 | 三一专用汽车有限责任公司 | 搅拌筒磨损量预测方法、装置和可读存储介质 |
CN112818611B (zh) * | 2021-01-28 | 2023-08-22 | 南京大学 | 一种单裂隙岩石水力压裂过程流固耦合的数值模拟方法 |
CN112765895A (zh) * | 2021-01-28 | 2021-05-07 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
CN112818611A (zh) * | 2021-01-28 | 2021-05-18 | 南京大学 | 一种单裂隙岩石水力压裂过程流固耦合的数值模拟方法 |
CN112765895B (zh) * | 2021-01-28 | 2023-10-17 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
CN113094792A (zh) * | 2021-04-19 | 2021-07-09 | 湘潭大学 | 一种基于离散元的钢纤维混凝土纤维定向分布的方法 |
CN113094792B (zh) * | 2021-04-19 | 2024-04-09 | 湘潭大学 | 一种基于离散元的钢纤维混凝土纤维定向分布的方法 |
CN113515878B (zh) * | 2021-07-07 | 2023-06-20 | 重庆交通大学 | 一种基于块石形状和破碎的堆积体离散元三维建模方法 |
CN113515878A (zh) * | 2021-07-07 | 2021-10-19 | 重庆交通大学 | 一种基于块石形状和破碎的堆积体离散元三维建模方法 |
CN113887092B (zh) * | 2021-08-31 | 2023-03-31 | 浙江工业大学 | 土工织物包裹泥浆加排水板抽水的离散元模型构建方法 |
CN113887092A (zh) * | 2021-08-31 | 2022-01-04 | 浙江工业大学 | 土工织物包裹泥浆加排水板抽水的离散元模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107644120B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107644120A (zh) | 一种用于岩土体三维离散元快速建模和模拟的通用数值模拟箱 | |
Munjiza et al. | A combined finite‐discrete element method in transient dynamics of fracturing solids | |
WO2021169336A1 (zh) | 隧道岩体破坏突涌水灾变模拟的近场动力学方法及系统 | |
Hopkins et al. | On the numerical simulation of the sea ice ridging process | |
US7966165B2 (en) | Soil-water coupled analyzer and soil-water coupled analysis method | |
Long et al. | Nonlinear study on the structure-soil-structure interaction of seismic response among high-rise buildings | |
CN106446402A (zh) | 一种土体失水开裂多场耦合离散元快速模拟建模方法 | |
Firoj et al. | Effect of CPRF on nonlinear seismic response of an NPP structure considering raft-pile-soil-structure-interaction | |
CN105404758A (zh) | 一种基于有限单元法的固体连续介质变形的数值模拟方法 | |
Jia et al. | Coupled three-dimensional discrete element–finite difference simulation of dynamic compaction | |
Gong et al. | A non-matching nodes interface model with radial interpolation function for simulating 2d soil–structure interface behaviors | |
Ding et al. | CFD/DEM coupled approach for the stability of caisson-type breakwater subjected to violent wave impact | |
Manna et al. | Effect of surcharge load on stability of slopes-testing and analysis | |
Gong et al. | The simulation of high compressive stress and extrusion phenomenon for concrete face slabs in CFRDs under strong seismic loads | |
Qu et al. | A flexible various-scale approach for soil-structure interaction and its application in seismic damage analysis of the underground structure of nuclear power plants | |
CN113761760A (zh) | 工程尺度岩体破裂全过程模拟的pd-fem数值计算方法及系统 | |
Xu et al. | Investigation on artificial boundary problem in the combined finite-discrete element method (FDEM) | |
Demir et al. | Numerical investigation of seismic performance of high modulus columns under earthquake loading | |
Salahudeen et al. | Numerical modelling of soil reinforcement using geogrids | |
CN114564899A (zh) | 一种砂土渗流破坏的模拟方法、装置及存储介质 | |
Sunbul et al. | A case study on 3d non-linear analysis of a clay core rockfill dam | |
Liu et al. | Circumventing volumetric locking in stabilized smoothed particle finite element method and its application to dynamic large deformation problems | |
Feng et al. | 2D particle contact-based meshfree method in CDEM and its application in geotechnical problems | |
Nguyen et al. | A multiphase approach for evaluating the horizontal and rocking impedances of pile group foundations | |
Gordan et al. | Soil saturated simulation in embankment during strong earthquake by effect of elasticity modulus |
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 |