FLAC程序基本命令

1、flac程序使用手冊flac 輸入命令 flac 的輸入和一般的數(shù)值模擬的程序不一樣, 它可以用交互的方式從鍵盤輸入各個命令, 也可以寫成命令文件, 類似于批處理, 由文件來驅動。flac 命令大小寫一樣。所有的命令可以附帶若干個關鍵詞和有關的數(shù)值。在下面的命令解釋中 , 只有大寫的字母起作用, 小寫的字母寫不寫、寫多少個都沒有崐關系。i,j,m 和 n 開始的變量要求整型數(shù), 否則要求實型數(shù)。?實型數(shù)的小數(shù)點可崐以忽略, 但是整型數(shù)不能帶小數(shù)點。數(shù)值間可以用空格隔開, 空格的數(shù)目不限,? 也可以用下面的分隔符隔開: ( ), / = 表示可選的參數(shù), 輸入時括號不用輸入; . 表示可以有任意

2、個參數(shù)。由 * 號開始到行末為注釋, flac 在執(zhí)行時不理會。下面的 flac 命令按字母排列。apply 關鍵詞 = 數(shù)值 可以有下面的關鍵詞: pressure 壓力 xforce x-方向的力 yforce y-方向的力attach 該命令可以將一條線上的結點和另一條線上的結點互相接合在一起, 用以形成復雜的網格形狀。call 文件名寫成的命令文件可以用 call 命令來調用 , 命令文件的最后一行必須是return, 以返回到交互方式。命令文件中不能有 call 命令本身。config 關鍵詞 flac 用以解平面應變問題, ? 但經過配置命令也可以用于解平面應力問題或軸對稱問題。需

3、要時應在形成網格之前發(fā)。關鍵詞有: p_str 平面應力問題 ax 軸對稱問題cyc n 該命令同 step, 為執(zhí)行 n 個時步的循環(huán)運算。fix x y x y 用此命令可以使 內結點的 x- 或 y- 方向的速度保持不變。 的格式可以是 i = i1,i2, j = j1,j2; i 和 j ? 何者先輸入沒有關系。如果要求位移固定, 則必須將速度初始化為零 ( 開始時速度的缺省值為零) 。fix 和 ini xv, yv ?可以聯(lián)合使用以提供一個剛性移動的邊界條件。如果有 mark 的關鍵詞 , 則只有在此范圍內被標記的結點將被固定。free x y x y 該命令與 fix 命令相反

4、 , 用于放松對結點的約束。gen x1,y1 x2,y2 x3,y3 x4,y4 用于在全域或局域中產生網格。(x1,y1)的點放在左下側, (x2,y2)、 (x3,y3) 和 (x4,y4) 為順時針排列。如果點和前一個 gen 命令中的相同則用same 來代替。如果網格不是均勻排列, 則用關鍵詞 ratio,? 在 i 和 j 方向的比例系數(shù)分別為 ri 和 rj. gen circle xc,yc rad arc xc,yc xb,yb theta line x1,y1 x2,y2 該命令產生園、弧或直線。 arc 弧的圓心為 (xc,yc), 起始點為 (xb,yb), 逆時針的角

5、為theta 度 circle 園的中心為 (xc,yc), 半徑為 rad line 直線的兩端為 (x1,y1) 和 (x2,y2). 用 arc 等關鍵詞所形成的邊界系被標記 ( 參見命令 mark 和 unmark) 過 , 標記結點所圍住的區(qū)域在以后用命令 inital, model 和 property 時可以用關鍵詞 region 來說明。gen adjust 用此命令來微調網格, 使之離散平滑。該命令可以連續(xù)使用以增進平滑效果的。grid icol jrow 此命令用以產生數(shù)目為 icol 列和 jrow 行的網格。對于 640k 內存的flac 版本, 如用摩爾 -庫侖的本構

6、模型約可分兩千個單元。help 幫助命令 , 可在屏幕上顯示命令表。his 每隔 nstep 時步 , 記錄一下關鍵詞所示項目在 i1 和 ji ?結點或單元的數(shù)值。nstep 的缺省值為 10. 用戶應記住所要求 his ( 歷史 )?的順序 , 因以后在打印或繪圖時要用。歷史值在 flac 停止運行時刪除, 因此如果要保留的話, 要用到命令 his write n ( 見下面 )。可以用的關鍵詞有: ang 單元 i,j 內的最小主應力和 x- 軸所形成的夾角 pp 單元 i,j 內的孔隙壓力 sig1 單元 i,j 的最大主應力 sig2 單元 i,j 的最小主應力 sxx 單元 i,j

