24小时热门版块排行榜    

查看: 1860  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见