首页/文章/ 详情

固体断裂(solid fracture)的前沿研究

8小时前浏览5

点击上方蓝字了解更多计算与STEM领域研究前沿



文一:

 

裂缝交叉口角度模式对非线性流动行为和局部流场的影响:可视化实验

摘要:

相贯裂隙角度形态的影响对于深入理解裂隙岩体的流动特性至关重要。采用水力学和可视化实验相结合的方法,研究了“一进两出”不同角度裂隙交叉口的宏观非线性流动行为和局部流场特性。角度模式对交叉口的形态有重要影响,导致交叉口非线性流增强,流量分布发生改变。无论是非线性系数(b)还是流量分布比率 ,都表现出与角度模式相对应的正态分布模式。通过观察不同角度条件下裂缝的交叉形态和局部流场,发现进出口裂缝引起的突起对进入裂缝的流体流动有显著影响,特别是在高雷诺数(Re)条件下。采用一种基于神经网络的机器学习算法,预测了不同 Re 值和角度方向图下交叉口裂缝出口处的流量分布。该模型对不同条件下交叉裂缝出口流量分布进行了准确预测。此外,通过分析交叉口内不同位置的流速和流场特性,阐明了交叉口内宏观流速分布与流场特性之间的关系。这些结果将有助于研究地下水流动行为、微流控芯片中的流体流动和渗流传热。

 

图:岩体中的天然裂隙和裂隙交叉。

 

图:实验系统和断裂模型示意图:a断裂流可视化实验系统,b相交断裂模型,c具有断裂相交的PMMA模型,d角度分布模式。

 

图:实验误差与图像处理过程:a流速误差,b压力梯度误差,c速度误差,d图像处理过程。

 

图:不同角度模式下的局部焦场特征。

文二:

 

慢三维断裂中的拓扑缺陷形成

摘要:

裂纹在三维(3D)材料中传播时会形成各种表面图案。名义拉伸(I型)断裂中的对称破坏拓扑缺陷以表面台阶的形式出现在缓慢(非惯性)区域。我们发现,最近为三维动态(惯性)拉伸断裂提供基本线索的相同相场框架也会产生裂纹表面台阶。台阶形成被证明涉及两个基本的物理成分:有限强度淬火无序和一个小的介观反平面剪切(模式III)加载分量(在主要的拉伸模式I加载分量之上)。我们量化了无序(包括其强度和空间相关长度)与介观I+III模式混合在控制台阶形成中的非线性相互作用。最后,我们表明,表面台阶是从小尺度的背景表面粗糙度中生长出来的,由两个重叠的裂纹段组成,由桥接裂纹连接,这与实验结果一致。

 

图:长杆中裂纹产生的断裂面(绿色)。

文三:

 

非线性断裂力学模型:从极限全寿命模型到循环载荷下基于亥姆霍兹自由能分离的石灰岩断裂演化

关键词:

岩石结构和荷载条件在岩石的变形和破坏过程中起着相当大的作用。为了探索脆性岩石的非线性力学机制,本研究提出了不等幅循环载荷下石灰岩试样的极限全寿命本构方程、断裂方程和岩石阻力断裂能(RBE)方程。该理论模型基于微积分方法和热力学原理,对复合系统作用后岩石系统的亥姆霍兹自由能进行非线性分离。基于非线性分离参数和断裂演化的解析解,推导出了石灰岩的极端全寿命本构方程,而基于残余断裂能和热力学原理推导出了断裂方程。我们给出了材料RBE的解析解,它反映了岩石在负载下抵抗断裂的能力。通过在不等幅循环荷载下具有不同预制裂缝的石灰岩样品验证了模型性能。本研究可为不同条件下非线性断裂模型的研究奠定基础。

 

图:复合热力学系统中的脆性岩石破裂:(a)复合热力学系统:该系统包括储热器、加载装置和材料体;(b) 实验曲线和本构目标曲线;以及(c)石灰石样品渐进破碎过程中的能量演化。

 

