本文主要是介绍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. 保存图片的路径
绘制结果如下：
二维模量投影图与三维一样：
说明：将二维投影图的数据导出，在其他画图软件绘制效果会好一点！
欢迎大家交流，我们会持续更新，但是微信上面更新的快，关注微信公众号：原子与分子模拟
这篇关于GULP案例3代码解释：计算弹性常数和绘制弹性模量二维、三维图的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
