 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1729 | 回復: 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 附件中是模型視圖 |
回復此樓» 猜你喜歡
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有7人回復
-
265求調(diào)劑
已經(jīng)有5人回復
-
279分求調(diào)劑 一志愿211
已經(jīng)有3人回復
-
求調(diào)劑
已經(jīng)有3人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有12人回復
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有3人回復
-
329求調(diào)劑
已經(jīng)有5人回復
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
已經(jīng)有4人回復
-
材料學碩318求調(diào)劑
已經(jīng)有5人回復
-
一志愿中國海洋大學，生物學，301分，求調(diào)劑
已經(jīng)有6人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 超低界面張力表面活性劑驅(qū)油中， 泡沫穩(wěn)定性和界面張力之間的平衡關(guān)系
已經(jīng)有7人回復
- 多粗的玻璃管（或鐵管）往管頂端注水水不以水柱形式出現(xiàn)，而會順利流到管底？
已經(jīng)有5人回復
- 求助comsol模擬熱傳問題
已經(jīng)有11人回復
- 不銹鋼304如何減小表面摩擦系數(shù)呢？
已經(jīng)有6人回復
- 高溫摩擦磨損試驗后，試樣表面磨痕呈一圈凹坑，不知道怎么解釋
已經(jīng)有4人回復
- 土壤上非離子表面活性劑的吸附：為什么吸附量一直增加，達不到平衡
已經(jīng)有12人回復
- 無機非金屬版資源索引
已經(jīng)有63人回復
- 注塑成型的幾個問題
已經(jīng)有18人回復
- 問個問題：工材下面分幾個處，如何分管，代碼是多少
已經(jīng)有7人回復
Ace_xt
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時
- 蟲號: 2991731
- 注冊: 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學
2樓2015-05-21 10:58:45
Ace_xt3樓2015-05-21 15:23:24
Ace_xt
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時
- 蟲號: 2991731
- 注冊: 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學
4樓2015-05-21 17:49:19
5樓2015-05-21 20:18:41
Ace_xt
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 168.4
- 帖子: 16
- 在線: 12.2小時
- 蟲號: 2991731
- 注冊: 2014-02-24
- 性別: GG
- 專業(yè): 污染生態(tài)化學
6樓2015-05-21 21:53:07
7樓2017-03-13 13:12:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Ace_xt 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工 求調(diào)劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿，福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|