首页/文章/ 详情

固体大变形+流固耦合方法——物质点法(MPM)

9月前浏览756


文一:

 

广义插值物质点法

摘要

使用变分形式和Petrov–Galerkin离散化格式推广了计算固体力学的材料点法(MPM)离散解程序,产生了一系列称为广义插值材料点(GIMP)方法的方法。该泛化允许识别其他基于点或节点的离散求解技术的各个方面,这些技术不使用固定网格,即“无网格方法”。注意到了相似之处,并确定了与其中一些方法相比的一些实际优势。示例用于演示和解释MPM计算中可能出现的数值伪影噪声。这种噪声导致材料点的非物理局部变化,从而评估本构响应。在一种情况下,它会破坏显式解决方案,而在另一种情况中,它会严重降级。由于噪声输入,依赖历史的非弹性本构定律可能会错误地发展,并报告不准确的应力状态。噪声是由于插值函数缺乏平滑性,并且是由于材料点跨越计算网格边界而产生的。GIMP方法中可用的下一个平滑度被证明能够消除小区交叉噪声。

 

图:由缓慢增加的体力b引起的连续棒准静态压缩的数值和分析解。MPM溶液材料点应力用开三角形表示,连续颗粒GIMP溶液材料点应由力用开圆表示。垂直线表示单元边界。实心对角线表示解析解。在增加柱重量的情况下获得解决方案。

 

图:针对不同粒子密度的MPM、GIMP和有限GIMP算法的收敛性研究结果。相对于精确解的误差是根据计算单元大小绘制的。计算的是柱在重力作用下的准静态压缩

文二:

 

颗粒材料的物质点法

摘要:

提出了一个颗粒材料模型,描述了每个颗粒的内部变形和颗粒之间的相互作用。该模型基于FLIP材料点、单元中粒子法,求解每个晶粒的连续本构模型。晶粒之间的相互作用是用接触算法计算的,该算法禁止相互渗透,但允许分离和摩擦滑动和滚动。细胞内颗粒法消除了单独接触检测步骤的需要。在接触模型中使用公共休息框架产生了计算成本与晶粒数量的线性比例。该模型的特性通过滑动和滚动接触的数值解来说明,颗粒材料的特性通过剪切计算来说明。数值计算结果表明,对于形状由计算网格解析的光滑颗粒,可以准确地对接触进行建模。

 

图:用于模拟倾斜平面上圆盘的几何图形。

 

图:可变形圆盘的质心位置(x分量)作为时间的函数(虚线),用于摩擦系数p=0.3和0.9以及倾角0=rr/3的模拟。为了进行比较,还显示了刚性圆盘(实线)的解析解。在计算解和分析解中,p=0.3的曲线高于p=0.9的曲线。

文三:

 

流固耦合问题的物质点法

摘要:

本文提出了一种新的材料点法,用于求解流体饱和土在大变形作用下的流体力学耦合问题。提出了一种基于Biot混合理论的土壤-孔隙流体耦合MPM算法,用于解决包括地下水位位置随时间变化在内的流体力学相互作用问题。通过将一维固结试验的模拟结果与相应的解析解进行比较,检验了所提出方法的准确性。

对所提出的方法中使用的MPM参数进行了灵敏度分析,以检查每网格粒子数和网格尺寸对求解精度的影响。为了证明该方法的能力,对一个大型堤防渗流破坏的物理模型实验进行了模拟。地下水位随时间变化的堤坝模型的行为与实验观测结果非常吻合。通过检查孔隙水压力以及模拟计算的有效应力,讨论了渗流引起的破坏机制。

 

图:典型计算网格和粒子示意图,以及分析流程:(a)典型计算网格与粒子示意图;(b) 分析流程示意图

 

图:提出方法的操作流程

 

图:MPM模型和固结模拟结果:(a) MPM模型;(b) 孔隙水压力分布;(c) 固结度与时间因素

 

图:(a)颗粒的位置和(b)堤防中的偏剪切应变(膨胀角5 0.0)

文四:

 

考虑有限变形的岩土材料脆韧性破坏转变的相场隐式粒子域插值材料点法

摘要

提出了一种基于对流粒子域插值的相场隐式材料点法(PF-ICPDI),用于建立压敏土工材料在有限变形条件下的脆韧性破坏转变模型。在该方法中,采用基于微力平衡定律和热力学第二定律的相场断裂公式作为非局部损伤函数进行弹塑性断裂分析。利用光滑的双屈服面塑性本构模型来评估压敏土工材料的脆韧性破坏转变行为。通过引入有效相场应力并将总储存能量分解为弹性和塑性部分,建立了相场断裂模型和帽塑性模型的耦合效应。针对准静态弹塑性断裂有限变形问题,提出了隐式材料点法,该方法可以避免网格严重变形,提高数值稳定性。

