FLAC3D命令流

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

1、实用标准文案命令流文档大全实用标准文案文档大全实用标准文案1genzonebrip0000p11000p20100p30010&p415150p501515p615010p7202020&size101010rat1.00.91.1groupbrick_1genzonebrip02000p1add1000p2add0200p3add0015&size101010rat1.00.91.1groupbrick_2genzonebrip04000edge10size101010rat1.00.91.1groupbrick_3

2、plotsur2-1newgenzonbrisize333modelelaspropbulk3e8shear1e8inidens2000fixzranz-.1.1fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solveappnstress-10e4ranz3x12y12histgpvel003histgpvel033plohist1redploaddhist2bluesolve3-1;-------------------------

3、------------工程信息;ProjectRecordTreeexport;Title:Simpletest;---------------------------------计算第一步;...STATE:STATE1....configgrid10,10modelelasticgroup'User:Soil'notnullmodelelasticnotnullgroup'User:Soil'propdensity=1500.0bulk=3E6shear=1E6notnullgroup'User:Soil'fix

4、xyj1fixxi1fixxi11setgravity=9.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.0rad_size=5endparm文档大全实用标准文案;===============================;建立八分之一球体外接立方体网格;===============================genzonepyramidp0rad00p1rad0radp2radrad0p3000&p4radradradsizerad_sizerad_sizerad_sizegroup1gen

7、zonepyramidp00rad0p1radrad0p20radradp3000&p4radradradsizerad_sizerad_sizerad_sizegroup2genzonepyramidp000radp10radradp2rad0radp3000&p4radradradsizerad_sizerad_sizerad_sizegroup3;==================================;利用FISH语言将内部立方体节点调整到球面;===========================

8、=======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)ifdist>0then;节点位置