Flac3D常用命令流

1、文檔供參考，可復(fù)制、編制，期待您的好評與關(guān)注！ FLAC3D 常用命令流  模型鏡像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一點）;沿X軸鏡像,通過對稱平面法線向量確定對稱面gen zone reflect normal 0 0 -1 ;沿z軸鏡像;繪圖控制pl <add> contour szz outline on;在模型中顯示位移-應(yīng)變曲線hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone

2、syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y軸,n代表x軸(不管m,n的正負(fù)); "-"表示對其值作"mirror"對模型進(jìn)行壓縮實驗的方法;即在模型兩側(cè)施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1;修改模型的坐標(biāo)值ini x add -100 y add -100 z add -100;顯示云圖的同時也顯示模型網(wǎng)格輪廓plot add cont disp outline on <gradient

3、/average>gradient更精確;輸入角度、弧度方法 pi=，90°為90.0*degraddef set_valsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint ptA ptB ptC;施加結(jié)構(gòu)單元方法sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &end2=(0.0,25.0,0.0) radius=24.5 notplot add sel

4、 geom black black cid on scale=0.03sel node init zpos add -25.0;如何顯示某一平面plot create name_planeplot set plane  origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom blackplot add axes redplot current name_planeplot show;如何施加法向應(yīng)力并顯示所施加的力apply nstress -3.0E7 range cyli

5、nder end1 0 -1 0 end2010 radius 39.9 not;柱體外表面施加法向力apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱體內(nèi)表面施加法向力plot grid green fap black;顯示施加的力;獲取某網(wǎng)格點的ID，坐標(biāo)值def aaget_id=gp_id(gp_near(3,0,0)endaaprint get_idprint gp position range id 3;保存命令流文件set log onset logfile hua.log;文件存放的路徑與調(diào)用

6、的文本路徑相同.set log off;fish函數(shù)的說明與解釋def find_add  ;定義fish函數(shù)find_addhead = null   ;給head賦值p_gp = gp_head  ;第一個網(wǎng)格結(jié)點的指針賦給p_gploop_while p_gp # null  ;當(dāng)p_gp值不為null時作循環(huán)x_pos = gp_xpos(p_gp)  ;將指針為p_gp的結(jié)點的x坐標(biāo)值賦給x_posif x_pos = 10.0 then  ;如果x_pos = 10.0 則(執(zhí)行)new = get_mem(2)

7、  ;從主內(nèi)存空間里得到2個fish變量對象并返回第一個對象的開始地址mem(new) = head  ;將head類型和數(shù)值置于地址為new的fish變量mem(new+1)=p_gp  ;第二個變量存儲符合條件的zone地址head = new       ;將new值賦給headendifp_gp = gp_next(p_gp)  ;將結(jié)點指針為p_gp的下一個結(jié)點的指針賦給p_gpendloop;結(jié)束循環(huán)end  ;結(jié)束fish函數(shù);實際上這個fish函數(shù)為滿足條件（x坐標(biāo)為1

8、0的）的zone的地址開辟一定的地址空間，各地址之間存在一定的聯(lián)系;找到第一個符合條件的zone地址后，用new = get_mem(2) 從主內(nèi)存空間里得到2個fish變量對象并返回第一個對象的開始地址，并用;下面的mem(new) = head  ，將第一次開辟的兩個變量的第一個變量存儲地址head（注意第一次headnull），第二個變量存儲第一;個符合條件的zone地址，并將第一個變量的地址賦予head（headnew），第一次循環(huán)結(jié)束；下次循環(huán)，同樣開辟兩個變量對象，第;一個變量對象記錄上次循環(huán)開辟的第一個變量的地址，第二個記錄第二個符合條件的zone地址，其余循環(huán)依次類推

9、，這樣子就建立了;一個符合條件的zone地址鏈條，方便以后使用。；FLAC后處理；輸出JPG圖片set plot jpgset plot quality 100plot hard file XX.jpg;輸出TXT文件hist write X file xx.txt  ；(x是歷史記錄ID）print con zd;輸出單元應(yīng)力和節(jié)點位移set log onset logfile xx.logprint zo stressprint gp disset log off；生成影片文件plot set rot 20 0 30plot con szz ou on magf 10plot a

10、dd hist 1set movie avi step 1 file xx.avimovie startsolve;計算間隔命令 在solve命令之前set hist_rep 1   ；(步數(shù)）；檢查模型賦值情況plot bcon prop bu;變形云圖加網(wǎng)格 放大20 倍plot con zd ou on magf 20;變形網(wǎng)格plot sk magf 20;變形矢量圖（加網(wǎng)格）plot sk dis scale 0.07  ；(scale是箭頭長度，默認(rèn)0.05）；如何觀察塑性區(qū)plot block stateplot block state sh_nplot block state sh_p;切片操作plot set plane ori 1 2 3 norm 0 0 -1plot con zd p