此外,在模拟大变形脆塑性失效过渡过程时,采用了对流粒子域插值技术来消除数值噪声,提高材料点跨越单元边界时的计算精度。采用交错增量迭代法求解耦合离散控制方程。

通过几个具有代表性的数值例子证明了所提出的PF-ICPDI方法的准确性和能力,并讨论了主要材料参数对岩土材料相场脆性-韧性破坏转变模型的影响。

 

图:在弹塑性断裂建模中,通过选择相应的返回映射算法来更新本构关系的过程。

 

图:不同围压下平面应变压缩试验的等效塑性应变(a–g)和相场(h–n)的等值线:(a,h)σc=5 MPa,(b,i)20 MPa,(c,j)40 MPa,(d,k)60 MPa,(e,l)100 MPa,(f,m)σc=150 MPa和(g,n)200 MPa。

文五:

 

冲击和穿透模拟中材料点法的接触算法

标准材料点法(MPM)中固有的无滑动接触约束产生了更大的穿透阻力。因此,标准MPM不能很好地处理涉及冲击和穿透的问题。为了克服这些不足,提出了两种MPM的接触方法,并在我们的三维显式MPM代码MPM3D中实现。在MPM中,在每个时间步长开始时,重新定义的规则网格上可能不满足不可穿透性条件,即使在最后一个时间步长结束时已将其施加在变形网格上。物体之间的不可穿透性条件在第一接触方法中仅施加在变形网格上,而在第二接触方法中,它同时施加在变形格栅和重新定义的规则格栅上。此外,提出了三种方法用于冲击和穿透模拟,以确定满足接触表面共线条件的表面法向量。通过模拟两个弹性环的碰撞和球体滚动问题验证了接触算法,并将其应用于粒子失效模型模拟钢球穿透和厚板穿孔。在弹性环碰撞的模拟中,第一接触算法在总能量中引入了显著的扰动,但第二接触算法可以通过使用更大的时间步长来获得稳定的解。这两种接触算法似乎都能很好地解决其他问题,如球体滚动和弹丸穿透

 

图:MPM中的材料离散化

 

图:通过第一接触方法(=0.25)获得的弹性环在不同时间步长的配置:(a)t=0ms;(b) t=0.2毫秒;(c) t=0.4毫秒;(d) t=0.8毫秒;(e) t=1.8ms;

(f) t=2.8毫秒;(g) t=3.8ms;以及(h)t=4.8ms。

 

图:在V0=575m/s的打击速度下的弹-目标相互作用。Vr表示射弹的剩余速度:(a)实验照片,Vr=455m/s;(b) 模拟标准MPM,Vr=334.1m/s;(c) 模拟第一次接触MPM,Vr=442.8m/s;以及(d)第二接触MPM的模拟,Vr=453.4m/s。


来源:STEM与计算机方法
断裂碰撞岩土求解技术理论材料数字孪生控制试验人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:9月前
江野
博士 等春风得意,等时间嘉许。
获赞 58粉丝 96文章 153课程 0
点赞
收藏
作者推荐

数据驱动框架下的材料科学与工程(material science and engineering)

