手册的中文整合和翻译

FISH语言

4FISH语言初学指南

4.1简介和概述

FISH是一种内嵌于PFC2D的编程语言，使顾客可以定义新的变量和函

数。这些函数可以扩大PFC2D日勺应用范围或增长顾客定义特性。例如，

可以绘制或打印新的变量，生成特殊粒子，将伺服控制用于数值试验，

指定性能日勺不均匀分布，以及自动进行参数研究。

FISH语言是针对那些想用Itasca软件做些事情日勺人开发的I,而用现存日勺

Itasca软件很难或者说不也许做到。它提供了一种内置语言，因此顾客可

以写自己需要的函数，而不是在原则代码中加入许多新的和专门的功能。

某些有用日勺FISH函数已经写出；这些函数日勺互相关联设置，称为扩张的

Fishtank,在PFC2D程序中给出（见FISHvolume第3节）。对于没有编

程经验的人，写某些简朴的函数或对既有的函数作某些简朴修改是可行

的。第4.2节为非程序员作了简介。不过FISH程序也可以变的复杂（在

任何编程语言中的代码都是这样）；更多详细资料参见FISHvolume第2

节。

与所有的编程任务同样，FISH函数应按一种渐进的方式编写，在编写愈

加复杂日勺代码之前每一步都要检查运行。FISH函数的错误检查不会比大

多数编译器少，因此所有的函数在用于实际应用之前，应进行简朴的数

据测试。

FISH程序简朴地嵌入一种正常的PFC2D数据文献一在DEFINE背面口勺

行是FISH函数；函数碰到END终止。FISH函数还可以调用其他函数。

定义函数口勺次序并不重要，只要它们都是先定义后使用（例如通过PFC2D

命令引用）。由于FISH函数的汇编形式是储存在PFC2D的存储空间，

SAVE命令保留FISH函数和有关变量口勺目前值。

FISH语言规则和固有功能的完整定义见FISHvolume第2节。其中包括

语法规则，数据类型，运算，变量和函数。所有FISH语言名称在FISH

volume第2节进行了描述，在命令和FISH参照概要中给出了这些名称

的I概要。

4.2新手指南和教程

本节是为已经运行PFC2D（至少是某些简朴问题）而没有使用FISH语

言的人准备日勺；假如没有编程经验。为了从例子中获得最大日勺好处，你

应当尝试用PFC2D运行这些例子。短的程序可以直接键入。在运行一种

例子之后，使用命令NEW“清除界面”，就可认为下一种例子作准备。

此外较长的程序可以先在文献中创立然后用CALL调用。

在PFC2D命令提醒背面输入例4.1中时各行，每输入一行便在其末尾按

venter〉键。（在此和背面时所有例子，第一行用了分号，这仅表达注释

不必输入。）

Example4.1Defining*aFISHfunction

;zn«ir>e:fishbl.DAT

defabc

abc=22*3+5

end

注意在输入第一行后（不是注释行），命令提醒符变为DEF>；当输入END

命令后，命令提醒符变为常用提醒符。假如你是在给PFC2D或FISH输

入行得到话，提醒符的变化会让你理解。一般状况下，DEFINE语句背

面的所有行都作为FISH函数定义日勺一部分（直到输入END语句）。不过

假如输入口勺行中有错误（例如将“二”输入为“+”），那么会使PFC2D提

醒符回到本来状态。在这种状况下，应输入NEW命令重新开始。由于输

入很轻易出错，因此FISH程序一般通过编辑器输入文献。它们被PFC2D

调用就像止常的PFC2D数据文献。这一过程将梢后描述；目前继续运行。

假如在输入下面行时没有错误，在PFC2D＞提醒符下，实现abc函数日勺值

（在前面例4.1中定义的），输入行：

printabc

将得到：

abc=71

通过定义符号abc（用DEFINE…END构造，加例4.1）,目前在使用PFC2D

的许多方面都波及到该构造。

例如PRINT命令使FISH函数的值显示出来；该值通过一系列的算术运

算：

abc=22*3+5

这是一种赋值语句。假如等号成立，等号右边的体现式值将赋给笔号左

边的变量。注意算术运算按照一般的法则；加、减、乘、除分别用符号+,

