




Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Prepara tus exámenes
Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Prepara tus exámenes con los documentos que comparten otros estudiantes como tú en Docsity
Encuentra los documentos específicos para los exámenes de tu universidad
Estudia con lecciones y exámenes resueltos basados en los programas académicos de las mejores universidades
Responde a preguntas de exámenes reales y pon a prueba tu preparación
Consigue puntos base para descargar
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Comunidad
Pide ayuda a la comunidad y resuelve tus dudas de estudio
Ebooks gratuitos
Descarga nuestras guías gratuitas sobre técnicas de estudio, métodos para controlar la ansiedad y consejos para la tesis preparadas por los tutores de Docsity
efecto corona en matlab en donde se crea el codigo y despues se procede a simular
Tipo: Apuntes
1 / 8
Esta página no es visible en la vista previa
¡No te pierdas las partes importantes!





mt=input('cual es el coeficiente meteorologico:'); b=input('cual es el radio del conductor en metros:'); r=b100; d=input('cual es el factor de correccion de la densidad del aire:'); %calculo de Vc variando DMG y r% disp('calculo de VC') disp('opciones') disp('1. con radio del conductor constante y DMG variable') disp('2. con DMG constante y radio del conductor variable') a=input('cual es su opcion:'); if a== Vc=84mcdmtrlog10(DMG/r); k=0; Vcc=[]; dm=[]; for DMG=DMG+100:100:DMG+ Vc=84mcdmtrlog10(DMG/r); k=k+1; Vcc(1,k)=Vc; dm(1,k)=DMG; end plot(dm,Vcc) title('Variación del Vc con DMG variable a radio constate'); xlabel('DMG(cm)'); ylabel('Vcc (kV)'); grid on; end if a== Vc=84mcdmtrlog10(DMG/r); m=0; Vcc1=[]; r1=[]; for r=r+0.1:0.1:r+ Vc=84mcdmtrlog10(DMG/r); m=m+1; Vcc1(1,m)=Vc; r1(1,m)=r; end plot(r1,Vcc1) title('Voltaje critico disruptivo con r variable y DMG constante'); xlabel('r(cm)'); ylabel('Vcc (kV)'); grid on; end if Vc < Vmax disp('Existe posibilidad de efecto corona'); disp ('su Vc es de:'); disp (Vc); disp('su Vmax es de:'); disp (Vmax); disp('y el cálculo de perdidas seria en kv/km'); P=(241/d)(f+25)(r/DMG)^(1/3)((Vmax-Vc)/1.732050808)^2*0. end if Vc > Vmax
disp('no existe posibilidad de efecto corona'); disp ('su Vc es de:'); disp (Vc); disp('su Vmax es de:'); disp (Vmax); disp('Por lo tanto no se realizan cálculos de perdidas'); end
Vcc(1,k)=Vc; dm(1,k)=DMG; end plot(dm,Vcc) title('Variación del Vc con DMG variable a radio constate'); xlabel('DMG(cm)'); ylabel('Vcc (kV)'); grid on; end if a== Vc=84mcdmtrlog10(DMG/r); m=0; Vcc1=[]; r1=[]; for r=r+0.1:0.1:r+ Vc=84mcdmtrlog10(DMG/r); m=m+1; Vcc1(1,m)=Vc; r1(1,m)=r; end plot(r1,Vcc1) title('Voltaje critico disruptivo con r variable y DMG constante'); xlabel('r(cm)'); ylabel('Vcc (kV)'); grid on; end if Vc < Vmax disp('Existe posibilidad de efecto corona'); disp ('su Vc es de:'); disp (Vc); disp('su Vmax es de:'); disp (Vmax); disp('y el cálculo de perdidas seria en kvkm'); P=(241/d)(f+25)(r/DMG)^(1/3)((Vmax-Vc)/1.732050808)^20. end if Vc > Vmax disp('no existe posibilidad de efecto corona'); disp ('su Vc es de:'); disp (Vc); disp('su Vmax es de:'); disp (Vmax); disp('Por lo tanto no se realizan cálculos de perdidas'); end disp('--------- el calculo para varias fases seria el siguiente:-------------'); n=input('el número de conductores por fase es:'); R=input('el radio del haz de conductores es de:'); h1=input('la altura de fase a en metros es de:'); h2=input('la altura de fase b en metros es de:'); h3=input('la altura de fase c en metros es de:'); Flecha=input('el valor de la flecha en metros es de:'); HMG=((h1((2/3)Flecha))(h2((2/3)Flecha))(h3((2/3)Flecha)))^(1/3); RMG=((r2.718281828^(-1/4))DabDca(r2.718281828^(- 1/4))DbcDab(r2.718281828^(-1/4))DcaDbc)^(1/9); disp('el voltaje critico disruptivo es el siguiente:');
Vcrit=84d^(2/3)mtmcnr(1-(0.07r))(1-((n-1)r)/R)log10(((DMG/ RMG)2HMG)/((4*HMG^2)+DMG^2)^(1/2)); disp(Vcrit);