| 查看: 2923 | 回復(fù): 10 | ||
2009云飛揚鐵桿木蟲 (著名寫手)
|
[求助]
lammps中某個區(qū)域壓強的求解
|
|
|
各位朋友,還是前面提到過的一個問題。我建立了一個立方體,中心挖去原子形成一個球孔,在球孔中填入了另一種原子A,想求解球孔區(qū)域A原子所形成的壓強是多大,目前知道應(yīng)該局部壓強采用compute stress/atom 來計算(單位是能量——壓強.體積),但是這樣就需要求解球孔區(qū)域A原子所形成的體積有多大,計算體積目前是采用的compute voronoi/atom 指令。但是在計算過程中總是覺得體積算的不對,相關(guān)的具體化指令如下,請知道的朋友幫忙指點指點,應(yīng)該怎么計算才能正確得到。先謝謝啦。 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) 體積計算部分是不是有問題呢?希望知道的朋友能夠詳細(xì)說說,先謝謝啦 |
木蟲 (正式寫手)
| 印象中compute voronoi/atom是唯一計算體系部分原子體積的方法,這個算法有些思路清奇,嚴(yán)格上說可能不太準(zhǔn)確,可我用過一次,覺得精度還行。但在用之前你總要知道這個命令輸出了什么東西吧。剛才簡單看了一下手冊,我有些記不清了,只能大致說一下當(dāng)時是怎么做的。1)voronoi體積不是真正的體積,所以你是不是要V(總體積)*voronoi(A原子)/voronoi(所有原子)來得到A原子的體積? 2)我記得這個命令輸出的不是單個的數(shù),而是數(shù)組,好像還是二維數(shù)組,所以你是不是應(yīng)該輸出到一個文本看看它的輸出格式,把其中一個數(shù)組所有的結(jié)果加起來,甚至為了結(jié)果準(zhǔn)確,還要做一個系綜平均 |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵蟲 (小有名氣)
新蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
木蟲 (正式寫手)
鐵桿木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | danyyyy 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
[考研] 求調(diào)劑,學(xué)校研究所都可以,材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩,總分302求調(diào)劑 +3 | 五維天空 2026-03-04 | 5/250 |
|
|
[考研] 中國科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑,有科研和獲獎經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +3 | T123 tt 2026-03-04 | 3/150 |
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里? +3 | 22735 2026-03-04 | 4/200 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研] 268求調(diào)劑 +6 | 好運連綿不絕 2026-03-02 | 6/300 |
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 化工專碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
[基金申請] 成果系統(tǒng)訪問量大,請一小時后再嘗試。---NSFC啥時候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | 材料專碩調(diào)劑; 2026-02-28 | 11/550 |
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|