CN109815599A - 一种离散元材料自动训练方法 - Google Patents
一种离散元材料自动训练方法 Download PDFInfo
- Publication number
- CN109815599A CN109815599A CN201910080559.2A CN201910080559A CN109815599A CN 109815599 A CN109815599 A CN 109815599A CN 201910080559 A CN201910080559 A CN 201910080559A CN 109815599 A CN109815599 A CN 109815599A
- Authority
- CN
- China
- Prior art keywords
- particle
- vector
- discrete element
- setting value
- actual value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000000463 material Substances 0.000 title claims abstract description 35
- 238000012549 training Methods 0.000 title claims abstract description 27
- 239000002245 particle Substances 0.000 claims abstract description 44
- 230000006835 compression Effects 0.000 claims abstract description 12
- 238000007906 compression Methods 0.000 claims abstract description 12
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 230000000452 restraining effect Effects 0.000 claims abstract 2
- 238000012795 verification Methods 0.000 claims abstract 2
- 230000005484 gravity Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 claims description 7
- 238000011439 discrete element method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000012856 packing Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000012669 compression test Methods 0.000 claims description 2
- 239000003292 glue Substances 0.000 claims 1
- 238000005303 weighing Methods 0.000 claims 1
- 238000004088 simulation Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种离散元材料自动训练方法，包括以下步骤：步骤A：建立颗粒随机堆积模型；步骤B：输入颗粒宏观力学向量的设定值K0，引入宏观力学向量的中间值K，令K＝K0；步骤C：通过中间值K转换得到颗粒微观力学参数，将微观力学参数赋予模型中的颗粒；步骤D：对赋值后的堆积模型进行拉伸和压缩测试，获得颗粒宏观力学向量的实际值K1；步骤E：判断实际值K1是否收敛于设定值K0；若不收敛，利用如下规则更新中间值K：K＝(K0/K1)·K0并重复步骤C～E，更新颗粒的微观力学参数，直到实际值K1收敛于设定值K0。本发明提供的离散元材料自动训练方法的优点在于：通过离散元堆积模型力学性质测试，调整颗粒间力学参数，实现自动训练具有特定力学性质的离散元材料。
Description
技术领域
本发明涉及离散元数值模拟的材料赋值技术领域，尤其涉及一种离散元材料自动训练方法。
背景技术
岩土工程、地质工程和能源工程领域的很多问题都涉及到大变形和破坏，离散元数值方法是模拟这些问题的一种有效方法。离散元法通过随机堆积不同大小的颗粒来表示岩土体材料，但是离散元定量化建模通常依赖于一系列人工的调节操作和数值模拟测试，以确定合适的颗粒力学参数，所以其建模过程非常繁琐且需要丰富的经验，极大地限制了离散元法的实际工程应用。因此如何在不需要大量经验的情况下确定颗粒参数成为离散元法广泛应用的制约条件。
发明内容
本发明所要解决的技术问题在于提供一种能够自动对颗粒进行训练获得具有特定力学性质的离散元材料的方法。
本发明是通过以下技术方案解决上述技术问题的：
一种离散元材料自动训练方法，包括以下步骤：
步骤A：建立颗粒随机堆积模型；
步骤B：输入颗粒宏观力学向量的设定值K0，引入宏观力学向量的中间值K，令K＝K0；
步骤C：通过中间值K转换得到颗粒微观力学参数，将微观力学参数赋予模型中的颗粒；
步骤D：对赋值后的堆积模型进行拉伸和压缩测试，获得颗粒宏观力学向量的实际值K1；
步骤E：判断实际值K1是否收敛于设定值K0；若不收敛，利用如下规则更新中间值K：
K＝(K0/K1)·K0
并重复步骤C～E，更新颗粒的微观力学参数，直到实际值K1收敛于设定值K0。
优选地，步骤A所述的建立颗粒随机堆积模型的方法为根据模型尺寸建立由六块压力板拼合构成的长方体模型箱，模型箱上方的压力板能够沿竖直方向运动，输入颗粒的平均半径和半径范围，在模型箱内随机生成单元，给单元施加随机初速度，对无序运动的单元赋予具有步骤C中微观参数的材料特征使其具有重力和边界范围，赋值后的单元在模型箱内随机堆积，然后移除模型箱侧面的压力板，平衡模型能量和受力。
优选地，所述单元为单个颗粒或由多个颗粒构成的clump团簇，单元在模型箱中的排列方式为将模型箱内部划分为规则排列的格点，相邻格点的距离等于颗粒的平均半径，在每个格点位置随机生成一个单元，每个格点处的单元的种类和尺寸随机，单元在模型箱内重力堆积后，对团簇单元进行胶结。
优选地，所述平衡模型能量和受力的方法为基于牛顿力学方法计算初始时间每个单元受到其他单元和模型箱边界的作用力以及重力的合力；基于合力计算加速度并结合当前速度求得ΔT时间后的初始速度，利用ΔT时间内的平均速度计算ΔT时间内的单元位移，以ΔT时刻的位置和初始速度重复上述过程，从而重复迭代实现离散元法动态模拟，直到系统稳定。
优选地，时间步ΔT满足
其中，T为单元振动周期，f为一个周期内的迭代频率，f＝50。
优选地，步骤B所述的宏观力学向量的设定值K0为：
K0＝[E0,v0,Tu0,Cu0,μi0]
其中，E0为杨氏模量设定值，v0为泊松比设定值，Tu0为抗拉强度设定值，Cu0为抗压强度设定值，μi0为内摩擦系数设定值。
优选地，步骤C所述的宏观力学性质向量与微观力学参数的转换方法为：
其中，d是四面体晶格单元的边缘长度，p是正压力；宏观力学参数包括杨氏模量E、泊松比v、抗拉强度Tu、抗压强度Cu和内摩擦系数μi；微观力学参数包括法向刚度Kn、切向刚度Ks、断裂位移Xb、抗剪力Fs和颗粒摩擦系数μp。
优选地，步骤D所述的获取颗粒宏观力学向量的实际值K1的方法为：
根据宏观力学向量的中间值K对应的杨氏模量E、泊松比v、抗压强度Cu确定对模型箱上方压力板的压力进行单轴压缩试验，通过应力应变曲线得到杨氏模量实际值E1，通过侧向应变和竖向应变的比值得到泊松比实际值v1，通过应力应变曲线峰值确定抗压强度实际值Cu1；
根据求出内摩擦系数的实际值μi1。
根据中间值K对应的抗拉强度Tu确定加压大小，对模型箱上方的压力板进行单轴拉伸试验，根据应力应变曲线确定模型破坏时的抗拉强度实际值Tu1。
优选地，步骤E中判断宏观力学向量的实际值K1是否收敛于初始宏观力学性质K0的方法为：引入收敛系数，
σE＝E0/E1
σv＝v0/v1
如果收敛系数在预设区间内，则判断实际值K1收敛于设定值K0，训练过程结束，得到具有接近宏观力学向量设定值的离散元材材料。
优选地，当所有宏观力学性质参数的收敛系数同时满足σ∈[0.9,1.1]时，认为宏观力学向量的实际值K1收敛于宏观力学向量的设定值K0。
本发明提供的离散元材料自动训练方法的优点在于：利用转换公式建立模型宏观力学性质和微观力学参数之间的关系，通过离散元堆积模型力学性质测试，自动调整颗粒间力学参数，实现自动训练具有特定力学性质的离散元材料，可显著提高离散元软件的建模效率，便于离散元法的推广应用。
附图说明
图1是本发明的实施例所提供的离散元材料自动训练方法的流程图；
图2是本发明的实施例所提供的离散元材料自动训练方法的模型箱示意图；
图3是本发明的实施例所提供的离散元材料自动训练方法的建模流程图；
图4是本发明的实施例所提供的离散元材料自动训练方法的平衡颗粒速度的方法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。
如图1所示，一种离散元材料自动训练方法，包括以下步骤：
步骤A：建立颗粒随机堆积模型；
结合图2和图3，根据模型尺寸建立由六块压力板拼合构成的长方体模型箱，模型箱上方的压力板能够沿竖直方向运动，输入颗粒的平均半径和半径范围，将模型箱内部划分为规则排列的格点，相邻格点的距离等于颗粒的平均半径，在每个格点位置随机生成一个单元，所述单元为单个颗粒或由多个颗粒构成的clump团簇，每个格点处的单元的种类和尺寸随机，给每个单元随机施加不同方向和大小的初速度，单元运动一段时间后对无序运动的单元赋予具有特定微观参数的材料特征使其具有重力和边界范围，赋值后的单元在模型箱内随机堆积，对团簇单元进行胶结从而在计算模型参数时忽略其内部颗粒的互相作用；然后移除模型箱侧面的压力板，平衡模型箱内部的能量和受力即可。
所述平衡模型内部能量和速度的方法为本领域现有技术，具体方案可参考论文(刘春,张晓宇,许强,朱晨光,汤强.三维离散元模型的滑坡能量守恒模拟研究[J].地下空间与工程学报,2017,13(S2):698-704)。
参考图4，本申请对具体过程进行简单的描述，基于牛顿力学方法计算初始时间每个单元受到其他单元和模型箱边界的作用力以及重力的合力；基于合力计算加速度并结合当前速度求得ΔT时间后的初始速度，利用ΔT时间内的平均速度计算ΔT时间内的单元位移，以ΔT时刻的位置和初始速度重复上述过程，从而重复迭代实现离散元法动态模拟，直到系统稳定。
所述时间步ΔT满足：
其中，T为单元振动周期，f为一个周期内的迭代频率，优选实施例中f＝50。
步骤B：输入颗的粒宏观力学向量的设定值K0，
K0＝[E0,v0,Tu0,Cu0,μi0]
其中，E0为杨氏模量设定值，v0为泊松比设定值，Tu0为抗拉强度设定值，Cu0为抗压强度设定值，μi0为内摩擦系数设定值；
引入宏观力学向量的中间值K，令K＝K0。
步骤C：通过中间值K转换得到颗粒微观力学参数，转换方法为：
其中，d是四面体晶格单元的边缘长度，p是正压力；宏观力学参数包括杨氏模量E、泊松比v、抗拉强度Tu、抗压强度Cu和内摩擦系数μi；微观力学参数包括法向刚度Kn、切向刚度Ks、断裂位移Xb、抗剪力Fs和颗粒摩擦系数μp；将转换得到的微观力学参数赋予步骤A建立的模型中的颗粒。
步骤D：对赋值后的堆积模型进行拉伸和压缩测试，获得颗粒宏观力学向量的实际值K1；具体方法为：
根据中间值K对应的杨氏模量E、泊松比v、抗压强度Cu确定对模型箱上方压力板的压力进行单轴压缩试验，通过应力应变曲线得到杨氏模量实际值E1，通过侧向应变和竖向应变的比值得到泊松比实际值v1，通过应力应变曲线峰值确定抗压强度实际值Cu1；
根据中间值K对应的抗拉强度Tu确定加压大小，对模型箱上方的压力板进行单轴拉伸试验，根据应力应变曲线确定模型破坏时的抗拉强度实际值Tu1；
根据求出内摩擦系数的实际值μi1。
步骤E：判断实际值K1是否收敛于设定值K0；引入收敛系数，
σE＝E0/E1
σv＝v0/v1
如果收敛系数在预设区间内，则判断实际值K1收敛于设定值K0，训练过程结束，得到具有接近宏观力学向量的设定值K0的离散元材料。
若收敛系数不在指定区间内，则认为实际值K1不收敛于设置值K0，利用如下规则更新中间值K：
K＝(K0/K1)·K0
具体为：
E＝(E0/E1)·E0
v＝(v0/v1)·v0
Tu＝(Tu0/Tu1)·Tu0
Cu＝(Cu0/Cu1)·Cu0
μi＝(μi0/μi1)·μi0
并重复步骤C～E，更新模型颗粒的微观力学参数，重复获取宏观力学参数的实际值，直到实际值K1收敛于设定值K0。
优选实施例中，当所有宏观力学性质参数的收敛系数同时满足σ∈[0.9,1.1]时，认为实际值K1收敛于设定值K0；本领域技术人员可以根据具体情况，确定收敛系数的取值范围。
上述步骤中的序号A～E仅用于区别各步骤，不用于限定各步骤的顺序。
Claims (10)
1.一种离散元材料自动训练方法，其特征在于：包括以下步骤：
步骤A：建立颗粒随机堆积模型；
步骤B：输入颗粒宏观力学向量的设定值K0，引入宏观力学向量的中间值K，令K＝K0；
步骤C：通过中间值K转换得到颗粒微观力学参数，将微观力学参数赋予模型中的颗粒；
步骤D：对赋值后的堆积模型进行拉伸和压缩测试，获得颗粒宏观力学向量的实际值K1；
步骤E：判断实际值K1是否收敛于设定值K0；若不收敛，利用如下规则更新中间值K：
K＝(K0/K1)·K0
并重复步骤C～E，更新颗粒的微观力学参数，直到实际值K1收敛于设定值K0。
2.根据权利要求1所述的一种离散元材料自动训练方法，其特征在于：步骤A所述的建立颗粒随机堆积模型的方法为根据模型尺寸建立由六块压力板拼合构成的长方体模型箱，模型箱上方的压力板能够沿竖直方向运动，输入颗粒的平均半径和半径范围，在模型箱内随机生成单元，给单元施加随机初速度，对无序运动的单元赋予具有步骤C中微观参数的材料特征使其具有重力和边界范围，赋值后的单元在模型箱内自然堆积，移除模型箱侧面的压力板；平衡模型能量和受力。
3.根据权利要求2所述的一种离散元材料自动训练方法，其特征在于：所述单元为单个颗粒或由多个颗粒构成的clump团簇，单元在模型箱中的排列方式为将模型箱内部划分为规则排列的格点，相邻格点的距离等于颗粒的平均半径，在每个格点位置随机生成一个单元，每个格点处的单元的种类和尺寸随机，单元在模型箱内重力堆积后，对团簇单元进行胶结。
4.根据权利要求2所述的一种离散元材料自动训练方法，其特征在于：所述平衡模型能量和受力的方法为基于牛顿力学方法计算初始时间每个单元受到其他单元和模型箱边界的作用力以及重力的合力；基于合力计算加速度并结合当前速度求得ΔT时间后的初始速度，利用ΔT时间内的平均速度计算ΔT时间内的单元位移，以ΔT时刻的位置和初始速度重复上述过程，从而重复迭代实现离散元法动态模拟，直到模型中的能力和受力达到相对平衡。
5.根据权利要求4所述的一种离散元材料自动训练方法，其特征在于：时间步ΔT满足
其中，T为单元振动周期，f为一个周期内的迭代频率，f＝50。
6.根据权利要求1所述的一种离散元材料自动训练方法，其特征在于：步骤B所述的宏观力学向量的设定值K0为：
K0＝[E0,v0,Tu0,Cu0,μi0]
其中，E0为杨氏模量设定值，v0为泊松比设定值，Tu0为抗拉强度设定值，Cu0为抗压强度设定值，μi0为内摩擦系数设定值。
7.根据权利要求6所述的一种离散元材料自动训练方法，其特征在于：步骤C所述的宏观力学性质向量与微观力学参数的转换方法为：
其中，d是四面体晶格单元的边缘长度，p是正压力；宏观力学参数包括杨氏模量E、泊松比v、抗拉强度Tu、抗压强度Cu和内摩擦系数μi；微观力学参数包括法向刚度Kn、切向刚度Ks、断裂位移Xb、抗剪力Fs和颗粒摩擦系数μp。
8.根据权利要求2所述的一种离散元材料自动训练方法，其特征在于：步骤D所述的获取颗粒宏观力学向量的实际值K1的方法为：
根据宏观力学向量的中间值K对应的杨氏模量E、泊松比v、抗压强度Cu确定对模型箱上方压力板的压力进行单轴压缩试验，通过应力应变曲线得到杨氏模量实际值E1，通过侧向应变和竖向应变的比值得到泊松比实际值v1，通过应力应变曲线峰值确定抗压强度实际值Cu1；
根据中间值K对应的抗拉强度Tu确定加压大小，对模型箱上方的压力板进行单轴拉伸试验，根据应力应变曲线确定模型破坏时的抗拉强度实际值Tu1；
根据求出内摩擦系数的实际值μi1。
9.根据权利要求1所述的一种离散元材料自动训练方法，其特征在于：步骤E中判断宏观力学向量的实际值K1是否收敛于宏观力学向量的设定值K0的方法为：引入收敛系数，
σE＝E0/E1
σv＝v0/v1
如果收敛系数均在预设区间内，则判断实际值K1收敛于设定值K0，训练过程结束，得到具有接近宏观力学向量设定值K0的离散元材料。
10.根据权利要求9所述的一种离散元材料自动训练方法，其特征在于：当所有宏观力学性质参数的收敛系数同时满足σ∈[0.9,1.1]时，认为宏观力学向量的实际值K1收敛于设定值K0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910080559.2A CN109815599B (zh) | 2019-01-28 | 2019-01-28 | 一种离散元材料自动训练方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910080559.2A CN109815599B (zh) | 2019-01-28 | 2019-01-28 | 一种离散元材料自动训练方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109815599A true CN109815599A (zh) | 2019-05-28 |
| CN109815599B CN109815599B (zh) | 2023-03-24 |
Family
ID=66605444
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910080559.2A Active CN109815599B (zh) | 2019-01-28 | 2019-01-28 | 一种离散元材料自动训练方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109815599B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111611695A (zh) * | 2020-05-11 | 2020-09-01 | 太原理工大学 | 一种模拟岩土材料时离散元线性刚度参数的自动标定方法 |
| CN111610091A (zh) * | 2020-05-11 | 2020-09-01 | 太原理工大学 | 一种模拟岩土材料时离散元赫兹接触参数自动标定方法 |
| CN112765895A (zh) * | 2021-01-28 | 2021-05-07 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120253756A1 (en) * | 2011-03-31 | 2012-10-04 | Dem Solutions Ltd. | Method and apparatus for discrete element modeling involving a bulk material |
| CN107391788A (zh) * | 2017-06-09 | 2017-11-24 | 东南大学 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
-
2019
- 2019-01-28 CN CN201910080559.2A patent/CN109815599B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120253756A1 (en) * | 2011-03-31 | 2012-10-04 | Dem Solutions Ltd. | Method and apparatus for discrete element modeling involving a bulk material |
| CN107391788A (zh) * | 2017-06-09 | 2017-11-24 | 东南大学 | 运用三维离散实体解决连续介质构件非线性力学问题的方法 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111611695A (zh) * | 2020-05-11 | 2020-09-01 | 太原理工大学 | 一种模拟岩土材料时离散元线性刚度参数的自动标定方法 |
| CN111610091A (zh) * | 2020-05-11 | 2020-09-01 | 太原理工大学 | 一种模拟岩土材料时离散元赫兹接触参数自动标定方法 |
| CN112765895A (zh) * | 2021-01-28 | 2021-05-07 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
| CN112765895B (zh) * | 2021-01-28 | 2023-10-17 | 南京大学 | 一种基于机器学习的岩土材料离散元自动建模方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109815599B (zh) | 2023-03-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109815599A (zh) | 一种离散元材料自动训练方法 | |
| Horton et al. | A meshless Total Lagrangian explicit dynamics algorithm for surgical simulation | |
| Ding et al. | A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models | |
| CN107644120B (zh) | 一种用于岩土体三维离散元快速建模的通用数值模拟箱的模拟方法 | |
| Stahl et al. | Discrete element simulation of geogrid-stabilised soil | |
| CN103425899B (zh) | 用于页岩气水力压裂的三维离散元建模和模拟方法 | |
| CN107463740A (zh) | 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法 | |
| JP6789274B2 (ja) | Ｂｐｍ理論に基づくメゾの力チェーン粒子モデルの構築方法 | |
| CN112131709A (zh) | 基于近场动力学本构模型的节理岩体力学仿真方法及系统 | |
| Estrada et al. | Yield surfaces and plastic potentials of cemented granular materials from discrete element simulations | |
| Place et al. | Simulation of the micro-physics of rocks using LSMearth | |
| Nguyen et al. | Discrete modeling of strain accumulation in granular soils under low amplitude cyclic loading | |
| Dabeet et al. | Discrete element modeling of direct simple shear response of granular soils and model validation using laboratory element tests | |
| CN110929412B (zh) | 一种基于dda理论的节理摩擦系数动态衰减计算方法 | |
| Masson et al. | Multiscale simulations of the mechanical behaviour of an ensiled granular material | |
| Xiang et al. | Modelling breakwater armour layers and the dynamic response of armour units | |
| Baudet et al. | Integrating tensile parameters in 3D mass-spring system | |
| Vermeulen et al. | Soil dynamics of the origination of soil tare during sugar beet lifting | |
| Faria | Seismic behaviour of concrete dams and its interaction with the rock foundation | |
| Khelalfa et al. | Deciphering Soil-Plant and Soil-Insect Interactions: A Mathematical Modeling Approach | |
| Duan et al. | Prediction of crack propagation by vector form intrinsic finite element method | |
| Shi et al. | Fitting method research of sieving parameters of vibration screen and sieving efficiency based on dem | |
| Wang et al. | Influence of particle gradation curve on granular packing characteristics | |
| Smolin et al. | 3D Modeling of the mechanical behavior of ceramics with pores of different size | |
| CN116956406A (zh) | 滑坡动态破裂运动全过程sph分析方法与系统 |
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 |