| 5 | 1/1 | 返回列表 |
| 查看: 3465 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
leoluode銀蟲 (小有名氣)
|
[求助]
Lammps拉伸應(yīng)變計(jì)算 已有1人參與
|
||
|
請(qǐng)問大家,如果輸入文件中timestep=0.001(ps),拉伸變形命令 fix 1 all deform 100 x erate 0.0001 units box remap x , run 40000步,大家是怎么計(jì)算最終的應(yīng)變的啊? 按目前我的看法是,工程應(yīng)變速率為0.0001/ps,總共運(yùn)行40000步,每步的時(shí)間為0.001ps,那么總應(yīng)變?yōu)?.001ps*40000*0.0001/ps=0.0004=0.4%!但是感覺核deform后面的100沒有啥關(guān)系啊,換成1000,總應(yīng)變依然是0.4%! 所以我想問問大家對(duì)deform后面的100是如何理解的? |
銅蟲 (小有名氣)
|
好的,首先感謝對(duì)我的信任。 我的理解是,咱們所說的變形速率(或者說論文時(shí)體現(xiàn)的變形速率)應(yīng)該是erate所規(guī)定的那個(gè)值。可能5樓解釋的不太準(zhǔn)確,這里N=1或=10000,只是告訴計(jì)算機(jī),你多長時(shí)間(或者多少步)計(jì)算一次,但計(jì)算的依據(jù)是根據(jù)你說的erate、timestep和run共同決定的。只是說,如果N越小,我們就能得到越多的值,規(guī)律可能就越明顯。舉個(gè)例子,我們要找變形的中間規(guī)律,如果我們?cè)O(shè)定N=10000的話,那么用30000步得到0.03變形量我們得到30個(gè)數(shù),可以畫出變形規(guī)律,但如果N=1的話,那將得到30000個(gè)值,這樣我們得到的規(guī)律就更準(zhǔn)確。 所以,我認(rèn)為,N值只是給計(jì)算機(jī)看的,他對(duì)我們研究的結(jié)論、屬性沒有什么關(guān)系。 當(dāng)然,僅為個(gè)人見解,可能會(huì)有錯(cuò)誤,希望多加討論。 |

銅蟲 (小有名氣)
|
fix ID group-ID deform N parameter args ... keyword value ... N是每N步實(shí)施一次變形,如果這個(gè)數(shù)越大,那么每N步變形就越大,如果這個(gè)數(shù)較小,比如說是1,那么每跑一步就變形一次,這樣模擬的更細(xì),精度更高。至于你說的百分之零點(diǎn)四,那是總變形吧,總變形是由后面的參數(shù)確定的,N影響不了總變形。 |

銀蟲 (小有名氣)
|
哈哈,多謝一直的回答。 關(guān)于您說的,我還是有點(diǎn)疑問! (1)deform后面的N是每N 步變形一次,那么這個(gè)N設(shè)為1,就是每1步變形一次、總共運(yùn)行40000次,那么模型就變形4萬次,總變形時(shí)間為40000*0.001ps=40ps,應(yīng)變速率為0.0001/ps,則總變形為0.004; 如果N后面設(shè)為100,就是每100步變形一次,總共運(yùn)行40000次,實(shí)際上變形只有400次,則總共變形時(shí)間為400*0.001ps=0.4ps,應(yīng)變速率為0.0001/ps,則總變形為0.00004。 所以N是會(huì)影響總變形的; (2)deform后面的N為1或100的區(qū)別在于,N為1時(shí),是每運(yùn)行1步,就變形1次;而N為100時(shí),就是第一步變形了,第2-100步?jīng)]有變形,第101步變形,第102-200步?jīng)]有變形。第201步變形。。。。。。。 N為1或100區(qū)別就是變形后的弛豫時(shí)間步相同 不知道上面的理解您認(rèn)為有問題沒? |
銀蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|