 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1369 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
mathsphy新蟲 (小有名氣)
|
[求助]
為什么兩個(gè)碰撞后能量發(fā)散了? 已有2人參與
|
||
|
大家好,我初學(xué)lammps,今天做兩個(gè)鐵原子碰撞時(shí)程序提示有原子丟失,所以我就使用了可變時(shí)間步長,結(jié)果發(fā)現(xiàn)它系統(tǒng)的能量不太對,所以想請教一下. 我這個(gè)系統(tǒng)就只建了一個(gè)5*5*5的晶胞,然后放了三個(gè)原子,坐標(biāo)分別是(0,0,0),(0,3,3),(0,3,0),晶格單位. 第一個(gè)原子是用作參考的,第二個(gè)原子碰撞第三個(gè)原子.系綜和設(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原子的動(dòng)能: 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 我不明白為啥會(huì)發(fā)散啊.....請問要怎么解決這個(gè)問題?...  |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(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
» 收錄本帖的淘帖專輯推薦
 分子動(dòng)力學(xué) |
» 猜你喜歡
一志愿武理材料工程348求調(diào)劑
已經(jīng)有6人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑
已經(jīng)有6人回復(fù)
-
323求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有5人回復(fù)
-
317求調(diào)劑
已經(jīng)有16人回復(fù)
-
287求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
308求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 量子光學(xué)為什么和經(jīng)典理論有矛盾
已經(jīng)有57人回復(fù)
- workbench顯示動(dòng)力學(xué)碰撞分析問題 ：撞碎和能量太大
已經(jīng)有11人回復(fù)
- 如果一個(gè)直徑100英尺（約30米）的鉆石流星以光速撞向地球，會(huì)發(fā)生什么？（轉(zhuǎn))
已經(jīng)有8人回復(fù)
- 分子中電荷變化種種
已經(jīng)有82人回復(fù)
- 畢淑敏：旅行的33個(gè)理由
已經(jīng)有20人回復(fù)
- 30個(gè)寫在人臉上的密碼，為了自身健康務(wù)必保管好密碼
已經(jīng)有16人回復(fù)
- 【討論】大家看看吧！宇宙只是幻覺，哲學(xué)要騰飛了
已經(jīng)有63人回復(fù)
feng0603
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2.9
- 散金: 12
- 帖子: 108
- 在線: 61.4小時(shí)
- 蟲號: 2973885
- 注冊: 2014-02-17
- 專業(yè): 微/納機(jī)械系統(tǒng)
4樓2015-12-10 11:02:16
luftf
木蟲 (正式寫手)
- 應(yīng)助: 51 (初中生)
- 金幣: 3257.7
- 散金: 685
- 紅花: 19
- 帖子: 541
- 在線: 200.2小時(shí)
- 蟲號: 821742
- 注冊: 2009-08-04
- 性別: GG
- 專業(yè): 金屬材料的力學(xué)行為
2樓2015-11-29 20:18:00
mathsphy
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時(shí)
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
3樓2015-11-30 10:08:57
mathsphy
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 419
- 帖子: 94
- 在線: 8.9小時(shí)
- 蟲號: 3485761
- 注冊: 2014-10-20
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
5樓2015-12-10 23:39:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[考研] 招08考數(shù)學(xué) +3 | laoshidan 2026-03-20 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|