JP2011248826A - 粒子欠損を考慮した個別要素法 - Google Patents
粒子欠損を考慮した個別要素法 Download PDFInfo
- Publication number
- JP2011248826A JP2011248826A JP2010124150A JP2010124150A JP2011248826A JP 2011248826 A JP2011248826 A JP 2011248826A JP 2010124150 A JP2010124150 A JP 2010124150A JP 2010124150 A JP2010124150 A JP 2010124150A JP 2011248826 A JP2011248826 A JP 2011248826A
- Authority
- JP
- Japan
- Prior art keywords
- particle
- element method
- individual
- individual element
- elements
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 64
- 230000007547 defect Effects 0.000 title claims abstract description 31
- 238000011439 discrete element method Methods 0.000 title abstract description 4
- 239000002689 soil Substances 0.000 claims abstract description 22
- 230000009471 action Effects 0.000 claims abstract description 7
- 230000009467 reduction Effects 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 45
- 230000007812 deficiency Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 22
- 238000004364 calculation method Methods 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 6
- 238000012669 compression test Methods 0.000 description 6
- 239000004575 stone Substances 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 230000003993 interaction Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
【解決手段】連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、前記個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが設定された限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させる(縮小後半径:R*i)粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素の半径縮小比率[R*i/Ri]を段階的に小さく設定する。
【選択図】図3
Description
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法が提供される。
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きく設定し、個別要素の半径縮小比率を段階的に小さく設定することにより、簡易に土粒子の破砕の影響を再現し得るようにしたものである。
(1)要素の運動方程式
質量mi、慣性モーメントIiのある要素iについての運動方程式は次式で表される。
要素間の相互作用が働く条件として、要素間の接触判定を行う。時刻tにおける要素iと要素jの中心の座標をそれぞれ([X]t i、[Y]t i)、([X]t j、[Y]t j)とし、回転変位量を([θ]t i、[θ]t j)、半径をそれぞれri,rjとした時,下式(2)を満たすと要素iと要素jは接触したものとする。
図12に示されるように、接触した要素間には相互作用力が働く相互作用力の粒子間に弾性バネ52と粘性ダッシュポッド53を仮定している。
図12に示されるように、弾性ばね52によって作用する接線方向の力は,接触した時点での値を0とし下式(6)に従い接線方向力を要素が接触している間,毎ステップ増減させることによって計算を行う。
本発明では、個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させるものである(縮小後半径:R*i)。この際、土粒子の破砕形態は、1回の欠損で破砕が完了するわけではなく、欠損が何度か繰り返されるとともに、この破砕は初期の小さな荷重では小さく欠損し、次の大きな荷重で大きく破砕する傾向を示すとの知見の下、前記粒子欠損操作を、接触力Fiが漸次大きくなる過程で複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素iの半径縮小比率[R*i/Ri]を段階的に小さく設定するものである。
実施例1では、砕石の欠損を個別要素法に考慮するための基礎実験として、砕石の破砕実験を行った。その結果を図5に示す。試験は、アムスラーで試験体を押し潰す方法で行った。砕石は初期の小さな荷重で小さく欠損し、次の大きな荷重で大きく欠損するものが多かった。図中、●は局所欠損を示し、○が大きな破壊を示す。
実施例2では、原地盤材料に対する三軸圧縮試験を対象に本発明のシミュレーション検討を行った。
実際の地震によって地盤沈下が生じた地盤を対象として、本発明の個別要素法を適用してシミュレーションを行った。対象は図9に示されるように、建屋に挟まれた35m区間の地盤である。また、入力地震波を図10に示した。更に、個別要素法における地盤構成要素諸元を下表3に示し、解析パラメータを下表4に示す。なお、初期最大接触力の平均値Foは、地盤を10mの深度毎に分割し、それぞれの領域で求めた。もちろん、全深度の平均値としてもよい。
Claims (3)
- 連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法。 - 前記粒子欠損操作を2段階とする請求項1記載の粒子欠損を考慮した個別要素法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010124150A JP5663197B2 (ja) | 2010-05-31 | 2010-05-31 | 粒子欠損を考慮した個別要素法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010124150A JP5663197B2 (ja) | 2010-05-31 | 2010-05-31 | 粒子欠損を考慮した個別要素法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011248826A true JP2011248826A (ja) | 2011-12-08 |
JP5663197B2 JP5663197B2 (ja) | 2015-02-04 |
Family
ID=45413964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010124150A Expired - Fee Related JP5663197B2 (ja) | 2010-05-31 | 2010-05-31 | 粒子欠損を考慮した個別要素法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5663197B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019106175A (ja) * | 2017-10-27 | 2019-06-27 | 江西理工大学 | Bpm理論に基づくメゾの力チェーン粒子モデルの構築方法 |
JP2019148116A (ja) * | 2018-02-27 | 2019-09-05 | 公益財団法人鉄道総合技術研究所 | バラスト軌道におけるバラスト沈下量の推定方法 |
CN111008472A (zh) * | 2019-12-02 | 2020-04-14 | 东华理工大学 | 一种基于离散元的裂谷盆地伸展过程模拟方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004310381A (ja) * | 2003-04-04 | 2004-11-04 | Taisei Corp | 個別要素法の計算方法 |
JP2008152455A (ja) * | 2006-12-15 | 2008-07-03 | Tokyo Electric Power Co Inc:The | 個別要素法における離散モデルのバネ定数設定方法 |
-
2010
- 2010-05-31 JP JP2010124150A patent/JP5663197B2/ja not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004310381A (ja) * | 2003-04-04 | 2004-11-04 | Taisei Corp | 個別要素法の計算方法 |
JP2008152455A (ja) * | 2006-12-15 | 2008-07-03 | Tokyo Electric Power Co Inc:The | 個別要素法における離散モデルのバネ定数設定方法 |
Non-Patent Citations (2)
Title |
---|
CSNH199800214001; 武正文夫、五味裕司、鈴木章彦: '"個別要素法による粒状体の流動挙動解析"' 石川島播磨技報 第36巻、第6号, 19961101, p.479-482, 石川島播磨重工業株式会社 * |
JPN6014018967; 武正文夫、五味裕司、鈴木章彦: '"個別要素法による粒状体の流動挙動解析"' 石川島播磨技報 第36巻、第6号, 19961101, p.479-482, 石川島播磨重工業株式会社 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019106175A (ja) * | 2017-10-27 | 2019-06-27 | 江西理工大学 | Bpm理論に基づくメゾの力チェーン粒子モデルの構築方法 |
JP2019148116A (ja) * | 2018-02-27 | 2019-09-05 | 公益財団法人鉄道総合技術研究所 | バラスト軌道におけるバラスト沈下量の推定方法 |
CN111008472A (zh) * | 2019-12-02 | 2020-04-14 | 东华理工大学 | 一种基于离散元的裂谷盆地伸展过程模拟方法 |
CN111008472B (zh) * | 2019-12-02 | 2023-11-21 | 东华理工大学 | 一种基于离散元的裂谷盆地伸展过程模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5663197B2 (ja) | 2015-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Numerical simulation of the rock fragmentation process induced by indenters | |
Gao et al. | The application of a modified Voronoi logic to brittle fracture modelling at the laboratory and field scale | |
Wong et al. | Application of the numerical manifold method to model progressive failure in rock slopes | |
Kh et al. | Numerical simulation of particle breakage of angular particles using combined DEM and FEM | |
Choi et al. | Three-dimensional numerical analysis of the rock-cutting behavior of a disc cutter using particle flow code | |
Galouei et al. | Size effect, material ductility and shape of fracture process zone in quasi-brittle materials | |
Wang et al. | Dynamic validation of a discrete element code in modeling rock fragmentation | |
Gao et al. | On the role of pre-existing discontinuities on the micromechanical behavior of confined rock samples: a numerical study | |
Wang et al. | A 3D Voronoi clump based model for simulating failure behavior of brittle rock | |
Das et al. | A mesh-free approach for fracture modelling of gravity dams under earthquake | |
Liu et al. | Hybrid finite-discrete element modelling of asperity degradation and gouge grinding during direct shearing of rough rock joints | |
Liu et al. | Simulation of rock fracture process based on GPU-accelerated discrete element method | |
Xia et al. | Simulation of rock deformation and mechanical characteristics using clump parallel-bond models | |
JP5663197B2 (ja) | 粒子欠損を考慮した個別要素法 | |
Kim et al. | Numerical simulation of concrete abrasion induced by unbreakable ice floes | |
Yu et al. | A dynamic FEM-DEM multiscale modeling approach for concrete structures | |
Leclerc | Effect of packing characteristics on the simulation of elasticity and brittle fracture by the Cohesive Discrete Element Method | |
Zhang et al. | Influence of intermediate principal stress on failure mechanism of hard rock with a pre-existing circular opening | |
Gao et al. | Microscale characterization of fracture growth in increasingly jointed rock samples | |
An et al. | Hybrid finite-discrete element modelling of rock fracture during conventional compressive and tensile strength tests under quasi-static and dynamic loading conditions | |
Zhang et al. | Numerical simulation of particle breakage of granular assemblies in discrete element analyses | |
Kalogeropoulos et al. | Numerical simulation of rock cutting using Yade | |
Kwok et al. | Numerical simulation of strength and deformation behavior of inherently anisotropic rocks | |
Goodman et al. | Kinematics of shear banding in 3D plane strain DEM | |
Fakhimi et al. | Dilation angle in bonded particle simulation of rock |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130530 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140404 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140512 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140707 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20141201 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141208 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5663197 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |