PFC2D培训课件课件.ppt

wall id=4 nodes ( 4.0, 4.0) (-2.0, 4.0) wall id=5 nodes (-2.0, 4.0) (-2.0,-2.0) wall id=1 kn=1e8 ks=1e8 wall id=2 kn=1e8 ks=1e8 wall id=3 kn=1e8 ks=1e8 wall id=4 kn=1e8 ks=1e8 wall id=5 kn=1e8 ks=1e8 plot create the_view plot add ball yellow plot add axes black plot add wall blue id=on plot show pause property rad mul 1.5 plot add cf green cycle 500 4.边界条件 PFC2D中有三种边界条件，分别是：墙体边界、颗粒体边界和混合边界。其中颗粒体边界又分为速度边界和受力边界。 1).墙体边界 建模过程中，墙体可作为颗粒体的生成范围约束，但同时也可以将墙体作为边界来施加约束。对于墙体，我们只能施加速度约束，而不能直接对其施加外力，因为运动定律对墙体是不适用的。其速度由以下三个参数控制:线速度、角速度和旋转中心。墙的运动是通过不断更新定义墙的基点的位置来描述。采用WALL命令设置，如： wall id=1 x=1.0 y=1.0 spin=10.0 (2)颗粒体边界 PFC2D中的模型可以将一连串的颗粒体作为边界条件。基本方法：在模型紧密压缩至平衡后，我们通过FISH函数将与墙体相接触的颗粒体逐个提取，将这一系列的颗粒体采用共同的边界条件限制，最后删除初始的限制墙体，即实现了以颗粒体代替墙体来作为边界条件。颗粒体边界条件分速度边界和外力边界，下例为速度边界程序实现。 restore expand.sav def bound bp = ball_head loop while bp # null section cp = b_clist(bp) loop while cp # null if c_nforce(cp) # 0.0 then b2 = c_ball2(cp) if pointer_type( b2 ) = 101 then ; b2 is a wall b_xfix(bp) = 1 ; fix original ball in x,y b_yfix(bp) = 1 b_color(bp) = 1 exit section ; all done for this ball end_if end_if if c_ball1(cp) = bp cp = c_b1clist(cp) else cp = c_b2clist(cp) 程序6： end_if end_loop end_section bp = b_next(bp) end_loop end bound ini xvel=0.0 grad=(-0.1, 0) ini yvel=0.0 del wall 1 del wall 2 del wall 3 del wall 4 plot create assembly plot add ball red green plot add vel black plot show plot add cf white cycle 200 放大 放大 颗粒体速度边界示意图 其前期操作与速度边界类似，差异处是先将边界颗粒体的速度初始化为零，再删除限制墙体，采用FISH函数反向施加颗粒体所受的不平衡力，运行一定的计算步至平衡态。在此过程中，区域内部的颗粒体会因轻微的扰动而导致其自身产生运动而脱离边界颗粒体，此时前面将边界颗粒体的速度固定为零。 颗粒体外力边界 ini xvel=0 yvel=0 spin=0 del wall 1 del wall 2 del wall 3 del wall 4 cycle 1 def replace ; Replace unbalanced forces with applied forces b