*,/表达。符号“人”表达指数。

我们目前输入一种稍微不一样日勺程序（用命令NEW清除旧的））：

Example4.2Usingavariable

;fnarr.e:fishb2.DAT

new

defabc

hh-22

abc=hh*3+5

end

这里我们引入变量hh,给它赋值为22,在背面的行用到。假如我们输入

命令PRINTabc,输出的成果与前面完全相似。不过我们这里用到两个

FISH符号；它们均有值，一种（abc）是作为函数。此外一种（hh）是

作为变量。它们的区别如下：

当用到FISH函数符号名时（例如在PRINT语句中），假如符号对应一种

函数，那么有关B勺函数将被执行。不过，假如符号并不是一种函数名称，

那么仅仅使用符号的目前值。

下面日勺试验可以协助澄清辨别变量和函数。在此试验之前，请注意，

PFC2DH勺SET命令可以用来设置任何顾客定义的FISH符号的值，符号

在FISH程序中的独立性已简介。目前，键入如下行不使用NEW命令，

由于我们但愿保留此前输入的程序。

Example4.3SET〃〃为variables

;fname:fishb3.DAT

setabc«Ohn=O

printhh

printabc

printhh

SET命令将abc和hh的值设为0。由于hh是一种变量，第一种PRINT

命令简朴地显示了hh口勺值为0o第二个PRINT命令使abc执行（由于

abc是函数名）；hh和abc时值被重新计算。因此，第三个PRINT命令

表明hh已经重置其本来口勺值。为了测试您与否已经理解，请输入稍作修

改的程序（例4.4）并指出为何显示的答案是不一样的。

Example4.4Testyourunderstandingoffunctionandvariablenames

;fname:fishb4.DAT

new

defabc

abc=hh*3+5

end

sethh=22

printabc

setabc=0hh=O

printhh

printabc

printhh

在这个阶段，直接列出波及到FISH变量或函数的IPFC2D最重要的命令

是很有好处日勺。（见表4.1,var代表变量或函数名）

Table4.1Commandsthatdirectly

refertoFISHnames

PRINTvar

SETvar=value

HISTORYvar

我们已经看到了前面的两个例子（参照例4.3和4.4）；第三种有用日勺状

况是，当在历史变量的---------------------o例4.5阐明了怎样做到这

一点。

例4.5显示了粒子的不平衡力是怎样储存在FISH变量中和通过

HISTORY命令使用的。

Example4.5CapturingthehistotyofaFISHvariable

;fname:fishi>5.DAT

new

id5knle8nodes-5,05,0

callid1x0y1.1rad1

propdens1000knle8

setgrav0-10

sezdtmax=2e-3

oefdown_force

down_force=b_yfcb(ball_he£d)

PHf!

histdown_force

eye500

plothis1

在这个例子中，一种球落到平面上，经反弹后到达平衡。b_yfob在函数

中是一种预定义的变量名一在这种状况下，为对应的在y方向的不平衡

力。Ball_head是此外一种预定义名，为对应的球的）“地址”。PFC2D

所有预定义的名称在FISHvolume时表2.1中给出。运行

结束后，与其他历史同样，我们简朴MJ绘制了yforce（历史1,y方向的

不平衡力)。以类似的方式，可以用FISH函数绘制任何数量我们想得到

欧I历史，不管用多么复杂的公式来描述它。

除了上述预定义的变量名，尚有许多其他提供应FISH程序的预定义对

象。它们可以分为儿类；有一类是是由标量变量构成，它们有单一的号

码一例^口：

clockclocktimeinhundredthsofasecond

pi7i

Stepcurrentstepnumber

unbalmaximumunbalancedforce

urandrandomnumberdrawnfromunifbnndistributionbenveen

0.0and1.0.

内置于对象有用的另一类是一系列的固有功能，这些固有功能可以使诸

如正弦和余弦之类日勺函数能在FISH中程序计算。完整的清单见FISH

volume中第节；一部分如下：

abs⑷absolutevalueofa

COS(〃)cosineofa(aisinradians)

log(〃)base-tenlogarithmofa

max(〃,8)returnsmaximumofa.h

sqrt(〃)squarerootofa

