Problemas ciencia de materiales en matlab

Universidad Autónoma de Ciudad Juárez
Instituto de Ingeniería y Tecnología

“Difusión”

Materiales para el diseño

1 de Marzo de 2011

5.12
Sedice que se forma un par de difusión cuando se permite que dos mételas de difundan uno en otro a elevada temperatura. En el caso de un bloque de metal puro A adyacente a un bloque de metal puro B, elperfil de concentración de A (en porcentaje atómico) tras la interdifusión esta dado por:
cx=50 1-erfx2Dt
Para un par difusion con D=10-14m2s, dibujese el perfil de concentración del metal A en unintervalo de 20 µm a ambos lados de la entrecara original (x=0) después de una hora.
%José Alejandro Sierra Gómez
%28 de febrero de 2011
%Problema 5.12
%PERFIL DE CONCENTRACION PARA METAL AADYACENTE A UN BLOQUE DE METAL PURO B,
%D=10E-14M^2/S, INTERVALO 20 uM A AMBOS LADOS DE LA ENTRECARA ORIGINAL
%A 1 HR

clc;
clf;
clear;

x = linspace(-20e-6, 20e-6, 40)';
D = 10e-14;t = 3600;
c= 50 * (1 - erf(x/(2*sqrt(D*t))));

plot(x,c, 'b')
axis([x(1) x(size(x,1)) 0 100]);
grid
xlabel('x, m')
ylabel('% concentracion')
title('Función de error:concentración en 1 hora')
hold on

5.13
-Utilícese la información del problema 5.12 para representar como progresa la interdifusión de dos metales, X e Y, cuando D=10-12m2s. Dibújese el perfil deconcentración del metal X en un intervalo de 300 µm a ambos lados de la entrecara inicial, después de 1, 2 y 3 horas.

%José Alejandro Sierra Gómez
%28 de febrero de 2011
%Problema 5.13
%PERFIL DECONCENTRACION DE UN METAL X, D=10E-12M^2/S, INTERVALO 20uM, A
%1,2,3 HORAS
clc;
clf;
clear;

x = linspace(-300e-6, 300e-6, 600)';
D = 10e-12; %coeficiente dedifusión

t(1)= 3600;
t(2)= 7200;
t(3) = 10800; % tiempo en segundos

for w=1:3
c=50 * (1 - erf(x/(2*sqrt(D*t(w)))))

figure(w)
plot(x,c, 'r')
axis([x(1)...