24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1776  |  回复: 1

torresqian

新虫 (初入文坛)

[求助] create_atom建立特殊形状边界的原子集合

使用lammps的create_atom命令外加关键词var跟set创建形状为三角函数的原子真空边界,可根据某一点坐标确定是否在该位置生成原子(The var and set keywords can be used together to provide a criterion for accepting or rejecting the addition of an individual atom, based on its coordinates)。我参照手册中的示例分别建立了二维跟三维的模型,但均出现了variable for create_atom is invalid style。原示例命令如下:
variable        x equal 100
variable        y equal 25
lattice           hex 0.8442
region          box block 0 $x 0 $y -0.5 0.5
create_box      1 box

variable        xx equal 0.0
variable        yy equal 0.0
variable        v equal "(0.2*v_y*ylat * cos(v_xx/xlat * 2.0*PI*4.0/v_x) + 0.5*v_y*ylat - v_yy) > 0.0"
create_atoms    1 box var v set x xx set y yy
我自己建模如下
3d:
variable        x equal 20
variable        y equal 10
variable        z equal 10
units            metal
atom_style  atomic
lattice          fcc 3.615
region         box block 0 $x 0 $y 0 $z
create_box  1 box
variable       xx equal 0.0
variable       yy equal 0.0
variable       zz equal 0.0
variable       v equal "(0.3*v_y*ylat-0.15*v_y*ylat*cos(v_xx/xlat*2.8*PI/v_x) &
-v_yy)>0.0"
create_atom 1 box var v set x xx set y yy set z zz
2d参照示例代码仍然提示variable for create_atom is invalid style。
不知道可有人做过类似的模型,能否指点一二?lammps小白,请大神们赐教
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TW九九

铜虫 (初入文坛)

我也是小白,现在在做粗糙度,需要建立特殊形状边界的原子集合,我也看了手册里关于creatatoms的介绍,你的示例的错误应该是在于creatatoms  set命令需要用internal变量,而你用的时equal变量,所以会出错,
关于我自己的问题我想问类似的命令可以建立锯齿状边界的原子吗,主要是手册中的那个示例我也没怎么看懂
2楼2021-03-18 13:50:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 torresqian 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求助调剂,跨调 +14 X十甫寸Y 2026-04-11 15/750 2026-04-14 00:36 by 王珺璞
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 085801电气专硕272求调剂 +12 电气李 2026-04-13 13/650 2026-04-13 21:12 by gretl
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +20 努力奋斗112 2026-04-08 20/1000 2026-04-13 14:24 by 张zhihao
[教师之家] 转长聘了 +3 简单化xn 2026-04-13 3/150 2026-04-13 14:18 by wwwkkk83
[考研] 211本科材料化工求调剂 +17 YHLAH 2026-04-11 20/1000 2026-04-13 08:46 by Sammy2
[考研] 299求调剂 +8 ZVVZ13 2026-04-08 8/400 2026-04-12 00:40 by 蓝云思雨
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
信息提示
请填处理意见