使用固有功能的例子稍后给出。目前我们必须讨论一种更深入的措施,

运用这种措施，PFC2D数据文献可以运用顾客定义的FISH名。

当想在PFC2D中输入一种数字时，可以用FISH变量或函数来替代。

这个简朴的申明对于FISH非常强大的功能来说是很重要口勺，它容许在

FISH函数中计算，并使用通过PFC2D以符号形式输入的诸如范围，作

用力和性能等。因此，参数日勺变化可以很轻易，而无需在输入文献中变

化许多号码。

例如4.6规定将墙的几何形状和刚度定于文献开始的某个位置，而不是

许多位置。这就减少了出现错误和数据义献紊乱的风险，以及在许多模

型中反复运用一块工作。

Example4.6UsingsymbolicconstantsinaFISHfunction

;fnsir；e:fishb6.DAT

cefbox_georr.ez.ry

n_stiff-2e8

=le8

xx-5.0

yy=3.0

end

box_geometry

racrowalllnodesrnodes(0,0)(xx,0)，

racrcwa112nodes'nodes(xx,0)(xx,yy)'

racrowa±13nodesrnodesxxyyoyy'

racrowa114nodesznodes0yyoO'

wallid=lks=s_stizfkn=n_stizfwalllnodes

wallid-2ks-s_stiffkn-n_stiffwall2nodes

wallid=3ks=s_stizfkn=n_stizzwall3nodes

wallid=4ks-sstiffkn-nstiffwall4nodes

plotwallredidon

printwallprop

printmacro

顺便请注意，为FISH变量和函数选择名称有很大的灵活性；字符

可以包括在名称中。名称必须以非数字开头，并且不能包括任何算术运

算符（+,・，*或人）。选定H勺名称不应与内置（或保留）名称相似；FISH

volume第节中表2.1包括所有需要防止的名称的清单，以及应当遵照日勺

某些规则。

在上面日勺例子中，我们检查了FISH变量的计算值，通过将它们的名称直

接给PRINT命令。此外，我们可以通过命令printfish列出目前所有口勺

变量和函数。

我们目前研究怎样在FISH程序中进行选择和循环运行。下面的FISH语

句容许程序日勺指定部分反复多次。

LOOPvar（exprl9expr2）

ENDLOOP

LOOP和ENDLOOP是FISH语句，字符var代表loop变量，exprl和expr2

代表体现式（或单个变量）。例4.7显示了循环的应用（或反复次序）,

得到10日勺阶乘。

Example4.7ControlledloopinFISH

;fnane:fishb7.DAT

new

defxxx

sum-0

prod=1

loopn(I,10)

sum=sum+n

prod=prod*n

endloop

end

xxx

printsum,prod

在这种状况下，循环变量n从1持续变化到10,循环里面的语句（在LOOP

和ENDLOOP之间的语句）对每个价都执行。如前所述，变量名称或算

术体现式可取代1到10之间日勺数字。

在PFC2D中一种实际使用的循环，将按次序选用模型中的每个球并访问

和/或变化某些属性。例4.8包括了FISH函数makecolors,阐明了（随

机，在这种状况下）怎样变化模型中每个球的颜色指数。其中。到2内

欧I整数值指定为颜色指数。

