24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2460  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 明天琳琳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 生物学调剂 +10 纸扇zhishan 2026-04-13 10/500 2026-04-18 21:19 by zqndavala
[考研] 化学070300 求调剂 +29 哈哈哈^_^ 2026-04-12 29/1450 2026-04-18 15:56 by Equinoxhua
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 化工学硕294分,求导师收留 +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 5/250 2026-04-17 10:02 by bobvan
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
信息提示
请填处理意见