PFC3D中文讲解.doc

（5）通过例子学习几种常见命令 例二：PFC3D目录下的Guide\Start\footing.dat，程序如下（注意：本程序与原程序不同，特加了plot set rotation (24.0,352.0,340.0)这句以使得模型转过一定角度，让用户看的更清楚，另外还将原程序一分为二，主要是为了更好地学习每个命令的作用），以下是程序代码。 ;fname: footing.DAT (tutorial example for PFC3D) new ; clear program state to begin new problem set random ; reset random-number generator title Tutorial Example wall id=1 face ( 0, 0, 0) (10, 0, 0) (10, 0 -5) ( 0, 0,-5) wall id=2 face ( 0, 0, 0) ( 0, 0,-5) ( 0, 2,-5) ( 0, 2, 0) wall id=3 face (10, 0,-5) (10, 0, 0) (10, 2, 0) (10, 2,-5) wall id=4 face (10, 2,-5) (10, 2, 0) ( 0, 2, 0) ( 0, 2,-5) wall id=5 face ( 0, 0,-5) (10, 0,-5) (10, 2,-5) ( 0, 2,-5) wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0) gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0 ; plot create Footing plot set title text Basic collection of particles in a box plot set rotation (24.0,352.0,340.0) plot add ball yellow plot add wall white plot add axes brown plot show pause ; 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 wall id=6 kn=1e8 ks=1e8 ; prop density 2000 kn 1e8 ks 1e8 ini rad mul 1.47 ; set hist_rep=5 hist ball zvel 3 1 0 hist diagnostic muf set dt dscale cycle 3000 以下是程序中逐条命令的解释： 1、我们首先定义了模型的边界，使用WALL命令，如下： wall id=1 face ( 0, 0, 0) (10, 0, 0) (10, 0 -5) ( 0, 0,-5) wall id=2 face ( 0, 0, 0) ( 0, 0,-5) ( 0, 2,-5) ( 0, 2, 0) wall id=3 face (10, 0,-5) (10, 0, 0) (10, 2, 0) (10, 2,-5) wall id=4 face (10, 2,-5) (10, 2, 0) ( 0, 2, 0) ( 0, 2,-5) wall id=5 face ( 0, 0,-5) (10, 0,-5) (10, 2,-5) ( 0, 2,-5) wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0) 这些命令创造了一个由六个墙面围成的箱子，每个墙面都赋予了一个独一无二的ID号。每个墙面四边形通过给出四边形的四个顶点坐标来定义，即四个顶点连成四条线，然后再围成一个四边形墙面。注意：顶点的顺序决定了墙面是否被激活，激活的墙面可以自动识别与小球的接触，不激活的话就会被小球穿透， 2、然后在箱子里填充小球，用Generate命令，也可简写，如下： gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0 这个命令规定了小球填充的空间范围，小球的大小以及填充的总个数。 Id=1,750的含义是：ID号从1到750其实就生成了750个小球，并给每个小球指定了各自的ID号。 Rad 0.15,0.20的含义是：小球半径从0.15m到0.20m随机选择 x=0,10 y=0,2