Example4.8M(uiipulaiinf»variablesinaloopconstruct

;fname:fishb8.DAT

new

generatex-0,1y-0,1rad-0.05z0.07id-1,20

oefrr.akecolors

bp=ball_head

loopwhilebp¥null

xx=int(urand*3)

bcolor(bp)-xx

cp=b_next(bp)

endLOOD•

end

makecolors

plotballredgreenblue

我们已经看到了几种FISH程序的例子，让我们简要地检查一下程序语法

和风格的问题。一种完整的FISH语句必须占一行；没有持续行。假如一

种公式太长而不能放在一行，那么必须用一种临时变量将公式分开。例

4.9阐明了这个问题。

Example4.9Splittinglines

;fnar.e:fishb9.DAT

new

deflongsum/exampleofasumofmanythings

tempi-vl+v2+v3+v4+v5+v6+v7+v8+v9+vlO

longsum-tempi+vll+vl2+v!3+v14+vl5

prri

在这种状况下15个变量的和被分为两部分。注意分号之后的函数名表达

注释。分号之后的任何字符都会被FISH编译器忽视，不过它们反应了日

志文献。带有信息的注释程序是非常好的。有些程序出现了缩进一一在

某些行日勺开头插入空格来阐明有关口勺一组语句。可以在变量名和算术运

算之间插入（随意地）任何数量日勺空格字符使程序更具可读性。通过缩

进来表达循环，条件语句等是很好的程序实践。空格在FISH中的“重要”

意义，即空格字符不可以插入变量和函数名中。

需要阐明的另一种问题是变量的类型。您也许已经注意到，当从多种程

序的例子打印出变量时，这数字没有小数点或“E格式”一一指数日勺次

数。在任何时候，FISH函数或变量名分为四种类型：整数，浮点数，字

符串和指针。这些类型可以根据状况发生动态变化，但临时日勺顾客不必

紧张变量类型，由于它是自动设置的。思索例4.10。

Example4.10Variabletypes

;fname:fishblO.DAT

new

dezhaveone

aa-2

bb=3.4

cc='Haveaniceday7

dd=aa*bb

ee=cc,oldchap'

end

haveone

printfish

printaabbccddee

与变量aa,bb,cc相对应分派给它们日勺数字（或字符串）分别转化为整数,

浮点数和字符串。整型是精确H勺数字（没有小数点），但范围有限；浮点

数精度有限（大概精确到小数点后15位），范围较大；字符串变量是任

意的字符序列；在PFC2D中指针用来处理内部变量。I四种类型存在一定

转换规则。例如，dd成为一种浮点数，由于它是一种浮点数和一种整数

的乘积；ee变变成一种字符串，由于它是两个字符串的和（串联）。这个

问题可以相称复杂，它的I充足解释见FISHvolume第节。

此外，FISH中尚有常用的语言要素——if语句。下面的三个语句通过

FISH程序可以做出选择。

IFexprltestexpr2THEN

ELSE

ENDIF

这些语句容许FISH程序有条件H勺执行，ELSE和THEN表达选择。测试

的项目包括下列其中一种符号或符号配对：

=#><>=<=

Themeaningsarestandardexceptfor#,whichmeans“notequal/'expr1和

expr2是任意有效H勺体现或单一变量。假如测试为真，立即执行IF后H勺

语句直到碰到ELSE或ENDIFo假如测试为假，则执行ELSE和ENDIF

之间的语句；否则程序跳到ENDIF后的第一行。这些语句的应用见例

4.11o

Example4.11ActionoftheIFELSEENDIFconstruct

;fnane:fishbl1.DAT

new

dezabc

ifxx>0then

abc=33

else

abc-11

end_if

end

setxx-1

printsbe

setxx=-1

printabc

例4.1中显示的abc时值取决于xx日勺值。可以用不一样的测试符号检查

（例如用“v”替代

到目前为止，我们己经调用了来自PFC2D的FISH程序，通过使用PRINT

命令，或在PFC2D输入窗口独立行中给定函数名。也有也许是相反n勺一

-就是在FISH函数中给？？？？？？？？？？？？？？？？？？？？

最有效的PFC2D命令是位于下列两FISH语句之间：

COMMAND

ENDCOMMAND

从FISH程序中发出PFC2D命令重要有两个原因。首先，它可以使用FISH

函数执行我们已经讨论过的预定义变量不能实现的操作。第二，我们可

以控制FISH在PFC2D中日勺整个运行。

Example4.12UseoftheCOMMANDENDCOMMANDconstruct

;zname:fishbl2.DAT

new

PFC2DVersion3.1

defmake_walls

command

wa«lia=_nodes=(0,0)(5,0)

wallid=2noces=(5,0)(5,3)

wailid=3nodes=(5,3)(0,3)

wallid=4nodes=(0,3)(0,0)

endcommand

end

make_;valls

plotwallid=onred

在例4.12中，我们通过函数makewalls创立了四个墙。为了从FISH函

数中执行这四个PFC2D墙命令，它们必须在置于COMMAND

ENDCOMMAND构造之间。

我们目前已经覆盖到了FISH语言的某些方面以及怎样与PFC2D互相作

用。FISH语言口勺完整指南见FISHvolume第2节。

墙W、J状态信息功能

墙的状态信息包括位置，速度，非平衡加载和固定。

Table2.21Wallstaleinformationfunctions

typehuictionnamemodifiableavailability

F

W_PTRwp,INTdof)

