2_Getting Started

离散单元法 PFC2D PFC3D 入门操作 ITASCA 武汉 咨询有限公司中国水电顾问集团 ITASCA研发中心 朱永生 一个简单的例子 模型概念图 一个简单的例子 Wall墙体生成wallid 1face 0 0 0 10 0 0 10 0 5 0 0 5 wallid 2face 0 0 0 0 0 5 0 2 5 0 2 0 wallid 3face 10 0 5 10 0 0 10 2 0 10 2 5 wallid 4face 10 2 5 10 2 0 0 2 0 0 2 5 wallid 5face 0 0 5 10 0 5 10 2 5 0 2 5 wallid 6face 0 0 0 0 2 0 10 2 0 10 0 0 一个简单的例子 球体颗粒生成genid 1 750rad0 15 0 20 x 0 10y 0 2z 5 0模型查看生成视图plotcreateFootingplotsettitletext Basiccollectionofparticlesinabox plotsetrot20 00 00 0plotsetbackwhiteplotaddwalllcyanplotaddballorangeplotaddaxesbrownplotshow 一个简单的例子 对象属性wall墙体属性wallid 1kn 1e8ks 1e8wallid 2kn 1e8ks 1e8wallid 3kn 1e8ks 1e8wallid 4kn 1e8ks 1e8wallid 5kn 1e8ks 1e8wallid 6kn 1e8ks 1e8ball颗粒属性propdensity2000kn1e8ks1e8接触状态生成及其计算控制膨胀法生成初始应力场iniradmul1 47 一个简单的例子 接触状态生成及其计算控制计算控制sethist rep5 模型响应监测迭代间隔histballyvel 3 1 0 监测某点y向速度histdiagnosticmuf 监测系统最大不平衡力setdtdscale 拟静力分析设置cyc3000模型结果查看printhistplothist1plothist2savefoot0 sav 一个简单的例子 建立真实初始应力场 重力场setgrav00 9 81propertyfric1 0cyc1000初始模型结果查看plotcreateContact Forceplotsettitletext ModelatEquilibrium plotsetbackwhiteplotaddwalllcyanplotaddballorangeplotaddcforceblackplotaddaxesbrownplotshow plotprintview plotprintitemSaveFoot1 sav 一个简单的例子 后期加载deletewall6wallid10face 0 0 0 0 2 0 2 2 0 2 0 0 wallid10kn1e8ks1e8fric1 0wallid10zvel 0 5e 3propertyn bond5e5s bond5e5地基承载力监测histwallzforceid10求解cycle2700 plothist3 Example Footing dat PFC数值分析一般性步骤 数值 网格 模型生成 ball generate 颗粒体的接触状态建立 一般采用膨胀法 接触模型设置及其参数赋值 边界条件和初始状态建立 后续工况分析 视图处理相关操作 数值 网格 模型建立 Creatingamixtureoftwodiscreteparticlesizes FNAME gen1 datnewsetrandommacrobig id1 20 macrosmall id51 100 genbigrad 8 8x 0 10y 0 10z 0 5gensmallrad 4 4x 0 10y 0 10z 0 5changec index0rangebigchangec index1rangesmallplotcreatefredplotsetbackwhiteplotaddballorangelgrayplotshow 2 Confiningparticlestoanangledregion FNAME gen2 datnewsetrandomwallid1orig000norm 101genno shadowgaussx05y010z05rad 2 22id 1 200delwall1wallid1orig000norm10 1genno shadowgaussx05y010z05rad 3 33id 201 300delwall1changec index0rangeid 1 200changec index1rangeid 201 300plotcreatethe viewplotsetbackgroundwhiteplotaddballorangelgrayplotaddaxesbrownplotshow Generatingahemisphericalregionofparticles FNAME gen3 datnewsetrandomgenx05y05z05rad 22 22 id 1 150annulus2 5 2 5 002 5plotcreatethe annulusplotsetbackwhiteplotaddballlgrayplotaddaxesmagentaplotshow 接触状态的建立 Example Expand dat 接触模型及其参数赋值 常用的两种接触强度准则Slipmodel frictionContactBond n bonds bondParallelBond pb radiuspb knpb kspb nstrengthpb sstrength 通过wall或者property命令赋值 wall knksfricBall knksfric 边界条件及其初始状态建立 荷载种类 均布荷载线性分布 边界条件建立 Example Raft apply dat 边界条件