首页/文章/ 详情

PDE求解、力学模量表征、微裂纹定位、网格生成和本构模型反演

7月前浏览8334


文一:

 

https://doi.org/10.1016/j.camwa.2023.03.021

由具有尖锐梯度或不连续性的数据驱动的偏微分方程的发现

摘要:

在偏微分方程的数据驱动发现中,以前的研究人员已经成功地使用了各种方法来从平滑数据中推导参数估计,但不是从具有尖锐梯度或不连续性的数据中推导。为了捕捉数据中的尖锐梯度/不连续部分,我们在高斯过程先验中引入了Sigmoid函数的非零均值函数。我们使用无噪声和有噪声的数据对Burgers方程、无粘性Burgers方程式和非线性波动系统(NLWS)的回归问题和反问题进行了测试,并验证了该方法的有效性和鲁棒性。

 

图:总结本工作中开发的GPNOM方法的流程图。

 

图:带有跳转的Forrester函数的回归问题。(a) 具有零均值函数的GP。(b) GPNOM。

 

图:Burgers方程:预测的相对误差𝜆 具有无噪声数据。(a) 驻波。(b) 移动的波浪。

文二:

 

DOI:10.1109/TBME.2023.3315846

基于数据驱动的低分辨率CT松质骨各向异性力学特性表征

摘要:

探讨松质骨的各向异性力学行为对体内骨生物力学分析至关重要。然而,由于缺乏微观结构信息,在低分辨率(LR)临床CT图像下表征各向异性力学行为是一个挑战。本文提出的数据驱动方法从LR临床CT图像中准确地表征了松质骨的各向异性力学特性。方法:利用绵羊骨小梁立方体获得高分辨率(HR)显微CT和LR临床CT图像数据集。首先,使用HR图像数据来训练自动编码器模型。编码器提取微结构特征。训练快速超分辨率(FSR)模型以将LR骨立方体映射到从相应HR样本中提取的特征。预训练的FSR模型用于将LR临床CT图像转换为编码的微观结构特征。这些特征后来被用于预测目标组织形态参数、各向异性弹性张量和基于全连接神经网络的织物张量。结果:数据驱动模型准确预测了小梁骨的弹性张量和织物张量,LR CT图像的空间分辨率为0.6mm/像素。验证了LR临床CT图像可以使用生成深度学习模型和上采样操作生成微观结构信息。意义:本研究证明松质骨的临床医学图像可以用于分析复杂的力学性能,使用一种数据驱动的方法,可用于实时骨骼骨缺损诊断与个性化假体设计在临床应用中。

 

图:用带环形槽的标准圆柱形容器进行模型配准。

 

图:由微CT图像及其弹性表面建模的三次RVE。

 

图:整个步骤分为3个步骤。

文三:

 

https://doi.org/10.1016/j.actamat.2023.119073

高温载荷下CMC内部微裂纹演化行为的原位定量跟踪:一种深度学习方法

摘要:

微裂纹在陶瓷基复合材料的失效中起着极其重要的作用。微观裂纹演化行为的原位定量跟踪仍然是一个巨大的挑战。在本研究中,开发了一种基于生成对抗性网络的深度学习微裂纹分割方法,用原位X射线计算机微断层扫描(μCT)定量表征CMC在高温拉伸载荷下的微裂纹演化行为。该方法实现了对μCT图像中微裂纹的精确而稳健的分割,具有低灰度对比度和由噪声和伪影引起的图像质量。得到了各微裂纹的开孔面积、开孔位移和裂纹体积等裂纹参数。对于大多数微裂纹,这些参数的值在加载过程中没有明显增加。在高载荷水平下,一些大型微裂纹的裂纹参数显著增加。还跟踪了每个微裂纹的演变,以进一步确定主导最终失效的关键损伤。表现出更显著体积演变的四个最大裂纹被捕获为主要裂纹。这些主裂纹都起源于预先存在的微裂纹,通常是由相邻的小微裂纹随着拉伸载荷的增加而聚结而形成的。最终断裂被证明发生在这些主要裂缝附近。

 

图:(a) SiC/SiC样品的尺寸和光学图像。(b) 现场μCT设备,带实验室X射线源,温度升高。

 

图:(a) CT图像切片,带有矩阵、空隙、纬线和经线,用箭头标记。(b) CT图像的灰度分布,基于两个阈值将其划分为三个区域,分别指示空隙、纤维和基质的灰度范围。(c) (3D)从具有由箭头标记的不同相位的样本中提取的感兴趣区域的渲染体积。(d) -(f)分别从(c)中提取的基质、纤维和空隙。

 

