| 5 | 1/1 | 返回列表 |
| 查看: 2434 | 回復(fù): 4 | ||
明天琳琳新蟲 (初入文壇)
|
[求助]
納米壓痕相關(guān)的lammps程序(新手上路) 已有2人參與
|
|
我剛剛學(xué)lammps沒多久,自己剛剛弄好了一個(gè)程序,怎么改都錯(cuò)了好多,特請(qǐng)教大神幫忙提點(diǎn)錯(cuò)誤…… 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 ![]() ![]() ![]() |
新蟲 (初入文壇)
新蟲 (初入文壇)
|
create_box命令用多了,應(yīng)該是只有一個(gè)box nvt系綜的溫度設(shè)置為0.01 在box1到3 的region里建立原子,如果分別在三個(gè)區(qū)域中建立原子,可能會(huì)產(chǎn)生原子距離過近的問題,所以可以用 region union的命令把三個(gè)region整合為一個(gè)再建立原子(region之間有交集關(guān)系并無影響) 暫時(shí)就想到這么多,不一定完全正確,僅做參考 fix indent這玩意我也不太懂,如果有什么進(jìn)展,還請(qǐng)不吝賜教! |
新蟲 (初入文壇)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +5 | shansx 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 招08考數(shù)學(xué) +3 | laoshidan 2026-03-20 | 7/350 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +4 | 皮皮皮123456 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|