离散单元法-PFC2D&PFC3D——入门操作Yongsheng(Freddie)ZhuItascaConsultingChinaLtd.ItascaSoftwareTrainingCourse一个简单的例子模型概念图一个简单的例子模型设置setdiskonWall墙体生成wallid=1nodes(0,-5)(10,-5)wallid=2nodes(10,-5)(10,0)wallid=3nodes(10,0)(0,0)wallid=4nodes(0,0)(0,-5)一个简单的例子球体颗粒生成genid=1,500rad0.08,0.13x=0,10y=-5,-0.5模型查看生成视图plot>createFootingsettitletext‘BaseCollectionofParticlesinABox’addballyellowaddwallblackaddaxesbrownshow视图调整setsize-1.01.0-2.02.0setcaptionright一个简单的例子对象属性wall墙体属性wallid=1kn=1e8ks=1e8wallid=2kn=1e8ks=1e8wallid=3kn=1e8ks=1e8wallid=4kn=1e8ks=1e8ball颗粒属性propdensity1000kn1e8ks1e8接触状态生成及其计算控制膨胀法生成初始应力场propradmul1.51一个简单的例子接触状态生成及其计算控制计算控制sethist_rep5-----模型响应监测迭代间隔histballyvel(3,0)-----监测某点y向速度histdiagnosticmuf-----监测系统最大不平衡力setdtdscale-----拟静力分析设置cyc4000模型结果查看printhistplothist0plothist1……一个简单的例子建立真实初始应力场---重力场setgrav0-9.81propertyfric1.0cyc1000初始模型结果查看plotcreateContact_Forceplotsettitletext‘ModelatEquilibrium’plotaddballlblueplotaddwallblackplotaddcforceblackplotaddaxesbrownplotshow;;;plotprintview,plotprintitemSaveFoot1.sav一个简单的例子后期加载deletewall3wallid10nodes(2,0)(0,0)wallid10kn1e8ks1e8fric1.0wallid10yvel-0.5e-3propertyn_bond5e5s_bond5e5地基承载力监测histwallyforceid10求解cycle2000;;;;plothist3Example:Footing.datPFC数值分析一般性步骤数值’网格‘模型生成:ball,generate;颗粒体的接触状态建立:一般采用膨胀法;接触模型设置及其参数赋值;边界条件和初始状态建立;后续工况分析;视图处理相关操作。数值’网格‘模型建立CreatingamixtureoftwodiscreteparticlesizesConfiningparticlestoanon-rectangularregionGeneratingasemi-circularregionofparticlesExample:Gen1.datExample:Gen2.datExample:Gen3.dat接触状态的建立iniradmul1.60;iniradmul1.66Example:Expand.dat接触模型及其参数赋值常用的两种接触模型ContactBond:n_bonds_bondParallelBond:pb_radiuspb_knpb_kspb_nstrengthpb_sstrength通过wall或者property命令赋值。边界条件及其初始状态建立荷载种类:均布荷载线性分布Example:Raft.dat边界条件建立边界条件及其初始状态建立初始状态建立Example:Tun1.dat后续工况分析Example:Tun2.dat后续工况分析Example:Tun_Unl.dat无支护隧洞后续工况分析Example:Tun_lin.dat支护隧洞视图处理Example:Trench1.dat