 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 802 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問題，作者驕傲的九零后將贈送您 50 個(gè)金幣 | |||
驕傲的九零后銀蟲 (正式寫手)
|
[求助]
求助 lammps上DPD模擬計(jì)算
|
||
|
各位蟲友好，最近在lammps上，用DPD計(jì)算模擬聚合物相分離，總是出現(xiàn)ERROR: Fix mvv/dpd requires atom attribute vest (src/DPD-MESO/fix_mvv_dpd.cpp:66)這句話。這里的vest是什么意思啊，找了不少資料也沒解決，在這里卡了好久了，哪位蟲友遇到過啊，下面是我的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é)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有100人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 驕傲的九零后 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 0856材料與化工353分求調(diào)劑 +10 | NIFFFfff 2026-03-09 | 10/500 |
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
|
[考研] 0703求調(diào)劑 +3 | jtyq001 2026-03-10 | 3/150 |
|
|
|
[考博] 找博導(dǎo) +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
|
[碩博家園] 2026級碩士研究生招生/調(diào)劑 +3 | 知足常樂的樂 2026-03-06 | 5/250 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +11 | 鶴遨予卿 2026-03-04 | 13/650 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +5 | 鋼鐵大炮 2026-03-04 | 5/250 |
|
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
|
[考研] 材料工程330分求調(diào)劑，一志愿985 +3 | 小材化本科 2026-03-07 | 3/150 |
|
|
|
[考研] 083000，總分284，求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
|
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國金及主持省級大創(chuàng)經(jīng)歷 +6 | 熊境喆 2026-03-04 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|