Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Karen Iter10

Descargar como xlsx, pdf o txt
Descargar como xlsx, pdf o txt
Está en la página 1de 32

Pagina Web de la asignatura:

http://umss-ie-sistemasdepotencia2.jimdo.com/

PRACTICA 1
FLUJO DE CARGA RESUELTO POR EL METODO GAUSS-SEIDEL

Determinar las tensiones de barra, considera un error maxima (Epsilon) igual a 0.0001

Datos de Barra

Barra Tension Nominal Tension (pu) Angulo Pgen (MW)


1 115 kV 1.02 0 ?
2 115 kV 1.05 ? 20
3 115 kV ? ? 0

Datos de Linea

Linea Z pu B100 Ysh/2 pu B100


LIN 1-2 0.008+j*0.08 0+j*0.0045
LIN 2-3 0.01+j*0.1 0+j*0.006

Datos de Demanda

P (MW) P (MW) P (MW) P (MW) P (MW) P (MW)


Barra grupo 1 grupo 2 grupo 3 grupo 4 grupo 5 grupo 6
1 10 15 20 25 30 35
2 60 55 50 45 40 35
3 20 20 20 20 20 20
Total 90 90 90 90 90 90

Factor de Potencia

Barra FP grupo 1 FP grupo 2 FP grupo 3 FP grupo 4 FP grupo 5 FP grupo 6


1 0.95 0.95 0.95 0.95 0.95 0.95
2 0.97 0.97 0.97 0.97 0.97 0.97
3 0.88 0.88 0.88 0.88 0.88 0.88

Grupos

Grupo 1 Grupo 2 Grupo 3

Jorge Ruiz Quispe Ricardo Toalina Cruz Jhussen Villegas


Froilan Luna Apaza Fernando Lizon Mamani RiveroJuan J. Martinez
Nilton Vargas Moller Diego Gomez Soruco Walter Santi
Jhimy Alanes Pena Marco A, Bautista Luna Rudy Romera
Torres T Wilder Romel Rojas Diaz Iver Cotari

Notlin32@hotmail.com odracir_tc@hotmail.com ivercotari.91@gmail.com


tecnodisco431@hotmail.com ro_mel23@hotmail.com j.j.mc.60026007@gmail.com
Wildr_fx@hotmail.com sal-vatore_009@hotmail.com
B2

B1 B3

P (MW) P (MW)
grupo 7 grupo 8
40 45
30 25
20 20
90 90

FP grupo 7 FP grupo 8
0.95 0.95
0.97 0.97
0.88 0.88

Grupo 4

Felix Araque Huarachi


Enrique Mendez Soto
Remberto Olivera Carballo
Richard Sanabria Alcocer
Grover Veizaga Sejas
David Ascui Huiza
e_mendez_631@hotmail.com
PRACTICA 1
FLUJO DE CARGA RESUELTO POR EL METODO GAUSS-SEIDEL

Determinar las tensiones de barra, considera un error maxima (Epsilon) igual a 0.0001

Datos de Barra

Barra Tension NominTension (pu) Angulo Pgen (MW) fi [rad] fi [deg]


1 115 kV 1,2 0 ? 0.3176 18.195
2 115 kV 1,05 ? 20 0.2456 14.070
3 115 kV ? ? 0 0.4949 28.358

Datos de Linea
CARGA
Linea Z pu B100 Ysh/2 pu B100 1 15+4,93026157768295i 0,15+0,0493026157768295i
LIN 1-2 0,008+0,08i 0,0045i 2 55+13,7842993394408i 0,55+0,137842993394408i
LIN 2-3 0,01+0,1i 0,006i 3 20+10,7948564427617i 0,2+0,107948564427617i

Datos de Demanda Factor de Potencia

P (MW)
Barra grupo 2 Q[MVAr] Barra FP grupo 2
1 15 4.930 1 0.95 alfa
2 55 13.784 2 0.97 1.20E+00
3 20 10.795 3 0.88
Total 90

