24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2444  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 明天琳琳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +15 ORCHID1 2026-04-10 16/800 2026-04-12 00:45 by 蓝云思雨
[考研] 0831一轮调剂失败求助 +6 小熊睿睿_s 2026-04-11 6/300 2026-04-11 23:55 by baobaoye
[考研] 267求调剂 +8 再忙也要吃饭啊 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 一志愿085502,267分求调剂 +11 再忙也要吃饭啊 2026-04-08 12/600 2026-04-11 11:08 by zhq0425
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 080500求调剂 +17 黄宇博 2026-04-06 17/850 2026-04-11 08:36 by zhq0425
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 生物与医药调剂 +5 十七sa 2026-04-05 5/250 2026-04-10 08:14 by kangsm
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 材料科学与工程320求调剂,080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 材料工程专业日语生求调剂 +9 111623 2026-04-07 9/450 2026-04-07 23:31 by 一只好果子?
[考研] 材料求调剂 +18 一样YWY 2026-04-05 18/900 2026-04-07 15:49 by dxlg
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 第一志愿东南大学物理313,有科研竞赛获奖经历,希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
信息提示
请填处理意见