Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clc
- open=input('preiona enter')
- //Variables de entrada//
- fp=0.95
- fprad=acos(fp)
- S=10 //Potencia aparente del sistema//
- V=500 //Voltaje de la línea//
- f=50 //Frecuencia del sistema//
- l=400 //Largo de la línea//
- cir=1 //Número de circuitos//
- seccion=765 //Sección del conductor en mm2//
- n=4 //Número de subconductores//
- G=0 //Conductancia de la línea//
- d=0.45 //Distancia entre conductores//
- D12=14 //Distancia entre conductores//
- D23=14 //Distancia entre conductores//
- D13=28 //Distancia entre conductores//
- RL=(0.0402/n) //Resistencia de la línea por subconductor//
- V2=complex(0.875,0) //Voltaje en pu del consumo//
- I2=complex(1.142857*cos(fprad),1.142857*sin(-fprad)) //Corriente en pu del consumo//
- r=(sqrt(seccion/%pi))*(0.001) //Radio del conductor//
- W=(2*%pi*f) //Omega//
- E0=8.8542*10^(-12) //Epsilon en vacío//
- XT1=complex(0.00266666,0.007999999)
- XT2=complex(0.03,0.12)
- //Condición dependiente del n //
- if n==1 then reql=(%e^-0.25)*r
- elseif n==2 then reql=sqrt((%e^-0.25)*r*d)
- elseif n==3 then reql=((%e^-0.25)*r*(d^2))^(1/3)
- elseif n==4 then reql=1.09*((((%e)^-0.25)*r*(d^3))^(1/4))
- else reql=r
- end
- if n==1 then reqc=r
- elseif n==2 then reqc=sqrt(r*d)
- elseif n==3 then reqc= (r*(d^2))^(1/3)
- elseif n==4 then reqc=1.09*((r*(d^3))^(1/4))
- else reqc=r
- end
- Deq=(D12*D13*D23)^(1/3) //Distancia media entre conductores//
- Zb=(V*V)/S //Cálculo de la impedancia base de la línea//
- //Ecuaciones de inductancia//
- L=(2*10^(-7))*(log(Deq/reql))*1000 //Inductancia en H/KM//
- XL=W*L //Reactancia en ohm/KM//
- Z=complex(RL,XL) //Impedancia de la línea//
- //Ecuaciones de conductancia//
- c=((2*%pi*E0)/(log(Deq/reqc)))*1000//Capacitancia en F/KM//
- B=W*c //Admitancia en S/KM//
- Y=complex(G,B) //Conductancia de la línea//
- //Ecuaciones para la matriz//
- Gamma=sqrt(Z*Y) //Constante de propagación//
- Zc=sqrt(Z/Y) //ni idea qué es//
- A=cosh(Gamma*l)
- B=((sinh(Gamma*l)*Zc)/Zb)
- C=((sinh(Gamma*l))/Zc)*Zb
- D=A
- T=[A,B;C,D]
- F=1/sqrt(real(A)^2+(imag(A)^2)) //Efecto ferranti en la línea//
- //Calculo de la matriz de entrada//
- MT1=[1,XT1;0,1]
- MT2=[1,XT2;0,1]
- MC=[V2;I2]
- MTG=MT1*T*MT2*MC
- V1=MTG(1.1)
- I1=MTG(2.1)
- //Cálculo de la potencia del generador//
- phi1=(atan((imag(V1))/(real(V1)))) //Ángulo phi del voltaje 1//
- phi2=(atan((imag(I1))/(real(I1)))) //Ángulo phi de la corriente 1//
- phi12=phi1-phi2 //Diferencia de ángulos entre V1 e I1//
- fp1=cos(phi12) //Factor de potencia del generador//
- I1r=real(I1) //Valor real de la corriente 1//
- I1i=imag(I1) //Valor imaginario de la corriente 1//
- P1=real(V1*(complex(I1r,-I1i)))
- //Cálculos de la potencia de salida//
- phi3=atan(imag(V2)/real(V2)) //Ángulo phi del voltaje 2//
- phi4=atan(imag(I2)/real(I2)) //Ángulo phi de la corriente 2//
- phi34=phi3-phi4 //Diferencia de ángulos entre V2 e I2//
- fp2=cos(phi34) //Factor de potencia del consumo//
- I2r=real(I2) //Valor real de la corriente 2//
- I2i=imag(I2) //Valor imaginario de la corriente //
- P2=real(V2*(complex(I2r,-I2i))) //Potencia real del consumo//
- Rend=P2/P1 //Rendimiento del sistema//
- Reg=(((sqrt(imag(V1)^2+real(V1)^2))-(sqrt(imag(V2)^2+real(V2)^2)))/(sqrt((imag(V1)^2)+(real(V1))^2)))*100
- S1=V1*I1
- disp(I1,'La corriente I1 es:')
- disp(V1,'El voltaje V1 es:')
- disp(Rend,'El rendimiento del sistema es:')
- disp(Reg,'La regulación del sistema es:')
- disp(S1,'La potencia aparente es:')
- disp (fp1,'factor de potencia')
- disp(r,'el radio')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement