| 查看: 1869 | 回复: 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 附件中是模型视图  |
回复此楼» 猜你喜欢
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
不应该看fileCode
已经有8人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
-
重要来源：本周末出结果
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
好奇怪的filecode
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 超低界面张力表面活性剂驱油中， 泡沫稳定性和界面张力之间的平衡关系
已经有7人回复
- 多粗的玻璃管（或铁管）往管顶端注水水不以水柱形式出现，而会顺利流到管底？
已经有5人回复
- 求助comsol模拟热传问题
已经有11人回复
- 不锈钢304如何减小表面摩擦系数呢？
已经有6人回复
- 高温摩擦磨损试验后，试样表面磨痕呈一圈凹坑，不知道怎么解释
已经有4人回复
- 土壤上非离子表面活性剂的吸附：为什么吸附量一直增加，达不到平衡
已经有12人回复
- 无机非金属版资源索引
已经有63人回复
- 注塑成型的几个问题
已经有18人回复
- 问个问题：工材下面分几个处，如何分管，代码是多少
已经有7人回复
7楼2017-03-13 13:12:55
2楼2015-05-21 10:58:45
Ace_xt3楼2015-05-21 15:23:24
4楼2015-05-21 17:49:19