图:加载下的客观理论曲线和实验数据:(a)石灰石在循环加载下的实验数据;(b)循环加载应力极限状态下的实验数据和理论曲线。

 

图:不同试样的实验数据:(a)弹性能与应变之间的关系;(b)断裂与应变之间的关系。

 

图:石灰岩样品的微观结构分析:(a)石灰岩样品的扫描电镜(SEM)图像;(b、c和f)局部放大图像;(g) 矿物成分分析结果;以及(h)不同数量和位置的预制裂缝样品及其相应的破坏状态。

 

图:大理石试样的理论模型验证:(a和c)大理石试样D和白色大理石试样在循环载荷下的实验应力-应变曲线;以及(b和d)理论曲线。

文四:

 

低温气冰压力下非饱和岩石裂纹的弹塑性裂纹萌生机理

摘要:

在分析寒冷地区岩体的冻胀和开裂时,通常假设裂缝完全饱和。然而,在实际工程实践中,岩体裂缝往往处于非饱和状态。在水分冻结时,冰的体积增加,气体的体积减少,气体压力增加,而裂纹壁受到气体压力、冻结压力和冰的摩擦力的影响,导致复杂的应力条件。本文研究了气体-冰压力影响下非饱和岩石裂纹萌生的机理。假设屈服范围较小,根据复变函数和弹塑性裂纹力学理论,推导出了冻结后气体压力、塑性区、应力强度因子、裂纹萌生角和裂纹萌生应力的计算公式。对冻结温度Ti、充水体积比sw和裂纹形状等关键参数进行了分析和讨论。结果表明,只有当sw>0.7时,冷冻后才会产生有效气压(>1atm);Ti越低,冷冻后的气体压力越小。此外,裂纹形状越倾向于圆形,冷冻后就越容易产生有效的气体压力,而Ti越低,裂纹形状就越倾向于椭圆形。此外,Ti越小,裂纹起始角越大,而sw和裂纹形状的变化对裂纹起始角几乎没有影响。

 

图:非饱和裂缝的应力模型。

 

图:不同冷冻温度下冷冻压力的变化。

 

图:计算模型:(a)模型几何;以及(b)断裂起始的描述。

 

图:交替迭代过程。

文五:

 

埋地爆炸引起的土壤破碎和坑化的耦合并行近场动力学-SPH建模与模拟

摘要:

在这项工作中,我们开发了一个耦合的近场动力学(PD)和平滑粒子流体动力学(SPH)模型,用于模拟地下爆炸下的土壤破碎、喷射和火山口形成。具体来说,基于非寻常状态的PD理论描述了土壤的动态响应,而SPH方法则控制着爆炸性气体产物的运动。采用了一种稳定高效的耦合算法,用于PD和SPH之间的数据传输。为了准确捕捉爆炸荷载下的土壤行为,我们提出了一种改进的Drucker-Prager塑性模型,其中加入了额外的高压压实状态方程(EoS)。此外,还为耦合的PD-SPH模型开发了一个使用OpenMP算法的并行代码,以处理大规模、长时间的土壤爆炸。PD-SPH预测与两组精心设计的离心模型试验进行了比较。模拟中考虑了重力效应,以复 制土壤中埋藏爆炸的真实物理过程。研究发现,耦合的PD-SPH模型成功地捕捉到了埋藏爆炸引起的土壤破碎、喷射和火山口形成的整个物理过程。预测的土壤喷射和火山口形态,以及喷射速度、喷射高度和火山口直径等定量结果与试验结果一致。此外,研究结果表明,重力效应对土壤喷射和坑化过程有显著影响。随着重力加速度的增加,土壤喷射高度和火山口尺寸显著减小。

 

图:PD和SPH的直接耦合方案。

 

图:NNPS中条带PIB算法的模式化。


