24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3380  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 086000生物与医药求调剂 +3 老天眷顾之人 2026-03-31 3/150 2026-04-05 22:24 by syh9288
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 301求调剂 +12 121. 2026-04-04 12/600 2026-04-05 09:00 by 来看流星雨10
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 08:15 by 544594351
[考研] 323分(计算机视觉和大模型项目)能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 一志愿南昌大学324求调剂 +9 hanamiko 2026-03-30 9/450 2026-04-04 11:04 by 猪会飞
[考研] 292求调剂 +21 是妍子也是研子 2026-03-30 22/1100 2026-04-03 21:44 by qlm5820
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 286求调剂 +7 Faune 2026-03-30 7/350 2026-04-03 10:14 by linyelide
[考研] 085602化工求调剂(331分) +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 一志愿中海洋材料357 +4 麦恩莉. 2026-03-30 4/200 2026-03-31 14:35 by 记事本2026
信息提示
请填处理意见