我算的一个纳米线，结果我没算完就让他停掉了，因为最大应力比人家文献里的要小，如果取sss的边界条件，再用velocity 或displace_atom。算出来的最大应力跟他文献里差不多，但初始应力（应变为0时的应力）有1点多GPA，因为非周期边界没法用fix nPt 降压，但纳米线真正的边界条件应该取ssp亚，在P方向拉伸，我是这样做的，再用fix deform ，但算出来的应力又太小了，我觉得可能是这个命令我没用对，请大家帮我看看代码有没有错误，特别是fix deform命令，谢谢了。我是参考附件也附上

boundary    p s s
units        metal
atom_style    atomic
neighbor    2.0 bin
neigh_modify    delay 1 check yes

# create geometry

lattice        fcc 3.61
region        box block 0 30 0 3 0 3
create_box    1 box
create_atoms    1 box

mass        1 63.546

# LJ potentials

pair_style    eam
pair_coeff    * * Cu_u3.eam

# define groups

region            1 block INF 1 INF INF INF INF  
group        left region 1
region        2 block  29 INF INF INF INF INF
group        right region 2
group        boundary union left right
group        mobile subtract all left

# initial velocities
velocity    left set 0.0 0.0 0.0

compute p all pressure thermo_temp
variable     pressx equal c_p[1]
variable     pressy equal c_p[2]
variable     pressz equal c_p[3]


thermo_style    custom step temp etotal press v_pressx  v_pressy v_pressz vol
dump    1 all atom 1000 dump.cupoly
minimize 1.0e-6 1.0e-6 1000 1000

fix        1 left setforce 0.0 NULL NULL
fix        2 all npt 1.0 1.0 1.0 aniso 0.0 0.0 NULL NULL NULL NULL 10.0 drag 1.0

thermo 1000
timestep    0.002
run        10000


unfix 2
fix  3  all nve
fix  4  all temp/rescale 100 1.0 1.0 5 1.0
fix  5  all deform 1000 x erate 0.002 units box

compute         s all stress/atom
compute         strea all reduce sum c_s[1]
variable        stressa equal c_strea/vol
compute         streb all reduce sum c_s[2]
variable        stressb equal c_streb/vol
compute         strec all reduce sum c_s[3]
variable        stressc equal c_strec/vol

thermo_style    custom step temp etotal press v_stressa  v_stressb v_stressc vol
run        300000

[ 本帖最后由 redream 于 2008-9-5 16:35 编辑 ]

我说两句吧。
1。fix 1感觉没必要，可能你是想让left边界处原子固定不动，但这样起不到效果，fix deform会让指定group里的原子都参与deform，同时，由于你计算应力时考虑的是group all，而group left中的原子受力被人为改变，这将会对计算应力造成一定影响
2。thermo 1000 热力学结果输出间隔感觉过长，不知道你得到的最后结果怎样，可否将应力-时间步和应力-应变曲线贴上来看看
3。用ssp边界条件是合理的，你说在sss边界条件下初始应力高达1.0GPa，可能是因为你没有让体系平衡，正常情况下sss边界条件对应于零应力状态
4。最大应力与各方面因素有关，体系温度，应变速率，模型尺寸，加载方式等，不知道你是怎么进行对比的？