24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2442  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 明天琳琳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工求调剂! +35 RichLi_ 2026-04-06 35/1750 2026-04-11 11:02 by zhq0425
[考研] 293求调剂 +6 勇远库爱314 2026-04-08 6/300 2026-04-11 10:08 by zhq0425
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 291分调剂 +3 上岸小莹加油 2026-04-09 4/200 2026-04-09 17:28 by 1753564080
[考研] 085404 293求调剂 +7 勇远库爱314 2026-04-08 7/350 2026-04-09 16:02 by 猪会飞
[考研] 262求调剂 +10 天下第一文 2026-04-04 13/650 2026-04-09 15:16 by 探123
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 280求调剂 +5 兮兮夜夜 2026-04-09 8/400 2026-04-09 11:15 by 兮兮夜夜
[考研] 一志愿北京2,材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
信息提示
请填处理意见