关于Lammps拉伸过程中应力的换算方法问题。

sasuke1000

# Initialization
units       real
boundary    p p p
atom_style  full  #molecular
read_data   2304265K.data
timestep    0.0005
# Dreiding potential information
neighbor    2 bin
neigh_modify every 1 delay 0 check yes
neigh_modify one 5000
#-------------------------应力应变保存到文件---------------------------
fix                 def3 all print 100 "${strainx} ${stressx} ${strainy} ${stressy} ${strainz} ${stressz}" &
                    file stress-strain.dat screen no
#-------------------三向拉伸--------
fix                 1 all nve
fix                 2 all temp/rescale 10 265 265 5 1.0
fix                 3 all deform 100 y erate 0.005  remap x units box
dump                1 all custom 100 water_tension.xyz id type x y z
run                 1200000

做一个材料的拉伸，部分拉伸代码如上，最后输出的stress y达到11GPa的数量级，明显不是这个材料的屈服强度。

我查阅了real unit的单位制，pressure的单位是atm，也就是1个大气压，这个11是已经除以了10000的，所以数值计算没有问题。
我又看了L-J势下关于pressure的定义， ,感觉并不能据此算出一个与系统结果不同的、有意义的p。
吊诡的是，有的国内的EI期刊做水泥混凝土的分子模拟的，拉伸强度达到GPa数量级也能发表，是我混凝土原理学错了吗，钢筋的抗拉强度也才几百兆帕，这混凝土早就拉开裂了。所以在这里真诚地询问大家，这些究竟是咋回事，我物理的底子薄，是工科生，可能很多统计力学、分子模拟原理不是很精通，想问下是否存在一个转换方式，把lammps的这个结果转换成接近真实的数字？抑或这是这就是我的模拟的真实结果，只是我的模拟的参数设置有误呢？PS：这个系统下材料的密度还是比较正常的。

sasuke1000

我更新了拉伸时的系综为NPT，然后根据材料的极限拉应变的数值调整了erate、timestep、run步数的关系，目前运行下来貌似正常了，没那么大了。

清风与鹿

请问real里force与Gpa的换算楼主知道吗？
我看metal里force单位可以直接和Gpa换算，想请教一下real里的换算

sasuke1000

清风与鹿 发表于 2023-4-14 15:34
请问real里force与Gpa的换算楼主知道吗？
我看metal里force单位可以直接和Gpa换算，想请教一下real里的换 ...

在手册里有，你用VSCODE写代码的时候放到real上面会自动显示手册里这块内容，lammps里面real单元的压力（press）的单位是atm，你输出的时候调PXX/10000，出来的单位大约就是GPa，得换算一下

uu4

想请问下楼主，我做高熵合金的模拟，做出来也是应力很大，有22 GPa这么大，实验测出来只有400 MPa左右。我尝试调整了模型大小、系综、erate、timestep、run这些参数，结果几乎没什么变化。我感到很迷惑，请问楼主您当时怎么调整计算的？

==============================================以下是我的in文件==========================
#--------------------初始化--------------------------
units                         metal
dimension                    3
boundary                  p p p
atom_style               atomic

#-------------------2--------------------------------
lattice                      fcc 3.592   #Cu晶格 a=b=c=3.592
region                       box1 block 0 20 0 15 0 70

create_box               4    box1  #4 种原子
create_atoms            1 region box1

#------------原子类型转化---------------------------
set type 1 type/ratio 2 0.25 87393
set type 1 type/ratio 3 0.33333 87393
set type 1 type/ratio 4 0.5 87393

#-----------原子类型信息---------------------------
mass 1 55.845           #Fe
mass 2 51.996           #Cr
mass 3 58.693           #Ni
mass 4 58.933194     #Co




#----------------势函数设置-------------------------
pair_style eam/alloy
pair_coeff   * * FeNiCrCoCu-with-ZBL.eam.alloy Fe Cr Ni Co



timestep 0.01
velocity all create 300 12345 dist gaussian mom yes rot yes

#----------------------设置输出格式-----------------
thermo 100
thermo_style custom step lx ly lz press pxx pyy pzz ke pe temp

min_style cg
minimize 1e-10 1e-12 10000 10000
variable tmp equal "lz"
variable L0 equal ${tmp}
variable strain equal "(lz-v_L0)/v_L0"
variable stressz equal "-pzz/10000"
variable srate equal 1.0e10
variable srate1 equal "v_srate / 1.0e12"
fix def3 all print 100 "${strain} ${stressz}" file stress-strain.dat screen no
fix     4  all print 100 "${stressz}" file stress.txt screen no
fix     5  all print 100 "${strain}" file strain.txt screen no

fix 1 all deform 10 z erate 0.005 remap x units box

fix 6 all npt  temp 300 300 0.1  x 0 0 1000 y 0 0 1000


dump 1 all custom 50 dump.min id type x y z ix iy iz vx vy vz
dump_modify                 1 sort id

thermo_style custom step time temp pe etotal press vol

run 14000

jackieiii

uu4 发表于 2024-9-21 21:36
想请问下楼主，我做高熵合金的模拟，做出来也是应力很大，有22 GPa这么大，实验测出来只有400 MPa左右。我 ...

你现在解决了吗，我也遇到相同的问题，但我看了一下你的in文件觉得有两个疑问：
fix 1 all deform 10 z erate 0.005 remap x units box
z轴拉伸为啥要remap x？
fix 6 all npt  temp 300 300 0.1  x 0 0 1000 y 0 0 1000
x, y轴方向的压力为何设置成0而不是1bar呢？