文一:
基于真实集料形态的离散元模型研究沥青混合料的开裂行为
摘要:
本研究旨在通过实验室试验和离散元模型研究沥青混合料的抗裂性。针对虚拟裂缝模型的构建,提出了基于各种软件的交互式二维建模方法和三维骨料的建立方法。在建立二维模型时,考虑到AC-20沥青混合料的细观结构特征,基于PFC2D软件和图像处理,建立了集料-沥青砂浆空隙三相结构。考虑到粗骨料的级配和形态,利用PFC3D建立了包含真实骨料形状的三维模型。基于该数值模型,对不同长度的预切狭缝试样进行了半圆形弯曲和拉伸试验,分析了低温条件下试样的应力水平、开裂过程和裂纹位置。以断裂性能参数为评价指标,通过与室内试验数据的比较,验证了数值模拟的可行性。在模型加载的早期阶段,应力分布具有压缩-拉伸区的特征。裂缝的形成经历了一个平稳膨胀和快速渗透的过程,压缩区的裂缝穿过骨料。粒子流程序给出的失效状态与试验获得的裂纹演化结果基本一致。此外,模拟的断裂参数与实验观察结果有很好的重叠,表明离散元模拟可以充分反映沥青混合料中的裂纹扩展过程。该方法对于深入有效地探索沥青混合料的开裂性能是有效的。
图:2D虚拟模型生成过程。
图:三维虚拟模型生成过程。
图:不同缺口长度试样的三维模拟失效结果。
图:(a) 具有不同缺口长度的试样中每个时间步长的断裂数量,以及(b)裂纹发展过程的3D模型。
文二:
基于 GPU 的实际非凸集料离散元模型:ITZ 体积分数与混凝土扩散系数的细观研究
摘要:
具有真实非凸形态的骨料及其周围的界面过渡区(ITZ)显著影响混凝土的扩散率。在此,我们最初提出了一种数学可控的参数化方法来设计三维(3D)非凸聚集体。然后,我们开发了一个基于GPU的离散元建模(DEM)来生成多分散非凸聚集体的随机紧密堆积。接下来,我们建立了包含非凸集料和ITZ的三相混凝土细观结构,并通过蒙特卡罗随机点采样方法获得ITZ体积分数。最后,我们提出了一个双概率布朗运动(DP-BM)方案来获得混凝土的有效扩散率。此外,还研究了骨料的体积分数、粒径多分散性和形状以及ITZ厚度对ITZ体积分数和混凝土扩散率的影响,以探索混凝土的严格组分-结构-性能关系。这项工作为混凝土的细观结构表征和扩散率评估提供了一个全面的数值框架,可作为混凝土耐久性评估的宝贵工具。
图:收缩范围 θr 对骨料形态的影响: Nc = 50,α = -0.2。
图:两个非凸集料的接触检测示意图:(a)主粒子和从粒子; (b)基于主粒子的水平集函数 φm 和从粒子的离散表面点的接触检测。
图:混凝土三相细观结构的可视化:(a)尺寸为h=0.05Deq的ITZ层周围的单个骨料及其横截面;(b) 三相细观结构和一个具有遵循富勒分布的聚集体的横截面,Deqmin=10mm,Deqmax=40mm,fp=0.6,h=0.05Deqmin。
文三:
道床颗粒离散元建模的标定方法
摘要:
离散元法(DEM)被广泛用于模拟静态和动态负载条件下的压载行为。由于道碴颗粒的形状不规则,在校准道碴颗粒DEM输入参数时,应考虑多个参数(即静摩擦系数和滚动摩擦系数、杨氏模量、泊松比和恢复系数)的同时影响,这是一个挑战。为了减少形状因素对压载样品体积行为的影响,首先创建并生成了一系列用于描述由17-24个球组成的颗粒形状的团块。在下文中,提出了一种有效的校准方法,包括四种不同的试验,包括空心圆柱体试验、承压试验、直剪试验和压载箱试验。由于摩擦系数的变化对所有测试结果的影响,执行测试的顺序很重要。由于道砟样品上没有外力,因此休止角试验受接触和阻尼参数的影响不大,因此,首先,通过几个休止角试验来校准颗粒摩擦系数。然后,使用直接剪切试验来校准颗粒的杨氏模量和泊松比。最后,利用循环载荷下压载箱试验的DEM模型对颗粒的恢复系数进行了标定。所提出的方法将有助于校准压载颗粒特性,这是模拟任何实验之前的必要阶段。
图:a粗糙表面的赫兹接触,b DEM建模的赫兹-明德林接触模型
图:a直接剪切试验装置,b剪切试验的DEM模拟,c,d初始和最终状态下的法向/切向力分布,e XZ平面模拟剪切试验的法向力玫瑰图
图:a休止角试验示意图,b道砟休止角试验的DEM模拟
文四:
烟叶拉伸行为的离散元模拟
摘要:
烟叶的力学性能对于烟叶收获机械的研究和设计具有重要意义。采用离散单元法(DEM)建立了烟叶试样拉伸行为的数值模型。模型由粘结在一起的球形颗粒构成。模型输出为叶片的抗拉强度(σ 宏)和杨氏模量(Emacro)。通过拉伸试验测定了不同取样方向和不同取样位置的烟草叶肉样品的拉伸性能。通过仿真结果与试验数据的比较,对模型进行了标定和验证。试验结果表明,取样位置和取样方向对烟叶的拉伸性能影响不大。平均拉伸强度为0.57 MPa,平均杨氏模量为3.50 MPa。用两个方程很好地描述了模型输出(σ 宏观和 Emacro)与两个最敏感的模型参数(粒子弹性模量(Emicro)和键的临界拉应力(σ 微观)之间的关系。利用拉伸强度和杨氏模量的方程和实测值,标定出 σ 值为8.02 MPa,E 值为23.07 MPa。校正后的模型能够预测烟叶的抗拉强度和杨氏模量,相对误差小于2% 。所提出的 DEM 模型和方法可以应用于其他植物叶片的拉伸行为和性质的模拟。
图:模型叶肉样品。
图:在微观水平上的模型断裂行为: (a)粘结力; (b)模型颗粒分布(左)与相关联的接触力链(右)。
图:模型微观参数和宏观参数之间的关系(a)颗粒的杨氏模量(Emicro)和杨氏模量(Emacro)之间的关系;(b) 颗粒的杨氏模量(Emicro)、键的临界拉伸应力(σmicro)和拉伸强度(σmacro)之间的关系。
文五:
X射线计算机断层扫描表征的锂离子阴极结构的离散元法和电化学建模
摘要:
电极微观结构对锂离子电池的性能有着深刻的影响。在这项工作中,使用具有键合颗粒模型的离散单元法(DEM)研究了压延工艺对电极微观结构的影响。对使用X射线计算机断层扫描(XCT)表征的真实电极结构和理想DEM结构进行了综合评估。研究了层析成像扫描和DEM结构的电极结构和传输特性,即孔隙率分布、比表面积和扭曲因子。在考虑碳结合域(CBD)相之后,进一步进行电化学分析。考虑到压延的影响,层析成像和DEM模拟的理想化结构之间实现了极好的一致性。采用电极压缩,压延后电池性能得到改善。本研究为今后利用DEM和电化学分析对电池性能进行定量评价提供了依据。
图:AM颗粒相的图像处理过程。
图:DEM模拟设置(a)颗粒间结合(b)压缩试验(c)AM粒度分布(d)DEM模拟的孔隙压力关系
图:(a)锂化状态(SOL)和过电位等值线的比较(b)在25%放电度下每个颗粒中心点的SOL和(c)在不同放电速率下基于断层扫描和理想化微观结构的放电曲线。
图:DEM生成的AM粒子结构的结构分析(a)沿z方向的最小路径长度图(b)通过欧几里得距离变换计算的粒子间间隙(c)通过Navier-Strokes方程计算的流场。