| 10 | 1/1 | 返回列表 |
| 查看: 3363 | 回復(fù): 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 現(xiàn)在我有FCC10x10x10.cfg 文件,問題是我不知道這句話“% voronoirize FCC10x10x10.cfg bcc 2 2 2 FCC10x10x10-voronoirized.cfg”怎么用, 哪位大俠知道啊,跪求。。 |
分子動力學(xué)第一性原理 | 材料 |
金蟲 (小有名氣)
金蟲 (初入文壇)
金蟲 (小有名氣)
|
voronoirize 輸入結(jié)構(gòu)文件 晶體結(jié)構(gòu) x方向晶粒數(shù)目 y方向晶粒數(shù)目 z方向晶粒數(shù)目 輸出結(jié)構(gòu)文件 對應(yīng)的FCC10x10x10.cfg是輸入文件,bcc是結(jié)構(gòu) 2 2 2表示將初始結(jié)構(gòu)隨機切成2*2*2個晶粒, FCC10x10x10-voronoirized.cfg是輸出結(jié)構(gòu)文件,下面的代碼可以將這個cfg直接轉(zhuǎn)成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位置 |
銅蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
銀蟲 (小有名氣)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 17/850 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|