 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 797 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者驕傲的九零后將贈(zèng)送您 50 個(gè)金幣 | ||
驕傲的九零后銀蟲(chóng) (正式寫(xiě)手)
|
[求助]
求助 lammps上DPD模擬計(jì)算
|
|
|
各位蟲(chóng)友好，最近在lammps上，用DPD計(jì)算模擬聚合物相分離，總是出現(xiàn)ERROR: Fix mvv/dpd requires atom attribute vest (src/DPD-MESO/fix_mvv_dpd.cpp:66)這句話。這里的vest是什么意思啊，找了不少資料也沒(méi)解決，在這里卡了好久了，哪位蟲(chóng)友遇到過(guò)啊，下面是我的in文件。 units lj dimension 3 boundary p p p timestep 0.005 ############## Atom and Bond Styles ############## atom_style full read_data Equil-2-5-80.data ############## Atom and Bond Styles ############## bond_style harmonic bond_coeff * 10.0 1.0 bond_coeff 4 100.0 1.0 bond_coeff 1 60.0 1.0 angle_style hybrid harmonic cosine angle_coeff 4 harmonic 1000 120 angle_coeff 1 cosine 1 angle_coeff 2 cosine 1 angle_coeff 3 cosine 1 angle_coeff 5 cosine 1 angle_coeff 6 cosine 1 angle_coeff 7 cosine 1 angle_coeff 8 cosine 1 angle_coeff 9 cosine 1 angle_coeff 10 cosine 1 ############## Groups and Fixes ############## #Set neighbor list radius to 10 neighbor 3 bin neigh_modify one 40000 neigh_modify page 400000 neigh_modify every 1 delay 0 check yes comm_modify vel yes # Fix beads of type 3 as rigid bodies group yellow type 3 4 fix 2 yellow rigid/small molecule langevin 1 1 0.5 456447 ############## Tem compute ############## group my type 1 2 5 region boundary block -24 28 -24 28 -28 24 compute mine my temp/region boundary ############## Output Settings ############## thermo 1000 thermo_style custom step temp time ebond eangle edihed pe ke etotal lx ly lz density ############## Minimization ################ velocity all create 1 34387 #minimize 1e-6 1e-6 1000 1000 #write_data min.data ############## Soft Equilibration ############## dump 1 all custom 1000 pe.lammpstrj id mass mol type xu yu zu dump_modify 1 flush yes ############## Pair Styles ############## # Soft potential for initial relaxation pair_style soft 2.5 pair_coeff * * 10 # Initial NvT equilibration fix 3 all nvt temp 1 1 0.5 run 100000 unfix 3 write_data SEquil-2-5-80.data ############## DPD Equilibration ################ pair_style hybrid dpd 1 1 34387 lj/cut 2.5 #pair_style dpd T cutoff seed #pair_style style args #lj/cut args = cutoff pair_coeff 1 1 dpd 25 4.5 pair_coeff 1 2 dpd 60 4.5 pair_coeff 1 3 dpd 60 4.5 pair_coeff 1 4 dpd 60 4.5 pair_coeff 1 5 dpd 60 4.5 pair_coeff 2 2 dpd 25 4.5 pair_coeff 2 3 dpd 80 4.5 pair_coeff 2 4 dpd 80 4.5 pair_coeff 2 5 dpd 25 4.5 pair_coeff 3 3 dpd 25 4.5 pair_coeff 3 4 dpd 25 4.5 pair_coeff 3 5 dpd 80 4.5 pair_coeff 4 4 dpd 25 4.5 pair_coeff 4 5 dpd 80 4.5 pair_coeff 5 5 dpd 25 4.5 pair_coeff 3 3 lj/cut 2 1.2 pair_coeff 4 4 lj/cut 2 1.2 fix 4 all langevin 1 1 0.5 12345 ####################################################################Deform to the density of 3 g/cm^3########### fix 11 all mvv/dpd 0.65 run 1000000 write_data deformed.data run 1000000 write_data 1deformed.data run 100000 unfix 11 write_data 2deformed.data |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有181人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有5人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有2人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 驕傲的九零后 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +5 | Seele_v5 2026-03-04 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|
|
[考研] 中國(guó)科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑，有科研和獲獎(jiǎng)經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
[考研] 材料與化工，291分，求調(diào)劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 3/150 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛(ài)讀書(shū)的小帥 2026-02-28 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 1 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|