ECUACION DEL CAMBIO DE ESTADO DEL CONDUCTOR

1. Consideraciones Previas

Tendido el conductor y en condiciones de servicio normal, éste se

ve sometido a los efectos de cambios en la presión de viento ó el peso

adicional de costra de hielo según donde se instale la línea.

Estas circunstancias hacen que el cable no mantenga

estática su ecuación, por tanto, su parámetro es cambiante.

Es posible deducir una ecuación que teniendo como dato un tiro

inicial (o esfuerzo inicial) en determinadas condiciones, calcular un tiro

final en otras condiciones.

Esta ecuación se denomina ECUACIÓN DE CAMBIO DE

ESTADO (ECE) del conductor.

2. ECE para Vanos de Nivel

En este caso la variación geométrica

Del conductor en su longitud, matemáticamente se expresa de la

siguiente manera (considerando que 𝐿1 ≅ 𝑎)

𝑇𝑜2 − 𝑇𝑜1
𝐿2 − 𝐿1 = 𝑎 ∝ (𝑡2 − 𝑡1 ) + 𝑎( )
𝑆. 𝐸.

Donde:

𝐿2 𝑦 𝐿1 : Longitud curvada de todo el vano del conductor en las condiciones

finales e iniciales (m).

∝: Coeficiente de dilatación lineal (1/°C)

𝑎: Longitud total del vano (m).

𝑡2 𝑦 𝑡1 : Temperatura final e inicial (°C).

𝑇𝑜2 𝑦 𝑇𝑜1: Tensión final e inicial en el punto más bajo (Kg).

𝑆: Sección del conductor (mm2)

𝐸: Módulo de elasticidad (Kg/mm2)

Antes de continuar, es necesario aclarar que en todos los modelos

matemáticos, los subíndices (1) significan condiciones iniciales, mientras

que los subíndices (2) significan condiciones finales.

Se tiene:

𝑎3 ∗ 𝑤 1 𝑣1
𝐿1 = 𝑎 +
24𝑇𝑜12

𝑎3 ∗ 𝑤 1 𝑣2
𝐿1 = 𝑎 +
24𝑇𝑜22

Donde:

𝑤 1 𝑣1 y 𝑤 1 𝑣2 Es la carga unitaria resultante en el conductor en las

condiciones finales e iniciales (Kg/m) y cuyo valor se determina partiendo

de la Figura:
 Fv

Pn

w
wr

Figura: Fenómenos naturales que afectan al conductor

𝑤𝑟 = √𝐹𝑣2 + (𝑤 + 𝑃ℎ )2

Donde:

𝑤𝑟 : Peso unitario resultante (Kg/m)

𝑤 : Peso unitario del conductor (Kg/m)

𝑃ℎ : Peso de la costra de hielo (Kg/m)

𝐹𝑣2 : Fuerza debido a la presión del viento (Kg/m)

La fuerza debido a la presión del viento y el peso de la costra de hielo de

pueden determinar con las siguientes relaciones.

𝐹𝑣 = 𝑃𝑣 (𝐷 + 2𝑒) ∗ 10−3

Donde:

𝑃𝑣 : Presión del viento sobre el conductor (Kg/m)

𝐷 : Diámetro exterior del conductor (mm)

𝑒 : Espesor de la costra de hielo (mm). Este valor está determinado en la

tabla 250 CNE.

Observando la Figura, se determina que la superficie neta de la costra de

hielo (𝑆ℎ ) es:

𝜋 𝜋
𝑆ℎ = ( ) (𝐷 + 2𝑒)2 − ( ) 𝐷2
4 4

Resolviendo y simplificando la ecuación, se tiene:

𝑆ℎ = (𝜋𝑒)(𝐷 + 𝑒) 𝑚𝑚2

Convirtiendo a m2 sería:

𝑆ℎ = (𝜋 ∗ 10−6 ∗ 𝑒)(𝐷 + 𝑒) 𝑚2

Entonces, considerando que la densidad del hielo es de 913 Kg/m3

(según el CNE) el peso de la costra de hielo (𝑃ℎ ), sería:

𝑃ℎ = (913𝜋 ∗ 10−6 ∗ 𝑒)(𝐷 + 𝑒)

Operando y redondeando la ecuación, finamente se tiene:

𝑃ℎ = (0.0029𝑒)(𝐷 + 𝑒)

El valor de la presión del viento sobre el conductor se puede calcular con

la siguiente relación:

𝑄𝑉2
𝑃𝑣 =
2𝑔

Donde:

𝑃𝑣 : Presión del viento sobre el conductor (kg/m2)

𝑄 : Peso del aire (1,225 Kg/m2)

𝑔 : Aceleración de la gravedad (9,81 m/s2)

