| 查看: 7430 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
heibeidedong木虫 (正式写手)
|
[求助]
关于金刚石纳米压痕的in文件
|
|
|
最近写了一个金刚石纳米压痕文件，压头在基体上方2A处，现在设定压头沿着Z轴负方向移动10A，速度为0.01A/PS,时间步长设为0.001PS,按照这样算的话，10/0.01=1000PS,1000/0.001=1000000步，按照这样的计算，需要运行1000000步压头才能到达设定的深度，但事实上运行1800步压头就压入基体了，但是1800*0.001*0.01=0.018A<2A,这样的话压头还没进入基体，所以不知道什么原因，求解，谢谢。以下是in文件。 #3d indenter simulation units metal dimension 3 boundary p p s atom_style atomic neighbor 2.0 bin neigh_modify delay 5.0 lattice diamond 3.567 read_data data.0 mass 1 12.0 #Potential Setting newton on pair_style tersoff pair_coeff * * C.tersoff C #define groups region 1 block O 142.68 O 142.68 0 6.62 units box region 2 block O 142.68 O 142.68 6.62 20 units box region 3 block O 142.68 O 142.68 22 42.68 units box region 4 block O 142.68 O 142.68 O 42.68 units box group 11 region 1 group 22 region 2 group 33 region 3 group 44 region 4 set group 44 type 1 #define velcity fix 1 44 nve fix 2 33 rigid/nve single fix 3 11 setforce 0.0 0.0 0.0 velocity 33 set 0.0 0.0 -0.01 units box fix dong 33 setforce 0.0 0.0 -0.01 timestep 0.001 thermo 100 thermo_style custom step pe ke dump 1 all atom 100 dump.atom dump_modify 1 scale no run 1000000 #run without indenter velocity 33 set 0.0 0.0 0.01 units box unfix dong run 1000000 Sample Text |
回复此楼» 猜你喜欢
请问社科规划办官网上的“项目查询模块”与“项目申报模块”的账号密码一样吗？
已经有3人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有9人回复
-
40+中年人失业后，还能找到好工作吗？
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
极限拉扯，如何抉择
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
11楼2015-12-19 11:26:06
2楼2013-04-26 08:53:52
heibeidedong
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1481
- 散金: 258
- 红花: 2
- 帖子: 503
- 在线: 216.4小时
- 虫号: 2197207
- 注册: 2012-12-20
- 性别: GG
- 专业: 碳素材料与超硬材料
3楼2013-04-29 17:38:54
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
【答案】应助回帖
★ ★ ★ ★ ★
heibeidedong: 金币+5, ★★★很有帮助 2013-05-03 14:27:41
heibeidedong: 金币+5, ★★★很有帮助 2013-05-03 14:27:41
|
region 3 block O 142.68 O 142.68 22 42.68 units box region 4 block O 142.68 O 142.68 O 42.68 units box 这2句：区域4包含区域3 fix 1 44 nve fix 2 33 rigid/nve single 这2句， 对3都有作用， 所以压头实际设定速度为0.02A/PS fix dong 33 setforce 0.0 0.0 -0.01 为什么要加这个力？ 这样压头就一直在加速， 当然失控了。 你看看dump输出的位置，速度，力， 肯定有问题。 |
4楼2013-05-02 20:36:01