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

《FLAC/FLAC3D基础与工程实例》全部命令流1genzonebrip0000p11000p20100p30010&p415150p501515p615010p7202020&size101010rat1.00.91.1groupbrick_1genzonebrip02000p1add1000p2add0200p3add0015&size101010rat1.00.91.1groupbrick_2genzonebrip04000edge10size101010rat1.00.91.1groupbrick_3plotsur2-1newgenzonbrisize333modelelaspropbulk3e8shear1e8inidens2000fixzranz-.1.1fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solveappnstress-10e4ranz3x12y12histgpvel003histgpvel033plohist1redploaddhist2bluesolve3-1;-------------------------------------工程信息;ProjectRecordTreeexport;Title:Simpletest;---------------------------------计算第一步;...STATE:STATE1....configgrid10,10modelelasticgroup'User:Soil'notnullmodelelasticnotnullgroup'User:Soil'propdensity=1500.0bulk=3E6shear=1E6notnullgroup'User:Soil'fixxyj1fixxi1fixxi11setgravity=9.81history999unbalancedsolvesavestate1.sav;----------------------------------计算第二步;...STATE:STATE2....initialxdisp0ydisp0initialxvel0yvel0modelnulli47j810group'null'i47j810groupdelete'null'history1xdispi=4,j=11solvesavestate2.sav;--------------------------------绘图命令;***plotcommands****;plotname:syyplotholdgridsyyfill;plotname:Unbalancedforceplotholdhistory999;plotname:gridplotholdgridmagnify20.0lredgriddisplacement;plotname:Xdis-Aplotholdhistory1line5-1new;===============================;定义球体半径和半径方向上单元网格数;===============================defparmrad=10.0rad_size=5endparm;===============================;建立八分之一球体外接立方体网格;===============================genzonepyramidp0rad00p1rad0radp2radrad0p3000&p4radradradsizerad_sizerad_sizerad_sizegroup1genzonepyramidp00rad0p1radrad0p20radradp3000&p4radradradsizerad_sizerad_sizerad_sizegroup2genzonepyramidp000radp10radradp2rad0radp3000&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=sqrt(px*px+py*py+pz*pz)ifdist0then;节点位置调整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xpos(p_gp)=k*pxgp_ypos(p_gp)=k*pygp_zpos(p_gp)=k*pzend_ifp_gp=gp_next(p_gp)end_loopendmake_sphere;===============================;利用镜像生成完整球体网格;===============================genzonerefgenzonerefdip90genzonerefdip90dd90;===============================;显示球体网格;===============================plotsurfplsetbackwhplblgr5-3ngenzonbrisize112groupsoilranz11grouprockranz01expgrid1.flac3d6-1genzonbrisize333modelmohrpropbu3e6sh1e6coh10e3fric15fixzranz-.1.1fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1inidens2000histunbalsetgrav10solveelasticsave6-1.sav6-2rest6-1.savinixd0yd0zd0xv0yv0zv0appnstress-100e3ranz2.93.1x12y12solvesave6-2.sav6-3rest6-1.savinixd0yd0zd0xv0yv0zv0appnstress-100e3ranz2.93.1x12y12histid=2gpzdis113histid=3gpzdis112histid=4gpxdis113histid=5gpxdis113histid=6zoneszz113histid=7zoneszz1.51.52.5histid=8zonesxz1.51.52.5solvesave6-3.sav6-4rest6-3.savsetlogonsetlogfile6-2.logprintzonestressprintgpdissetlogoff6-5rest6-1.savinixd0yd0zd0xv0yv0zv0appnstress-100e3ranz2.93.1x12y12plotsetrot20030plotconszzouonmagf10plotaddhist1setmovieavistep1file6-5.avimoviestartsolvemoviefinish7-1ngenzonbrisize112modelelaspropbulk3e7shear1e7fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000setgrav00-10solveploconsz7-2ngenzonbrisize112modelmohrpropbulk3e7shear1e7c1e10f15tension1e10fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000setgrav00-10solvepropbulk3e7shear1e7c10e3f15ten0solveploconsz7-3ngenzonebricksize112modelmohrpropbulk3e7shear1e7coh10e3fri15ten0fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000setgrav00-10solveelasploconsz7-4newgenzonebricksize112modelmohrpropbulk3e7shear1e7coh10e3fri15ten0fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000iniszz-40e3grad0020e3ranz02inisyy-20e3grad0010e3ranz02inisxx-20e3grad0010e3ranz02setgrav00-10solveploconsz7-5ngenzonbrisize112modelmpropbulk3e7shear1e7c10e10f15ten1e10fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000ranz01inidens1500ranz12iniszz-35e3grad0020e3ranz01inisyy-17.5e3grad0010e3ranz01inisxx-17.5e3grad0010e3ranz01iniszz-15e3grad0015e3ranz12inisyy-7.5e3grad007.5e3ranz12inisxx-7.5e3grad007.5e3ranz12inipp10e3grad00-10e3ranz01setgrav00-10solveploconsz7-6ngenzonbrisize112modelmpropbulk3e7shear1e7c10e10f15ten1e10fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000ranz02iniszz-50e3grad0020e3ranz01inisyy-30e3grad0010e3ranz01inisxx-30e3grad0010e3ranz01inipp30e3grad00-10e3ranz02appnstress-10e3ranz2setgrav00-10solveploconsz7-7newgenzonebrickp0000p16000p20600p30090&p460600p506090p6600150p76060150&size6610modelelasprobulk10e10she10e10iniden2500applysxx-1e9grad001.1111111e7rangex-.1.1applysxx-1e9grad006.6666666e6rangex59.960.1applysyy-1e9grad008.3333333e6rangey-.1.1applysyy-1e9grad008.3333333e6rangey59.960.1applyszz-1e8grad008.3333333e5ranz0120setgrav00-10step30000inixdisp0ydisp0zdisp0inixvel0yvel0zvel0plocontszz7-8newgenzonebrickp0000p16000p20600p30090&p460600p506090p6600150p76060150&size6610modelelasprobulk10e10she10e10iniden2500inisxx-1e9grad001.11111