首页/文章/ 详情

使用 Abaqus 模拟增材和减材工艺过程

1天前浏览15


对于建筑和制造中的许多增材和减材工艺,有限元分析 (FEA) 可有效地用于了解添加或去除材料的影响。例如混凝土、塑料、金属等。FEA 使用户能够虚拟评估传热效果、自重变化或与这些过程相关的其他负载场景,从而提高计算精度并减少原型设计需求,以评估由此产生的位移形状、应力状态、潜在裂纹、和部分变形。数字模型中的材料可以立即添加或删除。或者,根据逐渐的材料变化,在进行实物投资之前指导合适的材料选择和工艺优化。


建筑和制造中的材料添加和去除


材料的添加用于建筑和制造过程。滑模混凝土结构和碾压混凝土被放置在升降机中,材料依次添加到结构中。热量作为混凝土固化的副产品释放,引起热应力,从而导致开裂。新材料贡献的额外质量将由本身固化并获得强度的混凝土承载。该过程与增材制造的过程类似。尽管规模和材料有所不同,但评估材料添加影响的建模过程是一致的。焊接是其中的一个子集,涉及高温和渐进式材料添加。

与此相反的是在物理、化学或其他加工过程中去除材料。用于去除材料的技术会产生热量,热量会传递到剩余的材料上。由此产生的热膨胀可能会导致暂时和永久的变形。并且,要去除的材料量应考虑零件的最终收缩和可能的变形。


如何使用 Abaqus 虚拟评估和优化材料添加和去除


SIMULIA Abaqus 软件提供了多种选项来模拟材料的去除或添加;状态变化、材料场变量和渐进元素激活。接触激活/停用可以单独使用,也可以与可用的不同方法一起使用,在 Abaqus 2022 中,此功能已扩展到 Abaqus Standard 中的一般接触。在这篇博文中,我们将了解如何在 Abaqus/Standard 中应用状态和材料场变量的变化。

状态的改变允许元素在一个步骤内被停用或激活,而场变量可用于改变材料属性以模拟移除或添加的效果。当需要逐渐增加或减少刚度而不是状态的瞬时变化时,场变量可能是有利的。这些程序可以独立或串联应用以达到所需的结果。渐进式元素激活是一种保真度较高的方法,允许完全和部分元素激活。有关此方法的更多信息可以在 Abaqus 文档中找到。


 模型设置


将使用说明性示例来解释这些过程,因为将使用现场变量和模型更改方法来实现相同的结果。这些示例将重点关注应力/位移元素的背景;然而,这些方法可以得到更广泛的应用。矩形零件的一部分将被分段移除为八个单元。几何特性如表 1 所示,材料特性如表 2 所示。
 
表 1:几何特性  


 
表 2:材料特性  


初始和最终几何形状如图 1 所示。通过单独去除或降解零件内的八个单元,零件从初始配置转变为最终配置。



 

图 1:初始和最终几何形状  

无论选择哪种方法,此过程的第一步都是对几何结构进行分区,以确定要移除的单元格。零件几何形状应反映最大尺寸 - 如果材料被去除,这将反映初始几何形状;如果添加材料,这将反映最终的几何形状。通过划分部分,可以将各个组分配给要删除的单元。


对于此示例,零件几何形状被划分为 3 × 3 × 4 单元,如图 2 所示。然后根据零件级别的几何形状将集 合分配给各个单元。当需要重新网格化并且在装配体中以partname-instance.setname 的形式访问这些集时,此方法非常可靠。这里,这些组是根据它们在零件形式X1-Y1-Z1X3-Y3-Z4中的 x、y、z 位置来命名的 - 后面的零件在图 2 中突出显示。零件名称为Part-1,第一个零件为 Part-1。并且该部分的唯一实例被标记为Part-1-1 。随后,集 合X3-Y3-Z4在装配级别被识别为Part-1-1.X3-Y3-Z4 。



 
图 2:分区名称和集名称  


关键字编辑

场变量的空间变化用于引发局部材料降解,以及随后的材料去除或添加的效果。材料属性的编写依赖于指定的字段变量。温度依赖性是该通用框架的一个具体示例。通过在步骤内更改字段变量的值,属性会在模拟内的指定时间发生更改。振幅与所需组的场变量相关联,以管理空间变化,并随后控制哪些单元受到材料属性变化的影响。在此示例中,杨氏模量和泊松比降级到对模型内八个单元的零件刚度影响可以忽略不计的程度。


