FLAC3D 高手總結(jié).doc

;模型鏡像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一點(diǎn)）;沿X軸鏡像,通過(guò)對(duì)稱平面法線向量確定對(duì)稱面gen zone reflect normal 0 0 -1 ;沿z軸鏡像;繪圖控制pl contour szz outline on;在模型中顯示位移-應(yīng)變曲線hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y軸,n代表x軸(不管m,n的正負(fù)); -表示對(duì)其值作mirror;對(duì)模型進(jìn)行壓縮實(shí)驗(yà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;顯示云圖的同時(shí)也顯示模型網(wǎng)格輪廓plot add cont disp outline on ;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 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 cylinder 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)格點(diǎn)的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)用的文本路徑相同.set log off;fish函數(shù)的說(shuō)明與解釋def find_add ;定義fish函數(shù)find_addhead = null ;給head賦值p_gp = gp_head ;第一個(gè)網(wǎng)格結(jié)點(diǎn)的指針賦給p_gploop_while p_gp # null ;當(dāng)p_gp值不為null時(shí)作循環(huán)x_pos = gp_xpos(p_gp) ;將指針為p_gp的結(jié)點(diǎn)的x坐標(biāo)值賦給x_posif x_pos = 10.0 then ;如果x_pos = 10.0 則(執(zhí)行)new = get_mem(2) ;從主內(nèi)存空間里得到2個(gè)fish變量對(duì)象并返回第一個(gè)對(duì)象的開(kāi)始地址mem(new) = head ;將head類型和數(shù)值置于地址為new的fish變量mem(new+1)=p_gp ;第二個(gè)變量存儲(chǔ)符合條件的zone地址head = new ;將new值賦給headendifp_gp = gp_next(p_gp) ;將結(jié)點(diǎn)指針為p_gp的下一個(gè)結(jié)點(diǎn)的指針賦給p_gpendloop;結(jié)束循環(huán)end ;結(jié)束fish函數(shù);實(shí)際上這個(gè)fish函數(shù)為滿足條件（x坐標(biāo)為10的）的zone的地址開(kāi)辟一定的地址空間，各地址之間存在一定的聯(lián)系;找到第一個(gè)符合條件的zone地址后，用new = get_mem(2) 從主內(nèi)存空間里得到2個(gè)fish變量對(duì)象并返回第一個(gè)對(duì)象的開(kāi)始地址，并用;下面的mem(new) = head ，將第一次開(kāi)辟的兩個(gè)變量的第一個(gè)變量存儲(chǔ)地址head（注意第一次headnull），第二個(gè)變量存儲(chǔ)第一;個(gè)符合條件的zone地址，并將第一個(gè)變量的地址賦予head（headnew），第一次循環(huán)結(jié)束；下次循環(huán)，同樣開(kāi)辟兩個(gè)變量對(duì)象，第;一個(gè)變量對(duì)象記錄上次循環(huán)開(kāi)辟的第一個(gè)變量的地址，第二個(gè)記錄第二個(gè)符合條件的zone地址，其余循環(huán)依次類推，這樣子就建立了;一個(gè)符合條件的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é)點(diǎn)位移set log onset logfile xx.logprint zo stressprint gp disset log off；生成影片文件plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie avi step 1 file xx.avimovie startsolve;計(jì)算間隔命令 在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是箭頭長(zhǎng)度，默認(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 pla