图:荷载-位移曲线,插图显示了实验照片和荷载配置。

 

图:随着拉伸载荷的增加,SiC/SiC试样的微观结构演变。顶部:三维渲染体,包括用圆圈标记的裂纹扩展区域,底部:位于顶行虚线框位置的切片。

 

图:deeplab v3+网络的关键元素包括(a)巧妙的卷积,(b)深度可分离卷积和(c)不对称的空间金字塔池。(d) 修改后的deeplab v3+网络的详细架构。

 

图:采用不同方法得到的微裂纹分割结果。绿色表示被正确分割的像素,而红色表示被错误分割的像素。

 

图:(a) 以及(b)不同加载步骤下的裂纹张开位移场和裂纹张开面积分布。

 

图:基于μCT的断裂分析:(a)断裂形态,(b)不同加载步骤下微裂纹与断裂表面的距离。微裂纹的体积和开口面积由不同尺寸和颜色的符号表示,(c)主裂纹相对于断裂面的分布。

 

图:断口形貌的SEM图像。

文四:

 

https://doi.org/10.1016/j.cma.2023.116480

基于数据驱动的聚类降阶模型聚类相互作用矩阵的快速构建及弹塑性应力-应变曲线和屈服面的预测

摘要:

基于数据驱动的聚类降阶模型(CROM)方法在预测非均质材料的非线性有效性质方面非常有效。这些方法中的聚类交互矩阵(CIM)可以看作是一个聚类离散格林函数,是CROM方法的关键部分,越来越多的学者已经认识到这一点。然而,格林函数的显式表达式仅适用于均质材料。非均质材料的CIM需要通过数值方法进行计算,这面临着计算效率的挑战。因此,本文提出了一种构建非均质材料CIM的有效方法。通过在FCA控制方程下推导傅立叶空间中的格林函数,有效地计算了能提供自平衡应力场的均匀材料的CIM。然后,为了避免应力平均对簇的影响,我们在不引入额外计算成本的情况下,用非均质材料的自平衡应力基向量来扩充所提出的簇元相互作用矩阵。基于增广的簇元相互作用矩阵,改进了自平衡应力空间,并利用最小余能原理计算了高度非均质材料的CIM。所提出的算法显著提高了FCA的离线计算精度,能够有效预测复杂负载路径下的非线性有效特性和屈服面。数值算例验证了该方法的有效性和准确性。此外,我们还提供了FCA和SCA均质材料的CIMs之间的关系。基于本文提出的算法,可以对其他CROM方法有所帮助。

 

图:FCA的计算框架。

 

图:(算法1)基于FEM的CIM构造算法。

 

图:FCA示意图。

 

图:由算法2和算法1计算的同源CIMs的比较。(3D RVE)。

 

图:三维RVE模型和聚类模型。

 

图:当RVE最终产生时,由FCA和DNS计算的Mises应力分布。

文五:

 

https://doi.org/10.1016/j.neunet.2022.10.022

用于自动生成四边形网格的强化学习:温和的演员-评论家方法

摘要:

本文提出、实现并评估了一种基于强化学习(RL)的网格自动生成计算框架。网格生成在计算机辅助设计与工程(CAD/E)领域的数值模拟中发挥着重要作用。它被确定为美国国家航空航天局CFD 2030愿景研究的关键问题之一。现有的网格生成方法存在计算复杂度高、复杂几何形状中的网格质量低以及速度限制的问题。这些方法和工具,包括商业软件包,通常是半自动的,需要人类专家的输入或帮助。通过将网格生成公式化为马尔可夫决策过程(MDP)问题,我们能够使用最先进的强化学习(RL)算法,称为“软参与者-批评者”,从试验中自动学习网格生成的行动策略。这种用于网格生成的RL算法的实现使我们能够在没有人工干预和任何额外清理操作的情况下构建一个全自动网格生成系统,这填补了现有网格生成工具的空白。在与两个有代表性的商业软件包进行比较的实验中,我们的系统在可扩展性、可推广性和有效性方面表现出了良好的性能。

 

图:网格生成器为完成网格而采取的一系列操作。

 

图:用于自动网格生成的基于RL的计算框架。代理通过实现各种RL技术充当网格生成器。它在一种状态被感知后产生一个元素,并通过观察奖励来改善行动。环境对网格划分边界进行建模,并通过剪切通过动作生成的元素来更新边界。

 

图:每个规则类型的操作空间。子图(a)-(c)分别对应三种类型的行动。

 

