Facultad de Ingeniería de sistemas y eléctrica

Carrera de Ingeniería Electrica y de Potencia

Proyecto del Curso Turbomáquinas

Diseño de una microcentral eléctrica para la electrificación rural de 40 casas

en el distrito de Carhuamayo, Junín.

Autores:

Agelvis Vargas, Jose Alejandro (U20230741)

Sante Mancilla, Carlos Andre (U20203213)

Docente:

Villarreal, Luis

Lima, Perú

Noviembre; 2024
 ii

INDICE DE CONTENIDO

INDICE DE TABLAS ............................................................................................................. iv

INTRODUCCION ...................................................................... ¡Error! Marcador no definido.

CAPITULO 1: INTRODUCCION ............................................ ¡Error! Marcador no definido.

ASPECTOS GENERALES...................................................... ¡Error! Marcador no definido.

1.1. Definición del Problema .......................................¡Error! Marcador no definido.

1.1.1. Descripción del Problema ............................¡Error! Marcador no definido.

1.1.2. Planteamiento del problema ........................¡Error! Marcador no definido.

1.2. Definición de objetivos..........................................¡Error! Marcador no definido.

1.2.1. Objetivo general ............................................¡Error! Marcador no definido.

1.2.2. Objetivos específicos ....................................¡Error! Marcador no definido.

1.3. Alcances y limitaciones ........................................¡Error! Marcador no definido.

1.3.1. Alcances .........................................................¡Error! Marcador no definido.

1.3.2. Limitaciones ...................................................¡Error! Marcador no definido.

1.4. Justificación ............................................................¡Error! Marcador no definido.

1.5. Estudios de viabilidad ...........................................¡Error! Marcador no definido.

BIBLIOGRAFÍAS .................................................................................................................. 36

ANEXOS ................................................................................................................................. 38

ANEXO 1 Lista de dispositives........................................................................................ 38

 iii

INDICE DE FIGURAS

Figura 1. Producción acumulada durante el 2023 ...... ¡Error! Marcador no definido.

Figura 2. Grafica de sistema con seguimiento solar vs sistema fijo¡Error! Marcador

no definido.

Figura 3. Imagen termográfica de un punto caliente ocasionado por suciedad ¡Error!

Marcador no definido.
 iv

INDICE DE TABLAS

Tabla 2. Funciones del controlador ............................ ¡Error! Marcador no definido.

CAPITULO 1: INTRODUCCION

En la región de Junín se encuentra la localidad de Carhuamayo,

comprendida por 40 viviendas, entre ellas algunos talleres industriales como

metalurgia y panadería. Se busca establecer una micro central hidráulica para

satisfacer la demanda electrica. En la figura 1 se muestra la localidad en el

noroeste de Junin,

Figura 1. Localidad de Carhuamayo

Fuente: Google Earth

1.1. Definición del Problema

Ante la situación mencionada se presenta la siguiente problemática ¿Cómo

establecer una micro central hidráulica con los recursos disponibles para

satisfacer la demanda electrica de la población?

1.2. Definición de objetivos

1.2.1. Objetivo general

Establecer un suministro eléctrico mediante una micro central hidráulica

1.2.2. Objetivos específicos

Determinar la máxima demanda requerida por la población

Establecer los parámetros técnicos de la turbina

 2

Determinar el sistema de transformación y distribución electrica

Analizar el comportamiento del flujo de cargas

1.3. Alcances

El presente proyecto tiene como alcance determinar el sistema de generación

electrica mediante una micro central hidráulica para abastecer de energía

electrica a una población ubicada en Carhuamayo, Junín

1.4. Limitaciones

El presente proyecto tiene como limitaciones los cálculos eléctricos y

mecánicos para el desarrollo de la micro central hidráulica

 3

CAPITULO 2: MARCO TEORICO

2.1. Determinación de la máxima demanda

La máxima demanda para zonas rurales está definida por el Ministerio de

Energía y Minas a través de la norma “Calificación de los Sistemas de

Distribución Eléctrica en Sectores de Distribución Típicos” donde establecen

la potencia recomendada para el diseño de proyectos de electrificación rural,

teniendo en cuenta la demanda de la localidad, su ámbito geográfico y la

accesibilidad técnica – comercial.

2.2. Cálculo de la potencia de la turbina

Mediante la obtención de la electricidad se producen intercambios de energía

entre la turbina y el medio, ocasionando una reducción en la eficiencia de la

máquina. Por ello, al momento de establecer la potencia mecánica que debe

suministrar la turbina, será necesario considerar algunos factores de

eficiencia para obtener una potencia efectiva de generación como se muestra

en la siguiente ecuación (1)

𝑷 = 𝝆. 𝑯𝒏 . 𝑸. ɳ Ecu. (1)

Donde:

P Potencia mecánica producida en el eje de la turbina [W]

𝝆 Densidad de agua [1.000 kg/m^3]

𝑸 Caudal [m^3/s]

H Altura teórica [m]

ɳ Eficiencia manométrica
 4

De acuerdo con la altura teórica y caudal disponible es posible determinar la

turbina que proporcionará el mayor rendimiento como se muestra en la

siguiente figura 2

Figura 2. Seleccion de turbina en función de altura y caudal

Fuente: Journal of the Institute of Engineering

2.3. Cálculo de la velocidad especifica

La velocidad especifica es un parámetro adimensional que facilita la elección

de la maquina mas apropiada, se determina a partir de variables medibles

como la velocidad de rotación del eje, el caudal aplicado a la turbina y la altura

manométrica. Se representa con la siguiente ecuación 2

