PFC3D5.0球体规则排列操作实例代码

John.Lee

1. ; fname: ball_gen.p3dat
   
2. ; 问题描述：按照矩形排列生成9*9的颗粒平面；按照六边形排列生成9*9颗粒平面。
   
3. new
   
4. domain extent -40 40
   
5. 
   
6. def orthogon_arrange
   
7.     local ini_x = -21.0
   
8.     local ini_y = 0.0
   
9.     local ini_z = 5.0
   
10.     local ini_rad = 1.0
    
11.     loop local n(1,9)
    
12.         o_x = ini_x
    
13.         o_y = ini_y
    
14.         o_z = ini_z + 2.0*(n-1)*ini_rad
    
15.         loop local m(1,9)
    
16.             command
    
17.                 ball create group arrange_one ...
    
18.                      radius @ini_rad x @o_x y @o_y z @o_z
    
19.             ;pause key
    
20.             end_command
    
21.             o_x = ini_x + 2.0*ini_rad*m
    
22.         end_loop
    
23.     end_loop
    
24. end
    
25. @orthogon_arrange
    
26. 
    
27. def hexagon_arrange
    
28.     local ini_x = 5.0
    
29.     local ini_y = 0.0
    
30.     local ini_z = 5.0
    
31.     local ini_rad = 1.0
    
32.     loop local n(1,9)
    
33.         h_y = ini_y
    
34.         h_z = ini_z + (n-1)*ini_rad*math.sqrt(3.0)
    
35.         test_n = n/2
    
36.         if (test_n*2 - n) < 0 then
    
37.         h_x = ini_x
    
38.         else
    
39.         h_x = ini_x + ini_rad
    
40.         end_if
    
41.         loop local m (1,9)
    
42.             command
    
43.                 ball create group arrange_two ...
    
44.                      radius @ini_rad x @h_x y @h_y z @h_z
    
45.             end_command
    
46.             h_x += 2.0*ini_rad
    
47.         end_loop
    
48.     end_loop
    
49. end
    
50. @hexagon_arrange
    
51. 
    
52. ball generate group arrange_three ...
    
53.      box -21 -5 -0.5 0.5 -21 -5 ...
    
54.      radius 1.0 cubic
    
55. ball generate group arrange_four ...
    
56.      box 5 22 -0.5 0.5 -19 -5 ...
    
57.      radius 1.0 hexagon
    

复制代码

规则排列效果图

公子无心

谢谢分享！！！！

rmy

赞一个，楼主威武