 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 3462 | 回復(fù): 6 | |||
leoluode銀蟲 (小有名氣)
|
[求助]
Lammps拉伸應(yīng)變計算 已有1人參與
|
|
請問大家，如果輸入文件中timestep=0.001(ps)，拉伸變形命令 fix 1 all deform 100 x erate 0.0001 units box remap x , run 40000步，大家是怎么計算最終的應(yīng)變的�。� 按目前我的看法是，工程應(yīng)變速率為0.0001/ps,總共運(yùn)行40000步，每步的時間為0.001ps,那么總應(yīng)變?yōu)?.001ps*40000*0.0001/ps=0.0004=0.4%！但是感覺核deform后面的100沒有啥關(guān)系啊，換成1000，總應(yīng)變依然是0.4%! 所以我想問問大家對deform后面的100是如何理解的？ |
回復(fù)此樓» 猜你喜歡
308求調(diào)劑
已經(jīng)有3人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有8人回復(fù)
-
330求調(diào)劑
已經(jīng)有4人回復(fù)
-
304求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有6人回復(fù)
-
南京大學(xué)化學(xué)376求調(diào)劑
已經(jīng)有3人回復(fù)
-
考研調(diào)劑求學(xué)校推薦
已經(jīng)有5人回復(fù)
-
323求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- lammps碳納米管拉伸in文件報錯求助 急
已經(jīng)有1人回復(fù)
- lammps計算彈性模量應(yīng)力應(yīng)變
已經(jīng)有4人回復(fù)
- 為什么lammps中的應(yīng)變率這么大？
已經(jīng)有4人回復(fù)
- lammps做拉伸試驗(yàn)，求大神進(jìn)來指導(dǎo)下（烈日炎炎，有勞了`(*∩_∩*)′）
已經(jīng)有14人回復(fù)
- LAMMPS進(jìn)行拉伸模擬的一些問題
已經(jīng)有4人回復(fù)
- 關(guān)于lammps計算應(yīng)力——應(yīng)變的幾個問題
已經(jīng)有2人回復(fù)
- lammps做拉伸和壓縮的應(yīng)力問題
已經(jīng)有9人回復(fù)
- lammps里，關(guān)于應(yīng)力應(yīng)變的問題
已經(jīng)有0人回復(fù)
- LAMMPS模擬宏觀應(yīng)力可以實(shí)現(xiàn)嗎？
已經(jīng)有2人回復(fù)
- lammps實(shí)驗(yàn)拉伸Al納米線的in文件的錯誤提示好久找不出？
已經(jīng)有8人回復(fù)
- 關(guān)于lammps fix npt/deform做拉伸時溫度的變化
已經(jīng)有6人回復(fù)
- lammps 拉伸應(yīng)力應(yīng)變曲線
已經(jīng)有4人回復(fù)
- fix deform 指令拉伸求教
已經(jīng)有10人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
leoluode: 金幣+20, ★★★很有幫助 2016-07-27 11:43:20
感謝參與，應(yīng)助指數(shù) +1
leoluode: 金幣+20, ★★★很有幫助 2016-07-27 11:43:20
|
fix ID group-ID deform N parameter args ... keyword value ... N是每N步實(shí)施一次變形，如果這個數(shù)越大，那么每N步變形就越大，如果這個數(shù)較小，比如說是1，那么每跑一步就變形一次，這樣模擬的更細(xì)，精度更高。至于你說的百分之零點(diǎn)四，那是總變形吧，總變形是由后面的參數(shù)確定的，N影響不了總變形。 |
2樓2016-07-25 11:42:49
leoluode
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 346.5
- 帖子: 62
- 在線: 53.4小時
- 蟲號: 3487208
- 注冊: 2014-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
哈哈，多謝一直的回答��！ 關(guān)于您說的，我還是有點(diǎn)疑問！ （1）deform后面的N是每N 步變形一次，那么這個N設(shè)為1，就是每1步變形一次、總共運(yùn)行40000次，那么模型就變形4萬次，總變形時間為40000*0.001ps=40ps，應(yīng)變速率為0.0001/ps，則總變形為0.004； 如果N后面設(shè)為100，就是每100步變形一次，總共運(yùn)行40000次，實(shí)際上變形只有400次，則總共變形時間為400*0.001ps=0.4ps，應(yīng)變速率為0.0001/ps，則總變形為0.00004。 所以N是會影響總變形的； （2）deform后面的N為1或100的區(qū)別在于，N為1時，是每運(yùn)行1步，就變形1次；而N為100時，就是第一步變形了，第2-100步?jīng)]有變形，第101步變形，第102-200步?jīng)]有變形。第201步變形。。。。。。。 N為1或100區(qū)別就是變形后的弛豫時間步相同 不知道上面的理解您認(rèn)為有問題沒？ |
3樓2016-07-25 15:16:16
leoluode
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 346.5
- 帖子: 62
- 在線: 53.4小時
- 蟲號: 3487208
- 注冊: 2014-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
4樓2016-07-25 15:18:41
5樓2016-07-25 16:20:21
leoluode
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 346.5
- 帖子: 62
- 在線: 53.4小時
- 蟲號: 3487208
- 注冊: 2014-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
追夢2016 你好， 通過近一段時間的實(shí)際模擬測試，我發(fā)現(xiàn)關(guān)于拉伸應(yīng)變計算結(jié)果與你的解釋是一致的！ 對于timestep=0.001(ps)，拉伸變形命令 fix 1 all deform N x erate 0.0001 units box remap x , run 300000步， 總應(yīng)變是由erate=0.0001/ps、timestep=0.001ps、run 300000步三者決定，總應(yīng)變?yōu)?.03，與deform后的N值無關(guān)！ （1）N=1時，意味著每運(yùn)行一步，就變形一次，總共運(yùn)行300000步要達(dá)到變形量0.03，也就是說每次運(yùn)行的應(yīng)變?yōu)?0^-7，即每運(yùn)行一步，應(yīng)變增加10^-7，運(yùn)行到10000步時，總應(yīng)變?yōu)?0^-3;則在step=1~10000的過程中，總應(yīng)變依次為10^-7,2*10^-7,3*10^-7.........,10^-3 (2)N =10000時，意味著每運(yùn)行10000步，就變形一次，總共運(yùn)行300000步要達(dá)到變形量0.03，亦即在每個10000th變形一次，且應(yīng)變?yōu)?0^-3;;則在step=1~10000的過程中，總應(yīng)變依次為0,0 0,...，0,10^-3，step達(dá)到10000時，應(yīng)變一下子就由0增加大10^-3。 上述分析說明，實(shí)際上N=1和N=10000時，模型的實(shí)際變形速率確實(shí)是不一致的，N=1是每一步變形，變形速率為10^-4/ps; N=10000時在10000步之前都保持0應(yīng)變，而在每個10000 th時，直接變形達(dá)到規(guī)定值,變形速率為1/ps。N=10000時的實(shí)際變形速率是N=1的10000倍。 無論如何，erate后規(guī)定的應(yīng)變速率是一樣的（10^-4/ps），但是N改變時實(shí)際變形速率也改變了，只有當(dāng)N=1時，erate規(guī)定的變形速率才與實(shí)際變形速率是一樣的！那么變形速率我們在寫論文時到底說采用哪一個呢？ 還望“追夢2016”同學(xué)能繼續(xù)交流啊，能夠指出我理解上的錯誤 |
6樓2016-07-27 15:30:22
|
好的，首先感謝對我的信任。 我的理解是，咱們所說的變形速率（或者說論文時體現(xiàn)的變形速率）應(yīng)該是erate所規(guī)定的那個值�？赡�5樓解釋的不太準(zhǔn)確，這里N=1或=10000，只是告訴計算機(jī)，你多長時間（或者多少步）計算一次，但計算的依據(jù)是根據(jù)你說的erate、timestep和run共同決定的。只是說，如果N越小，我們就能得到越多的值，規(guī)律可能就越明顯。舉個例子，我們要找變形的中間規(guī)律，如果我們設(shè)定N=10000的話，那么用30000步得到0.03變形量我們得到30個數(shù)，可以畫出變形規(guī)律，但如果N=1的話，那將得到30000個值，這樣我們得到的規(guī)律就更準(zhǔn)確。 所以，我認(rèn)為，N值只是給計算機(jī)看的，他對我們研究的結(jié)論、屬性沒有什么關(guān)系。 當(dāng)然，僅為個人見解，可能會有錯誤，希望多加討論。 |
7樓2016-07-27 16:07:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 leoluode 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 279分求調(diào)劑 一志愿211 +10 | chaojifeixia 2026-03-19 | 11/550 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研]
11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|