| 查看: 1782 | 回復: 11 | |||
lchunmei001木蟲 (正式寫手)
|
[求助]
沖擊結果異常
|
|
請教大家,我沿x方向給粒子2KM/s的速度沖擊加載,具體命令如下,但得到的結果很奇怪,壓強很低的情況下,溫度就已達到幾萬K,顯然很不合理,但是在不知道是哪里出問題了,測試了幾種物質(zhì),均出現(xiàn)這種情況,應該是命令哪里沒對,還請高手指點一二,感激不盡! velocity all set 20 NULL NULL sum yes units box fix 1 all nve fix 2 all wall/reflect xlo xhi |
condensed matter physics |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
輸入和結果如下,請予以指點,謝謝! 輸入: 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 結果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文件運行了一下,4萬多個原子都一起向x方向運動嗎?安裝我的理解,所謂的“沖擊”應該是一個物體以較高的速度撞擊另一個物體,比如炮彈攻擊掩體的時候。我看過使用lammps模擬納米顆粒沖擊石墨烯的。等等這些。我沒弄明白你的模擬要干啥。 2、在你的in文件中,你能告訴我哪句命令是thermostat的?你設置的溫度是多少? PS 把dump命令改成 dump id all atom 10000 dump.lammpstrj 吧,這樣方便模擬完了使用VMD觀看軌跡動畫。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
![]() |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 考研292求調(diào)劑 +3 | YIXIAO789 2026-03-26 | 3/150 |
|
|
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 333求調(diào)劑 +7 | 87639 2026-03-21 | 12/600 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 蘆lty 2026-03-25 | 6/300 |
|
|
[考研] 0703化學求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 機械學碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 070300化學求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 一志愿南大,0703化學,分數(shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|