CN103234785A - 集料试样离散元虚拟压实方法 - Google Patents
集料试样离散元虚拟压实方法 Download PDFInfo
- Publication number
- CN103234785A CN103234785A CN2013101098407A CN201310109840A CN103234785A CN 103234785 A CN103234785 A CN 103234785A CN 2013101098407 A CN2013101098407 A CN 2013101098407A CN 201310109840 A CN201310109840 A CN 201310109840A CN 103234785 A CN103234785 A CN 103234785A
- Authority
- CN
- China
- Prior art keywords
- wall
- particle system
- aggregate particle
- aggregate
- deadweight
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种集料试样离散元虚拟压实方法,包括如下步骤:(1)集料生成:在给定空间内,生成随机分布的、互不重叠的集料颗粒体系;(2)自重压实:对生成的集料体系赋予一定高度,使集料颗粒体系在自重作用下自由下落,通过自重完成自重压实;(3)墙体加速压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v1,使集料颗粒体系在墙体压力和自重作用下加速压实;(4)墙体均匀压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v2,使集料颗粒体系在墙体压力和自重作用下加速压实。本发明提供的集料试样离散元虚拟压实方法,可以在未知空隙率情况下压实集料试样,不仅能保证集料体系被充分压实,而且能保证集料体系不会被过度压实。
Description
技术领域
本发明涉及集料力学性能离散元仿真技术,尤其涉及一种集料试样离散元虚拟压实方法。
背景技术
要模拟集料的力学性能,首先需要获得压实后的集料试件。在离散元建模过程中,通常是在已知空隙率的情况下准备试件。如果已知集料的空隙率,则可以根据体积关系计算集料试样的高度,然后将试样压实至这个高度即可。要知道集料压实状态的空隙率,需要一系列的实验室测试,这样既耗时又增加了成本。而且,为了特定研究目的,需要生成由理想化的集料颗粒组成的集料试样,这种试样的空隙率是无法通过实验室测量得到的。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种在未知空隙率情况下集料试样离散元虚拟压实方法,能够快速得到集料离散元试样。
技术方案:为实现上述目的,本发明采用的技术方案为:
集料试样离散元虚拟压实方法,包括如下步骤:
(1)集料生成:在给定空间内,生成随机分布的、互不重叠的集料颗粒体系;
(2)自重压实:对生成的集料体系赋予一定高度,使集料颗粒体系在自重作用下自由下落,通过自重完成自重压实;
(3)墙体加速压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v1,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于设定值f1时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率是否小于等于50%,若判断为否,则重复步骤(3),否则进入步骤(4);
(4)墙体均匀压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v2,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于0时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率相对于前一次是否发生变化,若判断为否,则重复步骤(4)直至集料颗粒体系的空隙率不再发生变化。
有益效果:本发明提供的集料试样离散元虚拟压实方法,可以在未知空隙率情况下压实集料试样,不仅能保证集料体系被充分压实,而且能保证集料体系不会被过度压实。
具体实施方式
下面结合实例对本发明作更进一步的说明。
集料试样离散元虚拟压实方法,包括如下步骤:
(1)集料生成:在给定空间内,生成随机分布的、互不重叠的集料颗粒体系;
(2)自重压实:对生成的集料体系赋予一定高度,使集料颗粒体系在自重作用下自由下落,通过自重完成自重压实;
(3)墙体加速压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v1,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于设定值f1时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率是否小于等于50%,若判断为否,则重复步骤(3),否则进入步骤(4);
(4)墙体均匀压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v2,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于0时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率相对于前一次是否发生变化,若判断为否,则重复步骤(4)直至集料颗粒体系的空隙率不再发生变化。
步骤(3)的具体过程为:①检查墙体上的向上的竖向力;②当墙体上的竖向力大于设定值f1时,令墙体的速度为零,执行一定数目的计算步,使集料体系在墙体压力f1和自重作用下压实;③当墙体上的竖向力小于设定值f1时,给墙体赋予一定的速度v1,运行一个计算步;④重复上述过程,当集料颗粒体系的空隙率小于等于50%时,停止执行这个程序。
步骤(4)的具体过程为:①检查墙体上的向上的竖向力;②当墙体上的竖向力大于0时,令墙体的速度为零,执行一定数目的计算步,使集料体系在墙体压力和自重作用下压实;③当墙体上的竖向力等于0时,给墙体赋予一定的速度v2,运行一个计算步;④重复上述过程,运行足够多的计算步,直到集料颗粒体系的空隙率不再变化为止。
关于集料颗粒体系的空隙率误差计算
在上述步骤中,只有墙体的最后一次运动才会产生过压实。设墙体的速度为vm/s,计算时间步长为ts,压实后试样的高度为hm,则由于过压实对空隙率造成的误差△P可由下式计算。
墙体的速度需要根据计算时间步长确定,保证由于过压实对空隙率造成的误差可以 忽略不计。例如,如果墙体最后一次向下运动0.1mm,压实后试样的高度为7cm,则空隙率的误差仅为0.143%。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.集料试样离散元虚拟压实方法,其特征在于:包括如下步骤:
(1)集料生成:在给定空间内,生成随机分布的、互不重叠的集料颗粒体系;
(2)自重压实:对生成的集料体系赋予一定高度,使集料颗粒体系在自重作用下自由下落,通过自重完成自重压实;
(3)墙体加速压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v1,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于设定值f1时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率是否小于等于50%,若判断为否,则重复步骤(3),否则进入步骤(4);
(4)墙体均匀压实:对集料颗粒体系上方的墙体赋予一定的向下运动速度v2,使集料颗粒体系在墙体压力和自重作用下加速压实,同时对墙体受到的向上的竖向力进行实时检测,当该竖向力大于0时,停止墙体运动;对集料颗粒体系的空隙率进行检测,判断集料颗粒体系的空隙率相对于前一次是否发生变化,若判断为否,则重复步骤(4)直至集料颗粒体系的空隙率不再发生变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310109840.7A CN103234785B (zh) | 2013-03-29 | 2013-03-29 | 集料试样离散元虚拟压实方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310109840.7A CN103234785B (zh) | 2013-03-29 | 2013-03-29 | 集料试样离散元虚拟压实方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234785A true CN103234785A (zh) | 2013-08-07 |
CN103234785B CN103234785B (zh) | 2015-01-21 |
Family
ID=48882836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310109840.7A Expired - Fee Related CN103234785B (zh) | 2013-03-29 | 2013-03-29 | 集料试样离散元虚拟压实方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234785B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983547A (zh) * | 2014-05-23 | 2014-08-13 | 重庆大学 | 一种细长金属管粉料夯实过程粉料颗粒尺度的划分方法 |
CN105512436A (zh) * | 2016-01-15 | 2016-04-20 | 东南大学 | 一种基于离散元的沥青混合料压实模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005047793A (ja) * | 2003-06-25 | 2005-02-24 | Kenji Izumo | ポーラスコンクリート内の空隙管理方法 |
CN101793638A (zh) * | 2009-11-19 | 2010-08-04 | 长安大学 | 模拟现场碾压效果的公路基层材料圆柱体试件制备方法 |
CN102720114A (zh) * | 2012-03-23 | 2012-10-10 | 云南省公路科学技术研究院 | 一种沥青混合料压实度的快速控制方法 |
-
2013
- 2013-03-29 CN CN201310109840.7A patent/CN103234785B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005047793A (ja) * | 2003-06-25 | 2005-02-24 | Kenji Izumo | ポーラスコンクリート内の空隙管理方法 |
CN101793638A (zh) * | 2009-11-19 | 2010-08-04 | 长安大学 | 模拟现场碾压效果的公路基层材料圆柱体试件制备方法 |
CN102720114A (zh) * | 2012-03-23 | 2012-10-10 | 云南省公路科学技术研究院 | 一种沥青混合料压实度的快速控制方法 |
Non-Patent Citations (5)
Title |
---|
CHRISTOPHE L.MARTIN 等: "Discrete Element Simulations of the Compaction of Aggregated Ceramic Powders", 《JOURNAL OF THE AMERICAN CERAMIC SOCIETY》 * |
D.HE 等: "Computer simulation of powder compaction of spherical particles", 《JOURNAL OF MATERIALS SCIENCE LETTERS》 * |
OLLE SKRINJAR 等: "On discrete element modelling of compaction of powders with size ratio", 《COMPUTATIONAL MATERIALS SCIENCE》 * |
刘丽萍 等: "山区土石料压实标准及其应用", 《土木工程学报》 * |
高涛涛 等: "计算机模拟热拌沥青混合料的压实", 《中外公路》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983547A (zh) * | 2014-05-23 | 2014-08-13 | 重庆大学 | 一种细长金属管粉料夯实过程粉料颗粒尺度的划分方法 |
CN103983547B (zh) * | 2014-05-23 | 2016-06-29 | 重庆大学 | 一种细长金属管粉料夯实过程粉料颗粒尺度的划分方法 |
CN105512436A (zh) * | 2016-01-15 | 2016-04-20 | 东南大学 | 一种基于离散元的沥青混合料压实模拟方法 |
CN105512436B (zh) * | 2016-01-15 | 2018-05-08 | 东南大学 | 一种基于离散元的沥青混合料压实模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103234785B (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102567575A (zh) | 航天器虚拟正弦振动试验方法 | |
CN102564855A (zh) | 一种级配碎石动三轴试验的数值方法 | |
CN102222126A (zh) | 一种级配碎石三轴试验的数值模拟方法 | |
CN102230870A (zh) | 一种级配碎石cbr数值试验方法 | |
CN103234785A (zh) | 集料试样离散元虚拟压实方法 | |
CN105547831A (zh) | 一种测定含复杂结构面岩体变形性质的模型试验方法 | |
CN106354016A (zh) | 一种金属管夯实高度的控制方法 | |
Seck et al. | Dilepton signature of a first-order phase transition | |
CN107966311B (zh) | 基于加速栅数据的极小子样离子推力器可靠度确定方法 | |
CN103163033A (zh) | 级配碎石回弹模量数值测定方法 | |
CN107894367A (zh) | 一种输电线路湿陷性黄土本构模拟方法 | |
CN105987823A (zh) | 通信装备基于振动应力变换的加速试验方法 | |
Caliez et al. | Viscoelastic plastic model and experimental validation for a granular energetic material | |
CN104062563B (zh) | 一种高电压虚拟实验平台的实现方法及系统 | |
Xu et al. | Hydro-mechanical triaxial behaviour of compacted earth at different temperatures | |
CN103308385A (zh) | 级配碎石动态模量数值测定方法 | |
Miao et al. | Microcrack growth behavior and cwfs criterion parameters optimization of granite with PFC | |
WO2022051030A3 (en) | Generative modeling of quantum hardware | |
Royo et al. | Comparison of laboratory direct shear test results with the numerical analysis | |
Baumketner et al. | Finite-size dependence of the bridge function extracted from molecular dynamics simulations | |
Fenton et al. | SHOCK‐LESS HIGH RATE COMPACTION OF POROUS BRITTLE MATERIALS | |
Bo et al. | Mechanical characteristics optimization of ballast bed variable on operating process for the dynamic track stabilizer | |
Sazzad et al. | Modeling the behavior of granular materials under plane strain compression by 3D DEM | |
Sand et al. | DEM modelling and simulation of coal compaction by stamping | |
CN106007484A (zh) | 一种软岩相似模型试验材料及其制备方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 210093 Nanjing University Science Park, 22 Hankou Road, Gulou District, Nanjing City, Jiangsu Province Patentee after: SOUTHEAST University Address before: 211103 No. 5 Runfa Road, Jiangning District, Nanjing City, Jiangsu Province Patentee before: Southeast University |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150121 |