FLTW_pos(

LTW_X(WJPTRwp)

F

LTW_y(WJPTRwp)

F

LTW_Z(WJPTRwp)3Donly

F

W_vvel(W_PTRup.INTdof)yes

FLT

LTw_xvel(WJPTRwp)yes

F

LTw_yvel(W_PTRup)yes

F

LTw_zvel(W_PTRup)yes3Donly

F

w_rvel(WJPTRup)yes2Donly

FLT

LTw_vrvel(W_PTRup.INTdof)yes3Donly

F

LTw_rxvel(WJPTRwp)yes3Donly

F

LTw_ryvel(W_PTRwp)yes3Donly

FLT

wrzvel(WPTRup)yes3Donly

F

W_vfob(W_PTRup.INTdof)

FLT

W_xfob(W_PTRwp)

FLT

LTw_yfob(W_PTRMR)

F

LTwzfob(WPTRup)3Donly

FLTw_mom(W_PTRwp)2Donly

FLTw_vmom(W_PTR［叩.INTdof)3Donly

FLTw_xmom(W_PTRwp)3Donly

FLTw_ymom(W_PTRup)3Donly

FLTw-zmom(W_PTRup)3Donly

F

w_radend1(W_PTR“p)yes3Donly

FLT

w_radend2(W_PTRwp)yes3Donly

FLT

LTw_radfob(W_PTRup)3Donly

F

LTw_radvel(W_PTRup)yes3Donly

INTW_flX(WJPTRup)yes2Donly

W_pos(wp,dof)旋转中心的I位置（dof为分量形式，其中

dofe{1,2,3})

W_x(wp)旋转中心的位置（x分量）

W_y(wp)旋转中心的位置（y分量）

W_z(wp)旋转中心出J位置（z分量）

W_vvel(wp,dof)速度（dof为分量形式，其中dofe{l,2,3}）

W_xvel(wp)速度（x分量）

W_yvel(wp)速度（y分量）

W_yvel(wp)速度（z分量）

W_rvel(wp)旋转中心周围的旋转速度

W_vrvel(wp,dof)旋转中心周围日勺旋转速度（dof为分量形式，

其中dofe{l,2,3（）

W_vxvel(wp)旋转中心周围日勺旋转速度（x分量）

W_vyvel(wp)旋转中心周围的旋转速度（y分量）

W_vzvel(wp)旋转中心周围日勺旋转速度（z分量）

W_vfob(wp,dof)考虑所有墙产生的非平衡力（dof为分量形式,

其中dofe（1,2,3））

W_xfob(wp)非平衡力（X分量）

W_yfob(wp)非平衡力（y分量）

W_zfob(wp)非平衡力（Z分量）

W_mom(wp)由所有发生接触口勺球产生的旋转中心口勺非平

衡力矩

W_vmom(wp,dof)旋转中心的非平衡力矩(dof为分量形式，

其中dofe{l,2,3})

W_xmom(wp)旋转中心的非平衡力矩(x分量)

W_ymom(wp)旋转中心的非平衡力矩(y分量)

W_zmom(wp)旋转中心的非平衡力矩(z分量)

W_radend1(wp)圆柱墙的lendl(endpointl)日勺半径

W_radend2(wp)圆柱墙E^Jend2(endpoint2)日勺半径

W_radfob(wp)圆柱墙日勺径向非平衡力

W_radvel(wp)expand或contract圆柱墙Ef、J径向速度

W_fix(wp)假如此值被设置成非零时，那么、墙节点速度

W、J更新是向内(inhibited)W、J。一般来说，这

个标志是在墙节点被顾客自定义一种功能的

时候直接使用。

2.5.1.5wall-sagment的功能

在PFC2D中，wall-sagment的J功能以wall-sagment数据构造

执行。由每个构成墙的Iwall-sagment被储存在一种连着的I序列

中，这个序列都可以由命令w_wlist(wp)得到。这个连着的序列

可以运用如下两个命令两面翻转(即可以以任一端为头)

ws_next(wsp)和ws_prev(wsp)o注意，由每个wall-sagment构

成内墙可以是封闭II勺，也可以是开放I付。

Table2.22Wall-segnientfunctions(2Donly)

typefimctionnamemodifiable

WS_PTRws_next(WSJPTRwsp)

WS_PTRws_prev(WS-PTRwsp)

FLTwsJengthfWS_PTRwsp)

