PFC數(shù)值軟件基礎(chǔ)

1、PFCPFC2D2D顆粒流軟件培訓(xùn)顆粒流軟件培訓(xùn)中仿科技有限公司中仿科技有限公司主要內(nèi)容主要內(nèi)容第一部分第一部分 PFC2D顆粒流程序簡(jiǎn)介顆粒流程序簡(jiǎn)介 第二部分第二部分 有限差分法基礎(chǔ)介紹有限差分法基礎(chǔ)介紹第三部分第三部分 離散元法基礎(chǔ)介紹離散元法基礎(chǔ)介紹第四部分第四部分 PFC2D的使用的使用第一部分第一部分 PFCPFC2D2D顆粒流程序簡(jiǎn)介顆粒流程序簡(jiǎn)介1 1 、理論背景、理論背景2 2、顆粒流方法的基本假設(shè)、顆粒流方法的基本假設(shè)3 3、顆粒流方法的特點(diǎn)、顆粒流方法的特點(diǎn)4 4、可選特性、可選特性5 5、應(yīng)用領(lǐng)域、應(yīng)用領(lǐng)域6 6、求解步驟、求解步驟 作為離散元的一種，二維顆粒流程序作為

2、離散元的一種，二維顆粒流程序(Particle (Particle Follow Code PFC2D)Follow Code PFC2D)數(shù)值模擬新技術(shù)，其理論基礎(chǔ)是數(shù)值模擬新技術(shù)，其理論基礎(chǔ)是Cundall 1979Cundall 1979提出的離散單元法，用于顆粒材料力提出的離散單元法，用于顆粒材料力學(xué)性態(tài)分析，如顆粒團(tuán)粒體的穩(wěn)定、變形及本構(gòu)關(guān)系，學(xué)性態(tài)分析，如顆粒團(tuán)粒體的穩(wěn)定、變形及本構(gòu)關(guān)系，專門用于模擬固體力學(xué)大變形問題。它通過圓形專門用于模擬固體力學(xué)大變形問題。它通過圓形( (或異或異型型) )離散單元來模擬顆粒介質(zhì)的運(yùn)動(dòng)及其相互作用。由離散單元來模擬顆粒介質(zhì)的運(yùn)動(dòng)及其相互作用。由

3、平面內(nèi)的平動(dòng)和轉(zhuǎn)動(dòng)運(yùn)動(dòng)方程來確定每一時(shí)刻顆粒的平面內(nèi)的平動(dòng)和轉(zhuǎn)動(dòng)運(yùn)動(dòng)方程來確定每一時(shí)刻顆粒的位置和速度。作為研究顆粒介質(zhì)特性的一種工具，它位置和速度。作為研究顆粒介質(zhì)特性的一種工具，它采用有代表性的數(shù)百個(gè)至上萬個(gè)顆粒單元，通過數(shù)值采用有代表性的數(shù)百個(gè)至上萬個(gè)顆粒單元，通過數(shù)值模擬實(shí)驗(yàn)可以得到顆粒介質(zhì)本構(gòu)模型。模擬實(shí)驗(yàn)可以得到顆粒介質(zhì)本構(gòu)模型。1.1. PFC2D (Particle Follow Code 2 Dimension)PFC2D (Particle Follow Code 2 Dimension)即二維顆粒流程序，是通過離散單元方法來模擬圓即二維顆粒流程序，是通過離散單元方法來模擬

4、圓形顆粒介質(zhì)的運(yùn)動(dòng)及其相互作用。最初，這種方法形顆粒介質(zhì)的運(yùn)動(dòng)及其相互作用。最初，這種方法是研究顆粒介質(zhì)特性的一種工具，它采用數(shù)值方法是研究顆粒介質(zhì)特性的一種工具，它采用數(shù)值方法將物體分為有代表性的數(shù)百個(gè)顆粒單元，期望利用將物體分為有代表性的數(shù)百個(gè)顆粒單元，期望利用這種局部的模擬結(jié)果來研究邊值間題連續(xù)計(jì)算的本這種局部的模擬結(jié)果來研究邊值間題連續(xù)計(jì)算的本構(gòu)模型。以下兩種因素促使構(gòu)模型。以下兩種因素促使PFC2DPFC2D方法產(chǎn)生變革與發(fā)方法產(chǎn)生變革與發(fā)展展:(1):(1)通過現(xiàn)場(chǎng)實(shí)驗(yàn)來得到顆粒介質(zhì)本構(gòu)模型相當(dāng)通過現(xiàn)場(chǎng)實(shí)驗(yàn)來得到顆粒介質(zhì)本構(gòu)模型相當(dāng)困難困難:(2):(2)隨著微機(jī)功能的逐步增強(qiáng)，用

5、顆粒模型模隨著微機(jī)功能的逐步增強(qiáng)，用顆粒模型模擬整個(gè)問題成為可能，一些本構(gòu)特性可以在模型中擬整個(gè)問題成為可能，一些本構(gòu)特性可以在模型中自動(dòng)形成。因此，自動(dòng)形成。因此，PFC2DPFC2D便成為用來模擬固體力學(xué)和便成為用來模擬固體力學(xué)和顆粒流問題的一種有效手段。顆粒流問題的一種有效手段。 2 2、顆粒流方法的基本假設(shè)、顆粒流方法的基本假設(shè) 顆粒流方法在模擬過程中作了如下假設(shè)顆粒流方法在模擬過程中作了如下假設(shè): :1)1)顆粒單元為剛性體顆粒單元為剛性體; ;2)2)接觸發(fā)生在很小的范圍內(nèi)，即點(diǎn)接觸接觸發(fā)生在很小的范圍內(nèi)，即點(diǎn)接觸; ;3)3)接觸特性為柔性接觸，接觸處允許有一定的接觸特性為柔性接

6、觸，接觸處允許有一定的“重重疊疊”量量; ;4) “4) “重疊重疊”量的大小與接觸力有關(guān)，與顆粒大小相量的大小與接觸力有關(guān)，與顆粒大小相比，比，“重疊重疊”量很小量很小; ;5)5)接觸處有特殊的連接強(qiáng)度接觸處有特殊的連接強(qiáng)度; ;6)6)顆粒單元為圓盤形顆粒單元為圓盤形( (或球形或球形) )。 其中，顆粒為剛性體的假設(shè)，對(duì)于模擬介質(zhì)運(yùn)其中，顆粒為剛性體的假設(shè)，對(duì)于模擬介質(zhì)運(yùn)動(dòng)為只沿相互接觸面的表面發(fā)生的問題非常重要，動(dòng)為只沿相互接觸面的表面發(fā)生的問題非常重要，比如象砂土或糧食這種顆粒組合體材料，利用這比如象砂土或糧食這種顆粒組合體材料，利用這種假設(shè)在總體上來講是比較恰當(dāng)?shù)?，因?yàn)檫@種材種假

