| 10 | 1/1 | 返回列表 |
| 查看: 3361 | 回復: 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”怎么用, 哪位大俠知道啊,跪求。! |
分子動力學第一性原理 | 材料 |
金蟲 (小有名氣)
金蟲 (初入文壇)
金蟲 (小有名氣)
|
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ā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 材料學碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 【考研調(diào)劑】化學專業(yè) 281分,一志愿四川大學,誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 化學求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|