文一: 由自然语言处理和信息抽取驱动的数据驱动的材料研究摘要:随着数据科学和机器学习在社会各个方面的出现,尤其是在科学领域,人们越来越重视获取数据。基于所探索的材料类别的显著范围和感兴趣的材料特性的多样性,材料科学中的数据尤其具有异质性。这导致了范围在许多数量级的数据,这些数据可能表现为数字文本或基于图像的信息,这需要定量解释。因此,通过改编自自然语言处理领域的技术,跨领域自动消费和编纂科学文献的能力具有巨大的潜力,可以解锁和生成数据科学和机器学习所需的丰富数据集。这篇综述侧重于材料科学文献的自然语言处理和文本挖掘的进展和实践,并强调了在文章的图表中提取文本之外的额外信息的机会。我们讨论并举例说明了对材料进行自然语言处理的几个原因,包括数据汇编、假设发展以及理解领域内和跨领域的趋势。详细介绍了当前和新兴的自然语言处理方法及其在材料科学中的应用。然后,我们讨论材料科学领域的自然语言处理和数据挑战,未来的方向可能会被证明是有价值的。 图:高能物理、天文学和天体物理学中大型集中式数据集与材料物理学中异构、分散数据的比较。与其他领域不同,材料科学缺乏足够的动力来实现数据集中,这不仅是因为数据如此多样化,还因为数据来自各种独立科学家和实验室来源。由于科学家研究的材料种类繁多,材料科学领域的数据尤其异质。数据显示为数字文本或基于图像的信息,这需要定量解释。 图:自然语言处理示意图,包括每个步骤的工具和模型示例。很明显,大多数自然语言提取遵循类似的方法: (1)获取相关文本资源,(2)将文本处理成单独的术语(也称为标记化) ,(3)文档分割和段落分类,(4)将标记识别为信息类,(5)实体关系提取,(6)命名实体链接 图:图像提取示意图,包括来自显微镜图像或分子结构的例子。科学文献中的数据来源于基于图像的计量学,主要来源于(扫描)透射电子显微镜和原子力显微镜。尽管这些图像包含了关于材料结构的大量定量数据,但是大多数图像在其周围的文本中都进行了定性的讨论。这些数据可以揭示特别重要的依赖于纳米技术或晶体学的研究问题。该图提供了从文本中提取此信息的路径。文二: 一种用于数据驱动材料发现的设备流水线设计摘要:世界需要新材料来刺 激我们经济的关键部门的化学工业:环境和可持续性、信息存储、光通信和催化。然而,几乎所有的功能材料仍然是通过“试错”发现的,其中缺乏可预测性是技术创新的主要材料瓶颈。目前,材料发现的平均“分子到市场”交付时间为20年。正如材料基因组计划所强调的那样,这对工业需求来说太长了,该计划的宏伟目标是将平均分子到市场的交付周期缩短4倍。只有采用一种全新的材料发现方法,才能现实地实现进展中如此巨大的变化。幸运的是,一种全新的材料发现方法 正在出现,通过这种方法,人工智能的数据科学提供了一种有前景的解决方案,可以加快这些平均分子到市场的交付周期。这种方法被称为数据驱动的材料发现。鉴于“大数据”、人工智能和高性能计算(HPC)的及时和重大进展,其广阔的前景直到最近才成为现实。政府监管的数据和文献开放访问要求刺 激了对海量数据集的访问。可以挖掘数据并在其中找到模式的自然语言处理(NLP)和机器学习(ML)工具正在成为主流。Exascale HPC功能可以帮助数据挖掘和模式识别,也可以通过计算生成自己的数据,现在我们已经掌握了。这些及时的进展为开发数据驱动的材料发现策略提供了一个理想的机会,以系统地设计和预测特定设备应用的新化学品。该账户展示了数据科学如何通过四步“设计到设备”的管道来提供材料发现,该管道包括(1)数据提取、(2)数据富集、(3)材料预测和(4)实验验证。同源化学和性质信息的大规模数据库首先由ChemDataExtractor等“化学感知”自然语言处理工具伪造,并使用机器学习方法和高通量量子化学计算进行丰富。然后,可以通过挖掘这些数据库来预测定制应用的新材料,这些数据库具有已知可提供功能材料的化学结构和物理性质之间关系的算法编码。这些可以采取分类、枚举或机器学习算法的形式。数据挖掘工作流程将这些预测简短地列出给少数几个主要的候选材料,然后进行实验验证。正在开发这种从设计到设备的方法,为加速发现功能应用的新化学品提供路线图。所提供的案例研究证明了其在光伏、光学和催化应用方面的实用性。虽然该账户侧重于物理科学的应用,但所讨论的通用管道很容易转移到生物学和医学等其他科学学科。 图:数据驱动材料发现的四步设计到设备方法的总体工作流程 图:(a) 化学和性质空间及其结构-功能关系的概念说明,为挖掘这些数据源提供了知识库。(b) 通用数据挖掘工作流程,将化学数据过滤到连续的短列表中,在那里它们遵循适合给定材料应用的结构-功能关系,直到只剩下少量材料;它们服从所有关系,因此代表领先的候选材料预测。 图:化学感知的自然语言处理文本挖掘工具 ChemData萃取器的操作工作流程 图:使用人工神经网络预测优化催化过程的机器学习算法示例。文三: 基于数据驱动的原子模拟材料模型摘要:经典分子动力学模拟材料的中心近似是用来计算原子受力的原子间势。用传统的方法构造可行的函数形式和找到准确的参数化势,需要付出巨大的努力和创造性。机器学习已经成为开发准确而强大的原子间势能的一种有效的替代方法。从一个非常普遍的模型形式开始,势能直接从电子结构计算的数据库中学习,因此可以被视为量子模拟和经典原子模拟之间的多尺度联系。不准确的外推风险存在于狭窄的时间和长度范围之外,这两种方法可以直接进行比较。在这项工作中,我们使用频谱邻域分析势(SNAP) ,并显示了如何可以产生最小的插值误差拟合,这也是鲁棒的外推超越训练。为了证明这种方法,我们已经开发了一个适用于所有二元化合物的钨铍电位。随后,在固体钨的(001)表面进行了大规模的分子动力学模拟。机器学习 W-Be 电位产生一个与缺陷形成能量的量子计算一致的注入结构的种群。预测了很浅(< 2nm)的 Be 平均注入深度,这可能解释了铍存在下 ITER 偏滤器的降解。 图:由31k 个原子组成的 NVE MD 模拟的模拟速率(ns/day)与 SNAP 电位中使用的描述符数目之间的关系。 图:拟合SNAP电位的工作流程。 图:(左)钨材料性能预测包括在 SNAP 优化回路表示为百分比误差的 DFT 预测。由于空位结合能(0.12 eV)的绝对值很小,EAM 的误差百分比比这里显示的 SNAP 电位大得多。(右)铍材料性能预测,也包括在优化程序。这里显示的两个 SNAP 电位都显著改进了现有的自组织形成能的 BOP 预测。文四: 加速材料设计的数据驱动策略摘要:近年来,随着机器学习和人工智能的出现,正在进行的自然科学革命引发了材料科学界的极大兴趣。可实现材料空间的固有高维性使得传统方法对于大规模探索无效。针对日益复杂的问题开发的现代数据科学和机器学习工具是一种有吸引力的替代方案。一场迫在眉睫的气候灾难要求在可能采取行动的几年内彻底改革现有技术,实现清洁能源转型。应对这场危机需要以前所未有的速度和规模开发新材料。例如,有机光伏有潜力在很大程度上取代现有的硅基材料,并开辟新的应用领域。近年来,有机发光二极管已成为数字屏幕和便携式设备的最先进技术,并通过柔性显示器实现了新的应用。网状框架允许原子精确合成纳米材料,并有望通过实现从气体储存、气体分离、电化学储能到纳米医学的多功能纳米颗粒的潜力,彻底改变该领域。在最近的十年里,模拟和机器学习在性质预测、性质优化和化学空间探索中的全面应用,以及计算能力和算法效率的显著进步,促进了所有这些领域的重大进展。在本报告中,我们回顾了我们团队在材料科学机器学习这一蓬勃发展的领域的最新贡献。我们首先总结了我们小组所涉及的最重要的材料类别,重点关注有机电子材料和晶体材料等小分子。具体而言,我们强调了我们用于加快发现和推导材料设计策略的数据驱动方法。随后,我们的重点放在我们团队开发和使用的数据驱动方法上,详细阐述了高通量虚拟筛选、反向分子设计、贝叶斯优化和监督学习。我们讨论了它们的总体思想、工作原理及其用例,并举例说明了在数据驱动的材料发现和设计工作中的成功实施。此外,我们详细阐述了这些方法的潜在陷阱和剩余挑战。最后,我们对该领域进行了简要展望,因为我们预计在材料发现和设计活动中会越来越多地采用和实施大规模数据驱动的方法。 图:热激活延迟荧光有机发射器的逆向设计工作流程,从选择碎片到设备集成和测试。 图:使用超分子变分自动编码器(SmVAE)的自动网状框架(RF)发现平台。 图:高通量虚拟筛选从大的候选空间开始(例如,如图所示,组合生成)。使用虚拟筛选,大多数候选人被淘汰,这样可以进行更少(更昂贵和耗时)的实验测试。 图:分子特性监督式学习的工作流程。已知(标记)数据集用于优化模型,该模型随后用于估计未知(未标记)数据集的分子特性。文五: 用于数据驱动材料设计的各种集成模拟(JARVIS)联合自动化存储库摘要:各种集成模拟联合自动化库(JARVIS)是一个集成基础设施,用于使用密度泛函理论(DFT)、经典力场(FF)和机器学习(ML)技术加速材料发现和设计。JARVIS的动机是材料基因组计划(MGI)原则,即开发开放存取数据库和工具,以减少材料发现、优化和部署的成本和开发时间。JARVIS的主要功能有:JARVIS-DFT、JARVIS FF、JARVIS-ML和JARVIS工具。到目前为止,JARVIS由大约40000种材料和大约100万个JARVIS-DFT计算特性、大约500种材料和≈110个力场以及大约25 ML JARVIS-ML材料特性预测模型组成,所有这些模型都在不断扩展。JARVIS工具提供用于运行和分析各种模拟的脚本和工作流。我们将计算数据与实验或高保真度计算方法进行比较,以评估预测中的误差/不确定性。除了现有的工作流程外,作为数据驱动材料设计范式的一部分,该基础设施还可以支持各种其他技术上重要的应用程序。 图:基于长度和时间尺度的计算材料设计技术。我们主要集中在计算方法的最低两个层次,DFT 和 MD,但我们结合其他模拟方法的具体应用。 图:JARVIS 基础设施概述。对于给定的材料性能指标,几个 JARVIS 组件可以一起工作,以设计优化的或全新的材料。 图:JARVIS工具包的三个主要组件及其功能。 图:JARVIS中的三个主要数据库及其内容摘要。来源:STEM与计算机方法

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