PFC2D学习笔记之边界条件与初始条件.pdfVIP

PFC2D 学习笔记之边界条件与初始条件 边界条件 墙边界 一般，颗粒组是在一组墙内创建并压缩。这些墙也可以当作边界约束，按一定速度运动 监测其反力；或保持反力一定，伺服控制其速度。但是不能在墙上直接施加力。 当墙为多段线时，要考虑到墙与球在线段转角处的接触。当转角处为凹时，两段墙可能会同 时与一个球接触，但 PFC2D 只允许每个墙与球的接触最多一个，因此要在此处将墙断开， 设置成两个墙；当转角为凸时，两段墙不肯能同时与一个球接触，因此不需要将墙断开。墙 也可以是圆形或弧或点。 颗粒边界 可以创建一个颗粒串，并使用这些颗粒串作为边界条件。 固定速度的颗粒边界 用 FISH 函数得到边界颗粒，并固定其平动自由度，然后删除墙，并施加加速度。从而可以 得到内部颗粒的速度。边界颗粒的速度保持初始值不变；该模拟对应与应变控制式试验。有 两种方法可以得到边界颗粒。一种是，遍历每个球的接触列表并探测与墙接触的球，作为边 界；另一种是，使用 RANGE ，将落在一个空间范围内的所有颗粒作为边界颗粒。 第一种方法的核心程序为： bp = ball_head loop while bp # null ; scan all balls section cp = b_clist(bp) loop while cp # null ; scan ball’s contacts 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 ; identify by color 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) end_if end_loop end_section bp = b_next(bp) end_loop 第二种方法参考 RANGE 的用法，较为简单，此处略去。 在颗粒边界上施加力 先固定所有边界颗粒，并移除墙，然后 CYCLE 一次，这样每个边界颗粒上的不平衡力就与 之前的墙反力相等或相反。接着，用 FISH 函数在边界颗粒上施加与不平衡力相反的力，这 样，边界颗粒受力平衡。然后撤除所有颗粒的约束，对边界颗粒施加加速度。 但是为了防止由于体力的未平衡导致一些边界颗粒移动位置，就需要使用各种稳定方法。比 如，固定边界颗粒的旋转自由度，防止颗粒从平衡位置旋转开。 也可以使用混合边界条件。 与固定速度的颗粒边界不同的是，这种边界条件要使用 FISH 函数施加力，平衡掉颗粒边界 上的不平衡力，其核心程序为： bp = ball_head loop while bp # null if b_xfix(bp) = 1 b_xfap(bp) = -b_xfob(bp) b_yfap(bp) = -b_yfob(bp) end_if bp = b_next(bp) end_loop 混合边界条件 垂直于 X 轴方向的两个边界固定住并施加速度边界；其他边界施加力。当几何边界改变时， 这种双轴试验就不能继续了，因为边界力将不可用。应该使用 FISH 定义一个可以随几何边 界更新的边界力，参考《双轴试验》部分。 核心程序与前面相似，不再重复。 应用命令 free x y spin range x=0.2,9.8，只释放上下边界和内部颗粒的自由度。 初始条件 为了得到所需要的初始应力状态，需要调整颗粒组状态，因此必须在一定程度上改变孔隙比， 其结果是，孔隙比的很小的变化都会引起应力的很大变化。 这里的应力是指测量区域内的平均应力，有两种计算方法：一种是，将墙的反力之和除以墙 的面积(长度与厚度之积) ；另一种是使用一个或多个测量圆。 获得各向同性应力状态 1、生成所需要的紧密颗粒组，用半径扩展法或颗粒排斥法。其核心程序参考《颗粒生成细 节》 2 、定义一个测量圆，用于测量区域内平均应力。测量圆的大小将会影响计算结果。 meas id=1 x 5 y 2.5 rad 2.0 3、用内置的MP=FIND_MEAS （1）函数获取测量圆 1 的内存地址。 4 、设置所要求的应力大小及其精度，并适当 CYCL