;
macro ’wall_fric’ ’0.5’ wall id 4 wall_props node pt1 pt2 ;bottom
new
gen id 1 400 x 0 5 y 0 5 rad 0.05 0.05 ini rad mul 1.4 ；定义一组范围的球体，名为Tunnel group Tunnel range x 1.5 3.5 y 0 5 plot create the_view plot add ball yellow plot show
pause
；删除先前定义的组
delete ball range group Tunnel
；只产生一个球，半径是必需的，如果ID未指定，则自动赋一个比现在ID大的数，如果x与y未指 定，则在坐标原点，指定参数hertz后，则质点将服从Hertz-Mindlin contact而不是默认的Linear law ，指定该法则后，则弹性体必需给出其属性，通过PROPERTY给出。产生的该质点可以与其它质 点重叠，而不像GENERATE命令，当执行循环后，重叠的质点将产生排斥力，除非给定边界限制 。
PFC2D学习初期总结[geotechman]
欢迎您： HBT2008 隐身 | 我的主页 | 控制面板 | 搜索 | 风格 | 展区 | 插件 | 退出 geotechman → 数值模拟 → PFC2D → PFC2D学习初期总结
左栏 右栏 我的收件箱 (1)
主题：PFC2D学习初期总结
总数 19 1 2 下一页 1 << 上一主题 下一主题 >> 共有2492人关注过本帖 树形 打印 收藏 修复 推荐
PFC2D学习初期总结[geotechman]
add ball lblue add wall black add axes brown show ；建立视点Footing2 create Footing2
add ball lblue
add wall black add axes brown set size -1.0 1.0 -2.0 2.0 set title bottom
hist 1 hist 2
；将显示当前的几个视点，如在plot 模式下。 print list ；如果在命令模式下，则如下： plot print list ；在命令模式下，如下的命令将显示所有的历史信息 print hist ；提供重力，加速度为9.81m/s2 set grav 0,-9.81 ；提供球之间的磨擦力,磨擦系数为1.0 property fric 1.0 cycle 1000
PFC2D学习初期总结[geotechman]
macro pt3 ’1.0 1.0’ macro pt4 ’0.0 1.0’ macro wall_props ’kn wall_n_stiff ks wall_s_stiff fric wall_fric’ wall id 1 wall_props node pt2 pt3 ;side wall id 2 wall_props node pt3 pt4 ;top
内存分配 当载入PFC程序后，程序将动态分配内存，随着模型的构建，内存相应地增大。Print memory 将显 示当前内存大小。 内存分配 当载入PFC程序后，程序将动态分配内存，随着模型的构建，内存相应地增大。Print memory 将显 示当前内存大小。 图形控制 在图形模式下，按键M为放大图形，缩小为Shift+M；箭头来移动；Ctrl+G为在彩图与黑白图形间 转换；Ctrl+C显示照相机对话框；Ctrl+L为cutting plane；Ctrl+R为恢复图形显示为默认值；Ctrl+Z 为用一矩形框来放大图形局部；F9为重画。 实例 new ；将球体看成一定厚度的盘子；Note that only mass properties are affected by the SET disk command set disk on ；下面为建立四个直线边壁，分别赋予1~4的编号。每一直线壁以两点连线建立。 wall id=1 nodes ( 0, -5) (10, -5) wall id=2 nodes (10, -5) (10, 0) wall id=3 nodes (10, 0) ( 0, 0) wall id=4 nodes ( 0, 0) ( 0, -5) ；下面为产生500个球，每个球的编号被赋予1~500，产生球的地方在0<x<10,-5<y<-0.5的范围内，
;球为白色但有影子 plot add ball white shade on
一． 粒子创建 两种创建粒子的命令：BALL and GENERATE BALL 创建单一的粒子，并指定其位置，其不受已存在粒子的制约，也就是说其创建的粒子
可以与已存在的粒子叠加。 GENERATE 创建一系列的粒子，其不可重叠。 有两种粒子集：规则的和不规则的。
property n_bond=5e5 s_bond=5e5 ;记录ID号为10的wall在Y方向的力 hist wall yforce id 10
cycle 2000 plot current 0 plot hist 3 ；拷贝Contact_ Force 到Foot_Load下，并定义标题为：Loaded footing plot copy Contact_ Force Foot_Load both current 4 set title text ’Loaded footing’ show ；显示所有的球 print ball ；显示所有的墙 print wall ；显示模型的信息 print info
wjwhappy
加好友 发短信
等级：岩土一级 帖子：32 积分：0 威望：0 精华：1 注册：8/24/2005 7:33:42 AM
个性首页 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC
小 大 1楼
PFC2D学习初期总结
Post By：3/11/2006 9:39:55 PM
最近学了一段时间PFC2D，由于做课题要用到它，时间又很紧，只是粗粗的学了一下，现将我的学习笔记抄录如下，比 较乱一些。以期对比我更初级的学习者有所帮助。
;查看当前视图的详细信息
plot print view
;查看当前视图的具体模型
plot print item
;改变当前视图为另一视图
plot current Footing
;显示任何已存在的视图
plot show <viewname>
；以下命令将壁3删除，并建立壁10，因密度缩放比例是开着的，所以其速度单位为m/step，如 果是动态分析的话，则速度单位为length/time delete wall 3 wall id 10 nodes (2,0) (0,0) wall id 10 ks 1e8 kn 1e8 fric 1.0 wall id 10 yvel = -0.5e-3 ；在接触面上提供法向强度normal strength和切向强度shear strength，即内聚力cohesion，单位为： N.当超过法向或切向强度后，bond就会失稳破坏，球之间的摩擦力即发生效果。
ball rad v id i x v y v <hertz> ；建立球的最大与最小半径相同；如果半径不同的话，则产生的半径是最小的或是最大的，随机分布的。如果指定gaus s后，则按高斯分布，即按(R1+R2)/2与 (R2-R1)/2确定半径。 gen big rad=.6,.6 x=0,10 y=0,10 gen small rad=.3,.3 x=0,10 y=0,10 ；将大球定义Index 为0，小球为1 change c_index 0 range big change c_index 1 range small ；定义的wall的光面是向上的，故而在定义球时的no_shadow即这光面上方。 wall id 1 nodes -10 -10 10 10 gen no_shadow gauss x 0 6 y 0 6 rad .2 .22 id=1,50 del wall 1 ；定义的wall的光面是向下的，故而在定义球时的no_shadow即这光面上方。 wall id 1 nodes 10 10 -10 -10 gen no_shadow gauss x 0 6 y 0 6 rad .3 .32 id=51 75
；给球赋属性，密度为1000kg/m3,刚度都为108N/m
prop density 1000 kn 1e8 ks 1e8
；给球的半径乘以一个系数，从而减小球与球间的空隙
prop rad mul 1.51
；每5步监测一下变量的变化，从而将变化的数据存入一数组，默认为每10步；两个变量被监测， 一个是圆点在（3，0）附近的球在Y方向的速度，一个是平均不不衡力mean unbalanced force(muf) 的诊断
set hist_rep = 5
history ball yvel (3, 0) history diagnostic muf ；因为这是一个静态分析，用密度缩放参数(density-scaling option)使问题最优化 set dt dscale cycle 4000 ；视图Footing中，图像每20步进行一次更新，更新间隔可以通过SET pinterval 命令改变 ；以下4条命令进入到视点0中，显示刚才设定的历史1的图形，在该例中是圆点在（3，0）附近的 球在Y方向的速度和hist 2的图形 plot current 0 ；当前视点
； 定义宏 macro wall_s_stiff ’1e6’ macro wall_n_stiff ’1e6’ macro wall_fric ’0.0’ macro pt1 ’0.0 0.0’ macro pt2 ’1.0 0.0’
/old/dispbbs.asp?boardid=49&Id=3606[2009-4-26 23:31:23]