| 5 | 1/1 | 返回列表 |
| 查看: 3635 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
張牧童新蟲 (正式寫手)
|
[求助]
in文件運(yùn)行失。╨ost atoms)求助!!
|
||
|
剛剛接觸lammps沒多久,寫個(gè)小的in文件試試手,運(yùn)行中丟失原子,無(wú)法輸出想要的結(jié)果~求指導(dǎo)!最好請(qǐng)大神運(yùn)行下,給點(diǎn)建議!謝謝~ #3d contact simulation units metal dimension 3 boundary p p s atom_style atomic neighbor 3.0 bin neigh_modify delay 5 # create geometry region box block 0 80 0 80 0 100 units box create_box 3 box mass 1 64 mass 2 12 #atom region region boundary block INF INF INF INF INF 35 units box region coating block INF INF INF INF 35 40 units box region tool cone z 40 40 0 18 40 60 units box lattice fcc 3.6149 create_atoms 1 region boundary units box lattice diamond 3.57 create_atoms 2 region coating units box lattice diamond 3.57 create_atoms 3 region tool units box # Eam Potentials pair_style eam pair_coeff * * Ni_u3.eam # define groups group boundary region boundary group mobile subtract all boundary group tool region tool set group boundary type 1 set group mobile type 2 #initial velocities compute new mobile temp/partial 0 0 1 velocity mobile create 0.1 482748 temp new velocity tool set 0 0 -5 sum yes # fixes fix 1 all nve fix 2 boundary setforce 0 0 0 fix 3 mobile temp/rescale 100 0.1 0.1 0.02 1.0 fix_modify 3 temp new # Run timestep 0.01 thermo 100 thermo_modify temp new dump 1 all atom 500 dump.yahen.lammpstrj run 2000 |
木蟲 (著名寫手)
|
dimension 3 boundary s s p newton on units metal atom_style atomic neighbor 2.0 bin neigh_modify delay 5 lattice fcc 3.52 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box1 block 0 30 0 15 -3.0 3.0 side in units lattice region 1 sphere 15.0 20.0 0.0 2.0 side in units lattice region box2 block 0 30 15 22 -3.0 3.0 units lattice region box3 block 0 0.5 15 22 -3.0 3.0 region box4 block 29.5 30 15 22 -3.0 3.0 region box union 5 box1 1 box2 box3 box4 units lattice create_box 3 box mass 1 58.71 mass 2 58.71 mass 3 12 region 2 block INF INF INF 0.5 INF INF units lattice region 3 block INF 0.5 INF INF INF INF units lattice region 4 block 29.5 INF INF INF INF INF units lattice create_atoms 1 region box1 units lattice create_atoms 2 region 2 units lattice create_atoms 2 region 3 units lattice create_atoms 2 region 4 units lattice delete_atoms region box3 delete_atoms region box4 lattice none lattice diamond 3.57 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 create_atoms 3 region 1 units lattice pair_style eam pair_coeff * * Ni_u3.eam # define groups group sphere region 1 group lower region 2 group leftboundary region 3 group rightboundary region 4 group boundary union lower rightboundary leftboundary group mobile subtract all boundary set group mobile type 1 set group lower type 2 set group leftboundary type 2 set group rightboundary type 2 set group sphere type 3 # initial velocities compute new mobile temp velocity mobile create 0.1E-4 482748 temp new fix 1 all nve fix 2 boundary setforce 0.0 0.0 0.0 fix 3 all temp/rescale 100 0.1E-4 0.1E-4 0.1E-7 1.0 # run with indenter timestep 0.001 variable y equal "0.8*ylat - step*dt*1.6*ylat" fix 4 sphere move variable NULL v_y NULL NULL NULL NULL thermo 100 thermo_modify temp new compute stress mobile stress/atom pair compute displ mobile displace/atom compute centro mobile centro/atom fcc compute pe mobile pe/atom dump 1 all custom 100 dump.*.trajectory id type xs ys zs c_stress[1] c_stress[2] c_stress[3] c_centro c_pe c_displ[1] c_displ[2] c_displ[3] c_displ[4] run 3000 # run without indenter variable dy equal "step*dt*1.6*ylat-20.0" fix 6 sphere move variable NULL v_dy NULL NULL NULL NULL run 1000 我的源程序! |
木蟲 (正式寫手)
木蟲 (著名寫手)
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|