PFC2D&PFC3D——入门操作

离散单元法-PFC2D&PFC3D——入门操作Yongsheng(Freddie)ZhuItascaConsultingChinaLtd.ItascaSoftwareTrainingCourse一个简单的例子模型概念图一个简单的例子•模型设置setdiskon•Wall墙体生成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•模型查看生成视图plotcreateFootingsettitletext‘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；•颗粒体的接触状态建立：一般采用膨胀法；•接触模型设置及其参数赋值；•边界条件和初始状态建立；•后续工况分析；•视图处理相关操作。数值’网格‘模型建立1.Creatingamixtureoftwodiscreteparticlesizes2.Confiningparticlestoanon-rectangularregion3.Generatingasemi-circularregionofparticlesExample:Gen1.datExample:Gen2.datExample:Gen3.dat接触状态的建立iniradmul1.60;iniradmul1.66Example:Expand.dat接触模型及其参数赋值常用的两种接触模型1.ContactBond：n_bonds_bond2.ParallelBond：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