24小时热门版块排行榜    

查看: 1340  |  回复: 1

Berryfoss

铁虫 (初入文坛)

[交流] 【求助】Lammps 中加外力的问题 已有1人参与

在非周期性边界方向,想加一个确定大小的外力,用什么命令呢“?用fix setforce时,run几步后总是出现atom lost的错误,而且怎么确定setforce的值也是个问题,如果用fix deform,是不是不能用在非周期性条件下阿?这种情况下应该用什么命令呢?希望大家给点建议。

in文件如下

units                        metal
boundary                p p s
atom_style                atomic
neighbor        2 bin
neigh_modify        every 5 delay 0 check yes

#create geometry
lattice         fcc 3.61
region L block 0 40 0 20 0 20 units lattice
lattice  fcc 3.52
region R block 42 80 0 20 0 20 units lattice
region W block 0 80 0 20 0 20 units lattice
create_box 2 W
create_atoms 1 region L
create_atoms 2 region R
mass 1 63.546
mass 2 58.7

# potentials
pair_style        eam/alloy
pair_coeff        * * /home/foss/Applications/lammps-16Sep10/examples/CuNi/CuNi_Zhou.eam.alloy Cu Ni

#define groups
group a region L
group b region R
group c region W
region  3  block 77 80 0 20 0 20
group right region 3

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

velocity c create 1000.0 4928459 rot yes dist gaussian

#initialize
fix 1 c nvt temp 1000.0 1000.0 0.05
fix 2 a addforce 3.1e5 0 0 region L
fix 3 right setforce 0.0 0.0 0.0

#dump
dump 1 c xyz 10 dump.nicu
timestep 0.002
thermo 10
run 10000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lootoo

禁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

2楼2020-08-20 21:36:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Berryfoss 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见