𝒏√𝑸 Ecu. (2)

𝒏𝒔 = 𝟑
(𝒈𝑯𝒎 )𝟒

Donde:

𝒈 Gravedad [9.81 m/s^2]

𝑸 Caudal [m^3/s]
 5

𝑯𝒎 Altura manometrica [m]

𝒏 Velocidad de rotacion [rpm]

2.4. Tipo de rodete por velocidad especifica

CAPITULO 3: DESARROLLO

3.1. Cálculo de cargas

De acuerdo con los datos proporcionados, se conoce que; existen 40

casas en el centro poblado ubicado en Carhuamayo, Junín de las cuales 2

funcionan como talleres (01 metalúrgico y 1 panadería). Se conoce que se

encuentran en cercanía del Rio Carhuamayo, con un caudal de 17.14 m3/s.

Con estos datos es posible determinar el sector tipico vigente para realizar el

cálculo de cargas como se muestra en la figura 1 de la Norma DGE

 6

“Calificacion electrica para la elaboracion de proyectos de subsistemas de

distribucion secundaria”

Figura 3. Sector típico vigente Carhuamayo, Junin

Fuente: Ministerio De Energía Y Minas

En este sentido, es posible determinar la demanda maxima de la

instalacion con el sector tipico y el tipo de habilitacion, en la figura 4 se

muestra la calificacion electrica para la elaboracion de proyectos de

subsistemas de distribucion

Figura 4. Calificación electrica por sector típico

Fuente: Ministerio De Energía Y Minas
En el caso de estudio, se evaluará el sector de distribucion tipico 2

correspondiente a Santa. Para las 38 casas se considera habilitaciones para

viviendas en vias de regularizacion (parcial o totalmente edificadas)

calificados como asentamientos humanos marginales o pueblos jovenes, con

suministro monofasico de 250 W cada vivienda. Mientras que, para las 2

casas que funcionan como taller se considerará 1000 W por predio,

resultando una maxima demanda mostrada en la siguiente tabla 1

 7

Tabla 1. Cálculo de cargas

Tipo de habilitación Potencia [W] Cantidad Total [W]

Vivienda 250 38 9,500

Vivienda taller 1,000 2 2,000

Alumbrado publico 10 50 500

Total 12,000

A la potencia calculada de 12,000 W se debe adicionar una 20% por factor

de seguridad, resultando una potencia necesaria de 14,400 W

3.2. Cálculo de la potencia de la turbina

Se debe verificar que la potencia disponible pueda satisfacer la demanda

electrica de la población, para ello se emplea la siguiente ecuación:

𝑷 = 𝒑𝒈𝑸𝑯

De acuerdo con un estudio de la municipalidad de Satipo (2008), el rio

Carhuamayo cuenta con un caudal de 17.14 m3/s. Una velocidad media de

1.04m/s. una profundidad de 0.25m y se encuentra ubicado a una altitud de

4.126msnm. Se tiene conocimiento que la altura neta disponible para la micro

central es de 5 metro.

𝒌𝒈 𝒎 𝒎𝟑
𝑷 = 𝟏𝟎𝟎𝟎 𝟑 × 𝟗. 𝟖𝟏 𝟐 × 𝟏𝟕. 𝟏𝟒 ×𝟏𝒎
𝒎 𝒔 𝒔

𝑷 = 𝟏𝟔𝟖. 𝟏𝟒𝐤𝐖

De acuerdo con el resultado obtenido podemos verificar que el rio puede

ofrecer una mayor potencia de la requerida de 14.4 kW. Teniendo en cuenta

que la región de Junin es afectada durante la época de estiaje, se puede

alcanzar un nivel mínimo de caudal de 1.47 m3/s para satisfacer la demanda

 8

prevista. Tomando como referencia la grafica de tablas de acuerdo con la

altura y el caudal disponible, se tomará la turbina kaplan debido a su

capacidad de operar en bajos saltos de altura y poder regularse ante

variaciones del caudal. En la siguiente figura se puede observar el diseño del

canal recomendado para esta aplicacion

Figura 5. Turbina kaplan

Fuente: Ciencia Latina

3.3. Cálculo de la velocidad especifica

Teniendo en cuenta que la turbina Kaplan tendrá una velocidad de rotacion de

360rpm considerando que las centrales hidráulicas operan a bajas revoluciones.

En este sentido, la velocidad especifica resulta:

𝒏√𝑸
𝒏𝒔 = 𝟑
(𝒈𝑯𝒎 )𝟒

𝟑𝟔𝟎𝒓𝒑𝒎√𝟏. 𝟒𝟕 m3 /𝑠
𝒏𝒔 = 𝟑
(𝟗. 𝟖𝟏m/𝑠 2 𝟏𝒎)𝟒

𝒏𝒔 = 𝟕𝟖. 𝟕𝟒𝒓𝒑𝒎
 9

3.4. Cálculo de la toma de agua

3.5. Determinación los materiales a emplear

De acuerdo con el ultimo estudio realizado por Source International en los

años 2009 y 2016 en la laguna Yanacocha (del cual nace el rio Carhuamayo),

se detectaron trazas de zinc, boro, bario y otros cuatro elementos en

concentraciones por debajo de los estandares nacionales e internacionales

de calidad. Sin embargo, la laguna de Quiulacocha representa una de las

principales fuentes de contaminacion por la presencia de agua acidas (pH de

2.4) con un elevado contenido de minerales disueltos, producto de

meteorizacion de los relaves. Son muy ricas en manganeso, zinc, hierro,

cadmio, arsenico y mercurio. En la siguiente figura 5 se aprecia la toma de

