首页/文章/ 详情

开源代码大放送,你想要的全都有:拓扑优化,图像识别,XFEM,进场动力学

20小时前浏览8

文一:

 

用于复杂3D 拓扑优化的简单而紧凑的 Python 代码

摘要:

本文提出了一个用于一般三维拓扑优化的100行Python代码。该代码采用Abaqus脚本接口,可方便地访问高级有限元分析(FEA)。它是使用双向进化结构优化(BESO)方法为体积约束下的柔度最小化而开发的。源代码由控制迭代过程的主程序和实现输入模型准备、有限元分析、网格无关滤波器和BESO算法的五个独立函数组成。该代码从模型数据库(.cae文件)中读取初始设计,该数据库可以是Abaqus/cae中生成的任意3D几何体,也可以是从各种广泛使用的CAD建模包转换而来的。这种结构良好的代码可以方便地扩展到各种其他拓扑优化问题。作为对代码进行简单修改的例子,给出了对多种负载情况和非线性的扩展。此代码对拓扑优化领域的研究人员和寻求自动化概念设计工具的实践工程师非常有用。通过进一步的扩展,该代码可以解决结构工程、机械工程和建筑实践中复杂的三维概念设计问题。完整的代码在附录部分给出也可以从网站上下载www.rmit.edu.au/research/cis m/。

 

图:三维悬臂的拓扑优化:(上)以初始完整设计为设计域;优化解决方案的(中间和底部)侧视图和透视图

 

图:具有初始空隙和不可设计部分的拓扑优化:(顶部)设计域(四分之一模型);全桥优化解决方案的(中间和底部)侧视图和透视图。

文二:

 

一种有效的多裂纹扩展可变节点 XFEM 模型—— Matlab 面向对象实现

摘要:

我们提出了一种用于脆性材料多裂纹扩展模拟的高效变节点扩展有限元方法(XFEM)的Matlab面向对象实现。为了提高该方法的精度,不仅采用了局部网格细化策略,还采用了可变节点单元直接连接不同尺度的单元,节省了计算成本。为了有效地进行聚结检测,还提出了一种改进方案。进行稳定性分析以确定主动裂纹尖端,并采用裂纹长度控制方案对多次裂纹扩展进行建模。网格通过直接的方式与裂纹区域的扩展同时进行局部细化。描述了用于可变节点XFEM的Matlab面向对象编程代码的通用计算机程序结构。通过线弹性断裂力学分支的五个数值例子,证明了本工作的所有这些理想特征。此外,这项工作还提供了一个由Matlab编程语言编写的用于计算裂缝模拟的高效计算机代码,供教育和学术使用。

 

图:当尖端接近另一个裂纹时,积分域的修正:(a)该域被另一个裂缝穿过;(b) 该域与即将到来的裂缝无关。

 

图:多尺度网格示意图。(a) 原始大尺度网格;(b) 细化后的多尺度网格

 

图:边裂板的尺寸、加载条件和网格: (a)拉伸状态下的边裂板,(b)初始粗网格,(c)初始时间步中的多尺度网格,(d)标准小尺度网格,(e)初始时间步中的增强策略放大视图。

文三:

 

均匀和非均匀材料中动态断裂的开源近场动力学代码

摘要:

我们提出了一种新的开源近场动力学软件,用于模拟均质和非均质固体中的二维和三维动态裂缝。该软件及其实现进行了详细描述,并在GitHub上提供,并获得了MIT的许可。通过包括波浪传播问题在内的几个数值例子证明了该软件的能力。此外,在细观尺度上对聚合物基复合材料的三维断裂模拟显示了该代码对复杂微观结构建模的能力,包括界面断裂和本体中断裂之间的相互作用。

 

图:Peridigm 开发团队报告的 Peridigm 架构及其可扩展区域

 

图:20 m/s恒定速度下的波传播(左半部分)和裂缝生长(右半部分)作为应用于板的边界条件

文四:

 

基于数字图像处理的异质材料建模开源软件包 DigiSim

摘要:

提出了一个基于数字图像的有限元和离散元数值模型生成软件包。支持2D位图和3D体积图像。该软件包包括用于矢量化位图图像、计算几何信息统计和生成数值模型的几何和分析例程。处理复杂连接区域的能力使具有复杂微观结构的异质材料能够建模。四个例子表明该代码健壮、高效且用户友好。作为一个开源的、随时可用的Matlab软件包,它在非均质材料的力学行为分析中具有广阔的应用前景。

 

图:DigiSim软件包的一般流程图

 

图:使用x射线显微成像的混凝土典型微观结构:(a)混凝土的3D CT图像,(b)预处理后的2D CT切片

 

图:多相图像分析的示例:(a)彩色图像,(b)灰色图像,(c)孔隙的二元图像,(d)聚集体的二元图,(e)两相的几何形状,(f)使用Abaqus的FEM网格

文五:

 

基于 MPI/Fortran 混合阵列多尺度 CAFE 框架的高性能混凝土非均质材料断裂模型

