1、一、建立模型筆記1鍵盤上的鍵的功能:X,Y,Z分別控制所建立的模型圍繞X,Y,Z三個軸做逆時針旋轉(zhuǎn).如果打開大小寫鎖定鍵,分邊按XYZ時建立的模型會做順時針旋轉(zhuǎn).+號的功能是加大模型旋轉(zhuǎn)的角度，-號的功能是減小模型旋轉(zhuǎn)的角度。2鍵盤上的上下左右四個鍵可以實現(xiàn)圖形的向上向下,向左和向右的移動.3當在一個區(qū)域建立多個模型的時候,有時候為了鏡像單個模塊兒,需要預(yù)先編組,編組的方法如下:gen zon tunint p0 0 0 0 p1 0 10 0 p2 5 5 0 p3 0 0 -5 p4 5 10 0 &p5 5 5 -5 p6 0 10 -5 p7 5 10 -5 dim sq2 sq2 1

2、 sq2 1 1 1 group(編組) z1(組名為:z1) &fill group tunnel4鏡像整個模型的方法:首先判斷即將鏡像的模型在原模型的哪個方向(即三個軸的方向),然后在軸所在的方向上去兩個點.其中一個是原點(origin(xyz).5鏡像部分模型的方法: gen zone reflect normal (x,y,z) origin (x,y,z) range group z1就可以了.6特別要注意的一點:交叉巷道的巷道充填和巷道的網(wǎng)格是兩個組,所以在鏡像的時候不要忘了給充填部分鏡像.7 group radcyl range group cylint not (編組的時候不包

3、括cylint)7 快捷鍵總結(jié)Ctrl+Z 局部縮放的功能。Ctrl+R的功能是恢復(fù)到初始狀態(tài)。8模型建立的方法:方法1:利用默認值生成網(wǎng)格,各個默認值:(1)p0的默認值是(0,0,0),網(wǎng)格的每邊的單元格數(shù)size默認值是10,網(wǎng)格每邊的長度默認值為size的值.方法2:利用4個點的坐標來生成矩形網(wǎng)格.p0p3 size方法3:利用edge來確定邊長生成矩形網(wǎng)格.方法4:利用參數(shù)ratio來確定單元體幾何邊畫邊生成矩形網(wǎng)格.方法5:利用參數(shù)add(相對坐標)來生成矩形網(wǎng)格.方法6:利用八個點的左邊來生成矩形網(wǎng)格.9、用戶自定義模型的方法我來貢獻一點自己的成果：; U j/ ?7 L+ a0

4、 uFLAC3D的二次開發(fā)環(huán)境提供了開放的用戶接口，在軟件安裝文件中包含了軟件自帶所有本構(gòu)模型的源代碼，且給出了Mohr-Coulomb模型和應(yīng)變軟化模型的編譯示例，因此可以方便地進行本構(gòu)模型的修改與開發(fā)。為了方便起見，下面的說明以建立UserModel模型為例。 c% b1 h. K6 Q: l(1) 在模型頭文件(usermodel.h)中進行新的本構(gòu)模型派生類的聲明，修改模型的ID(為避免與已有模型沖突，一般要求大于100)、名稱和版本，修改派生類的私有成員，主要包括模型的基本參數(shù)及程序執(zhí)行過程中主要的中間變量。( ? p6 u J5 Q3 y( a(2) 在程序C+文件(usermod

5、el.cpp)中修改模型結(jié)構(gòu)(UserModel:UserModel(bool bRegister): Constit -utive Model)的定義，這是一個空函數(shù)，主要功能是給(1)中定義的所有私有成員賦初值，一般均賦值為0.0。8 k6 x$ 8 N7 U! |6 Y(3) 修改const char *UserModel: roperties()函數(shù)，該函數(shù)包含了給定模型的參數(shù)名稱字符串，在FLAC3D的計算命令中需要用到這些字符串進行模型參數(shù)賦值。0 H; q3 _ ( X4 T- L, t% V. D(4) const char *UserModel:States()函數(shù)是單元在計

6、算過程中的狀態(tài)指示器，可以按照需要進行修改指示器的內(nèi)容。/ p! s$ c; h: + p- t(5) 按照派生類中定義的模型參數(shù)變量修改double UserModel:GetProperty()和void UserModel: SetProperty()函數(shù)，這兩個函數(shù)共同完成模型參數(shù)的賦值功能。 U; e G W Q# R4 q/ 9 G h(6) const char * UserModel:Initialize()函數(shù)在執(zhí)行CYCLE命令或大應(yīng)變模式下對于每個模型單元(zone)調(diào)用一次，主要執(zhí)行參數(shù)和狀態(tài)指示器的初始化，并對派生類聲明中定義的私有變量進行賦值。值得注意的是，Init

7、ialize()函數(shù)調(diào)用時沒有定義應(yīng)變分量，但可以調(diào)用應(yīng)力分量，但不能對應(yīng)力進行修改。9 n# e8 | c- B/ q, B B5 m(7) const char * UserModel:Run()是整個模型編制過程中最主要的函數(shù)，它對每一個字單元(sub-zone)在每次循環(huán)時均進行調(diào)用，由應(yīng)變增量計算得到應(yīng)力增量，從而獲得新的應(yīng)力。在計算過程中，要根據(jù)單元應(yīng)力情況對單元狀態(tài)指示器進行賦值。當進行塑性模型編制時，需對達到塑性的應(yīng)力狀態(tài)進行修正。2 / . U( A G. b3 6 . (8) 修改const char * UserModel:SaveRestore()中的變量，修改方法同(

8、2)和(5)，該函數(shù)的主要功能是對計算結(jié)果進行保存。; 8 ! Q! C. S# _5 (9) 程序的調(diào)試有兩種方法。在VC+的工程設(shè)置中將FLAC3D軟件中的EXE文件路徑加入到程序的調(diào)試范圍中，并將FLAC3D自帶的DLL文件加入到附加動態(tài)鏈接庫(Additional DLLs)中，然后在Initialize()或Run()函數(shù)中設(shè)置斷點，進行調(diào)試；在程序文件中加入return()語句，這樣可以將希望得到的變量值以錯誤提示的形式在FLAC3D窗口中得到。10、網(wǎng)格剖面的畫法方法：首先在命令行里輸入edit11、想知道有多少單元體和節(jié)點？方法：用命令 print info12、繪制塑形區(qū)域命

9、令是什么？方法： Plot block state shear ;獲得剪切屈服的單元，包括shear-n和shear-pPlot block state tension-p ;獲得過去拉伸屈服的單元Plot block state now ;獲得當前處于塑性狀態(tài)的單元，包括shear-n 和tension-nPlot block state past ;獲得過去處于塑性狀態(tài)的單元，包括shear-p 和tension-p13、step后面的數(shù)字是怎么確定的？答：14、開挖模型的命令Model null range group 組名Model null range group x x1,x2畫圈

10、的正負號是啥意思？答:正號表示10的正的幾次方.而負號表示負的幾次方.11、關(guān)于ratioratio代表單元體尺寸變化率，關(guān)于只需要記住一點：不同模型的兩個尺寸相等的面，采用相同的集合變化率，網(wǎng)格才會重合，ratio對于減少單元的數(shù)目有很好的效果。12、地震波的模擬命令流set dyn onset dyn multi onset largefree x y z range z -35.1 -34.9 ;去掉模型底部原有的靜力條件apply nquiet squiet dquiet ran z -35.1 -34.9 ;靜態(tài)邊界條件table 1 read crvel20.txtapply sz

11、z -1.e7 hist table 1 ran z -35.1 -34.9 ;加動力載荷apply syz -1.26984e7 hist table 1 ran z -35.1 -34.9apply ff ;設(shè)置阻尼set dyn damp rayleigh 0.05 7.14set dyn time=0/thread-1-1.html /thread-1-4.html 參考資料二、材料參數(shù)、邊界條件，初始條件1要特別注意fix x y z range z 與fix z range z 之間的差別，前者是x、y、z

12、三方向的位移全被限制，后者只限制了z方向的位移。fix x range (范圍)命令的功能是保持網(wǎng)格節(jié)點在x軸方向制定的參數(shù)(速度、壓力和溫度)不變，range為指定的范圍Fix x y z range x -0.1 0.1 表示x=0是一個平面，且是固定的Fix x range x -0.1 0.1 表示x=0是一個平面，只是x的方向不能發(fā)生運動。2、plot hist m vs n 其中m代表y軸，n代表x軸。3、sxx sxy syy這都是指那個方向:SXX是x方向的正應(yīng)力，SYY是y方向的正應(yīng)力, SXY是指法向指向x軸方向平面（也就是yz面）上的y方向剪應(yīng)力。4、history命令

13、Hist unbal 監(jiān)測不平衡力,并保留記錄 Hist gp zdisp 4,4,8 監(jiān)測網(wǎng)格坐標點(4,4,8)在z軸方向的記錄,并保留記錄. Plot his -2 -3 vs 1 其中的符號表示圖形的方向及位置的不同,1代表hist代碼中的第一條記錄,2代表hist代碼中的第二條記錄,3代表hist代碼中的第三條記錄.一次類推.補充：在記錄歷史記錄的時候可以自己定義記錄的編號：如 his id 12 gp zdisp 4,4,4 繪制圖形: 如 plot his 12 5、在fish函數(shù)中如：(1)p_gp=gp_head (意思是指向節(jié)點表) 那么后面肯定有一個： p_gp=gp_n

14、ext(p_gp) (意思是查找下一個節(jié)點的ID號) (2) pnt=zone_head ;zone_head指向單元體表 pnt=z_next(pnt) ;查找下一個單元體表 6、sel del link range id 23sel link id=23 23 target node tgt_num（代表目標節(jié)點的ID號） 1,7、想知道一個節(jié)點的x y z 三個分量則用它:;已知IDdefgetgpxyz gp1=find_gp(1) gp1x=gp_xpos(gp1) gp1y=gp_ypos(gp1) gp1z=gp_zpos(gp1) endgetgpxyzprintgp1xgp1

15、ygp1z 8、交集就是range group a group b指的是a和b的公共部分；并集就是range group a any group b any指的是a和b兩個部分；Not指的是出了什么什么之外。三、求解1、FLAC計算的過程中想終止pause，想要繼續(xù),continue2、收斂的準則：第一：相對收斂準則：一般而言，大多數(shù)問題可以采用FLAC3D默認的收斂標準（或稱相對收斂標準），即當體系最大不平衡力與典型內(nèi)力的比率R小于定值10-5（也可由用戶自定義該值，命令：SET mech ratio ）時，計算即行終止。這里，闡述兩個名詞的定義。所謂體系最大不平衡力，是指每一個計算循環(huán)(或

16、稱計算時步)中，外力通過網(wǎng)格節(jié)點傳遞分配到體系各節(jié)點時，所有節(jié)點的外力與內(nèi)力之差中的最大值；所謂典型內(nèi)力，則是指計算模型所有網(wǎng)格點力的平均值。圖 2-9為簡單示例中R默認為10-5時，計算終止時所花費的時步以及R值。由于R為無量綱的，所以它適用于不同的單位系統(tǒng)。有時，也可采用體系最大不平衡力小于某一臨界值（命令：SET mech force ）作為一個收斂標準，也稱絕對收斂標準。由于這一臨界值需用戶自行定義，并且沒有比較統(tǒng)一的取值范圍，若取較小值，對單元數(shù)眾多的復(fù)雜模型而言無疑是十分苛刻的，因而有較大的局限性。 第二：用戶自定義準則：FLAC3D還允許用戶根據(jù)實際分析的需要，自定義收斂標準，比

17、如以某關(guān)鍵點的位移或速度達到某一臨界值作為收斂標準等等。3、記錄最大不平衡力：his unbal 記錄最大不平衡力（命令流位置），plot his 1 繪制最大不平衡力曲線圖。注意：體系最大不平衡力隨著計算的進行，逐漸趨近于0，表明體系最終達到了力平衡狀態(tài)。不過，力平衡狀態(tài)僅表示所有網(wǎng)格節(jié)點的合力為0，并非表明體系處于真實的物理平衡狀態(tài)，因為在力平衡狀態(tài)下，體系也有可能正在發(fā)生穩(wěn)定的塑性流動。這時，就需要借助其它方式如觀察網(wǎng)格節(jié)點速度、塑性區(qū)標識等來進一步評估模型所處的狀態(tài)。（平衡狀態(tài)就是說靜止狀態(tài)和勻速直線運動狀態(tài)。）4、觀察網(wǎng)格節(jié)點速度的方法：（1）跟蹤記錄網(wǎng)格中一些關(guān)鍵點的速度（命令：H

18、IST gp vel ）并繪制其歷時曲線。 注意： 如果速度歷時曲線在最后階段顯示為水平線，表明體系達到了穩(wěn)定狀態(tài)。 如果速度歷時曲線最終收斂并趨近于零，如果速度歷時曲線最終收斂并趨近于零。 如果速度歷時曲線最終收斂并趨近于零，表明與記錄相應(yīng)的網(wǎng)格點進入塑性流動狀態(tài)； 如果一個或多個點的速度的歷時曲線都出現(xiàn)明顯的上下波動，表明系統(tǒng)此時可能處于瞬時調(diào)整狀態(tài)。（2）采用PLOT vel命令繪制完整的速度矢量場圖。注意：通常，體系在達到真實平衡狀態(tài)時，因網(wǎng)格節(jié)點力的變化已低于計算機的精度（約小數(shù)點后六位），速度場表現(xiàn)為大小和方向不確定的低振幅隨機場。如果速度場矢量是一致的，且量級很大，則體系有可能正

19、發(fā)生塑性流動或正進行彈性調(diào)整（如發(fā)生彈性阻尼振動）。在這種情況下，要確認塑性流動是否正在發(fā)生，則應(yīng)觀察塑性區(qū)標識來研判體系所處的狀態(tài)。5、塑性區(qū)標識命令：plot block state 某一區(qū)域的應(yīng)力正位于屈服面上，或者說正處于破壞狀態(tài)時，以shear-n 或tension-n 標識； 某一區(qū)域在計算過程中曾進入過屈服狀態(tài)，但現(xiàn)在已經(jīng)退出，以shear-p 或tension-p 標識。 判斷破壞機制是否起作用的步驟: 觀察連接兩個面的活性塑性區(qū)域（由shear-n 或tension-n 標識）是否存在交線（或有交集），若存在，表示破壞機制正在起作用； 觀察速度矢量圖，若速度矢量圖出現(xiàn)了與破壞

20、機制相一致的運動，則確認第（1） 步的判斷； 若活性塑性區(qū)域不存在交線或交集，則應(yīng)再執(zhí)行一些時步后比較前后兩種圖形：若塑性區(qū)域減少了，體系可能正向平衡發(fā)展；若增加了，體系最終有可能發(fā)生破壞。如果確認是連續(xù)塑性流動的情形，則有一個問題需注意，人工邊界的選取有可能影響塑性區(qū)的范圍。所謂“人工邊界”是指用來界定網(wǎng)格模型大?。ɑ蚍治鲇颍┑倪吔?。如果塑性流動沿著人工邊界發(fā)生，那么求解是不真實的，因為破壞機制受到了人工邊界的影響。因此，在實際分析時，應(yīng)做些簡單的人工邊界設(shè)置比較分析，將其設(shè)置在對分析所關(guān)注區(qū)域的應(yīng)力和位移影響較小的區(qū)域。6、計算塑性區(qū)域的面積def get_plastshearnow =

21、1 tensionnow = 2 shearpast = 4 tensionpast = 8 v_shear_now = 0v_tension_now = 0v_shear_past = 0v_tension_past = 0 ;flac默認的值;xx=z_state(find_zone(11779),0) ;通過本步可驗證 p_z = zone_headloop while p_z # nullif and(z_state(p_z,0),shearnow) = shearnow thenv_shear_now = v_shear_now + z_volume(p_z)endifif and(

22、z_state(p_z,0),tensionnow) = tensionnow thenv_tension_now = v_tension_now + z_volume(p_z)endifif and(z_state(p_z,0),shearpast) = shearpast thenv_shear_past = v_shear_past + z_volume(p_z)endifif and(z_state(p_z,0),tensionpast) = tensionpast thenv_tension_past = v_tension_past + z_volume(p_z)endifp_z

23、= z_next(p_z)endloopii = out(shear_now : + string(v_shear_now)ii = out(tension_now : + string(v_tension_now)ii = out(shear_past : + string(v_shear_past)ii = out(tension_past : + string(v_tension_past)endget_plast運行fish后，執(zhí)行pri fish就可以顯示shearnow（當前剪切破壞）、tensionnow（當前拉帳破壞）、shearpast（過去剪切破壞）、tensionpast（過去拉帳破壞） 的體積值。/thread-1-2.html有解釋四、基本命令1、怎么統(tǒng)計網(wǎng)格單元體的個數(shù)？答：在FLAC3D中輸入 prin information 2、怎么畫出網(wǎng)格的剖面？答：set pla dip 90 dd 90 ori 2 0 03、繪制網(wǎng)格線框命令？答：Plot grid4、應(yīng)力集中系數(shù)的求法。答：應(yīng)力集中系數(shù)是在巷道開挖過程中涉及到的名詞，在巷道沒有開挖之前某個水平有一個原始應(yīng)力值，開挖之后應(yīng)力的重新分布，之前的那個水平上的應(yīng)力值會變成成一個新的值。這個值呢可以通