| 查看: 1665 | 回复: 6 | |||
[求助]
表面摩擦做平衡的若干问题 已有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 附件中是模型视图  |
回复此楼» 猜你喜欢
限项规定
已经有6人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
救命帖
已经有6人回复
-
存款400万可以在学校里躺平吗
已经有37人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
最失望的一年
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 超低界面张力表面活性剂驱油中， 泡沫稳定性和界面张力之间的平衡关系
已经有7人回复
- 多粗的玻璃管（或铁管）往管顶端注水水不以水柱形式出现，而会顺利流到管底？
已经有5人回复
- 求助comsol模拟热传问题
已经有11人回复
- 不锈钢304如何减小表面摩擦系数呢？
已经有6人回复
- 高温摩擦磨损试验后，试样表面磨痕呈一圈凹坑，不知道怎么解释
已经有4人回复
- 土壤上非离子表面活性剂的吸附：为什么吸附量一直增加，达不到平衡
已经有12人回复
- 无机非金属版资源索引
已经有63人回复
- 注塑成型的几个问题
已经有18人回复
- 问个问题：工材下面分几个处，如何分管，代码是多少
已经有7人回复
2楼2015-05-21 10:58:45
Ace_xt3楼2015-05-21 15:23:24
4楼2015-05-21 17:49:19
5楼2015-05-21 20:18:41
6楼2015-05-21 21:53:07
7楼2017-03-13 13:12:55