PFC2D学习笔记之加载-求解-模拟-结果分析
5星 · 超过95%的资源 需积分: 15 106 浏览量
2009-08-19
17:02:45
上传
评论 8
收藏 96KB PDF 举报 
PFC2D学习笔记之加载-求解-模拟-结果分析
加载方法
加载分为被动加载和主动加载。
被动加载，是指模型内部的改变引起的模型内的加载，这种作用是自动发生的，如开挖，重力加载也属于被动加
载的一种，因为重力总是在模拟开始时就设置了的。需要注意的是，SET GRAVITY 命令只是根据颗粒的重度和
半径给颗粒施加一个体力，颗粒间的接触力不受该命令影响。
主动加载，由外部施加的速度和力组成。
墙的控制
可以给墙施加速度，包括平动和转到。墙运动时所受到的力和力矩可以用 PRINT WALL 显示，或用 HISTORY
WA LL 跟踪，或用 FISH 函数对其操作。还可以用 FISH 函数创建一个伺服控制机制，给 WAL L 施加速度以使力
或力矩保持为常量。也可以定义 WA LL 的运动方程，使其动态行为像一个惯性实体一样。
颗粒受力控制——应力边界
可以用一组平衡力作用在边界颗粒上，以取代 WA L L，且这组力可以放大或缩小以表示侧应力的变化。但是，应
力控制式边界在颗粒材料内不具有代表性；应变控制式边界更更真实地体现材料内部约束；施加的力边界对于为
破坏的粘结材料几乎已经足够了。
颗粒速度的控制
一组边界颗粒可以用命令 FIX 和 INIITIAL XVEL,YVEL 固定并施加速度。但是需要用一下 FISH 函数获取控制颗
粒的总反力：
fxsum_lh = 0.0
fxsum_rh = 0.0
bp = ball_head
loop while bp # null
if b_xfix(bp) = 1
if b_xvel(bp) > 0.0
fxsum_lh = fxsum_lh + b_xfob(bp)
end_if
if b_xvel(bp) < 0.0
fxsum_rh = fxsum_rh + b_xfob(bp)
end_if
end_if
bp = b_next(bp)
end_loop
周期性边界
当需要使所求解问题在给定区域循环进行时，使用周期性边界能够减少颗粒数目以计算时
间。设置周期性边界：set periodic on 0.0 width 0.0 height
求解程序
PFC2D 模型的求解即执行若干时间步（CYCLE 或 STEP）并观察结果行为，对于静态问题，是指系统达到平衡
资源评论
guochaoxu2014-09-03笔记不错，浅显易懂，解决了我很多的疑惑