7、 的 xx- 應力 syy 單元 i,j 的 yy- 應力 sxy 單元 i,j 的 xy- 應力 x 結點 i,j 的 x- 坐標 y 結點 i,j 的 y- 坐標 xdis 結點 i,j 的 x- 位移 ydis 結點 i,j 的 y- 位移 unbal 最大不平衡力 xvel 結點 i,j 的 x- 速度 yvel 結點 i,j 的 y- 速度 xxa 單元 i,j 三角形 a 的 xx- 應力 xya 單元 i,j 三角形 a 的 xy- 應力 yya 單元 i,j 三角形 a 的 yy- 應力 xxb 單元 i,j 三角形 b 的 xx- 應力 xyb 單元 i,j 三角形 b 的

8、xy- 應力 yyb 單元 i,j 三角形 b 的 yy- 應力 xxc 單元 i,j 三角形 c 的 xx- 應力 xyc 單元 i,j 三角形 c 的 xy- 應力 yyc 單元 i,j 三角形 c 的 yy- 應力 xxd 單元 i,j 三角形 d 的 xx- 應力 xyd 單元 i,j 三角形 d 的 xy- 應力 yyd 單元 i,j 三角形 d 的 yy- 應力his 關鍵詞關鍵詞可以是 : dump nhis 將第 nhis 的歷史寫屏 ; write nhis 將第 nhis 的歷史寫在文件 flac.his 上, ?該文件可在 flac 結束后打印出來。連續(xù)執(zhí)行 his wr

9、ite 命令可將結果順序寫在 flac.his 上, ? 但是首次執(zhí)行此命令會將以前盤上所存的同名文件沖掉; reset 所有的歷史都清除掉。initial 關鍵詞 = 值 某些結點值可以給初值: pp 空隙壓力 x x-坐標 y y-坐標 sxx xx-應力 syy yy-應力 sxy xy-應力 xdisp x-位移 ydisp y-位移 xvel x-速度 yvel y-速度還有幾個可選用的關鍵詞以協(xié)助實現(xiàn) initial 命令 : mark 只有標記的結點值將被初始化 region i j 為標記結點所圍住的區(qū)域將被初始化。i,j 為標記區(qū)域中的任一個單元 var xv, yv 在一定

10、范圍內的參數(shù)的變差值, xv 和 yv 為 x- 方向和 y-? 方向的變差值。范圍可以有形式 i=i1,j1 j=j1,j2, 其中 i 和 j 的先后次序無關。interface n 關鍵詞 n 關鍵詞 = 值 interface 為交界面命令。網格的一部分可以和網格的另一部分通過交界面而相互作用。界面的性質由其剛度, 粘結力和摩擦力來表征。用于本命令的關鍵詞有: aside i1,j1 bdide i1,j1 cohesion 粘結力值 friction 摩擦力值 glued kn 法向剛度值 ks 切向剛度值 tbond 抗拉強度 unglued 膠結命令 glue 將使交界面的上下兩

11、部分膠結在一起, 使之不能滑動和分離, 命令 unglue 則為命令 glue 之逆。mark 在給定范圍內的結點將被標記。計算的進行與結點的標記與否無關, 但標記的結點可以限定一個區(qū)域作為命令 initial, prop 和 model 的作用域。命令 gen 將自動對結點進行標記。 的形式可以是 i=?i1, i2, j=j1,j2, i 和 j 的順序無關。model 關鍵詞 該命令對于給定的區(qū)域或范圍賦予有關的本構模型。關鍵詞有: anisotropic 橫觀各向同性彈性模型 elastic 彈性 , 各同性本構模型 mohr-coulomb 摩爾 - 庫倫塑性模型 mull 零模型

12、, 用于開挖掉的單元 ss 應變軟化 ubiquitous 徹體節(jié)理模型new 該命令可以不用退出 flac 重新開始一個新問題。plot 關鍵詞 開關 . 用缺省值。可以在一行上連寫幾個關鍵詞以在一個圖上畫出幾個變量。在繪圖前先要用 model 命令賦予網格以一定的本構模型。?關鍵詞及其意義為: beam 畫出結構單元的幾何形狀 boundary 畫出網格的內外邊界 cable 畫出錨束的幾何形狀 disp 位移矢量 e_p 塑性應變的輪廓線 ( 只限于應變軟化材料) grid 畫出畫出紀錄在第 nhis 個歷史上的變量值 pp 孔隙壓力的等值線 rf 以矢量形式畫出的固定結點的反力 sta