图:使用FreeMesh RL对具有不同边界段密度的同一域进行网格划分。域1-3的边界形状相同,但边界上的顶点数量从低到高。通过将所有域与相同的训练模型进行网格划分来检查可伸缩性。

 

图:网格划分策略在四个域上的执行时间。网格划分策略主要由两种类型的操作组成:类型0和类型1。每个条上的数字表示每个动作类型的执行次数。

 

图:八种质量指标的网格化性能比较结果。BQ表示Blossom Quad方法;F-RL表示FreeMesh RL方法。L、 H分别指示是较低的值还是较高的值是优选的。


计算机技术在科学&技术&工程&数学中得到了广泛的应用,力学方面,计算机技术成为了科学的第四次革命性技术,现在基于计算机的数据科学已经逐步成为力学等其他科学发现的第四范式。人工智能、大数据、数字孪生等概念已经逐步成为当今时代的主题。智能制造、智能算法、数据驱动力学、大语言模型、自动驾驶在当今社会展现出巨大潜力,吸引了大量的研究人员。同时高性能显卡和多核中央处理器的出现为大规模数值模型的高性能计算提供了强大算力。然而因为该领域的论文较多,涉及内容较广,需要的知识量较大,不仅需要力学,数学,物理的知识,还需要计算机、数据科学、大数据分析的知识。

来源:STEM与计算机方法

附件

免费链接.txt
ACT断裂复合材料非线性化学动网格光学航空航天裂纹理论自动驾驶材料数字孪生
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-04-27
最近编辑:7月前
江野
博士 等春风得意,等时间嘉许。
获赞 48粉丝 47文章 310课程 0
点赞
收藏
作者推荐

计算力学前沿分享:颗粒复合材料、DEM冲击模拟、流固耦合、分子动力学