FLTws_x(WS_PTRwsp)

FLTws_y(WS_PTRwsp)

FLTws_xun(WSJPTRyvsp)

FLTws_yun(WS_PTR)vsp)

FLTws_xvel(WSJPTRwsp)yes

FLTws_yvel(WS.PTRwsp)yes

Ws_next(wsp)在上述序列中下一种wall-sagmentH勺地址，(当

墙是不闭合的时候，那么在该序列中，最终一

种sagment的Jws_next为null。相反，假如墙是

闭合区I时候，那么序列中最终一种segment口勺

ws_next则为序列日勺第一种segment)

Ws_prev（vvsp）在上述序列中前一种wall-sagment的J地址：,（当

墙是不闭合的时候，那么在该序列中，第一

种sagment日勺ws_prev为nulL相反,假如墙

是闭合I向时候，那么序列中第一种segment

日勺ws_prev则为序列的I最终一种segmento）

Ws_length(wsp)wallsegmentvvsp口勺长度

Ws_x(wsp)与在上述序列中此segment日勺第一种节点对

应的wallsegmentwsp日勺end-0的I位置坐标（x

分量）。（当墙是由n个segment构成日勺开放环

时，则此时会有n+1个wallsegment数据构

造，并且最终一种节点日勺位置坐标会储存在

最终一种segment的数据构造中。）

Ws_x(wsp)wallsegmentwspend-0的位置坐标（y分

量)

Ws_xun(wsp)wallsegmentwsp区J单位法向向量（乂分量）。

（这个单位法向向量指向墙的I有效面。假如

墙是由n个segment构成的开放环时，则会

有n+1个wallsegment数据构造，并且最终

代Isegment的单位法向向量不会被定义。)

Ws_yun(wsp)单位向量(y分量)

Ws_xvel(wsp)速度(x分量)

Ws_yvel(wsp)速度(y分量)

.6wall-face欧J功能

wall-face日勺功能的I执行基于PFC3D有限wall-face构造。每

个有限长度的wall构成的Wall-face的I信息储存在一种连接

的链当中，其中数据链日勺头由wall-support功能中日勺

w_flist(wp)命令调用。这个数据链可以用wf_next(wp)功能

截断。

Table2.23Wall-facefunctions(3Donly)

typefimctionnamemodifiable

WF_PTRwf_next(WS-PTRwfp)

FLTwf_xun(W_PTRwfp)

FLTwf_yun(WJPTRwfp)

FLTwLzun(W_PTRwfp)

Wf_next(wfp)由墙构成出JfacewallH勺数据链中下一种

facewall的1地址。

Wf_xun(wfp)wallfacewfp的J单位法向量(x分量)。(这个

单位法向量指向强的有效面的反向。)

Wf_yun(wfp)wallfacewfpI内单位法向量(y分量)

Wf_zun(wfp)wallfacewfp的J单位法向量(z分量)

.7(parallel-bond)平行键的功能

平行键功能日勺执行是在PFC2D日勺平行键数据构造上。平行

键的信息与每一种接触有关，储存地址可以通过c_pb命令

获得。这个功能被分为两类，每一种均有它自己日勺小节：

物理性质(seetable2.24)和陈说信息(seetable2.25).

平行键H勺语法的描述见理论与背景中的节。

平行键物理性质功能

Table2.24Parallelbondphysicalproperlyfunctions

typefimctionnamemodifiable

FLTpb_kn(PB_PTRpbp)yes

FLTpb_ks(PBJPTRpbp)yes

FLTpb_nstrength(PBJPTRpbp)yes

FLTpb_sstrength(PBJPTRpbp)yes

FLTpb_rad(PB_PTRpbp)yes

Pb_kn(pbp)法向硬度（应力/位移）

