帮我生成LAMMPS计算hydro_stress和von Mises stress 的脚本
时间: 2023-10-14 17:04:14 浏览: 129
以下是一个简单的 LAMMPS 脚本,用于计算 hydrostatic stress 和 von Mises stress:
```
# 初始化
units metal
dimension 3
boundary p p p
atom_style atomic
# 定义晶胞
lattice fcc 4.05
region box block 0 10 0 10 0 10
create_box 1 box
create_atoms 1 box
# 定义势能函数
pair_style eam/alloy
pair_coeff * * Al99.eam.alloy Al
# 定义温度和压力
velocity all create 300 87287
fix 1 all npt temp 300 300 0.1 iso 0.0 0.0 1.0
# 运行 MD
timestep 0.001
thermo 100
run 5000
# 计算应力
compute stress all stress/atom NULL
compute hs all reduce sum c_stress[1] c_stress[2] c_stress[3]
compute vs all reduce sum (c_stress[1]-c_stress[2])^2+(c_stress[2]-c_stress[3])^2+(c_stress[3]-c_stress[1])^2
# 输出结果
variable hs equal -(c_hs[1]+c_hs[2]+c_hs[3])/3.0
variable vs equal sqrt(3.0*c_vs/vol)
thermo_style custom step temp v_hs v_vs
thermo_modify flush yes
run 10000
```
这个脚本使用 EAM 势函数模拟铝的行为,并计算了每个原子的应力。然后,通过两个计算来得到 hydrostatic stress 和 von Mises stress。最后,将结果输出到屏幕上。
相关推荐
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)