| 查看: 2920 | 回復(fù): 10 | ||
2009云飛揚(yáng)鐵桿木蟲 (著名寫手)
|
[求助]
lammps中某個(gè)區(qū)域壓強(qiáng)的求解
|
|
|
各位朋友,還是前面提到過的一個(gè)問題。我建立了一個(gè)立方體,中心挖去原子形成一個(gè)球孔,在球孔中填入了另一種原子A,想求解球孔區(qū)域A原子所形成的壓強(qiáng)是多大,目前知道應(yīng)該局部壓強(qiáng)采用compute stress/atom 來計(jì)算(單位是能量——壓強(qiáng).體積),但是這樣就需要求解球孔區(qū)域A原子所形成的體積有多大,計(jì)算體積目前是采用的compute voronoi/atom 指令。但是在計(jì)算過程中總是覺得體積算的不對(duì),相關(guān)的具體化指令如下,請(qǐng)知道的朋友幫忙指點(diǎn)指點(diǎn),應(yīng)該怎么計(jì)算才能正確得到。先謝謝啦。 compute stress voidA stress/atom He_bubble compute stress1 voidA reduce sum c_stress[1] compute stress2 voidA reduce sum c_stress[2] compute stress3 voidA reduce sum c_stress[3] compute volum voidA voronoi/atom only_group variable p equal -(${c_stress1}+${c_stress2}+${c_stress3})/(3.0*c_volum) 體積計(jì)算部分是不是有問題呢?希望知道的朋友能夠詳細(xì)說說,先謝謝啦 |
木蟲 (正式寫手)
| 印象中compute voronoi/atom是唯一計(jì)算體系部分原子體積的方法,這個(gè)算法有些思路清奇,嚴(yán)格上說可能不太準(zhǔn)確,可我用過一次,覺得精度還行。但在用之前你總要知道這個(gè)命令輸出了什么東西吧。剛才簡(jiǎn)單看了一下手冊(cè),我有些記不清了,只能大致說一下當(dāng)時(shí)是怎么做的。1)voronoi體積不是真正的體積,所以你是不是要V(總體積)*voronoi(A原子)/voronoi(所有原子)來得到A原子的體積? 2)我記得這個(gè)命令輸出的不是單個(gè)的數(shù),而是數(shù)組,好像還是二維數(shù)組,所以你是不是應(yīng)該輸出到一個(gè)文本看看它的輸出格式,把其中一個(gè)數(shù)組所有的結(jié)果加起來,甚至為了結(jié)果準(zhǔn)確,還要做一個(gè)系綜平均 |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵蟲 (小有名氣)
新蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
木蟲 (正式寫手)
鐵桿木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[考研] 化學(xué)專業(yè)調(diào)劑 +3 | 好好好1233 2026-03-04 | 3/150 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +3 | T123 tt 2026-03-04 | 3/150 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
[考研] 一志愿鄭州大學(xué),學(xué)碩,物理化學(xué), 333求調(diào)劑 +5 | 李魔女斗篷 2026-03-04 | 5/250 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
[考研]
|
15779376950 2026-03-01 | 8/400 |
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
[考研] 0854總分272 +5 | 打小就是老實(shí)人 2026-03-02 | 6/300 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研] 化工專碩348,一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | 材料專碩調(diào)劑; 2026-02-28 | 11/550 |
|
|
[考研] 調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-02-28 | 3/150 |
|