本文摘要(由AI生成):这篇文章的摘要基于 SPH 法对返回舱入水载荷进行了模拟分析。文章首先通过模拟平板和楔形体两种算例模型的入水过程,并与相关文献的试验结果对比,验证了算法的有效性。然后,对返回舱的入水过程进行模拟,结果显示返回舱入水速度和倾角对砰击有较大影响,并且存在二次砰击现象。最后文章表明,SPH 气-液两相流算法能够较好地模拟返回舱入水过程。文一: DOI:10.7511jslx20220525002一种颗粒复合材料物理性能计算模型与工程应用摘要:结合傅里叶描述子法与格林函数法提出了任意形状夹杂物的Eshelby张量计算方法,基于分步夹杂法建立了一种具有任意形状颗粒的复合材料物理性能预测理论模型。通过傅里叶描述子法对扫描电镜(SEM)图像中颗粒形貌进行数学表征,构建由多个细观结构参数组合的几何模型;利用格林函数求解夹杂物Eshelby张量,考虑颗粒形态、大小和组分性能,研究细观结构对颗粒复合材料有效性能的影响。在此基础上,设计了一款颗粒复合材料物理性能预测软件。以铜铬合金为研究对象,分析重构的细观模型与预测的结果,并与有限元计算结果进行对比,验证了该计算模型的有效性和可靠性,以及在工程计算应用中的可行性。 图:软件人机交互界面 图:SEM和合金重构细观结构 图:弹性性能计算结果显示界面 图:电导率预测可靠性验证文二: DOI:10.7511/js1x20220315002不同重力场下颗粒冲击过程的离散元分析摘要:开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟,并与地球重力下冲击的试验结果进行对比验证。在此基础上,进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明,在所有重力加速度下,冲击物的穿透深度d与冲击速度v。的关系可以用Poncelet模型表达;d与冲击物下落的总高度H表现为d~H”的幂律关系,当H<10m时,d与H的幂率标度为0.322,而H>10m时,d与H的幂率标度下降到0.211。此外,穿透深度小于冲击物半径时,重力加速度对冲击物减速过程无影响。在所有的重力加速度下,当冲击速度大于5m/s时,冲击物的持续碰撞时间t。是恒定的,且与重力的一1/2次方呈正比。 图:H1=0.5cm时,冲击力变化过程的试验结果和模拟结果与本文离散元计算值对比。 图:不同颗粒厚度下试验结果和模拟结果与本文离散元计算获得的底板最大冲击力 图:v=15m/s时,不同阶段的颗粒床切片图 图:冲击物的最终穿透深度随重力加速度和冲击速度的变化,内插图为系数α和3随重力加速度的变化文三: DOI:10.7511js1x20220610002地震和风耦合作用下连体高层建筑联合易损性分析方法研究摘要:摘要:近年来,城市地标建筑的建设中,高层连体建筑越来越受到青睐,在其服役期内不可避免地受到风和地震等多种灾害作用。而单一灾害作用下的结构失效概率低估了结构损伤乃至倒塌的可能性,迫切需要多灾害作用下的结构安全性能评价方法。本文以苏州国际会议酒店的连体高层建筑为对象,开展了地震和风耦合作用下的联合易损性分析方法研究。首先,建立了地震和风耦合作用下的多灾害易损性函数和结构性能需求函数,形成了多灾害联合易损性分析方法。随后,选取了地震动记录,开展了基于数值风洞的结构风荷载模拟,并进行了单一灾害作用下的结构易损性分析。最后,基于均匀分布抽样了100组地震和风耦合作用的荷载组合,开展了结构多灾害联合易损性分析。结果表明,结构易损性随着PGA和风速的增大而增大;相比单一灾害作用,地震与风耦合作用下,结构各破坏状态的失效概率均有所增加;当风速达40/s时,罕遇地震作用下结构轻微破坏、中等破坏、严重破坏和倒塌破坏的概率相比仅地震作用时分别增加了1.13%,35.67%,84.00%和98.33%。本文研究为复杂高层建筑多灾害安全性能评价提供了可用工具和参考。 图:苏州市国际会议酒店项目 图:地震动记录 图:结构数值风洞模拟 图:地震与风耦合作用下联合易损性曲面与位移时程文四: DOI:10.7511/js1x20220615002分子动力学模拟含不同初始缺陷单晶镍的断裂行为与应力演化摘要:断裂是一个跨尺度复杂的物理过程,对宏观尺度的断裂行为已有深入的研究和发展,然而对微观尺度的断裂行为及断裂过程中应力场的变化缺乏深入的理解。本文通过分子动力学模拟,研究了具有不同初始缺陷(尖锐裂纹、钝裂纹和孔洞)的单晶镍的断裂行为和应力分布特征。结果表明,不同的初始缺陷导致了不同的断裂机制、断裂强度和抗断裂性能。含初始孔洞的单晶镍样品有最高的断裂强度和最强的抗断裂性能,这与孔洞扩展过程中堆积层错的形成密切相关。其次是含初始钝裂纹的样品,在裂纹扩展过程中出现由[100]超位错发射引起的裂尖钝化;含尖锐裂纹的样品表现为脆性断裂,裂尖原子没有出现微结构的变化,其强度和抗断裂性能最低。此外,不同的初始缺陷也会导致断裂过程中应力分布的变化,对含有尖锐裂纹的脆性断裂试样,高应力(拉伸应力、平均应力和米塞斯应力)总是出现在扩展裂纹的裂尖。而对于含有钝裂纹或孔洞的韧性断裂试样,高应力不仅分布在裂尖,也分布在位错发射和堆积层错形成的区域,在裂纹/孔洞扩展之前,应力随着加载时间的增加而迅速增加,而一旦裂纹或孔洞开始扩展,应力增加非常缓慢或几乎不增加,但拉伸应力值始终大于平均应力和米塞斯应力值。这表明,在I型加载条件下,纳米尺度下材料的断裂由拉应力(正应力)控制。 图:含不同初始缺陷样品的几何形状 图:含不同初始十陷样品的原子应力-应变曲线和断裂强度。 图:I型载荷作用下含初始缺陷的单晶镍拉伸应力场变形云图 图:三种不同样品的微观结构特征 图:t=310s处不同初始缺陷样品的裂纹/孔洞张开位移和不同初始缺陷样品的裂纹长度随加载时间的变化文五: DOI:10.7511jslx20220528002基于SPH法的返回舱入水载荷模拟分析摘要:返回舱水上着落前期面临较大砰击,为研究过程中受到的砰击载荷,基于光滑粒子流体动力学(SPH)气-液两相流算法,首先通过模拟平板和楔形体两种算例模型的入水过程,并与相关文献的试验结果进行对比,验证算法的有效性。在此基础上,对返回舱的入水过程进行模拟。结果显示,两种算例模型的计算结果与相关文献试验结果吻合良好。返回舱入水速度和倾角对砰击有较大影响并且过程中存在二次砰击现象。砰击载荷随入水速度增大而增加。第一次砰击载荷峰值随倾角增大而减小,第二次砰击载荷峰值随倾角增大先增后减。结果表明,SPH气-液两相流算法能够较好地模拟返回舱入水过程。 图:计算模型(单位:mm) 图:平板入水压力云图 图:楔形体入水过程中液面变化对比 图:返回舱模型 图:不同倾角下自由液面压力云图及姿态变化来源:STEM与计算机方法

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