 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1367 | 回復: 4 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
mathsphy新蟲 (小有名氣)
|
[求助]
為什么兩個碰撞后能量發(fā)散了? 已有2人參與
|
||
|
大家好,我初學lammps,今天做兩個鐵原子碰撞時程序提示有原子丟失,所以我就使用了可變時間步長,結(jié)果發(fā)現(xiàn)它系統(tǒng)的能量不太對,所以想請教一下. 我這個系統(tǒng)就只建了一個5*5*5的晶胞,然后放了三個原子,坐標分別是(0,0,0),(0,3,3),(0,3,0),晶格單位. 第一個原子是用作參考的,第二個原子碰撞第三個原子.系綜和設置如下: fix 1 all nve fix 2 all dt/reset 1 NULL 0.001 0.01 velocity 1 set 1000 0 0 velocity 2 set 0 0 -1000 velocity 3 set 0 0 0 輸出的信息如下,其中最右邊三列分別代表第1,2,3原子的動能: Step Dt Time Temp TotEng PotEng KinEng 4 5 6 0 1e-005 0 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 10 1e-005 0.0001 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 20 1e-005 0.0002 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 30 1e-005 0.0003 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 40 1e-005 0.0004 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 50 1e-005 0.0005 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 60 1e-005 0.0006 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 70 1e-005 0.0007 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 80 1e-005 0.0008 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 90 1e-005 0.0009 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 100 1e-005 0.001 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 110 1e-005 0.0011 1.8351477e+008 47442.293 0 47442.293 23721.146 23721.146 0 120 1e-005 0.0012 1.8351477e+008 47442.293 -2.2147594e-005 47442.293 23721.146 23721.146 5.3606842e-015 130 9.9999999e-006 0.0013 1.8351478e+008 47442.293 -0.00057027847 47442.293 23721.146 23721.147 3.4424435e-012 140 9.9999994e-006 0.0014 1.8351478e+008 47442.293 -0.002693678 47442.295 23721.146 23721.149 7.6714312e-011 150 9.9999977e-006 0.0015 1.8351481e+008 47442.293 -0.010243151 47442.303 23721.146 23721.157 1.108068e-009 160 9.9999835e-006 0.0015999999 1.8351504e+008 47442.293 -0.067726385 47442.361 23721.146 23721.215 4.9021806e-008 170 9.9998924e-006 0.0016999994 1.8351669e+008 47442.293 -0.49445878 47442.787 23721.146 23721.641 2.5772336e-006 180 9.9997992e-006 0.0017999979 1.8351836e+008 47442.293 -0.92702466 47443.22 23721.146 23722.074 9.0716298e-006 190 9.9996775e-006 0.0018999954 1.8352057e+008 47442.293 -1.497937 47443.791 23721.146 23722.645 2.3667644e-005 200 9.9995503e-006 0.0019999916 1.8352292e+008 47442.293 -2.1048515 47444.398 23721.146 23723.252 4.6713654e-005 210 9.9994425e-006 0.0020999866 1.8352492e+008 47442.293 -2.6212737 47444.914 23721.146 23723.768 7.2437174e-005 220 9.9993716e-006 0.0021999806 1.8352627e+008 47442.293 -2.9705227 47445.263 23721.146 23724.117 9.3009162e-005 230 9.9993898e-006 0.0022999743 1.8352597e+008 47442.293 -2.8944083 47445.188 23721.146 23724.041 8.8322947e-005 240 9.9997456e-006 0.0023999692 1.8351944e+008 47442.287 -1.2123728 47443.5 23721.146 23722.353 1.5352359e-005 250 1e-005 0.0024999687 1.8347824e+008 47442.269 9.4212164 47432.847 23721.146 23711.7 0.00094058064 260 1e-005 0.0025999687 1.8330903e+008 47442.22 53.116248 47389.104 23721.146 23667.928 0.029882795 270 1e-005 0.0026999687 1.8275428e+008 47442.011 196.32033 47245.691 23721.146 23524.134 0.41077332 280 1e-005 0.0027999687 1.8045156e+008 47440.57 790.17944 46650.39 23721.146 22922.405 6.8394779 290 8.1905831e-012 0.0028299688 2.7355272e+020 7.0718929e+016 165.60334 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 300 8.1905831e-012 0.0028299689 2.7355272e+020 7.0718929e+016 6.8461037 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 310 8.1905831e-012 0.0028299689 2.7355272e+020 7.0718929e+016 -2.9300923 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 320 8.1905831e-012 0.002829969 2.7355272e+020 7.0718929e+016 -2.556017 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 330 8.1905831e-012 0.0028299691 2.7355272e+020 7.0718929e+016 -1.4125551 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 340 8.1905831e-012 0.0028299692 2.7355272e+020 7.0718929e+016 -0.44704995 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 350 8.1905831e-012 0.0028299693 2.7355272e+020 7.0718929e+016 -0.0086311757 7.0718929e+016 23721.146 3.5359435e+016 3.5359493e+016 我不明白為啥會發(fā)散啊.....請問要怎么解決這個問題?...  |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 1.lammps
2015-11-28 23:12:36, 71.03 K
» 收錄本帖的淘帖專輯推薦
 分子動力學 |
» 猜你喜歡
070300化學求調(diào)劑
已經(jīng)有3人回復
-
311求調(diào)劑
已經(jīng)有3人回復
-
生物學調(diào)劑
已經(jīng)有5人回復
-
269專碩求調(diào)劑
已經(jīng)有6人回復
-
08工科 320總分 求調(diào)劑
已經(jīng)有8人回復
-
尋找調(diào)劑
已經(jīng)有3人回復
-
291 求調(diào)劑
已經(jīng)有3人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有8人回復
-
一志愿東華大學化學070300，求調(diào)劑
已經(jīng)有6人回復
-
318求調(diào)劑
已經(jīng)有7人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 量子光學為什么和經(jīng)典理論有矛盾
已經(jīng)有57人回復
- workbench顯示動力學碰撞分析問題 ：撞碎和能量太大
已經(jīng)有11人回復
- 如果一個直徑100英尺（約30米）的鉆石流星以光速撞向地球，會發(fā)生什么？（轉(zhuǎn))
已經(jīng)有8人回復
- 分子中電荷變化種種
已經(jīng)有82人回復
- 畢淑敏：旅行的33個理由
已經(jīng)有20人回復
- 30個寫在人臉上的密碼，為了自身健康務必保管好密碼
已經(jīng)有16人回復
- 【討論】大家看看吧！宇宙只是幻覺，哲學要騰飛了
已經(jīng)有63人回復
mathsphy
新蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
5樓2015-12-10 23:39:21
luftf
木蟲 (正式寫手)
- 應助: 51 (初中生)
- 金幣: 3257.7
- 散金: 685
- 紅花: 19
- 帖子: 541
- 在線: 200.2小時
- 蟲號: 821742
- 注冊: 2009-08-04
- 性別: GG
- 專業(yè): 金屬材料的力學行為
2樓2015-11-29 20:18:00
mathsphy
新蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
3樓2015-11-30 10:08:57
feng0603
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 2.9
- 散金: 12
- 帖子: 108
- 在線: 61.4小時
- 蟲號: 2973885
- 注冊: 2014-02-17
- 專業(yè): 微/納機械系統(tǒng)
4樓2015-12-10 11:02:16
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿南大，0703化學，分數(shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 材料學碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 北科281學碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|