第9章ADAMS用户子程序-(陈立平)机械系统动力学分析及ADAMS应用

第9章ADAMS用户子程序本章对ADAMS用户子程序做了简要介绍，着重介绍了CONSUB、GFOSUB和REQSUB的使用方法，以及在用户子程序中两个最常用的功能子程序SYSARY和SYSFNC的使用情况。

通过本章的学习，读者将具备基本的开发用户子程序的能力。

9.1 ADAMS用户子程序简介一般情况下，ADAMS的大部分功能可以通过函数表达式完成，函数表达式很容易操作，因为不必编译或连接程序，而且ADAMS/Solver还会实时地提供函数表达式。

但函数表达式提供的只是有限的编程结构，因此有些复杂的情况，特别是涉及到一些逻辑表达，用函数表达式则很难表达出来。

因此在需要采用一些ADAMS没有提供的特殊函数时，可以采用用户子程序。

用户子程序更具有通用性，可以利用编程语言来定义模型元素或者特定的输出。

用户可以将函数表达式写成子程序的形式并将其与ADAMS/View连接，它具有函数表达式所没有的通用性和灵活性。

子程序利用通用程序设计语言（FORTRAN或C）的功能来定义ADAMS/View不能提供的函数，并使之按照需要而量身设计。

通过连接用户子程序，不会失去ADAMS/View的任何功效，也不会降低仿真速度。

当出现下列情况时，通常会使用到用户子程序：(1)数学函数很难表达。

(2)需要定义多用户使用的函数。

(3)GSE和UCON声明时，需要用户子程序。

(4)需要控制复杂仿真运行时，以及需要作决策逻辑时。

使用用户子程序时要小心，因为不正确的用户子程序会很难调试。

9.1.1用户子程序的种类根据用户子程序的功能，可以将它们分为三类：(1)Driver Subroutine——驱动子程序(2)Evaluation Subroutines——计算子程序(3)Restart subroutines——重启子程序1．驱动子程序（Driver Subroutine）驱动子程序（Driver Subroutine）－－CONSUB－－用来驱动ADAMS/Solver。

这样ADAMS/Solver在仿真时可更进一步地完成对模型的分析和修改的自动化工作。

CONSUB发出request(请求)，regulate(调整)，和vary simulations(改变仿真)的命令，并通过作出的交互式命令来修改数据，或者进行不同的分析。

大多数运行ADAMS/Solver的情况都是交互式，或者是以批处理的形式，一次只涉及到单个的分析。

使用CONSUB，可以完成对数据的动态分析，改变前提条件，并重复动态分析。

2．计算子程序（Evaluation Subroutines）计算子程序（Evaluation Subroutines）为用户自定义。

ADAMS/Solver可写子程序来对某些非标准的或复杂的计算式计算其数值。

计算子程序（Evaluation Subroutine）的程序名和功能如表9-1所示。

表9-1 计算子程序（Evaluation Subroutine）的程序名和功能3．重启子程序（Restart Subroutines）重启子程序（Restart subroutines）由一对子程序组成，SA VESUB和RELSUB。

它们可以保存和重新加载子程序运行过程中的某些状态变量，这些状态变量是那些相关的用户子程序在重新启动或进行某些变量初始化时所必须的数据。

如果在用户子程序之间相互调用时，需要使用它们在子程序内部保存数据。

如果没有使用RELSUB和SA VSUB，在保存仿真文件时，内部数据就会丢失。

虽然SA VE命令可将模型和仿真数据保存到文件中去，而且可以通过使用RELOAD命令重新调入保存文件，并可以从该保存点重新启动仿真。

但SA VE和RELOAD不可能对用户子程序的内部数据进行特殊操作，在RELOAD命令之后，ADAMS/Solver会从仿真流程中断点继续调用用户子程序，但用户子程序中某些感兴趣的内部变量却没有得到保存。

例如，假设SFOSUB子程序中力的系数是一个内部状态变量，如果打算将第一次调用SFOSUB子程序时所计算的力的系数保存以作今后使用，若用“SA VE”和“RELOAD”时会产生数据丢失。