来源:STEM与计算机方法
ACT断裂非线性芯片裂纹理论自动驾驶爆炸材料数字孪生控制试验人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-12-04
最近编辑:8小时前
江野
博士 等春风得意,等时间嘉许。
获赞 48粉丝 50文章 315课程 0
点赞
收藏
作者推荐

颗粒力学(particle mechanics)的前沿研究分享

文一: 欧拉和拉格朗日方法在运动列车产生的微粒浓度数值研究中的比较摘要:最近的研究表明,由于颗粒物的连续产生和积累,地铁系统内的空气质量较差。其中许多颗粒是由刹车片磨损产生的含铁空气磨损颗粒。本研究调查了列车制动过程中列车制动片排放的磨损颗粒的浓度分布。由于列车在地铁系统的大几何形状内运动引起的非定常三维湍流,这种数值模拟需要很高的计算成本,特别是在使用跟踪每个粒子的拉格朗日方法时。因此,采用将颗粒相视为连续体的成本较低的欧拉方法更为可取。然而,应研究使用欧拉方法代替拉格朗日方法的可行性,特别是对于列车诱导湍流流场内运动物体发射的大颗粒。因此,在本工作中,对通过这两种方法获得的颗粒浓度预测进行了比较。欧拉方法的预测更加连续和均匀,而基于单个粒子累积的拉格朗日方法的预测则更加分散。结果表明,欧拉方法可以在保持计算效率的同时提供合理的预测。 图:a隧道和车站计算域。b列车(初始位置)和隧道模型示意图。c显示计算域横截面几何形状和隧道沿线网格界面的示意图。 图:a列车模型的速度。b列车乘客车厢示意图,重点是刹车片位置。 图:位于左右平台人类呼吸区的表面示意图。比较左表面b右表面归一化平均颗粒浓度时变的欧拉和拉格朗日方法的结果。 图:比较欧拉和拉格朗日方法对车站隧道连接表面归一化平均颗粒浓度时变的结果:a隧道入口1和b隧道入口2。 图:比较通过欧拉和拉格朗日方法获得的人类呼吸表面归一化粒子浓度轮廓的结果。 图:比较通过欧拉和拉格朗日方法获得的距离站入口74米处人类呼吸表面b a表面t=102秒归一化粒子浓度轮廓的结果。文二: 颗粒性质异质性对砂土材料力学性质的影响摘要:颗粒组成的随机排列或颗粒固有性质的异质性是导致砂土力学参数变化的原因。将离散元法与随机场理论相结合,改进了一种随机离散耦合建模方法,其中使用滚动阻力接触模型来考虑砂土颗粒的滚动摩擦特性。考虑到两个随机场因素,即变异系数和波动范围,为典型工作情况开发了四个随机离散元模型(RDEM)。采用蒙特卡洛方法模拟了三轴RDEM,分析了应力-应变关系、强度和变形参数。研究了变异系数(Cov)和波动范围(δ)对力学特性变异性的影响。结果表明,随着Cov的增加,随机非均质样品的平均剪切模量G、峰值摩擦角φp和临界摩擦角φc下降,而平均泊松比μ增加。随着δ的增加,随机非均质样品的变形和强度参数的平均值最终接近均质样品的平均值。 图:滚动阻力模型图。 图:随机离散耦合过程原理。 图:三轴试验模型。 图:不同Cov和δ粒子的典型随机场模型。 图:非均质土样直剪试验结果比较。 图:代表性均匀和非均匀试样的平均法向接触力。文三: 颗粒形状对颗粒材料剪切行为影响的DEM研究摘要:本研究使用二维离散元法模拟研究了颗粒凸度、球形度和纵横比(AR)对剪切颗粒材料行为的影响。制备了具有不同颗粒形状的各向同性、致密和松散的组件,并通过双轴压缩进行排水剪切,直至达到临界状态。介绍了强度和剪胀性等宏观特征。然后通过考虑配位数、织物各向异性、颗粒力矩、接触处的摩擦运动和颗粒旋转来研究宏观行为背后的因素。对于这里考虑的形状范围,数据表明,剪切强度随着颗粒凸度和球度的增加而降低,而剪切强度随着AR的增加而增加。剪切强度和凸度之间的相关性较弱,但AR和强度之间的相关性更强。随着AR的增加,大应变下的体积应变趋于增加。临界状态强度与临界状态配位数和临界状态机械孔隙比之间的相关性比临界状态孔隙比与临界状态强度之间的相关性更强。接触织物各向异性、颗粒传递的力矩大小和接触处产生的摩擦力是影响强度的重要因素。临界状态强度随着平均颗粒力矩和平均移动摩擦的增加而增加。颗粒旋转分析为颗粒材料对剪切的响应提供了见解。 图:DEM研究中使用的粒子类型及其球度S值。 图:双轴剪切下的样品图示。 图:颗粒形状对堆积的影响:a凸度b球度c纵横比(“2D-DEM”表示当前的研究。) 图:三维团簇致密样品的粒子取向彩色图。文四: 通过光滑粒子流体力学(SPH)模拟和实验研究研究横向速度对Al-16Si合金摩擦堆焊的影响摘要:移动速度对Al-16Si合金摩擦片显微组织、力学性能及耐磨性能的影响 使用平滑粒子流体动力学 (SPH) 模拟和实验评估 AA1050 合金表面 结果显示,随着移动速度的增加,涂层的高度和宽度分别减少了 54% 和 20% 速度从 75 毫米/分钟增加到 115 毫米/分钟。此外,观察到界面处的未粘合区域相应增加 模拟结果显示了以横向速度涂覆的样品在涂层/基材界面处的最大剪切应力 75、95 和 115 mm/min 分别为 83、95 和 112 MPa。随着移动速度从 75 增加到 115 mm/min, 摩擦堆焊所需的预测扭矩和垂直力分别增加了92%和22%。当移动速度从 75 至 115 毫米/分钟。移动速度从 75 至 95 毫米/分钟增加导致晶粒尺寸减小 11% 和 13%, 与 AA1050 基材相比,硬度、强度和耐磨性分别提高 12% 和 8%。 图:显微组织检查、SPT试验、磨损试验和试样提取位置推压的示意图。突出显示的区域显示了检查的位置。 图:不同样品的横截面;a样品C-75,b样品C-95,C样品C-115。 图:不同样品涂层过程中的预测扭矩和垂直力。 图:A390消耗棒和样品C-75不同区域的光学显微镜图像。 图:不同放大倍数下A390自耗棒的SEM图像。 图:磨损表面的SEM图像;a AA1050基板,b A390耗材棒,c样品c-75,d样品c-95,e样品c-115。文五: 基于SPH的流体和结构力学修正输运速度公式摘要:为了提高光滑粒子流体力学(SPH)方法的计算精度,采用了粒子移动技术(PST)。移动确保了粒子在空间中的均匀分布。这可以通过使用输运速度移动粒子来实现。在本文中,我们提出了一类扩展的输运速度公式(TVF)方法。我们以一致的方式导出了方程,并且证明了还有一些额外的项显著地提高了方法的精度。特别地,我们将其应用于熵阻尼人工压缩 SPH 方法。我们用一种简单的方法来识别自由表面的粒子及其法线,从而使这种方法适用于自由表面问题。我们展示了新方法如何应用于弹性动力学问题。我们考虑了一组基准问题,包括流体力学和固体力学,以证明该方法的准确性和适用性。实现是开源的,手稿是完全可重复的。 图:流体圆环自由表面粒子的识别。描绘流体粒子的法线,b流体粒子的边界粒子。 图:溃坝模拟中流体法线的识别。向我们显示所有流体粒子的法线。 图:雷诺数为1000的泰勒-格林涡旋的粒子图,分辨率为150×150。颜色代表压力。 图:Re=100的空腔的粒子图,粒子排列为150×150,左侧有校正项,右侧没有校正项。来源:STEM与计算机方法

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