PFC常用命令

New清除所有颗粒setrandom/随机产生颗粒setdiskon;treatballsasdisksofunitthicknesstitle’Idealizedmineblock-cavingprocess’为整个程序命名，在窗体上显示为“Jobtitle。。。”wallid1ks1e9kn1e9fric0.5nodes-20,020,0/生成墙wallid墙编号；ks切向刚度；kn法向刚度；fric摩擦系数；nodesx，y产生墙的起始及结束坐标genx-20,20y0,3id=20000,22550rad0.05,0.1/生成颗粒xa1,b2yb1,b2,产生颗粒范围；id颗粒编号；radr1,r2产生颗粒半径范围；iniradmul1.55rangeid=20000,22550ini初始条件设置；radmuln颗粒半径扩大倍数；id颗粒编号；plotcreatecave1产生新的窗体plotaddballlbluerange（id）添加显示颗粒，range确定范围，xy；id确定编号；lblue颗粒颜色；（颗粒颜色：black，blue，green，cyan蓝绿色，red，magenta紫红色，orange，brown，lgray灰色，dgray，lblue，lgreen，lcyan，lred，lmagenta，lorange，yellow，white）plotaddclump显示簇plotaddcfblack添加接触应力wallid6xvel5设置墙体速度，x正方向数值为正，x负方向数值为负。wallid6spin5xa1yb1墙体以（a1,b1）为中心，以5rad/s的角速度旋转，正值为顺时针旋转，负值为逆时针旋转histid120wallpowerid6监测墙体与球接触力的大小histid=1diagmcf监测平均接触力histid=2diagmuf监测平均不平衡力measid1x-14y6rad6设置测圆，measure；id，编号；xy，测圆中心坐标；rad，测圆半径histid=3measureporosityid=1监测测圆孔隙度，id=3监测编号；id=1测圆编号propks1e9kn1e9dens2000设置颗粒材料参数，kskn，切向法向刚度；dens材料密度单位kg/m3setgravity0,-10设置重力加速度cyc6000循环次数，或者step命令savecave0.savsavecave0.dat保存文件Callabc.txt调用记事本编辑的程序Restoreabc.savRestoreabc.dat调用保存的数据文件deleteballsrangeid=20000,22550删除颗粒range确定范围histid=101ballyposition04监测单个颗粒y方向位移，此颗粒质心坐标距给定坐标（0，4）最近setlogonprinthistsetlogoff记录文件及其数据存取的过程，并在与原文件相同的文件夹里产生pfc2d.log记事本文件rangenamemeicengx0500y2030自定义名称代替范围macroraft'id=5000,5079'propertypb_kn1e10pb_ks1e10c_index1rangeraft自定义任意字符，代替常用或重复使用的字符串Solve自动运行至合适的结果deletewall3删除墙propertypb_rad1.0pb_s1e20pb_n1e20rangeraft设置平行粘结半径，切向法向刚度inixv0rangeRaftiniyforce-4e6.98rangeRaftx方向速度，y方向应力setechooffsetx0=-19.75y0=12.25radius=0.25/确定初始颗粒质心坐标及颗粒半径setid_start=5000n_col=80n_row=1/所产生颗粒的起始id编号，每排颗粒数及排数chanshengsetechoon产生规则排列的颗粒。defchanshengxc=x0yc=y0rc=radiusidc=id_startr2=2.0*radiusyinc=radius*sqrt(3.0)looprow(1,n_row)loopcol(1,n_col)commandballid=idcx=xcy=ycrad=rcend_commandidc=idc+1xc=xc+r2end_loopyc=yc+yincxc=x0+radius*(row-(row/2)*2)end_loopend产生规则颗粒的fish程序ballx0y2rad1产生单个颗粒，质心坐标，半径ballx0y2rad1ballx2y2rad1ballx4y2rad1clumpid=1rangex=-15y=1,3产生一系列单个颗粒，将这些颗粒连接成簇，xy分别限定簇起作用的范围deftraverse_clumpsclp=clump_headloopwhileclp#nulloo=out('Clump,ID='+string(cl_id(clp))+'contains...')bp=cl_list(clp)loopwhilebp#nulloo=out('Ball'+string(b_id(bp)))bp=b_cllist(bp)endLoopclp=cl_next(clp)endLoopend未知propxdisp0.0ydisp0.0将xy位移归零，进行初始化plotadddispred显示位移向量为红色plotaddmeasred显示测圆为红色histid911measporosityid1/孔隙率histid912meass11id1/x-xstresshistid913meass12id1/x-ystresshistid914meass22id1/y-ystressjsetid1origin030dip45spacing10number100rangeguanxia产生节理，originxy初始点坐标；dip节理方向角度；spacing节理间间距；number产生节理数目；range产生范围plotaddcontactyellowrangejset12显示节理颜色，jset编号propkn3e8ks3e8dens2000pb_kn1e9pb_ks1e9pb_nstr1e15pb_sstr1e15rangeid1100设置材料参数