欢迎来到天天文库
浏览记录
ID:58846114
大小:145.51 KB
页数:7页
时间:2020-09-23
《渗碳的计算机模拟及程序.docx》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、20钢在碳含量为1.2%的渗碳气氛中渗碳情况程序:co=0.2;cp=1.2;d=0.162*exp((-)/(8.314*1203));x=0.0437*exp((-79953)/(8.314*1203));t=3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;t=4*3600;a=0:0
2、.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;t=8*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)h
3、oldon;t=20*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;co=0.2;cp=1.2;d=0.162*exp((-)/(8.314*1203));x=0.0473*2*exp((-79953)/(8.314*1203));t=3600;a=0:0.01:0.5;c=co+(cp
4、-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c,'-.')holdon;t=4*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c,'-.')holdon;
5、t=8*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c,'-.')holdon;t=20*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t
6、/d)^0.5));plot(a,c,'-.')holdon;其中虚线为界面传递系数增大一倍时的渗碳情况。从左至右依次是1、4、8、20小时的渗碳曲线。c/%x/cm碳浓度随渗碳层厚度的变化20钢在碳含量为1.2%的渗碳气氛中渗碳情况程序:co=0.2;cp=1.2;d=0.162*exp((-)/(8.314*1203));x=0.0437*exp((-79953)/(8.314*1203));t=3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0
7、.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;t=4*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;t=8*3600;a=0:0.01:0.5;c=co+(cp-co).*
8、(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5));plot(a,c)holdon;t=20*3600;a=0:0.01:0.5;c=co+(cp-co).*(erfc(a/(2.*(d.*t)^0.5))-exp((x.*a+x^2.*t)/d).*erfc(a/(2.*(d.*t)^0.5)+x.*(t/d)^0.5))
此文档下载收益归作者所有