Abaqus显示动力学载荷子程序VDLOAD的应用

1 VDLOAD子程序接口

   首先，还是去abaqus帮助文档找到VDLOAD的接口位置：

2 VDLOAD形参

找到之后，需要做的就是了解子程序形参所表示的含义，这一步的工作，只需要老老实实把帮助文档读一读就行了，其实很多子程序的形参都有共同的地方，大家一开始写的时候，可能会比较困难，等你写得多了，就会大彻大悟了，什么样的子程序拿到手，都能很快上手：

大家可以把这些形参读一读，我挑几个比较重要的讲：

（1）nBlock

   这是个变量，它有个值，这个值表示的是调用子程序的点的数量。什么意思呢，比如在我们的VDLOAD里。由于是节点调用子程序，所以这个值就表示调用vdload子程序的节点数量。而像VUMAT那种等等，由于它是单元积分点调用，这个变量，就表示的是单元积分点数量。那问题就来了，我们怎么知道是积分点还是节点还是别的呢。其实，它已经说的明明白白的：

   每个Explicit子程序都有这么一个变量。这个有什么用，大家要明白，与Standard不同的是，Explicit求解过程是大增量步去代替迭代的过程，所以你会看到Explicit的增量步都非常的多，如果每个增量步的每个积分点都单独调用子程序，那时间将会非常的长。本来explicit的求解效率就很低了，时间再长，根本顶不住。所以，每个增量步只调用一次，这一次就把所有的积分点处理了，于是乎，就有了子程序接口开始的这么一句话，有个循环。而这个循环，就表示，一次调用把所有积分点的数据都处理了。这就是它与Standard的区别。standard的每个积分点都是独立调用子程序的。

   搞清楚这一点之后，其实接下来的事，就和写Standard子程序，没有任何区别了。大家只要搞清楚子程序的形参之后。记住待定义的变量之后，就没有任何问题了。

（2）ndim

   这个好理解，如果你的模型是二维的话，这个值就是2，如果三维模型，这个值就是3

（3）curCoords(nblock,ndim)

   这个变量也是常用的，是个二维数组。储存节点的坐标

   剩下的形参，有些比较简单，有些不怎么用得到。如果有不理解的，把它写到文件里，然后看一下值是什么，有输入有输出，多练习，大家就明白是什么了。

3 待定义的变量（重要）

   没有什么比这个更重要了。这个是最最最最重要的，刻在心里。任何子程序，需要你定义的变量，都是非常重要的。

   这个变量就是在VDLOAD里更新的，就是我经常说的，要牢记在心里，这就是的目的。你写程序的目标，其他形参都是你的工具，你要用那些工具，去做各种各样的操作，然后得到你想要的值，给它。这就是我们写子程序的初心。就是为了定义这个值的。

4 模型

   说完子程序的三个要素之后，解释一下我们今天的模型。100*100*2mm的板，材料随机。大家注意一下我这里的坐标原点在哪，你们做的时候，保证和我的位置要一样，因为我写的子程序是依照这个坐标原点来的。

   step分析步建两个，时间都是0.001s，explicit的时间能给小一点就给小一点。算的快一点。

   载荷的话，给pressure。选user-defined，得启动一下vdload。下面全固定就行了。

5 写VDLOAD

   现在就是写了，比如我今天想实现一个，step-1力以一个圆向外扩散，然后step-2力再以一个圆返回中心。大家先记住，等会看了结果就明白了。

   算法的实现方式，当然是通过子程序来做了。因为我们这里有两个分析步，形参里却没有关系分析步的变量，所有只能通过时间来区别哪个是哪个分析步了。我们的两个分析步时间都是0.001s，所以，如果总时间再[0,0.001]之间的，就是Step-1，如果总时间在(0.001,0.002]之间，是Step-2。

6 结果

   把它丢进abaqus算一下，看一下结果，没问题啦：

   其他的杂七杂八的颜色，应该是显示动力学的应力波互相干涉产生的。

GIF

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删