首页/文章/ 详情

基于无网格(mesh-free)策略实现单积分点几何必须为错(GND)的计算

11月前浏览564
参考文献:《Physically based crystal plasticity FEM including geometrically necessary dislocations: Numerical implementation and applications in micro-forming》
GND 演化方程依赖依赖于剪切应变率的梯度或者塑性变形梯度的旋度,而标准FEM/VUMAT 只告诉你每个积分点本身的   、Fp,不会直接给梯度。以往广泛应用的数值方案通常是:先把 积分点的数据外推到节点,再用线性形函数求梯度,然而这类方案只能用特定单元(如 C3D8),对自适应网格、复杂接触不友好。
该文章提出的一个mesh-free的方案,该方案的主要优势是不改单元、不加 DOF,只在材料子程序内部,用邻近积分点的数据做一次局部重构,就算出梯度,该策略对某个积分点   ,附近有一团“邻居积分点” ,作者把它们当成 mesh-free 的“节点”,对每个场变量 (可以是 Fp 的分量)做 MLS 拟合,如下图所示:
权函数使用立方样条,有紧支撑,距离越近权越大:
在实现上作者提到,立方支撑三个方向尺寸约为5个单元尺寸,最多取最邻近60个(3D)或者30个(2D)积分点,作者指出:当邻域尺寸比网格尺寸还小的时候,这个非局域模型就自然退化为传统的局域模型。也就是说,邻域尺寸本身就扮演了“材料内在长度标度”的角色。
为了提高计算效率秒作者使用了一个“时间滞后 + 公共块”的策略对GND进行更新。
作者使用的方案对于显示大变形分析计算效率非常高,使用标准的C3D8R单积分点即可正常运行,并将所提出的数值模型应用于铜箔拉伸和杯冲过程中的尺寸效应分析,模拟效果如下:
作者的研究证明:通过 MLS 在 VUMAT 里计算 GND,可以在 ABAQUS 中完整重现微成形的尺寸效应,并清晰揭示“GND 在晶界和局部剪切带聚集”是强化的主要来源,同时保证数值方法可扩展、可工程化。详细的数值实现策略可以参考原始文献。
使用文章提到的策略,尝试进行数值显示,首先在umat隐式中进行实现,并在后续中修改为vumat即可。实现策略验证使用包含200个晶粒的二维模型拉伸验证。分别使用CPE6单元(二维多积分点,使用传统的GND计算方案),CPE3单元(mesh-free策略),模型共包含27119,SSD计算使用经典的KM模型,流动方程使用唯象的幂律模型,取向随机分配给不同晶粒
初始多晶模型和网格如下:
拉伸变形10%后应力分布:
传统方案:
MLS方案:
拉伸变形10%后累计剪切分布:

传统方案:

MLS方案:

拉伸变形10%后SSD分布:

传统方案:

MLS方案:

拉伸变形10%后GND分布:

传统方案:

MLS方案:
总的计算时间
传统方案:60分钟15秒

MLS方案:36分钟12秒

可以看到MLS方案的计算策略计算效率通常会显著高于传统的GND计算策略,并且数值实现显然更适合在vumat框架中(计算效率高,单元类型适应性强),而且可以抑制局部网格噪声,因此在大变形,接触分析,成型计算中是一种非常合适的选择。感兴趣的可以进行数值尝试,也可以加入知识星球讨论交流

来源:我的博士日记
AbaqusUM材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-28
最近编辑:11月前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 134文章 164课程 0
点赞
收藏
作者推荐

圆锥曲线和天体运动,以及微分方程

