| 查看: 1888 | 回复: 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 附件中是模型视图  |
回复此楼» 猜你喜欢
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有11人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
面上函评意见出来了，像什么等级？
已经有14人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 超低界面张力表面活性剂驱油中， 泡沫稳定性和界面张力之间的平衡关系
已经有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
20
投票:
