PFC2D&PFC3D――命令简介

离散单元法-PFC2D&PFC3D——命令简介Yongsheng(Freddie)ZhuItascaConsultingChinaLtd.ItascaSoftwareTrainingCourseRangekeyword……•anuuluscenterxyradiusr1r2•circlecenterxyradiusr•coloriliu•fishfname•groupgname•idiliu•jsetiliu•linekeyword:above,below,dip,distance,normal,origin•namername(rangename)•radiusrlru:ballswithradiusbetweenrlandruwillbeselected.•xxlxu•yylyuPFC数值分析常用命令数值分析一般步骤建立数值’网格’模型本构选择及其材料属性定义边界条件，初始条件工况设置后处理建立数值“网格”模型创建数值模型单元•创建PFC数值模型的命令：•Ball-----创建单个球体单元•Generate-----创建单个球体单元BallBALLradiusrkeyword...Keyword:density密度；hertz创建Hertz接触模型；id球体的标识号；x球体形心x坐标(default:x=0.0)y球体形心y坐标(default:y=0.0)该方法生成的球体单元可以在空间上有重叠。GenerateGeneratekeyword…xxlxuyylyuradiusrlruidiliuannulusxcycr1r2Keyword:filterfname-----通过fish函数确定在空间上生成颗粒的位置；guass-----颗粒半径服从高斯分布，数学期望E=(rl+ru)/2，方差sigma=(ru-rl)/2hertz-----指定接触模型为hertz模型；minrmin-----颗粒半径guass分布时，颗粒最小半径；no_shadow-----以wall作为辅助工具，颗粒生成定位；tries-----随机数生成控制。Example:Filter.datJsetJsetididkeywordkeyword……range…Keyword:area_ratiomddipmd结构面倾角，以+x为起始，顺时针为正,m为数学期望，d为方差；guss结构面倾角为高斯分布；number生成几条结构面；origin参考点radiusmd结构面小段长度一半；spacingmd结构面间距。MeasureMeasureididxxcyycradiusr设置后处理测量圈本构选择及其材料属性定义PropertyPropertykeywordaddormultiplyvgradientgxgykeywordv…range…BallPropertiescolor----ball的颜色damp----局部阻尼系数density----密度[质量/体积]friction----球体表面的摩擦系数(注：不是摩擦角！)kn----法向刚度(线性接触模型)[力/位移]ks----剪切刚度(线性接触模型)[力/位移]moment----施加力矩[力*距离]poiss-----泊松比(Hertz接触模型radius-----球体半径PropertyBallPropertiesshear----剪切模量(Hertz接触模型)spin----转动角速度xdisplacement----累计x方向位移xforce----施加于球体形心的x方向作用力xvelocity----x向速度ydisplacement----累计y方向位移yforce----施加于球体形心的y方向作用力yvelocity----y向速度PropertyContact-BondPropertiesn_bond----法向接触强度s_bond----切向接触强度Parallel-BondPropertiespb_kn---法向刚度[应力/位移]pb_ks---剪切刚度[应力/位移]pb_nstrength----法向接触强度[应力]pb_radius----接触约束单元半径pb_sstrength---切向接触强度[应力]Wall见边界条件与初始条件章节wall命令介绍。边界条件及其初始条件WallWallkeyword……Keyword:frictionf----摩擦系数；idid----标识号；knkn----法向刚度；ksks----剪切刚度；maxfricon/off------ball-wall接触时，摩擦系数取两者摩擦属性的最大值/最小值，defualt=off；nodesx1,y1x2,y2…xn,yncloseWallkeyword……Keyword:propetykeywordvaluekeywordvalue…frictionf-----摩擦系数；knkn-----法向刚度；ksks-----剪切刚度；maxfricon/offspinrv-----转动角速度；temperature----温度;xx-----转动参考点x坐标;xvelocityxv-----平动速度…..spinrv-----转动角速度；Wallkeyword……Keyword:typekeywordcirclekeywordanglea1a2------该类型wall的起始角和终止角；centerxy------圆心坐标；radr-------半径；pointkeywordcenterpxpyFix功能：控制ball的自由度（速度和角速度，而不是位移）Fixkeword…range…Keyword:spin----转动角速度；x----x向平动速度；y----y向平动速度。fixxrangeid5propertyxvel=1.5rangeid=5Initial功能：该命令实现的功能与property命令一致。SetGravity功能：设置重力加速度方向。Setgravitygxgy