01
在Abaqus材料子程序(如UMAT、VUMAT、UMATHT等)开发中,经常需要获取积分点的空间位置,用于判断材料属性随位置的变化(例如梯度材料、焊接残余应力场),或者计算与坐标相关的边界条件。最直接的方式就是利用子程序自带的COORD数组。
以UMAT为例,其接口中包含一个COORD参数。该数组存放的是当前积分点的坐标值,维度与模型的空间维数一致(平面问题为2,三维问题为3)。使用时只需读取COORD(1)、COORD(2)、COORD(3)即可得到x、y、z坐标。
这里有一个关键细节:COORD返回的是“当前构形”坐标还是“初始构形”坐标?
它取决于你是否打开了几何非线性(NLGEOM选项)。当NLGEOM = ON时,COORD返回的是当前构形下的坐标,即在大变形、大转动分析中,每个增量步的积分点位置会随网格移动而更新。当NLGEOM = OFF(小变形假设)时,COORD返回的是初始构形坐标(即参考构形),整个分析过程中保持不变。
实操提示:如果你需要在大变形分析中追踪某材料点的原始位置(例如用于计算对数应变),应使用初始坐标。但初始坐标并不直接作为参数传入UMAT,可以通过公用块或状态变量在第一步时保存COORD的初始值。另外,在GETVRM实用子程序中也可以请求坐标输出,但UMAT自带的COORD是最高效的。记住:打开几何非线性后,每次调用UMAT时COORD都会更新,千万别想当然地认为它始终是初始坐标。

02
通过history output输出积分点坐标
Abaqus可以直接在历程变量history output中输出积分点坐标。直接在history output中勾选COORD选项,但是这里需要注意的是,Domain中的Set集合如果是node set,这里输出来的是节点坐标,当这里是element set的时候,输出来的才是积分点坐标。


免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删
武汉格发信息技术有限公司,格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求,再低成本合规性管理软件许可,帮助贵司提高软件投资回报率,为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks 等。