 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1757 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
Ace_xt銅蟲 (初入文壇)
|
[求助]
表面摩擦做平衡的若干問題 已有1人參與
|
|
|
本人最近在用LAMMPS做摩擦，模型為z方向上下兩個(gè)表面，x、y、z方向都賦予周期性邊界條件，中間加水。我想在上下表面分別沿x軸正負(fù)方向運(yùn)動(dòng)的同時(shí)能在z方向做npt，使盒子在z方向有一定的伸縮以保證中間用于摩擦的水的壓強(qiáng)不變。但是在做平衡時(shí)，賦予所有原子初速度并在z方向npt，為上下表面賦予定時(shí)位移后模型沒能正常運(yùn)動(dòng)。請(qǐng)問可能是什么原因？另外腳本有無其它問題？下面是我做平衡的腳本，新手上路，請(qǐng)大家多多批評(píng)指正！ ########################## # 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ù)此樓» 猜你喜歡
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有23人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有12人回復(fù)
-
334分 一志愿武理 材料求調(diào)劑
已經(jīng)有16人回復(fù)
-
0703化學(xué)321分求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有7人回復(fù)
-
297求調(diào)劑
已經(jīng)有15人回復(fù)
-
11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè)
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有9人回復(fù)
-
化工專碩求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 超低界面張力表面活性劑驅(qū)油中， 泡沫穩(wěn)定性和界面張力之間的平衡關(guān)系
已經(jīng)有7人回復(fù)
- 多粗的玻璃管（或鐵管）往管頂端注水水不以水柱形式出現(xiàn)，而會(huì)順利流到管底？
已經(jīng)有5人回復(fù)
- 求助comsol模擬熱傳問題
已經(jīng)有11人回復(fù)
- 不銹鋼304如何減小表面摩擦系數(shù)呢？
已經(jīng)有6人回復(fù)
- 高溫摩擦磨損試驗(yàn)后，試樣表面磨痕呈一圈凹坑，不知道怎么解釋
已經(jīng)有4人回復(fù)
- 土壤上非離子表面活性劑的吸附：為什么吸附量一直增加，達(dá)不到平衡
已經(jīng)有12人回復(fù)
- 無機(jī)非金屬版資源索引
已經(jīng)有63人回復(fù)
- 注塑成型的幾個(gè)問題
已經(jīng)有18人回復(fù)
- 問個(gè)問題：工材下面分幾個(gè)處，如何分管，代碼是多少
已經(jīng)有7人回復(fù)
Ace_xt3樓2015-05-21 15:23:24
Ace_xt
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時(shí)
- 蟲號(hào): 2991731
- 注冊(cè): 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學(xué)
2樓2015-05-21 10:58:45
Ace_xt
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時(shí)
- 蟲號(hào): 2991731
- 注冊(cè): 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學(xué)
4樓2015-05-21 17:49:19
5樓2015-05-21 20:18:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +11 | 匹克i 2026-03-23 | 11/550 |
|
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|