FLAC程序基本命令

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

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

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

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

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

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

7、 的 xx- 應(yīng)力 syy 單元 i,j 的 yy- 應(yīng)力 sxy 單元 i,j 的 xy- 應(yīng)力 x 結(jié)點(diǎn) i,j 的 x- 坐標(biāo) y 結(jié)點(diǎn) i,j 的 y- 坐標(biāo) xdis 結(jié)點(diǎn) i,j 的 x- 位移 ydis 結(jié)點(diǎn) i,j 的 y- 位移 unbal 最大不平衡力 xvel 結(jié)點(diǎn) i,j 的 x- 速度 yvel 結(jié)點(diǎn) i,j 的 y- 速度 xxa 單元 i,j 三角形 a 的 xx- 應(yīng)力 xya 單元 i,j 三角形 a 的 xy- 應(yīng)力 yya 單元 i,j 三角形 a 的 yy- 應(yīng)力 xxb 單元 i,j 三角形 b 的 xx- 應(yīng)力 xyb 單元 i,j 三角形 b 的

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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