24小时热门版块排行榜    

查看: 2643  |  回复: 4

明天琳琳

新虫 (初入文坛)

[求助] 纳米压痕相关的lammps程序(新手上路) 已有2人参与

我刚刚学lammps没多久,自己刚刚弄好了一个程序,怎么改都错了好多,特请教大神帮忙提点错误……

units             metal
dimension         3
boundary          p p m

atom_style        atomic
neighbor          2.0   bin
neigh_modify      delay 3
lattice           fcc   3.615

#define groups
region            box1 block 0 20 0 20 0  10 units box
region            box2 block 0 20 0 20 10 20 units box
region            box3 block 0 20 0 20 25 31 units box
create_box        1 box1
create_box        1 box2
create_box        1 box3
create_atoms      1 region box1
create_atoms      1 region box2
create_atoms      1 region box3
group             11 region box1
group             22 region box2
group             33 region box3

#Potential Setting
pair_style        eam
pair_coeff        * * Cu_u3.eam Cu
mass              1 63.55

#indention
fix               1 11 setforce 0.0 0.0 0.0
variable          z equal vdisplace(28,-0.1)
fix               2 33 indent 100 sphere 10 10 v_z 3.0

compute           cna all cna/atom 4.363
compute           str all stress/atom
fix               3   all nvt temp 0 0 0.1

variable          Temp equal temp
variable          EnerPot equal 20

thermo            ${ThermoDump}
thermo_style      custom step temp pe
dump              5 all cfg 100 dump.*.cfg mass type xs ys zs id c_cna c_str[2]

fix               6 all print 500 "${Step} ${Temp} ${EnerPot} ${pe}" file Cu.dat title "Temp pe" screen no
minimize          1e-10 1e-12 2000 2000

timestep          0.001
run               1000000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明天琳琳

新虫 (初入文坛)

对了,忘记说,我的程序的目的是为了看如何去金属材料的力学性质
2楼2017-01-08 21:07:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MD小师弟

新虫 (初入文坛)

【答案】应助回帖

create_box命令用多了,应该是只有一个box
nvt系综的温度设置为0.01
在box1到3 的region里建立原子,如果分别在三个区域中建立原子,可能会产生原子距离过近的问题,所以可以用 region union的命令把三个region整合为一个再建立原子(region之间有交集关系并无影响)

暂时就想到这么多,不一定完全正确,仅做参考
fix indent这玩意我也不太懂,如果有什么进展,还请不吝赐教!
3楼2017-01-13 12:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明天琳琳

新虫 (初入文坛)

对不起,当时有段时间没回复,我就没有再看小木虫,结果金币没法给了,抱歉……并且感谢你的回复……
4楼2018-03-28 09:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

回头看我

新虫 (小有名气)

【答案】应助回帖

楼主,你这样的压头是单向运动的,怎么返回压头得到卸载曲线呢?
5楼2020-06-13 21:04:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 明天琳琳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
信息提示
请填处理意见