24小时热门版块排行榜

>论坛更新日志 (2822)
>虫友互识 (321)
>考博 (157)
>导师招生 (149)
>文献求助 (122)
>基金申请 (94)
>休闲灌水 (91)
>仿真模拟 (45)
>文学芳草园 (26)
>论文投稿 (21)
>硕博家园 (18)
>催化 (17)
>博后之家 (16)
>公派出国 (16)
>考研 (16)
>找工作 (14)

返回列表

cmhrbeu

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 123.5
帖子: 9
在线: 4.1小时
虫号: 2753315
注册: 2013-10-25
专业: 核技术及其应用

[求助] lammps弛豫阶段一直无法在设定的温度上下波动，材料的辐照损伤

# 3d UO2 crash simulation

units metal
dimension 3
boundary p p p
atom_style charge
read_data uo2.lammps

# hybrid/overlay born/coul/long potentials

group u type 1
group o type 2
set group o charge -1.1104
set group u charge 2.2208
pair_style hybrid/overlay born/coul/long 10.4 morse 10.4
pair_coeff 1 1 born/coul/long 187.03 0.3422 0 0 0
pair_coeff 1 2 born/coul/long 432.18 0.3422 0 0 0
pair_coeff 2 2 born/coul/long 883.12 0.3422 0 3.996 0

pair_coeff 1 2 morse 0.5055 1.864 2.378

kspace_style pppm 1e-4

# fixes

min_style sd
minimize 1.0e-6 1.0e-6 10000 10000
velocity all create 500.0 50000 rot yes dist gaussian
fix 1 all nvt temp 300.0 300.0 3.0 drag 1

# run

timestep 0.001
thermo 100
thermo_style custom step temp pe ke

#neigh_modify every 1 delay 0 check yes page 10000

compute 1 all pe/atom
compute 2 all ke/atom
dump 1 all custom 1000 cm.crash id type x y z c_1 c_2

run 50000

unfix 1
write_restart cm.rs

以上是in文件，设定的温度是300K
小弟初学lammps，毕业设计进行中，望大神给力。

lammps.jpg

回复此楼