在区域内创建大小不同的小球，使大球固定不动，小球施加重力下落
模型创建new ; 释放当前内存，开始新的任务分析 title 'test 1' ; 设置模型名称用于图像显示 ; 设置计算区域，必须在ball、wall等实体部分建立前设置 domain extent -100.0 100.0 condition periodic ; 设置周期性边界，当颗粒、wall、clump等碰到domain边界时，处理方式为destroy （删除），stop （停止运动），reflect （速度反向，弹回），periodic（从domain）相对面重新出现，常用于均匀化方法。 set random 10000 ;设置随机种子 ; 设置模型的边界wall wall generate box -25.0 25.0 -50.0 50.0 ; 生成一个矩形wall，box后的参数为xmin, xmax, ymin, ymax ; 创建颗粒体系（ball、clump、cluster等）并分组 ball generate radius 1.0 1.5 box -25.0 25.0 -50.0 50.0 number 450 ; 在box内随机生成半径为1.0~1.5的450个球 ball group small_balls range radius 1.0 1.3 ; 将半径处于1.0~1.3的球编为一组，名为small_balls ball group big_balls range radius 1.3 1.5 ; 将半径处于1.3~1.5的球编为一组，名为big_balls
运行结果：
显示已经创建完毕，但是右方绘图区尚未激活，所以看不到创建的实体。
激活方式，首先点击一下绘图区，再右上角处添加绘图选项。
添加之后，图形显示正常：
; 设置球的实体属性（密度、速度、阻尼等） ball attribute density 100.0 ; 设置密度 ball fix velocity range group big_balls ; 将big_balls组球的速度设为固定 ball attribute radius multiply 1.2 ; 半径扩大 ball attribute damp 0.7 ; 设置阻尼 cmat default model linear property kn 1.0e8 fric 1.0 ; 设置接触（contact、cmat等） ball property kn 2e8 ks 1e8 fric 1.0 ; 属性继承 set gravity 10.0 ; 添加外力求解计算
set timestep auto ; 设定时间步长 set timestep max 5e-3 wall history id 1 contactforce id 1 ; 创建数据记录 ; 计算求解 step 1000 cycle 2002 solve time 10.0; solve fishhalt @stop_me ; 输出数据 history write 1 file testresult
求解完成后的图形：
输出结果会自动放到文件所在位置：
