| 查看: 1744 | 回复: 8 | ||
ldj168168金虫 (正式写手)
|
[求助]
dump 问题
|
|
做了一个crack 的例子，源代码如下： # 2d LJ crack simulation(问题的基本初始化) dimension 2 #2维分子动力学模拟 boundary s s p #x,y方向不存在周期边界条件，z方向周期边界条件。 atom_style atomic #原子类型（金属） neighbor 0.3 bin #建近邻表参数（rc+0.3）bin表示为近邻表类型。 neigh_modify delay 5 #间隔多少载荷步重新形成近邻表 # create geometry 创建初始几何构形 lattice hex 0.93 #指定晶格类型（二维hex）和晶格常数 region box block 0 100 0 40 -0.25 0.25 #定义一个区域 create_box 5 box #在指定区域建立一个simulation box,5表示原子类型的种类数 create_atoms 1 box #在simulation box中创建类型为1的原子（原子位置初始化） mass 1 1.0 mass 2 1.0 mass 3 1.0 mass 4 1.0 mass 5 1.0 #指定单个原子的质量 # LJ potentials(指定原子作用势) pair_style lj/cut 2.5 #指定lj势，截断半径为2.5 pair_coeff * * 1.0 1.0 2.5 #指定lj势参数 # define groups（便于加载） region 1 block INF INF INF 1.25 INF INF group lower region 1 #定义lower组（便于施加外加速度） region 2 block INF INF 38.75 INF INF INF group upper region 2 #定义upper组（便于施加外加速度） group boundary union lower upper #定义总边界组 group mobile subtract all boundary #定义可动原子组（便于统计温度） region leftupper block INF 20 20 INF INF INF region leftlower block INF 20 INF 20 INF INF group leftupper region leftupper group leftlower region leftlower #定义左上、左下原子组（便于指定裂纹的存在） set group leftupper type 2 set group leftlower type 3 set group lower type 4 set group upper type 5 #指定原子类型（便于指定裂纹的存在） # initial velocities初始化速度 compute new mobile temp #定义温度的计算（可动区域内统计平均） #compute new2 mobile stress/atom #定义原子应力的计算（整个区域） velocity mobile create 0.01 887723 temp new #按指定的温度（0.01）计算方法，初始化原子的速度 velocity upper set 0.0 0.3 0.0 #upper原子组y方向的速度为0.3 velocity mobile ramp vy 0.0 0.3 y 1.25 38.75 sum yes #mobile原子的速初始度从0到0.3线性变化 # fixes施加约束 fix 1 all nve #nve系综的积分算法 fix 2 boundary setforce NULL 0.0 0.0 #边界boundary上力条件，钢化原子，便于加载！！ # run运行计算 timestep 0.003 #时间间隔步 thermo 200 #每200步输出热动力学统计量 thermo_modify temp new #计算温度通过new指示的方法计算 neigh_modify exclude type 2 3 #原子2，3之间作用取消（也就是通过不使他们在近邻表中出现实现） dump 1 all atom 500 dump.crack #每隔500步将原子信息写入文件dump.crack dump 2 mobile custom 500 dump2.crack tag x y z c_new2[2] run 5000 #进行5000步的模拟 结果总是输出ERROR：Invalid attribute in dump custom command (dump_custom.cpp:100) 这里的错误是tag 的问题吗？我试过多次，还是不行，哪位好心人指导下。 |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有9人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
Tanya_snail
木虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 3386.2
- 红花: 3
- 帖子: 298
- 在线: 333.3小时
- 虫号: 1656666
- 注册: 2012-03-01
- 性别: MM
- 专业: 金属材料的微观结构
2楼2013-09-03 19:02:53
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
ldj168168: 金币+4, ★★★★★最佳答案 2013-09-04 21:51:56
感谢参与，应助指数 +1
ldj168168: 金币+4, ★★★★★最佳答案 2013-09-04 21:51:56
|
请你参考lammps手册，dump命令里有“tag ”这个关键词吗？ http://lammps.sandia.gov/doc/dump.html "我试过多次" ，你再试一次，把tag 去掉 |
3楼2013-09-04 14:49:15
Tanya_snail
木虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 3386.2
- 红花: 3
- 帖子: 298
- 在线: 333.3小时
- 虫号: 1656666
- 注册: 2012-03-01
- 性别: MM
- 专业: 金属材料的微观结构
4楼2013-09-04 19:14:42
5楼2013-09-04 20:00:58
Tanya_snail
木虫 (小有名气)
- 应助: 25 (小学生)
- 金币: 3386.2
- 红花: 3
- 帖子: 298
- 在线: 333.3小时
- 虫号: 1656666
- 注册: 2012-03-01
- 性别: MM
- 专业: 金属材料的微观结构
6楼2013-09-04 20:07:11
ldj168168
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3425.1
- 散金: 9219
- 红花: 43
- 帖子: 625
- 在线: 136.2小时
- 虫号: 1958429
- 注册: 2012-08-27
- 专业: 半导体微纳机电器件与系统
7楼2013-09-04 21:18:16
ldj168168
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3425.1
- 散金: 9219
- 红花: 43
- 帖子: 625
- 在线: 136.2小时
- 虫号: 1958429
- 注册: 2012-08-27
- 专业: 半导体微纳机电器件与系统
8楼2013-09-04 21:19:56
ldj168168
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3425.1
- 散金: 9219
- 红花: 43
- 帖子: 625
- 在线: 136.2小时
- 虫号: 1958429
- 注册: 2012-08-27
- 专业: 半导体微纳机电器件与系统
9楼2013-09-04 21:54:17
5