 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1368 | 回復(fù): 4 | |||
mathsphy新蟲 (小有名氣)
|
[求助]
為什么兩個碰撞后能量發(fā)散了? 已有2人參與
|
|
大家好,我初學(xué)lammps,今天做兩個鐵原子碰撞時程序提示有原子丟失,所以我就使用了可變時間步長,結(jié)果發(fā)現(xiàn)它系統(tǒng)的能量不太對,所以想請教一下. 我這個系統(tǒng)就只建了一個5*5*5的晶胞,然后放了三個原子,坐標(biāo)分別是(0,0,0),(0,3,3),(0,3,0),晶格單位. 第一個原子是用作參考的,第二個原子碰撞第三個原子.系綜和設(shè)置如下: 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ā)散啊.....請問要怎么解決這個問題?...  |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 1.lammps
2015-11-28 23:12:36, 71.03 K
» 收錄本帖的淘帖專輯推薦
 分子動力學(xué) |
» 猜你喜歡
298求調(diào)劑一志愿211
已經(jīng)有3人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有8人回復(fù)
-
298求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有3人回復(fù)
-
275求調(diào)劑
已經(jīng)有8人回復(fù)
-
305分求調(diào)劑（食品工程）
已經(jīng)有6人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有17人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 量子光學(xué)為什么和經(jīng)典理論有矛盾
已經(jīng)有57人回復(fù)
- workbench顯示動力學(xué)碰撞分析問題 ：撞碎和能量太大
已經(jīng)有11人回復(fù)
- 如果一個直徑100英尺（約30米）的鉆石流星以光速撞向地球，會發(fā)生什么？（轉(zhuǎn))
已經(jīng)有8人回復(fù)
- 分子中電荷變化種種
已經(jīng)有82人回復(fù)
- 畢淑敏：旅行的33個理由
已經(jīng)有20人回復(fù)
- 30個寫在人臉上的密碼，為了自身健康務(wù)必保管好密碼
已經(jīng)有16人回復(fù)
- 【討論】大家看看吧！宇宙只是幻覺，哲學(xué)要騰飛了
已經(jīng)有63人回復(fù)
luftf
木蟲 (正式寫手)
- 應(yīng)助: 51 (初中生)
- 金幣: 3257.7
- 散金: 685
- 紅花: 19
- 帖子: 541
- 在線: 200.2小時
- 蟲號: 821742
- 注冊: 2009-08-04
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
2樓2015-11-29 20:18:00
mathsphy
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
3樓2015-11-30 10:08:57
feng0603
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2.9
- 散金: 12
- 帖子: 108
- 在線: 61.4小時
- 蟲號: 2973885
- 注冊: 2014-02-17
- 專業(yè): 微/納機械系統(tǒng)
4樓2015-12-10 11:02:16
mathsphy
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
5樓2015-12-10 23:39:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 mathsphy 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +8 | 梨花珞晚風(fēng) 2026-03-17 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分數(shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|