主程序p2voxel.py，读取PYline3Data.mat内含有各类微结构拓扑信息，由matlab生成。之后会整理matlab方面代码。

test.py为读入模型，使用p2voxel生成体素，输出DH的代码示例。

## 参数说明

```python
in_parameter_names, out_parameter_names = showPara(mtype)
```

得到类似输出

```
in_parameters ['R1', 'R2', 'R3', 'F12_1', 'F12_2', 'F13_1', 'F13_2']
out_parameters ['E5', 'E6', 'E9']
```

in_ 表示结构内部参数

out_ 表示接触面上参数。不同微结构类型，如果有着同样的接触面参数与取值，则可以相连。

```python
init_parameters = {
    'R1':0.036, 'R2':0.036, 'R3':0.036,
    'E5':0.1,'E6':0.1,'E7':0.1,'E8':0.1,'E9':0.1,'E10':0.1,
    'F11_1':0.2,'F12_1':0.2,'F13_1':0.2,'F14_1':0.2,
    'F11_2':0.2,'F12_2':0.2,'F13_2':0.2,'F14_2':0.2,
    'T1':0.1,'T2':0.1,'T3':0.1
}
```

生成的voxel单元，视为单位长度1\*1\*1的体素。

"R1", "R2", "R3"为对应杆单位长度半径

"E{num}"为对应六条边上点的位移比例，范围0<E{num}<1

"F{num}"为对应四个面上点的位移比例，范围0<F{num}\_1+F{num}\_2 < 1

"T1", "T2", "T3"为体内点位移比例，范围0<T1+T2+T3<1

且所有参数均需要满足 大于0 小于1

## 其他说明

目前主要处理的体素分辨率为40\*40\*40

```
voxel = p2voxel(mtype, parameters, resolution=40)
```