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