Vusdfld\vumat等显示动力学子程序获取单元编号的一种方法

在编写vusdfld\vumat等子程序的时候，有时需要知道单元编号，以便对单元直接赋值，不同于usdfld，umat子程序，他们的单元编号信息会传入子程序中。但是显示动力学对应的子程序传入的是积分点编号，并没有单元编号。可以有一种取巧的方法来实现积分点和单元编号的对应。具体思路如下：

以边长为50mm立方体为例，网格种子密度为5，共1000个单元，单元类型C3R8R。而在显示动力学中，使用单核运算时，abaqus会一次性传入136个积分点，1000个单元会传入8次，最后一次为48个积分点，可以通过write(*,*) ‘nblock’,nblock输出来看传入积分点的变化情况。因此可以在程序中设置一个全局变量，每8次循环后就归0，重新循环。利用状态变量statenew(k,1)来存储单元。

在后处理中显示单元编号，可以看出云图从1到1000，正好对应所有单元的个数，利用查询值功能，可以看到单元ID和状态变量的值是相同的。

知道单元编号之后就可以对单元进行赋值操作，比如读取一个随机材料场文件，然后将数据赋值给单元。还可以利用该功能，实现类似“生死单元”的功能。一开始就根据单元编号抑制一部分，在之后的分析步在激活。比如第一个分析步杀死一些单元，然后下一个分析步在激活，比如删除编号为501的单元。

可以看到编号501的状态变量为值为0，单元应力为0，说明该单元被删除，不参与运算过程。如果在场输出中把STATUS勾选上，可以视觉上显示出来单元删除的效果。结果如下图所示：

需要注意的是如果一个单元有多个积分点，则循环中需要跳过一部分积分点。该方法需要单核计算才能使用。

最后有子程序开发等相关需求，欢迎联系我们。