v : Velocidad del viento (m/s)

Reemplazando valores, en la ecuación, se tiene:

𝑃𝑣 = 0,062 436 29 𝑣 2

Convirtiendo la velocidad del viento (v) de (m/s) a (km/h) y considerando

que 1 m/s=3,6 Km/h; se tiene:

0.062 436 29 𝑣 2
𝑃𝑣 =
3,62

Finalmente, simplificando se tiene:

𝑃𝑣 = 0,004 82 𝑣 2

Donde:

: Velocidad del viento (Km/h)

Pero según el CNE, la presión del viento sobre el conductor se determina

con la siguiente relación:

𝑃𝑣 = 𝐾 𝑣 2 𝑆𝑓

Donde:

𝑃𝑣 : Presión del viento sobre el conductor (N/m2)

𝑣 : Velocidad del viento (m/s)

𝐾 : Constante de presión que depende de la altitud

K= 0,613 para las elevaciones hasta la 3 000 msnm

K= 0,455 para las elevaciones mayores de 3 000 msnm

𝑆𝑓 : Factor de forma (se considera sd=1)

Entonces las fórmulas que recomienda el CNE, son:

Para elevaciones menores a 3 000 msnm, convirtiendo la velocidad del

viento en (Km/h) y la presión del viento sobre el conductor en Kg/m2),

sería:

𝑃𝑣 = 0,004 82 𝑣 2

Mientras que para elevaciones mayores a 3 000 msnm, convirtiendo la

velocidad del viento en (km/h) y la presión del viento sobre el conductor

en (kg/m2), sería:

0,455
𝑃𝑣 = 𝑣 2 = 0,003 58 𝑣 2
3,62 ∗ 9,81

Si observamos la ecuación es parecida a una ecuación anterior, lo que

indica que el CNE, acepta dicha ecuación, pero para altitudes menores de

3 000 msnm.

Por otro lado, con respecto a la velocidad del viento, espesor del hielo y

temperatura; en el CNE, se ha identificado las zonas de carga en todo el

territorio peruano, el cual se ha dividido en tres zonas:

 Zona A – zona de carga ligera

 Zona B – zona de carga regular

 Zona C – zona de carga fuerte

Además, de acuerdo a la altitud sobre el nivel del mar, se ha dividido en

cuatro áreas:

 Área 0 altitud menor a 3 000 msnm.

 Área 1 altitud de 3 001 a 4 000 msnm.

 Área 2 altitud de 4 001 a 4 500 msnm.

  Área 3 altitud mayor a 3 500 msnm.

Las zonas de carga A, B y C; incluyen cargas de viento y las áreas 1,1 2

y; incluyen cargas de viento incrementadas con las cargas de hielo. En el

cuadro se muestra las velocidades del viento y las temperaturas que se

deben usar según la zona de carga.

A B C
Zonas de Carga
Ligera Ligera Ligera

19,5 m/s 22,2 m/s 25 m/s

Velocidad horizontal del viento
(70 Km/h) (80 km/h) (90 km/h)
Temperatura 20 °C 15 °C 10 °C
También, en el cuadro se muestra la velocidad del viento, el espesor

radial del hielo y la temperatura que se deben utilizar para calcular las

cargas según las zonas de carga, según el cual se presenta tres casos:

caso de sólo viento, caso de sólo hielo o caso de combinación de viento

y hielo.

Zonas de carga Área 0 Área 1 Área 2 Área 3

Caso de sólo viento 26m/s 29m/s 31,5 m/s 33,5 m/s
Velocidad horizontal del viento (94 km/h) (104 km/h) (113 km/h) (120 km/h)
Temperatura 10°C 5°C 0°C -5°C
Caso de sólo hielo
no hay 6 mm 25 mm 50 mm
Espesor radial del hielo
0°C 0°C 0°C 0°C
Temperatura
Caso de hielo y viento No hay No hay 12 mm 25 mm
Espesor radial del hielo 14 m/s 14 m/s 15,5 m/s 17 m/s
Velocidad horizontal del viento (50 km/h) (50 km/h) (56 km/h) (61 km/h)
Temperatura 5°C 5°C -5°C -10°C

Por otro lado, continuado con el desarrollo de los modelos matemáticos;

reemplazando las ecuaciones, se tiene:

2 2
𝑎2 . 𝑤𝑟2 𝑎2 . 𝑤𝑟1 𝑇𝑜2 − 𝑇𝑜1
𝑎+ + 𝑎 − = 𝑎 ∝ (𝑡2 − 𝑡1 ) + 𝑎( )
24𝑇𝑜22 24𝑇𝑜12 𝑆. 𝐸.

Simplificando la ecuación, se tiene:

