Vuamp子程序用于显示动力学分析中，目前经常用于金属轧制成型。

Vuamp基本功能：

(1)可以通过vuamp定义当前幅值函数的值；

(2)可以用于模型的逻辑控制，一般通过设置的传感器的值；

(3)可以预定义状态变量；

(4)可以计算幅值函数的倒数和积分。

Vuamp基本格式如下：

SUBROUTINE VUAMP(

* ampName, time, ampValueOld, dt, nprops, props,

nSvars,

* svars, lFlagsInfo, nSensor, sensorValues,

sensorNames,

*

jSensorLookUpTable,

* AmpValueNew,

* lFlagsDefine,

* AmpDerivative, AmpSecDerivative, AmpIncIntegral)

INCLUDE 'VABA_PARAM.INC'

C time indices

parameter

(iStepTime = 1,

* iTotalTime = 2,

* nTime = 2)

C flags passed in for information

parameter (iInitialization =

1,

* iRegularInc = 2,

* ikStep = 3,

* nFlagsInfo = 3)

C optional flags to be defined

parameter

(iComputeDeriv = 1,

* iComputeSecDeriv = 2,

* iComputeInteg = 3,

* iStopAnalysis = 4,

* iConcludeStep = 5,

* nFlagsDefine = 5)

dimension time(nTime), lFlagsInfo(nFlagsInfo),

* lFlagsDefine(nFlagsDefine),

* sensorValues(nSensor),

* props(nprops),

* sVars(nSvars)

character*80 sensorNames(nSensor)

character*80 ampName

dimension jSensorLookUpTable(*)

user coding to define

AmpValueNew, and

optionally lFlagsDefine, AmpDerivative,

AmpSecDerivative,

AmpIncIntegral

RETURN

END

前面是变量定义，后面需要根据需要定义AmpValueNew等变量的值。

下面对变量进行解释。

需要定义的变量值：

AmpValueNew：当前幅值的值。

可以更新的变量：

lFlagsDefine(iComputeDeriv)：默认条件下为0， Abaqus将自动计算幅值对时间的一阶导数， 如果设置为1， 则需要设置幅值对时间的一阶导数。

lFlagsDefine(iComputeSecDeriv)：默认条件下为0， Abaqus将自动计算幅值对时间的二阶导数， 如果设置为1， 则需要设置幅值对时间的二阶导数。

lFlagsDefine(iComputeInteg)：默认条件下为0， Abaqus将自动计算幅值的积分， 如果设置为1， 则需要设置幅值的增量积分。

lFlagsDefine(iStopAnalysis)：默认条件下为0， Abaqus在出现错误信息时，将继续进行计算;若设置为1， 在出现错误信息后， 将停止计算。

lFlagsDefine(iConcludeStep)：默认条件下为0， 在运算结束后自动停止计算;设置为1， 如果还有下一步命令， 将继续进行下一步计算。

Svars：一组随着幅值变化而改变的阵列， 这个阵列的变量值数量就是nsvars，这个数组需要用户自己定义。在vuamp程序调用过程中， 阵列的信息传递给包含当前增量起始值的子程序中， 大多数时候， 用户需要自己更新增量的最后一组数据。

AmpDerivative：当前幅值的一阶导数。

AmpSecDerivative：当前幅值的二阶导数。

AmpIncIntegral：当前幅值的增量积分。

用于子程序和主程序之间传递信息的变量:

ampName：用户自定义的幅值名。

time(iStepTime)：当前步数的时间。

time(iTotalTime)：截止当前的总时间。

ampValueOld：前一幅值时间。

Dt：当前时间增量。

Nprops：相关的属性数量。

props(nprops)：用户定义的幅值属性。

nSvars：用户自定义的状态变量数目。

lFlagsInfo：和vuamp子函数相关的整数矩阵。

lFlagsInfo(iInitialization)：如果子函数调用从初始第一步开始， 其值为1， 否则为0.

lFlagsInfo(iRegularInc)：子函数时间增量为常数， 其值为1， 否则为0。

nSensor：模型中总的传感器数量。

sensorValues：当前的前一时间， 传感器的矩阵数值。每一个传感器矩阵数值对应一个历史输出历史输出数据库需要自定一个传感器。

sensorNames：用户自定义数组的传感器名。每一个传感器名对应一个数值， 这个数值提供输出请求， 如果系统中有小写格式或者大写、小写混合格式，则要求每个传感器名必须转化成大写格式。

jSensorLookUpTable：传递给功能函数IVGETSENSORID和VGETSENSORVALUE的值。

功能函数释义：

(1)IVGETSENSORID：获取传感器sensor的ID，返回值是该传感器的ID。格式：

IVGETSENSORID(mySensorName, jSensorLookUpTable

)

(2)VGETSENSORVALUE：获取传感器sensor的值，返回值是该传感器的数值。格式：

VGETSENSORVALUE(mySensorName, jSensorLookUpTable, sensorValues )

其中变量：

mySensorName：大写，用户定义的字符，传感器名称。

jSensorLookUpTable：指针。

下面是Abaqus帮助文档上的一个例子。

后面的逻辑代码主要是用来定义新的幅值的值，也即是变量ampValueNew的值，在这个过程中也利用状态变量svar来帮助更新。