7、設(shè)在總體上來講是比較恰當(dāng)?shù)模驗(yàn)檫@種材料的變形是來自于顆粒剛性體間的滑動(dòng)和轉(zhuǎn)動(dòng)以料的變形是來自于顆粒剛性體間的滑動(dòng)和轉(zhuǎn)動(dòng)以及接觸面處的張開和閉鎖，而不是來自于每個(gè)剛及接觸面處的張開和閉鎖，而不是來自于每個(gè)剛性顆粒本身的變形，對(duì)于這種特殊材料，沒有必性顆粒本身的變形，對(duì)于這種特殊材料，沒有必要采用非常精確的數(shù)值模型，來得到對(duì)材料特性要采用非常精確的數(shù)值模型，來得到對(duì)材料特性的近似。的近似。3 3、顆粒流方法的特點(diǎn)、顆粒流方法的特點(diǎn) PFC2DPFC2D可以直接模擬圓形顆粒的運(yùn)動(dòng)和相互作用可以直接模擬圓形顆粒的運(yùn)動(dòng)和相互作用問題。顆料可以代表材料中的個(gè)別顆粒，例如砂粒，問題。顆料可以代表材料中的個(gè)

8、別顆粒，例如砂粒，也可以代表粘結(jié)在一起的固體材料，例如混凝土或也可以代表粘結(jié)在一起的固體材料，例如混凝土或巖石。當(dāng)粘結(jié)以漸進(jìn)的方式破壞時(shí)，它能夠破裂。巖石。當(dāng)粘結(jié)以漸進(jìn)的方式破壞時(shí)，它能夠破裂。粘結(jié)在一起的集合體可以是各向同性，也可以被分粘結(jié)在一起的集合體可以是各向同性，也可以被分成一些離散的區(qū)域或塊體。這類物理系統(tǒng)可以用處成一些離散的區(qū)域或塊體。這類物理系統(tǒng)可以用處理角狀塊體的離散單元程序理角狀塊體的離散單元程序UDECUDEC和和3DEC3DEC來模擬。來模擬。 PFC2DPFC2D有三個(gè)優(yōu)點(diǎn)：有三個(gè)優(yōu)點(diǎn)： 第一、它有潛在的高效率。因?yàn)閳A形物體間的接第一、它有潛在的高效率。因?yàn)閳A形物體間的

9、接觸探測(cè)比角狀物體間的更簡(jiǎn)單。觸探測(cè)比角狀物體間的更簡(jiǎn)單。 第二、對(duì)可以模擬的位移大小實(shí)質(zhì)上沒有限制。第二、對(duì)可以模擬的位移大小實(shí)質(zhì)上沒有限制。 第三、由于它們是由粘結(jié)的粒子組成，塊體可以第三、由于它們是由粘結(jié)的粒子組成，塊體可以破裂，不象破裂，不象UDECUDEC和和3DEC3DEC模擬的塊體不能破裂。模擬的塊體不能破裂。 用用PFC2DPFC2D模擬塊體化系統(tǒng)的缺點(diǎn)是，塊體的邊模擬塊體化系統(tǒng)的缺點(diǎn)是，塊體的邊界不是平的，用戶必須接受不平的邊界以換取界不是平的，用戶必須接受不平的邊界以換取PFC2DPFC2D提供的優(yōu)點(diǎn)。提供的優(yōu)點(diǎn)。 PFC2D PFC2D中幾何特征、物理特性和解題條件的說明

10、中幾何特征、物理特性和解題條件的說明不如不如FLACFLAC和和UDECUDEC程序那樣直截了當(dāng)。程序那樣直截了當(dāng)。 例如用連續(xù)介質(zhì)程序，創(chuàng)建網(wǎng)格、設(shè)置初始?jí)毫?、例如用連續(xù)介質(zhì)程序，創(chuàng)建網(wǎng)格、設(shè)置初始?jí)毫?、設(shè)置固定或自由邊界。在象設(shè)置固定或自由邊界。在象PFC2DPFC2D這樣的顆粒程序這樣的顆粒程序中，由于沒有唯一的方法在一個(gè)指定的空間內(nèi)組合中，由于沒有唯一的方法在一個(gè)指定的空間內(nèi)組合大量的粒子，粒子緊密結(jié)合的狀態(tài)一般不能預(yù)先指大量的粒子，粒子緊密結(jié)合的狀態(tài)一般不能預(yù)先指定。必須跟蹤類似于物體壓實(shí)的過程，直到獲得要定。必須跟蹤類似于物體壓實(shí)的過程，直到獲得要求的孔隙率。求的孔隙率。 由于顆粒

11、相對(duì)位置變化產(chǎn)生接觸力，初始應(yīng)力狀由于顆粒相對(duì)位置變化產(chǎn)生接觸力，初始應(yīng)力狀態(tài)的確定與初始?jí)好苡嘘P(guān)。由于邊界不是由平面組態(tài)的確定與初始?jí)好苡嘘P(guān)。由于邊界不是由平面組成，邊界條件的設(shè)定比連續(xù)介質(zhì)程序更復(fù)雜。成，邊界條件的設(shè)定比連續(xù)介質(zhì)程序更復(fù)雜。 當(dāng)要求滿足有實(shí)驗(yàn)室實(shí)際測(cè)試的模擬物體的力當(dāng)要求滿足有實(shí)驗(yàn)室實(shí)際測(cè)試的模擬物體的力學(xué)特性時(shí)，出現(xiàn)了更大的困難。在某種程度上，這學(xué)特性時(shí)，出現(xiàn)了更大的困難。在某種程度上，這是一個(gè)反復(fù)試驗(yàn)的過程，因?yàn)槟壳斑€沒有完善的理是一個(gè)反復(fù)試驗(yàn)的過程，因?yàn)槟壳斑€沒有完善的理論可以根據(jù)微觀特性來預(yù)見宏觀特性。論可以根據(jù)微觀特性來預(yù)見宏觀特性。 然而，給出一些準(zhǔn)則應(yīng)該有助于

