24小时热门版块排行榜

返回列表

Sebastiansxy

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 253.3
红花: 1
帖子: 37
在线: 13.1小时
虫号: 5023625
注册: 2016-09-17
性别: GG
专业: 机械动力学

[求助] 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