CN107784177A - 基于bpm理论构建细观力链颗粒模型的方法 - Google Patents

基于bpm理论构建细观力链颗粒模型的方法 Download PDF

Info

Publication number
CN107784177A
CN107784177A CN201711061238.5A CN201711061238A CN107784177A CN 107784177 A CN107784177 A CN 107784177A CN 201711061238 A CN201711061238 A CN 201711061238A CN 107784177 A CN107784177 A CN 107784177A
Authority
CN
China
Prior art keywords
mover
mrow
msup
msub
overbar
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.)
Pending
Application number
CN201711061238.5A
Other languages
English (en)
Inventor
蔡改贫
姜志宏
郭晋
占鹏飞
郭进山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangxi University of Science and Technology
Original Assignee
Jiangxi University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangxi University of Science and Technology filed Critical Jiangxi University of Science and Technology
Priority to CN201711061238.5A priority Critical patent/CN107784177A/zh
Publication of CN107784177A publication Critical patent/CN107784177A/zh
Priority to JP2018202287A priority patent/JP6789274B2/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

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)
  • Disintegrating Or Milling (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于BPM理论构建细观力链颗粒模型的方法,第一步,将物体的几何模型在空间进行网格划分,并获得每个网格中心的三维坐标;第二步,利用所得到的坐标信息按照网格的数量生成相应数量的刚性基础颗粒单元,并赋予颗粒固有属性;第三步,通过BPM理论将各个刚性基础颗粒进行粘结,相邻颗粒间形成一种粘结键,并给粘结键赋予一定的力学特性;第四步,给所构成的粘结模型施加外力,力沿着粘结键向内部进行传递,粘结颗粒模型受到外力作用时,力沿着粘结键向模型内部传递。根据本方法构建的颗粒模型可以体现出外力加载时,力在物体内部的传递路径演化情况和内部力的传递规律,为从微细观角度研究物体破碎机理提供一定的参考。

Description

基于BPM理论构建细观力链颗粒模型的方法
技术领域
本发明属于岩土工程技术领域,具体涉及一种基于BPM理论构建细观力链颗粒模型的方法,尤其是涉及模拟物质破碎过程的颗粒细观结构模型构建方法。
背景技术
颗粒材料的破碎是一个复杂的变形过程,伴随着微裂纹的滋生、扩展以及交汇贯通,这就需要运用接触理论、细观力学、统计力学等多方面的知识去研究破碎的微观、细观或宏观变化过程,揭示破碎的本质。
在细观力学研究方面,当颗粒材料受到外部力作用时,力在其内部的传递路径被称之为力链。国外学者R.H.Ewoldt通过光弹试验测试法对散体颗粒系统内部颗粒间的接触力进行了检测,国内学者孙其诚等人对散体颗粒体系内部颗粒间的接触力分布特性以及接触力形成的网络进行了研究,给出了颗粒体系受力条件下力链网络的形态演化;杨涵等人通过散体颗粒介质变形局部宏观和细观机制进行研究,将力链结构分为两种:准直性较好的力链和准直性较差的力链;洪俊等人通过对散粒体系统进行动力学分析,得到了散粒系统的变形行为和外部载荷之间的关系;刘恩龙等人通过构建细观力学模型,给出了颗粒集合体破碎过程中的细观应力-应变本构方程;以上细观力学的研究都是针对散体系统颗粒进行的,而实际颗粒材料是由不同的颗粒和胶结物质组成(如混凝土、矿岩等),因此在从细观方面对颗粒破碎进行的研究较少。
发明内容
本发明的目的在于克服现有技术的不足,适应现实的需要,提供一种基于BPM理论构建细观力链颗粒模型的方法。
为了实现本发明的目的,本发明采用的技术方案为:
一种基于BPM理论构建细观力链颗粒模型的方法,包括如下步骤:
第一步,将物体的几何模型在空间进行网格划分,并获得每个网格中心的三维坐标;
第二步,利用所得到的坐标信息按照网格的数量生成相应数量的刚性基础颗粒单元,并赋予颗粒固有属性;
第三步,通过BPM理论将各个刚性基础颗粒进行粘结,相邻颗粒间形成一种粘结键,并给粘结键赋予一定的力学特性,其中的力学关系如下:
分别为法向、切向的力、力矩;ni、ti为接触平面的单位矢量;
随着外力的作用,粘结键在每时步产生相应弹性力和力矩增量为:
式中:A、I、J分别为平行键的初始面积、惯性矩和极矩;
作用在平行键上的最大拉应力和剪应力为:
当最大拉应力超过抗拉强度或是最大剪应力大于剪切强度时,平行键发生破裂失效;
第四步,给所构成的粘结模型施加外力,力沿着粘结键向内部进行传递,根据上述公式结合粘结颗粒模型进行判断,粘结颗粒模型受到外力作用时,力沿着粘结键向模型内部传递,形成的传递路径称之为力链,力链有强弱之分,同时力在内部进行传递的过程中表现出一定的传递规律。
所述第四步中,强弱力链的判断方法如下:
式中:Fi表示粘结键上受到的力(N),为粘结键上受到的平均力(N)。
力链传递过程中的规律统计方法为:
定义无量纲键上受力为
通过统计学的方法,对颗粒系统内部粘结键上受力分布进行统计分析,其概率计算方法下:
式中,Dj为f在区间[fj,fj+1]内的数量,Df为总体样本数。
所述第一步中,粘结键上受到的平均力计算方法为:
式中:N为颗粒体系内部粘结键的总数量。
本发明的有益效果在于:
根据本方法构建的颗粒模型可以体现出外力加载时,力在物体内部的传递路径演化情况和内部力的传递规律,同时能从微细观角度描述不同加载阶段下的裂纹变化,完成物体破碎过程中微细观方面力及裂纹的变化,为从微细观角度研究物体破碎机理提供一定的参考。
附图说明
图1为本案采用的巴西圆盘模型示意图;
图2为图1中巴西圆盘几何体网格划分示意图;
图3圆盘颗粒模型;
图4粘结键空间骨架;
图5钨矿石动态断破碎时粘结键断裂过程;
图6颗粒模型加载过程中力链传递形态;
图7不同应变阶段下粘结键力分布概率
具体实施方式
下面结合附图和实施例对本发明进一步说明:
实施例:参见图1——图7。
以一种钨矿石巴西圆盘颗粒模型的构建方法为例,来详细说明本案的技术方案,本发明公开了一种基于BPM理论构建细观力链颗粒模型的方法,包括如下步骤:
第一步,将物体的几何模型在空间进行网格划分,并获得每个网格中心的三维坐标;如附图1建立巴西圆盘三维几何模型,并按照图2所示对图1中的巴西圆盘模型进行网格划分,并提取每个网格单元中心点的三维坐标;
第二步,利用所得到的坐标信息按照网格的数量生成相应数量的刚性基础颗粒单元,并赋予颗粒固有属性;如图3所示,将每个刚性的基础颗粒单元赋予钨矿石属性,其基本属性如下表1所示:
表1为刚性基础颗粒固有属性
第三步,通过BPM理论将各个刚性基础颗粒进行粘结,相邻颗粒间形成一种粘结键,图4即为各粘结键形成圆盘模型的空间骨架结构示意图,并给粘结键赋予一定的力学特性,其中的力学关系如下:
分别为法向、切向的力、力矩;ni、ti为接触平面的单位矢量;
随着外力的作用,粘结键在每时步产生相应弹性力和力矩增量为:
式中:A、I、J分别为平行键的初始面积、惯性矩和极矩;
作用在平行键上的最大拉应力和剪应力为:
当最大拉应力超过抗拉强度或是最大剪应力大于剪切强度时,平行键发生破裂失效;
本实施例中粘结键表现为线弹性,其中粘结键的相应参数如表2所示:
表2粘结键相应粘结参数
图5表示钨矿石动态断破碎时粘结键断裂过程。图中表示巴西圆盘受到外部力时,内部粘结键空间结构网络的断裂变化,由此体现出不同加载阶段下,微裂纹滋生、扩展以及宏观断裂的演变;从力的角度也可以体现出,力在内部空间结构中的传递路径。
第四步,根据上述公式结合粘结颗粒模型进行判断,粘结颗粒模型受到外力作用时,力沿着粘结键向模型内部传递,形成的传递路径称之为力链,力链有强弱之分,同时力在内部进行传递的过程中表现出一定的传递规律。
强弱力链的判断方法如下:
式中:Fi表示粘结键上受到的力(N),为粘结键上受到的平均力(N)。
力链传递过程中的规律统计方法为:
定义无量纲键上受力为
通过统计学的方法,对颗粒系统内部粘结键上受力分布进行统计分析,其概率计算方法下:
式中,Dj为f在区间[fj,fj+1]内的数量,Df为总体样本数。
粘结键上受到的平均力计算方法为:
式中:N为颗粒体系内部粘结键的总数量。
图6表示颗粒模型加载过程中力链传递形态。图中表示巴西圆盘受到外部力时,力在内部空间结构中的传递路径,力的分布特征以及力链网络的演化规律。
图7表示为巴西圆盘受力时不同应变阶段下力分布特征。图中表示巴西圆盘模型在加载过程中,力的分布概率呈指数衰减特征。
图7中力分布衰减函数如下所示:
P(f)=ae-bf
式中,a、b为拟合参数,f粘结键所受合力与平均合力的比值。
对图7中的数据进行拟合得出拟合参数,如表3所示:
表3指数分布拟合结果
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (3)

1.一种基于BPM理论构建细观力链颗粒模型的方法,包括如下步骤:
第一步,将物体的几何模型在空间进行网格划分,并获得每个网格中心的三维坐标;
第二步,利用所得到的坐标信息按照网格的数量生成相应数量的刚性基础颗粒单元,并赋予颗粒固有属性;
第三步,通过BPM理论将各个刚性基础颗粒进行粘结,相邻颗粒间形成一种粘结键,并给粘结键赋予一定的力学特性,其中的力学关系如下:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <msub> <mi>n</mi> <mi>i</mi> </msub> <mo>+</mo> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <msub> <mi>n</mi> <mi>i</mi> </msub> <mo>+</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
分别为法向、切向的力、力矩;ni、ti为接触平面的单位矢量;
随着外力的作用,粘结键在每时步产生相应弹性力和力矩增量为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <mo>=</mo> <msup> <mover> <mi>k</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <msup> <mi>A&amp;Delta;U</mi> <mi>n</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> <mo>=</mo> <msup> <mover> <mi>K</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <msup> <mi>A&amp;Delta;U</mi> <mi>s</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <mo>=</mo> <mo>-</mo> <msup> <mi>k</mi> <mi>s</mi> </msup> <msup> <mi>J&amp;Delta;&amp;theta;</mi> <mi>n</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> <mo>=</mo> <mo>-</mo> <msup> <mi>k</mi> <mi>n</mi> </msup> <msup> <mi>I&amp;Delta;&amp;theta;</mi> <mi>s</mi> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
式中:A、I、J分别为平行键的初始面积、惯性矩和极矩;
作用在平行键上的最大拉应力和剪应力为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mover> <mi>&amp;sigma;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> </mrow> <mi>A</mi> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>|</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> <mo>|</mo> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <mi>I</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>&amp;tau;</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msup> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msup> </mrow> <mi>A</mi> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>|</mo> <msup> <mover> <mi>M</mi> <mo>&amp;OverBar;</mo> </mover> <mi>n</mi> </msup> <mo>|</mo> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <mi>J</mi> </mfrac> </mtd> </mtr> </mtable> </mfenced>
当最大拉应力超过抗拉强度或是最大剪应力大于剪切强度时,平行键发生破裂失效;
第四步,给所构成的粘结模型施加外力,力沿着粘结键向内部进行传递,根据上述公式结合粘结颗粒模型进行判断,粘结颗粒模型受到外力作用时,力沿着粘结键向模型内部传递,形成的传递路径称之为力链,力链有强弱之分,同时力在内部进行传递的过程中表现出一定的传递规律。
2.根据权利要求1所述的基于BPM理论构建细观力链颗粒模型的方法,其特征在于:所述第四步中,强弱力链的判断方法如下:
式中:Fi表示粘结键上受到的力(N),为粘结键上受到的平均力(N)。
力链传递过程中的规律统计方法为:
定义无量纲键上受力为
<mrow> <mi>f</mi> <mo>=</mo> <mfrac> <msub> <mi>F</mi> <mi>i</mi> </msub> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> </mfrac> </mrow>
通过统计学的方法,对颗粒系统内部粘结键上受力分布进行统计分析,其概率计算方法下:
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>D</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>&lt;</mo> <mi>f</mi> <mo>&lt;</mo> <msub> <mi>f</mi> <mrow> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <msub> <mi>D</mi> <mi>f</mi> </msub> </mfrac> </mrow>
式中,Dj为f在区间[fj,fj+1]内的数量,Df为总体样本数。
3.根据权利要求1所述的基于BPM理论构建细观力链颗粒模型的方法,其特征在于:所述第一步中,粘结键上受到的平均力计算方法为:
<mrow> <mover> <mi>F</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </msubsup> <msub> <mi>F</mi> <mi>i</mi> </msub> </mrow> <mi>N</mi> </mfrac> </mrow>
式中:N为颗粒体系内部粘结键的总数量。
CN201711061238.5A 2017-10-27 2017-10-27 基于bpm理论构建细观力链颗粒模型的方法 Pending CN107784177A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201711061238.5A CN107784177A (zh) 2017-10-27 2017-10-27 基于bpm理论构建细观力链颗粒模型的方法
JP2018202287A JP6789274B2 (ja) 2017-10-27 2018-10-26 Bpm理論に基づくメゾの力チェーン粒子モデルの構築方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711061238.5A CN107784177A (zh) 2017-10-27 2017-10-27 基于bpm理论构建细观力链颗粒模型的方法

Publications (1)

Publication Number Publication Date
CN107784177A true CN107784177A (zh) 2018-03-09

Family

ID=61432683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711061238.5A Pending CN107784177A (zh) 2017-10-27 2017-10-27 基于bpm理论构建细观力链颗粒模型的方法

Country Status (2)

Country Link
JP (1) JP6789274B2 (zh)
CN (1) CN107784177A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362096A (zh) * 2023-03-10 2023-06-30 同济大学 一种颗粒力链动力失效侦测方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111259593B (zh) * 2020-02-10 2023-05-23 江西理工大学 基于edem实现岩石破碎过程可视化的方法
CN111914321B (zh) * 2020-06-09 2024-01-16 西安理工大学 一种堆石混凝土三相细观模型的建立方法
CN114216820B (zh) * 2021-12-13 2024-03-19 青岛理工大学 一种识别三维颗粒体系内部力链的装置及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104615482A (zh) * 2015-02-10 2015-05-13 长安大学 一种冷再生混合料在压实过程中破碎的模拟方法
CN106813973A (zh) * 2016-12-15 2017-06-09 长江水利委员会长江科学院 岩体幂函数型细观时效破裂三维模型的构建方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5663197B2 (ja) * 2010-05-31 2015-02-04 東京電力株式会社 粒子欠損を考慮した個別要素法
US9020784B2 (en) * 2012-04-27 2015-04-28 Livermore Software Technology Corp. Methods for providing a bonded-particle model in computer aided engineering system
JP5632433B2 (ja) * 2012-09-27 2014-11-26 東光株式会社 シミュレーション装置、シミュレーション方法およびプログラム
JP2017059444A (ja) * 2015-09-17 2017-03-23 トヨタ自動車株式会社 全固体電池用電極シミュレーション方法及び装置、並びに全固体電池用電極の製造方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104615482A (zh) * 2015-02-10 2015-05-13 长安大学 一种冷再生混合料在压实过程中破碎的模拟方法
CN106813973A (zh) * 2016-12-15 2017-06-09 长江水利委员会长江科学院 岩体幂函数型细观时效破裂三维模型的构建方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
孙其诚 等: "《静态堆积颗粒中的力链分布》", 《物理学报》 *
孙其诚 等: "《颗粒体系中的骨架及力链网络》", 《岩土力学》 *
林龙飞: "《基于多尺度颗粒模型的预磨机破碎能耗分析》", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 *
郭进山 等: "《一种多尺度内聚颗粒模型的构建方法及数值模拟研究》", 《科学技术与工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362096A (zh) * 2023-03-10 2023-06-30 同济大学 一种颗粒力链动力失效侦测方法
CN116362096B (zh) * 2023-03-10 2023-11-10 同济大学 一种颗粒力链动力失效侦测方法

Also Published As

Publication number Publication date
JP6789274B2 (ja) 2020-11-25
JP2019106175A (ja) 2019-06-27

Similar Documents

Publication Publication Date Title
CN107784177A (zh) 基于bpm理论构建细观力链颗粒模型的方法
Whittles et al. Laboratory and numerical investigation into the characteristics of rock fragmentation
Rong et al. Effect of particle shape on mechanical behaviors of rocks: a numerical study using clumped particle model
Tromans Mineral comminution: energy efficiency considerations
Martin et al. Discrete element simulations of the compaction of aggregated ceramic powders
Metzger et al. Numerical investigation of the breakage of bonded agglomerates during impact
Wang et al. Crashworthiness index of honeycomb sandwich structures under low-speed oblique impact
CN106840877A (zh) 一种基于应力的多轴小裂纹全寿命预测方法
Xia et al. Particle simulation of the failure process of brittle rock under triaxial compression
Li et al. Multi-scale modeling and numerical analysis of sandwich beams with FG auxetic 3D lattice cores and GRC face sheets subjected to drop-weight impacts
CN114510760B (zh) 一种基于离散元的大孔隙混凝土抗压试件模型构建方法
CN110705087A (zh) 一种含封闭应力的岩石离散元试样构建方法
Chen et al. Improving dust resistance of mine tailings using green biopolymer
So et al. Simulation of fabrication and degradation of all-solid-state batteries with ductile particles
Kim et al. Random vibration fatigue analysis of a multi-material battery pack structure for an electric vehicle
Lundkvist et al. A discrete element analysis of the mechanical behaviour of a lithium-ion battery electrode active layer
Li et al. Simulation of flow in sandstone with fluid coupled particle model
CN110750871A (zh) 一种基于离散元的摩擦型颗粒材料蠕变试验模拟方法
Gong et al. Shear failure and mechanical behaviors of granite with discontinuous joints under dynamic disturbance: laboratory tests and numerical simulation
CN114492155A (zh) 一种运用多目标遗传算法优化卸胎机法兰盘结构的方法
Prabhu et al. Effect of particle size on high-strain rate response of sand
Wei et al. Research on aluminum honeycomb buffer device for soft landing on the lunar surface
Chen et al. Mechanical degradation behavior of single crystal LiNixMnyCozO2 cathode in li-ion battery by indentation analysis
Hong et al. Numerical simulation of the collision breakage process between the agglomerate and hammer in a hammer crusher using DEM
CN111157373A (zh) 评估汽车电池多孔结构吸能垫性能的落锤冲击测试法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180309

RJ01 Rejection of invention patent application after publication