图 3 显示了具有场变量依赖性的材料属性。场变量数设置为 1。当场 1 等于 1 时,杨氏模量值为 200 000 MPa,泊松比为 0.3。当 Field 1 值增加到 2 时,这些属性分别降低到 2 MPa 和 0.01,可以忽略不计。


 
图 3:材料特性  


确定一个或多个幅度以引入场变量的空间变化并管理在时间步长上元素的移除。为要去除的八个细胞中的每一个设置唯一的幅度,并且通过改变表格值来交错去除。在图 4 中,显示了与X3-Y3-Z4组相关的表格振幅。与幅度相关的场变量将以 4 的步长时间从 1 增加到 2,并保持等于 2,直到步长结束。这是该步骤中要降级的第四组,总共使用 8 秒的步骤时间来单独降级每组


 
 图 4:幅度  


通过在 Model -> Edit keywords -> ModelName中找到的关键字编辑器中对模型进行修改,将一个或多个 amThe 振幅与唯一集关联起来。
在步骤 1 之前,将字段变量设置为等于 1 的初始条件,如图 5 所示。


 
图 5:设置初始条件的关键字  

在所需的步骤中,根据基于时间的幅度更改场变量。这里,幅度函数Amp-X3-Y3-Z4与装配级别和场变量 1 上的Part-1-1.X3-Y3-Z4集相链接,如图 6 所示。


 
图 6:将幅度链接到集 合的关键字  


通过定义具有场变量依赖性的材料,编写多个振幅函数,并将每个振幅函数与场变量和模型中的唯一集相关联。通过分析可以改变不同位置的材料特性。我们快速降低属性以模拟材料去除的效果。但是,对于材料添加来说,相反的情况也是可能的。此外,如果建立了材料强度曲线,则可以使用振幅函数来定义强度随时间逐渐增加或减少。弹性体的老化是随着时间的推移在负载、温度下发生降解。使用此技术时,应注意应用边界条件和载荷。这是为了避免材料性能下降时发生元件变形。

 型号变更

模型更改类型交互是在步骤期间启动模型区域的添加或删除。移除后,Abaqus/Standard 会存储该区域在模型剩余部分边界上的节点处施加的力/通量。在移除步骤中,这些力逐渐降至零。这保证了模型的平滑效果,不妨碍收敛。如果模型更改是为了添加材料,则初始步骤将是有效的最终几何形状。应在第一步中使用模型更改交互删除将添加到模型的单元,然后可以在后续步骤中根据需要添加单个单元。为应力/位移单元提供两种类型的重新激活:无应变重新激活和有应变重新激活。有关这些选项的更多信息,请参阅 Abaqus 文档。
第二个模型显示了使用模型更改交互来移除细胞。在图 5 中,模型更改了X3-Y3-Z4的交互集。正如对字段变量的使用所做的那样,该集 合是我们从模型中删除的第四个集 合。为了实现顺序删除,设置X3-Y3-Z4随后在步骤 4 中停用。为了同时移除多个单元格,可以定义包含所有所需单元格的备用组。在“编辑交互”对话框中,为名为Part-1-1.X3-Y3-Z4的集按集选择几何区域。我们可以选择删除/添加的其他区域类型是皮肤、纵梁或元素。默认情况下,“停用”选项处于启用状态,因为在模型更改分析中,停用必须在重新激活之前发生。
 
图 7:模型变更交互  


材料去除/添加可能会导致数值问题,需要特别注意。为了避免这些问题和其他常见问题,请遵循以下准则:


  • 边界条件:静态应力分析中材料的去除。确保模型的其余部分受到充分约束以防止出现无约束刚体模式非常重要。
  • 载荷:在元件被移除或重新激活的区域施加的分布式和集中载荷可能需要修改。
  • 元件:目前不适用于刚性元件、粘性元件、垫圈元件和压电元件。可以删除和重新激活 Abaqus/Standard 中的所有其他元素类型。
  • 接触:如果连接到接触对的元件被移除,则在启动模型更改时应移除或停用接触对。
  • 约束:如果要删除附加到受多点约束或线性约束方程约束的节点的所有元素,则该节点应为多点约束或线性约束方程的从属节点。
  • 警告:在某些情况下,元素删除可能会导致 Abaqus/Standard 在消息文件中报告额外的未连接区域。忽略这些消息是安全的。