12、模型與原型的然而，給出一些準(zhǔn)則應(yīng)該有助于模型與原型的匹配，如哪些因素對(duì)力學(xué)行為的某些方面產(chǎn)生影響，匹配，如哪些因素對(duì)力學(xué)行為的某些方面產(chǎn)生影響，哪些將不產(chǎn)生影響。應(yīng)該意識(shí)到，由于受現(xiàn)有知識(shí)哪些將不產(chǎn)生影響。應(yīng)該意識(shí)到，由于受現(xiàn)有知識(shí)的限制，這樣的模擬很難。然而，用的限制，這樣的模擬很難。然而，用PFC2DPFC2D進(jìn)行試進(jìn)行試驗(yàn)，對(duì)固體力學(xué)，特別是對(duì)斷裂力學(xué)和損傷力學(xué)，驗(yàn)，對(duì)固體力學(xué)，特別是對(duì)斷裂力學(xué)和損傷力學(xué)，可以獲得一些基本認(rèn)識(shí)?？梢垣@得一些基本認(rèn)識(shí)。 PFC2D PFC2D能模擬任意大小圓形粒子集合體的動(dòng)態(tài)力能模擬任意大小圓形粒子集合體的動(dòng)態(tài)力學(xué)行為。學(xué)行為。 粒子生成器根據(jù)粒子的指定分

13、布規(guī)律自動(dòng)概粒子生成器根據(jù)粒子的指定分布規(guī)律自動(dòng)概率地生成。粒子半徑按均勻分布或按高斯分布規(guī)率地生成。粒子半徑按均勻分布或按高斯分布規(guī)律分布。律分布。 初始孔隙度一般比較高，但通過控制粒子半初始孔隙度一般比較高，但通過控制粒子半徑的擴(kuò)大可以獲得密度壓實(shí)。在任何階段任何因徑的擴(kuò)大可以獲得密度壓實(shí)。在任何階段任何因素都可以改變半徑。所以不需反復(fù)試驗(yàn)就可以獲素都可以改變半徑。所以不需反復(fù)試驗(yàn)就可以獲得指定孔隙度的壓實(shí)狀態(tài)。得指定孔隙度的壓實(shí)狀態(tài)。 屬性與各個(gè)粒子或接觸有關(guān)，而不是與屬性與各個(gè)粒子或接觸有關(guān)，而不是與“類型號(hào)類型號(hào)”有關(guān)。有關(guān)。 因此，可以指定屬性和半徑的連續(xù)變化梯度。因此，可以指定屬

14、性和半徑的連續(xù)變化梯度?！肮?jié)理生成器節(jié)理生成器”用來修改沿指定軌跡線的接觸特性。用來修改沿指定軌跡線的接觸特性。假定這些線疊加在顆粒集合體上。用這種方法，模型假定這些線疊加在顆粒集合體上。用這種方法，模型可以被成組的弱面，如巖石節(jié)理切割?？梢员怀山M的弱面，如巖石節(jié)理切割。 粒子顏色也是一種屬性，用戶可以指定各種標(biāo)記方粒子顏色也是一種屬性，用戶可以指定各種標(biāo)記方案。案。 PFC2DPFC2D模型中為了保證數(shù)據(jù)長(zhǎng)期不漂移，用雙精度模型中為了保證數(shù)據(jù)長(zhǎng)期不漂移，用雙精度數(shù)據(jù)存儲(chǔ)坐標(biāo)和半徑。接觸的相對(duì)位移直接根據(jù)坐數(shù)據(jù)存儲(chǔ)坐標(biāo)和半徑。接觸的相對(duì)位移直接根據(jù)坐標(biāo)而不是位移增量計(jì)算。接觸性質(zhì)由下列單元組成

15、：標(biāo)而不是位移增量計(jì)算。接觸性質(zhì)由下列單元組成： 1 1）線性彈簧或簡(jiǎn)化的）線性彈簧或簡(jiǎn)化的Hertz-MindlinHertz-Mindlin準(zhǔn)則；準(zhǔn)則； 2 2）庫侖滑塊；）庫侖滑塊； 3 3）粘結(jié)類型：粘結(jié)接觸可承受拉力，粘結(jié)存在）粘結(jié)類型：粘結(jié)接觸可承受拉力，粘結(jié)存在有限的抗拉和抗剪強(qiáng)度。有限的抗拉和抗剪強(qiáng)度。 可設(shè)定兩種類型的粘結(jié)，接觸粘結(jié)和平行粘結(jié)?？稍O(shè)定兩種類型的粘結(jié)，接觸粘結(jié)和平行粘結(jié)。這兩種類型粘結(jié)對(duì)應(yīng)兩種可能的物理接觸：接觸這兩種類型粘結(jié)對(duì)應(yīng)兩種可能的物理接觸：接觸粘結(jié)再現(xiàn)了作用在接觸點(diǎn)一個(gè)很小區(qū)域上的附著作粘結(jié)再現(xiàn)了作用在接觸點(diǎn)一個(gè)很小區(qū)域上的附著作用；平行粘結(jié)再現(xiàn)了粒子

16、接觸后澆注其它材料的用；平行粘結(jié)再現(xiàn)了粒子接觸后澆注其它材料的作用（如水泥灌漿）。平行粘結(jié)中附加材料的有效作用（如水泥灌漿）。平行粘結(jié)中附加材料的有效剛度具有接觸點(diǎn)的剛度。剛度具有接觸點(diǎn)的剛度。 塊體邏輯支持附屬粒子組或塊體的創(chuàng)建，促進(jìn)了塊體邏輯支持附屬粒子組或塊體的創(chuàng)建，促進(jìn)了程序的推廣普及。塊體內(nèi)粒子可以任意程度的重疊，程序的推廣普及。塊體內(nèi)粒子可以任意程度的重疊，作為剛性體具有可變形邊界的每一個(gè)塊體，可作為一作為剛性體具有可變形邊界的每一個(gè)塊體，可作為一般形狀的超級(jí)粒子。通過指定墻的速度、混合的粒子般形狀的超級(jí)粒子。通過指定墻的速度、混合的粒子速度、施加外力和重力來給系統(tǒng)加載。速度、施加