2 2
𝑎2 . 𝑤𝑟2 𝑎2 . 𝑤𝑟1 𝑇𝑜2 − 𝑇𝑜1
2 − 2 =∝
(𝑡2 − 𝑡1 ) + ( )
24𝑇𝑜2 24𝑇𝑜1 𝑆. 𝐸.

De la ecuación, usualmente sólo es desconocido 𝑇𝑜2 , por lo que

tratando de despejar divo valor se obtiene:

2 2
𝑎2 . 𝑤𝑟1 . 𝑆. 𝐸 𝑎2 . 𝑤𝑟2 . 𝑆. 𝐸
𝑇𝑜22 [𝑇𝑜2 + 𝑆. 𝐸. ∝. (𝑡2 − 𝑡1 ) + − 𝑇𝑜1 ] =
24𝑇𝑜12 24
 Si emplemos el concepto de los esfuerzos unitarios en el conductor
𝑇𝑜1 𝑇𝑜2
definidos por 𝜎01 = y 𝜎02 = que están expresados en /kg/mm2) y
𝑆 𝑆

lo reemplazamos en la ecuación, se tiene:

2 2
2
𝑎2 . 𝑤𝑟1 .𝐸 𝑎2 . 𝑤𝑟2 . 𝑆. 𝐸
𝜎02 [𝜎01 +∝. 𝐸(𝑡2 − 𝑡1 ) + − 𝜎01 ] =
24𝑆 2 𝜎01 2 24𝑆 2

A esta ecuación se conoce como la ECUACIÓN DE CAMBIO DE

ESTADO para vanos nivelados donde el único valor desconocido es el

esfuerzo unitario en las condiciones finales (𝜎02 ) y cuyo valor se obtiene

resolviendo la ecuación cúbica corresponiente.

3. Ecuación de cambio de estado para vanos desnivelados por el método

de la parábola.

Cuando se presenta el caso de vanos desnivelados se puede demostrar

que introduciendo el valor de 𝑐𝑜𝑠𝛿 en la ecuación anterior se obtiene la

ecuación de cambio de estado para vanos desnivelados, entonces la

ecuación quedaría como:

2 2
2
𝑎2 . 𝑤𝑟1 . 𝐸. 𝑐𝑜𝑠 3 𝛿 𝑎2 . 𝑤𝑟2 . 𝐸𝑐𝑜𝑠 3 𝛿
𝜎02 [𝜎01 +∝. 𝐸(𝑡2 − 𝑡1 )𝑐𝑜𝑠𝛿 + − 𝜎01 ] =
24𝑆 2 𝜎01 2 24𝑆 2

Observando la Figura, que corresponde a vanos desnivelados, se puede

determinar el valor de δ de la siguiente manera:

ℎ
𝛿 = 𝑡𝑎𝑛−1 ( )
𝑎

Donde:

: Es el desnivel de una estructura con la otra (m).

: Vano (m)
 Para resolver la ecuación con más facilidad se debe decir que todos los

datos son conocidos, excepto el esfuerzo unitario en la condición final

(𝜎02 ), entonces para facilitar su solución se puede hacer los siguientes

cambios:

𝑋 = 𝜎02

2
𝑎2 . 𝑤𝑟1 . 𝐸. 𝑐𝑜𝑠 3 𝛿
𝐴 =∝. 𝐸(𝑡2 − 𝑡1 )𝑐𝑜𝑠𝛿 + − 𝜎01
24𝑆 2 𝜎01 2

2
𝑎2 . 𝑤𝑟2 . 𝐸𝑐𝑜𝑠 3 𝛿
𝐵=
24𝑆 2

Entonces, reemplazando las ecuaciones, se tiene:

𝑥 2 (𝑥 + 𝐴) = 𝐵

Resolviendo se tiene:

𝑥 3 + 𝐴𝑥 2 − 𝐵 = 0

Es una ecuación cúbica y por tanto tendrá tres soluciones para “x”, pero

solo una de ellas será un valor real positivo; para dicha solución existen

varios métodos que se estudia en el curso de métodos numéricos, pero

el más recomendable es el método de Newton-Raphson, que se

recomienda al lector buscar en internet donde hay mucha información al

respecto.

4. Ecuación de cambio de estado para vanos desnivelados por el método

de la catenaria.

El análisis con respecto a la ecuación de cambio de estado para vanos

novelados y desnivelados, corresponden al método de la parábola, ya

que todos los modelos matemáticos utilizados son de dicho método.

Ahora, haremos un análisis, pero por el método de la catenaria; psrs ello

partiremos de la siguiente relación:

 𝑇𝑜2 − 𝑇𝑜1
𝐿2 − 𝐿1 =∝ (𝑡2 − 𝑡1 )𝐿1 + ( )𝐿1
𝑆. 𝐸.

