24小时热门版块排行榜    

查看: 3630  |  回复: 9

xiaorongwei

金虫 (初入文坛)

[求助] AtomEye下面的voronoirize的使用

AtomEye下面可以用voronoi方法建立多晶,http://li.mit.edu/Archive/Graphics/A/utils.html#voronoirize
作者给的例子中,用了下面的命令:
% voronoirize FCC10x10x10.cfg bcc 2 2 2 FCC10x10x10-voronoirized.cfg
现在我有FCC10x10x10.cfg 文件,问题是我不知道这句话“% voronoirize FCC10x10x10.cfg bcc 2 2 2 FCC10x10x10-voronoirized.cfg”怎么用,
哪位大侠知道啊,跪求!!!
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学第一性原理 材料

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

rush7

金虫 (小有名气)

我不太懂,不过MDBBS上岂不是有更好的资源?atomeye不是可视化么,LiJu的意思是在那里面建模。。?另外LZ要是有voronoirize高版本的可否分享一下~~~?
2楼2012-06-25 16:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaorongwei

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by rush7 at 2012-06-25 16:47:14
我不太懂,不过MDBBS上岂不是有更好的资源?atomeye不是可视化么,LiJu的意思是在那里面建模。。?另外LZ要是有voronoirize高版本的可否分享一下~~~?

是可视化的软件,不过在atomeye网页上,voronoirize是属于 “Utilities”下面的,所以我认为是能voronoirize建模的。你说的我有“voronoirize高版本” 我没有哦,难道你有关于voronoirize的版本?能不能法我参考一下?
3楼2012-06-25 16:57:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icic2012

金虫 (著名写手)

感谢参与,应助指数 +1
jiaoyixiong: 应助指数-1 2012-06-26 07:32:15
现在党 员没有党 原证。除非组织证明,否则还真不好办。
从来不知道当有什么秘密,也不知道有什么几率
4楼2012-06-25 17:52:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rush7

金虫 (小有名气)

引用回帖:
3楼: Originally posted by xiaorongwei at 2012-06-25 16:57:11
是可视化的软件,不过在atomeye网页上,voronoirize是属于 “Utilities”下面的,所以我认为是能voronoirize建模的。你说的我有“voronoirize高版本” 我没有哦,难道你有关于voronoirize的版本?能不能法我参考一 ...

没呀 只记得MDBBS上有人说过不止一个版本……
5楼2012-06-25 22:23:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chen_kg

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
chaizhm: 金币+2, 谢谢~ 2012-08-07 19:36:58
xiaorongwei: 金币+5, ★★★很有帮助 2012-09-18 15:48:14
voronoirize 输入结构文件 晶体结构 x方向晶粒数目 y方向晶粒数目 z方向晶粒数目 输出结构文件
对应的FCC10x10x10.cfg是输入文件,bcc是结构 2 2 2表示将初始结构随机切成2*2*2个晶粒, FCC10x10x10-voronoirized.cfg是输出结构文件,下面的代码可以将这个cfg直接转成lammps的data文件
#! /bin/gawk -f
BEGIN {
OFMT = "%.16f";
CONVFMT ="%.16f";
count = 0;
H11 = 0.;
H22 = 0.;
H33 = 0.;
entry_count = 0;
N = 0.;
}
function abs(y){
    y =(y >= 0) ? y : 0 - y;
    return(y);
}


{}{
    if ($1 == "Number" && $2 == "of" && $3 == "particles" && NR == 1){
        N = $5;
        printf ("Input file\n\n";
        printf ("%d  atoms\n", $5);
        printf ("1  atom types\n\n";
    }
    if ( $1 == "H0(1,1)" {H11 = $3; printf ("0.0000 %.16f xlo xhi\n", H11);}
    if ( $1 == "H0(2,2)" {H22 = $3; printf ("0.0000 %.16f ylo yhi\n", H22);}
    if ( $1 == "H0(3,3)" {
        H33 = $3;
        printf ("0.0000 %.16f zlo zhi\n\n", H33);
        printf ("Masses\n\n";
        printf ("1 63.546000\n\n\n";
        printf ("Atoms\n\n";
    }
    if ( $0 == ".NO_VELOCITY." type = "extended";
    if ( type == "extended" && $1 == "entry_count" entry_count = $3;
    if ( type == "extended" && NF == entry_count && entry_count > 0 && NR > 27) {
        printf("%d 1 %.16f %.16f %.16f %.12f %.12f\n",count+1,$1*H11,$2*H22,$3*H33, 1./sqrt($4),$5);count++;
    }

    if ( $0 == "R = 1.0 [ns^-1]" type = "standard";
    if ( type == "standard" entry_count = 8;
    if ( type == "standard" && NF == entry_count && entry_count > 0 && NR > 57 && NR < (59+N) ){
        printf("%d 1 %.16f %.16f %.16f %.8f %.8f\n",count+1,$3*H11,$4*H22,$5*H33, 1./sqrt(47.084529), 1.0);count++;
    }  
   
   
}
  END{}



另外,有谁知道这个小程序有没有哪个版本可以使得不同尺度的初始样品voroinize后的晶粒具相同的形状么?即可以保存和读入seed位置
6楼2012-08-06 20:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gswyyyx

铜虫 (小有名气)

【答案】应助回帖

引用回帖:
6楼: Originally posted by chen_kg at 2012-08-06 20:38:08
voronoirize 输入结构文件 晶体结构 x方向晶粒数目 y方向晶粒数目 z方向晶粒数目 输出结构文件
对应的FCC10x10x10.cfg是输入文件,bcc是结构 2 2 2表示将初始结构随机切成2*2*2个晶粒, FCC10x10x10-voronoirized. ...

请问有没有把cfg文件或者lammps文件输入到MS格式的代码?你这个代码是哪儿下载的啊?
7楼2013-08-27 18:34:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啾啾哈喽

新虫 (小有名气)

引用回帖:
6楼: Originally posted by chen_kg at 2012-08-06 20:38:08
voronoirize 输入结构文件 晶体结构 x方向晶粒数目 y方向晶粒数目 z方向晶粒数目 输出结构文件
对应的FCC10x10x10.cfg是输入文件,bcc是结构 2 2 2表示将初始结构随机切成2*2*2个晶粒, FCC10x10x10-voronoirized. ...

请问 这个代码怎么用啊
8楼2015-11-03 12:15:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hello806

新虫 (小有名气)

引用回帖:
8楼: Originally posted by 啾啾哈喽 at 2015-11-03 12:15:35
请问 这个代码怎么用啊...

找到这个代码怎么用的了吗
9楼2016-01-30 07:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1073269969

银虫 (小有名气)

引用回帖:
6楼: Originally posted by chen_kg at 2012-08-06 20:38:08
voronoirize 输入结构文件 晶体结构 x方向晶粒数目 y方向晶粒数目 z方向晶粒数目 输出结构文件
对应的FCC10x10x10.cfg是输入文件,bcc是结构 2 2 2表示将初始结构随机切成2*2*2个晶粒, FCC10x10x10-voronoirized. ...

我觉得您说的不对,您看它运行完的说明,他是给他分成了16个晶粒 我试过1 1 1 这样是分成了两个晶粒
10楼2018-01-15 09:38:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaorongwei 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +4 ziyangfang 2026-08-25 4/200 2026-08-25 14:41 by ZJTJZ
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +6 yuleib84 2026-08-25 7/350 2026-08-25 14:40 by 六两废铜
[基金申请] 什么时候开奖? +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见