1、PFC3D中文讲解例题通过例子学习几种常见命令例二:PFC3D目录下的GuideStartfooting.dat,程序如下(注意:本程序与 原程序不同,特加了 plot set rotation (24.0,352.0,340.0)这句以使得模型转过一定 角度,让用户看的更清楚, 另外还将原程序一分为二, 主要是为了更好地学习每 个命令的作用),以下是程序代码。;fname: footing.DAT (tutorial example for PFC3D) new ; clear program state to begin new problemset random ; reset rand
2、om-number generator title Tutorial Examplewall 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,
3、-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,0Jplot create Footingplot set title text Basic collection of particles in a boxplot set rotation (24.0,352.0,340.0)plot add ball yellowplot add wall white plot
4、 add axes brown plot show pausewall id=1kn=1e8ks=1e8wall id=2kn=1e8ks=1e8wall id=3kn=1e8ks=1e8wall id=4kn=1e8ks=1e8wall id=5kn=1e8ks=1e8wall id=6kn=1e8ks=1e8prop density 2000 kn 1e8 ks 1e8ini rad mul 1.47 set hist_rep=5 hist ball zvel 3 1 0 hist diagnostic muf set dt dscalecycle 3000以下是程序中逐条命令的解释:1、
5、我们首先定义了模型的边界,使用 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,
6、 2,-5) ( 0, 2,-5)wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0)这些命令创造了一个由六个墙面围成的箱子, 每个墙面都赋予了一个独一无二的ID号。每个墙面四边形通过给出四边形的四个顶点坐标来定义,即四个顶 点连成四条线,然后再围成一个四边形墙面。注意:顶点的顺序决定了墙面是否 被激活,激活的墙面可以自动识别与小球的接触,不激活的话就会被小球穿透,2、 然后在箱子里填充小球,用 Gen erate命令,也可简写,如下:gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0这个
7、命令规定了小球填充的空间范围,小球的大小以及填充的总个数。Id=1,750的含义是:ID号从1到750其实就生成了 750个小球,并给每个小 球指定了各自的ID号。Rad 0.15,0.20的含义是:小球半径从 0.15m到0.20m随机选择x=0,10 y=0,2 z=-5,0 的含义是:小球填充的空间是 0vx10, 0y2,-5zSJiny-compLiniii of jng jlar velocirj of uill abau匚 Loiter of rcuiionXVQlOvityyklochyjf*jllEl Hntircf rcijnMpimomponc nc AAguliF ml
8、 ill About iNcntcraf raw ionf wiillBill CandiiianRXIfiwed S ball畑inrii:saijrielociT fixsdibrban、尸wH曲fix詞 r tll艸in-ai01IJT lMK) fiwdftflltll 叫通 lot i 寸 fi suf J fbr Ki 1c&piht.-znsu lj 1. cIcciLv fixed: for hdtlINHUU-rtnpoflcm of icnce 祁* ied v billMmm曲ml i-f rnormiu unl-Usmldcityl of vcloc: JLppli&I
9、QtMlIy-Oorfipor1iil of fdhit appl 圧d ballyinoffij -wnipeinanl 址m g_?nlLjj 口 billyvelocimy-campariEfil ofty applied, to bjlldorCGi-compontnc of fiirre appl ifd w hallETT0TTfnc amping eic of 口口rngm口卫1區4 中 bjl llzvfil&c t|fcumporienc 住f verity applied re bil 1SETgnvttyjrutiiljzt gmity1、给墙面施加速度和角速度Wall id 1 xvel 1e-3给ID号为1的墙面施加x方向速度,大小为1e-3Wall id 2 yspin 1e-3 x=2 y=2 z=0给ID号为2的墙面施加绕y轴转动的角速度,旋转中心坐标为( 2,2,0)2、给球施加恒定不变的速度或角速度为球施加速度或角速度,可用命令 initial,change, property,但这种定义只能对球的速度进行初始化,不能保证以后球一直以恒定的速度运动下去,因此需要用 到fix命令。程序如下:Ini xvel 5e-3 zsp in 1e-3 range id
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1