17、外力和重力來給系統(tǒng)加載。“擴(kuò)展的擴(kuò)展的FISHFISH庫庫”提供了在集合體內(nèi)設(shè)置指定應(yīng)力場(chǎng)或施加應(yīng)力邊提供了在集合體內(nèi)設(shè)置指定應(yīng)力場(chǎng)或施加應(yīng)力邊界條件的函數(shù)。時(shí)步計(jì)算是自動(dòng)的，包括因?yàn)榻鐥l件的函數(shù)。時(shí)步計(jì)算是自動(dòng)的，包括因?yàn)镠ertzHertz接接觸模型剛度變化的影響。模擬過程中，根據(jù)每個(gè)粒子觸模型剛度變化的影響。模擬過程中，根據(jù)每個(gè)粒子周圍接觸數(shù)目和瞬間剛度值，時(shí)步也在變化?；诠乐車佑|數(shù)目和瞬間剛度值，時(shí)步也在變化。基于估計(jì)的粒子數(shù)，單元映射策略采用最佳的單元數(shù)目，自計(jì)的粒子數(shù)，單元映射策略采用最佳的單元數(shù)目，自動(dòng)調(diào)整單元的外部尺寸來適應(yīng)粒子缺失和指定的新對(duì)動(dòng)調(diào)整單元的外部尺寸來適應(yīng)粒子

18、缺失和指定的新對(duì)象。單元映射方案支持接觸探測(cè)算法以保證求解時(shí)間象。單元映射方案支持接觸探測(cè)算法以保證求解時(shí)間隨粒子數(shù)目線性增加，而不是二次方增加。隨粒子數(shù)目線性增加，而不是二次方增加。 類似于類似于FLACFLAC，PFCPFC提供了局部無粘性阻尼。這種提供了局部無粘性阻尼。這種阻尼形式有以下優(yōu)點(diǎn)：阻尼形式有以下優(yōu)點(diǎn)： 1 1）對(duì)于勻速運(yùn)動(dòng)，體力接近于零，只有加速運(yùn)）對(duì)于勻速運(yùn)動(dòng)，體力接近于零，只有加速運(yùn)動(dòng)時(shí)才有阻尼；動(dòng)時(shí)才有阻尼； 2 2）阻尼系數(shù)是無因次的；）阻尼系數(shù)是無因次的； 3 3）因阻尼系數(shù)不隨頻率變化，集合體中具有不）因阻尼系數(shù)不隨頻率變化，集合體中具有不同自然周期的區(qū)域被同等阻

19、尼，采用同樣的阻尼系同自然周期的區(qū)域被同等阻尼，采用同樣的阻尼系數(shù)。數(shù)。 PFC2DPFC2D可以在半靜態(tài)模式下運(yùn)行以保證迅速收斂可以在半靜態(tài)模式下運(yùn)行以保證迅速收斂到靜態(tài)解，或者在完全動(dòng)態(tài)模式下運(yùn)行。到靜態(tài)解，或者在完全動(dòng)態(tài)模式下運(yùn)行。 PFC2DPFC2D包含功能強(qiáng)大的內(nèi)嵌式程序語言包含功能強(qiáng)大的內(nèi)嵌式程序語言FISHFISH，允，允許用戶定義新的變量和函數(shù)使數(shù)值模型適合用戶的許用戶定義新的變量和函數(shù)使數(shù)值模型適合用戶的特殊需求。例如，用戶可以定義特殊材料的模型和特殊需求。例如，用戶可以定義特殊材料的模型和性質(zhì)、加載方式、實(shí)驗(yàn)條件的伺服控制、模擬的順性質(zhì)、加載方式、實(shí)驗(yàn)條件的伺服控制、模擬

20、的順序以及繪圖和打印用戶定義的變量等。序以及繪圖和打印用戶定義的變量等。 4 4、可選特性、可選特性 1 1）熱學(xué)分析）熱學(xué)分析2 2）并行處理技術(shù)）并行處理技術(shù)3 3）能寫用戶定義接觸模型）能寫用戶定義接觸模型4 4）用戶寫）用戶寫C+C+程序的程序的C+C+編程。編程。 熱學(xué)選項(xiàng)熱學(xué)選項(xiàng)用來模擬材料內(nèi)熱量的瞬間流動(dòng)和熱誘用來模擬材料內(nèi)熱量的瞬間流動(dòng)和熱誘導(dǎo)位移和力的順序發(fā)展。熱學(xué)模型可以獨(dú)立運(yùn)行或?qū)灰坪土Φ捻樞虬l(fā)展。熱學(xué)模型可以獨(dú)立運(yùn)行或耦合到力學(xué)模型。通過修改粒子半徑和平行粘結(jié)承耦合到力學(xué)模型。通過修改粒子半徑和平行粘結(jié)承受的力，產(chǎn)生熱應(yīng)變來解釋粒子和粘結(jié)材料的受熱。受的力，產(chǎn)生熱應(yīng)變

21、來解釋粒子和粘結(jié)材料的受熱。 用戶定義的接觸本構(gòu)模型用戶定義的接觸本構(gòu)模型可以用可以用C+C+語言來編寫，語言來編寫，并編譯成動(dòng)態(tài)鏈接庫文件，一旦需要就可以加載。并編譯成動(dòng)態(tài)鏈接庫文件，一旦需要就可以加載。 用戶寫的用戶寫的C+C+程序程序選項(xiàng)允許用戶用選項(xiàng)允許用戶用C+C+語言寫自己語言寫自己的程序，創(chuàng)建可執(zhí)行的的程序，創(chuàng)建可執(zhí)行的PFC2DPFC2D個(gè)人版本。這個(gè)選項(xiàng)可個(gè)人版本。這個(gè)選項(xiàng)可以用來代替以用來代替FISHFISH函數(shù)，大大提高運(yùn)行的速度。函數(shù)，大大提高運(yùn)行的速度。 并行處理技術(shù)并行處理技術(shù)允許將一個(gè)允許將一個(gè)PFC2DPFC2D模型分成幾個(gè)部模型分成幾個(gè)部分，每個(gè)部分可以在單獨(dú)

22、的處理器上平行運(yùn)行。與分，每個(gè)部分可以在單獨(dú)的處理器上平行運(yùn)行。與一個(gè)一個(gè)PFC2DPFC2D模型在一個(gè)處理器上運(yùn)行相比，平行處理模型在一個(gè)處理器上運(yùn)行相比，平行處理在內(nèi)存容量和計(jì)算速度方面得到大大提高。在內(nèi)存容量和計(jì)算速度方面得到大大提高。5 5、應(yīng)用領(lǐng)域、應(yīng)用領(lǐng)域 PFC2D PFC2D既可解決靜態(tài)問題也可解決動(dòng)態(tài)問題，既可解決靜態(tài)問題也可解決動(dòng)態(tài)問題，既可用于參數(shù)預(yù)測(cè)，也可用于在原始資料詳細(xì)情況既可用于參數(shù)預(yù)測(cè)，也可用于在原始資料詳細(xì)情況下的實(shí)際模擬。下的實(shí)際模擬。PFC2D PFC2D 模擬試驗(yàn)可以代替室內(nèi)試驗(yàn)。模擬試驗(yàn)可以代替室內(nèi)試驗(yàn)。在巖石與土體中開挖問題的研究與設(shè)計(jì)方面，實(shí)測(cè)在巖