muestras por parte de los investigadores

Figura 6. Toma de muestras para análisis fisicoquimicos

Fuente: Asociación Civil Centro De Cultura Popular Labor
 10

Por este motivo se considera una tuberia y componentes de acero inoxidable

304 por su resistencia a agente nocivos (ya que el mercurio disuelve el

aluminio). Ademas, se conoce que el rio Carhuamayo presenta un fondo

pedregoso y color turbio. En la figura 6 se puede observar la reaccion quimica

que ocasiona el mercurio sobre el aluminio.

Figura 7. Efecto del aluminio sobre el aluminio

Fuente: Ciencia de sofa
 11

3.6. Diseño de la bocatoma Convencional

Para el diseño de la bocatoma debemos determinar lo siguiente:

1. Ancho de la captación (B1)

2. Número de barrotes (N) y espaciamiento (b´)
3. Ancho total de la captación (B)
4. Ancho del vertedero de desripiador (B2)
5. Longitud del desripiador (Ld)
6. Longitud de transición (Lt)
7. Ancho de compuerta de lavado (Bc)

Hallando el ancho de la captación (B1)

𝑄 = (𝐵1 )(𝐾 )(𝑆)(𝑀)(𝐻 3)

Donde:
Q = Caudal
B1 = Ancho de la captación
K = Coeficiente de pérdida del vertedero
S = Coeficiente de sumersión
 12

M = Coeficiente del vertedero

H = Carga sobre la cresta de captación
3⁄ 0.385
ℎ𝑛 2
𝑆 = [1 + ( ) ]
𝑌2

Donde:
S = Coeficiente de sumersión
hn = Altura de flujo que pasa sobre el muro
Y2 = Altura de agua acumulada
2
0.0133 𝐻
𝑀 = [1.794 + ] [1 + 0.55 ( ) ] [√2𝑔]
𝐻 + 𝑌1 𝐻 + 𝑌1

Donde:
M = Coeficiente de vertedero
H = Carga sobre la cresta
Y1 = Altura de la ventana
g = Gravedad

3⁄ 0.385
0.5 2
𝑆 = [1 + ( ) ]
1.5

𝑆 = 1.07

2
0.0133 1
𝑀 = [1.794 + ] [1 + 0.55 ( ) ] [√(2)(9.81)]
1 + 1.2 1 + 1.2
 13

𝑀 = 8.88

𝑄
𝐵1 =
(𝐾 )(𝑆)(𝑀)(𝐻 3)

𝐵1 = 2. 𝑚

Hallando el número de barrotes (N) y espaciamiento entre ellos (b´)

𝐵1
𝑁= −1=𝑛−1
𝑏´
Donde:

N = Número de barrotes
b´ = Separación entre los barrotes
B1 = Ancho de la captación
2
𝑁= −1= 𝑛−1
0.2

𝑁=9

Hallando el ancho total de la captación (B)

𝐵 = 𝐵1 + (𝑁)(𝑎´)

Donde:
B = Ancho total de la captación
B1 = Ancho de la captación
N = Número de barrotes
a´ = Ancho de los barrotes

𝐵 = 2 + (9)(0.05)

𝐵 = 2.45 𝑚

Hallando el ancho del vertedero de disripiador (B2)

 14

𝑄
𝐵2 = 3⁄
(𝑆)(𝑀) (𝐻 2)

Donde:
Q = Caudal
S = Coeficiente de sumersión
M = Coeficiente del vertedero
H = Carga sobre la cresta
B2 = Ancho de vertedero de disripiador
17.14
𝐵2 = 3⁄
(1.07)(8.88) (1 2)

𝐵2 = 1.8 𝑚

Hallando la longitud de desripiador (Ld)

𝐵 − 𝐵2
𝐿𝑑 =
(2)(𝐭𝐚𝐧(12.5°))

Donde:

B = Ancho total de la captación

B2 = Ancho de vertedero de disripiador

 15

Ld = Longitud de disripiador

2.45 − 1.8
𝐿𝑑 =
(2)(𝐭𝐚𝐧(12.5°))

𝐿𝑑 = 1.46 𝑚

Hallando la longitud de transición (Lt)

𝐵2 − 𝐵3
𝐿𝑡 =
(2)(tan(12.5°))

Donde:

B2 = Ancho de vertedero de disripiador

B3 = Ancho del canal de conducción

Lt = Longitud de transición

1.8 − 1
𝐿𝑡 =
(2)(𝐭𝐚𝐧(12.5°))

𝐿𝑡 = 1.8 𝑚
 16

Hallando el ancho de la compuerta de lavado (Bc)

𝑣2
𝑄 = (𝐾𝑐 )(𝑒)(𝑎)(𝐵𝑐 ) (√(2𝑔) (𝐻𝑑 + − (𝑒)(𝑎)))
2𝑔

Donde:

Q = Caudal

Kc = Coeficiente entre 0.95 y 0.97

e = Coeficiente de relación a/Hd

a = Apertura de la compuerta

Bc = Ancho de la compuerta

Hd = Altura del nivel de agua en el desripiador

v = Velocidad del caudal

g = gravedad

Utilizamos el siguiente gráfico para hallar el valor de “e”

 17

𝑎 𝑎 1.7
= =
𝐻𝑑 𝑌2 + ℎ𝑛 2

17.14
𝐵𝑐 =
22
(0.96)(0.745)(1.7) (√((2)(9.81)) (2 + ( )( )
(2)(9.81) − 0.745 1.7 ))

𝐵𝑐 = 3.28 𝑚
 18

