GULP案例3代码解释：计算弹性常数和绘制弹性模量二维、三维图

GULP计算弹性常数和弹性模量的输入文件如下：

说明：输入文件与前面相比，唯一变动的地方是strain命令。

opti prop conp strain
name ZnO-WZTemperature 0.0 K
Pressure 15.0 GPacell
3.242  3.242  5.188  90.0  90.0  120.0
frac
Zn core 0.333333   0.666667   0.000000  
O  core 0.333333   0.666667   0.381900
space
186species
Zn   1.14
O   -1.14buck 
Zn  Zn  78.91   0.5177  0.00     0.0 12.0
Zn  O   257600  0.1396  0.00     0.0 12.0Morse
O   O   0.1567  1.164   3.405    0.0 12.0 
dump every 1 ZnO-WZ.dump

弹性各向异性

ElasticPOST 弹性常数各向异性后处理软件（带图形界面）

主要功能

本软件(elasticpost)由matlab编写而成，主要是对计算得到的弹性刚度矩阵进行后处理，主要功能如下：

1. 对任意对称性的三维材料可以体弹模量、杨氏模量、剪切模量、泊松比等与弹性刚度矩阵相关量进行三维各向异性绘图；

2. 对任意对称性的二维材料可以对杨氏模量、剪切模量、泊松比进行弹性常数二维各向异性图的绘制；

3. 可以按照一定格式读取 excel、 txt 以及 mat 文件，以达到批量处理的目的；

4. 可以根据 voigt-reuss-hill 近似对上述量进行计算， 并可以将数据进行保存并与 matlab 工作空间交互；

5. 对三维图形提供多种图形属性的调控操作，如从最高位置看，切片等；

6. 对图片的输出提供.tif、 .jpg、 .emf 等多种位图格式以及.eps、.pdf、.fig 等多种矢量格式。

有编译后的版本(见bin目录)，该版本不需要matlab，但是需要matlab runtime（2019b version），安装会自动下载（大约2.5 G）

下载地址

elasticpost github：https://github.com/hitliaomq/elasticpost#opennewwindow

参考链接：http://muchong.com/t-14302716-1-pid-9

（转载）参考链接：http://muchong.com/t-14302716-1-pid-9

此链接中有详细的安装过程。

操作过程：

安装好之后，打开界面如下：

以ZnO-WZ结构为例，进行弹性模量的绘制

第一步：选择晶系（WZ结构属于六方晶系）

第二步：输入完ZnO-WZ结构所有的弹性常数：

1. 选择二维或者三维模量投影图

2. 选择弹性模量：图中为体弹模量B，此外还可选择杨氏模量E、剪切模量G等

3. 输入文件名称：图中为Diamond

4. 选择保存的数据类型

5. 保存

体弹模量B绘制结果如下

第三步：右边菜单栏介绍：

1. 选择晶向

2. 选择切片方向

3. 选择添加坐标轴

4. 保存图片的路径

绘制结果如下：

二维模量投影图与三维一样：

说明：将二维投影图的数据导出，在其他画图软件绘制效果会好一点！

欢迎大家交流，我们会持续更新，但是微信上面更新的快，关注微信公众号：原子与分子模拟