23、石與土體中開挖問題的研究與設(shè)計(jì)方面，實(shí)測(cè)資料相對(duì)較少，關(guān)于初始應(yīng)力、不連續(xù)性等問題也資料相對(duì)較少，關(guān)于初始應(yīng)力、不連續(xù)性等問題也只能部分了解。而在松散介質(zhì)流動(dòng)問題中，影響流只能部分了解。而在松散介質(zhì)流動(dòng)問題中，影響流動(dòng)介質(zhì)不規(guī)律分布的影響因素很難定量描述。因此，動(dòng)介質(zhì)不規(guī)律分布的影響因素很難定量描述。因此，應(yīng)用應(yīng)用PFC2D PFC2D 初步研究影響整個(gè)系統(tǒng)的一些參數(shù)的特初步研究影響整個(gè)系統(tǒng)的一些參數(shù)的特性，對(duì)整個(gè)系統(tǒng)的特性有所了解后，就可以方便地性，對(duì)整個(gè)系統(tǒng)的特性有所了解后，就可以方便地設(shè)計(jì)模型模擬整個(gè)過程。設(shè)計(jì)模型模擬整個(gè)過程。 PFC2DPFC2D可以模擬顆粒間的相互作用問題、大變形問

24、可以模擬顆粒間的相互作用問題、大變形問題、斷裂問題等，適用于以下領(lǐng)域：題、斷裂問題等，適用于以下領(lǐng)域： （1 1）在槽、管、料斗、筒倉中松散物體的流動(dòng)問）在槽、管、料斗、筒倉中松散物體的流動(dòng)問題；題； （2 2）礦山冒落法開采中的巖體斷裂、坍塌、破碎）礦山冒落法開采中的巖體斷裂、坍塌、破碎和巖塊的流動(dòng)問題；和巖塊的流動(dòng)問題； （3 3）鑄模中粉料的壓實(shí)問題；）鑄模中粉料的壓實(shí)問題； （4 4）由粘結(jié)粒子組成物體的碰撞及其動(dòng)態(tài)破壞；）由粘結(jié)粒子組成物體的碰撞及其動(dòng)態(tài)破壞； （5 5）梁結(jié)構(gòu)的地震響應(yīng)及垮塌；）梁結(jié)構(gòu)的地震響應(yīng)及垮塌； （6 6）顆粒材料的基本特性研究，如屈服、流動(dòng)、）顆粒材料的基

25、本特性研究，如屈服、流動(dòng)、體積變化等；體積變化等； （7 7）固體的基本特性研究，如累積破壞、斷裂。）固體的基本特性研究，如累積破壞、斷裂。 6 6、求解步驟、求解步驟 1) 1)定義模擬對(duì)象定義模擬對(duì)象 根據(jù)模擬意圖定義模型的詳細(xì)程序，假如只對(duì)某根據(jù)模擬意圖定義模型的詳細(xì)程序，假如只對(duì)某一力學(xué)機(jī)制的不同解釋作出判斷時(shí)，可以建立一個(gè)一力學(xué)機(jī)制的不同解釋作出判斷時(shí)，可以建立一個(gè)比較粗略的模型，只要在模型中能體現(xiàn)要解釋的機(jī)比較粗略的模型，只要在模型中能體現(xiàn)要解釋的機(jī)制即可，對(duì)所模擬問題影響不大的特性可以忽略。制即可，對(duì)所模擬問題影響不大的特性可以忽略。2)2)建立力學(xué)模型的基本概念建立力學(xué)模型的基

26、本概念 首先對(duì)分析對(duì)象在一定初始特性形成初步概念。首先對(duì)分析對(duì)象在一定初始特性形成初步概念。為此，應(yīng)先提出一些問題，如系統(tǒng)是否將變?yōu)椴环€(wěn)為此，應(yīng)先提出一些問題，如系統(tǒng)是否將變?yōu)椴环€(wěn)定系統(tǒng)、問題變形的大小、主要力學(xué)特性是否非線定系統(tǒng)、問題變形的大小、主要力學(xué)特性是否非線性、是否需要定義介質(zhì)的不連續(xù)性、系統(tǒng)邊界是實(shí)性、是否需要定義介質(zhì)的不連續(xù)性、系統(tǒng)邊界是實(shí)際邊界還是無限邊界、系統(tǒng)結(jié)構(gòu)有無對(duì)稱性等。際邊界還是無限邊界、系統(tǒng)結(jié)構(gòu)有無對(duì)稱性等。 綜合以上內(nèi)容來描述模型的大致特征，包括顆綜合以上內(nèi)容來描述模型的大致特征，包括顆粒單元的設(shè)計(jì)、接觸類型的選擇、邊界條件的確定粒單元的設(shè)計(jì)、接觸類型的選擇、邊界

27、條件的確定以及初始平衡狀態(tài)的分析。以及初始平衡狀態(tài)的分析。 3) 3)構(gòu)造并運(yùn)行簡(jiǎn)化模型構(gòu)造并運(yùn)行簡(jiǎn)化模型 在建立實(shí)際工程模型之前，先構(gòu)造并運(yùn)行一在建立實(shí)際工程模型之前，先構(gòu)造并運(yùn)行一系列簡(jiǎn)化的測(cè)試模型，可以提高解題效率。通過系列簡(jiǎn)化的測(cè)試模型，可以提高解題效率。通過這種前期簡(jiǎn)化模型的運(yùn)行，可對(duì)力學(xué)系統(tǒng)的概念這種前期簡(jiǎn)化模型的運(yùn)行，可對(duì)力學(xué)系統(tǒng)的概念有更深入的了解，有時(shí)在分析簡(jiǎn)化模型的結(jié)果后有更深入的了解，有時(shí)在分析簡(jiǎn)化模型的結(jié)果后( (例如所選的接觸類型是否有代表性、邊界條件對(duì)例如所選的接觸類型是否有代表性、邊界條件對(duì)模型結(jié)果的影響程度等模型結(jié)果的影響程度等) )，還需將第二步加以修改。，還