3.7. Diseño del canal de distribución

Se diseñarán 2 canales, uno desde la bocatoma hasta el desarenador y otro

desde el desarenador hasta la cámara de carga, asumiendo una longitud de

100m entre cada uno. A su vez, la variación de las alturas entre cada uno

será de 1m.

Para el canal de distribución se utilizará hormigón afinado, por lo que

tendremos los siguientes datos:

Coeficiente de Manning 0.012

Pendiente del canal (i) 0.0004 < i < 0.0005

Velocidad del agua (v´) 0.65 < v´ < 2.5

Para el diseño del canal 1 (Bocatoma - Desarenador) tomaremos el valor de

la pendiente de i1 = 0.0005.

Mientras que para el diseño del canal 2 (Desarenador – Cámara de carga)

usaremos la siguiente ecuación

1
𝑖2 = = 0.01
100

Diseñando el canal 1:
 19

Longitud del canal 100

Caudal 17.14

Ancho 10

Coeficiente de Manning 0.012

Pendiente 0.0005

Profundidad de control de aguas abajo 2.2

Talud 0

Para hallar la altura normal y crítica usamos el programa Hcanales V3

Altura normal: 1.0247 m

 20

Altura crítica: 0.6690 m

Resultados del canal 1

Pendiente 0.0005

Material Hormigón

Tipo de sección Rectangular

Ancho 10

Área 10.25

Altura normal 1.0247

Altura crítica 0.6690

Profundidad 2.5

Longitud del canal 100

 21

3.8. Diseño de la cámara de carga

Hallando el caudal del estanque:

(0.693)(𝑄 2 )
𝑉𝑡 =
(𝐴)(𝑖 )(𝑔)

Donde:
Vt = Volumen del estanque
Q = Caudal
A = Área del canal de conducción
i = Pendiente del canal de llegada
g = Gravedad
(0.693)(17.142 )
𝑉𝑡 =
(10.25)(0.01)(9.81)

𝑉𝑡 = 202.47 𝑚3

Hallando las dimensiones del estanque:

𝑉𝑡
𝐻= +𝑘
3

𝑉𝑡
𝐵=√
𝐻

Donde:
H = Profundidad del estanque
B = Ancho del estanque
Vt = Volumen del estanque
k = Constante de capacidad del estanque (2)
 22

202.47
𝐻= +2
3

𝐻 = 69.49 𝑚

202.47
𝐵=√
69.49

𝐵 = 1.71 𝑚

Hallando la transición de entrada:

𝑇1 − 𝑇2
𝐿𝑇 =
(2)(𝐭𝐚𝐧(22.5°))

Donde:

T1 = Ancho del canal

T2 = Ancho de la cámara

LT = Largo de transición

10 − 1.71
𝐿𝑇 =
(2)(𝐭𝐚𝐧(22.5°))

𝐿 𝑇 = 10 𝑚

Hallando la altura mínima

𝐻𝑚𝑖𝑛 = (𝐶 )(𝑉)(√𝐷)

Donde:
V = Velocidad del flujo
C = Constante (0.3 – 0.4)
Hmin = Altura mínima
D = Diámetro de la tubería forzada
 23

𝐻𝑚𝑖𝑛 = (0.35)(1.67)(√1)

𝐻𝑚𝑖𝑛 = 0.5845 𝑚

3.9. Diseño de la tubería de presión

𝟑
𝑸 = 𝟏𝟕. 𝟏𝟒 𝒎 ⁄𝒔

𝑳 = 𝟏𝟎𝟎 𝒎

𝑯𝒏 = 𝟏𝟎 𝒎

𝑯𝒇 = 𝟑. 𝟕𝟔 (𝒑é𝒓𝒅𝒊𝒅𝒂 𝒅𝒆 𝒄𝒂𝒓𝒈𝒂 𝒅𝒆𝒍 𝟓%)

𝒏𝒎 = 𝟎. 𝟎𝟏𝟐

Diámetro de la tubería

𝟕 (𝟓. 𝟐)(𝑸𝟑 )
𝑫= √
𝑯𝒏

𝟕 (𝟓. 𝟐)(𝟏𝟕. 𝟏𝟒𝟑 )

𝑫= √
𝟏𝟎

𝑫 = 𝟑. 𝟎𝟕𝟖 𝒎

Cálculo de la velocidad de la tubería

 24

(𝟒)(𝑸)
𝑽=
(𝝅)(𝑫𝟐 )

(𝟒)(𝟏𝟕. 𝟏𝟒)
𝑽=
(𝝅)(𝟑. 𝟎𝟕𝟖𝟐 )

𝑽 = 𝟐. 𝟑𝟎𝟑 𝒎/𝒔

Cálculo de la altura bruta

𝑯𝒃 = 𝟏𝟎 𝒎
 25

𝟕 (𝟓. 𝟐)(𝑸𝟑 )
𝑫=√
𝑯𝒏

𝟕 (𝟓. 𝟐)(𝟏𝟕. 𝟏𝟒𝟑 )

𝑫=√
𝟏𝟎

𝑫 = 𝟑. 𝟎𝟕𝟖 𝒎

3.3 Determinar el tipo de turbina a emplear

Se debe verificar que la potencia disponible pueda satisfacer la demanda

electrica de la población, para ello se emplea la siguiente ecuación.

𝑷 = 𝒑𝒈𝑸𝑯

Donde:

P Potencia mecánica producida en el eje de la turbina [W]

p Densidad de agua [1.000 kg/m^3]

g Gravedad (9.81 m/s^2)

Q Caudal

H Altura teórica [m]

