Optimización de un dispositivo de ultrasonidos para monitorizar el fraguado de un cemento

Memoria Proyecto de Fin de Carrera APÉNDICES Apéndice 1: Programa «load_bcp.m» pág. 60 Apéndice 2: Programa «cement_analyse.m» (Método-1:FFT) pág. 62 Apéndice 3: Programa «make_plots.m» pág. 66 Apéndice 4: Programas de soporte-I pág. 69 Apéndice 5: Programas de soporte-II pág. 70 Apéndice 6: Programas de soporte-III pág. 71 Apéndice 7: Programa con aire–agua pág. 72 Apéndice 8: Nuevo programa «load_bcp.m» pág. 73 Apéndice 9: Programa «cement_analyse.m» (Método-2:FFT) pág. 75 Apéndice 10: Impedancia acústica del CSC versus L/P pág. 77 Apéndice 11: Coeficiente de reflexión del CSC versus L/P pág. 80 Apéndice 12: Programa para calcular la velocidad del sonido pág. 83 Apéndice 13: Pulsos registrados en los ensayos con CFC pág. 90 Apéndice 14: Funciones Ac(t), Rp,c(t) y Zc(t) en los ensayos con CFC pág. 98 PFC Lorenzo González Pág. 59 59 PFC Lorenzo González Pág. 60 60 Apéndices Apéndice 1: Programa «load_bcp.m» Function out = load_bcp(filename) % out = load_bcp(filename); % ======================================================================== % 'filename'≡'nombre de fichero'. % 'bcp'≡'Programa de adquisición de datos'; genera el fichero 'desc.txt'. % 'desc.txt'≡'contiene condiciones de contorno' y ficheros '*.bin'. % '*.bin'≡'cada fichero contiene un registro de los pulsos P1 y P5'. % ======================================================================== % out = load_bcp(filename); % ======================================================================== % 'out=load_bcp(filename)'≡'matriz de salida del experimento'. % 'Columna No. 1'≡'t1:Timescale of pulses (in microsec)'. % 't1'≡'registro temporal generado a partir de 'Ts:Sampling time % (in nanosec)' del fichero 'desc.txt'. % 'Columna No. i'≡'i-fichero '*.bin'' según orden de registro de 't2:Time % scale (in sec)'. % ======================================================================== % Created: November 23, 1998, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: November 23, 1998, by Dr. J. Carlson (EISLAB) % ======================================================================== % Esta parte busca el fichero. Es una parte de control. % ======================================================================== if nargin < 1, disp('No filename specified, aborting'); return; end; if ~exist([filename, '_desc.txt']), disp(['File not found: ', filename]); return; end; % ======================================================================== % Esta parte lee la información de 'desc.txt' y la guarda en variables. % ======================================================================== desc = fopen([filename '_desc.txt'],'r'); fgetl(desc); out.Ts = str2num(fgetl(desc)); fgetl(desc); blocksize = str2num(fgetl(desc)); fgetl(desc); n_blocks = str2num(fgetl(desc)); fgetl(desc); Tstart = str2num(fgetl(desc)); fgetl(desc); ChLev = str2num(fgetl(desc)); t1 = (1:blocksize)'*out.Ts*1e-3+Tstart; fgetl(desc); fgetl(desc); fgetl(desc); noMeas = str2num(fgetl(desc)); fgetl(desc); dt2 = str2num(fgetl(desc)); PFC Lorenzo González Pág. 61 61 out.t2 = (0:noMeas-1)*dt2; fclose(desc); % ======================================================================== % Esta parte lee los ficheros '*.bin' y genera la matriz de los pulsos. % ======================================================================== P = zeros(blocksize,n_blocks*noMeas); for n = 1:noMeas, fid = fopen([filename, '_',num2str(n),'.bin'],'r'); % ======================================================================== % Esta parte filtra de ceros la señal de los pulsos. % ======================================================================== for k = 1:n_blocks, tmp = fread(fid,blocksize*2,'int8'); if length(find(tmp(1:2:end)==0)) < length(find(tmp(2:2:end)==0)), P(:,n*k) = tmp(1:2:end); else P(:,n*k) = tmp(2:2:end); end; end; fclose(fid); end; % ======================================================================== % Esta parte escala los pulsos según el valor del 'Channel Level' del OSC. % ======================================================================== out.pulses = detrend(P(4:end,:),0)*ChLev/25; out.t1 = t1(4:end); % ======================================================================== PFC Lorenzo González Pág. 62 62 Apéndice 2: Programa «cement_analyse.m» (Método-1:FFT) VERSION No. 11 % out=cement_analyse(pw,pc) % ===================================================================== % Input1≡out=load_bcp(MC) % Input2≡out=load_bcp(ME) % This program calculates Aw(t), Ac(t), Rp,w(t), Rp,c(t), zw(t) and zc(t). % ===================================================================== % Created: September 16, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: September 16, 2003, by Dr. J. Carlson (EISLAB) % ===================================================================== % ===================================================================== % Todos los cálculos se realizan en el campo de frecuencias. % Por este motivo, se aplica la FFT a los pulsos (ficheros '*.bin'). % ===================================================================== % Esta parte especifica variables de cálculo. % ===================================================================== dx = 30e-3; % Espesor del PMMA (MC) pw.temp = mean(pw.temp); % Temperatura de calibración del agua % ===================================================================== % Esta parte: 1) Ejecuta la FFT sobre el agua (MC). %2) Determina el valor máximo de la FFT. %3) Determina el índice de ese máximo. % ===================================================================== Nfft = 2^nextpow2(size(pc.pulses,1)); % Define el ancho de la FFT N = size(pc.pulses,2); % Calcula el número de pulsos del ME Aw = abs(fft(pw.pulses,Nfft)); % Calcula las Aiw de los i-pulsos [mx, maxidx] = max(Aw(1:(Nfft/2))); % Selecciona (Ajw,j)|Aiw:máximo:i=j Aw = Aw(maxidx).*ones(N,1); % Construye la matriz de los Ajw % ===================================================================== % Esta parte: 1) Repite los cálculos anteriores para el cemento (ME). % ===================================================================== Ap = abs(fft(pc.pulses,Nfft)); Ap = Ap(maxidx,:)'; % ===================================================================== % Esta parte: 1) Fija la temperatura del MC y del ME. % ===================================================================== pw.temp = 19; % Temperatura del agua pc.temp = 21; % Temperatura del cemento % ===================================================================== % Esta parte: 1) Calcula variables a las temperaturas especificadas. % ===================================================================== zw = z_water(pw.temp); % Calcula zw zp = z_PMMA(pc.temp); % Calcula zp Rpw = ((zw - zp) ./(zw + zp)); % Calcula Rp,w Rpc = Ap./Aw.* Rpw; % Calcula Rp,c zc = (1+Rpc)./(1-Rpc).*zp; % Calcula zc % ===================================================================== 1 Versión cedida por el Dr. J. Carlson (EISLAB-Sweden). Se modificó convenientemente. PFC Lorenzo González Pág. 63 63 % Esta parte: 1) Recoge todas las variables de salida. % ===================================================================== out.t = pc.t2/60; out.Rpc = Rpc; out.Rpw = Rpw; out.zc = zc; out.zp = zp; out.zw = zw; %out.dT = dT; out.Aw = Aw; out.Ac = Ap; % ===================================================================== VERSION No. 22 % out=cement_analyse(pw,pc) % ===================================================================== % Input1≡out=load_bcp(MC) % Input2≡out=load_bcp(ME) % This program calculates Aw(t), Ac(t), Rp,w(t), Rp,c(t), zw(t) and zc(t). % ===================================================================== % Created: September 16, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: February, 2004, by Dr. J. López López (GRICCA-UPC) % ===================================================================== % ===================================================================== % Todos los cálculos se realizan en el campo de frecuencias. % Por este motivo, se aplica la FFT a los pulsos (ficheros '*.bin'). % ===================================================================== % Esta parte especifica variables de cálculo. % ===================================================================== dx = 30e-3; % Espesor del PMMA (MC) pw.temp = mean(pw.temp); % Temperatura de calibración del agua % ===================================================================== % Esta parte: 1) Ejecuta la FFT sobre el agua (MC). %2) Determina el valor máximo de la FFT. %3) Determina el índice de ese máximo. % ===================================================================== Nfft = 2^nextpow2(size(cemento.pulses,1)); % Define el ancho de la FFT N = size(agua.pulses,2); % Calcula el número de pulsos del MC NN = size(cemento.pulses,2); % Calcula el número de pulsos del ME Aw = abs(fft(agua.pulses,Nfft)); % Calcula las Aiw de los i-pulsos [mx, maxidx] = max(Aw(1:(Nfft/2))); % Selecciona (Ajw,j)|Aiw:máximo:i=j Aw = Aw(maxidx).*ones(NN,1); % Construye la matriz de los Ajw % ===================================================================== % Esta parte: 1) Repite los cálculos anteriores para el cemento (ME). % ===================================================================== Ap = abs(fft(cemento.pulses,Nfft)); Ap = Ap(maxidx,:)'; % ===================================================================== % Esta parte: 1) Fija la temperatura del MC y del ME. % ===================================================================== 2 Se utilizó para obtener los resultados con el Método-1 de cálculo de la FFT en las muestras de CFC. PFC Lorenzo González Pág. 64 64 pw.temp = 19; % Temperatura del agua pc.temp = 21; % Temperatura del cemento % ===================================================================== % Esta parte: 1) Calcula variables a las temperaturas especificadas. % ===================================================================== zw = z_water(pw.temp); % Calcula zw zp = z_PMMA(pc.temp); % Calcula zp Rpw = ((zw - zp)./(zw + zp)); % Calcula Rp,w Rpc = Ap./Aw.* Rpw; % Calcula Rp,c alfa = Rpc.^2; % Calcula el parámetro alfa zc1 = (1 + alfa - 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 1ªSol. zc2 = (1 + alfa + 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 2ªSol. % ===================================================================== % Esta parte: 1) Recoge todas las variables de salida. % ===================================================================== out.t = cemento.t2/60; out.Rpc = Rpc; out.Rpw = Rpw; %out.zc = zc; out.zc1 = zc1; out.zc2 = zc2; out.zp = zp; out.zw = zw; %out.dT = dT; out.Aw = Aw; out.Ac = Ap; % ===================================================================== VERSION No. 33 % out=cement_analyse(pw,pc) % ===================================================================== % Input1≡out=load_bcp(MC) % Input2≡out=load_bcp(ME) % This program calculates Aw(t), Ac(t), Rp,w(t), Rp,c(t), zw(t) and zc(t). % ===================================================================== % Created: September 16, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: February, 2004, by Dr. J. López López (GRICCA-UPC) % ===================================================================== % ===================================================================== % Todos los cálculos se realizan en el campo de frecuencias. % Por este motivo, se aplica la FFT a los pulsos (ficheros '*.bin'). % ===================================================================== % Esta parte especifica variables de cálculo. % ===================================================================== dx = 30e-3; % Espesor del PMMA (MC) pw.temp = mean(pw.temp); % Temperatura de calibración del agua % ===================================================================== % Esta parte: 1) Ejecuta la FFT sobre el agua (MC). %2) Determina el valor máximo de la FFT. %3) Determina el índice de ese máximo. 3 Se utilizó para obtener los resultados con el Método-2 de cálculo de la FFT en las muestras de CFC. PFC Lorenzo González Pág. 65 65 % ===================================================================== Nfft = 2^nextpow2(size(cemento.pulses,1)); % Define el ancho de la FFT N = size(agua.pulses,2); % Calcula el número de pulsos del MC NN = size(cemento.pulses,2); % Calcula el número de pulsos del ME Aw = abs(fft(agua.pulses,Nfft/2)); % Calcula las Aiw de los i-pulsos for k=1:N [mx, maxidx] = max(Aw(1:(Nfft/4),k)); % Encuentra (Ajw,j)|Aiw:máximo:i=j maximo_ag(k) = mx; indice_ag(k) = maxidx; end medio_ag=mean(maximo_ag); % Calcula el valor medio de Ajw Aw = medio_ag.*ones(NN,1); % Construye la matriz de los Ajw % ===================================================================== % Esta parte: 1) Repite los cálculos anteriores para el cemento (ME). % ===================================================================== Ap = abs(fft(cemento.pulses,Nfft)); for n=1:NN [mx, maxidx] = max(Ap(1:(Nfft/2),n)); maximo(n) = mx; indice(n) = maxidx; end Ap = flipud(rot90(maximo)); % Pasa los valores de fila a columna % ===================================================================== % Esta parte: 1) Fija la temperatura del MC y del ME. % ===================================================================== pw.temp = 19; % Temperatura del agua pc.temp = 21; % Temperatura del cemento % ===================================================================== % Esta parte: 1) Calcula variables a las temperaturas especificadas. % ===================================================================== zw = z_water(pw.temp); % Calcula zw zp = z_PMMA(pc.temp); % Calcula zp Rpw = ((zw - zp)./(zw + zp)); % Calcula Rp,w Rpc = Ap./Aw.* Rpw; % Calcula Rp,c alfa = Rpc.^2; % Calcula el parámetro alfa zc1 = (1 + alfa - 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 1ªSol. zc2 = (1 + alfa + 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 2ªSol. % ===================================================================== % Esta parte: 1) Recoge todas las variables de salida. % ===================================================================== out.t = cemento.t2/60; out.Rpc = Rpc; out.Rpw = Rpw; %out.zc = zc; out.zc1 = zc1; out.zc2 = zc2; out.zp = zp; out.zw = zw; %out.dT = dT; out.Aw = Aw; out.Ac = Ap; % ===================================================================== PFC Lorenzo González Pág. 66 66 Apéndice 3: Programa «make_plots.m» VERSION No. 14 % make_plots(instruct) % ============================================================== % Plot acoustic impedance zc(t) and reflection coefficient % Rp,c(t) for cement sample. The data should be passed as an % instruct delivered by cement_analyze.m. % ============================================================== % Created: September 18, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: September 18, 2003, by Dr. J. Carlson (EISLAB) % ============================================================== % close all; set(0,'defaultaxesfontname','times'); set(0,'defaultaxesfontsize',12); figure (1); plot(instruct.t, instruct.zc*1e-6,'black'); grid on; xlabel('Setting time (minutes)'); ylabel('Acoustic impedance z_c (MPa\cdots/m)'); set(get(gca,'children'),'linewidth',2); figure (2); plot(instruct.t, instruct.Rpc,'black'); grid on; xlabel('Setting time (minutes)'); ylabel('Reflection coefficient, R_{p,c}'); set(get(gca,'children'),'linewidth',2); VERSION No. 25 % make_plots(instruct) % ============================================================== % Plot acoustic impedance zc(t), reflection coefficient Rp,c(t) % and Fourier amplitude Ac(t) for cement sample. The data % should be passed as an instruct delivered by cement_analyze.m. % ============================================================== % Created: September 18, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: February, 2004, by Dr. J. López (GRICCA-UPC) % ============================================================== % close all; set(0,'defaultaxesfontname','times'); set(0,'defaultaxesfontsize',12); figure (1); plot(resultados.t, resultados.zc1*1e-6,'black'); grid on; xlabel('Setting time (minutes)'); ylabel('Acoustic impedance z_c (MPa\cdots/m)'); set(get(gca,'children'),'linewidth',2); 4 Versión cedida por el Dr. J. Carlson (EISLAB-Sweden). Se modificó convenientemente. 5 Se utilizó para dibujar los resultados con el Método-1 de cálculo de la FFT en las muestras de CFC. PFC Lorenzo González Pág. 73 73 Apéndice 8: Nuevo programa «load_bcp.m» Function out = load_bcp(filename) % out = load_bcp(filename); % ======================================================================== % 'filename'≡'nombre de fichero'. % 'bcp'≡'Programa de adquisición de datos'; genera el fichero 'desc.txt'. % 'desc.txt'≡'contiene condiciones de contorno' y ficheros '*.bin'. % '*.bin'≡'cada fichero contiene un registro de los pulsos P1 y P5'. % ======================================================================== % out = load_bcp(filename); % ======================================================================== % 'out=load_bcp(filename)'≡'matriz de salida del experimento'. % 'Columna No. 1'≡'t1:Timescale of pulses (in microsec)'. % 't1'≡'registro temporal generado a partir de 'Ts:Sampling time % (in nanosec)' del fichero 'desc.txt'. % 'Columna No. i'≡'i-fichero '*.bin'' según orden de registro de 't2:Time % scale (in sec)'. % ======================================================================== % Created: November 23, 1998, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: February, 2004, by Dr. J. López (GRICCA-UPC) % ======================================================================== % Esta parte busca el fichero. Es una parte de control. % ======================================================================== if nargin < 1, disp('No filename specified, aborting'); return; end; if ~exist([filename, '_desc.txt']), disp(['File not found: ', filename]); return; end; % ======================================================================== % Esta parte lee la información de 'desc.txt' y la guarda en variables. % ======================================================================== desc = fopen([filename '_desc.txt'],'r'); fgetl(desc); out.Ts = str2num(fgetl(desc)); fgetl(desc); blocksize = str2num(fgetl(desc)); fgetl(desc); n_blocks = str2num(fgetl(desc)); fgetl(desc); Tstart = str2num(fgetl(desc)); fgetl(desc); ChLev = str2num(fgetl(desc)); t1 = (1:blocksize)'*out.Ts*1e-3+Tstart; fgetl(desc); fgetl(desc); fgetl(desc); noMeas = str2num(fgetl(desc)); fgetl(desc); dt2 = str2num(fgetl(desc)); out.t2 = (0:noMeas-1)*dt2; fclose(desc); PFC Lorenzo González Pág. 74 74 % ======================================================================== % Esta parte es nueva con respecto al programa load_bcp.m (Apéndice 1) % ======================================================================== Espacio=0.0083; Dt=0.000000002; v=zeros(noMeas,2); % ======================================================================== % Esta parte lee los ficheros '*.bin' y genera la matriz de los pulsos. % ======================================================================== P = zeros(blocksize,n_blocks*noMeas); for n = 1:noMeas, fid = fopen([filename, '_',num2str(n),'.bin'],'r'); % ======================================================================== % Esta parte filtra de ceros la señal de los pulsos. % ======================================================================== for k = 1:n_blocks, tmp = fread(fid,blocksize*2,'int8'); if length(find(tmp(1:2:end)==0)) < length(find(tmp(2:2:end)==0)), P(:,n*k) = tmp(1:2:end); else P(:,n*k) = tmp(2:2:end); end; end; fclose(fid); end; % ======================================================================== % Esta parte es nueva con respecto al programa load_bcp.m (Apéndice 1). % Esta parte permite calcular la velocidad del sonido en el ME. % ======================================================================== for n = 1:noMeas, fid = fopen([filename, '_',num2str(n),'.bin'],'r'); Tt = fread(fid,'int8'); L= size(Tt); [mx1, maxidx1] = max(Tt(L/100:(L/2))); [mx2, maxidx2] = max(Tt((L/2):L)); v(n,1) = Espacio/((maxidx2-maxidx1)*Dt)/1000; % Calcula velocidad en m/s v(n,2) = n/4; % Ahora se calcula el tiempo fclose(fid); end % ======================================================================== % Esta parte escala los pulsos según el valor del 'Channel Level' del OSC. % ======================================================================== out.pulses = detrend(P(4:end,:),0)*ChLev/25; out.t1 = t1(4:end); % ======================================================================== % Esta parte recoge la variable de salida 'velocidad'. % ===================================================================== out.v = v % ======================================================================== PFC Lorenzo González Pág. 75 75 Apéndice 9: Programa «cement_analyse.m» (Método-2:FFT) Function out = cement_analyze(pw,pc); % out=cement_analyse(pw,pc) % ===================================================================== % Input1≡out=load_bcp(MC) % Input2≡out=load_bcp(ME) % This program calculates Aw(t), Ac(t), Rp,w(t), Rp,c(t), zw(t) and zc(t). % ===================================================================== % Created: September 16, 2003, by Dr. J. Carlson (EISLAB-Sweden) % Last modified: February, 2004, by Dr. J. López López (GRICCA-UPC) % ===================================================================== % ===================================================================== % Todos los cálculos se realizan en el campo de frecuencias. % Por este motivo, se aplica la FFT a los pulsos (ficheros '*.bin'). % ===================================================================== % Esta parte especifica variables de cálculo. % ===================================================================== dx = 30e-3; % Espesor del PMMA (MC) pw.temp = mean(pw.temp); % Temperatura de calibración del agua % ===================================================================== % Esta parte: 1) Ejecuta la FFT sobre el agua (MC). %2) Determina el valor máximo de la FFT. %3) Determina el índice de ese máximo. % ===================================================================== Nfft = 2^nextpow2(size(cemento.pulses,1)); % Define el ancho de la FFT N = size(agua.pulses,2); % Calcula el número de pulsos del MC NN = size(cemento.pulses,2); % Calcula el número de pulsos del ME Aw = abs(fft(agua.pulses,Nfft/2)); % Calcula las Aiw de los i-pulsos for k=1:N [mx, maxidx] = max(Aw(1:(Nfft/4),k)); % Encuentra (Ajw,j)|Aiw:máximo:i=j maximo_ag(k) = mx; indice_ag(k) = maxidx; end medio_ag=mean(maximo_ag); % Calcula el valor medio de Ajw Aw = medio_ag.*ones(NN,1); % Construye la matriz de los Ajw % ===================================================================== % Esta parte: 1) Repite los cálculos anteriores para el cemento (ME). % ===================================================================== Ap = abs(fft(cemento.pulses,Nfft)); for n=1:NN [mx, maxidx] = max(Ap(1:(Nfft/2),n)); maximo(n) = mx; indice(n) = maxidx; end Ap = flipud(rot90(maximo)); % Pasa los valores de fila a columna % ===================================================================== % Esta parte: 1) Fija la temperatura del MC y del ME. % ===================================================================== pw.temp = 19; % Temperatura del agua PFC Lorenzo González Pág. 76 76 pc.temp = 21; % Temperatura del cemento % ===================================================================== % Esta parte: 1) Calcula variables a las temperaturas especificadas. % ===================================================================== zw = z_water(pw.temp); % Calcula zw zp = z_PMMA(pc.temp); % Calcula zp % zp = 3.2626*10^6; Rpw = ((zw - zp)./(zw + zp)); % Calcula Rp,w Rpc = Ap./Aw.* Rpw; % Calcula Rp,c alfa = Rpc.^2; % Calcula el parámetro alfa zc1 = (1 + alfa - 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 1ªSol. zc2 = (1 + alfa + 2 * sqrt(alfa))./(1 - alfa)*zp; % Calcula la 2ªSol. % ===================================================================== % Esta parte: 1) Recoge todas las variables de salida. % ===================================================================== out.t = cemento.t2/60; out.Rpc = Rpc; out.Rpw = Rpw; %out.zc = zc; out.zc1 = zc1; out.zc2 = zc2; out.zp = zp; out.zw = zw; %out.dT = dT; out.Aw = Aw; out.Ac = Ap; % ===================================================================== PFC Lorenzo González Pág. 77 77 Apéndice 10: Impedancia acústica del CSC versus L/P L/P = 3 ml/g L/P = 3 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) L/P = 3 ml/g L/P = 2.5 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) L/P = 2ml/g L/P = 1.5 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 78 78 L/P = 1.0ml/g L/P = 1.0 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) L/P = 0.8ml/g L/P = 0.8 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) L/P = 0.6ml/g L/P = 0.6 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) L/P = 0.4ml/g L/P = 0.4 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 79 79 L/P = 0.2ml/g L/P = 0.2 ml/g 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) polvo de CaSO4.½H2O polvo de CaSO4.½H2O 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) polvo de CaSO4.2H2O polvo de CaSO4.2H2O 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 0 10 20 30 40 50 60 0 0.5 1 1.5 2 2.5 3 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 80 80 Apéndice 11: Coeficiente de reflexión del CSC versus L/P L/P = 3 ml/g L/P = 3 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c L/P = 3 ml/g L/P = 2.5 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c L/P = 2ml/g L/P = 1.5 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c PFC Lorenzo González Pág. 81 81 L/P = 1.0ml/g L/P = 1.0 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c L/P = 0.8ml/g L/P = 0.8 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c L/P = 0.6ml/g L/P = 0.6 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c L/P = 0.4ml/g L/P = 0.4 ml/g 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c PFC Lorenzo González Pág. 82 82 L/P = 0.2ml/g L/P = 0.2 ml/g 0 10 20 30 40 50 60 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c polvo de CaSO4.½H2O polvo de CaSO4.½H2O 0 10 20 30 40 50 60 -1 -0.8 -0.6 -0.4 -0.2 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -1 -0.8 -0.6 -0.4 -0.2 0 Setting tim e (minutes) Reflection coefficien t , R p,c polvo de CaSO4.2H2O polvo de CaSO4.2H2O 0 10 20 30 40 50 60 -1 -0.8 -0.6 -0.4 -0.2 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 10 20 30 40 50 60 -1 -0.8 -0.6 -0.4 -0.2 0 Setting tim e (minutes) Reflection coefficien t , R p,c PFC Lorenzo González Pág. 89 89 if cc(k)>0; v(k) = 2* esp /(cc(k)*dt); k=k+1; end x2=x2+1; x1=x1+1; end vv = mean(v); else vv = 0; end end % ================================================================== % Variable de salida % ================================================================== out.vv = vv; % ================================================================== Apéndice 13: Pulsos registrados en los ensayos de CFC. M1-01 M2-01 M1-02 M2-02 M1-03 M2-03 PFC Lorenzo González Pág. 91 91 M1-04 M2-04 M1-05 (*) M2-05 M1-06 M2-06 (*) Se ha sustituido el primer pulso por el segundo, porque salía una singularidad. PFC Lorenzo González Pág. 92 92 M1-07 M2-07 M1-08 M2-08 M1-09 M2-09 PFC Lorenzo González Pág. 93 93 M1-10 (**) M2-10 M1-11 M2-11 M1-12 M2-12 (**) Se ha sustituido el pulso 432 por el 433, porque aparecía una singularidad. PFC Lorenzo González Pág. 94 94 M1-13 M2-13 M1-14 M2-14 (***) M1-15 M2-15 (***) Por error no se realiza correctamente el ensayo M2-14. PFC Lorenzo González Pág. 95 95 M1-16 M2-16 M1-17 M2-17 M1-18 M2-18 PFC Lorenzo González Pág. 96 96 M1-19 M2-19 M1-20 M2-20 M1-21 M2-21 PFC Lorenzo González Pág. 97 97 M1-22 M2-22 M1-23 M2-23 M1-24 M2-24 PFC Lorenzo González Pág. 98 Apéndice 14: Funciones A c (t), R p,c (t) y Z c (t) en los ensayos con CFC M1-01 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 105 M2-04 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 106 M1-05 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 107 M2-05 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 108 M1-06 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 109 M2-06 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 110 M1-07 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 111 M2-07 1º RESULTADO. 0 20 40 60 80 100 120 0 5 10 15 20 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 2.6 2.8 3 3.2 3.4 3.6 3.8 4 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 2º RESULTADO. 0 20 40 60 80 100 120 0 10 20 30 40 50 60 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 -1 0 1 2 3 4 5 6 7 8 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 112 M1-08 1º RESULTADO. 0 20 40 60 80 100 120 0 5 10 15 20 25 30 35 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 3.05 3.1 3.15 3.2 3.25 3.3 3.35 3.4 3.45 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 2º RESULTADO. 0 20 40 60 80 100 120 0 5 10 15 20 25 30 35 40 45 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 1 2 3 4 5 6 7 8 9 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 113 M2-08 1º RESULTADO. 0 20 40 60 80 100 120 0 5 10 15 20 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 2.9 3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) 2º RESULTADO. 0 20 40 60 80 100 120 0 5 10 15 20 25 30 35 40 45 Setting tim e (minutes) Fourier Transform coefficien t , Ac 0 20 40 60 80 100 120 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 Setting tim e (minutes) Reflection coefficien t , R p,c 0 20 40 60 80 100 120 0 5 10 15 20 25 30 Setting tim e (minutes) Acoustic im p edan ce z c (MPa ⋅ s/m) PFC Lorenzo González Pág. 114 M1-09 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 121 M2-12 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 122 M1-13 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 123 M2-13 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 124 M1-14 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 125 M1-15 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 126 M2-15 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 127 M1-16 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 128 M2-16 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 129 M1-17 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 130 M2-17 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 137 M1-21 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 138 M2-21 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 139 M1-22 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 140 M2-22 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 141 M1-23 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 142 M2-23 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 143 M1-24 1º RESULTADO. 2º RESULTADO. PFC Lorenzo González Pág. 144 M2-24 1º RESULTADO. 2º RESULTADO.