PFC3D中文講解例題

1、.通過例子學(xué)習(xí)幾種常見命令例二：pfc3d目錄下的guidestartfooting.dat，程序如下（注意：本程序與原程序不同，特加了plot set rotation (24.0,352.0,340.0)這句以使得模型轉(zhuǎn)過一定角度，讓用戶看的更清楚，另外還將原程序一分為二，主要是為了更好地學(xué)習(xí)每個命令的作用），以下是程序代碼。;fname: footing.dat (tutorial example for pfc3d)new ; clear program state to begin new problemset random ; reset random-number generat

2、ortitle 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,-5) (10, 0,-5) (10, 2,

3、-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 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 whiteplot add axes brownplot showp

4、ause;wall id=1 kn=1e8 ks=1e8wall id=2 kn=1e8 ks=1e8wall id=3 kn=1e8 ks=1e8wall id=4 kn=1e8 ks=1e8wall id=5 kn=1e8 ks=1e8wall id=6 kn=1e8 ks=1e8;prop density 2000 kn 1e8 ks 1e8ini rad mul 1.47;set hist_rep=5hist ball zvel 3 1 0hist diagnostic mufset dt dscalecycle 3000以下是程序中逐條命令的解釋：1、我們首先定義了模型的邊界，使用w

5、all命令，如下：精品.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,-

6、5)wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0)這些命令創(chuàng)造了一個由六個墻面圍成的箱子，每個墻面都賦予了一個獨一無二的id號。每個墻面四邊形通過給出四邊形的四個頂點坐標(biāo)來定義，即四個頂點連成四條線，然后再圍成一個四邊形墻面。注意：頂點的順序決定了墻面是否被激活，激活的墻面可以自動識別與小球的接觸，不激活的話就會被小球穿透，2、然后在箱子里填充小球，用generate命令，也可簡寫，如下：gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0這個命令規(guī)定了小球填充的空間范圍，小球的大小以及填充的總個數(shù)。id=1,750的含義是：id號從1到750其實就生成了750個小球，并給每個小球指定了各自的id號。rad 0.15,0.20的含義是：小球半徑從0.15m到0.20m隨機(jī)選擇x=0,10 y=0,2 z=-5,0的含義是：小球填充的空間是0x10，0y2，-5z0 thenabc=11elseabc=2