lammps数据后处理：python绘制应力应变曲线 附程序代码

大家好，我是小马老师。

从本文开始，公众号推出一个新的专题：lammps数据后处理，主要介绍使用python编程处理lammps模拟的相关数据。

本文先介绍一个python绘制应力应变曲线的脚本程序。

绘制应力应变的方法有很多，常规的做法是把数据文件拖入到origin绘图，还有一个简单的方法是使用python脚本，在模拟完成后，直接运行一下脚本就能得到应力应变曲线，可以快速的观察运行结果。

基本的流程：

（1）读取数据

应力应变数据使用fix print命令保存到strain_stress.txt中
fix def1 all print 100 “${p1} ${p2}” file strain_stress.txt screen no
使用open命令读取文件

filename = "strain_stress.txt"
f = open(filename,'r')

（2）分析数据

st