1. Lectrura de datos
2. Conversión a por unidad
3. Formacion del YBUS

Ylin12= 1,23762376237624-12,3762376237624i
YLlin23= 0,99009900990099-9,9009900990099i

1,23762376237624-12,3717376237624i -1,23762376237624+12,3762376237624i

YBUS= -1,23762376237624+12,3762376237624i2,22772277227723-22,2667277227723i -0,99009900990099+9,9009

0 -0,99009900990099+9,9009900990099i 0,99009900990099-9,89499

4. Definir tipos de nodos

1 SL
2 PV
3 PQ
5. Calculo de las potencias especificadas

1 SL Sesp1= Pgen1+Qgen1-0,15-0,0493026157768295i
2 PV Sesp2= 0,2+Qgen2 -0,55-0,137842993394408i Sinyesp2= -0,35+Qiny2'
3 PQ Sesp3= 0 -0,2-0,107948564427617i

6. Definir el vector de estado

V1= 1,2
VE= V2= 1,05
V3= 1

7. PROCESO ITERATIVO 1

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,33910891089109-23,3800641089109i
Y23*V3= -0,99009900990099+9,9009900990099i
Y21*V1+Y22*V2+Y23*V3= -0,13613861386139+1,3724111386139i
V2*= 1,05
V2*{ }= -0,14294554455446+1,44103169554459i
Qiny2= 0,14294554455446-1,44103169554459i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,20705445544554-1,44103169554459i
Sesp2*= -0,20705445544554+1,44103169554459i
V2*= 1,05
Sesp2*/V2*= -0,197194719471943+1,3724111386139i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,33910891089109-23,3800641089109i
Y23*V3= -0,99009900990099+9,9009900990099i
[-Y21V1-Y22V2-Y23V3]= -2,47524752475248+24,7524752475248i
{-[ ]}= 2,47524752475248-24,7524752475248i
[ + ]= 2,27805280528054-23,3800641089109i
V2= 1,04972838608404-0,00271485895266421i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
1-alfa)*V2^k= -0.21
alfa*V2^k+1= 1,25967406330085-0,00325783074319705i
V2acel= 1,04967406330085-0,00325783074319705i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049679119
|V2esp|/|V2acel|= 1,00030569448043
V2correg= 1,04999494286825-0,00325882664407342i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1
Sesp3*/V3*= -0,0018348623853211+0,000990353802088229i
Y31xV1= 0
Y32xV2= -1,00733334299754+10,3992160943827i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 1,00733334299754-10,3992160943827i
[ + ]= 1,00549848061222-10,3982257405806i
V3= 1,05050765360405-0,00349753832679661i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)*V3^k -0.2
alfa*V3^k+1= 1,26060918432486-0,00419704599215593i
V3acel= 1,06060918432486-0,00419704599215593i
|V3acel|= 1.060617489

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -5,05713175003031E-06-0,00325882664407342i
(|V3^(k+1) - V3^k| ) ( 0,06060918432486-0,00419704599215593i

7. PROCESO ITERATIVO 2

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,29950495049505+12,9950495049505i
Y22*V2= 2,26653423942423-23,3872112654611i
Y23*V3= -1,00855319247852+10,5052365240007i
Y21*V1+Y22*V2+Y23*V3= -0,0415239035493402+0,113074763490101i
V2*= 1,04999494286825+0,00325882664407342i
V2*{ }= -0,04396837978699+0,118592610627377i
Qiny2= 0,04396837978699-0,118592610627377i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,30603162021301-0,118592610627377i
Sesp2*= -0,30603162021301+0,118592610627377i
V2*= 1,04999494286825+0,00325882664407342i
Sesp2*/V2*= -0,291106739974726+0,113849383599246i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,26653423942423-23,3872112654611i
Y23*V3= -1,00855319247852+10,5052365240007i
[-Y21V1-Y22V2-Y23V3]= -2,49370170733001+25,3567216725156i
{-[ ]}= 2,49370170733001-25,3567216725156i
[ + ]= 2,20259496735528-25,2428722889164i
V2= 1,13222244651204-0,0143569708132774i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,20999898857365+0,000651765328814684i
alfa*V2^k+1= 1,35866693581445-0,0172283649759329i
V2acel= 1,1486679472408-0,0165765996471182i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.148787551
|V2esp|/|V2acel|= 0,914007119373709
V2correg= 1,04989068157448-0,0151511300924737i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,06060918432486+0,00419704599215593i
Sesp3*/V3*= -0,188165163378337+0,102524383042804i
Y31xV1= 0
Y32xV2= -0,889484535296775+10,4099583622151i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,889484535296775-10,4099583622151i
[ + ]= 0,701319371918438-10,3074339791723i
V3= 1,03837763674727-0,0330245199700664i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k-0,212121836864972+0,000839409198431186i
alfa*V3^k+1= 1,24605316409672-0,0396294239640797i
V3acel= 1,03393132723175-0,0387900147656485i
|V3acel|= 1.034658714

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,000104261293770058-0,0118923034484003i
(|V3^(k+1) - V3^k| ) ( -0,0266778570931099-0,0345929687734926i

7. PROCESO ITERATIVO 3

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,00149929128372-23,4113824628275i
Y23*V3= -0,639634831262639+10,2753497891912i
Y21*V1+Y22*V2+Y23*V3= -0,123284054830409+1,7154524748786i
V2*= 1,04989068157448+0,0151511300924737i
V2*{ }= -0,155425823967505+1,79916967530586i
Qiny2= 0,155425823967505-1,79916967530586i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,194574176032495-1,79916967530586i
Sesp2*= -0,194574176032495+1,79916967530586i
V2*= 1,04989068157448+0,0151511300924737i
Sesp2*/V2*= -0,160564317898009+1,71599047195347i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,00149929128372-23,4113824628275i
Y23*V3= -0,639634831262639+10,2753497891912i
[Y21V1+Y22V2+Y23V3]= -2,12478334611413+25,1268349377061i
{-[ ]}= 2,12478334611413-25,1268349377061i
[ + ]= 1,96421902821612-23,4108444657526i
V2= 1,04970091473136-0,0168064031801335i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209978136314896+0,00303022601849474i
alfa*V2^k+1= 1,25964109767763-0,0201676838161602i
V2acel= 1,04966296136273-0,0171374577976655i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04980285
|V2esp|/|V2acel|= 1,00018779669018
V2correg= 1,04986008459268-0,017140676155518i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03393132723175+0,0387900147656485i
Sesp3*/V3*= -0,189253067336274+0,111506142301261i
Y31xV1= 0
Y32xV2= -0,869755765383663+10,4116252693884i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,869755765383663-10,4116252693884i
[ + ]= 0,680502698047389-10,3001191270871i
V3= 1,0374373000179-0,0350341882170638i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k-0,20678626544635+0,0077580029531297i
alfa*V3^k+1= 1,24492476002148-0,0420410258604766i
V3acel= 1,03813849457513-0,0342830229073469i
|V3acel|= 1.038704414

7.5 Verificar la convergencia


(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,0000305969817999685-0,0019895460630443i
(|V3^(k+1) - V3^k| ) ( 0,00420716734337989+0,0045069918583016i

7. PROCESO ITERATIVO 4

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 1,95713044921288-23,4151333252358i
Y23*V3= -0,688424025249169+10,3125425432264i
Y21*V1+Y22*V2+Y23*V3= -0,216442090887779+1,7488943665055i
V2*= 1,04986008459268+0,017140676155518i
V2*{ }= -0,257211143815341+1,8323844237768i
Qiny2= 0,257211143815341-1,8323844237768i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,092788856184659-1,8323844237768i
Sesp2*= -0,092788856184659+1,8323844237768i
V2*= 1,04986008459268+0,017140676155518i
Sesp2*/V2*= -0,0598703024062671+1,74633807699528i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 1,95713044921288-23,4151333252358i
Y23*V3= -0,688424025249169+10,3125425432264i
[Y21V1+Y22V2+Y23V3]= -2,17357254010066+25,1640276917413i
{-[ ]}= 2,17357254010066-25,1640276917413i
[ + ]= 2,11370223769439-23,417689614746i
V2= 1,05067027452481-0,0101900854873814i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209972016918536+0,0034281352311036i
alfa*V2^k+1= 1,26080432942977-0,0122281025848577i
V2acel= 1,05083231251123-0,0087999673537541i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.050869159
|V2esp|/|V2acel|= 0,999172914565602
V2correg= 1,04996318441156-0,00879268902893263i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03813849457513+0,0342830229073469i
Sesp3*/V3*= -0,189012512060244+0,110224681298601i
Y31xV1= 0
Y32xV2= -0,952511182299241+10,4043807258857i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,952511182299241-10,4043807258857i
[ + ]= 0,763498670238997-10,2941560445871i
V3= 1,03767159392011-0,0266699556910215i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,207627698915026+0,00685660458146938i
alfa*V3^k+1= 1,24520591270413-0,0320039468292258i
V3acel= 1,0375782137891-0,0251473422477564i
|V3acel|= 1.037882912

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 0,000103099818880015+0,00834798712658537i
(|V3^(k+1) - V3^k| ) ( -0,000560280786029876+0,0091356806595905i

7. PROCESO ITERATIVO 5

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,1432424834081-23,3988320198065i
Y23*V3= -0,778321575555976+10,2979499803354i
Y21*V1+Y22*V2+Y23*V3= -0,120227606999366+1,7506031090438i
V2*= 1,04996318441156+0,00879268902893263i
V2*{ }= -0,141627069850141+1,83701169105137i
Qiny2= 0,141627069850141-1,83701169105137i

7.2 Calcular la correcion al vector de estado (V.E.)


V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,208372930149859-1,83701169105137i
Sesp2*= -0,208372930149859+1,83701169105137i
V2*= 1,04996318441156+0,00879268902893263i
Sesp2*/V2*= -0,183792864166332+1,75113542250737i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,1432424834081-23,3988320198065i
Y23*V3= -0,778321575555976+10,2979499803354i
[Y21V1+Y22V2+Y23V3]= -2,26347009040747+25,1494351288503i
{-[ ]}= 2,26347009040747-25,1494351288503i
[ + ]= 2,07967722624114-23,3982997063429i
V2= 1,04965673898398-0,0116167493260894i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209992636882312+0,00175853780578653i
alfa*V2^k+1= 1,25958808678078-0,0139400991913073i
V2acel= 1,04959544989847-0,0121815613855208i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049666137
|V2esp|/|V2acel|= 1,00031806602238
V2correg= 1,04992929054833-0,0121854359262971i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,0375782137891+0,0251473422477564i
Sesp3*/V3*= -0,190123324194575+0,108646904139154i
Y31xV1= 0
Y32xV2= -0,918886070579563+10,4074042984253i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,918886070579563-10,4074042984253i
[ + ]= 0,728762746384988-10,2987573942861i
V3= 1,03778422592714-0,0301916813667382i

Acelerar el valor corregido Vi


V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,20751564275782+0,00502946844955128i
alfa*V3^k+1= 1,24534107111257-0,0362300176400858i
V3acel= 1,03782542835475-0,0312005491905345i
|V3acel|= 1.03829432

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,0000338938632300501-0,00339274689736447i
(|V3^(k+1) - V3^k| ) ( 0,0002472145656498-0,0060532069427781i

7. PROCESO ITERATIVO 6

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06762160588125-23,4056354139063i
Y23*V3= -0,718633600444955+10,306390923503i
Y21*V1+Y22*V2+Y23*V3= -0,136160509415195+1,7522406581116i
V2*= 1,04992929054833+0,0121854359262971i
V2*{ }= -0,164310723317866+1,83806961587788i
Qiny2= 0,164310723317866-1,83806961587788i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,185689276682134-1,83806961587788i
Sesp2*= -0,185689276682134+1,83806961587788i
V2*= 1,04992929054833+0,0121854359262971i
Sesp2*/V2*= -0,156519665303353+1,75247694563289i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06762160588125-23,4056354139063i
Y23*V3= -0,718633600444955+10,306390923503i
[Y21V1+Y22V2+Y23V3]= -2,20378211529644+25,1578760720179i
{-[ ]}= 2,20378211529644-25,1578760720179i
[ + ]= 2,04726244999309-23,405399126385i
V2= 1,0498282143777-0,0130896543840686i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209985858109666+0,00243708718525942i
alfa*V2^k+1= 1,25979385725324-0,0157075852608823i
V2acel= 1,04980799914357-0,0132704980756229i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049891871
|V2esp|/|V2acel|= 1,00010299042765
V2correg= 1,04991611931835-0,0132718648098948i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03782542835475+0,0312005491905345i
Sesp3*/V3*= -0,189412426425283+0,109708562774469i
Y31xV1= 0
Y32xV2= -0,908116308138022+10,4083495623697i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,908116308138022-10,4083495623697i
[ + ]= 0,718703881712739-10,2986409995952i
V3= 1,03767187007868-0,0311970003370854i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,20756508567095+0,0062401098381069i
alfa*V3^k+1= 1,24520624409442-0,0374364004045025i
V3acel= 1,03764115842347-0,0311962905663956i
|V3acel|= 1.038110005

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,0000131712299800135-0,0010864288835977i
(|V3^(k+1) - V3^k| ) ( -0,000184269931279957+4,25862413890032E-06i
7. PROCESO ITERATIVO 7

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,04340104789105-23,407762396079i
Y23*V3= -0,718493319563875+10,3045622522783i
Y21*V1+Y22*V2+Y23*V3=-0,160240786524315+1,7482850047142i
V2*= 1,04991611931835+0,0132718648098948i
V2*{ }= -0,191442386975862+1,83342591355621i
Qiny2= 0,191442386975862-1,83342591355621i

5+0,0493026157768295i 7.2 Calcular la correcion al vector de estado (V.E.)


5+0,137842993394408i
+0,107948564427617i V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,158557613024138-1,83342591355621i
Sesp2*= -0,158557613024138+1,83342591355621i

V2*= 1,04991611931835+0,0132718648098948i
Sesp2*/V2*= -0,128924456136964+1,74788913870586i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,04340104789105-23,407762396079i
Y23*V3= -0,718493319563875+10,3045622522783i
[Y21V1+Y22V2+Y23V3]= -2,20364183441536+25,1560474007932i
{-[ ]}= 2,20364183441536-25,1560474007932i
[ + ]= 2,0747173782784-23,4081582620873i
V2= 1,0500730351814-0,0118811456287278i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
0 (1-alfa)*V2^k= -0,20998322386367+0,00265437296197896i
alfa*V2^k+1= 1,26008764221768-0,0142573747544734i
99009900990099+9,9009900990099i V2acel= 1,05010441835401-0,0116030017924944i

9009900990099-9,8949900990099i 7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.05016852
|V2esp|/|V2acel|= 0,999839531106601
V2correg= 1,04993590926004-0,0116011398716367i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03764115842347+0,0311962905663956i
Sesp3*/V3*= -0,189445920926263+0,109728275039168i
Y31xV1= 0
Y32xV2= -0,924677733211557+10,4068913192792i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,924677733211557-10,4068913192792i
[ + ]= 0,735231812285294-10,29716304424i
V3= 1,03768946451991-0,0295284276383251i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207528231684694+0,00623925811327912i
alfa*V3^k+1= 1,24522735742389-0,0354341131659901i
V3acel= 1,0376991257392-0,029194855052711i
|V3acel|= 1.03810973

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 0,0000197899416900515+0,00
(|V3^(k+1) - V3^k| ) ( 0,0000579673157301563+0,00

7. PROCESO ITERATIVO 8

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,08064671169466-23,4044811413311i
Y23*V3= -0,738366906150584+10,3031545667769i
Y21*V1+Y22*V2+Y23*V3=-0,142868709307414+1,7501585739607i
V2*= 1,04993590926004+0,0116011398716367i
V2*{ }= -0,17030682262555+1,83589689382073i
Qiny2= 0,17030682262555-1,83589689382073i

7.2 Calcular la correcion al vector de estado (V.E.)


V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,17969317737445-1,83589689382073i
Sesp2*= -0,17969317737445+1,83589689382073i
V2*= 1,04993590926004+0,0116011398716367i
Sesp2*/V2*= -0,151807549133194+1,75025734259048i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,08064671169466-23,4044811413311i
Y23*V3= -0,738366906150584+10,3031545667769i
[Y21V1+Y22V2+Y23V3]= -2,22351542100207+25,1546397152918i
{-[ ]}= 2,22351542100207-25,1546397152918i
[ + ]= 2,07170787186888-23,4043823727013i
V2= 1,04989175223166-0,0119981658881299i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,209987181852008+0,00232022797432734i
alfa*V2^k+1= 1,25987010267799-0,0143977990657559i
V2acel= 1,04988292082598-0,0120775710914286i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04995239
|V2esp|/|V2acel|= 1,00004534768056
V2correg= 1,0499305305813-0,0120781187812644i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,0376991257392+0,029194855052711i
) Sesp3*/V3*= -0,189657258906037+0,10936270234422i
) Y31xV1= 0
Y32xV2= -0,919949844325402+10,4073103213805i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,919949844325402-10,4073103213805i
[ + ]= 0,730292585419365-10,2979476190363i
V3= 1,0377185170381-0,0300304990588845i

Acelerar el valor corregido Vi


V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,20753982514784+0,0058389710105422i
alfa*V3^k+1= 1,24526222044572-0,0360365988706614i
V3acel= 1,03772239529788-0,0301976278601192i
|V3acel|= 1.03816168

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -5,37867873995346E-06-0,0004
(|V3^(k+1) - V3^k| ) ( 0,0000232695586799281-0,001

7. PROCESO ITERATIVO 9

7.1 Calcular Qiny

Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,07001396997936-23,405423952535i
Y23*V3= -0,728461501679889+10,3043778028108i
Y21*V1+Y22*V2+Y23*V3=-0,143596046552019+1,7504389987907i
V2*= 1,0499305305813+0,0120781187812644i
V2*{ }= -0,17190788349249+1,83610497664374i
Qiny2= 0,17190788349249-1,83610497664374i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,17809211650751-1,83610497664374i
Sesp2*= -0,17809211650751+1,83610497664374i
V2*= 1,0499305305813+0,0120781187812644i
Sesp2*/V2*= -0,149485402607074+1,7505067483607i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,07001396997936-23,405423952535i
Y23*V3= -0,728461501679889+10,3043778028108i
[Y21V1+Y22V2+Y23V3]= -2,21361001653138+25,1558629513257i
{-[ ]}= 2,21361001653138-25,1558629513257i
[ + ]= 2,06412461392431-23,405356202965i
V2= 1,04990131873476-0,0123396875322646i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,20998610611626+0,00241562375625288i
alfa*V2^k+1= 1,25988158248171-0,0148076250387175i
V2acel= 1,04989547636545-0,0123920012824646i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04996861
|V2esp|/|V2acel|= 1,00002990021246
V2correg= 1,04992686846325-0,0123923718059357i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,03772239529788+0,0301976278601192i
) Sesp3*/V3*= -0,189542162322847+0,109540170496768i
) Y31xV1= 0
Y32xV2= -0,916834802380092+10,4075852043945i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,916834802380092-10,4075852043945i
[ + ]= 0,727292640057245-10,2980450338977i
V3= 1,03769822882881-0,0303316472153118i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207544479059576+0,00603952557202384i
alfa*V3^k+1= 1,24523787459457-0,0363979766583742i
V3acel= 1,03769339553499-0,0303584510863504i
|V3acel|= 1.03813738

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -3,66211805014061E-06-0,0003
(|V3^(k+1) - V3^k| ) ( -0,0000289997628899474-0,00

7. PROCESO ITERATIVO 10

7.1 Calcular Qiny


Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06300842525917-23,4060424777688i
Y23*V3= -0,726840479872758+10,3042499073626i
Y21*V1+Y22*V2+Y23*V3=-0,148980569465078+1,7496925781087i
V2*= 1,04992686846325+0,0123923718059357i
V2*{ }= -0,17810154373435+1,83520302669839i
Qiny2= 0,17810154373435-1,83520302669839i

7.2 Calcular la correcion al vector de estado (V.E.)

V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,17189845626565-1,83520302669839i
Sesp2*= -0,17189845626565+1,83520302669839i
V2*= 1,04992686846325+0,0123923718059357i
Sesp2*/V2*= -0,143073278579984+1,74962285387627i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06300842525917-23,4060424777688i
Y23*V3= -0,726840479872758+10,3042499073626i
[Y21V1+Y22V2+Y23V3]= -2,21198899472425+25,1557350558775i
{-[ ]}= 2,21198899472425-25,1557350558775i
[ + ]= 2,06891571614427-23,4061122020012i
V2= 1,04995624789793-0,0121300143674224i

7.3 Acelerar el valor corregido Vi

V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,20998537369265+0,00247847436118714i
alfa*V2^k+1= 1,25994749747752-0,0145560172409069i
V2acel= 1,04996212378487-0,0120775428797198i

7.4 Calcular el valor de la tension corregida en nodos PV

V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.05003158
|V2esp|/|V2acel|= 0,999969920438937
V2correg= 1,04993054138505-0,0120771795925313i
|V2correg=| 1.05000000

Estimar el valor de V3

V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,03769339553499+0,0303584510863504i
) Sesp3*/V3*= -0,189529544482794+0,109572238122028i
) Y31xV1= 0
Y32xV2= -0,919959153920532+10,4073094984584i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,919959153920532-10,4073094984584i
[ + ]= 0,730429609437738-10,2977372603364i
V3= 1,03769884055141-0,0300146823994633i

Acelerar el valor corregido Vi

V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)

(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207538679106998+0,00607169021727008i
alfa*V3^k+1= 1,24523860866169-0,036017618879356i
V3acel= 1,03769992955469-0,0299459286620859i
|V3acel|= 1.03813193

7.5 Verificar la convergencia

(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 3,67292180003531E-06+0,0003
(|V3^(k+1) - V3^k| ) ( 6,53401969996636E-06+0,0004
)
)
)
)
)
)
)
)
)
. PROCESO ITERATIVO 7

2850047142i

0474007932i
8913192792i

)
,0000197899416900515+0,0016707249382581i )
,0000579673157301563+0,0020014355136846i )

. PROCESO ITERATIVO 8

1585739607i
6397152918i

3103213805i
)
5,37867873995346E-06-0,000476978909627699i )
,0000232695586799281-0,0010027728074082i )

. PROCESO ITERATIVO 9

4389987907i

8629513257i
5852043945i

)
3,66211805014061E-06-0,0003142530246713i )
0,0000289997628899474-0,000160823226231201i )

. PROCESO ITERATIVO 10
6925781087i

7350558775i
3094984584i

)
,67292180003531E-06+0,000315192213404399i ) < 0,00001
,53401969996636E-06+0,000412522424264501i )

También podría gustarte