PFC3d手冊(cè)命令中文詳解.docx

PFC3D 學(xué)習(xí)記錄5 COMMAND REFERENCE通用規(guī)定和特征命令語(yǔ)法默認(rèn)情況下，命令行忽略大小寫(xiě)，然而可以通過(guò)SET case命令激活區(qū)別大小寫(xiě)情況。斜體字代表數(shù)字，以i, j, m n 為開(kāi)頭的字母代表整數(shù)，其他情況為實(shí)數(shù)命令詞，關(guān)鍵字和數(shù)值可用空格，以及（），分隔。在分號(hào)（;） 后為注釋一行最多可包含80個(gè)字符。而&表示續(xù)行Range 語(yǔ)法定義一定限制范圍的物體可運(yùn)用Range的命令：Change/ initial/ propertyClumpFix/ FreeGroupJsetModelPlot ( 運(yùn)用于不同輸出項(xiàng)目)PrintRange( 創(chuàng)建一個(gè)命名的范圍)一個(gè)Range 定義一系列物體。由許多范圍元素構(gòu)成。如果使用多個(gè)范圍元素，最終物體幾何為不同元素集合的交集。但是可通過(guò)關(guān)鍵詞any改變選擇特性。如果使用關(guān)鍵詞any選項(xiàng)，任何包含在范圍內(nèi)的單元可被認(rèn)為最終物體集合的一部分。即求合集除此之外，通過(guò)在定義范圍單元后接關(guān)鍵詞not確定所有不在指定范圍的單元通常還可以用定義group 方式定義復(fù)雜的range。內(nèi)置范圍單元：Annulus center ( x, y, z) radius = r1, r2球形空心區(qū)域中心在（x, y, z ） 內(nèi)徑r1, 外徑r2Color il 球體或墻體在間隔（il, iu ）內(nèi)顏色索引。若無(wú)iu 則iu = il球體顏色索引列表通過(guò)property命令建立。墻體顏色索引列表必須用FISH 內(nèi)置的w_color設(shè)置Cylinder end1( x1, y1, z1 ) end2( x2, y2, z2 ) radius = r圓柱范圍。圓柱中心軸過(guò)end1( x1, y1, z1 ) end2( x2, y2, z2 ) 半徑為r.Fish fname調(diào)用用戶定義的范圍單元。通過(guò)FISH函數(shù)fname 執(zhí)行。而此函數(shù)中必須包含一個(gè)range_element指定范圍單元函數(shù)。此時(shí)fname就會(huì)調(diào)用所有在范圍內(nèi)的各個(gè)物體。物體的地址將通過(guò)fc_arg( 0 )傳入函數(shù)，而坐標(biāo)x-, y-, 和z-分量將通過(guò)fc_arg( 1 ), fc_arg( 2 ), fc_arg( 3 )分別傳入函數(shù)。物體類型由FISH內(nèi)置的pointer_type指定。如果物體不在指定范圍內(nèi)，則fname的值必須在函數(shù)內(nèi)設(shè)置為零，如果未指定，默認(rèn)值為。group gname定義屬于group名字gname的顆粒球體，（粒子組通過(guò)group命令定義）id il id編號(hào)在il，iu范圍內(nèi)的物體，如果iu沒(méi)有指定，則默認(rèn)il = iujset il 節(jié)理編號(hào)在在il，iu范圍內(nèi)的接觸，如果iu沒(méi)有指定，則默認(rèn)il = iu節(jié)理通過(guò)JSET命令創(chuàng)建name rname創(chuàng)建一個(gè)命名范圍rname，通過(guò)在命令行中范圍元素集合定義，通常一個(gè)范圍用range命令創(chuàng)建。plane keyword無(wú)限平面包圍區(qū)域Keyword:Origin x, y, z 平面上一點(diǎn)（x, y, z ）dd dd 方向角x-y平面上方向角，從y軸正向順時(shí)針?lè)较驕y(cè)量dip dip 俯角從x-y平面，向z軸負(fù)向。Distance d 在距離d范圍內(nèi)區(qū)域Above 平面上方below 平面下方normal nx, ny, nz 平面的單位法向矢量radius rl 半徑在（rl, ru ）范圍內(nèi)球體。ru未指定，則范圍（rl t, ru t ）t 未指定，則t = rl 10-6. 注意誤差范圍只有在未指定ru 條件下指定。rname 命名范圍rnamesphere center( x, y, z ) radius = r 球形區(qū)域，中心（x, y, z ） 半徑rx = xl, , y = yl, , z = zl, , 指定x, y, z 坐標(biāo)在指定范圍內(nèi)的物體。其他同radius 條重點(diǎn)理解fish定義range部分。在線幫助Command ? 適用command的各個(gè)關(guān)鍵詞？/ help 顯示所有可用命令命令分類程序控制命令這些命令支持批處理執(zhí)行模式，允許改變程序狀態(tài)。其與用程序建立實(shí)際數(shù)值模型無(wú)關(guān)。Cycle & Step 執(zhí)行指定時(shí)間步Solve執(zhí)行計(jì)算，直至滿足收斂條件Build 控制AC/DC 操作 call 調(diào)用data文件，執(zhí)行命令，批處理模式 return 從批處理模式中返回到控制臺(tái)?；蛟诙嗉?jí)調(diào)用中返回上一級(jí)調(diào)用文件continue- pause重新/ 中止調(diào)用data文件 Fluid設(shè)置流體分析條件 / MPI 控制MPI 操作/ PARALLEL控制并行計(jì)算程序操作NEW 重啟新模型Quit & Stop / Save Restore Set 改變?nèi)殖绦蛟O(shè)置/ Thermal 設(shè)置熱分析/ Title 設(shè)置模型名稱，這些將保存在隨后的輸出圖并記錄在save文件中/ Ucode注冊(cè)用戶編譯的C+ 類特殊計(jì)算模式可選計(jì)算模式，通過(guò)config命令調(diào)用。Config 可在任何分析階段給出，但是必須在調(diào)用時(shí)指定計(jì)算模式。Cppudm C+用戶定義接觸模型Cppuwc C+用戶編寫(xiě)代碼fluid 流體分析thermal 熱分析模型屬性部分模型建立和修改Ball 建立新球體Clump 建立新團(tuán)塊，或修改新團(tuán)塊屬性Delete 刪除球體，團(tuán)塊，墻體和時(shí)程記錄Generate 生成球體集合Jset 生成一系列節(jié)理Wall 生成墻體或修改墻體屬性Fix 對(duì)球體設(shè)置固定速度Free 將球體固定速度清零Model 對(duì)指定接觸指定用戶定義接觸模型Property 修改已有球體屬性（材料屬性和施加外力和速度），修改屬性，粘結(jié)屬性，接觸屬性模型監(jiān)控History 指定執(zhí)行過(guò)程變量（預(yù)定或用戶定義fish變量Measure 建立一個(gè)測(cè)量球或改變測(cè)量球?qū)傩訫ovie 對(duì)屏幕輸出進(jìn)行控制抓圖，從而實(shí)現(xiàn)動(dòng)畫(huà)輸出Plot 控制模型幾何和反應(yīng)屬性輸出Print 生成輸出變量和模型信息Trace 追蹤模型中能量及功的輸出命令工具Define 定義fish函數(shù)Group 定義和修改組集合Macro 命令宏Range 定義范圍Table 控制表的創(chuàng)建和修改。成對(duì)數(shù)值索引命令集合建立球體Ball radius = r Radius :球體半徑density : 密度（m/V）Hertz: 指定Hertz接觸模型否則為線性接觸模型Id: 編號(hào)，唯一性，默認(rèn)值大于當(dāng)前編號(hào)最大值的整數(shù)x, y, z ：x-, y- z-坐標(biāo)，默認(rèn)為0.0用ball 建立球體允許重疊，球體屬性由property命令設(shè)置Build 控制執(zhí)行AC/DC 操作Build assembly xa, ya, 通過(guò)復(fù)制粒子和包含在pbrick中的力（必須已經(jīng)存在）建立粒子集合，而初始粒子由Call 調(diào)用data文件，默認(rèn)為pfc3d.dat, 而若擴(kuò)展名未指定，等于*.dat文件調(diào)用允許多級(jí)調(diào)用，但不允許相互迭代。Change 同property命令一樣創(chuàng)建團(tuán)塊Clump full = on / off permanent release update_cycle = n創(chuàng)建一個(gè)團(tuán)塊，若id等于已有團(tuán)塊編號(hào)，則為修改團(tuán)塊屬性Range 為要建立團(tuán)塊的范圍，范圍內(nèi)原來(lái)屬于其他團(tuán)塊的粒子轉(zhuǎn)化到當(dāng)前團(tuán)塊，必須指定range防止所有粒子變成一個(gè)團(tuán)塊Print clump 輸出團(tuán)塊信息，plot clump 輸出團(tuán)塊1、Full on 在計(jì)算團(tuán)塊運(yùn)動(dòng)時(shí)，使用全部慣性向量力矩，這會(huì)耗費(fèi)更多計(jì)算時(shí)間。因?yàn)閏ross-coupling term使用，必須進(jìn)行矩陣轉(zhuǎn)換。當(dāng)旋轉(zhuǎn)速度較大和陀螺效應(yīng)影響大時(shí)使用。Full off 計(jì)算團(tuán)塊運(yùn)動(dòng)時(shí)，只有使用diagonal inertia，而忽略旋轉(zhuǎn)速度乘積。只有3D計(jì)算受full關(guān)鍵詞影響。當(dāng)full on 時(shí)，應(yīng)設(shè)置update_cycle = 12、Permanent 通過(guò)阻止構(gòu)成團(tuán)塊粒子間生成接觸使得構(gòu)成團(tuán)塊的粒子在隨后模擬過(guò)程保持接觸。這樣不需要存儲(chǔ)接觸信息，并且避免經(jīng)常性更新接觸信息及其對(duì)應(yīng)的矢量。當(dāng)在模擬過(guò)程中團(tuán)塊有個(gè)形成和分解過(guò)程，就不能指定這個(gè)關(guān)鍵詞然而，指定顆粒間永久接觸有利于節(jié)約內(nèi)存和計(jì)算時(shí)間。當(dāng)一個(gè)團(tuán)塊有其permanent屬性設(shè)置，其設(shè)置為永久的然而永久團(tuán)塊可被刪除且團(tuán)塊內(nèi)球體可被釋放。Release關(guān)鍵詞用于從團(tuán)塊提取球體，釋放時(shí)，球體將不存在接觸3、release 范圍內(nèi)所有屬于同一團(tuán)塊將被釋放。如果id未指定，所有范圍內(nèi)屬于任何團(tuán)塊的粒子將被釋放。團(tuán)塊內(nèi)所有粒子被釋放時(shí)，團(tuán)塊本身不被刪除，形成“空?qǐng)F(tuán)塊”, 這種情況在計(jì)算過(guò)程中忽略，而在粒子重新加入之后，被激活。Delete clump 刪除團(tuán)塊4、update_cycle n控制數(shù)據(jù)更新速度，默認(rèn)值n= 10。N越大，計(jì)算時(shí)間小?？杀桓滤俣扔蠱oment of inertia tensor 慣性張量力矩（3D）contact coordinates 接觸坐標(biāo)contact normal vector 接觸法向矢量接觸剪力矢量(3D) 平行粘結(jié)力矩（3D） 平行粘結(jié)剪力矢量(3D )注意團(tuán)塊阻尼系數(shù)為構(gòu)成團(tuán)塊粒子的阻尼系數(shù)的算術(shù)平均值當(dāng)模型中存在團(tuán)塊時(shí)，在本版本中就不能使用密度放大系數(shù)Fix/free命令用于團(tuán)塊內(nèi)任何粒子時(shí)，此工況用于團(tuán)塊重心。Print clump 可看到各個(gè)團(tuán)塊的fix/free狀況。Config keyword 指定可選計(jì)算模式Keyword:Cppudm C+用戶定義模型Cppuwc C+用戶編寫(xiě)代碼Fluid 流體分析Thermal 熱分析Continue 重啟由pause中止的data文件過(guò)程Cycle/Step n執(zhí)行n 時(shí)間步對(duì)于熱分析和熱力分析計(jì)算為熱計(jì)算步，否則為力學(xué)步計(jì)算過(guò)程輸入Esc，中止計(jì)算過(guò)程。計(jì)算過(guò)程輸入空格，繼續(xù)中止當(dāng)前cycle.當(dāng)前計(jì)算過(guò)程由set mechanical 和set thermal命令指定。對(duì)于熱力耦合計(jì)算，每個(gè)熱步后跟著最多n個(gè)力學(xué)計(jì)算子步，命令：thermal set substep.Damp keyword 設(shè)置局部阻尼的阻尼屬性（球體）和粘性阻尼的阻尼屬性（接觸） 局部阻尼：施加一個(gè)與各個(gè)球體自身不平衡力成比例大小的阻尼力。粘性阻尼：在每個(gè)接觸施加一個(gè)法向和切向粘壺，這些粘壺與原有接觸作用平行，從而生成與兩接觸實(shí)體（球球或球墻）相對(duì)速度差值成比力的粘滯阻力兩種阻尼可以同時(shí)存在當(dāng)未指定阻尼系數(shù)時(shí)，局部阻尼發(fā)生作用，且所有球體阻尼系數(shù)為0.7，而粘性阻尼不發(fā)生作用。Keyword1、 default 指定新建球的阻尼參數(shù)default local = v 局部阻尼系數(shù)默認(rèn)v= 0.7 viscous normal = v 法向臨界阻尼系數(shù)，默認(rèn)v = 0 notension off/ on 無(wú)拉力標(biāo)記on 無(wú)張力默認(rèn)off shear = v 切向臨界阻尼系數(shù)，默認(rèn)v = 02、 local = v 范圍內(nèi)所有球體局部阻尼系數(shù)為v3、 viscous normal = v notension = off/on shear = v設(shè)定范圍內(nèi)粘滯阻尼系數(shù)Define functionnameEnd定義fish函數(shù)所有define end 內(nèi)的語(yǔ)句被編譯并以壓縮方式存儲(chǔ)以便后續(xù)過(guò)程使用Delete keyword刪除球體，團(tuán)塊，時(shí)程關(guān)系，墻體Keyword:1、Balls 刪除id號(hào)球體 刪除質(zhì)心在范圍內(nèi)球體未指定任何參數(shù)時(shí)，模型內(nèi)所有球體被刪除2、Clump id刪除編號(hào)為id的團(tuán)塊，不能使用range參數(shù)。團(tuán)塊內(nèi)所有球體被釋放，但是球體未刪除。4、 histories 刪除所有時(shí)程記錄。時(shí)程記錄內(nèi)容可以用history reset命令刪除。Walls id1 刪除編號(hào)為id1,id2的墻體Id1 to id2 刪除編號(hào)在id1,id2間墻體Fix keyword 固定球體指定自由度速度值。即賦值速度常量。速度屬性可以由property命令設(shè)置。關(guān)鍵詞：x，y, z x-,y-,z-方向速度xspin, yspin, zspin : x-, y-, z- 方向角速度例1， 固定球體（id = 5 ）x-方向速度Fix x range id = 5Property xvel = 1.5 range id = 5Fluid keyword 流體分析關(guān)鍵詞:1、 model xl, xh, yl, yh, zl, zh size nx, ny, nz設(shè)置x-, y-, z-方向模型尺寸和單元數(shù)目。在xl, xh yl, yh zl, zh 范圍內(nèi)建立總數(shù)為nx* ny * nz數(shù)目流體單元，每個(gè)流體單元x-, y-,z-方向尺寸為（xh-xl）/ nx, (yh yl ) / ny, ( zh zl ) / nz，注意模型總單元數(shù)為（nx + 2）*( ny + 2) * (nz+2), 邊界上附加一個(gè)單元為邊界條件2、 property density d 流體密度d Viscosity 流體粘滯系數(shù)v3、boundary nonslip keyword 在keyword指定位置設(shè)置無(wú)滑動(dòng)邊界Slip keyword 在keyword指定位置設(shè)置滑動(dòng)邊界 Pressure pv keyword 在keyword指定位置設(shè)置壓力為pv Velocity xv, yv, zv keyword 在keyword指定位置指定apparent速度為xv, yv, zv注：apparent速度定義：對(duì)于進(jìn)入網(wǎng)格的流體將flow rate除以單元長(zhǎng)度（或單元面積）邊界位置指定：xl : x 方向下邊界xu: x 方向上邊界yl : y 方向下邊界yu: y 方向上邊界zl : z 方向下邊界zu: z 方向上邊界at xp,yp,zp 在點(diǎn)（xp, yp, zp ）4、set keyword 控制流體計(jì)算循環(huán)參數(shù)buoterm off( 默認(rèn) )/ on on: 打開(kāi)浮力作為外部體力（流體計(jì)劃中 grterm off(默認(rèn))/on off: 在Navier-Strokes方程中忽略重力on: 在流體計(jì)算中考慮重力影響 set gravity visterm off（默認(rèn)）/on off: 忽略粘滯影響下面參數(shù)控制SIMPLE項(xiàng)目，注意如果迭代不收斂，pfc3D發(fā)出一個(gè)警告消息而不是結(jié)束計(jì)算 Crit cr 指定SIMPLE 計(jì)劃收斂標(biāo)準(zhǔn)，默認(rèn)10-6 Intval intv 如果在tdel auto 時(shí)間步模式下指定。則流體計(jì)算將在每intv計(jì)算步調(diào)用，而流體計(jì)算步有模型狀況指定Iteration it 指定SIMPLE scheme迭代限制，默認(rèn)it = 5000 P_relax pr 在壓力調(diào)整計(jì)算中指定松弛因子默認(rèn)pr = 0.8Por_relax por 在孔隙率中指定松弛因子por = 1.0 Tdel auto 流體分析計(jì)算中時(shí)間步為自動(dòng)計(jì)算 Tdel tdel 指定流體時(shí)間步tdel V_relax vr 指定速度計(jì)算的松弛因子, 默認(rèn)vr= 0.5Free keyword 消除指定范圍內(nèi)指定自由度的固定速度約束關(guān)鍵詞：x/y/z/ xspin/yspin/zspinGenerate x = xl, xu y = yl, yu z = zl, zu radius = rl, ru id = il, iu 生成顆粒集合，指定粒子數(shù)目（編號(hào)起始和終止），粒子生成區(qū)域和粒子半徑分布。粒子位置和半徑以隨機(jī)分布。因此粒子集合最終狀態(tài)有隨機(jī)種子生成器（SET random ）影響。1、粒子生成區(qū)域定義：(1) x, y, z 關(guān)鍵詞決定長(zhǎng)方體區(qū)域（xl, xu ） ( yl, yu ) ( zl, zu )。當(dāng)xl = xu/ yl = yu / zl = zu 生成一排粒子集合。(2) annulus 關(guān)鍵詞定義一圓形空心區(qū)域中心（xc, yc, zc ）內(nèi)徑r1外徑r2。在這種情況下可忽略x, y, z 關(guān)鍵詞。如果同時(shí)指定。生成圓形空心區(qū)域和長(zhǎng)方體區(qū)域交集(3) 生成區(qū)域可用用戶定義的fish 函數(shù)定義（filter ）2、 半徑分布Radius 關(guān)鍵詞指定，默認(rèn)為均勻分布，高斯正態(tài)分布用 gauss 關(guān)鍵詞粒子數(shù)目粒子數(shù)目指定 id = il, iu 粒子數(shù)目 iu il + 1 。粒子生成時(shí)與已存在的粒子不能重疊。但空間不能足夠容納粒子，則可用tries關(guān)鍵詞。定義更多嘗試次數(shù)。如果不能生成足夠粒子，則可發(fā)生錯(cuò)誤提示，Set gen_error 用于將錯(cuò)誤提示改為警告提示，則命令將繼續(xù)。球體屬性將用隨后property 命令指定可選關(guān)鍵詞：1、 filter fname調(diào)用用戶定義生成粒子過(guò)濾器。在每個(gè)試驗(yàn)粒子生成過(guò)程中，F(xiàn)ISH 函數(shù)fname將被調(diào)用。試驗(yàn)粒子半徑將通過(guò)fc_arg( 0 ),而粒子坐標(biāo)三個(gè)分量將通過(guò)fc_arg( 1 ), fc_arg( 2 ), fc_arg( 3 )分別導(dǎo)入函數(shù)。如果符合條件（即粒子落入指定區(qū)域）則fname必須被設(shè)定為0，當(dāng)調(diào)用時(shí)fname設(shè)置為1。2、 gauss半徑為高斯分布，平均半徑（rl + ru）/ 2, 標(biāo)準(zhǔn)差為（ru - rl）/ 23. hertz新建球?yàn)镠ertz接觸模型，未指定為線性接觸模型4Local只有在并行計(jì)算時(shí)有用。通常只有在局部處理器，這樣就處理器之間沒(méi)有共享信息，必須保證沒(méi)有重復(fù)的id號(hào)5. Min rmin只有在gauss關(guān)鍵詞時(shí)定義，此時(shí)高斯分布中最小半徑為rmin, 默認(rèn)值為rl / 106. no_shadow 防止球體在墻體陰影中（即墻體從非活動(dòng)面的垂直投影）生成。默認(rèn)球體將在所有墻體活動(dòng)面和非活動(dòng)面任何一面生成。7. tries tmax設(shè)置粒子生成嘗試次數(shù)。默認(rèn)tmax = 20000。GROUP keyword定義群組，定義屬于球體的群組（不論位置和屬性改變）。可生成多級(jí)嵌套群組，比如一個(gè)球體可屬于多個(gè)群組Print group 顯示當(dāng)前存在群組Print ball info 顯示球體信息，包括其屬于的群組關(guān)鍵詞1. delete gname群組gname將刪除2. gname 建立群組gname. 群組將包括所有指定范圍內(nèi)球體。若無(wú)指定，將包括模型內(nèi)所有球體， 如果gname已經(jīng)存在，則位于指定范圍的球體將加入群組gname.Owner 關(guān)鍵詞， gname0被指定為gnames上一級(jí)（即gname作為群組gname0的從對(duì)象。如果上一級(jí)組gname0用range 單元定義，則gname0包含所有指定范圍內(nèi)球體和所有從組內(nèi)球體。而從組gname作為從組內(nèi)單元，則從組之包含從組內(nèi)球體。注意:range 定義中為疊加過(guò)程。History keyword 存儲(chǔ)和采樣一系列模型計(jì)算過(guò)程的變量。這些變量可以繪制與步數(shù)或其他相關(guān)時(shí)程變量（plot history）的圖形。時(shí)程變量可輸出到文件（write關(guān)鍵詞），每個(gè)History命令中只能有一個(gè)變量，history時(shí)程變量可在任何時(shí)間添加。所有時(shí)程變量?jī)?nèi)的具體內(nèi)容可用purge關(guān)鍵詞刪除而所有時(shí)程變量可用delete關(guān)鍵詞刪除Print history 顯示所有時(shí)程變量項(xiàng)目每個(gè)時(shí)程變量有唯一的ID號(hào)，默認(rèn)按history命令順序生成。可通過(guò)id關(guān)鍵詞指定。所有時(shí)程變量用一個(gè)采樣間隔。默認(rèn)采樣間隔為10步，采樣間隔可用nstep關(guān)鍵詞修改（更好的方法是set hist_rep命令指定） 不同時(shí)程變量不能使用不同的nstep。除此之外還可將FISH符號(hào)作為關(guān)鍵詞。關(guān)鍵詞1. 內(nèi)置時(shí)程變量（1）. 球體History ball 關(guān)鍵詞 x, y, z / id = n監(jiān)控特定球變量: 距離坐標(biāo)為（x, y, z ）最近的球/ id = n的球。關(guān)鍵詞有:s11：對(duì)應(yīng)應(yīng)力xx，s12 (xy), s13 (xz)。s33 Temperature: 溫度, 若熱學(xué)分析未激活，返回0 Xposition/ xspin/ xvelocity：質(zhì)心坐標(biāo),角速度，平移速度的x分量。（ y/z方向類似）（2） 診斷標(biāo)志History Diagnostic keyword mcf: 平均接觸力（所有非零法向力接觸的法向力的平均絕對(duì)值。 muf: 平均不平衡力（所有球體不平衡力平均絕對(duì)值 smr: 力學(xué)計(jì)算過(guò)程的比例限制 str: 熱學(xué)計(jì)算過(guò)程的比例限制（3）能量和功History Energy keyword ( 必須通過(guò)trace energy激活此命令) Body: 體力功 Bond: 粘結(jié)能Boundary: 邊界功（外力功）Frictional: 摩擦功Kinetic: 動(dòng)能Strain: 應(yīng)變能（3）測(cè)量球內(nèi)時(shí)程變量History Measure keyword id = id 編號(hào)為id 測(cè)量球內(nèi)數(shù)據(jù)變量 Coord 坐標(biāo)數(shù)目（每個(gè)球的平均接觸數(shù)目） ed11/ ed12 ed33 應(yīng)變率 1-x 2-y z-3 porosity 孔隙率 sliding_fraction: 滑動(dòng)比例（接觸中發(fā)生滑動(dòng)的比例） s11/s12/ s33 應(yīng)力（5）墻體History wall keyword id = idID = id 墻體的參數(shù) Power: 由于接觸球體所有力矩和外力施加在墻體上的功率 xforce/yforce/zforce: 由所有與之接觸球體產(chǎn)生不平衡力（x/y/z分量） xmoment/ymoment/zmoment: 由與之接觸球體作用在旋轉(zhuǎn)中心的不平衡力矩。2、時(shí)程英語(yǔ)語(yǔ)法支持（1） 刪除所有時(shí)程數(shù)據(jù)History delete 等效于：reset關(guān)鍵詞 + delete histories(2) 輸出與計(jì)算步對(duì)應(yīng)的編號(hào)為id1 到 idn 時(shí)程變量數(shù)據(jù)到屏幕History dump id1 關(guān)鍵詞用于指定輸出計(jì)算步范圍。 Begin nb: 所有計(jì)算步大于或等于nb的時(shí)程數(shù)值將被輸出。 End ne : 所有計(jì)算步小于或等于ne的時(shí)程數(shù)值將被輸出。 Skip ns : 每隔ns樣本采樣間隔的數(shù)據(jù)輸出 Vs id0 : 編號(hào)為id1 到idn 時(shí)程數(shù)據(jù)值與編號(hào)為id0s時(shí)程數(shù)據(jù)輸出 Xmaximum xmax: x分量值小于等于xmax的數(shù)據(jù)將被輸出，默認(rèn)x分量為計(jì)算步。然而如果vs id0關(guān)鍵詞，則x分量為id0的時(shí)程數(shù)據(jù)。 Xminimum xmin: x分量值大于等于xmin的數(shù)據(jù)將被輸出，默認(rèn)x分量為計(jì)算步。然而如果vs id0關(guān)鍵詞，則x分量為id0的時(shí)程數(shù)據(jù)。 ymaximum ymax: y分量值小于等于ymax的數(shù)據(jù)將被輸出 yminimum ymax: y分量值大于等于ymin的數(shù)據(jù)將被輸出(3) history limits: 輸出所有時(shí)程最大最小計(jì)算步數(shù)和數(shù)值(4). History print: 輸出所有時(shí)程數(shù)據(jù)信息。 print history(5). history purge: 所有時(shí)程項(xiàng)目?jī)?nèi)數(shù)據(jù)刪除，保存時(shí)程項(xiàng)目(6). History rename id = hname: 設(shè)置編號(hào)為hname的時(shí)程項(xiàng)目名稱，而時(shí)程項(xiàng)目的id號(hào)不變。(7). History reset: 刪除所有時(shí)程項(xiàng)目(8). History write id1 所有編號(hào)為id1 到idn 時(shí)程項(xiàng)目數(shù)據(jù)與計(jì)算步關(guān)系將被輸出至文件（file關(guān)鍵詞）或表格（table關(guān)鍵詞）關(guān)鍵詞中定義作用范圍：Begin nb 同前面End ne 同前面File 以ASCII格式輸出至文件fnameSkip ns 同前面Table n: 輸出結(jié)果至編號(hào)id =n 的表格中，如果未指定n, 則將創(chuàng)建。若指定n已存在，輸出數(shù)據(jù)將加入已存在表格末尾，注意表格使用方便性Vs id0 同上Xmaximum xmax 同上Xminimum xmin 同上Ymaximum ymax 同上Yminimum ymin 同上Initialize keyword 命令等同于 property命令JSET id = id keyword 生成節(jié)理，指定所有位于節(jié)理平面兩個(gè)面上的接觸為節(jié)理編號(hào)為id。節(jié)理平面可以為無(wú)限平面或包含許多有限圓盤(pán)。 只有粘結(jié)接觸或法向力非零的接觸受影響。Property命令（用jset范圍單元）用于指定鄰近指定節(jié)理的接觸屬性。關(guān)鍵詞有：Area_ratio m 節(jié)理平面被盤(pán)占用的比例，m為平均值，而d為標(biāo)準(zhǔn)偏差, m = 0 1 默認(rèn)m = 0, d = 0dd m 方位角（全局xy-平面從y軸順時(shí)針?lè)较颍?m、d同上dip m 俯沖角（從xy平面向下）m, d 同上gauss 對(duì)m,d服從高斯分布number n 節(jié)理面數(shù)目為norigin ( x, y, z ) 節(jié)理面原點(diǎn)或起始點(diǎn)radius m 圓形節(jié)理盤(pán)的半徑。 m, d 同上spacing m 節(jié)理間距. m, d 同上JSET命令必須在球體生成后使用，生成指定數(shù)目節(jié)理，在原點(diǎn)開(kāi)始，然后改變每邊，如果未指定radius和area_ratio. 節(jié)理面為無(wú)限平面。JSET中的range用于選擇要改變的接觸，從而創(chuàng)建有限節(jié)理平面。當(dāng)指定radius和area_ratio, 則沿著每隔節(jié)理平面生成許多圓盤(pán)。如果接觸投影在圓盤(pán)內(nèi)，則穿過(guò)節(jié)理平面的球間接觸屬性將改變。圓盤(pán)將隨機(jī)分布在節(jié)理平面的方形區(qū)域內(nèi)。圓盤(pán)生成直到所有圓盤(pán)面積占所有節(jié)理面面積比例達(dá)到area_ratio. 不考慮圓盤(pán)重合，故此參數(shù)為近似值。Macro string1 string2宏替換 string2為當(dāng)引號(hào)引用，嵌套時(shí)，仍然需要單引號(hào)。嵌套允許10級(jí)別，允許迭代，當(dāng)不能正確執(zhí)行。 不允許多行Measure id = id x = xc y = yc z = zc radius = r建立一測(cè)量球指定id, 若測(cè)量球id已經(jīng)存在，則指定測(cè)量球的屬性將被修改。若不指定，則比當(dāng)前最大值大1。Radius 不能為零Print measure 顯示測(cè)量球星系History measure 測(cè)量球時(shí)程變量。Model mname load filename指定范圍內(nèi)接觸為用戶定義接觸模型。如果沒(méi)有指定范圍，則為整個(gè)模型。當(dāng)使用此命令后，新接觸模型將被賦予內(nèi)置接觸模型。如果需要在新形成接觸上定義新的接觸模型，則必須通過(guò)合適的fishcall動(dòng)態(tài)賦值。關(guān)鍵詞load 為加載dll版本接觸模型，并且其用filename指定。一旦加載，則model命令range可加載dll接觸模型。而model load只有在config cppudm調(diào)用情況下可以加載。如下接觸模型可以使用，并且通過(guò)mname指定Ductile 簡(jiǎn)單的韌性材料與內(nèi)置線性接觸模型類似，他有用戶定義軟化過(guò)程，而不是象內(nèi)置材料的脆性破壞。Softening 通用位移軟化模型Viscous 簡(jiǎn)單粘彈性模型：其剪力特性中包含一個(gè)彈簧連著粘壺。下列接觸模型通過(guò)model load指定Burger 綜合Kelvin模型和Maxwell模型模擬蠕變： Model load burwrv.dllHysdamp: 通過(guò)摩擦滑動(dòng)將遲滯阻尼將能量消散引入線性接觸模型。 Model load hyswrv.dllVisdamp: 通過(guò)摩擦滑動(dòng)將粘滯阻尼將能量消散引入線性接觸模型。 Model load viswrv.dllMOVIE keyword 控制屏幕繪圖捕捉并且允許快速回放生成動(dòng)畫(huà)。在使用movie命令使用set plot命令控制文件格式和規(guī)格。如果未指定，將使用默認(rèn)avi格式文件并且使用默認(rèn)關(guān)鍵詞。關(guān)鍵詞有(1 ).Avi_open/ avi_close Avi_open 打開(kāi)一個(gè)新的avi文件，如果文件已存在覆蓋。未指定文件默認(rèn)文件為pfc3d.avi, 同時(shí)可打開(kāi)10個(gè)avi文件Avi_close： 關(guān)閉avi文件。其余同上注意：每次一個(gè)新的輸出窗口時(shí)，一個(gè)dcx格式文件自動(dòng)打開(kāi)和關(guān)閉。(2). Step n 編號(hào)為viewid輸出視口(未指定viewid時(shí)為當(dāng)前視口) 中內(nèi)容每n計(jì)算步存儲(chǔ)于file文件指定文件。默認(rèn)文件pfc3d.avi/ pfc3d.dcx.(3) snap 編號(hào)為viewid輸出視口(未指定viewid時(shí)為當(dāng)前視口) 中內(nèi)容以不規(guī)則或用戶定義間隔計(jì)算步存儲(chǔ)于file文件指定文件。默認(rèn)文件pfc3d.avi/ pfc3d.dcx.MPI keyword控制MPI并行操作。關(guān)鍵詞：(1). Root 指定當(dāng)前線程為根過(guò)程， 這個(gè)線程可執(zhí)行新的命令行，所有其他線程為從線程(2) transmit on使得根線程處理的命令（各種內(nèi)容：關(guān)鍵詞，F(xiàn)ISH代碼和data文件）發(fā)送至其他從線程。在從線程執(zhí)行命令避免其接受鍵盤(pán)輸入的命令，并且內(nèi)置call命令被忽略（ 因?yàn)楦€程讀取的內(nèi)容被傳遞至所有從線程）New清除所有程序狀態(tài)信息，在不離開(kāi)PFC3d情況下，允許開(kāi)始新任務(wù)。此時(shí)文件PFC3D.INI重新被所有起始命令調(diào)用。日志文件，回應(yīng)模式，以及隨機(jī)種子生成器（set log/ echo/ random ） 不受new命令影響。而其他程序狀態(tài)信息包括fish函數(shù)和變量，時(shí)程變量，表格以及輸出視圖，將清零。這些信息可通過(guò)save命令和restore命令實(shí)現(xiàn)保存和調(diào)用，另外fish函數(shù)可以寫(xiě)入一個(gè)文件保存。PARALLEL keyword 控制并行線程程序操作和初始說(shuō)明連接布局（topology ）關(guān)鍵詞：(1). end 當(dāng)此命令被主線程接受，所有線程將轉(zhuǎn)化為孤立模式(2). master mname在從線程中用此命令指定稱為mname網(wǎng)絡(luò)為其主線程(3). numprocess n主線程用此命令指定并行計(jì)算線程總數(shù)目為n(4). process pn每個(gè)線程用此命令給自身指定一個(gè)id。Pn范圍為 0, n-1 ，其中n為線程總數(shù)。而0為主線程保留(5). Start 對(duì)于主線程，此命令使得所有從鍵盤(pán)和文件輸入到PFC3D中的內(nèi)容不僅對(duì)主線程適用，同樣對(duì)從線程適用。對(duì)于從線程，此命令不允許所有鍵盤(pán)輸入，且命令所有PFC3D接受主線程特征輸入。FISH并行操作，以及FISHS生成的命令可在環(huán)境中解釋。多個(gè)其他命令使得數(shù)據(jù)交換在底層。Start命令是在假設(shè)process, numprocess和master已經(jīng)處理情況下運(yùn)行在并行計(jì)算中應(yīng)注意問(wèn)題：(1). 如果有n個(gè)處理器，編碼0，1，, n-1, 每個(gè)對(duì)于一個(gè)空間區(qū)域，處理器0位于最左區(qū)域等等。直到處理器n-1為最右區(qū)域，處理器0被指定為主線程。(2). 與并行命令parallel有關(guān)的參數(shù)proc, num 和master, 用來(lái)設(shè)定數(shù)據(jù)交換和處理器個(gè)數(shù)，并行處理通過(guò)參數(shù)start啟動(dòng)。通常將par命令放置在pfc3d.ini中。(3). 參數(shù)num只能在主線程中指定，用于定義并行處理處理器個(gè)數(shù)。參數(shù)proc用一個(gè)不同數(shù)字( procnum)確定每個(gè)線程，其中0為主線程， 每個(gè)從線程必須在通過(guò)參數(shù)master提供一個(gè)網(wǎng)絡(luò)名稱作為主線程。而對(duì)于主線程為可選項(xiàng)。這個(gè)名稱在一開(kāi)始用指定名稱為處理器尋找網(wǎng)絡(luò)。參數(shù)start啟動(dòng)并行操作，而參數(shù)end為反向操作即結(jié)束并行操作。如果并行操作未啟動(dòng)，PFC3d就像在幾個(gè)獨(dú)立處理器獨(dú)立運(yùn)行。PAUSE 中止data文件執(zhí)行。可通過(guò)下列三種方式調(diào)用：如果為指定參數(shù)，則PFC3d在遇到pause命令時(shí)中止，此時(shí)程序進(jìn)入交互模式，可以在鍵盤(pán)輸入命令，當(dāng)鍵入continue, 程序繼續(xù)運(yùn)行如果指定可選關(guān)鍵字key，則PFC3d在遇到pause命令時(shí)中止，擊打任意除ESC鍵外任何鍵，程序繼續(xù)運(yùn)行。輸入Esc鍵，程序中止，且返回交互模式。若指定變量t, 則PFC3d在遇到pause命令時(shí)中止，且等待t秒后繼續(xù)運(yùn)行程序。Plot keyword switch 屏幕繪圖或?qū)⑵聊粌?nèi)容輸出至硬拷貝設(shè)備或文件。繪圖輸出語(yǔ)法建立在視口基礎(chǔ)上，一個(gè)特定視口包括視口設(shè)置參數(shù)（背景，大小等）， 繪圖項(xiàng)目（模型表面，矢量等）。視圖輸出用一個(gè)默認(rèn)視圖初始化視口，其視口編號(hào)為0，名稱為“Base”。可創(chuàng)建和保存多個(gè)視口，可通過(guò)定義活動(dòng)視口（plot current）切換已保存視口。其中viewid 可以為整數(shù)（id）,也可為字符串（名稱）Plot create/ current * plot print 每個(gè)視口可以存儲(chǔ)許多繪圖項(xiàng)目（這些是特定繪圖項(xiàng)目，例如模型，以及速度矢量），繪圖項(xiàng)目增加（add）減少（subtract）修改（modify） 重新編碼（move ）顯示（ show ） 顯示視口輸出項(xiàng)目（ print item ）繪圖控制可分為四類：(1).視口控制 定義視口，輸出條件(2).視口設(shè)置控制 設(shè)置背景，前景色，視口位置設(shè)置，標(biāo)題，名稱(3).視口輸出項(xiàng)目控制在一個(gè)視口內(nèi)建立（加，減，修改）繪圖項(xiàng)目(4) 交互控制 在繪圖模式下，允許特定鍵盤(pán)輸入使得用戶控制視口。四類列表如下1. 視口控制(1). plot close 關(guān)閉當(dāng)前或指定視口(2). Plot viewid1 viewid2 將編號(hào)為viewid1視口復(fù)制至編號(hào)為viewid2視口。Settings: 視口設(shè)置 items: 視口輸出項(xiàng)目 both: 以上兩項(xiàng)若編號(hào)為viewid2視口為新視口，則使其為活動(dòng)視口。(3). Plot create viewid創(chuàng)建視口viewid，并使之為活動(dòng)視口(4) plot current viewid使得視口視口viewid為活動(dòng)視口(5) plot destroy viewid刪除視口viewid, 當(dāng)前視口和Base視口不能刪除(6) plot export 將生成當(dāng)前視口或指定視口的設(shè)置所需PFC3D的命令顯示在屏幕上，如果指定，則輸入至文件filename, 這個(gè)文件可通過(guò)call命令調(diào)用。(7) plot hardcopy 視口viewid或當(dāng)前視口對(duì)于當(dāng)前硬拷貝設(shè)備( set plot )通用。如果設(shè)備為一個(gè)生成文件，這文件名為PFC3D.*, 其中*決定當(dāng)前硬拷貝設(shè)備名稱，另外，可以直接用關(guān)鍵詞file直接指定文件名。當(dāng)從圖形界面創(chuàng)建一個(gè)硬拷貝輸出時(shí)，輸出窗口必須為全屏，這個(gè)輸出界面與出現(xiàn)在硬拷貝內(nèi)圖像相似，但圖像比例可能與視圖窗口尺寸不一樣。(8). Plot print view item 關(guān)鍵詞print = print list即輸出所有視圖視口和名稱表格View 輸出視口具體設(shè)置，可通過(guò)增加viewid增加一個(gè)指定視口Information 輸出當(dāng)前顯示和硬拷貝設(shè)置Item 輸出視口中具體項(xiàng)目列表，如果指定viewid，則輸出具體項(xiàng)目設(shè)置(9)plot quit退出繪圖模式，返回命令模式(10) plot rename viewid1 viewid2視口更名，當(dāng)base視口不能更名(11)plot show 顯示指定視口2. 視口設(shè)置控制(1) plot reset所有視口設(shè)置恢復(fù)為默認(rèn)狀態(tài)。(2) plot set keyword1).Plot set angle a設(shè)置視口角度，控制視圖編寫(xiě)，a值越大，變形越大。故要調(diào)整合適視角2).plot set animate off/on 視口圖像內(nèi)容改變是否，通常隨命令執(zhí)行改變或計(jì)算過(guò)程默認(rèn)為20步改變（可通過(guò)set pinterval改變），off 視口圖像不改變，默認(rèn)為on3). Plot set background color設(shè)置視口背景4). Plot set caption keyword 設(shè)置標(biāo)題位置 left 視口左側(cè) 默認(rèn)值 off ： 關(guān)閉標(biāo)題顯示，允許整個(gè)屏幕繪圖 on : 打開(kāi)標(biāo)題顯示（默認(rèn)） right: 視口右側(cè) size n: 設(shè)置標(biāo)題占視口大小百分比n% n = 10,50 默認(rèn)為n= 354). Plot set center ( x, y, z ) 設(shè)置模型視口中心為（x, y, z ）Auto 模型幾何中心5). Plot set color off/onOff: 不使用顏色，灰度模式 on: 打開(kāi)彩色模式，用ctrl G 打開(kāi)或關(guān)閉彩色模式6). Plot set eyedistance d 通過(guò)透視計(jì)算指定觀察點(diǎn)到屏幕屏幕距離 Auto 自動(dòng)計(jì)算7) plot set foreground color設(shè)置屏幕前景色為color8)plot set magnification m設(shè)置視口放大比例為m。9). Plot set mode keywordMode toggles使得用戶用更靈活方式觀察模型，因?yàn)閙odel 模式更方便觀察模型從空間一個(gè)靜止點(diǎn)。模式可以在繪圖窗口以交互方式改變。注意mode toggle 本身不代參數(shù)合適