13、te 畫出單元中心當前的屈服狀態(tài) stress 主應力矢量 sxx xx-應力的等值線 ( 全應力 ) syy yy-應力的等值線 ( 全應力 ) sxy xy-應力的等值線 velocity 以箭頭畫出的速度矢量 water 水位線 xdisp x-位移的等值線 ydisp y-位移的等值線 xvel x-速度 yvel y-速度所謂開關其本身也是關鍵詞, 用于設置圖形的的某些特征, 有: color 如紅為 red, 綠為 green, 黃為 yellow 等。 interval = c 將等值線的間距置為 c.? ?在 ?plot ? 命令前發(fā)掃描線命令 sclin 可以畫出等值線的值。

14、 max = v 在畫矢量時 , 置箭頭的最大長度為 v. noh 不畫標題 zero 不畫零值線如果要放大或縮小圖形, 可以在 plot 命令前用 window 命令來設置窗口的大小。print 關鍵詞 . 和繪圖命令 plot 一樣 , 結點變量的值只有在給定的材料模型和性質后才可以打印出來。關鍵詞有: apply 打印出所施加力或壓力的大小和范圍 fix 打印出固定 x 或 y 的結點 limits 對 solve 命令的限制值 mark 打印出標記結點 mem 占用的內存 struct 打印出結構單元上有關結點的力, 力矩和位移 interface 打印出交界面的數(shù)據, 包括結點力和單

15、位法向矢量主要的網格關鍵詞 bulk 體積模量 cohesion 粘結力 den 質量密度 dilation 剪脹角 e_p 全塑性應變 ( 應變軟化模型) friction 摩擦系數(shù) shear 剪切模量 x x-坐標 y y-坐標 pp 孔隙壓力 jfric 節(jié)理摩擦力 jcoh 節(jié)理粘結力 jangle 節(jié)理角度 nuyx y-x 泊松比 nuzx z-x 泊松比 sig1 最大主應力 sig2 最小主應力 sxx xx-應力 syy yy-應力 sxy xy-應力 theta 最小主應力與 x- 軸的夾角 xdis x-位移 ydis y-位移 xmod x-模量 ymod y-模量

16、xvel x-速度 yvel y-速度 asxx xx-應力 ( 三角形 a) bsxx xx-應力 ( 三角形 b) csxx xx-應力 ( 三角形 c) dsxx xx-應力 ( 三角形 d) asyy yy-應力 ( 三角形 a) bsyy yy-應力 ( 三角形 b) csyy yy-應力 ( 三角形 c) dsyy yy-應力 ( 三角形 d) asxy xy-應力 ( 三角形 a) bsxy xy-應力 ( 三角形 b) csxy xy-應力 ( 三角形 c) dsxy xy-應力 ( 三角形 d) state 塑性狀態(tài) 0 彈性 1 正在屈服中 2 曾經屈服 , 現(xiàn)為彈性 3

17、已超過單軸抗張力 4 屈服并超過單軸抗張力 5 已經超過抗張力 6 徹體節(jié)理正在屈服 7 徹體節(jié)理過去屈服過, 現(xiàn)為彈性狀態(tài) tables 打印所存貯的表格 xreaction x-反力 yreaction y-反力prop 關鍵詞 = 值 本命令為 model 命令賦予材料的性質。下面的關鍵詞對 property 命令加以補充: region i,j 所有為標記結點所連續(xù)包圍的區(qū)域將賦予該性質。i,j 單元可為標記區(qū)域中的任一單元 var xv yv 性質可以有變差, xv 和 yv 分別為 x- 方向和 y-? 方向的變差值 可為 i=i1,i2, j=j1,j2, i 和 j 孰者在前無

