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

2、為注釋, FLAC 在執(zhí)行時不理會。 下面的 FLAC 命令按字母排列。Apply 關鍵詞 = 數(shù)值 可以有下面的關鍵詞: Pressure 壓力 XForce X-方向的力 YForce Y-方向的力ATtach 該命令可以將一條線上的結點和另一條線上的結點互相接合在一起, 用以形成復雜的網格形狀。Call 文件名 寫成的命令文件可以用 Call 命令來調用, 命令文件的最后一行必須是RETURN, 以返回到交互方式。命令文件中不能有 CALL 命令本身。Config 關鍵詞 FLAC 用以解平面應變問題, 但經過配置命令也可以用于解平面應力問題或軸對稱問題。 需要時應在形成網格之前發(fā)。關鍵

3、詞有: 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 命令相反, 用于放松對結點的約束。GEn x1,

4、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), 逆時針的角為theta 度 CIRCLE 園的中心為 (

5、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 版本, 如用摩爾-庫侖的本構模型約可分兩千個單元。Help 幫助命令, 可在屏幕上顯

6、示命令表。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 的 xx-應力 SYY 單元 i,j 的 yy-應力 SXY 單元 i

7、,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 的 xy-應力 YYB 單元 i,j 三角形 b 的 yy-應力 XXC 單元 i,j 三角形 c 的

8、 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 WRITE 命令可將結果順序寫在 FLAC.HIS 上, 但是首次執(zhí)行此命令會將以前盤上所存的同名文件沖掉; Reset

9、所有的歷史都清除掉。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 在一定范圍內的參數(shù)的變差值, xv 和 yv 為 x-方向和 y-方向的變差值。 范圍可以有形式 I=i1,j1 J=j1,j2

10、, 其中 i 和 j 的先后次序無關。INTerface n 關鍵詞 n 關鍵詞 = 值 INTERFACE 為交界面命令。網格的一部分可以和網格的另一部分通過交界面而相互作用。界面的性質由其剛度, 粘結力和摩擦力來表征。用于本命令的關鍵詞有: Aside i1,j1 Bdide i1,j1 Cohesion 粘結力值 Friction 摩擦力值 Glued KN 法向剛度值 KS 切向剛度值 TBond 抗拉強度 Unglued 膠結命令 GLUE 將使交界面的上下兩部分膠結在一起, 使之不能滑動和分離, 命令 UNGLUE 則為命令 GLUE 之逆。MArk 在給定范圍內的結點將被標記。計

11、算的進行與結點的標記與否無關, 但標記的結點可以限定一個區(qū)域作為命令 INITIAL, PROP 和 MODEL 的作用域。命令 GEN 將自動對結點進行標記。的形式可以是 I=i1, i2, J=j1,j2, i 和 j 的順序無關。Model 關鍵詞 該命令對于給定的區(qū)域或范圍賦予有關的本構模型。關鍵詞有: ANisotropic 橫觀各向同性彈性模型 Elastic 彈性, 各同性本構模型 Mohr-Coulomb 摩爾-庫倫塑性模型 Mull 零模型, 用于開挖掉的單元 SS 應變軟化 Ubiquitous 徹體節(jié)理模型NEW 該命令可以不用退出 FLAC 重新開始一個新問題。Plot

12、 關鍵詞 開關 . 用缺省值。可以在一行上連寫幾個關鍵詞以在一個圖上畫出幾個變量。在繪圖前先要用 MODEL 命令賦予網格以一定的本構模型。關鍵詞及其意義為: BEam 畫出結構單元的幾何形狀 Boundary 畫出網格的內外邊界 Cable 畫出錨束的幾何形狀 Disp 位移矢量 E_p 塑性應變的輪廓線 (只限于應變軟化材料) Grid 畫出 畫出紀錄在第 nhis 個歷史上的變量值 PP 孔隙壓力的等值線 RF 以矢量形式畫出的固定結點的反力 STAte 畫出單元中心當前的屈服狀態(tài) STress 主應力矢量 SXX xx-應力的等值線 (全應力) SYY yy-應力的等值線 (全應力)

13、SXY xy-應力的等值線 Velocity 以箭頭畫出的速度矢量 WAter 水位線 XDisp x-位移的等值線 YDisp y-位移的等值線 XVel x-速度 YVel y-速度 所謂開關其本身也是關鍵詞, 用于設置圖形的的某些特征, 有: color 如紅為 Red, 綠為 GREen, 黃為 Yellow 等。 Interval = c 將等值線的間距置為 c. 在 PLOT 命令前發(fā)掃描線命令 SCLIN 可以畫出等值線的值。 Max = v 在畫矢量時, 置箭頭的最大長度為 v. Noh 不畫標題 Zero 不畫零值線 如果要放大或縮小圖形, 可以在 PLOT 命令前用 WIN

14、DOW 命令來設置窗口的大小。Print 關鍵詞 . 和繪圖命令 PLOT 一樣, 結點變量的值只有在給定的材料模型和性質后才可以打印出來。關鍵詞有: Apply 打印出所施加力或壓力的大小和范圍 Fix 打印出固定 x 或 y 的結點 Limits 對 SOLVE 命令的限制值 MArk 打印出標記結點 MEm 占用的內存 STruct 打印出結構單元上有關結點的力, 力矩和位移 Interface 打印出交界面的數(shù)據, 包括結點力和單位法向矢量 主要的網格關鍵詞 Bulk 體積模量 COhesion 粘結力 Den 質量密度 Dilation 剪脹角 E_p 全塑性應變 (應變軟化模型)

15、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-模量 XVel X-速度 YVel Y-速度 ASXX XX-應力 (三角形 A) BSXX XX-應力 (三角形 B) CSXX XX-應力 (三角形 C) DSX

16、X 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 已超過單軸抗張力 4 屈服并超過單軸抗張力 5 已經超過抗張力 6 徹體節(jié)理正在屈服 7 徹體節(jié)理過去屈服過, 現(xiàn)為彈性狀態(tài) Tables 打印所存貯的表格 XReaction X-反力

17、 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 孰者在前無妨, 但 必需在輸入行的最后。 各種模型所需輸入的性質為: 彈性 (1) 剪切模量 (2) 體積模量 (3) 密度 摩爾-庫倫 (1) 剪切模量 (2) 體積模量 (3) 密度 (4) 摩擦角

18、 (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) 整體材料的剪脹角 (任選) 應變-硬化/軟化 (1) 剪切模量 (2) 體積模量 (3) 密度 (4) 初始摩擦角 (任選) (5) 初始粘結力 (任選) (6) 摩擦角和塑性應變的關系表的表號 (7) 粘結力和塑性應變的關系表的表號 (8) 剪脹角和塑

19、性應變的關系表的表號 如果表號給 0, 則取關鍵詞 COHESION, DILATION 或 FRICTION 所給出的值。Quit 該命令同 STOP, FLAC 終止運行。REstore 文件名 把以前用 SAVE 命令所存的文件讀入內存恢復現(xiàn)場。RETurn 該命令應是輸入數(shù)據文件中的最后一個命令。SAve 文件名 將內存中問題的現(xiàn)場存入文件中, 如果已經有同名的文件, 則該文件將被復蓋。SClin n x1,y1 x2,y2 該命令在屏幕上作一個掃描線, 以切割等值線, 交點以 A 到 Z 的字母來表示。在一張圖上最多可以有 5 條掃描線。掃描線的參數(shù)為: n = 掃描線的編號 (必需

20、為 1,2,3,4 或 5) x1,y1 = 掃描線的起點坐標 x2,y2 = 掃描線的終點坐標SEt 該命令設置整體條件, 其關鍵詞有: Aspect a 調整屏幕的高寬比, 以使輸出不變形 Col n n 為打印輸出的最大列數(shù)。用戶必需在運行 FLAC 前設置好打印機的列數(shù)。 EGA 設置 640x350 的 EGA 圖形模式 Force f 設置 SOLVE 命令時的不平衡力的限值。 Gravity g 重力加速度 g 及其與負 y-軸的夾角 Large 大變形 (坐標更新) Log ON 在當前盤上打開名為 FLAC.LOG 的文件, 以紀錄 FLAC 的運行進程。如果 FLAC.LO

21、G 文件已經存在, 則將被復蓋。 OFF 關閉 FLAC.LOG 文件。如果在稍后再 SET LOG ON, 則屏幕顯示將繼續(xù)紀錄在 FLAC.LOG 文件上。 Small 小變形 (坐標不更新) STep n SOLVE 命令進行計算循環(huán)的限制數(shù)目。 Time t SOLVE 命令進行計算的限制時間的分鐘數(shù)。SOLve 用本命令進行計算循環(huán), 關鍵詞為: Step 時步數(shù)目 Time 運行時間 (分鐘數(shù)) Force 不平衡力 缺省值為: S = 500 時步 T = 5 分鐘 F = 100 最長的時限為 1440 分鐘 (24 小時), 如果要求運算超過一天, 則可以用幾個 SOLVE

22、的命令。FLAC 在運行中如要中止, 可以按退出鍵 Esc退出循環(huán), 將控制轉給用戶, 進行交互操作STEp n 執(zhí)行 n 個時步的循環(huán)運算。STop FLAC 終止運行。STRuct 關鍵詞 . . . STRUCT 命令用于定義結構單元的幾何形狀, 性質等條件。STRUCT 命令要求說明結構單元的類型和形狀以及它與 FLAC 網格的連接情況。定義結構單元的關鍵詞有: (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) 結構單元的性質類型關鍵