计算材料学2总结.ppt

下载文档 降价啦

35
0
约1.11万字
约 55页
2016-12-04 发布于湖北
举报
版权申诉
保障服务

计算材料学2总结.ppt

1、本文档共55页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

　如图所示，第i个分子在经过N步随机行走后距原点距离为R，对n个分子每步的位移平方求和后取平均值就得到了所有分子距原点的方均距离R2：　　！Monte Carlo Simulation of One Dimensional Diffusion INTEGER X,XX(1:1000,1:1000) REAL XXM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XX(J,I)：X*X ,J:第n个原子,I:第几步跳跃 ! XXM(I): THE MEAN OF XX WRITE(*,*) “原子个数JMAX，跳动次数 IMAX READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第n个原子跳跃 X=0 DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN0.5)THEN X=X+1 ELSE X=X-1 END IF XX(J,I)=X*X END DO END DO OPEN(1,FILE=“f:\DIF1.DAT) DO I=1,IMAX XXM=0.0 XXM(I)=1.0*SUM(XX(1:JMAX,I))/JMAX !! WRITE(1,*) I, XXM(I) END DO CLOSE(1) END ! Monte Carlo Simulation of Two Dimensional Diffusion INTEGER X,Y,XY(1:1000,1:1000) REAL XYM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XY(J,I)：X*Y ,J:第几个原子,I:第几步跳跃 ! XYM(I): THE MEAN OF XY WRITE(*,*) “原子个数JMAX，跳跃次数 IMAX READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几个原子 X=0 !!! Y=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN.LT.0.25)THEN x=x y=y-1 END IF IF(RN.LT.0.5.AND.RN.GE.0.25)THEN x=x y=y+1 END IF IF(RN.LT.0.75.AND.RN.GE.0.5)THEN x=x-1 y=y END IF IF(RN.GE.0.75)THEN x=x+1 y=y END IF XY(J,I)=X*X+Y*Y END DO END DO OPEN(1,FILE=“f:\DIF2.DAT) DO I=1,IMAX XYM=0.0 XYM(I)=1.0*SUM(XY(1:JMAX,I))/JMAX !! WRITE(1,*) I, XYM(I) END DO CLOSE(1) END ! Monte Carlo Simulation of Two Dimensional Diffusion INTEGER X,XY(1:1000,1:1000),y,XN(1:4),YN(1:4),RN REAL XYM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XY(J,I)：X*Y ,J:第n个实验,I:第几步跳跃 ! XYM(I): THE MEAN OF XY WRITE(*,*) “原子个数JMAX，跳动次数 IMAX READ(*,*) JMAX,IMAX XN=(/0,0,-1,1/) YN=(/-1,1,0,0/) ISEED=RTC() DO J=1,JMAX !第n个原子 X=0 !!! Y=0 !!! DO I=1,IMAX !第几步跳跃 RN=4*RAN(ISEED)+1 X=X+XN(RN) Y=Y+YN(RN) XY(J,I)=X*X+Y*Y END DO END DO OPEN(1,FILE=C:\DIF2.DAT) DO I=1,IMAX XYM=0.0 XYM(I)=1.0*SUM