28、需將第二步加以修改。 4) 4)補(bǔ)充模擬問題的數(shù)據(jù)資料補(bǔ)充模擬問題的數(shù)據(jù)資料 模擬實(shí)際工程問題需要大量簡(jiǎn)化模型運(yùn)行的模擬實(shí)際工程問題需要大量簡(jiǎn)化模型運(yùn)行的結(jié)果，對(duì)于地質(zhì)力學(xué)來說包括結(jié)果，對(duì)于地質(zhì)力學(xué)來說包括: : a) a)幾何特性，如地下開挖酮室的形狀、地形地幾何特性，如地下開挖酮室的形狀、地形地貌、壩體形狀、巖土結(jié)構(gòu)等；貌、壩體形狀、巖土結(jié)構(gòu)等； b)b)地質(zhì)構(gòu)造位置，如斷層、節(jié)理、層面等；地質(zhì)構(gòu)造位置，如斷層、節(jié)理、層面等； c)c)材料特性，如彈材料特性，如彈/ /塑性、后破壞特性等；塑性、后破壞特性等； d)d)初始條件，如原位應(yīng)力狀態(tài)、孔隙壓力、飽初始條件，如原位應(yīng)力狀態(tài)、孔隙壓力

29、、飽和度等；和度等； e)e)外荷載，如沖擊荷載、開挖應(yīng)力等。外荷載，如沖擊荷載、開挖應(yīng)力等。 因?yàn)橐恍?shí)際工程性質(zhì)的不確定性因?yàn)橐恍?shí)際工程性質(zhì)的不確定性( (特別是應(yīng)力特別是應(yīng)力狀態(tài)、變形和強(qiáng)度特性狀態(tài)、變形和強(qiáng)度特性) )，所以必須選擇合理的參，所以必須選擇合理的參數(shù)研究范圍。第三步簡(jiǎn)化模型的運(yùn)行有助于這項(xiàng)數(shù)研究范圍。第三步簡(jiǎn)化模型的運(yùn)行有助于這項(xiàng)選擇，從而為更進(jìn)一步的試驗(yàn)提供資料。選擇，從而為更進(jìn)一步的試驗(yàn)提供資料。 5) 5)模擬運(yùn)行的進(jìn)一步準(zhǔn)備模擬運(yùn)行的進(jìn)一步準(zhǔn)備 a)a)合理確定每一時(shí)步所需時(shí)間，若運(yùn)行時(shí)間過合理確定每一時(shí)步所需時(shí)間，若運(yùn)行時(shí)間過長(zhǎng)，很難得到有意義的結(jié)論，所以應(yīng)該

30、考慮在多臺(tái)長(zhǎng)，很難得到有意義的結(jié)論，所以應(yīng)該考慮在多臺(tái)計(jì)算機(jī)上同時(shí)運(yùn)行。計(jì)算機(jī)上同時(shí)運(yùn)行。 b)b)模型的運(yùn)行狀態(tài)應(yīng)及時(shí)保存，以便在后續(xù)運(yùn)模型的運(yùn)行狀態(tài)應(yīng)及時(shí)保存，以便在后續(xù)運(yùn)行中調(diào)用其結(jié)果。例如如果分析中有多次加卸荷過行中調(diào)用其結(jié)果。例如如果分析中有多次加卸荷過程，要能方便地退回到每一過程，并改變參數(shù)后可程，要能方便地退回到每一過程，并改變參數(shù)后可以繼續(xù)運(yùn)行。以繼續(xù)運(yùn)行。 c)c)在程序中應(yīng)設(shè)有足夠的監(jiān)控點(diǎn)在程序中應(yīng)設(shè)有足夠的監(jiān)控點(diǎn)( (如參數(shù)變化如參數(shù)變化 處、不平衡等處、不平衡等) )，對(duì)中間模擬結(jié)果隨時(shí)作出比較分，對(duì)中間模擬結(jié)果隨時(shí)作出比較分析，并分析顆粒流動(dòng)狀態(tài)。析，并分析顆粒流動(dòng)狀

31、態(tài)。 6) 6)運(yùn)行計(jì)算模型運(yùn)行計(jì)算模型 在模型正式運(yùn)行之前先運(yùn)行一些檢驗(yàn)?zāi)Ｐ?，然在模型正式運(yùn)行之前先運(yùn)行一些檢驗(yàn)?zāi)Ｐ?，然后暫停，根?jù)一些特性參數(shù)的試驗(yàn)或理論計(jì)算結(jié)果后暫停，根據(jù)一些特性參數(shù)的試驗(yàn)或理論計(jì)算結(jié)果來檢查模擬結(jié)果是否合理，當(dāng)確定模型運(yùn)行正確無來檢查模擬結(jié)果是否合理，當(dāng)確定模型運(yùn)行正確無誤時(shí)，連接所有的數(shù)據(jù)文件進(jìn)行計(jì)算。誤時(shí)，連接所有的數(shù)據(jù)文件進(jìn)行計(jì)算。 7)7)解釋結(jié)果解釋結(jié)果 計(jì)算結(jié)果與實(shí)測(cè)結(jié)果進(jìn)行分析比較。圖形應(yīng)集中計(jì)算結(jié)果與實(shí)測(cè)結(jié)果進(jìn)行分析比較。圖形應(yīng)集中反應(yīng)要分析的區(qū)域如應(yīng)力集中區(qū)，各種計(jì)算結(jié)果應(yīng)反應(yīng)要分析的區(qū)域如應(yīng)力集中區(qū)，各種計(jì)算結(jié)果應(yīng)能方便地輸出分析。能方便地輸出分析

32、。第二部分第二部分 有限差分法基礎(chǔ)介紹有限差分法基礎(chǔ)介紹 連續(xù)介質(zhì)三維快速拉格朗日有限差分計(jì)算方法連續(xù)介質(zhì)三維快速拉格朗日有限差分計(jì)算方法( FLAC( FLAC3D3D) ) 是近是近2020年來逐步成熟完善起來的一種新型年來逐步成熟完善起來的一種新型數(shù)值計(jì)算方法數(shù)值計(jì)算方法, ,它基于顯式差分法來求解運(yùn)動(dòng)方程和它基于顯式差分法來求解運(yùn)動(dòng)方程和動(dòng)力方程動(dòng)力方程, ,可模擬巖土或其他材料的三維力學(xué)行為?？赡M巖土或其他材料的三維力學(xué)行為。其求解時(shí)首先將計(jì)算區(qū)域離散化其求解時(shí)首先將計(jì)算區(qū)域離散化, ,分成若干三維單元分成若干三維單元, ,單元之間由節(jié)點(diǎn)聯(lián)結(jié)單元之間由節(jié)點(diǎn)聯(lián)結(jié), ,節(jié)點(diǎn)受荷載作用后