摘要:

描述了一个用于非均质材料断裂模拟的三维多尺度元胞自动机有限元(CAFE)框架。该框架是在一个混合 MPI/Fortran 协同阵列代码中实现的,用于在 HPC 平台上进行高效的并行执行。使用了现代 Fortran 开发的两个开源 BSD 许可库: CGPACK,使用 Fortran 协阵列实现细胞自动机(CA) ; ParaFEM,使用 MPI 实现有限元素(FE)。该框架实现了结构层次(FE)和微结构层次(CA)之间的双向并发层次信息交换。描述了 MPI 到协阵列的接口和数据结构。利用 CAFE 框架对多晶铁圆棒在拉伸状态下的超晶解理传播进行了预测。通过这个 CAFE 框架实现的新结果包括模拟通过单个晶粒和跨晶粒边界的渐进解理传播,以及在单个晶体的优先取向解理平面上裂纹合并时出现宏观裂纹。CGPACK 和 ParaFEM 在 Cray XE6的前期工作中分别证明了几乎理想的至少可以扩展到数万个核的强扩展。在这项工作中使用了 Cray XC30和 XC40平台以及 CrayPAT 分析。最初 CGPACK/ParaFEM 混合 CAFE 模型的强缩放极限是2000核。在用最近邻算法取代所有到所有的通信模式后,Cray XC30的强缩放限制增加到7000核。TAU 对非 Cray 系统的分析发现了英特尔 Fortran 16在远程协同阵列操作优化方面的缺陷。最后,讨论了 CA 中线程并行化面临的挑战和机遇。

 

图:具有1百万晶粒的模拟等轴微观结构的示例

 

图:CAFE的实现1显示了晶界(顶部图像)和宏观裂纹(底部图像)

来源:STEM与计算机方法
HPCAbaqus断裂复合材料非线性拓扑优化通用建筑MATLABpython通信离散元裂纹理论多尺度控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:20小时前
江野
博士 等春风得意,等时间嘉许。
获赞 56粉丝 83文章 153课程 0
点赞
收藏
作者推荐

AIforScience的交叉研究

