24小时热门版块排行榜

返回列表

honghaier90

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 208.6
散金: 10
红花: 1
帖子: 52
在线: 46小时
虫号: 2581277
注册: 2013-08-05
性别: MM
专业: 高分子合成化学

[求助] lammps里，关于应力应变的问题

我在用lammps做材料的单轴拉伸和压时，遇到以下几个问题：
1，第0步输出的应力不为零，第0步材料未发生变形，为什么会出现应力，按理说应力也应该为零啊。拉压之前，均做过能量最小化，但是还是会出现这种情况，这个该如何处理？
2，对于同种材料分别做拉伸压缩，在得出应力应变曲线后，通过直线拟合，发现两条直线斜率不相等，也就是弹性模量不相等，这是为什么？
以下是我做拉伸时的in文件，恳请大家给于指点。
dimension       2
boundary       s s p

atom_style       angle
neighbor       0.3 bin
neigh_modify       every 10 delay 0 check yes
lattice       hex 0.91648

read_data          data.new1

pair_style       lj/cut  1.5
pair_coeff       1 1 1.0 1.0  1.5

bond_style    harmonic
bond_coeff    1 100 1.12246

angle_style harmonic
angle_coeff 1 1 60

region             left block INF 10    INF  INF  INF  INF units lattice
region             right  block 190 INF INF  INF  INF  INF units lattice
group             left region  left
group             right  region  right
group             boundary  union  left right
group             mobile  subtract  all boundary

thermo_style  custom lx

#min_style    sd
#minimize       1.0e-12 1.0e-12 10000 10000

variable       tmp  equal  lx
variable       L equal  ${tmp}
print          "initial  length , L:${L}"
#deformation
reset_timestep  0
velocity       left set 0.0 NULL 0.0
compute          new mobile temp
variable       pressx equal c_p[1]
compute       s all stress/atom
compute       strex all reduce sum c_s[1]
variable       x equal c_strex/vol
variable       strain1 equal "(lx - 224.6726)/224.6726"
#variable       strain2 equal "(-ly + 69.107548)/69.107548"
thermo_style custom step temp vol lx ly lz pxx v_x c_strex v_strain1
thermo_modify    temp new

velocity mobile create 0.01 887723 temp new
velocity right  set 0.3 0.0 0.0
velocity mobile ramp vx 0.0 0.3 x 10 190 sum yes



fix 1 all nve
fix 2  boundary setforce 0.0 NULL 0.0
fix 3 mobile temp/rescale 100  0  0.05  0.002  0.1

timestep 0.001
thermo 100

dump 1 all atom 100 dump.pull1.lammpstrj
run 80000

回复此楼