3.3.1 Condiciones de la fuente hídrica

 26

De acuerdo con un estudio de la municipalidad de Satipo (2008), el rio

Carhuamayo cuenta con un caudal de 17.14 m3/s. Una velocidad media

de 1.04m/s. una profundidad de 0.25m y se encuentra ubicado a una

altitud de 4.126msnm. Se tiene conocimiento que la altura neta es de 10

metros. Por lo que, el tipo de turbina a utilizar será de tipo Kaplan, ya.

que la Altura Neta es de 10 metros y el Caudal será de 17.14 m3/s

Elección del tipo de Turbina

La elección de la turbina a utilizar, se determinará con el caudal y la altura

que se tiene de dato, con ello se podrá determinar el tipo de turbina según

el siguiente gráfico.
 27

Después de analizar el gráfico tenemos que la turbina a utilizar será una

Kaplan.

Cálculo de la Potencia Hidráulica

Para hallar la potencia hidráulica utilizaremos la siguiente fórmula:

𝑷𝑯𝑻 = 𝜼𝒈 . 𝜼𝒕 . 𝝆. 𝒈. 𝑸. 𝑯𝒏

𝑷𝑯 = (𝟎. 𝟗)(𝟎. 𝟗𝟑)(𝟏𝟎𝟎𝟎)(𝟗. 𝟖𝟏)(𝟏𝟕. 𝟏𝟒)(𝟏𝟎)

𝑷𝑯 = 𝟏. 𝟓𝟏𝟑 𝑴𝑾

La eficiencia de las turbinas Kaplan son un 93% aproximadamente, por lo

que:

Escriba aquí la ecuación.

 28

Cálculo de la velocidad específica

Para hallar la velocidad específica utilizaremos el siguiente gráfico:

Después de ver el gráfico se tiene una velocidad específica (ns) de 700

rpm

Cálculo de la velocidad síncrona

Para hallar la velocidad síncrona utilizamos la siguiente fórmula

𝒏 = 𝒏𝒔 . 𝑷−𝟎.𝟓 . 𝑯𝒏

𝒏 = (𝟕𝟎𝟎)(𝟏𝟓𝟏𝟑)−𝟎.𝟓(𝟏𝟎)𝟏.𝟐𝟓

𝒏 = 𝟑𝟐𝟎. 𝟎𝟐 𝒓𝒑𝒎

Cálculo del coeficiente periférico

Para hallar el coeficiente periférico utilizamos la siguiente fórmula

 29

𝝓 = 𝟎. 𝟕𝟗 + 𝟏. 𝟔𝟏(𝟏𝟎−𝟑 )𝒏𝒔

𝝓 = 𝟎. 𝟕𝟗 + 𝟏. 𝟔𝟏(𝟏𝟎−𝟑 )(𝟕𝟎𝟎)

𝝓 = 𝟏. 𝟗𝟏𝟕

Cálculo del diámetro máximo de la turbina

Para hallar el diámetro máximo de la turbina utilizamos la siguiente

fórmula

(𝟖𝟒. 𝟓)(𝝓)(√𝑯𝒏 )
𝑫𝑴𝒂𝒙 =
𝒏

(𝟖𝟒. 𝟓)(𝟏. 𝟗𝟏𝟕)(√𝟏𝟎)

𝑫𝑴𝒂𝒙 =
𝟑𝟐𝟎. 𝟎𝟐

𝑫𝑴𝒂𝒙 = 𝟏. 𝟔𝟎 𝒎

Dimensionamiento general de la turbina

 30

Para dimensionar la turbina, utilizaremos las siguientes ecuaciones

𝑫𝒎 𝟎. 𝟐𝟓 + 𝟗𝟒. 𝟔𝟒
=
𝑫𝑴𝒂𝒙 𝒏𝒔

𝑫𝒎 = 𝟎. 𝟐𝟐 𝒎

𝑯𝒎
= 𝟔. 𝟗𝟒(𝒏𝒔 )−𝟎.𝟒𝟎𝟑
𝑫𝑴𝒂𝒙

𝑯𝒎 = 𝟎. 𝟕𝟗 𝒎

𝑯𝟏
= 𝟎. 𝟑𝟖 + 𝟓. 𝟏𝟕(𝟏𝟎−𝟓 )(𝒏𝒔)
𝑫𝑴𝒂𝒙

𝑯𝟏 = 𝟎. 𝟔𝟒 𝒎

𝑯𝒕
= 𝟎. 𝟐𝟒 + 𝟕. 𝟖𝟐(𝟏𝟎−𝟓 )(𝒏𝒔)
𝑫𝑴𝒂𝒙

𝑯𝒕 = 𝟎. 𝟒𝟒 𝒎

Dimensionamiento de la carcasa de acero (justificar pq no usamos de

cemento)

Para dimensionar la carcasa de acero usamos las siguientes formulas

 31

𝐴1
= 0.40(𝑛𝑠 )0.20
𝐷𝑀𝑎𝑥

𝐴1 = 2.37 𝑚

