| 5 | 1/1 | 返回列表 |
| 查看: 2922 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
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ì)算過程中總是覺得體積算的不對,相關(guān)的具體化指令如下,請知道的朋友幫忙指點(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è)命令輸出了什么東西吧。剛才簡單看了一下手冊,我有些記不清了,只能大致說一下當(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è)系綜平均 |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)