18、妨, 但 必需在輸入行的最后。各種模型所需輸入的性質為: 彈性 (1) 剪切模量 (2) 體積模量 (3) 密度摩爾 - 庫倫 (1) 剪切模量 (2) 體積模量 (3) 密度 (4) 摩擦角 (5) 粘結力 (6) 剪脹角 ( 任選 ) 橫觀各向同性 (1) 剪切模量 (2) x-模量 (3) y-模量 (4) 密度 (5) nuyx (6) nuyz 徹體節(jié)理 (1) 剪切模量 (2) 體積模量 (3) 密度 (4) 粘結力 ( 整體材料 ) (5) 摩擦力 ( 整體材料 ) (6) 節(jié)理粘結力 (7) 節(jié)理摩擦力 (8) 節(jié)理角 (9) 整體材料的剪脹角 ( 任選 ) 應變 - 硬化 /

19、 軟化 (1) 剪切模量 (2) 體積模量 (3) 密度 (4) 初始摩擦角 ( 任選 ) (5) 初始粘結力 ( 任選 ) (6) 摩擦角和塑性應變的關系表的表號 (7) 粘結力和塑性應變的關系表的表號 (8) 剪脹角和塑性應變的關系表的表號如果表號給 0, 則取關鍵詞 cohesion, dilation 或 friction 所給出的值。quit 該命令同 stop, flac 終止運行。restore 文件名把以前用 save 命令所存的文件讀入內存恢復現(xiàn)場。return 該命令應是輸入數(shù)據文件中的最后一個命令。save 文件名將內存中問題的現(xiàn)場存入文件中, 如果已經有同名的文件, 則

20、該文件將被復蓋。sclin n x1,y1 x2,y2 該命令在屏幕上作一個掃描線, 以切割等值線, 交點以 a 到 z 的字母來表示。在一張圖上最多可以有 5 條掃描線。掃描線的參數(shù)為: n = 掃描線的編號 ( 必需為 1,2,3,4 或 5) x1,y1 = 掃描線的起點坐標 x2,y2 = 掃描線的終點坐標set 該命令設置整體條件, 其關鍵詞有 : aspect a 調整屏幕的高寬比, 以使輸出不變形 col n n 為打印輸出的最大列數(shù)。用戶必需在運行 flac 前設置好打印機的列數(shù)。 ega 設置 640 x350 的 ega 圖形模式 force f 設置 solve 命令時的

21、不平衡力的限值。 gravity g 重力加速度 g 及其與負 y- 軸的夾角 large 大變形 ( 坐標更新 ) log on 在當前盤上打開名為 flac.log 的文件 , 以紀錄 flac 的運行進程。如果 flac.log 文件已經存在, 則將被復蓋。 off 關閉 flac.log 文件。如果在稍后再 set log on, ? 則屏幕顯示將繼續(xù)紀錄在flac.log 文件上。 small 小變形 ( 坐標不更新 ) step n solve 命令進行計算循環(huán)的限制數(shù)目。 time t solve 命令進行計算的限制時間的分鐘數(shù)。solve 用本命令進行計算循環(huán), 關鍵詞為 :

22、step 時步數(shù)目 time 運行時間 ( 分鐘數(shù) ) force 不平衡力缺省值為 : s = 500 時步 t = 5 分鐘 f = 100 最長的時限為 1440 分鐘 (24 小時 ), 如果要求運算超過一天, 則可以用幾個solve 的命令。 flac 在運行中如要中止, 可以按退出鍵 esc 退出循環(huán) , 將控制轉給用戶, 進行交互操作step n 執(zhí)行 n 個時步的循環(huán)運算。stop flac 終止運行。struct 關鍵詞 . . . struct 命令用于定義結構單元的幾何形狀, 性質等條件。struct 命令要求說明結構單元的類型和形狀以及它與 flac 網格的連接情況。定

23、義結構單元的關鍵詞有: (1) 結構類型關鍵詞梁一個梁單元系在關鍵詞 beam 后用梁的幾何形狀和性質的關鍵詞及其值來表征。錨束一個錨束單元系在關鍵詞 cable ?后以其幾何形狀和性質的關鍵詞及其值來表征。 (2) 結構單元的幾何形狀及其結點聯(lián)結 begin 關鍵詞關鍵詞 begin 表征梁或錨束的起始點 end 關鍵詞 end 表征梁或錨束的終了點結構單元的端點可以在 begin 和 end 關鍵詞后直接用下面三個之一的關鍵詞來表征 : grid i,j 結構單元的起點或終點為結點 i,j node n flac ?將自動按照構件單元輸入的順序為結構單元的結點編號, 此處n 為所編的第 n 號結點。 x,y 直接給端點以 x 和 y 坐標。 (3) 結構單元的性質類型關鍵詞 prop n 結構單元