圆锥曲线和天体运动,以及二阶线性偏微分方程什么是圆锥曲线 如图所示,圆锥曲线,是由一个圆锥被一个平面切割所产生的各种曲线。 圆锥是由一条与轴线保持固定夹角 α 的直线绕该轴线旋转一周而形成的立体图形。 当用不同角度的平面去切割这个圆锥,横截面与圆锥轴线之间的夹角 β,直接决定了我们得到哪种曲:当 0° ≤ β < α 时:平面与圆锥的两叶都相交,其交线是一条开放且有两支的双曲线(如图中橙色面所示)。当 β = α 时:平面恰好与圆锥的一条母线平行,此时得到的是一条开放的抛物线(如图中蓝色面所示)。当 α < β ≤ 90° 时:平面只与圆锥的一叶相交,交线是一个封闭的椭圆(如图中绿色面所示)。当 β = 90° 时:这是一个非常特殊的情况,截面与圆锥的轴线垂直,此时得到的便是一个完美的圆(如图中红色面所示)。我们可以将圆视为椭圆的一个特例。因此,双曲线、抛物线、椭圆和圆,这四种看似形态各异的曲线,实际上本质是统一的,它们都是平面与圆锥相交的所产生的曲线,统称为“圆锥曲线”。关于x和y的一般二元二次方程一般形式可以写作: 将该方程写成矩阵形式: 即: 其中 是实对称矩阵,他的行列式为: 将 作为判别式当 时,方程变为椭圆形(退化情况有点椭圆虚椭圆)当 时,方程变为抛物线(退化情况有两条平行直线)当时,方程变为双曲线(退化情况有两条相交直线)天体运行物理规律当天体的动能不足以克服引力源的引力势能时(总机械能<0)天体在一个封闭的,周期性的轨道上做椭圆运动(条件极其苛刻的情况下做圆周运动几乎不存在)例如地球围绕太阳,月球围绕地球,哈雷彗星光顾地球当天体的动能恰好够克服引力源的引力势能时(总机械能=0)由于条件极其苛刻因此几乎不存在,但是在地球上理想情况下的抛物运动总是抛物线轨迹当天体的动能大于引力源的引力势能时(总机械能>0)天体从某个方向飞来,然后在中心天体的引力作用下发生偏转(引力弹弓效应)然后朝另一个方向飞去做双曲线运动例如星际彗星‘奥陌陌’,最近的‘阿特拉斯’以及人类发射的星际探测器二阶线性偏微分方程一般形式可以写作: 类似于一般的二元二次方程,其判别式: 当 时,对应的是一个椭圆方程,他是封闭、有限的图形。椭圆上没有特定的“方向性”,它是匀称的。 最典型的例子是拉普拉斯方程 椭圆型方程没有实的特征线。这意味着信息传播速度是无穷大的。其描述的是平衡态、稳态问题,如稳定的温度场、静电场。某一点的值受区域内所有点的影响,同时也影响所有点,没有时间维度和因果关系。这类似于椭圆的“封闭性”和“整体性”。当 时,对应的是一个抛物线方程,他是开放、无限的,但它只有一个开放方向,具有独特的对称轴,指明了主导方向。 最典型的例子是热传导方程 有一族实的特征线描述的是扩散过程,如热传导、粘性流体的扩散。它描述的是随时间单向演化的系统。信息理论上可以瞬间传遍全域,但强度急剧衰减(比如在远处加热,近处只会感受到微弱的温度变化)。这类似于抛物线的“单向开放性”和“主导方向”。当时,对应的是一个双曲型方程,有两支,是开放、无限的,并且有两个明确的渐近方向。 最典型的例子是波动方程 有两族实的特征线描述的是波动现象,如声波、振动。信息以有限的速度 c沿特征线传播,具有清晰的因果律。某一点的扰动只影响其未来的一个锥形区域(影响域)。这完美地类比了双曲线的两支,代表了波前向前和向后传播的方向。总结圆锥曲线基于几何形状(不变量)分类,偏微分方程基于解的特性(不变量)分类,却具有可以类比的性质,相同的判别式,说明截然不同的物理现象,完全可能共享相同的数学本质后续将分享神经网络解微分方程的文章,本账号做分享的全部文章代码后续都会放在知识星球,如有需要请关注我吧~~~🤞🤞🤞来源:我的博士日记

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