𝐵1
= 1.26 + 3.79(10)−4 (𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐵1 = 2.44 𝑚

𝐶1
= 1.46 + 3.24(10)−4 (𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐶1 = 2.69 𝑚
 32

𝐷1
= 1.59 + 5.74(10)−4 (𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐷1 = 3.18 𝑚

𝐸1
= 1.21 + 2.71(10)−4 (𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐸1 = 2.24 𝑚

𝐹1 1.45 + 72.17
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝐹1 = 0.16 𝑚

𝐺1 1.29 + 41.63
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝐺1 = 0.09 𝑚

𝐻1 1.13 + 31.86
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝐻1 = 0.07 𝑚

𝐼1 0.45 − 31.80
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝐼1 = 0.07 𝑚

𝐿1
= 0.74 + 8.7(10)−4 (𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐿1 = 2.15 𝑚

𝑀1 1
=
𝐷𝑀𝑎𝑥 2.06 − 1.20(10)−3 (𝑛𝑠 )

𝑀1 = 1.31 𝑚
 33

Dimensionamiento del tubo de aspiración

Para dimensionar el tubo de aspiración usaremos las siguientes formulas

𝐻𝑡
= 0.24 + 7.82(10−5 )(𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝐻𝑡 = 0.47 𝑚

𝑁
= 2.00 − 2.14(10−6 )(𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝑁 = 3.19 𝑚

𝑂
= 1.40 − 1.67(10−5 )(𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝑂 = 2.22 𝑚

𝑃 1.26 − 16.35
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝑃 = 0.03 𝑚

𝑄 0.66 − 18.40
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝑄 = 0.04 𝑚

𝑅
= 1.25 − 7.98(10−5 )(𝑛𝑠 )
𝐷𝑀𝑎𝑥

𝑅 = 1.91 𝑚

𝑆 4.26 + 201.51
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝑆 = 0.47 𝑚

𝑇
= 1.20 + 5.12(10−4 )(𝑛𝑠 )
𝐷𝑀𝑎𝑥
 34

𝑇 = 2.49 𝑚

𝑍 2.58 + 102.66
=
𝐷𝑀𝑎𝑥 𝑛𝑠

𝑍 = 0.24 𝑚

𝑈
= 0.15
𝐷𝑀𝑎𝑥

𝑈 = 0.24 𝑚

𝑉
= 2.91
𝐷𝑀𝑎𝑥

𝑉 = 4.65 𝑚
35
 36

BIBLIOGRAFÍAS

[1].Bajracharya, Tri & Shrestha, Rajendra & Timilsina, Ashesh. (2019). A

Methodology for Modelling of Steady State Flow in Pelton Turbine

Injectors. Journal of the Institute of Engineering. 15. 246-255.

10.3126/jie.v15i2.27674. Comité de Operación Económica del Sistema.

(2023) Producción acumulada 2023.

https://www.coes.org.pe/Portal/home/

[2].Camasca, R. (2022). Mejoramiento de la eficiencia energética con sistema

On-Grid de paneles solares para una empresa textil en Villa El Salvador

– Lima. [Trabajo de suficiencia profesional de pregrado, Universidad

Nacional Mayor de San Marcos, Facultad de Ingeniería Electrónica y

Eléctrica, Escuela Profesional de Ingeniería Eléctrica]. Repositorio

institucional Cybertesis UNMSM.

[3].Tobajas, C (2020) Energía Solar Fotovoltaica 2da ed. Ecoe Ediciones.

Colombia

[4].Chero, J y Tesen, A (2021) Implementación de un servomecanismo de

posición a un panel solar utilizando lógica difusa para aumentar la

eficiencia en la captación de energía solar. [Tesis de grado, Universidad

Nacional Pedro Ruiz Gallo]. Repositorio Institucional de la Universidad

Nacional Pedro Ruiz Gallo https://hdl.handle.net/20.500.12893/9383

[5].Lozano Mestra, M. (2022). Estimación de la radiación solar incidente

sobre paneles fotovoltaicos para sistemas estáticos y sistemas dotados

con seguimiento solar

 37

[6].Cobian Lufin, E. (2023). Optimización de limpieza de paneles solares en

plantas fotovoltaicas. Disponible en

https://repositorio.uchile.cl/handle/2250/193421

[7].Sánchez Plúas, J. R., & Figueroa Sánchez, M. Á. (2022). Tesis.

Recuperado a partir de http://repositorio.ug.edu.ec/handle/redug/59832

[8].Pascual, D. L. (2023). Sistema de mejora de la eficiencia de paneles

solares fotovoltaicos mediante refrigeración geotérmica. Universidad de

Alcalá.

[9].Milosavljević, D., Kevkić, T. & Jovanović, S. (2022). Review and validation

of photovoltaic solar simulation tools/software based on case study. Open

Physics, 20(1), 431-451. https://doi.org/10.1515/phys-2022-0042

[10]. Goudelis, G., Lazaridis, P. I., & Dhimish, M. (2022). A Review of

Models for Photovoltaic Crack and Hotspot Prediction. Energies, 15(12),

4303. MDPI AG. Retrieved from http://dx.doi.org/10.3390/en15124303

[11]. Ahmed, M., Harbi, I., Kennel, R., Rodríguez, J., & Abdelrahem, M.

(2022). Maximum Power Point Tracking-Based Model Predictive Control

for Photovoltaic Systems: Investigation and New Perspective. Sensors,

22(8), 3069. MDPI AG. Retrieved from

http://dx.doi.org/10.3390/s22083069

[12]. Mateo, V. (2022) Integración de sensores y actuadores inteligentes

en ... - mincotur.gob.es, Generalitat Valenciana. Disponible en:

https://www.mincotur.gob.es/Publicaciones/Publicacionesperiodicas/Eco

nomiaIndustrial/RevistaEconomiaIndustrial/427/MASCAROS%20MATEO

.pdf
 38

ANEXOS

ANEXO 1 Lista de dispositives

Fuentes de alimentación:

Transformadores:

Designaci Tipo Sn ukr Uprim Pk P0 Grup Δu_t Aire Referencia

ón [kVA] [%] [kV]/ [kW] [kW] o de ransf forzad

Usec vect orm o

[V] ores er

[%]

Transform 250 4 20/ 3.2 1.1 Dyn 0.07 No 4GB54643C

ador 2.1 380 5 5

La Δu_del transformador [%] es independiente del punto de referencia que se

tome para el cálculo de la caída de tensión de la línea.

Renovables:

Designación Sn [kVA] Un [V] S_LF [kVA] Ikmax [kA] cos φ LF Tipo de

acometida

Renovables 2.10 60 380 48 0.109 1 Inductiva

Interruptores / Fusibles:

Mediante interruptores automáticos:

Sitio Designación Tipo Referencia de In [A] Cantidad

Catálogo Reles

MT 2.1 MT-IA 2.1 Interruptor CB-f NAR 7SJ8011 630 1

 39

Interruptores Automáticos / Termomagnéticos (PIA):

Sitio Designación Referencia In [A] Icu/Icn Icu/Icn Tipo de Cant

[kA] [kA] disparad idad

necesari or /

a caracter

ísticas

BTDP 2.1 IA 2.1B 3VA23405HN32 400 55 9.675 ETU350 1

0AA0

BTDP 2.10 IA 2.10B 3VA20105HN36 100 55 9.675 ETU350 1

0AA0

L 2.8 CONTROL 5SP43806 80 10 9.712 B 1

Combinación de arranque p/ motor:

Sitio Designació Clase de In Referencia Dispositivo Pmec Tipo Cant

n arranqu moto ánica de idad

e r [A] [kW] coordi

nación

/ Perfil

sobrec

argado

BTDP 2 MC-IA-RS Inversor 3.01 3RV23111D Interruptor 1.492 2 1

2.9A 2 C10 Relé de

3RB30161P sobrecarga

B0 Contactor
 40

Sitio Designació Clase de In Referencia Dispositivo Pmec Tipo Cant

n arranqu moto ánica de idad

e r [A] [kW] coordi

nación

/ Perfil

sobrec

argado

3RT20151A (sentido

P01 agujas reloj)

3RT20151A Contactor

P01 (sentido

contrario

agujas reloj)

Conexiones y Línea de distribución:

Cable / Conductor de media tensión:

Designaci Tipo / Punto de Ib [A] Materi Long Aislació Tipo de ftot θΔu [°C] /

ón Perfil inicio / Iz [A] al . [m] n instalaci θIkmax

[AWG/kc Punto de ón [°C] /

mil] llegada θIkmin [°C]

MT-C/L Ej: MV- Red 0.74 Cu 10 XLPE- Aire 1 20

2.1 90, MV- BTDP 7 cable 20

105 193 80

3x3
 41

Cable/Conductor de baja tensión:

Design Tipo / Punto Ib [A] Iz Materi Long. Aislac Tipo de Δu [%] θΔu [°C] / Núme

ación Perfil de inicio [A] al [m] ión instala / ∑ Δu θIkmax ro de

[AWG/k / Punto ción [%] [°C] / cable

cmil] de θIkmin s

llegada [°C]

C/L 2.9 ej: XHH, BTDP 2 3.012 CU 10 XLPE C 99.87 55 1

XHHW, Bomba 40 1 0.062 20

USE-2 s de 0.056 80

3x1x11/ agua

-/11

C/L 2.8 ej: XHH, BTDP 2 40 CU 50 XLPE C 98.78 55 1

XHHW, L 2.8 96 1 1.148 20

USE-2 1.142 80

3x1x50
 42

Sistema de electroductos/canalizaciones eléctricas:

Designación Tipo / Punto de Ib [A] Long. Con Grado ftot Δu θΔu [°C] Canti

Mater inicio / Punto In [A] [m] fig. de [%] / / θIkmax dad

ial de llegada de protec ∑ Δu [°C] /

elec ción [%] θIkmin

trod [°C]

ucto

s/ca

nali

zaci

one

eléc

trica

B 2.1 LDA1 BTDP 2.1 39.33 10 L1, IP54 1 99.9 55 1

625 BTDP 2 900 L2, 2 20

Al L3, 0.00 80

N, 6

PE 0

Cargas:
 43

Cargas estacionarias:

Designación Sitio Pn In [A] Un cos φ ai Cone Tipo de Cant

[kW] [V] xión carga idad

de

fase

L 2.8 Interior 23.69 40 380 0.9 1 L1- Inductiva 1

4 L2-

L3-N

Motores:

Designación Pmecá Un [V] ai Conexión Ia/In Clase Tipo de Ikre η Can

nica cos φ de fase de arranque tida

[kW] disparo d

In [A]

Bombas de 1.492 380 1 L1-L2-L3 5 Clase Combinaci 1 0.8 1

agua 3.012 0.84 10 ón de 96

arranque

para motor

ANEXO 2: Leyenda de símbolos

Símbolo [Unidad] Descripción

ai Factor de carga

cos φ Factor de potencia

 44

ftot Factor de reducción

Ia/In Relación de corrientes de inserción

Ib [A] Iz [A] Corriente de operación / Capacidad de carga permitida

Icu(fusible) [kA] Máxima capacidad de cortocircuito estimada - fusible

Icu [kA] Máxima capacidad de cortocircuito estimada acorde con IEC 60947-2

Icn [kA] Capacidad de cortocircuto estimada acorde con IEC 60898-1

Icu/Icn [kA] necesaria Necesario dispositivo en la instalación con capacidad de protección

Icw 1s [kA] Corriente nominal de tiempo corto admisible 1s

IΔn [mA] Corriente diferencial - Interruptor diferencial

Ik1máx Máxima corriente de cortocircuito monofásica

Ik1mín Mínima corriente de cortocircuito monofásica

Ik3máx Máxima corriente de cortocircuito trifásica

Ik3mín Mínima corriente de cortocircuito trifásica

Ik1D [kA] Corriente de cortocircuito monofásica ininterrumpida

Ik3D [kA] Corriente de cortocircuito trifásica ininterrumpida

Ikmáx/Ikmín Relación de máxima y mínima corriente de cortocircuito

Ikre Factor de recuperación de corriente de cortocircuito

In [A] Corriente nominal

P0 [kW] Pérdidas en vacío

 45

Pk [kW] Pérdidas en cortocircuito

Pmecánica [kW] mecánica. Potencia

Pn [kW] Potencia real nominal

R0 N [mΩ] Resistencia de Secuencia Homopolar entre Fase y Neutro

R0 PE(N) [mΩ] Resistencia de Secuencia homopolar entre fase y PE (N)

Relación de resistencia del sistema en secuencia de fase positiva /

R0/R1
secuencia de fase cero

R1 [%] Valor relativo de Resistencia fase-secuencia

R1 [mΩ] Resistencia de Secuencia de fase positiva

Sn [kVA] Potencia aparente nominal

ukr [%] Tensión de cortocircuito

Un [V] Tensión nominal

Uprim [kV] Sobretensión

Usec [V] Baja tensión

X0 N [mΩ] Reactancia de Secuencia Homopolar entre fase y neutro

X0 PE(N) [mΩ] Reactancia de Secuencia Homopolar entre fase y PE (N)

Relación de reactancias del sistema en secuencia de fase positiva /

X0/X1
secuencia de fase cero

X1 [mΩ] Reactancia de Secuencia de fase positiva

 46

xd'' [%] Reactancia subtransitoria

Z1 máx Máxima impedancia del sistema en secuencia de fase positiva

Z1 mín Mínima impedancia del sistema en secuencia de fase positiva

ZS Impedancia del circuito en falla

Zs máx Máxima impedancia del circuito en falla

Zs mín Mínima impedancia del circuito en falla

Tensión nominal relativa / Caída de tensión desde el principio al final del

Δu [%] / ∑ Δu [%] circuito / Caída de tensión acumulada desde bornes primario/ secundario

transformador hasta punto definido

Temperatura del conductor del cable de media tensión: / Temperatura del

θΔu [°C] / θIkmax [°C] /
conductor del cable de baja tensión:
θIkmin [°C]
para caída de tensión / con Ik máx. / para desconexión

η Grado de eficiencia

φ [°] Diferencia de fases

φ1 min/max [°] Fase ángulo Ik1 min/max

φ3 min/max [°] Fase ángulo at Ik3 min/max

ANEXO 3: Normativas para cálculos

Designación IEC HD EN DIN VDE

60364- 0100 –
Instalaciones eléctricas de baja tensión. * 384
1…6 100...710
 47

Cálculo de corrientes de cortocircuitos

60909 60909 0102
trifásicos/monofásicos en corriente alterna.

Cálculo de los efectos de las corrientes de

cortocircuitos 60865 60865 0103

Método de cálculo y definiciones

Aparellaje de maniobra y control de baja tensión -

60947-2 60947-2 0660 – 101
Interruptores

Aparellaje de maniobra y control de baja tensión en

61439 61439 0660 – 600
armarios/tableros eléctricos - Combinaciones

Efectos del aumento de la temperatura en

60890+
combinaciones de aparellaje de maniobra y control 528 S2 0660 – 507
C
mediante ensayo por extrapolación

Instalación Eléctrica de Edificios. Parte 5-52: 60364-

384 0298 – 4
Selección e instalación de equipamiento Eléctrico. 5-52

Instalaciones eléctrica de baja tensión – Selección de

equipamiento eléctrico – Parte 520: Sistemas de

0100-520
cableado – Suplemento 3:Capacidad de transporte de
Suplemento
corriente de los conductores en circuitos de
3
distribución trifásicos de corrientes de carga con

armón
 48

Accesorios eléctricos - Interruptores para protección

de sobrecorrientes para viviendas e instalaciones 60898-1 60898-1 0641 – 11

similares - Interruptores para corriente alterna

Aparellaje de maniobra y control de media tensión -

Combinación de seccionador con fusibles para 62271 62271 0671 – 105

corriente alterna

Instalaciones Eléctricas de Edificios, Selección e

60364- 60364-
instalación de equipamiento eléctrico, aislación, 0100-534
5-53 5-534
disparo y control

Instalaciones eléctricas de baja tensión. Protección

60364- 60364-
para seguridad, protección contra sobretensiones y 0100-443
4-44 4-443
disturbios electromagnéticos

62305-
Protección contra descargas atmosféricas-Parte 1…4 0185 – 1…4
1…4

Descargador para baja tensión. Descargador

61643-
conectado para sistemas en baja tensión; 0675-6-11
11
requerimientos y pruebas

Ensayos para cables eléctricos bajo condiciones de 60331- 50200 0472-814

fuego - integridad del circuito 11, 21 0482-200

 49

Compartamiento ante el fuego de materiales y 4102-12 :

componentes para la construcción - Parte 12: 1998-11

Mantenimiento de la integridad del circuito de

sistemas de cables eléctricos, requisitos y ensayos

Equipamento eléctrico de vehículos eléctricos - 61851 61851

Systema de carga de vehículos eléctricos

*) Consideraciones locales especiales diferentes a la norma IEC 60364-4-41:

2005 no están integradas y es necesario considerarlas!