FLAC3D命令流.doc

《FLAC3D命令流.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.命令流....1genzonebrip0000p11000p20100p30010&p415150p501515p615010p7202020&size101010rat1.00.91.1groupbrick_1genzonebrip02000p1add1000p2add0200p3add0015&size101010rat1.00.91.1groupbrick_2genzonebrip04000edge10size101010rat1.00.91.1groupbrick_3plotsur2-1newgenzonbrisi

2、ze333modelelaspropbulk3e8shear1e8inidens2000fixzranz-.1.1fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solveappnstress-10e4ranz3x12y12histgpvel003histgpvel033plohist1redploaddhist2bluesolve3-1;-------------------------------------工程信息;ProjectRec

3、ordTreeexport;Title:Simpletest;---------------------------------计算第一步;...STATE:STATE1....configgrid10,10modelelasticgroup'User:Soil'notnullmodelelasticnotnullgroup'User:Soil'propdensity=1500.0bulk=3E6shear=1E6notnullgroup'User:Soil'fixxyj1fixxi1fixxi11setgravity=9.

4、81history999unbalancedsolvesavestate1.sav;----------------------------------计算第二步;...STATE:STATE2....initialxdisp0ydisp0initialxvel0yvel0modelnulli47j810group'null'i47j810groupdelete'null'history1xdispi=4,j=11solvesavestate2.sav;--------------------------------绘图命令

5、;***plotcommands****;plotname:syyplotholdgridsyyfill;plotname:Unbalancedforceplotholdhistory999;plotname:gridplotholdgridmagnify20.0lredgriddisplacement;plotname:Xdis-Aplotholdhistory1line5-1new;===============================;定义球体半径和半径方向上单元网格数;====================

6、===========defparmrad=10.0..rad_size=5endparm;===============================;建立八分之一球体外接立方体网格;===============================genzonepyramidp0rad00p1rad0radp2radrad0p3000&p4radradradsizerad_sizerad_sizerad_sizegroup1genzonepyramidp00rad0p1radrad0p20radradp3000&p4rad

7、radradsizerad_sizerad_sizerad_sizegroup2genzonepyramidp000radp10radradp2rad0radp3000&p4radradradsizerad_sizerad_sizerad_sizegroup3;==================================;利用FISH语言将内部立方体节点调整到球面;==================================defmake_spherep_gp=gp_headloopwhilep_gp#nul

8、l;获取节点点坐标值：P=(px,py,pz)px=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)ifdist>0then;节点位置调整maxp=max(px,max(py,pz))k=(m