 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1758 | 回復(fù): 6 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
Ace_xt銅蟲 (初入文壇)
|
[求助]
表面摩擦做平衡的若干問題 已有1人參與
|
|
|
本人最近在用LAMMPS做摩擦，模型為z方向上下兩個表面，x、y、z方向都賦予周期性邊界條件，中間加水。我想在上下表面分別沿x軸正負方向運動的同時能在z方向做npt，使盒子在z方向有一定的伸縮以保證中間用于摩擦的水的壓強不變。但是在做平衡時，賦予所有原子初速度并在z方向npt，為上下表面賦予定時位移后模型沒能正常運動。請問可能是什么原因？另外腳本有無其它問題？下面是我做平衡的腳本，新手上路，請大家多多批評指正！ ########################## # Initialization #------------------------- boundary p p p units real #------------------------------ # Force field specifications #------------------------------ atom_style full bond_style harmonic angle_style charmm dihedral_style charmm pair_style lj/charmm/coul/charmm 10 12 pair_modify mix arithmetic neigh_modify every 2 delay 10 check yes special_bonds charmm #---------------- # Read geometry #---------------- read_data solved_d10.data #--------------------- # Define Groups #--------------------- group water type 1 10 group wato type 10 group allsulfur type 7 group uppersurf molecule <> 49 96 group lowersurf molecule <> 1 48 group uppersulfur intersect allsulfur uppersurf group lowersulfur intersect allsulfur lowersurf group uppernosulfur subtract uppersurf uppersulfur group lowernosulfur subtract lowersurf lowersulfur group pControl union lowernosulfur uppernosulfur water group twosurfs union lowersurf uppersurf #--------------------- # Fix all sulfur atoms #--------------------- fix fixsulfur allsulfur setforce 0.0 0.0 0.0 #-------------- # Minimization #-------------- min_style sd minimize 0.0 1.0e-8 1000 100000 #----------------------- # Initialize velocities #----------------------- velocity all create 10.0 4928459 rot yes mom yes dist gaussian #--------------------- # Shake all H atoms #--------------------- fix shakeH all shake 0.0001 20 0 m 1.0 a 11 #--------------------- ## Set Shear Velocity ##--------------------- variable xu equal vdisplace(0.0,0.0001) variable xl equal vdisplace(0.0,-0.0001) fix vu uppersurf move variable v_xu NULL NULL NULL NULL NULL fix vl lowersurf move variable v_xl NULL NULL NULL NULL NULL #-------------------------------- # Define the integration method #-------------------------------- fix 1 all npt temp 10.0 300.0 100.0 z 1.0 1.0 1000.0 #------------------------------- # Customize output infomration #------------------------------- thermo 10000 thermo_style multi dump mydump all atom 10000 trjequil.lammpstrj restart 500000 res.equil 附件中是模型視圖 |
回復(fù)此樓» 猜你喜歡
290分調(diào)劑求助
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有12人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有10人回復(fù)
-
本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分
已經(jīng)有5人回復(fù)
-
313求調(diào)劑
已經(jīng)有5人回復(fù)
-
362求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有27人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 超低界面張力表面活性劑驅(qū)油中， 泡沫穩(wěn)定性和界面張力之間的平衡關(guān)系
已經(jīng)有7人回復(fù)
- 多粗的玻璃管（或鐵管）往管頂端注水水不以水柱形式出現(xiàn)，而會順利流到管底？
已經(jīng)有5人回復(fù)
- 求助comsol模擬熱傳問題
已經(jīng)有11人回復(fù)
- 不銹鋼304如何減小表面摩擦系數(shù)呢？
已經(jīng)有6人回復(fù)
- 高溫摩擦磨損試驗后，試樣表面磨痕呈一圈凹坑，不知道怎么解釋
已經(jīng)有4人回復(fù)
- 土壤上非離子表面活性劑的吸附：為什么吸附量一直增加，達不到平衡
已經(jīng)有12人回復(fù)
- 無機非金屬版資源索引
已經(jīng)有63人回復(fù)
- 注塑成型的幾個問題
已經(jīng)有18人回復(fù)
- 問個問題：工材下面分幾個處，如何分管，代碼是多少
已經(jīng)有7人回復(fù)
Ace_xt
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時
- 蟲號: 2991731
- 注冊: 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學(xué)
送紅花一朵 
4樓2015-05-21 17:49:19
Ace_xt
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時
- 蟲號: 2991731
- 注冊: 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學(xué)
2樓2015-05-21 10:58:45
3樓2015-05-21 15:23:24
5樓2015-05-21 20:18:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | 圖鑒212 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +4 | zzll03 2026-03-31 | 4/200 |
|
|
[考研] 276求調(diào)劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +7 | 藍月亮亮 2026-03-30 | 7/350 |
|
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|