| 查看: 5135 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ayuancrystal铜虫 (初入文坛)
|
[求助]
新手求助！！lammps摩擦，摩擦力计算问题 已有3人参与
|
|
|
lammps小白，写了一个简单的单凸体接触模型，探针在基底上滑动，从网上找到一个摩擦例子，参照写了一个，可是发现摩擦力结果有问题，求助大神们 附上脚本： #initialization units metal boundary p p p dimension 3 atom_style atomic pair_style hybrid/overlay tersoff tersoff lj/cut 9 #Atom definition read_data Position.dat mass 1 28.0 mass 2 28.0 mass 3 28.0 mass 4 28.0 mass 5 28.0 mass 6 28.0 mass 7 28.0 mass 8 28.0 #Settings neighbor 1.0 bin neigh_modify delay 0 every 1 check yes group fixed-sub type 1 group temper-sub1 type 2 group temper-sub2 type 3 group interface-sub type 4 group bottom-tip type 5 group temper-tip1 type 6 group temper-tip2 type 7 group fixed-tip type 8 group low-tip union bottom-tip temper-tip1 temper-tip2 group hi-sub union temper-sub1 temper-sub2 interface-sub group boundary union fixed-tip fixed-sub group mobile subtract all boundary #Force field coefficients pair_coeff * * tersoff 1 Si.tersoff Si Si Si Si NULL NULL NULL NULL pair_coeff * * tersoff 2 Si.tersoff NULL NULL NULL NULL Si Si Si Si pair_coeff 1*4 5*8 lj/cut 0.0174 3.826 #Initial velocities compute tip_sub mobile temp/partial 0 1 1 velocity mobile create 298 34239 temp tip_sub compute 1 temper-sub1 temp compute 2 temper-sub2 temp compute 3 interface-sub temp compute 4 bottom-tip temp/partial 0 1 1 compute 5 temper-tip1 temp/partial 0 1 1 compute 6 temper-tip2 temp/partial 0 1 1 compute 7 hi-sub temp compute 8 low-tip temp/partial 0 1 1 #Fixes fix 1 driven nve fix 2 fixed-tip setforce 0.0 0.0 0.0 fix 3 fixed-sub setforce 0.0 0.0 0.0 fix 4 temper-sub1 langevin 298.0 298.0 0.00239 48279 fix 5 temper-tip1 langevin 298.0 298.0 0.00239 48279 #Run timestep 0.001 thermo 100 thermo_style custom step temp c_7 c_1 c_2 c_3 c_4 c_5 c_6 c_8 pe ke etotal vol f_2[1] f_2[2] f_2[3] f_3[1] f_3[2] f_3[3] thermo_modify format 1 %07d format 2 %7.3f format 3 %7.3f format 4 %7.3f format 5 %7.3f format 6 %7.3f format 7 %7.3f format 8 %7.3f format 9 %7.3f format 10 %7.3f format 11 %10.3f format 12 %10.3f format 13 %10.3f format 14 %10.3f format 15 %11.6f format 16 %11.6f format 17 %11.6f format 18 %11.6f format 19 %11.6f format 20 %11.6f dump 1 all xyz 1000 dump.tipsub.xyz dump_modify 1 element Si Si Si Si Si Si Si Si run 200000 fix 6 fixed-tip move linear 0.036 0.0 0.0 units box run 600000 参看摩擦例子，摩擦力取得是探针固定层原子受到的力。不知道这样计算摩擦力对不对，大神们是怎么计算摩擦力的，用的哪个命令啊? 还有处理的摩擦力如附图，摩擦力情况也不对啊！ 真心求指导啊！！！  摩擦力.png@oxox6085 |
回复此楼» 猜你喜欢
0703本科郑州大学求调剂
已经有3人回复
-
311求调剂
已经有4人回复
-
一志愿华理，数一英一285求A区调剂
已经有11人回复
-
调剂310
已经有5人回复
-
321求调剂
已经有6人回复
-
求调剂
已经有7人回复
-
349求调剂
已经有6人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有4人回复
-
材料求调剂
已经有9人回复
-
299求调剂
已经有7人回复
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
ayuancrystal: 金币+5, ★有帮助 2016-07-22 08:55:48
感谢参与，应助指数 +1
ayuancrystal: 金币+5, ★有帮助 2016-07-22 08:55:48
| fix 6命令放在thermo_style前面试试。另外，position_dat数据没问题吧？还有，在试验阶段，如果不知道程序对不对，建议别跑那么多步，少跑几步，如果没问题不出错了，再多跑点就行。 |
ayuancrystal
4楼2016-07-21 09:51:09
ayuancrystal
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160
- 帖子: 10
- 在线: 11.5小时
- 虫号: 3727254
- 注册: 2015-03-11
- 性别: MM
- 专业: 微/纳机械系统
3楼2016-07-21 09:04:30
ayuancrystal
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160
- 帖子: 10
- 在线: 11.5小时
- 虫号: 3727254
- 注册: 2015-03-11
- 性别: MM
- 专业: 微/纳机械系统
，
5楼2016-07-21 14:43:44
6楼2017-03-01 11:35:26
