《flac3d基础与工程实例》全部命令流

《《flac3d基础与工程实例》全部命令流》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《FLAC/FLAC3D基础与工程实例》全部命令流1genzonebrip0000p11000p20100p30010&p415150p501515p615010p7202020&size101010rat1.00.91.1groupbrick_1genzonebrip02000p1add1000p2add0200p3add0015&size101010rat1.00.91.1groupbrick_2genzonebrip04000edge10size101010rat1.00.91.1groupbrick_3plotsur2-1newgenzonb

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

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

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

5、otholdgridsyyfill;plotname:Unbalancedforceplotholdhistory999;plotname:gridplotholdgridmagnify20.0lredgriddisplacement;plotname:Xdis-Aplotholdhistory1line5-1new;===============================;定义球体半径和半径方向上单元网格数;===============================defparmrad=10.0rad_size=5endparm;====

6、===========================;建立八分之一球体外接立方体网格;===============================genzonepyramidp0rad00p1rad0radp2radrad0p3000&p4radradradsizerad_sizerad_sizerad_sizegroup1genzonepyramidp00rad0p1radrad0p20radradp3000&p4radradradsizerad_sizerad_sizerad_sizegroup2genzonepyramidp000radp1

7、0radradp2rad0radp3000&p4radradradsizerad_sizerad_sizerad_sizegroup3;==================================;利用FISH语言将内部立方体节点调整到球面;==================================defmake_spherep_gp=gp_headloopwhilep_gp#null;获取节点点坐标值：P=(px,py,pz)px=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=

8、sqrt(px*px+py*py+pz*pz)ifdist>0then;节点位置调整maxp=max(px,