8.8. 后处理

8.8.1. 数据

用户可以使用PyRunner引擎来保存 SudoDEM 运行中的数据:

1def savedata():
2        global saving_num
3        saving_num += 1
4        O.save('shear'+str(saving_num)+'.xml.bz2')
5# 假如已经载入了固结数据，下述代码可以省去，因为调用 'savedata()' 的PyRunner已定义
6O.engines = O.engines + [PyRunner( command = 'savedata()', iterPeriod = 1000 )]

Python 包 utilspost 提供了一系列的后处理函数用来协助用户后处理，包括应力-力-组构计算以及用于调用其他后处理软件的接口函数 (例如，颗粒形状的vtk文件保存函数, 颗粒信息3D直方图的vtk文件保存函数)。基于如上包用户可以非常容易调用或自定义适合自己的后处理函数。

1import sudodem.utilspost as utilspost
2#如下函数可以帮助用户生成接触力发挥系数计算文件
3for i in ['1.0','1.2','1.4','1.6']:
4    utilspost.calc_avgFricMobilIndex(i,steps_num)

Python 包 _superquadrics_utils 提供如下后处理辅助函数,

  • outputParticles(filename) : 该函数接收一个文件名，并将保存每个颗粒的形状信息到该文件中，包括形状参数 \(r_x\), \(r_y\), \(r_z\), \(\epsilon_1\), \(\epsilon_2\), 位置 \(x, y, z\), 和朝向 \(q0, q1, q2, q3\) (等价于 \(q_w, q_x, q_y, q_z\))。

  • outputParticlesbyIds(filename, ids) : 该函数接收一个文件名以及颗粒索引，该函数将具有该索引的颗粒的信息保存到该文件中，保存信息与 outputParticles(filename) 相同。

  • outputWalls(filename) : 该函数接收一个文件名, 并将墙体的位置信息保存到该文件中。

  • outputPOV(filename) : 该函数接收一个文件名，并将颗粒的信息以POV-Ray文件格式保存，该POV-Ray文件可被用于POV-Ray进行颗粒的渲染。

  • outputVTK(filename,*slices*) : 该函数接收一个文件名，并将颗粒的信息以VTK文件格式保存，用户可以使用Paraview等VTK软件对颗粒进行渲染。

8.8.2. 颗粒模拟可视化

8.8.2.1. SudoDEM3D

我们强烈推荐两款开源且免费的可视化软件用于 SudoDEM 的后处理可视化 Paraview and POV-Ray ：

../_images/paraview.png

图 8.21 Paraview 的主程序窗口

高精度的颗粒信息可视化相较于简单的截图可以更清晰表征模拟结果。这里我们将为大家演示如何利用vtk或者pov文件进行高精度颗粒信息可视化， 我们假设计算结果已经被保存至文件 "shearend.xml.bz2" 。

1O.load("shearend.xml.bz2")
2import sudodem._superquadrics_utils as superutils
3superutils.outputVTK("shearend.vtk",15)

打开 Paraview, 导入vtk文件 "shearend.vtk" 然后点击按钮 "Apply" 从而得到如 图 8.22 所示的可视化结果。 需要说明的是函数 outputVTK 将颗粒的表面视为多面体的切片组合(切片的精度这里被设置为 15) , 依此法可能会得到较大的后处理文件 (该例子文件大约180 Mb)。基于此用户可以选择多种多样的方式来减少后处理文件大小，比如用户可以通过 Paraview 提供的快速可视化方法，或者可以利用 Paraview 内置接口 Superquadric 来降低可视化文件大小。

../_images/shearendvtk.png

图 8.22 Paraview 渲染的颗粒

用户还可以通过使用*POV-Ray* 来降低可视化文件大小。

1O.load("shearend.xml.bz2")
2import sudodem._superquadrics_utils as superutils
3superutils.outputPOV("shearend.pov")

用户可以通过修改POV文件的相机设置来达到自定义的可视化效果, i.e., "shearend.pov" here,

1camera{ //orthographic angle 45
2        location <0.6,0.6,0.5>*50e-3
3        sky z
4        look_at  <0,0,0.05>*10e-3
5        right x*image_width/image_height
6        translate <-0.04,0.01,-0.1>*50e-3
7    }

这里给出了如何自定义光源的设置。

1// White background
2background{rgb 1}
3// Two lights with slightly different colors
4light_source{<4,8,5>*10e-3 color rgb <1,1,1>}
5light_source{<12,-6>*10e-3 color rgb <1,1,1>}

在命令行中运行如下命令来可视化颗粒:

1povray shearend.pov

如果要修改可视化精度，用户可以添加额外的命令参数来达到这一目的,

1povray shearend.pov +A0.01 +W1600 +H1200

图像处理软件, 例如, GIMP , 可用来修改图片大小或者制作动图。图 8.23 给出了使用 POV-Ray 渲染并利用 GIMP 裁剪的颗粒渲染图片。