在此示例中,模型更改交互的使用是从装配体中单独移除每个组。对于顺序移除,每组都处于单独的时间步长中,并注意确保随后调整受移除影响的载荷和边界条件,避免刚体运动。在添加材料的情况下,所添加的材料将首先被删除,然后在适当的时候添加回模型。您可以根据模拟需求将字段变量依赖项添加到包含模型更改交互的模型。



来源:ABAQUS仿真世界
Abaqus化学通用建筑增材焊接裂纹减材材料控制装配
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-05-30
最近编辑:1天前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 169粉丝 306文章 373课程 0
点赞
收藏
作者推荐

《Mechanics of Solid Polymers》5.3.7Mooney-Rivlin模型

5.3.7 Mooney-Rivlin模型 Mooney-Rivlin (MR) 模型是NH模型的扩展,通过在单位参考体积的亥姆霍兹自由能中加入对I2*的线性依赖,试图提高预测精度:如方程所示,MR模型的可压缩版本需要三个材料参数:C10、C01和κ。使用方程(5.41)可以证明,Mooney-Rivlin模型的柯西应力表达式为:对于MR模型的不可压缩版本(κ → ∞),单轴、平面和等双轴变形的柯西应力由以下表达式给出:Mooney-Rivlin模型预测弹性体行为的准确性在图5.15中通过与Treloar[16]的硫化天然橡胶数据对比得到了展示。图5.15 Treloar[16]的实验数据与不可压缩Mooney-Rivlin材料模型预测的对比。 这个图表明,Mooney-Rivlin模型能够改善NH模型的预测。然而,这种改进可能会付出代价——具有负C01项虽然通常能改善一种加载模式下的预测,但同时可能使模型在另一种有限变形下变得不稳定。这种情况的一个例子如图5.16所示。对于不可压缩单轴加载,Mooney-Rivlin材料模型可以用以下Python代码实现:Python代码:"MR_incompressible_uniaxial.py"from pylab import *defMR(strain, params):"""Mooney-Rivlin超弹性模型。 不可压缩单轴加载。 此函数使用真实应力和应变""" C10 = params[0] C01 = params[1] lam = exp(strain)return2 * (lam*lam - 1/lam) * (C10 + C01/lam)strain = linspace(0, 0.8)params = [1.0, 0.1]stress = MR(strain, params)plot(strain, stress, 'r-')show()以下代码示例展示了一种实现可压缩单轴加载的Mooney-Rivlin材料模型的方法。这段代码展示了Mooney-Rivlin模型的3D实现和可压缩单轴加载的应用。让我解释一下这两部分代码:"Polymer_Mechanics_Chap05.py"中的附加Python代码:这个函数MR_3D实现了3D Mooney-Rivlin模型:defMR_3D(stretch, param):"""Mooney-Rivlin。3D加载由拉伸指定。 param: [C10, C01, kappa]""" L1, L2, L3 = stretch[0], stretch[1], stretch[2] F = array([[L1,0,0], [0,L2,0], [0,0,L3]]) J = det(F) bstar = J**(-2.0/3.0) * dot(F, F.T) bstar2 = dot(bstar, bstar) I1s = trace(bstar) I2s = 0.5 * (I1s**2 - trace(bstar2)) C10, C01, kappa = param[0], param[1], param[2]return2/J*(C10+C01*I1s)*bstar - 2*C01/J*bstar2 + \ (kappa*(J-1) - 2*I1s*C10/(3*J) - 4*I2s*C01/(3*J))*eye(3)"MR_compressible_uniaxial.py":这个脚本使用上面定义的MR_3D函数来计算可压缩单轴加载下的应力-应变关系:from pylab import *from Polymer_Mechanics_Chap05 import *trueStrain = linspace(0, 0.8, 100)trueStress = uniaxial_stress(MR_3D, trueStrain, [1.0, 0.1, 100])plot(trueStrain, trueStress, 'r-')xlabel('真实应变')ylabel('真实应力 (MPa)')grid('on')show()这段代码计算了一系列真实应变对应的真实应力,并绘制了应力-应变曲线。它使用了之前定义的uniaxial_stress函数(在前面的代码中定义)来计算单轴应力。这个实现展示了Mooney-Rivlin模型如何处理可压缩材料的单轴加载情况。它考虑了体积变化(通过参数kappa),这使得模型更加通用,可以应用于更广泛的材料和加载情况。来源:ABAQUS仿真世界

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