计算化学公社

标题: 关于Lammps拉伸过程中应力的换算方法问题。 [打印本页]

作者
Author: sasuke1000 时间: 2023-4-11 23:56
标题: 关于Lammps拉伸过程中应力的换算方法问题。
本帖最后由 sasuke1000 于 2023-4-12 00:00 编辑

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

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

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

作者
Author: sasuke1000 时间: 2023-4-18 12:56

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

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

作者
Author: uu4 时间: 2024-9-21 21:36
本帖最后由 uu4 于 2024-9-21 21:38 编辑

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

作者
Author: jackieiii 时间: 2025-1-24 15:15

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)