33、節(jié)點(diǎn)受荷載作用后, ,其平衡方其平衡方程程( (運(yùn)動(dòng)方程運(yùn)動(dòng)方程) ) 可以寫成時(shí)間步長(zhǎng)為可以寫成時(shí)間步長(zhǎng)為t t 的有限差分的有限差分形式形式, ,由于采用動(dòng)態(tài)應(yīng)力松弛顯式差分求解技術(shù)由于采用動(dòng)態(tài)應(yīng)力松弛顯式差分求解技術(shù), , 在在某一微小的時(shí)段內(nèi)某一微小的時(shí)段內(nèi), , 作用于該節(jié)點(diǎn)的荷載只對(duì)周圍作用于該節(jié)點(diǎn)的荷載只對(duì)周圍若干節(jié)點(diǎn)有影響。若干節(jié)點(diǎn)有影響。 根據(jù)單元節(jié)點(diǎn)的速度變化和時(shí)段根據(jù)單元節(jié)點(diǎn)的速度變化和時(shí)段t ,t ,可求出單可求出單元之間的相對(duì)位移元之間的相對(duì)位移, ,進(jìn)而求出單元應(yīng)變進(jìn)而求出單元應(yīng)變, ,利用單元材利用單元材料的本構(gòu)關(guān)系即可求出單元應(yīng)力。在此基礎(chǔ)上料的本構(gòu)關(guān)系即可求出

34、單元應(yīng)力。在此基礎(chǔ)上, ,求出求出單元之間的不平衡力單元之間的不平衡力, ,將此不平衡力重新作用到節(jié)點(diǎn)將此不平衡力重新作用到節(jié)點(diǎn)上上, ,再進(jìn)行下一步的迭代過程再進(jìn)行下一步的迭代過程, ,直到整個(gè)系統(tǒng)不平衡直到整個(gè)系統(tǒng)不平衡力足夠小或節(jié)點(diǎn)位移趨于平衡為止。力足夠小或節(jié)點(diǎn)位移趨于平衡為止。 FLACFLAC3D3D可以解決諸多的有限元程序難以模擬的復(fù)可以解決諸多的有限元程序難以模擬的復(fù)雜的工程問題，例如分布開挖、大變形、非線性及雜的工程問題，例如分布開挖、大變形、非線性及非穩(wěn)定系統(tǒng)（甚至大面積屈服非穩(wěn)定系統(tǒng)（甚至大面積屈服/ /失穩(wěn)或完全塌方）。失穩(wěn)或完全塌方）。第三部分第三部分 離散元法基礎(chǔ)介

35、紹離散元法基礎(chǔ)介紹 離散單元法是一種模擬非連續(xù)介質(zhì)的計(jì)算方法，離散單元法是一種模擬非連續(xù)介質(zhì)的計(jì)算方法，自自CundallCundall在在7070年代提出以來，在巖石力學(xué)、土力學(xué)、年代提出以來，在巖石力學(xué)、土力學(xué)、結(jié)構(gòu)分析等領(lǐng)域的數(shù)值模擬中得到廣泛應(yīng)用，是一種結(jié)構(gòu)分析等領(lǐng)域的數(shù)值模擬中得到廣泛應(yīng)用，是一種新興的非連續(xù)體分析方法。離散單元法允許單元間的新興的非連續(xù)體分析方法。離散單元法允許單元間的相對(duì)運(yùn)動(dòng)，不一定滿足位移連續(xù)和變形協(xié)調(diào)條件，計(jì)相對(duì)運(yùn)動(dòng)，不一定滿足位移連續(xù)和變形協(xié)調(diào)條件，計(jì)算速度快，所需存儲(chǔ)空間小，特別適用于節(jié)理巖體的算速度快，所需存儲(chǔ)空間小，特別適用于節(jié)理巖體的大位移，大變形分

36、析。大位移，大變形分析。 離散單元法自問世以來有了長(zhǎng)足的發(fā)展，已經(jīng)成離散單元法自問世以來有了長(zhǎng)足的發(fā)展，已經(jīng)成為解決巖石力學(xué)問題的一種重要的數(shù)值方法，因?yàn)楣榻鉀Q巖石力學(xué)問題的一種重要的數(shù)值方法，因?yàn)楣こ讨兴姷降膸r體其形態(tài)呈非連續(xù)結(jié)構(gòu)，所形成的巖程中所見到的巖體其形態(tài)呈非連續(xù)結(jié)構(gòu)，所形成的巖石塊體運(yùn)動(dòng)和受力情況多是幾乎或材料非線性問題，石塊體運(yùn)動(dòng)和受力情況多是幾乎或材料非線性問題，所以很難用解決連續(xù)介質(zhì)力學(xué)問題的有限單元法或邊所以很難用解決連續(xù)介質(zhì)力學(xué)問題的有限單元法或邊界單元法等。界單元法等。 數(shù)值方法來進(jìn)行求解，而離散單元法正是充分?jǐn)?shù)值方法來進(jìn)行求解，而離散單元法正是充分考慮到巖體結(jié)構(gòu)的

37、不連續(xù)性，適用于解決節(jié)理巖石考慮到巖體結(jié)構(gòu)的不連續(xù)性，適用于解決節(jié)理巖石力學(xué)問題。力學(xué)問題。 近年來，離散元法的應(yīng)用領(lǐng)域又?jǐn)U展到求解連近年來，離散元法的應(yīng)用領(lǐng)域又?jǐn)U展到求解連續(xù)介質(zhì)向非連續(xù)介質(zhì)轉(zhuǎn)化的力學(xué)問題。混凝土等脆續(xù)介質(zhì)向非連續(xù)介質(zhì)轉(zhuǎn)化的力學(xué)問題?；炷恋却嘈圆牧显跊_擊、侵徹等動(dòng)荷載作用下產(chǎn)生的損傷和性材料在沖擊、侵徹等動(dòng)荷載作用下產(chǎn)生的損傷和破壞，其實(shí)質(zhì)是力學(xué)模型從連續(xù)體到非連續(xù)體的轉(zhuǎn)破壞，其實(shí)質(zhì)是力學(xué)模型從連續(xù)體到非連續(xù)體的轉(zhuǎn)變過程。建立在傳統(tǒng)的連續(xù)介質(zhì)力學(xué)基礎(chǔ)上的有限變過程。建立在傳統(tǒng)的連續(xù)介質(zhì)力學(xué)基礎(chǔ)上的有限元法等數(shù)值計(jì)算方法難以直接用于計(jì)算和模擬材料元法等數(shù)值計(jì)算方法難以直接用于

