| 查看: 1491 | 回复: 1 | ||
18717718552新虫 (小有名气)
|
[求助]
分子动力学纳米压痕模拟数据处理问题 已有1人参与
|
|
这是我用来做纳米压痕的in文件 运行没有问题 但是数据处理的时候发现 位移增大 载荷反而减小 大神们帮忙看下是什么原因 echo both #creates matrix units metal boundary p p s atom_style atomic dimension 3 read_data 550ns0.1v823k.data1 mass 1 55.845 mass 2 63.546 mass 3 58.6934 timestep 0.002 pair_style eam/alloy pair_coeff * * FeCuNi.eam.alloy Fe Cu Ni #Energy minimization minimize 1.0e-3 1.0e-6 1000 1000 #define groups region reg block INF INF INF INF INF INF units box group reg region reg region reg1 block 9 90.282 9 90.282 9 50 units box group reg1 region reg1 region newton_layer block 20 79.282 20 79.282 20 100.282 units box group newton_layer region newton_layer group thermostat_layer subtract reg1 newton_layer group boundary_layer subtract reg reg1 dump 1 all atom 1000 dump123.lammpstrj log log123.data #intial velocities compute newtemp thermostat_layer temp velocity thermostat_layer create 300 8728007 temp newtemp #fix fix 1 boundary_layer setforce 0.0 0.0 0.0 #固定边界层 fix 2 thermostat_layer nvt temp 300.0 300.0 0.05 # 恒温层温度设置 #assigning velocity to the indenter in z direction loading variable z equal "101.282-step*dt*0.05" #0.05 loading rate print "z is $z" #indenter position and radius at onset of loading fix 3 newton_layer indent 1000.0 sphere 50 50 v_z 30.0 units box thermo 100 thermo_style custom step temp c_newtemp f_3[1] f_3[2] f_3[3] run 80000 |
回复此楼» 猜你喜欢
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有51人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有20人回复
-
考博自荐
已经有1人回复
-
山东大学第二批博士研究生招生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2020-06-13 21:20:47
5