JP5515955B2 - 粒子充填構造シミュレーション方法 - Google Patents
粒子充填構造シミュレーション方法 Download PDFInfo
- Publication number
- JP5515955B2 JP5515955B2 JP2010076937A JP2010076937A JP5515955B2 JP 5515955 B2 JP5515955 B2 JP 5515955B2 JP 2010076937 A JP2010076937 A JP 2010076937A JP 2010076937 A JP2010076937 A JP 2010076937A JP 5515955 B2 JP5515955 B2 JP 5515955B2
- Authority
- JP
- Japan
- Prior art keywords
- particles
- particle
- space
- simulation method
- ellipsoid
- 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.)
- Expired - Fee Related
Links
- 239000002245 particle Substances 0.000 title claims description 80
- 238000000034 method Methods 0.000 title claims description 26
- 238000004088 simulation Methods 0.000 title claims description 25
- 238000012856 packing Methods 0.000 title description 6
- 238000007906 compression Methods 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 8
- 239000011149 active material Substances 0.000 description 11
- 230000008569 process Effects 0.000 description 7
- 239000007787 solid Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 239000000843 powder Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 239000003792 electrolyte Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000000329 molecular dynamics simulation Methods 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 239000007784 solid electrolyte Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Battery Electrode And Active Subsutance (AREA)
Description
本発明は、ある空間内に充填された粒子の充填構造に係るシミュレーションを行う粒子充填構造シミュレーション方法の技術分野に関する。
この種の方法として、例えば、粒子集合体を所定の大きさの空間内に配置し、重心が低い粒子から自由落下させ、その落下プロセスで、水平方向の自由度、又は既設置粒子を透過する条件として粒子透過係数を与え、最下部に粒子を設置して粒子充填構造をシミュレートする方法が提案されている（特許文献１参照）。或いは、粒子間接触、並びに静電的及び磁気的相互作用を考慮した個別要素法に基づいた紛体挙動計算方法が提案されている（特許文献２参照）。
ところで、例えばハイブリッド自動車、電気自動車等の車両に搭載される電池の安全性の向上を図るため、例えば固体電解質等を用いた全固体電池の開発が進められている。該全固体電池は、例えば活物質と電解質粒子とが型に入れられた後に、押し固められることによって作成される。このような紛体や粒体等の粒子を取り扱う分野では、粒子の挙動を把握することが重要な課題の一つであり、該粒子の挙動の把握には、上述したようなシミュレーションが利用されることが多い。
しかしながら、上記特許文献１では、粒子間に働く反発力については考慮されていないという技術的問題点がある。また、上記特許文献２では、粒子間に働く相互作用を計算する際に粒子を点として扱っているため、十分な精度を得られない可能性があるという技術的問題点がある。
本発明は、例えば上記問題点に鑑みてなされたものであり、全固体電池を構成する粒子の充填構造を再現可能な粒子充填構造シミュレーション方法を提案することを課題とする。
本発明の粒子充填構造シミュレーション方法は、上記課題を解決するために、複数の粒子各々の表面に複数の質点を配置する質点配置工程と、前記複数の粒子を相互に重ならないように空間に夫々配置する粒子配置工程と、前記複数の粒子各々における前記複数の質点の相互間の相対位置が固定された状態で、前記複数の粒子の相互間に働く反発力を設定しつつ、前記空間の圧縮を行う圧縮工程とを備える。
本発明の粒子充填構造シミュレーション方法によれば、質点配置工程において、複数の粒子各々の表面に複数の質点が配置される。複数の質点の各々は、粒子の表面に任意に配置されてもよいが、等間隔に配置されることが望ましい。また、粒子の表面に配置される質点の個数は、後述する圧縮工程において、粒子同士の重なりを回避するためにある程度多いことが望ましい。
粒子配置工程において、複数の粒子が相互に重ならないように空間に夫々配置される。ここで、「空間」は、複数の粒子各々の体積の合計である総体積よりも十分大きな容積を有するように設定される。尚、複数の粒子が相互に重ならないようにするには、例えば、複数の粒子各々の直径（粒子が楕円体である場合は、該楕円体の長軸）のうち最も長い直径よりも長い間隔で、複数の粒子を等間隔で空間に夫々配置すればよい。
圧縮工程において、複数の粒子各々における複数の質点の相互間の相対位置が固定された状態で、複数の粒子の相互間に働く反発力が設定されつつ、前記空間の圧縮が行なわれる。ここで、「複数の粒子各々における複数の質点」とは、同一粒子の表面に配置された複数の質点を意味する。また、「複数の粒子各々における複数の質点の相互間の相対位置が固定された状態」とは、粒子が変形しないこと（即ち、粒子を剛体として扱うこと）を意味する。
尚、「複数の粒子の相互間に働く反発力」は、例えば、一の粒子の表面に配置された一の質点と、他の粒子の表面に配置された一の質点との間に、一の粒子と他の粒子とが相互に近づくと反発する２体力間ポテンシャルを設定すればよい。また、空間の圧縮は、例えば分子動力学法等により行えばよい。
本願発明者の研究によれば、以下の事項が判明している。即ち、圧粉全固体電池は、活物質と電解質粒子とを型に入れ圧力を加え押し固め、所定の密度とすることによって作成される。圧粉全固体電池が、適切に充放電するためには、活物質と電極との間に導電性が確保されている必要があり、活物質粒子間の結合性が重要となる。活物質粒子間の結合性は、粒子形状や密度に影響を受ける。この影響を実験で解析するためには膨大な時間と費用がかかるため、該影響の解析にはシミュレーションが用いられることが多い。しかしながら、圧縮後の空間に所定の密度まで粒子を配置しようとすると、粒子同士の重なりの発生等により、実際の密度を再現することが困難である。
しかるに本発明では、実際の圧粉全固体電池の製造工程を模擬して、粒子配置工程において、複数の粒子が相互に重ならないように空間に夫々配置され、その後、圧縮工程において、複数の粒子各々における複数の質点の相互間の相対位置が固定された状態で、複数の粒子の相互間に働く反発力が設定されつつ、前記空間の圧縮が行われる。ここで、質点配置工程において、予め、複数の粒子各々の表面に複数の質点が配置されているので、空間の圧縮が行われる際に、粒子同士が相互に重なることを回避することができる。このように、本発明では、実際の圧粉全固体電池の製造工程を模擬すると共に、粒子の重なりを回避した上で粒子間反発力を考慮したシミュレーションを行っているので、全固体電池を構成する粒子の充填構造を精度よく再現することができる。
本発明の作用及び他の利得は次に説明する実施するための形態から明らかにされる。
以下、本発明に係る粒子充填構造シミュレーション方法の実施形態を、図面に基づいて説明する。
先ず、本実施形態に係る粒子充填構造シミュレーション方法の全体の処理の流れについて、図１及び図２を参照して説明する。図１は、本実施形態に係るシミュレーション方法の処理を示すフローチャートであり、図２は、本実施形態に係るシミュレーション方法の主要な処理の概念を示す概念図である。
図１において、先ず、圧粉全固体電池を構成する活物質を模擬した、本発明に係る「粒子」の一例としての楕円体１００（図２参照）が複数個発生させられる（ステップＳ１０１）。ここで、楕円体１００は、例えば、楕円体１００の長軸半径、該長軸半径の標準偏差、楕円体１００のアスペクト比、該アスペクト比の標準偏差、楕円体１００の個数等が指定されることにより発生させられる。尚、楕円体１００の長軸半径等の各種パラメータの値は、実際の活物質を測定して求めてもよいし、解析者（即ち、ユーザ）が恣意的に決定してもよい。
次に、図２（ａ）に示すように、質点配置工程において、楕円体１００の表面に、例えば等間隔で複数の質点２００が夫々配置される（ステップＳ１０２）。配置された質点２００は、後述する圧縮工程において使用される。尚、圧縮時の楕円体１００同士の重なりを回避するためには、質点２００の相互間の間隔は狭いほうがよいが、その分質点２００の個数が増加するので、圧縮工程に費やされる時間が長くなる。このため、質点２００の個数は、例えば、要求されるシミュレーション結果の精度、解析時間等に応じて調整すればよい。
次に、図２（ｂ）に示すように、粒子配置工程において、複数の楕円体１００各々が、空間１内に、相互に重ならないように配置される（ステップＳ１０３）。尚、楕円体１００同士の重なりの発生を回避するためには、例えば、ステップＳ１０１で発生した複数の楕円体１００各々の長軸半径のうち、最も長い長軸半径よりも長い間隔で空間１内に等間隔で、複数の楕円体１００を夫々配置すればよい。
次に、図２（ｃ）に示すように、圧縮工程において、複数の楕円体１００各々の質点２００の相対位置が固定された状態で（即ち、楕円体１００を剛体とする）、質点２００間に反発の２点間ポテンシャルが設定され、分子動力学法で、空間１が所定の密度となる空間１´まで圧縮される（ステップＳ１０４）。ここで、圧縮後の空間１´の容積は、例えば、複数の楕円体１００各々の体積の合計である楕円体総体積と、予め設定された体積密度とにより決定される。尚、圧縮工程では、楕円体１００の運動速度が低下され、楕円体１００同士の重なりを回避するために系の温度が冷却される。
上述したステップＳ１０１乃至Ｓ１０４の処理の結果、例えば図３に示すような活物質モデルが得られる。図３は、シミュレーション結果の一例を示す図である。
次に、本実施形態に係る粒子充填構造シミュレーションの有効性について、図４を参照して説明する。図４は、粒子の体積密度と結合粒子の割合との関係の一例を示す図である。
図４に示すように、結合粒子の割合が、粒子の体積密度が３０〜４５％の間で急激に変化していることがわかる。ここで、粒子の結合理論によれば、結合性は、体積密度を増加させた場合、ある閾値で急激に増加すると言われている（パーコレーションの閾値）。図４では、その現象が再現されていることから、本実施形態に係る粒子充填構造シミュレーションは、有効に機能していることがわかる。
尚、活物質モデルの結合性は、例えば次のように評価すればよい。即ち、先ず、空間１´全体をメッシュで分割する。次に、各メッシュ内に楕円体１００が共存していた場合、その楕円体１００は結合していると判定する。そして、全メッシュで結合判定を行えば、結合している楕円体１００群を抽出することができる。メッシュが楕円体１００に対して十分小さければ、このような方法で結合性の評価が可能であることが、本願発明者の研究により判明している。
結合した活物質群を模擬した楕円体１００群が、例えば電極を空間１´の上面と仮定して、該上面と接していれば有効な活物質と言える。上記方法で電極（即ち、空間１´の上面）と接触している楕円体１００のみを抽出すると、例えば図５に示すようになる。図５は、シミュレーション結果に基づいて評価された電極と接触している粒子群の一例を、体積密度毎に示す図である。
本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う粒子充填構造シミュレーション方法もまた本発明の技術的範囲に含まれるものである。
１、１´…空間、１００…楕円体、２００…質点
Claims (1)
- 複数の粒子各々の表面に複数の質点を配置する質点配置工程と、
前記複数の粒子を相互に重ならないように空間に夫々配置する粒子配置工程と、
前記複数の粒子各々における前記複数の質点の相互間の相対位置が固定された状態で、前記複数の粒子の相互間に働く反発力を設定しつつ、前記空間の圧縮を行う圧縮工程と
を備えることを特徴とする粒子充填構造シミュレーション方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010076937A JP5515955B2 (ja) | 2010-03-30 | 2010-03-30 | 粒子充填構造シミュレーション方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010076937A JP5515955B2 (ja) | 2010-03-30 | 2010-03-30 | 粒子充填構造シミュレーション方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2011210526A JP2011210526A (ja) | 2011-10-20 |
| JP5515955B2 true JP5515955B2 (ja) | 2014-06-11 |
Family
ID=44941361
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010076937A Expired - Fee Related JP5515955B2 (ja) | 2010-03-30 | 2010-03-30 | 粒子充填構造シミュレーション方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5515955B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104089864A (zh) * | 2014-06-30 | 2014-10-08 | 河海大学 | 一种计算凸形多面体粒子堆积体系孔径分布的方法 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015520927A (ja) * | 2012-05-25 | 2015-07-23 | バシウム・カナダ・インコーポレーテッド | リチウム電気化学的セルのための電極材料 |
| JP7291531B2 (ja) * | 2019-04-26 | 2023-06-15 | Fdk株式会社 | 粒子設計プログラム、設計支援プログラム、検査方法および極板の製造方法 |
-
2010
- 2010-03-30 JP JP2010076937A patent/JP5515955B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104089864A (zh) * | 2014-06-30 | 2014-10-08 | 河海大学 | 一种计算凸形多面体粒子堆积体系孔径分布的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2011210526A (ja) | 2011-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Miehe et al. | A phase‐field model for chemo‐mechanical induced fracture in lithium‐ion battery electrode particles | |
| Parteli | DEM simulation of particles of complex shapes using the multisphere method: application for additive manufacturing | |
| JP5515955B2 (ja) | 粒子充填構造シミュレーション方法 | |
| Wang et al. | A multi-sphere based modelling method for maize grain assemblies | |
| Wang et al. | Modelling and simulation of electron-rich effect on Li diffusion in group IVA elements (Si, Ge and Sn) for Li ion batteries | |
| Gustafsson et al. | Multi-particle finite element modelling of the compression of iron ore pellets with statistically distributed geometric and material data | |
| CN102298660B (zh) | 一种离散元法边界建模的通用方法 | |
| Siiriä et al. | Particle packing simulations based on Newtonian mechanics | |
| Li et al. | Fluid-solid interaction simulation for particles and walls of arbitrary polygonal shapes with a coupled LBM-IMB-DEM method | |
| JP2017059444A (ja) | 全固体電池用電極シミュレーション方法及び装置、並びに全固体電池用電極の製造方法 | |
| JP2020184536A (ja) | リチウムイオン電池正極材料導電性能シミュレーション生成方法 | |
| Xi et al. | Mechanical response analysis of battery modules under mechanical load: experimental investigation and simulation analysis | |
| JP2020194706A5 (ja) | ||
| Cho et al. | Modeling of stresses and strains during (de) lithiation of Ni3Sn2-coated nickel inverse-opal anodes | |
| Shishvan et al. | Cracking and associated volumetric expansion of NMC811 secondary particles | |
| Xiong et al. | Effect of the rolling friction coefficient on particles’ deposition morphology on single fibre | |
| Kumbhar et al. | Mesoscale analysis of Li‐ion battery microstructure using sequential coupling of discrete element and finite element method | |
| Boivin et al. | Breathing of a Silicon-Based Anode: Mechanical Discrete Approach Using DEM | |
| JP2013054418A (ja) | 粒子充填構造シミュレーション方法 | |
| Thakur et al. | Scaling of discrete element model parameters in uniaxial test simulation | |
| Dutt et al. | Granular dynamics simulations of the effect of grain size dispersity on uniaxially compacted powder blends | |
| Gregorski | High green density metal parts by vibrational compaction of dry powder in the three-dimensional printing process. | |
| JP2014016724A (ja) | 充填構造シミュレーション方法、充填構造シミュレーション装置、及びコンピュータプログラム | |
| JP6116315B2 (ja) | 全固体電池用材料及び全固体電池用材料の製造方法 | |
| Akhmetova et al. | Microstructural features of ternary powder compacts |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121227 |
|
| 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: 20140304 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140317 |
|
| LAPS | Cancellation because of no payment of annual fees |