38、計(jì)算和模擬材料具體的破壞形式和破壞的整個(gè)過程，而離散元法在具體的破壞形式和破壞的整個(gè)過程，而離散元法在這一方面顯示出巨大的生命力。這一方面顯示出巨大的生命力。第四部分第四部分 PFCPFC2D2D的使用的使用1.1.對(duì)對(duì)PFCPFC軟件的使用界面、菜單功能及作用進(jìn)行介軟件的使用界面、菜單功能及作用進(jìn)行介紹；紹；2.FISH2.FISH語言簡(jiǎn)介語言簡(jiǎn)介3.PFC2D3.PFC2D分析模型的生成方法分析模型的生成方法4.4.邊界條件的設(shè)置方法邊界條件的設(shè)置方法5.5.初始條件的設(shè)置初始條件的設(shè)置6.6.接觸本構(gòu)模型：接觸剛度模型滑動(dòng)模型連接觸本構(gòu)模型：接觸剛度模型滑動(dòng)模型連接模型接模型7.7.賦予

39、材料屬性：相關(guān)命令的使用方法介紹賦予材料屬性：相關(guān)命令的使用方法介紹第四部分第四部分 PFCPFC2D2D的使用的使用8.8.節(jié)理面的生成及屬性設(shè)置節(jié)理面的生成及屬性設(shè)置9.9.加載方法：主動(dòng)荷載和被動(dòng)荷載加載方法：主動(dòng)荷載和被動(dòng)荷載; ;10.10.求解過程：靜力求解、動(dòng)力求解求解過程：靜力求解、動(dòng)力求解; ;11.11.流體與熱分析簡(jiǎn)介流體與熱分析簡(jiǎn)介; ;12.12.介紹介紹PFCPFC2D2D軟件的用戶自定義本構(gòu)模塊的相關(guān)軟件的用戶自定義本構(gòu)模塊的相關(guān)功能、操作等功能、操作等; ;13.13.常用命令使用方法及相關(guān)的重要概念常用命令使用方法及相關(guān)的重要概念; ;14.14.講述講述PF

40、CPFC2D2D工程應(yīng)用的實(shí)例工程應(yīng)用的實(shí)例1.使用界面、菜單功能介紹2.FISH語言簡(jiǎn)介 利用FISH語言進(jìn)行編程，應(yīng)該首先編一些簡(jiǎn)單的函數(shù)，然后仔細(xì)檢查函數(shù)的功能，測(cè)試是否有錯(cuò)誤。如果沒有發(fā)現(xiàn)錯(cuò)誤，再逐漸增加其功能，增加一項(xiàng)功能檢查一下，直至發(fā)展到最后比較復(fù)雜的程序。這是因?yàn)殡m然FISH是一種編譯型語言，但它沒有自己獨(dú)立的編譯器，不象VC+或VB能夠?qū)崟r(shí)全面地檢查錯(cuò)誤，F(xiàn)ISH檢查錯(cuò)誤的能力很差，因此在使用他們到真實(shí)的應(yīng)用之前，一定要用一些簡(jiǎn)單的數(shù)據(jù)（假如可能的話）來檢查所有定義的函數(shù)。 FISH函數(shù)內(nèi)置于標(biāo)準(zhǔn)的Itasca軟件的數(shù)據(jù)文件中，函數(shù)的格式必須以DEFINE開始，以END結(jié)束。

41、函 數(shù)可以嵌套調(diào)用，但定義函數(shù)的次序沒有關(guān)系，只要在使用之前全部定義就行。由于FISH函數(shù)的編譯格式儲(chǔ)存在Itasca軟件的內(nèi)存中，因此可以用SAVE命令保存函數(shù)以及相關(guān)變量的當(dāng)前值。 FISH也可以用來改進(jìn)用戶寫的本構(gòu)模型，如例1： DEF abc abc=22*3+5 END Print abc 對(duì)上例子稍作改進(jìn)(例2)： new def abc hh=22 abc=hh*3+5 end1).稍有編程常識(shí)的人可以看出，執(zhí)行上面的例子（PRINT abc），其結(jié)果與例1相同：abc=71.在這個(gè)函數(shù)中，我們首先把22賦值給變量hh，然后把這個(gè)變量帶入abc的表達(dá)式中，因此二者的結(jié)果相同。2)

42、.FISH的執(zhí)行過程如下：當(dāng)在程序命令中使用一個(gè)FISH符號(hào)名時(shí)（例如執(zhí)行PRINT 符號(hào)名），如果符號(hào)名也是一個(gè)函數(shù)名，那么執(zhí)行這個(gè)函數(shù)（例如abc）；如果符號(hào)名不是函數(shù)名，那么使用符號(hào)目前的值(例如hh)。3).在輸入完例2的各行后，如果我們執(zhí)行命令：PRINT hh，此時(shí)hh=0,因?yàn)樵谶@個(gè)時(shí)候沒有執(zhí)行FISH函數(shù)，因此hh的初始值為0；我們接著執(zhí)行PRINT abc，結(jié)果顯示abc=71;再次執(zhí)行PRINT hh,此時(shí)結(jié)果為hh=22，這是因?yàn)槲覀兪紫冗\(yùn)行了abc函數(shù)，在這個(gè)過程中hh已被賦值。4). 下面的試驗(yàn)將進(jìn)一步解釋函數(shù)與變量之間的差別。注意：Itasca軟件的SET命令可以用

43、來設(shè)置任何用戶定義的FISH符號(hào)的值，與在FISH中使用的符號(hào)無關(guān)。下面的例3建立在例2的基礎(chǔ)之上，我們不使用NEW命令來清除內(nèi)存中的值，因?yàn)槲覀兿肜^續(xù)使用那些值: set abc=0 hh=0print hh print abcprint hh5).在這個(gè)例子中，我們首先把a(bǔ)bc和hh都賦值為0，由于hh是一個(gè)變量，第一個(gè)Print命令顯示當(dāng)前hh的值，hh=0;第二個(gè)Print命令由于abc是一個(gè)函數(shù)名，因此執(zhí)行abc函數(shù)，先前定義的abc=0不起作用，重新計(jì)算了hh和abc的值，因此第三個(gè)Print命令顯示的值是它在abc函數(shù)內(nèi)指定的值，即hh=22,例4是這個(gè)試驗(yàn)完整的命令。 現(xiàn)在我們總結(jié)一下：Itasca軟件的三個(gè)重要的命令PRINT,SET,HISTORY可以直接操作簡(jiǎn)單的FISH變量或函數(shù)。如下圖所示，其中var代表變量名或函數(shù)名。 HISTORY命令的