| 5 | 1/1 | 返回列表 |
| 查看: 1789 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
lchunmei001木蟲 (正式寫手)
|
[求助]
沖擊結(jié)果異常
|
||
|
請教大家,我沿x方向給粒子2KM/s的速度沖擊加載,具體命令如下,但得到的結(jié)果很奇怪,壓強(qiáng)很低的情況下,溫度就已達(dá)到幾萬K,顯然很不合理,但是在不知道是哪里出問題了,測試了幾種物質(zhì),均出現(xiàn)這種情況,應(yīng)該是命令哪里沒對,還請高手指點一二,感激不盡! velocity all set 20 NULL NULL sum yes units box fix 1 all nve fix 2 all wall/reflect xlo xhi |
condensed matter physics |
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
|
輸入和結(jié)果如下,請予以指點,謝謝! 輸入: units metal boundary s p p lattice fcc 4.07 atom_style atomic region box block 0 100 0 10 0 10 create_box 1 box create_atoms 1 box pair_style eam pair_coeff * * Pt_u3.eam timestep 1.0e-3 neighbor 1.0 bin neigh_modify every 1 delay 10 check yes thermo 100 thermo_style custom step press vol temp pxx pyy pzz lx ly lz velocity all set 20 NULL NULL sum yes units box fix 1 all nve fix 2 all wall/reflect xlo xhi dump id all atom 10000 dump.melt run 50000 結(jié)果log文件: Step Press Volume Temp Pxx Pyy Pzz Lx Ly Lz 0 19146.265 674326.27 31285.986 535234.08 -238897.64 -238897.64 407.0814 40.7 40.7 100 23424.423 673363.58 31311.738 540451.5 -235089.11 -235089.11 406.50024 40.7 40.7 200 27653.106 671960.25 31331.714 546177.73 -231609.21 -231609.21 405.65307 40.7 40.7 300 31059.835 671084.06 31347.923 550838.79 -228829.64 -228829.64 405.12413 40.7 40.7 400 36428.574 669579.59 31362.412 557393.91 -224054.09 -224054.09 404.21589 40.7 40.7 500 39253.628 668750.35 31377.061 561471.39 -221855.25 -221855.25 403.7153 40.7 40.7 600 44925.431 667277.7 31394.038 568443.51 -216833.61 -216833.61 402.82628 40.7 40.7 700 47685.276 666310.63 31413.841 572717.7 -214830.93 -214830.93 402.24247 40.7 40.7 。。。。。。 |
|
樓主,我有點兒不懂,請教你哈: 1、你是要模擬一個什么樣的沖擊呢?我把你的in文件運(yùn)行了一下,4萬多個原子都一起向x方向運(yùn)動嗎?安裝我的理解,所謂的“沖擊”應(yīng)該是一個物體以較高的速度撞擊另一個物體,比如炮彈攻擊掩體的時候。我看過使用lammps模擬納米顆粒沖擊石墨烯的。等等這些。我沒弄明白你的模擬要干啥。 2、在你的in文件中,你能告訴我哪句命令是thermostat的?你設(shè)置的溫度是多少? PS 把dump命令改成 dump id all atom 10000 dump.lammpstrj 吧,這樣方便模擬完了使用VMD觀看軌跡動畫。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +10 | zhubinhao 2026-03-27 | 10/500 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|