Hola buenas noches. Si se puede utilizar para cualquier sistema en por unidad. En caso de que tenga elementos shunt, ver la modificación en las respuestas a Juan Sebastian Ladino Cruz. Saludos
Hola ! Muchas gracias por su tiempo para compartir su conocimiento. Me puede ayudar ya que en su código me marca un error en la línea 16. Por favor. Gracias. Saludos.
Buenas tardes, entiendo la sumatoria de la relación de los números imaginarios dependiendo de la relación que haya en el bus, pero lo que no tengo entendido es el cambio de los signos, me ayudaría mucho a tender mejor el tema
Tutoriales Ing, una herramienta muy util para el analisis S.E.P, mi consulta es sobre cuando tengo mi Xg del generador al momento de calcular la diagonal me sale un error me podria ayudar, gracias
ingeniero espero este bien, analizando el programa cuando tenemos un modelo de con Yshunt a ambos costados de la rama modelo tipo PI se presentan inconvenientes.
Hola Juan Sebastian. Efectivamente como tu mencionas el programa presentado en este video no considera el efecto capacitivo de la línea en ambos lados. Sin embargo, sólo requiere un pequeño cambio en el codigo agregando una columna adicional a la tabla de Líneas. Luego en el mismo ciclo for que recorre las líneas (calculo de la diagonal) le sumas este valor a la diagonal principal y ya esta. Espero que se haya aclarado tu duda. Saludos. El código Modificado es este: clc % limpiar pantalla clear % limpiar variables %% Cargar los datos % k m rkm xkm gm bm Yshunt/2 Lineas = [1 2 0 0.125 0 0 0 1 3 0 0.250 0 0 0.85 2 3 0 0.250 0 -0.80 0 1 4 0 0.400 0 0 0.65 2 4 0 0.200 0 -0.80 0]; B = max(max(Lineas(:,1:2))); Ybus = zeros(B,B); % Fuera de la diagonal y elementos shunt for i = 1:size(Lineas,1) k = Lineas(i,1); m = Lineas(i,2); zkm = Lineas(i,3) + 1j*Lineas(i,4); Ybus(k,m) = -1/zkm; Ybus(m,k) = Ybus(k,m); end % Elementos de la diagonal principal Ybus = Ybus - diag(sum(Ybus,2)); for i = 1:size(Lineas,1) k = Lineas(i,1); m = Lineas(i,2); gm = Lineas(i,5); bm = Lineas(i,6); ys = Lineas(i,7); Ybus(m,m) = gm + 1j*bm + Ybus(m,m) + 1j*ys; Ybus(k,k) = Ybus(k,k) + 1j*ys; end
@@tutorialesing.electrica7337 a listo ingeniero muchas gracias, lo voy a implementar y le comento como me fue, por cierto ingeniero nos conocemos de un intersemestral de física 3 en la utp ahora estoy viendo analisis de sistemas de potencia con el ingeniero granada.
muy buena explicación
me ayudo mucho a entender de mejor manera Stevenson ya que ahí lo explica en un lenguaje un tanto matemático.
éxitos
Muchas gracias, un trabajo muy didáctico.
Muchas gracias. Me alegra que le haya gustado.
Con mucho gusto. La idea es ir nutriendo el canal con ese tipo de material.
Buen dia el codigo se puede utiliar para cualquier ejercicio? gracias
Hola buenas noches. Si se puede utilizar para cualquier sistema en por unidad. En caso de que tenga elementos shunt, ver la modificación en las respuestas a Juan Sebastian Ladino Cruz. Saludos
no puedes dejar el código en algun archivo txt?? gracias
Hola! muy buen video, quería saber si es posible agregar en su codigo, una forma de calcular la corriente de cortocircuito?
Hola ! Muchas gracias por su tiempo para compartir su conocimiento. Me puede ayudar ya que en su código me marca un error en la línea 16. Por favor. Gracias. Saludos.
Buenas tardes, entiendo la sumatoria de la relación de los números imaginarios dependiendo de la relación que haya en el bus, pero lo que no tengo entendido es el cambio de los signos, me ayudaría mucho a tender mejor el tema
Tutoriales Ing, una herramienta muy util para el analisis S.E.P, mi consulta es sobre cuando tengo mi Xg del generador al momento de calcular la diagonal me sale un error me podria ayudar, gracias
Hola buenas noches. Si el generador tiene un elemento Xg. Se debe escribir este como un elemento shunt en el nodo donde esté el generador.
Buenas tarde muy buena la explicación como puedo obtener el libro
Me alegra que le haya gustado el material. (www.engbookspdf.com/download/Power-System/Power-System-Analysis-by-Grainger-and-Stevenson)
ingeniero espero este bien, analizando el programa cuando tenemos un modelo de con Yshunt a ambos costados de la rama modelo tipo PI se presentan inconvenientes.
Hola Juan Sebastian. Efectivamente como tu mencionas el programa presentado en este video no considera el efecto capacitivo de la línea en ambos lados. Sin embargo, sólo requiere un pequeño cambio en el codigo agregando una columna adicional a la tabla de Líneas. Luego en el mismo ciclo for que recorre las líneas (calculo de la diagonal) le sumas este valor a la diagonal principal y ya esta. Espero que se haya aclarado tu duda. Saludos. El código Modificado es este:
clc % limpiar pantalla
clear % limpiar variables
%% Cargar los datos
% k m rkm xkm gm bm Yshunt/2
Lineas = [1 2 0 0.125 0 0 0
1 3 0 0.250 0 0 0.85
2 3 0 0.250 0 -0.80 0
1 4 0 0.400 0 0 0.65
2 4 0 0.200 0 -0.80 0];
B = max(max(Lineas(:,1:2)));
Ybus = zeros(B,B);
% Fuera de la diagonal y elementos shunt
for i = 1:size(Lineas,1)
k = Lineas(i,1); m = Lineas(i,2);
zkm = Lineas(i,3) + 1j*Lineas(i,4);
Ybus(k,m) = -1/zkm; Ybus(m,k) = Ybus(k,m);
end
% Elementos de la diagonal principal
Ybus = Ybus - diag(sum(Ybus,2));
for i = 1:size(Lineas,1)
k = Lineas(i,1); m = Lineas(i,2);
gm = Lineas(i,5); bm = Lineas(i,6); ys = Lineas(i,7);
Ybus(m,m) = gm + 1j*bm + Ybus(m,m) + 1j*ys;
Ybus(k,k) = Ybus(k,k) + 1j*ys;
end
@@tutorialesing.electrica7337 a listo ingeniero muchas gracias, lo voy a implementar y le comento como me fue, por cierto ingeniero nos conocemos de un intersemestral de física 3 en la utp ahora estoy viendo analisis de sistemas de potencia con el ingeniero granada.