| 查看: 2458 | 回复: 3 | ||
Sebastiansxy银虫 (初入文坛)
|
[求助]
lammps中shear例子，关于划定区域region的一些疑问
|
|
求助各位大神，最近想做剪切模拟，看了lammps上的shear例子，但是对于里面的区域划定感觉到很困惑，求大神指点，下面是代码： units metal boundary s s p atom_style atomic lattice fcc 3.52 region box block 0 16.0 0 10.0 0 2.828427 create_box 3 box lattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 & origin 0.5 0 0 create_atoms 1 box pair_style eam pair_coeff * * Ni_u3.eam neighbor 0.3 bin neigh_modify delay 5 region lower block INF INF INF 0.9 INF INF region upper block INF INF 6.1 INF INF INF #这里不太懂，前面划分区域y方向是0-10，但是这里划分的不对称啊，但是最后生成的区域是对称的 group lower region lower group upper region upper group boundary union lower upper group mobile subtract all boundary set group lower type 2 set group upper type 3 compute new3d mobile temp compute new2d mobile temp/partial 0 1 1 # equilibrate velocity mobile create 300.0 5812775 temp new3d fix 1 all nve fix 2 boundary setforce 0.0 0.0 0.0 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new3d thermo 25 thermo_modify temp new3d timestep 0.001 run 100 # shear velocity upper set 1.0 0 0 velocity mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes #这里赋予区域均匀变化的速度，但是区域y方向是1.4-8.6，按照前面的思路，不应该是0.9-6.1嘛，这个不是很懂 unfix 3 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new2d dump 1 all custom 100 dump.* id type xs ys zs ix iy iz thermo 100 thermo_modify temp new2d reset_timestep 0 run 3000 对于这个问题，考虑到了lattice中坐标轴方向的问题，但怎么算都觉得不对，望恳请大神指点 下面是划分好区域的图，和运行到最后的图  1.png  2.png |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有255人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
nicoleguozi
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5848.2
- 散金: 200
- 红花: 4
- 帖子: 532
- 在线: 217.7小时
- 虫号: 1253837
- 注册: 2011-04-02
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2018-07-31 18:18:28
Sebastiansxy
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 253.3
- 红花: 1
- 帖子: 37
- 在线: 13.1小时
- 虫号: 5023625
- 注册: 2016-09-17
- 性别: GG
- 专业: 机械动力学
3楼2018-11-13 13:34:13
4楼2020-05-25 00:07:06