文一: 人工智能时代的科学发现摘要:人工智能(AI)正越来越多地融入科学发现,以增强和加速研究,帮助科学家生成假设、设计实验、收集和解释大型数据集,并获得仅靠传统科学方法可能无法获得的见解。在这里,我们考察了过去十年的突破,包括自我监督学习,它允许在大量未标记的数据上训练模型,以及几何深度学习,它利用有关科学数据结构的知识来提高模型的准确性和效率。生成人工智能方法可以通过分析包括图像和序列在内的各种数据模式来创建设计,如小分子药物和蛋白质。我们讨论了这些方法如何在整个科学过程中帮助科学家,以及尽管取得了这些进步,但仍然存在的核心问题。人工智能工具的开发人员和用户都需要更好地了解这些方法何时需要改进,而糟糕的数据质量和管理带来的挑战仍然存在。这些问题跨越了科学学科,需要开发有助于科学理解或自主获取的基础算法方法,使其成为人工智能创新的关键关注领域。 图:人工智能时代的科学。科学发现是一个多方面的过程,涉及几个相互关联的阶段,包括假设形成、实验设计、数据收集和分析。人工智能准备通过在这一过程的每个阶段加强和加速研究来重塑科学发现。这里展示的原理和说明性研究突出了对增强科学理解和发现的贡献。 图:学习科学数据的有意义的表示 图:人工智能引导科学假设的产生。 图:人工智能与科学实验和模拟的集成。文二: 应用统计物理学对抗人工智能增强分子动力学的陷阱摘要:基于人工智能的方法通过从原始数据中提取相关信息的能力,在整个科学领域产生了无可争辩的影响。最近,人工智能也被用于提高分子模拟的效率,其中人工智能衍生的慢模式被用于以有针对性的方式加速模拟。然而,尽管使用人工智能的典型领域的特点是数据过多,但分子模拟在每个结构中都受到有限的采样,因此数据有限。因此,在分子模拟中使用人工智能可能会遇到一种危险的情况,即人工智能优化可能会陷入虚假状态,导致手头问题的反应坐标(RC)的错误表征。当这样一个不正确的RC被用来进行额外的模拟时,人们可能会开始逐渐偏离基本事实。为了解决虚假人工智能解决方案的问题,我们报告了一种新颖的自动化算法,该算法使用了统计力学的思想。它基于这样一种理念,即更可靠的人工智能解决方案将是最大限度地实现慢过程和快过程之间的时间间隔。为了从有限的数据中学习这种时间尺度的分离,我们使用了一个基于最大口径的框架。我们展示了这种自动方案对三个经典基准问题的适用性,即模型肽的构象动力学、配体与蛋白质的脱键以及蛋白质G的C末端结构域的折叠/展开能量景观。我们相信,我们的工作将使可信赖的人工智能在复杂系统的分子模拟中得到更多和稳健的使用。 图:描述丙氨酸九肽构象动力学的RC的伪AI解。 图:流程图说明了我们筛选人工智能解决方案的新颖且计算高效的协议。从短期无偏MD模拟开始,我们的协议自动筛选在基于AI的方法中获得的虚假解,并学习最优RC。在这项工作中,我们证明了我们的协议在RAVE环境中的适用性,并且通过基于路径熵的过程来实现对虚假解决方案的筛选。 图:捕获丙氨酸二肽中的假AI溶液 图:从无偏MD中捕获GB1-C16的动力学文三: 利用人工智能辅助分子动力学发现蛋白质构象的灵活性摘要:蛋白质对不同于其晶体结构的各种构象进行取样。这些结构、它们的倾向以及它们之间的移动路径包含了大量关于蛋白质功能的信息,这些信息从纯粹的结构角度来看是隐藏的。分子动力学模拟可以揭示这些替代构象,但通常需要高昂的计算成本。在这里,我们应用我们最近的统计力学和基于人工智能的分子动力学框架来增强蛋白质环的采样。我们通过对经典测试片蛋白T4溶菌酶的三个突变体的研究来举例说明这种方法。我们能够根据它们激发态的稳定性对它们进行正确的排序。通过分析反应坐标,我们还获得了至关重要的见解,了解为什么序列空间中的这些特定扰动会导致构象灵活性的巨大变化。因此,我们的框架允许在最小的先前人类偏差的情况下准确比较环构象群体,并应直接适用于生物学、化学及其他领域的一系列大分子。 图:单突变体 L99A 的晶体结构与腔相邻环(残基100-120)显示为蓝色 图:每个突变体的两个RC组分中的每个组分中三个最高加权OP的视觉表示。 图:在500ns的元动力学模拟之后,使用每个状态的两个复 制品来确定激发态和基态之间的自由能差。文四: 用于分析和增强分子动力学模拟的机器学习方法摘要:由于计算能力和软件可用性的提高,分子动力学(MD)已成为研究生物物理系统的强大工具。尽管MD为更好地理解这些复杂的生物物理系统做出了许多贡献,但仍有方法上的困难需要克服。首先,如何使运行中产生的大量数据,即使是一微秒长的MD模拟,也能让人理解。其次,如何有效地对下面的自由能表面和动力学进行采样。在这个简短的视角中,我们总结了基于机器学习的思想,这些思想解决了这两个局限性,重点是它们的关键理论基础和剩余的挑战。 图:示意图说明了使用机器学习来分析和增强MD模拟的一些方法的典型工作流程。描述系统在配置空间中的时间演化的高维数据被用作人工神经网络的输入。人工神经网络被训练以将输入投影到低维空间。根据神经网络和目标函数的结构,低维表示捕获了被认为是重要的不同特征,例如慢模式。在一些方法中,所学习的特征用于进一步增强MD模拟的采样,如底部箭头所示。文五: 通过基于简化的对特异性从头算力场的分子动力学模拟,研究了由甲烷、二氧化碳和丙烷组成的二元流体混合物的菲克扩散系数摘要:为了用分子动力学(MD)模拟计算流体的热物理性质,经常使用针对实验汽液平衡进行优化的有效力场(FF)。该FF类别的示例有TraPPE、OPLS和AMBER。另一种选择是简化的基于从头算的配对特异性FF(AI FF),它是从零密度极限下的量子化学计算中结合气体动力学理论得出的,因此具有预测性。在本研究中,我们表明,在选定的二元流体混合物的帮助下,MD模拟在计算Fick扩散系数时的预测能力可以使用基于纯物质相应FF的简化对特异性AI FF来提高。为了评估新型AI FF与TraPPE FF的性能,研究了由甲烷、二氧化碳和丙烷组成的二元混合物,从过热蒸汽到气态,从超临界区域到压缩液态。为了确定在(0.1和12)MPa之间的压力、(293和355)K之间的温度以及0.05和0.95之间的摩尔分数下的Fick扩散系数,考虑到系统尺寸效应,对MaxwelleTefan扩散系数和热力学因子的分析进行了单独的模拟。除了压缩液体状态和两相边界附近的区域,其中TraPPE FF通常是优越的,AI FF显示出对Fick扩散系数的改进预测,平均扩展统计不确定性为12%。这可以通过将模拟结果与理论从头计算和可用的实验数据进行比较来证明,导致AI FF和TraPPE FF的平均绝对偏差分别为7%和13%。 图:KirkwoodeBuff系数Gij是由CH4(物种1)和CO2(物种2)组成的等摩尔二元混合物在295K和10MPa下R的倒数的函数,基于具有AI FF的20ns的模拟运行 图:不同对CH4eCH4(实线)、CO2eCO2(虚线)和CH4eCO2(点线)的集体 位移Dri$Drj的乘积(上部),以及由CH4和CO2组成的等摩尔二元混合物在295K和0.1MPa下的相应b系数(下部),通过使用我们的AI FF的1ms的五次独立MD模拟运行(彩色线)计算得出 来源:STEM与计算机方法

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