../_images/shearend.png

图 8.23 POV-Ray 渲染的颗粒

包 snapshot 提供了扩展超椭球的可视化函数。例如，我们导入堆积结果完成时的颗粒状态文件并利用 snapshot 进行POV文件生成：

1from sudodem import snapshot
2snapshot.outPov("packingpolysuper",withbox=False)

此处我们得到两个文件分别是 packingpolysuper.inc 和 packingpolysuper.pov。 第一个文件 packingpolysuper.inc 包含了一些扩展超椭球的宏定义以及颗粒形状参数信息，例如颗粒形状参数，颗粒位置，颗粒朝向等

 1#macro myremain(sign1,sign2,sign3)
 2    clipped_by{plane{-sign1*x,0}}
 3    clipped_by{plane{-sign2*y,0}}
 4    clipped_by{plane{-sign3*z,0}}
 5#end
 6#macro randomColor(r,g,b)
 7texture{ pigment{ color rgb<r,g,b> transmit para_trans} finish { phong para_phong }}
 8#end
 9#macro octantSuper(ep1,ep2,rx,ry,rz,sign1,sign2,sign3)
10    superellipsoid{ <ep1,ep2>
11#if (singleColor = false)
12randomColor(rand(Random_r),rand(Random_g),rand(Random_b))
13#end
14myremain(sign1,sign2,sign3) scale <rx,ry,rz>}
15#end
16#macro polySuperEllipsoid(ep1,ep2,a1,a2,b1,b2,c1,c2,t1,t2,t3,r1,r2,r3,r,g,b)
17union{
18octantSuper(ep1,ep2,a2,b2,c2,-1,-1,-1)
19octantSuper(ep1,ep2,a1,b2,c2,1,-1,-1)
20octantSuper(ep1,ep2,a2,b1,c2,-1,1,-1)
21octantSuper(ep1,ep2,a1,b1,c2,1,1,-1)
22octantSuper(ep1,ep2,a2,b2,c1,-1,-1,1)
23octantSuper(ep1,ep2,a1,b2,c1,1,-1,1)
24octantSuper(ep1,ep2,a2,b1,c1,-1,1,1)
25octantSuper(ep1,ep2,a1,b1,c1,1,1,1)
26rotate<r1,r2,r3>
27translate<t1,t2,t3>
28#if (singleColor = true)
29randomColor(r,g,b)
30#end
31}
32#end
33polySuperEllipsoid(1.0000000e+00,1.0000000e+00,1.0000000e-01,
34                   5.0000000e-02,8.0000000e-02,1.5000000e-01,
35                   2.0000000e-02, 4.0000000e-02,-3.2218276e-01,
36                   1.8390308e-01,3.1674861e-01, -167.727379,
37                   16.417442,-104.719219,0.783099,0.394383,0.840188)
38...
39...

我们不建议用户直接修改第一个文件，但是用户可以通过修改第二个文件来达到自定义渲染效果:

 1//using the command: povray **.pov +A0.01 +W1600 +H1200
 2#include "colors.inc"
 3#declare singleColor =  true;
 4#declare Random_r = seed (1432);
 5#declare Random_g = seed (7242);
 6#declare Random_b = seed (9912);
 7#declare para_trans =  0 ;
 8#declare para_phong =  1.0 ;
 9camera {
10location < 1.39244459003 , 0.950562440499 , 0.605226784362 >
11sky z
12right -x*image_width/image_height
13look_at < 0.612872382451 , 0.349255160637 , 0.430021966053 >
14}
15#include "packingellipses.inc"
16light_source{<4,8,5>*10 color rgb <1,1,1>}
17light_source{<12,-6>*10 color rgb <1,1,1>}
18light_source { <0, 2, 10> White }
19background{rgb 1}
20//you may comment out this line to remove the plane
21plane { z, -5 pigment { checker Green White }}

命令行输入以下命令来完成渲染

1povray packingpolysuper.pov +A0.01 +W1600 +H1200

渲染完毕后用户将会得到如 图 8.24 所示的图片。

../_images/packingpolysuperpov.png

图 8.24 POV-Ray 渲染 扩展超椭球的堆积模拟 中的扩展超椭球堆积状态

8.8.2.2. SudoDEM2D

包 _superellipse_utils 提供了名为 drawSVG 的函数用于绘制颗粒力链的矢量图。

1from sudodem import _superellipse_utils
2O.load("your data file")
3drawSVG("test.svg", draw_forcechain=True,
4                force_line_color=(1.0,0,0),
5                solo_color=(93.0/255,152.0/255,208.0/255))

生成的矢量图文件可以通过文本编辑器或者矢量图编辑器 Inkscape (推荐)来修改。图 8.25 给出了椭球颗粒在周期边界条件下堆积过程结束后的力链矢量图示例，用户可通过 drawSVG 函数来配置对应的输出图片格式。

../_images/bed90.png

图 8.25 drawSVG 导出的颗粒和力链