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呢？

想请问下楼主，我做高熵合金的模拟，做出来也是应力很大，有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

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

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

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

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