| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2939 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
2009云飛揚(yáng)鐵桿木蟲(chóng) (著名寫(xiě)手)
|
[求助]
lammps中某個(gè)區(qū)域壓強(qiáng)的求解
|
||
|
各位朋友,還是前面提到過(guò)的一個(gè)問(wèn)題。我建立了一個(gè)立方體,中心挖去原子形成一個(gè)球孔,在球孔中填入了另一種原子A,想求解球孔區(qū)域A原子所形成的壓強(qiáng)是多大,目前知道應(yīng)該局部壓強(qiáng)采用compute stress/atom 來(lái)計(jì)算(單位是能量——壓強(qiáng).體積),但是這樣就需要求解球孔區(qū)域A原子所形成的體積有多大,計(jì)算體積目前是采用的compute voronoi/atom 指令。但是在計(jì)算過(guò)程中總是覺(jué)得體積算的不對(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ì)算部分是不是有問(wèn)題呢?希望知道的朋友能夠詳細(xì)說(shuō)說(shuō),先謝謝啦 |
木蟲(chóng) (正式寫(xiě)手)
| 印象中compute voronoi/atom是唯一計(jì)算體系部分原子體積的方法,這個(gè)算法有些思路清奇,嚴(yán)格上說(shuō)可能不太準(zhǔn)確,可我用過(guò)一次,覺(jué)得精度還行。但在用之前你總要知道這個(gè)命令輸出了什么東西吧。剛才簡(jiǎn)單看了一下手冊(cè),我有些記不清了,只能大致說(shuō)一下當(dāng)時(shí)是怎么做的。1)voronoi體積不是真正的體積,所以你是不是要V(總體積)*voronoi(A原子)/voronoi(所有原子)來(lái)得到A原子的體積? 2)我記得這個(gè)命令輸出的不是單個(gè)的數(shù),而是數(shù)組,好像還是二維數(shù)組,所以你是不是應(yīng)該輸出到一個(gè)文本看看它的輸出格式,把其中一個(gè)數(shù)組所有的結(jié)果加起來(lái),甚至為了結(jié)果準(zhǔn)確,還要做一個(gè)系綜平均 |
鐵桿木蟲(chóng) (著名寫(xiě)手)
鐵桿木蟲(chóng) (著名寫(xiě)手)
鐵桿木蟲(chóng) (著名寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +8 | mingmingry 2026-03-09 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +5 | 邱gl 2026-03-10 | 8/400 |
|
|
[考研] 材料工程307,求調(diào)劑 +7 | 我要燃燒你的夢(mèng) 2026-03-08 | 7/350 |
|
|
[考研] 327求調(diào)劑 +3 | Ffff03 2026-03-10 | 3/150 |
|
|
[考研] 0856材料與化工309分求調(diào)劑 +4 | ZyZy…… 2026-03-10 | 4/200 |
|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué),0856化工323分 | 本科應(yīng)化 | 有專利/競(jìng)賽/科研助手經(jīng)歷 | +7 | 橙子cyx 2026-03-06 | 9/450 |
|
|
[碩博家園] 2026級(jí)碩士研究生招生/調(diào)劑 +3 | 知足常樂(lè)的樂(lè) 2026-03-06 | 5/250 |
|
|
[考研] 一志愿山東大學(xué),總分327,英語(yǔ)二79,有論文,有競(jìng)賽,已過(guò)四六級(jí) +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
[考博] 26年博士申請(qǐng) +4 | 科研狗111 2026-03-07 | 4/200 |
|
|
[基金申請(qǐng)] 進(jìn)入個(gè)人成果庫(kù)好難,一下午都沒(méi)進(jìn)去 +6 | mi_dilee 2026-03-05 | 6/300 |
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過(guò)六級(jí)有競(jìng)賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) 求調(diào)劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 286 +6 | ksncj 2026-03-04 | 6/300 |
|