Pb_ks(pbp)切向硬度（应力/位移）

Pb_nstrength(pbp)法向载荷（应力）

Pb_sstrength(pbp)切向载荷（应力）

Pb_rad(pbp)半径相加，因此平行键半径为Pb_rad

与两接触小球中最小半径口勺乘积。

Table2.25Parallelbondstaleinformationfimcliom

typetunctionnamemodifiableavailability

F

FLTpb_nforce(PB.PTRpbp)yes

LTpb_sforce(PB_PTRpbp)yes2Donly

FLT

pb_vsforce(PB.PTRpbp、INTdof)yes3Donly

FLT

LT

Fpb_xsforce(PB_PTRpbp)yes3Donly

LTpb_ysforce(PB_PTRpbp)

Fyes3Donly

pb^zsforcefPBJPTRpbp)yes3Donly

r

FLTpb_mom(PB_PTRpbp)yes2Donly

LTpb_tmom(PB_PTRpbp)yes3Donly

FLT

pb_vmom(PB_PTRpbp,INTdof)yes3Donly

FLT

LTpb_xbmom(PB_PTRpbp)yes3Donly

FLT

Fpb_ybmom(PB_PTRpbp)yes3Donly

pb^zbmom(PBPTRpbp)yes3Donly

FLTpb_nstress(PB.PTRpbp)yes

FLTpbsstress(PBPTRpbp)yes

Pb_nforce(pbp)键材料承载的法向力；正值代表压缩;

方向由c_vun（cp）给定

Pb_sforce（pbp）在接触面上键材料承载日勺切向力，接

触面由c_vun（cp）给定

Pb_vsforce（pbp,dof）W向力向量（dof为向量的I分量,dof

属于｛1,2,3))

Pb_xsforce(pbp)切向力向量（x分量）

Pb_ysforce(pbp)切向力向量（y分量）

Pb_zsforce(pbp)切向力向量（z分量）

Pb_mom(pbp)键材料承载的弯曲力矩

Pb_tmom(pbp)键材料承载日勺扭曲力矩（这个力矩就

在两个圆心的连线上）

Pb_vmom(pbp,dof)键材料承载时弯曲力矩向量（dof为

向量的分量，dof属于｛1,2,3））

Pb_xbmom(pbp)键材料承载的弯曲力矩向量（x分量）

Pb_ybmom(pbp)键材料承载H勺弯曲力矩向量（y分量）

Pb_zbmom(pbp)键材料承载的弯曲力矩向量（z分量）

Pb_nstress(pbp)作用在键外围的最大法向力

Pb_sstress(pbp)作用在键外围的最大切向力

.8measurement-circle的功能

measurement-circle欧J功能H勺执行基于

measurement-circle数据构造。这些功能被分为如下几

类：辅助(seetable2.26)和陈说信息(seetable2.27)。

measurement-circle的I信息被储存在一种首尾相连链

中，可以用circ—head命令调用他的头。储存地址也可

以由辅助命令find-meas见table2.5.

measurement-circle的I语法日勺描述见理论与背景中的3.4

节。

measurement-circle辅助功能

Table2.26Measurementcirclesupportfunctions

typefimctionnamemodifiableavailability

INTmJd(M_PTRmp)

M_PTRm_next(M_PTRmp)

FLTm_rad(MJPTRnip)yes

FLTm_vpos(MJPTRmp.INTdof)yes

FLTm_x(MJPTRnip)yes

FLTm_y(MJPTRnip)yes

FLTm_z(M_PTRmp)yes3Donly

M_id(mp)measurement-circlempR\JID号码

M_next(mp)数据链中下一种measurement-circle的I位

置。

M_rad(mp)measurement-circlemp日勺半径

measurement-circle陈说信息功能

注意到应力或张力率日勺目前值在measure()功能被执行了

后来才能使用，调和数(详细见下文)、空隙率和滑动比则

没有此规定。

Table2.27Measurementcirclestateinformationfunctions

typemnctionnameavailability

FLTm_coord(M_PTRinp)

FLTm_poros(M_PTRinp)

FLTmsfrac(MJPTRinp)

INTmeasure(MJPTRmp,INTcode)

F

FLTm_ed11(M_P