可用SA VSUB将计算出来的系数保存在文件中并用RELSUB重新调入它们，或者力仅仅用RELSUB 来重新计算系数。

重启子程序（Restart Subroutine）的程序名和功能如表9-2所示。

表9-2 重启子程序（Restart Subroutine）的程序名和功能9.1.2 子程序的使用相对于ADAMS/View来讲，用户定义子程序有两个任务：输入和输出。

输入就是在相应的命令或状态，或在表明当前系统状态的计算结果中，提供给ADAMS/View的相关信息。

输出就是从ADAMS/View中获得的相关信息。

1. 声明浮点变量为了保证数值在子程序之间或者ADAMS/Solver和子程序之间能够正确地传递，需要声明浮点变量，就像在ADAMS/Solver中作的声明一样。

在UNIX和Windows机器上，ADAMS/Solver使用双精度浮点变量。

因为未声明的浮点变量通常都是单精度的(如REAL)，所以应当明确声明所有的浮点变量都是双精度的。

2. 参数的使用ADAMS/View会从FUNCTION=USER()的表达式中提供给用户子程序参数值。

例如，如果你使用FUNTION=(4.0，5.0，6.0)，ADAMS/View通过一个包含4.0，5.0，6.0的PAR数组参数，以及通过一个值为3的NPAR参数传递这些数值。

通过使用这些功能使得用户子程序更加灵活。

例如，如果计算需要marker101和marker102的距离，可以将这些标志符直接写进用户子程序。

3. 避免不连续在运行表达式中，当使用子程序来定义运动，传感器，力，或者模型中的状态平衡方程时，应当保证这些函数是连续的。

不连续的函数很难处理，因为大多数数学理论都是假定系统方程是连续的。

如果违反了该规则，ADAMS/Solver就会得出出人意料的结果。

当使用IF声明，将分段函数连接起来时，常常会产生不连续函数。

为了保证函数是连续的，在使用该类型的函数时，一定要小心。

4. 创建用户子程序为了创建自己的用户子程序，ADAMS提供了子程序模板并将它们放在安装目录install_dir/adams/solver/usersubs中。

子程序不能以C语言为模板。

如果需要一个C语言的子程序，必须通过修改FORTRAN文件来创建它，或者使用FORTRAN模板作为向导，创建自己的C语言格式子程序。

创建用户子程序的步骤：(1) 决定需要创建什么样的用户子程序。

例如，如果有GFORCE的声明，可以创建GFOSUB的用户子程序。

(2) 决定要计算什么，需要什么输入参数。

(3) 将合适的用户子程序模板从安装目录拷贝到当前目录。

模板目录中包含了要使用的所有子程序模板。

(4) 使用文字编辑器来察看用户子程序模板。

(5) 按需要修改用户子程序。

从用户子程序中调用功能子程序，或者其他FORTRAN或C语言子程序。

(6) 保存用户子程序。

因为可能对不同的模型会使用相同的子程序，所以最好给文件起一个容易区分的名字。

5. 调用用户子程序可以通过使用声明或者命令来调用用户子程序。

通过使用前面定义的用户子程序的名字来让ADAMS/Solver知道调用哪一个子程序。

也可以使用子程序来定义一些模型元素。

能调用用户子程序的声明或命令如表9-3所示。

表9-3 调用用户子程序的声明或命令可以从用户子程序中调用功能子程序。

这些子程序的功能就像是在用户的数据和用户使用的用户子程序的接口一样。

例如，在driver user——written subroutine(驱动用户子程序)中，可以使用control，access功能子程序。

从evaluation user——written subroutines(计算用户子程序)中，也可以调用access，function，以及某些情况下setup 函数。

换句话说，可以从上面任何一个用户子程序中调用access和function函数，但用户只能从driver subroutine中调用control函数，从某个evaluation subroutines中调用setup函数。

