重大计算机在材料工程中的应用课件03相图的计算机分析与计算

第3章相图的计算机分析与计算重庆大学材料科学与工程学院

汤爱涛3.1概述CALPHAD计算机在相图研究中应用的主要相图计算与绘制相图库系统计算机相图解析计算3.2相图计算基本方法热力学模型基本概念理想溶液模型规则溶液模型二元完全固溶相图的计算3.3相图计算软件简介CALPHAD技术Thermo-Calc(瑞典)F*A*C*T(加拿大)Lukas(德国)MTDATA(美国)上机实验5

相图的模拟计算实践实验目的了解相图计算基本原理了解常用相图计算软件的特点掌握采用理想溶液模型计算绘制二元相图的方法实验内容NiO-MgO二元相图的编程上机计算绘制实验要求:实验报告要求:NiO-MgO二元相图计算程序#include<stdio.h>#include<conio.h>#include<dos.h>#include<math.h>#include<graphics.h>#include<stdlib.h>#include<stdarg.h>externint

directvideo=0;

intTA=1960+273,TB=2800+273;

intn;floatR=8.310,DHA=52300,DHB=77400;floatdg(float

dh,int

t,int

tb){if(TA<=t&&t<=TB)return(dh*(tb-t)/tb);elseprintf("Thetemperatureisoutofrange\n");return0;}floatxb2(intk)

{floatdgb,dga,xbl;

dgb=dg(DHB,k,TB);

dga=dg(DHA,k,TA);

xbl=((exp(-dga/(R*k))-1)*exp(-dgb/(R*k)))/(exp(-dga/(R*k))-exp(-dgb/(R*k)));

returnxbl;

}

floatxb1(intk)

{floatdgb,dga,xbs;

dgb=dg(DHB,k,TB);

xbs=xb1(k)*exp(-dgb/(R*k));

returnxbs;

}initgrapher(){int

gdriver=DETECT,graphmode,errorcode;initgraph(&gdriver,&graphmode,"D:\\TC30\\bgi");

errorcode=graphresult();if(errorcode!=grOk){printf("graphics

error:%s\n",grapherrormsg(errorcode));

printf("pressanykeyhalt:");

getch();exit(1);}return0;}makemap1(){initgrapher();setviewport(50,50,getmaxx()-50,getmaxy()-50,1);

textcolor(BLUE);settextstyle(1,0,4);

outtextxy(70,100,"softwarefor");outtextxy(60,140,"twophasediagram\n");delay(5000);

closegraph();return0;}rectangle(X1,Y1,X2,Y2);setlinestyle(0,1,1);for(i=1;i<=9;i++){line(X1+i*xstep,Y1,X1+i*xstep,Y2);line(X1+i*xstep,Y1,X1+i*xstep,Y2);}makemap(float*solius,float*liquidus,int*p){intX1,X2,Y1,Y2,i,u,v,w,xstep,ystep;charmy[5],mx[5];

initgrapher();X1=getmaxx()/2-230;X2=getmaxx()/2+220;Y1=getmaxy()/2-150;Y2=getmaxy()/2+150;

xstep=(X2-X1)/10;ystep=(Y2-Y1)/10;

setbkcolor(WHITE);setcolor(1);setlinestyle(0,1,3);

for(i=1;i<=n+1;i++){p[i]=p[i]-273;u=((p[i]-1900)/1000.0*(Y2-Y1));v=(int)(solius[i]*(X2-X1));w=(int)(liquidus[i]*(X2-X1));putpixel(X1+w,Y2-u,BLUE);putpixel(X1+v,Y2-u,BLUE);}

for(i=1;i<=6;i++){gcvt(3.0-2.0*i,2,my);gcvt(0.2*i,2,mx);outtextxy(X1-30,Y1+2*i*ystep-40,my);outtextxy(X1+2*i*xstep-10,Y2+10,mx);}

settextstyle(1,1,1);outtextxy(30,140,"temperatue(1000)");settextstyle(1,0,2);outtextxy(200,50,"twophasediagram");outtextxy(190,