 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 3461 | 回復: 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,總共運行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是如何理解的？ |
回復此樓» 猜你喜歡
287求調(diào)劑
已經(jīng)有8人回復
-
一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有4人回復
-
求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分
已經(jīng)有4人回復
-
085600材料與化工
已經(jīng)有8人回復
-
材料與化工 322求調(diào)劑
已經(jīng)有4人回復
-
藥學383 求調(diào)劑
已經(jīng)有5人回復
-
一志愿蘇州大學材料求調(diào)劑，總分315（英一）
已經(jīng)有5人回復
-
324求調(diào)劑
已經(jīng)有4人回復
-
一志愿南昌大學，327分，材料與化工085600
已經(jīng)有8人回復
-
350求調(diào)劑
已經(jīng)有5人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- lammps碳納米管拉伸in文件報錯求助 急
已經(jīng)有1人回復
- lammps計算彈性模量應(yīng)力應(yīng)變
已經(jīng)有4人回復
- 為什么lammps中的應(yīng)變率這么大？
已經(jīng)有4人回復
- lammps做拉伸試驗，求大神進來指導下（烈日炎炎，有勞了`(*∩_∩*)′）
已經(jīng)有14人回復
- LAMMPS進行拉伸模擬的一些問題
已經(jīng)有4人回復
- 關(guān)于lammps計算應(yīng)力——應(yīng)變的幾個問題
已經(jīng)有2人回復
- lammps做拉伸和壓縮的應(yīng)力問題
已經(jīng)有9人回復
- lammps里，關(guān)于應(yīng)力應(yīng)變的問題
已經(jīng)有0人回復
- LAMMPS模擬宏觀應(yīng)力可以實現(xiàn)嗎？
已經(jīng)有2人回復
- lammps實驗拉伸Al納米線的in文件的錯誤提示好久找不出？
已經(jīng)有8人回復
- 關(guān)于lammps fix npt/deform做拉伸時溫度的變化
已經(jīng)有6人回復
- lammps 拉伸應(yīng)力應(yīng)變曲線
已經(jīng)有4人回復
- fix deform 指令拉伸求教
已經(jīng)有10人回復
【答案】應(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ù)越大，那么每N步變形就越大，如果這個數(shù)較小，比如說是1，那么每跑一步就變形一次，這樣模擬的更細，精度更高。至于你說的百分之零點四，那是總變形吧，總變形是由后面的參數(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)于您說的，我還是有點疑問！ （1）deform后面的N是每N 步變形一次，那么這個N設(shè)為1，就是每1步變形一次、總共運行40000次，那么模型就變形4萬次，總變形時間為40000*0.001ps=40ps，應(yīng)變速率為0.0001/ps，則總變形為0.004； 如果N后面設(shè)為100，就是每100步變形一次，總共運行40000次，實際上變形只有400次，則總共變形時間為400*0.001ps=0.4ps，應(yīng)變速率為0.0001/ps，則總變形為0.00004。 所以N是會影響總變形的； （2）deform后面的N為1或100的區(qū)別在于，N為1時，是每運行1步，就變形1次；而N為100時，就是第一步變形了，第2-100步?jīng)]有變形，第101步變形，第102-200步?jīng)]有變形。第201步變形。。。。。。。 N為1或100區(qū)別就是變形后的弛豫時間步相同 不知道上面的理解您認為有問題沒？ |
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 你好， 通過近一段時間的實際模擬測試，我發(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時，意味著每運行一步，就變形一次，總共運行300000步要達到變形量0.03，也就是說每次運行的應(yīng)變?yōu)?0^-7，即每運行一步，應(yīng)變增加10^-7，運行到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時，意味著每運行10000步，就變形一次，總共運行300000步要達到變形量0.03，亦即在每個10000th變形一次，且應(yīng)變?yōu)?0^-3;;則在step=1~10000的過程中，總應(yīng)變依次為0,0 0,...，0,10^-3，step達到10000時，應(yīng)變一下子就由0增加大10^-3。 上述分析說明，實際上N=1和N=10000時，模型的實際變形速率確實是不一致的，N=1是每一步變形，變形速率為10^-4/ps; N=10000時在10000步之前都保持0應(yīng)變，而在每個10000 th時，直接變形達到規(guī)定值,變形速率為1/ps。N=10000時的實際變形速率是N=1的10000倍。 無論如何，erate后規(guī)定的應(yīng)變速率是一樣的（10^-4/ps），但是N改變時實際變形速率也改變了，只有當N=1時，erate規(guī)定的變形速率才與實際變形速率是一樣的！那么變形速率我們在寫論文時到底說采用哪一個呢？ 還望“追夢2016”同學能繼續(xù)交流啊，能夠指出我理解上的錯誤 |
6樓2016-07-27 15:30:22
|
好的，首先感謝對我的信任。 我的理解是，咱們所說的變形速率（或者說論文時體現(xiàn)的變形速率）應(yīng)該是erate所規(guī)定的那個值�？赡�5樓解釋的不太準確，這里N=1或=10000，只是告訴計算機，你多長時間（或者多少步）計算一次，但計算的依據(jù)是根據(jù)你說的erate、timestep和run共同決定的。只是說，如果N越小，我們就能得到越多的值，規(guī)律可能就越明顯。舉個例子，我們要找變形的中間規(guī)律，如果我們設(shè)定N=10000的話，那么用30000步得到0.03變形量我們得到30個數(shù)，可以畫出變形規(guī)律，但如果N=1的話，那將得到30000個值，這樣我們得到的規(guī)律就更準確。 所以，我認為，N值只是給計算機看的，他對我們研究的結(jié)論、屬性沒有什么關(guān)系。 當然，僅為個人見解，可能會有錯誤，希望多加討論。 |
7樓2016-07-27 16:07:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 leoluode 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 本人考085602 化學工程 專碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 化學求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿福大288有機化學，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿，福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 復試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|