6. 编译，连接和运行库在写完定制的子程序后，必须要将它们编译并连接到ADAMS以创建到ADAMS/Solver的用户库文件中去。

然后，可以连同该库与ADAMS/Solver 一同运行。

必须有一个FORTRAN编译器，因为ADAMS本身不提供编译器，由用户子程序生成动态连接库的详细过程参见12.4.5节后半部分的介绍。

9．2 常用ADAMS用户子程序简介9.2.1 使用GFOSUB用户子程序实例为了使读者对用户子程序有一个直观印象，本节首先给出声明和调用GFOSUB子程序的例子，通过与GFORCE函数表达式的对比，说明使用用户子程序的必要性。

有关GFOSUB 用户子程序的更详细内容，将在9.2.2小节进行介绍。

另外，在12.4.4节的结尾部分提供了GFOSUB 用户子程序的具体工程应用实例。

现假定有如下一个GFORCE 函数表达式的声明，其函数表达式定义了小球和斜面之间的某种相互作用力关系（此处相互作用力不一定有确切的物理含义，仅为了说明GFORCE 的定义过程），包括力矢量和力矩矢量。

如下图所示，为了在GFORCE 中描述这种相互作用力关系，在小球球心上定义了标架Mark_I （ID 为2001），在斜面上定义了标架Mark_R （ID 为2003），同时为了确定小球对斜面的反作用力点的位置，生成与标架Mark_I 时刻重叠的标架Mark_J （ID 为2002），Mark_J 隶属于斜面。

于是可以使用GFORCE 函数表达式的声明，如下面所示（该GFORCE自身的ID 为100）：GFORCE/,100 I=2001, JFLOAT=2002, RM=2003,FX = － 10.0 * VX(2001，2003，2003) FY = － 10.0 * VY(2001，2003，2003) FZ = － 10.0 * VZ(2001，2003，2003) TX = － 100.0 * WX(2001，2003，2003) TY = － 100.0 * WY(2001，2003，2003) TZ = － 100.0 * WZ(2001，2003，2003)如果力和力矩的表达比较冗长，则可以使用GFOSUB 用户子程序来代替GFORCE 声明。

以及，用下面语句调用已创建好的GFOSUB 子程序,并向其传递参数。

GFORCE/,100 I=2001, JFLOAT=2002, RM=2003, FUNCTION=USER(10.0, 100.0, 2001, 2003, 2003)Part_I Mark_J ：：ID ＝2003Mark_I ：图9-1 小球与斜面的坐标定义及相对关系GFOSUB用户子程序的程序体创建如下：SUBROUTINE GFOSUB(ID，TIME，PAR，NPAR，DFLAG，& IFLAG，RESULT)======================= 定义与声明===================—————————————外部变量定义——————————INTEGER IDDOUBLE PRECISION TIMEDOUBLE PRECISION PAR( * )INTEGER NPARLOGICAL DFLAGLOGICAL IFLAGDOUBLE PRECISION RESULT(6)ID 标志号，在内部程序中调用TIME 仿真时间PAR 不定维的数组，用以向程序传递参数，本例中维数为5，分别为：PAR（1）＝10.0PAR（2）＝100.0PAR（3）＝2001PAR（4）＝2002PAR（5）＝2003NPAR PAR数组的维数，本例中NPAR ＝5DFLAG 积分标志位，IFLAG 初始标志位RESULT 返回的程序值，每一次调用可以返回六个值，本例中：RESULT（1）为FX的值RESULT（2）为FY的值RESULT（3）为FZ的值RESULT（4）为TX的值RESULT（5）为TY的值RESULT（6）为TZ的值———————————本地变量和参数定义—————————DOUBLE PRECISION FCT，TCR，VEL(6)INTEGER IPAR(3)，IM，JM，RM，NSTATESLOGICAL ERRFLG========================可执行代码=====================FCT = PAR(1)TCR = PAR(2)IM = PAR(3)JM = PAR(4)RM = PAR(5)……调用SYSARY子程序获得计算信息：如果变量IFLAG的值为true，则将子程序初始化。