| 9 | 1/1 | 返回列表 |
| 查看: 1813 | 回復: 8 | |||
[求助]
用DL_POLY模擬了硅酸鹽熔體后,如何計算硅氧四面體Qn結構?
|
| 模擬了CaO,Al2O3,SiO2的體系,然后想計算其中硅氧四面體類型分布,即計算與所有硅配位的氧原子是否橋氧, |
動力學模擬 |
金蟲 (著名寫手)
|
這是統(tǒng)計程序,matlab的! 供你參考一下! clear Ob=0;Onb=0;Of=0; A=xlsread('CS1'); [m,n]=size(A); Si=0;Ca=0;O=0; for i=2:m if A(i,4)==1 Si=Si+1; elseif A(i,4)==2 Ca=Ca+1; else O=O+1; end end for i=Si+Ca+2:Si+Ca+O+1 %O的個數 for j=2:Si+Ca+1 %Si和Ca的個數 B(i-Si-Ca-1,j-1)=sqrt((A(i,1)-A(j,1))^2+(A(i,2)-A(j,2))^2+(A(i,3)-A(j,3))^2); end [P,Q]=sort(B(i-Si-Ca-1,1:Si+Ca)); if Q(1)<Si+2&&Q(2)<Si+2 Ob=Ob+1; end if Q(1)<Si+2&&Q(2)>Si+1 Onb=Onb+1; end if Q(1)>Si+1&&Q(2)<Si+2 Onb=Onb+1; end if Q(1)>Si+1&&Q(2)>Si+1 Of=Of+1; end end Ob=Ob/O*100 Onb=Onb/O*100 Of=Of/O*100 發(fā)自小木蟲IOS客戶端 |
金蟲 (著名寫手)
金蟲 (著名寫手)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 326求調劑 +5 | 諾貝爾化學獎覬?/a> 2026-03-15 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 材料學碩333求調劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 307求調劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 265求調劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 初始318分求調劑(有工作經驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 材料專碩英一數二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 295材料求調劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 289求調劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
[考研] 085410人工智能專碩317求調劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|