La ecuación representa el cambio de longituf, que experimenta un

conductor soportado entre dos puntos, debido a los cambios de la

temperatura y cambios en las condiciones climáticas de la zona donde

se encuentra instalada la línea de transmisión.

𝑇𝑜1 𝑇𝑜2
Considerando los esfuerzos unitarios 𝜎01 = y 𝜎02 = y factorizando
𝑆 𝑆

la ecuación, se tiene:
𝜎02 𝜎01
𝐿2 − 𝐿1 = 𝐿1 [∝ (𝑡2 − 𝑡1 ) + ( − )]
𝑆 𝑆

Simplificando, se tiene:

𝐿2 𝜎02 𝜎01
− 1 =∝ (𝑡2 − 𝑡1 ) + ( − )
𝐿1 𝑆 𝑆

En la ecuación, los valores de 𝐿2 y 𝜎02 son desconocidos, por lo tanto lo

pasaremos al primer miembro, y el resto al segundo miembro, entonces

quedaría como:

𝐿2 𝜎02 𝜎01
− =∝ (𝑡2 − 𝑡1 ) − +1
𝐿1 𝐸 𝐸

Multiplicando por E a toda la ecuación y simplificando, se tiene:

𝐿2 𝐸
− 𝜎02 = 𝐸 ∝ (𝑡2 − 𝑡1 ) − 𝜎01 + E
𝐿1

Factorizando E se tiene:

𝐿2 𝐸
− 𝜎02 = 𝐸[∝ (𝑡2 − 𝑡1 ) + 1] − 𝜎01
𝐿1

Para que la ecuación de cambio de estado no sea muy larga, hacemos

los siguientes cambios de variable:

𝐵 = 𝐸[∝ (𝑡2 − 𝑡1 ) + 1] − 𝜎01

 𝐸
𝐴=
𝐿1

Reemplazando las ecuaciones, se tiene:

𝐴𝐿2 − 𝜎02 = 𝐵

Ahora, desarrollamos y calcularemos los valores de 𝐿1 , para ello

partiremos de la ecuación, que corresponde a la longitud curvada de un

conductor para vanos desnivelados, por lo tanto:

𝑎𝑟 2
𝐿1 = √[2𝑐1 . 𝑠𝑒𝑛ℎ ( )] + ℎ2
2𝑐1

Pero según la ecuación, el parámetro de la catenaria ( ), sería:

𝑇01 𝜎01 𝑆
𝑐1 = =
𝑤𝑟1 𝑤𝑟1

Reemplazando la ecuación, se tiene:

𝜎01 𝑆 𝑤𝑟1 𝑎𝑟 2
𝐿1 = √[2 . 𝑠𝑒𝑛ℎ ( )] + ℎ2
𝑤𝑟1 2𝜎01 𝑆

En la ecuación, todas las variables son conocidas y por tanto el valor de

𝐿1 se puede hallar sin ningún problema; en forma análoga, se pude

determinar 𝐿2 , sólo cambiando los subíndices ( 1 ) por subíndices ( 2 ).

𝜎02 𝑆 𝑤𝑟2 𝑎𝑟 2
𝐿2 = √[2 . 𝑠𝑒𝑛ℎ ( )] + ℎ2
𝑤𝑟2 2𝜎02 𝑆

Reemplazando la ecuación, se tiene:

𝐸
𝐴=
𝜎01 𝑆 𝑤𝑟1 𝑎𝑟 2 2
√[2
𝑤𝑟1 . 𝑠𝑒𝑛ℎ (2𝜎01 𝑆 )] + ℎ

Reemplazando la ecuación, se tiene:

 𝜎02 𝑆 𝑤𝑟2 𝑎𝑟 2
𝐴 ∗ √[2 . 𝑠𝑒𝑛ℎ ( )] + ℎ2 − 𝜎02 = 𝐵
𝑤𝑟2 2𝜎02 𝑆

Si hacemos:

2𝑆
𝑘=
𝑤𝑟2

Finalmente, reemplazando la ecuación y ordenando adecuadamente, se

tiene:

𝑎𝑟 2
𝐴 ∗ √[𝑘𝜎02 . 𝑠𝑒𝑛ℎ ( )] + ℎ2 − 𝜎02 − 𝐵 = 0
𝑘𝜎02

La ecuación corresponde a la ecuación de cambio de estado por el

método de la catenaria para vanos desnivelados, pero también sirve

para vanos nivelados, pues en ese caso el desnivel sería h=0. En esta

ecuación la única incógnita sería el esfuerzo unitario (𝜎02 ) que se tendría

que resolver mediante métodos numéricos, y que para solucionar e

sugiere el método de iteración de punto fijo.