Campo Elec en Aisladores
Campo Elec en Aisladores
Campo Elec en Aisladores
FACULTAD DE INGENIERÍA
T E S I S
QUE PARA OPTAR POR EL GRADO DE:
MAESTRO EN INGENIERÍA
INGENIERIA ELÉCTRICA – SISTEMAS ELÉCTRICOS DE POTENCIA
P R E S E N T A:
TUTOR:
2012
JURADO ASIGNADO:
TUTOR DE TESIS:
DR. RAFAEL ESCARELA PÉREZ
______________________________
FIRMA
AGRADECIMIENTOS
A la Universidad Nacional Autónoma de México por todas las oportunidades y apoyo que
he recibido por parte de ella.
A los ingenieros del área de líneas de transmisión del Grupo de Apoyo del Instituto de
Investigaciones Eléctricas para la Coordinación de Proyectos de Transmisión y Transformación de
la Comisión Federal de Electricidad, por su contribución en mi formación profesional y personal.
DEDICATORIA
A mi familia.
Rolling Stones
RESUMEN
En este trabajo se realiza el análisis del campo eléctrico en cadenas de aisladores estándar
de vidrio templado utilizadas en líneas eléctricas de transmisión de Sistemas Eléctricos de Potencia
aplicando el Método de los Elementos Finitos. El análisis se realiza para evaluar el
comportamiento del campo eléctrico en las cadenas de aisladores bajo los efectos de
desplazamientos transversales de las mismas, respecto a las estructuras de transmisión, y que son
producidos por la presencia de viento. Se consideran las configuraciones geométricas de
estructuras de transmisión empleadas en México.
i
ABSTRACT
This work is about the electric field analysis on standard toughed glass insulator strings
used in electric transmission lines of electrical power systems applying the Finite Elements
Method. Analysis is done in order to evaluate the electric field behavior in insulator string under
the effects of transversal displacements, related to transmission structures, and which are
produced by wind presence. Geometrical configurations of transmission structures used in Mexico
are considered.
Such as part of the work are considered electromagnetic theory topics involved, also
procedures and standards that define the behavior of insulators string displacement in presence of
wind mechanical loads on transmission lines.
Tridimensional numerical models are developed for structure configurations and insulator
strings of 115, 230, and 400 kV voltage levels, used by Comisión Federal de Electricidad in Mexican
transmission lines, and with different swing angles. This way, results of electric field and potential
distribution along strings were obtained, showing maximum electrical efforts in each study case
developed.
ii
CONTENIDO
RESUMEN i
ABSTRACT ii
CONTENIDO iii
FIGURAS v
TABLAS viii
NOMENCLATURA ix
1 INTRODUCCIÓN 1
1.1 Introducción 1
1.2 Antecedentes 3
1.2.1 Aportaciones realizadas al análisis del campo eléctrico y distribución de 4
potencial en cadenas de aisladores con métodos numéricos
1.3 Objetivo de la tesis 7
1.4 Alcance de la tesis 7
1.5 Estructura de la tesis 8
2 EL CAMPO ELÉCTRICO EN CADENAS DE AISLADORES 9
2.1 Introducción 9
2.2 Estructura y diseño de cadenas de aisladores en líneas de transmisión 10
2.3 Campo eléctrico y distribución de potencial en cadenas de aisladores 12
2.3.1 Métodos experimentales para el análisis de campo eléctrico y 14
distribución de potencial en aisladores
2.3.2 Métodos numéricos para el análisis de campo eléctrico y distribución de 15
potencial en aisladores
2.4 Fundamentos de teoría electromagnética 18
2.5 El método de los elementos finitos 20
2.6 Programas de cómputo aplicados al análisis de campo eléctrico por métodos 21
numéricos
3 ÁNGULO DE BALANCEO EN CADENAS DE AISLADORES DE LÍNEAS DE TRANSMISIÓN 23
3.1 Introducción 23
3.2 Normas y especificaciones 23
3.2.1 NOM-001-SEDE-2005 24
3.2.2 CIGRE 348 y 72 25
3.2.2.1 Libramientos 33
3.2.3 Comisión Federal de Electricidad 35
4 CASOS DE ESTUDIO EN EL ANÁLISIS DEL CAMPO ELÉCTRICO EN CADENAS DE 36
AISLADORES DE LÍNEAS DE TRANSMISIÓN
4.1 Introducción 36
iii
4.2 Casos de estudio 37
4.2.1 Características 38
4.3 Descripción geométrica de las líneas de transmisión de 115 kV 39
4.3.1 Campo eléctrico y distribución de potencial en cadena de aisladores en 41
disposición vertical para 115 kV
4.3.2 Campo eléctrico y distribución de potencial en cadena de aisladores con 49
ángulo de balanceo para 115 kV
4.4 Descripción geométrica de las líneas de transmisión de 230 kV 54
4.4.1 Campo eléctrico y distribución de potencial en cadena de aisladores en 56
disposición vertical para 230 kV
4.4.2 Campo eléctrico y distribución de potencial en cadena de aisladores con 64
ángulo de balanceo para 230 kV
4.5 Descripción geométrica de las líneas de transmisión de 400 kV 69
4.5.1 Campo eléctrico y distribución de potencial en cadena de aisladores en 70
disposición vertical para 400 kV
4.6 Modelo geométrico de cadena de aisladores con cambios en longitud del cable 74
conductor
4.7 Modelo geométrico de cadena de aisladores con cambios en dimensión de 76
fronteras
4.8 Resultados y análisis 79
5 CONCLUSIONES 81
5.1 Conclusiones 81
5.2 Trabajos futuros 82
BIBLIOGRAFÍA 84
APENDICE A 90
iv
FIGURAS
vi
balanceo para torre 2B2.
Figura 4.35. Líneas de flujo de campo eléctrico en cadena de aisladores de 230 kV con ángulo 68
de balanceo para torre 2B2.
Figura 4.36. Campo eléctrico en cadena de aisladores de 230 kV con ángulo de balanceo para 68
torre 2B2.
Figura 4.37. Dimensiones generales de estructura de 400 kV de la CFE. 70
Figura 4.38. Distribución de potencial en cadena de aisladores de 400 kV en configuración 71
“V” para torre EA4A22MA.
Figura 4.39. Distribución de potencial en cadena de aisladores de 400 kV en configuración 72
“V” para torre EA4A22MA.
Figura 4.40. Líneas de flujo de campo eléctrico en cadena de aisladores de 400 kV en 72
configuración “V” para torre EA4A22MA.
Figura 4.41. Campo eléctrico en cadena de aisladores de 400 kV en configuración “V” para 74
torre EA4A22MA.
Figura 4.42. Comparativa de la distribución de potencial en cadena de aisladores de 115 kV 75
para poste 1210SMP con longitudes de conductor diferentes.
Figura 4.43. Comparativa del campo eléctrico en cadena de aisladores de 115 kV para poste 75
1210SMP con longitudes de conductor diferentes.
Figura 4.44. Modelo geométrico establecido de sección de poste 1210SMP de 115 kV con 77
fronteras externas de 20 metros (distribución de potencial).
Figura 4.45. Campo eléctrico en modelo establecido de cadena de aisladores de 115 kV con 77
cambios en dimensiones de fronteras externas.
vii
TABLAS
viii
NOMENCLATURA
V [V] Voltaje
Número total de aisladores en la cadena o conductores agrupados
n
en haz
k = Cg/C Relación de capacitancias
v [m³] Volumen
S [m²] Superficie
V Gradiente de V
̅ [C/m²] Densidad de flujo eléctrico
[C/m³] Densidad volumétrica de carga
[C] Flujo eléctrico
̅ Divergencia de ̅
[F/m] Permitividad
̅ [V/m] Campo eléctrico
[S/m] Conductividad
[°] Ángulo de balanceo de aislador
[kg/m³] Densidad del aire
[m/s] Velocidad del viento
Coeficiente de arrastre de los conductores (1,0)
d [m] Diámetro del conductor
[m] Claro de viento
[N] Carga de viento en la cadena de aisladores
[N] Peso del conductor
[N] Peso del conjunto de aislamiento
[N] Carga de viento en conductores
Presión del viento correspondiente a la velocidad promedio del
( ) [N/m²]
viento
Factor de ráfaga del conductor
Factor de claro
[m/s] Velocidad de viento de referencia
[m] Altura de referencia de la velocidad del viento
Factor de ajuste acorde al tipo de terreno
ix
kW Parámetro
K Matriz de orden N por N
i,j,m Números enteros positivos
UR [kV] Voltaje nominal
US [kV] Voltaje máximo del sistema
Uff [kV] Voltaje resistido por frente de onda rápido (descarga)
Umax-sf [kV] Voltaje resistido por frente de onda lento (maniobra)
Dpe [m] Distancia mínima fase a fase
Dpp [m] Distancia mínima fase a tierra
ABREVIACIONES DESCRIPCIÓN
x
1
INTRODUCCIÓN
1.1. Introducción
1
Las cadenas de aisladores operan bajo esfuerzos mecánicos, eléctricos y ambientales
constantes que deben ser soportados para que la continuidad y seguridad del sistema de
transmisión no se vea afectada. Dentro de los fenómenos físicos que se presentan en las cadenas
de aisladores se pueden mencionar:
Las corrientes que se presentan a lo largo de los aisladores no son uniformes debido a: la
disposición de sus partes; el espaciamiento entre unidades; condiciones ambientales y la
resistencia eléctrica de las superficies de los aisladores, de tal forma que la distribución capacitiva
y consecuentemente las corrientes y voltajes son mayores en algunas partes de una cadena de
aisladores, lo que puede contribuir a un daño en algún punto particular de un aislador. Por lo
anterior, es importante que se distribuyan de manera uniforme (en la medida de lo posible)los
esfuerzos eléctricos a lo largo de la cadena, con la finalidad de que alguna unidad no presente un
sobreesfuerzo continuo y llegue a fallar con el paso del tiempo. Las cadenas de aisladores que
operan con niveles de voltaje mayores a 230 kV son más susceptibles de presentar una
2
distribución no uniforme de voltaje importante, afectándose así a los aisladores, de tal forma que
por medio de cambios en su forma, materiales y la adición de ciertos dispositivos como los anillos
equipotenciales se establece una mejor distribución de los esfuerzos eléctricos en las cadenas de
aisladores, previniendo un envejecimiento y falla prematura [Ilhan, 2007].
Con la finalidad de lograr un diseño tanto técnica como económicamente adecuado de los
aislamientos es necesario contar con información asociada a su comportamiento eléctrico.
Consecuentemente, las mediciones y los cálculos de los esfuerzos eléctricos en las cadenas de
aisladores representan las bases de un correcto diseño, siendo el campo eléctrico y la distribución
del potencial eléctrico escalar dos de los parámetros más importantes a ser considerados.
1.2. Antecedentes
3
sintéticas) mediante la aplicación de métodos numéricos permite el diseño del aislador, así como
estudiar el dimensionamiento de las estructuras de transmisión (torres y postes) y la configuración
de las líneas de transmisión, influyendo directamente en la optimización de costos y mejoras en el
desempeño de los sistemas de transmisión [Zhao, 2000].
Las aportaciones y trabajos que se han efectuado para el análisis del campo eléctrico en
los aislamientos de líneas de transmisión, ya sean de vidrio templado, porcelana ó sintéticas
aplicando métodos numéricos se describen a continuación.
4
equipotencial se aproxima a su parte aterrizada y e) el diámetro del tubo del anillo equipotencial
presenta poca influencia en el campo eléctrico a lo largo del aislador, pero es importante para el
control del efecto corona.
5
Izgi define un modelo simplificado capacitivo-resistivo equivalente de las cadenas de
aisladores en el que evalúa distribución de potencial y la eficiencia de la cadena bajo condiciones
de viento y contaminación, observando una mejora en la distribución de potencial conforme la
acción del viento aleja la cadena de aisladores del cuerpo de la estructura y viceversa conforme se
acerca la cadena a la estructura [Izgi, 2008].
Subba y Kumar analizan la distribución del campo para aisladores cerámicos empleando el
método numérico de simulación de carga superficial (CSM, por sus siglas en inglés) [Subba, 2011].
6
Desarrollan un electrodo conformado de acero inoxidable y aluminio que se coloca en la región del
perno ó pin del aislador para reducir el campo eléctrico máximo.
7
comportamiento ante diversos fenómenos como el causado por el desplazamiento transversal de
las cadenas ante la presión del viento para determinados valores de ángulo de balanceo.
8
2
EL CAMPO ELÉCTRICO EN CADENAS DE AISLADORES
2.1. Introducción
9
2.2. Estructura y diseño de cadenas de aisladores en líneas de transmisión
10
“Coordinación de Aislamiento – Parte 2 Guía de Aplicación” [NMX, 2004] al igual que
especificaciones de la Comisión Federal de Electricidad [CFE, 1991 & 1995],donde se establecen
los requisitos para las separaciones mínimas de los conductores de líneas aéreas con sus soportes
y conductores adyacentes. La longitud de las cadenas de aisladores varía acorde a ciertas
características (nivel de tensión, grado de contaminación, altitud sobre nivel del mar, distancias
de paso y de fuga del aislador empleado) que deben ser establecidas para obtener valores
seguros de aislamiento para las condiciones de operación a las que son sometidas.
Acotaciones en
Manguito mm sin escala
de zinc
146
18
254
11
2.3. Campo eléctrico y distribución de potencial en cadenas de aisladores
V1
Cg C
I1
Cuerpo recto de la i1
estructura a potencial Cg C V2
de tierra. I2 V
i2
Cg C V3
I3
i3
Cg C V4
I4
i4
Conductor
Figura 2.4. Circuito equivalente de cadena de aisladores [Still, 1927].
12
Los voltajes y corrientes obtenidos de una cadena de aisladores presentan un
comportamiento no lineal de tal forma que las unidades más cercanas a la parte conductora ó
energizada se encuentran más esforzadas eléctricamente que las alejadas de la misma. Una forma
tradicionalmente usada para determinar el grado de esfuerzo eléctrico presente en la cadena es
mediante el concepto de “eficiencia de la cadena” que se encuentra definido (Figura 2.4) mediante
la expresión [Still, 1927]:
( ) (2.1)
Distribución 100
de Voltaje
[%]
80
60
k↓↓
k↓
40
k
20
Número de
Aislador
0 1 2 3 4 5
Figura 2.5. Distribución del voltaje entre las unidades de una cadena de aisladores [Still, 1927].
13
El modelo del circuito equivalente planteado en la Figura 2.4 presenta muchas
simplificaciones como por ejemplo, no considera una resistencia para la parte dieléctrica del
aislador (vidrio/porcelana y el cemento) ó especifica una capacitancia constante tanto mutua
como a tierra. Otros factores no considerados por el modelo son la contaminación presente en los
aisladores ó la posición de la cadena cuando es desplazada por viento. Por lo anterior, los
resultados que arroja este modelo analítico simplificado del comportamiento eléctrico de las
cadenas no son correctos y se presenta la necesidad de desarrollar otros modelos y métodos que
definan el comportamiento eléctrico de la cadena de forma más real.
Los métodos experimentales y los análisis por métodos numéricos empleados para el
estudio de la distribución del campo eléctrico y potencial en las cadenas de aisladores presentan
ventajas ya que permiten considerar los fenómenos y características que los analíticos no pueden
debido a la complejidad de su definición dentro del modelo.
14
Un estudio comparativo publicado por R. Parraud [Que, 2002] de diferentes métodos de
medición del campo eléctrico señala dentro de los sensores empleados: a) medidores de
potenciales AC, b) sensores esféricos capacitivos con enlaces de datos ópticos y c) sensores
electro-ópticos de cuarzo cúbicos. Parraud concluyó que los primeros son poco confiables debido
a la influencia de los sensores de medición en la distribución del potencial, los segundos
presentan buenas correlaciones en la región de los electrodos pero diferencias significativas cerca
de los electrodos conectados a tierra, mientras que los últimos son una buena elección respecto a
los valores medidos y los calculados del campo eléctrico, con la excepción de la cercanía a los
electrodos debido a la distorsión del campo por la ubicación de los sensores de medición.
15
Los métodos de análisis del campo eléctrico se dividen en dos categorías [Zhou, 1993;
Que, 2002]: los métodos de frontera, de momentos ó distribución de fuente -fronteras abiertas-
(CSM y BEM) y los métodos de dominio espacial, de dominio del tiempo y aproximación de campo
-fronteras cerradas- (FEM y FDM).
El concepto básico del CSM es remplazar la distribución de carga de los conductores y las
cargas de polarización en las interfaces dieléctricas por una gran cantidad de cargas discretas
ficticias. Las magnitudes de estas cargas deben ser calculadas tal que su efecto integrado satisfaga
las condiciones de frontera. El potencial debido a cargas superficiales desconocidas puede ser
aproximado por tres formas de arreglos de cargas ficticias concentradas: línea, anillo y cargas
puntuales. Estas cargas pueden ser ubicadas en posiciones apropiadas, que usualmente se
encuentran dentro de la superficie del conductor ó electrodo, fuera del espacio en el cual el
campo eléctrico será determinado [Khalifa, 1990]. Una combinación adecuada de estas tres
formas de cargas puede efectuarse para simular casi cualquier sistema práctico de electrodos. El
CSM puede emplearse para resolver problemas de fronteras abiertas y es aplicado a problemas de
campo eléctrico tridimensionales sin simetría axisimétrica. Sin embargo, un problema es la
ubicación subjetiva de las cargas de simulación. También, es difícil o imposible calcular el campo
eléctrico cerca de electrodos muy delgados debido a que las cargas ficticias que aproximan el
campo deben estar dentro de los electrodos.
16
El FEM es un método numérico de solución de las ecuaciones de Maxwell en su forma
diferencial [Que, 2002]. La característica básica del método FEM es dividir el espacio del problema
entero, incluyendo las regiones aledañas, en un número de subregiones no separadas llamadas
“elementos finitos” (ver Figura 2.6 para la discretización en elementos finitos de una cadena de
aisladores). Este proceso de discretización se llama mallado. Dentro de cada elemento
geométrico, la variable desconocida es aproximada comúnmente por polinomios con valores
nodales como coeficientes. El potencial eléctrico dentro de cada elemento es una interpolación
de los potenciales en sus vértices. Mediante el uso de métodos de aproximaciones residuales, las
ecuaciones diferenciales parciales son reducidas a una ecuación matricial definida positiva,
simétrica y dispersa. Debido a que el tamaño y forma de los elementos es arbitrario, el método
resulta flexible y funciona para problemas de geometría complicada. El FEM es más efectivo para
problemas con fronteras cerradas. Si el problema es muy grande un número mayor de elementos
finitos es necesario y los cálculos son mayores [Haddad, 2004].
17
2.4. Fundamentos de teoría electromagnética
Dentro de las ecuaciones de Maxwell se encuentra la que expresa la ley de Gauss para los
campos eléctricos, la que establece que el flujo eléctrico total que pasa a través de una
superficie cerrada es igual a la carga que encierra la superficie:
∮ ̅ ̅ ∮ (2.2)
∮ ̅ ̅
̅
La ecuación (2.2) es la primera de las cuatro ecuaciones de Maxwell que aplica para campos
electrostáticos, señalando que la cantidad de flujo eléctrico por unidad de volumen dejando una
unidad de volumen es igual a la densidad de carga en el volumen.
̅ ̅ (2.4)
donde es la permitividad. En el caso de campos eléctricos estáticos se puede demostrar que las
ecuaciones de Maxwell llevan a una ecuación diferencial parcial conocida como la “ecuación de
Poisson”, con la cual se puede determinar de manera única al campo, si se conoce las
18
distribuciones de carga y las condiciones de frontera de la región analizada. Ante la ausencia de
fuentes de carga la ecuación de Poisson se simplifica a la ecuación de Laplace.
Para la determinación del campo eléctrico en los aislamientos de sistemas de alta tensión
en corriente alterna, los métodos de solución normalmente consideran que las corrientes en los
aislamientos son despreciables, ya que la permitividad es mucho mayor a la conductividad
[Asenjo, 1997]. El campo electrostático en un sistema puede determinarse mediante la
distribución del potencial dentro de una región. La ley de Gauss para campos electrostáticos
establece que (ver ecuaciones (2.3) y (2.4)):
̅ (2.5)
El valor de la magnitud de intensidad del campo eléctrico, que está dado por el máximo
valor de la relación de cambio del potencial con la distancia, es normal a las superficies
equipotenciales generadas por el potencial. Debido a que el campo eléctrico puede expresarse en
términos del potencial escalar mediante su gradiente [Hyatt, 2001]:
̅ (2.6)
(2.7)
ó,
(2.8)
si la permitividad del sistema es constante. Para sistemas más generales (como en cadenas de
aisladores) la permitividad no es constante y es necesario resolver la ecuación (2.7). La ecuación
(2.7) se conoce como la ecuación de Poisson y aplica para medios no homogéneos donde la
permitividad no es constante. En el caso de que , es decir cuando la densidad de carga es
nula, se tiene que [Hyatt, 2001]:
(2.9)
19
ó,
(2.10)
El método de los elementos finitos es una técnica numérica utilizada para resolver
ecuaciones diferenciales parciales y ordinarias. El FEM transforma problemas de región continua
en una subdivisión de elementos discretos (normalmente triangulares en modelos
bidimensionales) representando a la región de estudio. La solución que se obtiene del problema
es a través del potencial aproximado de forma polinomial en los nodos correspondientes a los
vértices de cada elemento y que puede presentar distintas propiedades [Haddad, 2004].
∑ (2.11)
j
Figura 2.7. Elemento triangular.
20
Aplicando el principio de residuales ponderados se deriva la ecuación matricial para los
elementos, que para elementos con tres nodos queda representada como,
[ ] [ ] (2.12)
Siendo μ un vector de orden N (número de nodos) que contiene todos los valores nodales
del dominio. La matriz K es de orden N x Ny se caracteriza por ser dispersa, simétrica y positiva en
la mayoría de los casos. La matriz B contiene los términos fuente incluidos en las ecuaciones de
Poisson y las condiciones de frontera conocidas.
Dentro de las opciones disponibles para el análisis del campo eléctrico se encuentran:
[Coulomb, 2012; ANSYS/Emag, 2012; Maxwell, 2012; COMSOL Multiphysics, 2012; MATLAB, 2012]
mediante la caja de herramientas de ecuaciones diferenciales parciales, [QuickField, 2012; Opera,
2012; Flux3D, 2012]. Coulomb fue programado utilizando el BEM y ha sido empleado para el
análisis de aisladores compuestos considerando modelos tridimensionales [Que, 2002]; por su
parte [ANSYS/Emag, 2012] está basado en el FEM y su desarrollo ha sido principalmente
encaminado a problemas de ingeniería estructural y térmica y en menor medida a la ingeniería
electromagnética. [Maxwell, 2012] también se basa en el FEM y su orientación es hacia la
21
ingeniería electromagnética.[Flux3D, 2012]tiene un enfoque total sobre la parte de estudios
electromagnéticos tridimensionales. Finalmente [COMSOL Multiphysics, 2012] utiliza el FEM para
la solución de problemas en varias áreas del conocimiento siendo la electromagnética
tridimensional una de ellas.
22
3
ÁNGULO DE BALANCEO EN CADENAS DE AISLADORES
DE LÍNEAS DE TRANSMISIÓN
3.1. Introducción
Uno de los fenómenos que se presenta en cadenas de aisladores, y en todas las partes que
constituyen una línea de transmisión, es el viento. Su presencia genera cargas mecánicas que
afectan la operación mecánica y eléctrica de las cadenas al modificar suposición con respecto a la
estructura que las sostienen. El cambio de la posición angular de la cadena de aisladores a causa
del viento con respecto a su posición original (sin viento) es conocido como “ángulo de balanceo”.
El ángulo de balanceo modifica las distancias dieléctricas del cable conductor de energía, que se
encuentra enganchado a la cadena de aisladores, respecto a la estructura, afectando los campos
eléctricos en la línea de transmisión y sus elementos. Estos cambios ocasionan esfuerzos eléctricos
que afectan el desempeño de las cadenas de aisladores y pueden ocasionar envejecimiento
prematuro, flameo, etcétera [EPRI, 2005].
23
desplazamiento de cables y aisladores, obteniendo así las distancias dieléctricas mínimas con
respecto al cuerpo recto de las estructuras y los obstáculos como arboles ó construcciones a lo
largo de la trayectoria.
Los documentos considerados en este trabajo para la evaluación del ángulo de balanceo
en cadenas de aisladores son: la norma oficial mexicana NOM-001-SEDE-2005, [2006],
“Instalaciones Eléctricas (utilización)”; la especificación CPTT-DDLT-001/02 de la CFE [2009a] y los
reportes técnicos “Tower top geometry and midspan clearances” [CIGRE 348, 2008] y “Guidelines
for the evaluation of the dielectric strength of external insulation” [CIGRE 72, 1992].
3.2.1. NOM-001-SEDE-2005
24
protección de personas contra los efectos eléctricos y térmicos. Es importante mencionar que la
NOM-001-SEDE-2005 no es una guía de diseño.
Otras metodologías [CIGRE 348, 2008; CIGRE 72, 1992; CFE, 2009a] definen de forma más
precisa el ángulo de balanceo por lo que son también incluidas en el presente trabajo para
obtener un mejor entendimiento del fenómeno y su determinación.
25
posición de los elementos y definir así los claros, además de ser preponderante en el cálculo
mecánico y estructural de las líneas de transmisión.
En el documento CIGRE 348 [2008] se establece que bajo las condiciones de un viento
estable ó moderado causante de un ángulo de balanceo pequeño y frecuente, los libramientos
deben ser suficientes para soportar las sobretensiones causadas por descargas atmosféricas y
maniobras con una probabilidad mayor al 90% en la línea de transmisión. Por otro lado, los vientos
extremos ó máximos son eventos raros de tal forma que su coincidencia con sobretensiones por
descargas atmosféricas ó por maniobra son poco probables. Consecuentemente, los diseños
deben cumplir que bajo condiciones de viento extremo los libramientos establecidos soporten la
tensión nominal del sistema.
Dentro del documento CIGRE 348 [2008] se encuentra información recabada de prácticas
internacionales en diversos países para la determinación de ángulo de balanceo de cadenas de
aisladores. En este documento, se observan que las condiciones atmosféricas preferentemente
evaluadas por diversos países para la determinación del ángulo de balanceo y los libramientos
asociados, con respecto a lo definido en la norma europea EN 50341-1[2009], son:
Viento estable ó moderado con máxima temperatura del conductor para esfuerzos de
descargas atmosféricas y maniobra (76% de aplicación).
Viento extremo con periodo de retorno de 50 años para esfuerzos de tensión nominal del
sistema (68% de aplicación).
26
Para la determinación del ángulo de balanceo en cadenas de aisladores, los principales
parámetros a considerar son la velocidad, dirección y periodo de tiempo del viento. Acorde a
resultados de varias investigaciones (Figura 3.1) se ha observado que los valores determinados
teóricamente para el ángulo de balanceo son mayores a los presentados en la realidad CIGRE 348
[2008].
[°]
50
Ángulo de balanceo [grados]
Valores
calculados
25
Valores
medidos
0
0 5 10 15 20 25 30 35 [m/s]
Velocidad del viento [m/s]
Figura 3.1.Valores calculados y medidos del ángulo de balanceo respecto a la
velocidad del viento [CIGRE 348, 2008].
Dentro de las expresiones consideradas como parte del documento CIGRE 348 [2008],
para determinar el ángulo de balanceo ( ) de cadenas de aisladores (Figura 3.1), se tiene la
expresión (3.1) que considera las cargas mecánicas de viento y peso en la cadena de aislamiento
(Figura 3.2):
[ ( ) ( )]
[ ( )]
donde:
27
C XC = coeficiente de arrastre de los conductores (1,0).
d = diámetro del conductor [m].
LW = claro de viento [m].
FWins = carga de viento en la cadena de aisladores [N].
WC = peso del conductor [N].
Wins =peso del conjunto de aislamiento [N].
Cadena de aisladores
Figura 3.2. Cargas mecánicas por viento y peso en cadena de aisladores para determinación de ángulo de
balanceo.
La carga por viento en los conductores, actuando de forma perpendicular a los mismos,
queda determinada mediante,
donde:
28
qZ VZ2 = presión del viento correspondiente a una altura determinada (Z)
2
sobre el terreno [N/m²].
Siendo,
Z
VZ VR =velocidad de viento promedio [m/s].
ZR
ZR= la altura de referencia (10 metros sobre el nivel del terreno ya que es la altura
de ubicación de los medidores de velocidad).
VR= velocidad de viento de referencia [m/s].
= es un factor acorde al tipo de terreno.
GC = factor de ráfaga del conductor dependiente de la categoría de terreno y
altura sobre el terreno.
GL = factor de claro.
( )
[ ]
( )
donde:
29
LC = claro de peso del conductor [m].
Figura 3.3.Registro de ángulos de balanceo [°] respecto a velocidad del viento [m/s] para 4
conductores ACSR 240/40 sin hielo [CIGRE 348, 2008].
Respecto a las mediciones del ángulo de balanceo realizadas en Canadá [Houle, 2000] se
tienen relaciones entre éstos y los valores calculados con la expresión (3.3), considerando
velocidades de viento de 5 minutos promedio. En la Figura 3.4 se considera un haz de cuatro
conductores y un claro de 320 metros, donde la correspondencia entre los valores calculados y los
30
medidos es muy similar. Se observa una tendencia en la que los valores calculados son mayores a
los medidos.
4 conductores en haz
Ángulo medido (grados)
(ACSR Bersfort)
Claro de 320 m
Ajuste lineal
Figura 3.4.Ángulos de balanceo medidos contra calculados en Canadá [CIGRE 348, 2008].
Con base en los resultados experimentales obtenidos, el Electric Power Research Institute
[EPRI, 1982] propuso la determinación del ángulo de balanceo acorde a un diagrama que
considera la relación entre la velocidad del viento y un parámetro kW definido mediante (Ver
Figura 3.5),
donde,
31
[°]
60
k w=35
32
40 25
20
30
15
20
10
7.5
10
Figura 3.5.Determinación del ángulo de balanceo en aisladores para valores de kw basados en mediciones
[EPRI, 1987]. Para una velocidad del viento de 120 km/h y un parámetro de kw=20 el ángulo de balanceo del
aislador es 20° [CIGRE 348, 2008].
De forma global se tiene que las condiciones de análisis para la determinación del ángulo
de balanceo, y consecuentemente la geometría superior de la estructura, acorde al documento
CIGRE 348[2008], se encuentran definidas en la Tabla 3.1.
Tabla 3.1. Coordinación de posición de conductores y cadena de aisladores respecto a esfuerzos eléctricos
[CIGRE 348, 2008].
32
3.2.2.1. Libramientos
En relación a los libramientos necesarios de fase a tierra se establece en los estándares IEC
60071-1[1993]e IEC 60071-2 [1996] que los voltajes mayores de los equipos se encuentran
divididos en dos rangos. El rango 1 abarca desde 1kV hasta los 245kV, mientras que el rango 2 es
para voltajes mayores de 245kV. Dentro del rango 1, las descargas atmosféricas definen los niveles
de aislamiento mientras que para el rango 2 son definidos mediante las sobretensiones por
maniobra. La Tabla 3.2 muestra algunos niveles de tensión y sus libramientos acorde a los frentes
de onda definidos para una probabilidad de 90% de soporte con coeficientes de variación de 3%
para descargas atmosféricas y 6% para maniobras.
donde,
Las condiciones consideradas en la Tabla 3.2 son: temperatura de 20°C, presión del aire de
1013 mbar y humedad de 11 g/m³ para altitudes de hasta 1000 msnm.
33
considerada por IEC 60071-2[1996], una opción para determinar el libramiento es lo establecido
en el trabajo de CIGRE 72 [1992], donde se establecen las siguientes expresiones:
[ ]
donde,
√
[ ]
donde,
Tabla 3.3. Libramientos para el nivel de voltaje a la frecuencia del sistema [CIGRE 348, 2008].
34
4.3. Comisión Federal de Electricidad
Las hipótesis de carga H3 “viento máximo” y H6 “Balanceo de cables” son parecidas a las
descritas previamente en el documento CIGRE 348[2008], por lo cual son viables para la
determinación de las condiciones climatológicas necesarias para el cálculo del ángulo de balanceo
en cadenas de aisladores. La presión del viento definida para la hipótesis de carga H6 es de 285 Pa
para zonas expuestas y rurales ó de 187 Pa para zonas urbanas y protegidas. Al observar estos
valores se presenta una similitud respecto a los valores establecidos en la norma NOM-001-SEDE-
2005[2006], pero que considera unidades de kg/m² en lugar de Pa. En el caso de la hipótesis de
carga H3, la determinación de la presión del viento en cables que afecta el ángulo de balanceo en
los conductores, y que consecuentemente sirve como parte de la información para su
determinación, se considera lo establecido en las especificaciones “Torres para líneas de
subtransmisión y transmisión” [CFE, 2006] y “Postes metálicos para líneas de transmisión y
subtransmisión” [CFE, 2008], cuyo Apéndice A define las metodologías para el cálculo de presiones
de viento a partir de velocidades regionales del viento en México.
35
4
CASOS DE ESTUDIO EN EL ANÁLISIS DEL CAMPO ELÉCTRICO EN
CADENAS DE AISLADORES DE LÍNEAS DE TRANSMISIÓN
4.1. Introducción
Las estructuras comúnmente empleadas dentro de las líneas de transmisión son torres
autosoportadas y postes troncocónicos. En este trabajo de tesis se considera el uso de modelos
36
representativos de ambos tipos de estructuras. Las torres autosoportadas son comúnmente
empleadas en proyectos rurales donde los derechos de vía no son tan restrictivos, mientras que
los postes troncocónicos se aplican principalmente en zonas urbanas donde se dispone de poco
espacio para la ubicación de los conductores y presentan un menor impacto estético para el
ambiente.
Las zonas rurales ó con terrenos abiertos presentan velocidades de viento mayores, que
causan desplazamientos grandes de conductores y cadenas de aisladores, por lo que se requieren
dimensiones mayores en las estructuras para garantizar que las distancias dieléctricas sean
apropiadas. Por otro lado, en las zonas urbanas debido a la cantidad de obstáculos que se
encuentran en el terreno se presenta una velocidad de viento menor y mayor turbulencia,
resultando en dimensiones menores de los postes troncocónicos.
Para todos los casos de estudio del presente trabajo se consideran estructuras tipo
suspensión ya que son susceptibles de exhibir un mayor cambio en su ángulo de balanceo debido
a las presiones transversales del viento. De igual forma, los modelos establecidos se asocian a la
cruceta o brazo medio de las estructuras con la finalidad de observar y analizar la interacción del
campo eléctrico y el potencial entre las estructuras y las cadenas de aisladores.
Los ángulos de balanceo considerados para todos los casos de estudio son 0°, 15° y 30°con
la finalidad de observar los cambios en los esfuerzos eléctricos en las cadenas de aisladores con
respecto a la modificación en su distancia al cuerpo de las estructuras ocasionado por la acción del
viento. El valor máximo de ángulo de balanceo considerado es de 30° ya que está establecido en la
NOM-001-SEDE-2005.Un valor intermedio de 15° se evalúa para observar los cambios intermedios
entre las condiciones máximas (deflexión máxima) y mínimas (a plomo).
37
Las cadenas de aisladores de 115 y 230 kV presentan los patrones de campo eléctrico y
potencial mostrados en las Figuras 4.3 a 4.36 para los casos de disposición vertical y bajo un
ángulo de balanceo. Los valores obtenidos son considerando el comportamiento interno a lo largo
de las cadenas. En los casos de disposición vertical se obtienen los patrones hasta alcanzar el brazo
inmediato inferior de la estructura de transmisión. En el caso de la estructura de 400 kV se
considera la condición de cadenas de aisladores en disposición en “V” para este trabajo.
4.2.1 Características
Las condiciones y características eléctricas definidas para los elementos que conforman los
modelos geométricos elaborados para los casos de estudio son descritas a continuación.
Las propiedades dieléctricas de cada una de las partes que integran un aislador se
observan en la Tabla 4.1. En la Figura 4.1 se observa un corte de sección transversal del aislador
estándar 25SVC111 de vidrio templado considerado en el presente trabajo. El perfil del modelo del
aislador es definido con base en información proporcionada por el fabricante NGK [2007]. Al
encontrarse enganchadas directamente las partes metálicas del aislador en una cadena de
aisladores (calavera y perno),éstas conforman en el modelo desarrollado un solo elemento ya que
sus propiedades dieléctricas son parecidas.
Las propiedades dieléctricas asignadas a la fuente de excitación (cable conductor) son las
del material aluminio, las secciones de las estructuras de transmisión tienen asignadas las
propiedades dieléctricas del acero. Finalmente, los materiales de la cadena de aisladores que
presentan menor conductividad eléctrica son el vidrio templado y el cemento (ver Tabla 4.1).
38
Tabla4.1. Propiedades dieléctricas de elementos utilizados dentro de los casos de estudio.
Permitividad relativa (r) Conductividad eléctrica () Voltaje de fuente (V0)
Material
[ ] [S/m] [kV]
Aire 1 0,0
6
Acero 1 4,032e 0,0
7
Hierro 1 1,12e
66,395 (para 115 kV)
7
Aluminio 1 3,774e 132,791 (para 230 kV)
230,940 (para 400 kV)
-14
Vidrio 2,09 1e
-12
Cemento 12,1 1e
Aire
Calavera y perno
(hierro).
Esbozo (vidrio).
Aislante (cemento).
Calavera y perno
(hierro).
Figura 4.1. Sección transversal de aislador definida dentro de los casos de estudio.
Las estructuras normalizadas de la CFE consideradas para especificar las dimensiones del
modelo son (Figura 4.2):
39
Poste troncocónico 1210SMP, tipo suspensión de 2 circuitos.
Torre autosoportada E71A21, tipo suspensión de 2 circuitos.
El calibre del cable conductor en este nivel de tensión es 795 kCM (28,11 mm de diámetro
exterior) con un conductor por fase.
5,1
2,8
4,85
3,2
1,4 1,4
2,3
3,18
28,4
26,33
1210SMP E71A21
40
El nivel de voltaje ó potencial eléctrico escalar de la fuente de excitación (cable
conductor), respecto a tierra, establecido en el modelo es de 66,395 kV. Un nivel de potencial
0,0kV se considera para toda la sección de las estructuras establecidas.
Los modelos tridimensionales desarrollados para análisis por el método de los elementos
finitos del poste 1210SMP y la torre E71A21 se observan en las Figuras4.3 a 4.18. Como parte
delos modelos desarrollados se establece una sección de conductor 795 kCM de 2 metros de
longitud. El modelo considerado para la cadena de aisladores está constituido de secciones
continuas de unidades de aisladores, cada una de las cuales se encuentra conformada por su pin o
perno, cemento, esbozo de vidrio templado y calavera acorde a lo observado en la Figura 4.1 y
considerando las propiedades de la Tabla 4.1.
Figura 4.3. Distribución de potencial en cadena de aisladores de 115 kV con disposición vertical en poste
1210SMP.
41
[V]
c) vista isométrica
Figura 4.3. Distribución de potencial en cadena de aisladores de 115 kV con disposición vertical en poste
1210SMP.Continuación.
42
Potencial eléctrico
Potencial [V]
eléctrico [V]
Longitud [m]
Figura 4.4. Distribución de potencial a lo largo de cadena de aisladores de 115 kV en disposición vertical para
poste 1210SMP (distancia brazo medio a brazo inferior de la estructura).
c) vista isométrica
Figura 4.5. Líneas de flujo de campo eléctrico en cadena de aisladores de 115 kV con disposición vertical en
poste 1210SMP.Continuación.
Respecto al campo eléctrico en el aire, la Figura 4.6 muestra que la intensidad mayor del
campo eléctrico es de aproximadamente 480,0 kV/m. Se observa en la Figura 4.6 que la parte
izquierda de la gráfica no presenta las variaciones del campo eléctrico que se encuentran en la
parte derecha, debido a que es un medio con propiedades dieléctricas homogéneas (aire). Los
cambios de las propiedades dieléctricas de los elementos constitutivos de los aisladores causan la
variación en el campo eléctrico que se observa en la parte derecha de la Figura 4.6, al pasar de un
elemento a otro. Los valores pico en la Figura 4.6 corresponden al interior de secciones metálicas
(hierro) de la cadena mientras que los valores mínimos y en decremento corresponden a las partes
dieléctricas de menor conductividad eléctrica (vidrio templado y cemento).
Mediante las Figuras 4.4 y 4.6 es posible observar los comportamientos de la distribución
de potencial y campo eléctrico de forma clara para una trayectoria lineal interna determinada. Sin
embargo, la Figura 4.7 ayuda a observar de cerca detalles relevantes en la trayectoria de los
fenómenos descritos entre las diferentes fronteras de los elementos.
44
Campo eléctrico [V/m]
Longitud [m]
Figura 4.6. Campo eléctrico a lo largo de cadena de aisladores de 115 kV en disposición vertical
para poste 1210SMP (distancia brazo medio a brazo inferior de la estructura).
[V/m]
Acorde a la Figura 4.8, los valores más grandes de potencial se encuentran en los primeros
cuatro aisladores de la cadena y presentan un cambio más pequeño en el resto de los aisladores
45
de la cadena, observándose un comportamiento similar del potencial al del modelo realizado para
la cadena de aisladores de la estructura tipo poste 1210SMP.
[V] [V]
c) vista isométrica
Figura 4.8. Distribución de potencial en cadena de aisladores de 115 kV con disposición vertical en torre
E71A21.
46
observan las protuberancias de la curva de potencial eléctrico causadas por la presencia de los
aisladores y las diferentes propiedades dieléctricas de los elementos que los conforman.
Longitud [m]
Figura 4.9. Distribución de potencial a lo largo de cadena de aisladores de 115 kV en disposición vertical para
torre E71A21 (distancia cruceta media a cruceta inferior de la estructura).
La vista frontal de la Figura 4.10 muestra que las líneas de flujo del campo eléctrico tienen
una interacción mayor respecto al cuerpo recto para la torre E71A21, comparado con la misma
vista de la estructura tipo poste 1210SMP (Figura 4.5). Esto es causado por las diferencias entre la
forma y las dimensiones de las crucetas y brazos de la torre E71A21 y el poste 1210SMP,
respectivamente. Consecuentemente, la distribución del campo eléctrico en el aire es más
homogénea, mientras se conserve constante la distancia de la fuente de excitación con respecto a
las fronteras que lo rodean en todas las direcciones.
47
[V/m] [V/m]
c) vista isométrica
Figura 4.10. Líneas de flujo de campo eléctrico en cadena de aisladores de 115 kV con disposición vertical
para torre E71A21.
48
de aisladores afectando el valor del campo eléctrico. El campo eléctrico decrece a lo largo de la
cadena, pero en el último aislador se observa un incremento antes de alcanzar la estructura, en la
cual el valor de campo eléctrico es despreciable.
Longitud [m]
Figura 4.11. Campo eléctrico a lo largo de cadena de aisladores de 115 kV en disposición vertical para torre
E71A21 (distancia cruceta media a cruceta inferior de la estructura).
49
La mayor parte del esfuerzo eléctrico se encuentra presente en los primeros cuatro elementos de
la cadena.
[V] [V]
50
En la Figura 4.14 se observan las líneas del campo eléctrico para la cadena de aisladores
con un ángulo de balanceo de 30°. En comparación con el caso de disposición vertical, la
interacción de las líneas de flujo del campo eléctrico, entorno a la cadena de aisladores, con la
estructura 1210SMP es mayor cuando se presenta deflexión de la cadena por viento como se
puede observar en la vista frontal de las Figuras 4.5 y 4.14.
[V/m] [V/m]
La Figura 4.15 muestra el campo eléctrico a lo largo de la cadena de aisladores para los
casos de 30 y 15° de deflexión. A diferencia del caso en el que se presenta una disposición vertical
en la cadena de aisladores (Figura 4.6), se tiene que la interacción del campo eléctrico con la
estructura es mayor, debido a que se acortan las distancias entre la fuente de excitación (cable
conductor) y la zona de potencial cero (brazo y cuerpo recto de la estructura). Consecuentemente,
los esfuerzos eléctricos máximos en la cadena se ven incrementados pasando de un valor de 700,0
kV/m en disposición vertical a valores de 970,0 y 960,0 kV/m para deflexiones de 15 y 30°
respectivamente.
La Figura 4.15 no muestra para los casos de ángulo de balanceo de 30 y 15° la parte del
comportamiento del campo eléctrico en el aire a diferencia de las Figuras 4.6 y 4.11 previas debido
a que el trabajo se encuentra enfocado en el comportamiento interno de la cadena de aisladores.
La conformación de superficies ó curvas de nivel para la distribución del potencial y el campo
eléctrico representan otras formas útiles de visualización del comportamiento de variables
eléctricas.
51
Campo eléctrico [V/m]
Campo eléctrico [V/m]
52
Potencial eléctrico [V]
Potencial eléctrico [V]
La Figura 4.19 muestra las curvas de comportamiento del campo eléctrico a lo largo de la
cadena de aisladores para la estructura E71A21 bajo un ángulo de balanceo. En el caso de campo
eléctrico se observa un valor del segundo aislador de la cadena similar en ambos casos 15 y 30°
53
(270,0 kV/m). El primer aislador de la cadena presenta un valor máximo próximo de 810,0 kV/m
para 15° y de 970,0 kV/m para 30° de ángulo de deflexión ó balanceo.
Campo eléctrico [V/m]
Las Figuras 4.12 a 4.15 muestran el comportamiento del campo eléctrico y la distribución
de potencial a lo largo de la cadena de aisladores para la estructura tipo poste troncocónico
1210SMP para los ángulos de balanceo de 15° y 30°, donde se observa que las diferencias del
potencial son pequeñas mientras que el campo eléctrico presenta un cambio de magnitud visible
siendo menor para el caso de 30° (960,0 kV/m) con respecto al de 15° (970,0 kV/m) de ángulo de
balanceo. Las Figuras 4.16 a 4.19 muestran el comportamiento del campo eléctrico y distribución
de potencial a lo largo de la cadena de aisladores para la estructura tipo torre autosoportada
E71A21, donde se observa un incremento del campo eléctrico conforme aumenta el ángulo de
balanceo.
54
Torre autosoportada 2B2, tipo suspensión de 2 circuitos.
El calibre del cable conductor en este nivel de tensión es 900 kCM (29,51 mm de diámetro
exterior) con un conductor por fase.
7,28
5,9
7,5
5,5
2,7 2,7
3,8
4,2
38,96
36,7
AS2 2B2
55
Todos los elementos previos se encuentran inmersos dentro de fronteras exteriores de
dimensiones cúbicas de 8 por 8 metros para el caso del poste AS2 y de 8 por 8 por 11,52 metros de
alto para la torre 2B2.
Los modelos tridimensionales desarrollados para análisis por el método de los elementos
finitos para el poste AS2 y la torre 2B2 se observan en las Figuras 4.21a 4.36. Como parte de los
modelos se establece una sección de conductor 900 kCM de 2 metros de longitud.
La Figura 4.21 muestra la distribución del potencial en la cadena de aisladores desde tres
vistas diferentes. Los dos primeros aisladores de la cadena se encuentra por encima de los 40,0 kV,
que representa el 30% del valor máximo (132,791 kV) de potencial. Los restantes dieciséis
aisladores de la cadena presentan una distribución de potencial más homogénea y por debajo de
los 50,0 kV. Como resultado, solamente el 11% de la longitud de la cadena se encuentra esforzada
con valores superiores al 30% del potencial máximo y el restante 89% de la longitud de la cadena
se encuentra esforzada por debajo del 30%. El comportamiento de la distribución del potencial se
observa en forma de escala de color acorde a la barra lateral donde el valor máximo es 1,328e5 V
(132,791 kV) y el mínimo 0,0 kV.
[V] [V]
56
[V]
c) vista isométrica
Figura 4.21. Distribución de potencial en cadena de aisladores de 230 kV con disposición vertical en poste
AS2.Continuación.
57
Potencial eléctrico [V]
Longitud [m]
Figura 4.22. Distribución de potencial a lo largo de cadena de aisladores de 230 kV en disposición vertical
para poste AS2 (distancia brazo medio a brazo inferior de la estructura).
Los valores máximos y mínimos del campo eléctrico presentados en la escala de la Figura
2.23 son respectivamente 2050,0 kV/m y 4,348 V/m respectivamente. Con relación al valor
máximo obtenido en la Figura 2.24 (520,0 kV/m) es mayor aproximadamente cuatro veces. Las
diferencias son importantes de tal forma que es significativo considerar los valores máximos del
modelo numérico probando los diferentes tipos de visualización de los resultados. Aun así, las
gráficas de comportamiento como las Figuras 2.23 y 2.24 son relevantes en el análisis de los
fenómenos de esfuerzos eléctricos.
58
[V/m] [V/m]
[V/m]
c) vista isométrica
Figura 4.23. Líneas de flujo de campo eléctrico en cadena de aisladores de 230 kV con disposición vertical en
poste AS2.
La Figura 4.24 muestra el comportamiento del campo eléctrico. En la curva se observa que
en la parte izquierda, aproximadamente a 1 metro de distancia de la fuente de excitación (punto
1,7 del eje vertical), el nivel de campo eléctrico disminuye y se mantiene hasta alcanzar el último
aislador de la cadena, donde aumenta la intensidad un poco, y es mínimo en el brazo medio del
poste AS2 (de 0 a 0,25 m aproximadamente del eje horizontal). En la parte derecha de la curva, las
propiedades dieléctricas continuas del aire influyen en el comportamiento del campo eléctrico a lo
largo de la cadena.
59
Campo eléctrico [V/m]
Longitud [m]
Figura 4.24. Campo eléctrico a lo largo de cadena de aisladores de 230 kV en disposición vertical para poste
AS2 (distancia brazo medio a brazo inferior de la estructura).
c) vista isométrica
Figura 4.25. Distribución de potencial en cadena de aisladores de 230 kV con disposición vertical en torre
2B2.Continuación.
La Figura 4.27 presenta tres vistas (frontal, lateral e isométrica) del campo eléctrico a lo
largo de la cadena de aisladores y en el entorno mediante sus líneas de flujo de campo. La
representación geométrica utilizada para la sección de la torre 2B2 muestra una sección del
cuerpo recto, cruceta media y cruceta inferior. La zona interna y alrededor de la cadena de
aisladores es el área en la cual las líneas del campo eléctrico se encuentran más concentradas,
indicando intensidades mayores. En comparación con los casos de 115 kV, se observa que las
líneas de flujo del campo para 230 kV no se ven influenciadas tanto con las fronteras externas del
modelo numérico tridimensional que se incrementaron en 2 metros pasando de 6 a 8 metros
respectivamente.
61
Potencial eléctrico [V]
Longitud [m]
Figura 4.26. Distribución de potencial a lo largo de cadena de aisladores de 230 kV en disposición vertical
para torre 2B2 (distancia cruceta media a cruceta inferior de la estructura).
El número de líneas de flujo del campo eléctrico definidas para representar la visualización
en la Figura 4.27 dentro del programa de cómputo son cincuenta. Este tipo de representación
facilita observar la trayectoria del campo eléctrico desde la fuente de excitación (cable conductor)
y su interacción con el entorno que son: aire, cadena de aisladores, estructura y fronteras
externas.
[V/m] [V/m]
c) vista isométrica
Figura 4.27. Líneas de flujo de campo eléctrico en cadena de aisladores de 230 kV con disposición vertical
para torre 2B2.Continuación.
El campo eléctrico a lo largo de una trayectoria recta de 7,6 metros definida entre los
extremos finales de las crucetas media e inferior de la torre 2B2 se presenta en la Figura 4.28. Se
observa en la curva los cambios entre las propiedades dieléctricas de los medios, primero el
comportamiento a lo largo de la cadena de aisladores, después la fuente de excitación (cable
conductor y herrajes), y finalmente en el aire.
De las estructuras tipo AS2 y 2B2 para 230 kV se presenta una mayor intensidad del campo
eléctrico a lo largo del centro de la cadena de aisladores para la torre 2B2 (Figura 4.28). Los valores
máximos son 520,0 y 1270,0 kV/m respectivamente en ambas estructuras. Sin embargo,
considerando los valores máximos mostrados en las Figuras 4.23 y 4.27 los valores máximos son
2050,0 y 1901,0 kV/m respectivamente. Debido a que las curvas 4.24 y 4.28 determinan el campo
eléctrico en una trayectoria lineal para la parte central de la cadena de aisladores, el valor máximo
del campo eléctrico es menor. Esto es porque la conformación geométrica y de las propiedades de
los materiales dieléctricos influye estableciendo un valor mayor en las secciones con vértices más
pronunciados. Consecuentemente, el perímetro del cemento, esbozo y perno presentan valores
mayores del campo eléctrico. Con el fin de mejorar el diseño de un elemento aislante, para un
mejor comportamiento respecto a la distribución de los esfuerzos eléctricos, las geometrías curvas
y homogéneas son recomendables.
63
Campo eléctrico [V/m]
Longitud [m]
Figura 4.28. Campo eléctrico a lo largo de cadena de aisladores de 230 kV en disposición vertical para torre
2B2 (distancia cruceta media a cruceta inferior de la estructura).
Los modelos numéricos tridimensionales desarrollados para análisis por el método de los
elementos finitos para el poste troncocónico AS2 y la torre autosoportada 2B2 consideran ahora
un ángulo de balanceo de 15 y 30°de la cadena de aisladores como se observan en las Figuras 4.29
a 4.36. Valores de ángulo de balanceo en sentido externo a la estructura no se consideran ya que
los esfuerzos eléctricos máximos en la cadena se presentan al acercarse la misma a la estructura
[Mehdi, 2010].
64
[V] [V]
65
[V/m] [V/m]
66
[V] [V]
eléctrico [V]
Potencial eléctrico
Potencial [V]
La Figura 4.35 muestra las líneas del campo eléctrico presentes en el entorno de los
aisladores con un ángulo de balanceo de 30°. En la parte final de la cadena de aisladores se
observa que el campo eléctrico interactúa más con la cruceta media de la torre 2B2 que para el
caso de deflexión vertical (Figura 4.27). Esto es consecuencia del cambio de posición de la cadena.
67
[V/m] [V/m]
Para la estructura AS2 se observa que el potencial (Figuras 4.29 y 4.30) presenta cambios
pequeños en las curvas de comportamiento considerando ángulos de balanceo de 15 y 30°.El
campo eléctrico (Figuras 4.31 y 4.32) si muestra diferencias visibles en las gráficas siendo mayor el
campo eléctrico para el caso de 30° de ángulo de balanceo en la cadena de aisladores.
68
Para la estructura 2B2 se muestran cambios en el potencial (Figuras 4.33 y 4.34), y
consecuentemente la eficiencia de la cadena de aisladores. El campo eléctrico presenta diferencias
visibles, siendo mayor el campo para el caso de 15° de ángulo de balanceo (Figura 4.36) en la
cadena de aisladores.
La estructura de la CFE considerada para especificar las dimensiones del modelo es (Figura
4.37):
El calibre del cable conductor en este nivel de tensión es 1113 kCM (31,96 mm de
diámetro exterior) con dos conductores por fase.
Como parte del análisis de cadenas de aisladores de 400 kV se establecen modelos de las
cadenas de aisladores y herrajes simplificados debido a limitantes del equipo de cómputo
empleado. Se consideran 2 aisladores en la parte terminal del conductor y un aislador en la
sección de enganche con la estructura para las cadenas en configuración “V”. Todos los elementos
69
previos se encuentran inmersos dentro de una frontera externa de dimensiones cúbicas de 15 por
15 metros.
10,4
10,3
4,5
45
9,43
EA4A22MA
La Figura 4.38 muestras tres vistas (frontal, lateral e isométrica) del modelo geométrico
tridimensional establecido de la torre EA4A22MA y la distribución del potencial a lo largo de las
cadenas de aisladores. El esfuerzo eléctrico no lineal es visible ya que prácticamente la mitad de la
70
longitud de la cadena presenta un potencial debajo de 50,0 [kV], que representa solamente el 21%
del valor máximo del potencial.
[V] [V]
c) vista isométrica
Figura 4.38. Distribución de potencial en cadena de aisladores de 400 kV en configuración “V” para torre
EA4A22MA.
71
Potencial eléctrico [V]
Longitud [m]
Figura 4.39. Distribución de potencial en cadena de aisladores de 400 kV en configuración “V” para torre
EA4A22MA.
En la Figura 4.40 se observa el comportamiento del campo eléctrico cerca de las cadenas
de aisladores mediante las líneas de campo. El campo eléctrico se distribuye entre las dos
secciones de la cadena de aislamiento. Se presenta una mayor interacción de las líneas de flujo de
campo eléctrico de la cadena derecha del conjunto de aislamiento con la estructura (ver Figura
4.40, vista frontal), esto es debido a un mayor grado de detalle definido en la cadena derecha.
[V/m] [V/m]
c) vista isométrica
Figura 4.40. Líneas de flujo de campo eléctrico en cadena de aisladores de 400 kV en configuración “V” para
torre EA4A22MA.Continuación.
La Figura 4.41 muestra la gráfica del comportamiento del campo eléctrico a lo largo de una
de las cadenas de aisladores. El valor máximo presente del campo es de aproximadamente 3800,0
kV/m que es el mayor de todos los casos de estudio para 115, 230 y 400 kV. Comparando este
valor con respecto al necesario para que el aire se vuelva un conductor que es de
aproximadamente 3000,0 kV/m [Hewitt, 1995] observamos un incremento del 26%. Como
consecuencia, los efectos de los fenómenos que ocasionan pérdidas y daños en los elementos en
la línea de transmisión (corona, térmicos, envejecimiento prematuro, descargas parciales, etc.) son
mayores.
73
Campo eléctrico [V/m]
Longitud [m]
Figura 4.41. Campo eléctrico encadena de aisladores de 400 kV en configuración “V” para torre EA4A22MA.
4.6. Modelo geométrico de cadena de aisladores con cambios en longitud del cable
conductor
74
Cambio en
de potencial
a lo largo de
la cadena de
aisladores
a) 2 metros b) 4 metros
Figura 4.43. Comparativa del campo eléctrico en cadena de aisladores de 115 kV para poste 1210SMP con
longitudes de conductor diferentes.
75
A partir de los resultados mostrados en las Figuras 4.42 y 4.43 se observa que la
evaluación de modelos axisimétricos de cadenas de aisladores, donde la fuente de excitación es
considerada como una carga puntual, no representa de manera correcta el comportamiento de la
distribución e intensidad del campo eléctrico y el potencial. Esto concuerda con los resultados
obtenidos por [Zhao, 2000] que establece: a) para la determinación del campo eléctrico y la
distribución de potencial en los aislamientos, los modelos desarrollados a través de técnicas
numéricas son una opción viable; b) para evaluar el campo eléctrico y la distribución de potencial
en los conjuntos de aislamiento instalados y en operación, se requiere un modelo geométrico
tridimensional que considere a la estructura y los conductores para evitar errores en la estimación
de los valores; y c) se pueden realizar simplificaciones en el modelo geométrico tridimensional de
los conjuntos de aislamiento sin perder precisión, siempre y cuando, elementos como la
estructura y los conductores se encuentren representados dentro del mismo.
76
[V]
Figura 4.44. Modelo geométrico establecido de sección de poste 1210SMP de 115 kV con fronteras externas
de 20 metros (distribución de potencial).
La Figura 4.45 muestra el comportamiento del campo eléctrico a lo largo del centro de la
cadena de aisladores de 115 kV y en una sección del aire hasta alcanzar el brazo inferior del poste
1210SMP. Se observa dentro de los resultados que, conforme disminuyen las dimensiones del
cubo que define las fronteras externas del modelo geométrico, el valor del campo eléctrico
aumenta y se presentan variaciones en la distribución del potencial. Esto es consecuencia de que
el campo eléctrico es inversamente proporcional a la distancia con respecto a la carga eléctrica y
las capacitancias en el modelo se ven modificadas.
En los resultados de la Figura 4.45se observa que aproximadamente en el punto 1,5 del eje
horizontal, correspondiente a la unión entre la fuente de excitación y el aislamiento, el campo
eléctrico presenta un valor entre 300,0 y 320,0 kV/m aproximadamente una variación del 6%.
Campo eléctrico [V/m]
c) frontera de 12,8 metros (300,0 [kV/m]) d) frontera de 10,24 metros (305,0 [kV/m])
Campo eléctrico [V/m]
Longitud [m]
Los valores presentes del campo eléctrico en el aire cambian en la Figura 4.45sobre un
rango entre 300,0 y 430,0 kV/m aproximadamente, lo cual representa un cambio del 43%.Esto es
resultado de los cambios en las dimensiones de las fronteras externas y la conformación de la
malla de elementos finitos.
78
4.8. Resultados y análisis
La Tabla 4.2 sintetiza los resultados obtenidos para los casos de estudio que se manejaron
en relación al ángulo de balanceo en cadenas de aisladores. Es evidente que el cambio de posición
de la fuente de excitación (cable conductor) con respecto a las estructuras influye en el
comportamiento del campo eléctrico en la cadena de aisladores, ya que se modifican las
capacitancias en los modelos geométricos al variar las distancias entre los elementos.
Tabla 4.2. Campo eléctrico en cadenas de aisladores ante diferentes estructuras y ángulos de balanceo.
Ángulo de balanceo [°]
Nivel de tensión
Tipo de estructura 0 15 30
[kV]
Campo eléctrico [kV/m]
1210SMP 700,0 970,0 960,0
115
E71A21 560,0 810,0 970,0
AS2 520,0 1320,0 1860,0
230
2B2 1270,0 1670,0 1410,0
400 EA4A22MA 3800,0
Para las estructuras de 115 kV, el valor el campo eléctrico para la cadena de aisladores en
el poste troncocónico 1210SMP, sin cambio de posición debido al viento (0°), es de 700,0 kV/m.
Un valor del campo eléctrico de 970,0 kV/m se presenta ante un ángulo de balanceo de 15°.
Consecuentemente, la intensidad del campo eléctrico en la cadena se incrementa un 38,6%. Por
otra parte, considerando la estructura tipo torre autosoportada E71A21 se observa un incremento
de 560,0 kV/m a 970,0 kV/m en el campo eléctrico de la cadena de aisladores al presentarse un
desplazamiento angular de 30° en la cadena, esto es un 73,2%, como resultado de la disminución
de la distancia con respecto a la estructura. Los resultados obtenidos coinciden con [Izgi, 2008] y
[Mehdi, 2010], respecto a un incremento del esfuerzo eléctrico en las cadenas de aisladores,
cuando se presenta desplazamiento angular que acerca la cadena a la estructura. Mientras [Izgi,
2008] y [Mehdi, 2010] evalúan la distribución del potencial, en el presente trabajo se analiza la
intensidad del campo eléctrico en la cadena que es una variable más fácil de asociar en relación a
gradientes disruptivos.
Respecto a las estructuras de 230 kV, se tiene que el campo eléctrico a lo largo de la
cadena de aisladores presenta un incremento del intensidad de 520,0 kV/m en el poste
troncocónico AS2 en disposición vertical a un valor aproximado de 1860,0 kV/m ante un ángulo de
balanceo de 30°, esto es un 357,7%. En la torre autosoportada 2B2 también se incrementa la
intensidad del campo eléctrico en la cadena de aisladores debido a que las distancias con respecto
79
a la estructura se reducen. Consecuentemente, la intensidad del campo eléctrico en la cadena de
aisladores de 1270,0 kV/m en condición de disposición vertical aumenta a 1670,0 kV/m cuando se
encuentra bajo un ángulo de balanceo de 15°, esto es un 31,5%. Los resultados obtenidos coincide
con [Mehdi, 2010] respecto a un aumento en los esfuerzos eléctricos en las cadenas. Las
intensidades de campo eléctrico en estos casos comparadas con el umbral para inicio de corona
(1500,0 kV/m) [Tanushri, 2010] son mayores. Consecuentemente, la implementación de anillos
equipotenciales en cadenas de aisladores para 230 kV se recomienda como una opción para
disminuir los esfuerzos del campo eléctrico.
80
5
CONCLUSIONES
5.1. Conclusiones
El análisis del campo eléctrico y la distribución del potencial eléctrico escalaren cadenas de
aisladores, mediante la aplicación del método de los elementos finitos, permite conocer su
comportamiento ante diversas configuraciones físicas de la cadena en relación a su posición con
respecto a la estructura que la sostiene y bajo diversas condiciones de operación.
Los resultados obtenidos en el presente trabajo de tesis coinciden con [Izgi, 2008] y
[Mehdi, 2010], respecto a un incremento del esfuerzo eléctrico en las cadenas de aisladores,
cuando se presenta desplazamiento angular que acerca la cadena a la estructura. Mientras [Izgi,
2008] y [Mehdi, 2010] evalúan la distribución del potencial, en el presente trabajo se analiza la
intensidad del campo eléctrico en la cadena que es una variable más fácil de asociar en relación a
gradientes disruptivos y aparición de fenómenos que afectan a los aislamientos en general como
las descargas parciales. Adicionalmente, los resultados obtenidos permiten conocer las partes
internas del aislador que presentan esfuerzos eléctricos mayores, siendo en el caso de las cadenas
de aisladores el primer elemento de la cadena en la parte de interacción física entre el perno, el
cemento y el esbozo.
81
del ángulo de balanceo. La NOM-SEDE-001-2005[2006] establece condiciones para la
determinación de la presión de viento en las líneas de transmisión y sus elementos, sin embargo,
el artículo 922-12 señala un ángulo de balanceo en las cadenas de aisladores de 30°. Por su parte
la especificación CPTT-DDLT-001/02 de [CFE, 2009a] define siete hipótesis de carga a considerar en
el diseño electromecánico de la línea de transmisión, siendo las hipótesis H3 “viento máximo” y H5
“Balanceo de cables”, como aplicables para determinar el ángulo de balanceo. Debido a que las
Normas Oficiales Mexicanas tienen mayor injerencia respecto a las especificaciones internas de
empresas. El análisis bajo la condición ángulo de balanceo de 30° en cadenas de aisladores aplica
para la CFE, independientemente de otros valores que se determinen del uso de sus
especificaciones internas y que sean aplicados ante condiciones operativas particulares de un
proyecto específico.
Una restricción importante en el análisis del campo eléctrico en las cadenas de aisladores
mediante el FEM es que, desarrollar modelos geométricos tridimensionales detallados, requiere
de equipos de cómputo con importantes capacidades de procesamiento. De tal forma que, la
simplificación de los modelos geométricos numéricos y la aplicación de otros métodos numéricos
como por ejemplo el BEM [Que, 2002; Zhao, 2000] resultan ser recomendables, debido a su
característica de menor complejidad en la definición del modelo, aunque se incrementa el tiempo
de solución computacional.
82
Un trabajo futuro involucra el desarrollo de modelos completos para las líneas de
transmisión trifásicas de uno y dos circuitos para diversos niveles de tensión, con la finalidad de
evaluar y analizar el comportamiento de los campos eléctricos y potenciales en las cadenas de
aisladores ante la presencia de las otras fases de la línea.
83
BIBLIOGRAFÍA
4. Asenjo, E., Morales, N., and Valdenegro, A. Solution of low frequency complex fields in
polluted insulators by means of the Finite Element Method. IEEE Transactions on Dielectrics
and Electrical Insulation, Vol. 4, No. 1, pp. 10-16, February 1997.
5. Binns, K.J., Lawrenson, P.J., Trowbridge, C.W. The Analytical and Numerical Solution of
Electric and Magnetic Fields. John Wiley & Sons Ltd. 1992.ISBN: 978-0471-92460-9.
8. CFE. Plano de diseño cuerpo básico. 2B2. Torres Mexicanas S.A. de C.V. Feb/87. Comisión
Federal de Electricidad. México. Febrero 1987.
11. CFE. Plano de diseño y montaje poste. AS2. Tomexsa. MAR/99. Comisión Federal de
Electricidad. México. Marzo 1999.
12. CFE. Plano de cuerpo básico y diseño. Torre E71A21. MAY/00. Comisión Federal de
Electricidad. México. Mayo 2000.
13. CFE. Plano de diseño. Torre tipo EA4A22MA. OCT/01. Comisión Federal de Electricidad.
México. Octubre 2001.
15. CFE. NRF-018-CFE. Aisladores tipo suspensión de porcelana o de vidrio templado. Comisión
Federal de Electricidad. México. 2004a.
16. CFE. Poste normalizado tipo 1210SMP. 09/03/04. Comisión Federal de Electricidad. México.
Marzo 2004b.
84
17. CFE. Especificación J1000-50.Torres para líneas de subtransmisión y transmisión. Comisión
Federal de Electricidad. Enero 2006.
23. CIGRE 72. Guidelines for the Evaluation of the Dielectric Strength of External Insulation.
Working Group 07 of Study Committee 33. Conseil International des Grands Réseaux
Électriques. 1992.
24. CIGRE 348. Tower Top Geometry and Mid Span Clearances. Working Group B2.06. Conseil
International des Grands Réseaux Électriques. June 2008.
25. Chakravorti, S., and Steinbigler, H. Boundary-element studies on insulator shape and electric
field around HV insulators with or without pollution. IEEE Transactions on Dielectrics and
Electrical Insulation, Vol. 7, No. 2, April 2000, pp 169-176.
26. Edison Electric Institute. EHV Transmission Line Reference Book. General Electric Company.
1968.
29. EPRI. Transmission Line Reference Book. 345 kV and Above, 2nd Edition. Electric Power
Research Institute, 1982.
30. EPRI. EPRI AC Transmission Line Reference Book – 200 kV and Above, 3rd Edition. Electric
Power Research Institute, December 2005.
31. Espino-Cortes, F.P. y Ramírez Vázquez, I. Uso de materiales con conductividad no-lineal para
atenuar la intensidad del campo eléctrico en aisladores no-cerámicos. IEEE Sección México.
Reunión de Verano de Potencia. RVP-AI/2009-TRA-14-P225. Acapulco, Gro. 5-11 julio 2009.
32. Feser, K., and Pfaff, W. A Potential Free Spherical Sensor for the Measurement of Transient
Electric Fields. IEEE Transactions on Power Apparatus and Systems. Vol. PAS-103, No. 10,
October 1984. pp 2904-2911.
85
33. Fink, D.G., y Beaty, H.W. Manual de Ingeniería Eléctrica, Tomo I. Decimotercera Edición.
McGraw-Hill. México. 1996.ISBN: 970-10-0872-3.
36. Guile, A. E. and Paterson,W. Electrical Power Systems. Volume One. 2nd Edition in SI/Metric
Units. Pergamon International Library. UK. 1977. ISBN: 978-0080-21729-1.
37. Haddad, A., and Warne, D.F. Advances in High Voltage Engineering. The Institution of
Electrical Engineers. UK. 2004. ISBN: 978-0852-96158-2.
38. Hewitt, P. Física Conceptual. Segunda Edición. Addison-Wesley Iberoamericana S.A. USA.
1995.ISBN: 0-201-62595-4.
39. Hileman, Andrew R. Insulation Coordination for Power Systems. Marcel Dekker, Inc. USA.
1999. ISBN: 0-8247-9957-7.
40. Houle, S.; Ghannoum, E.; Hardy, C. Static and dynamic testing of transmission lines
subjected to real wind conditions. CIGRE Symposium: Compact overhead transmission lines.
Leningrad 1991, Report 2000-2.
41. Hyatt, W.H., and Buck, J. Engineering Electromagnetics. Sixth Edition. McGraw-Hill. 2001.
ISBN: 978-0071-20229-9.
42. IEC Std. 60071-1. Insulation co-ordination – Part 1 Definitions, Principles, and Rules. Seventh
Edition. International Electrotechnical Commission. 1993.
43. IEC Std. 60071-2. Insulation co-ordination – Part 2Application Guide. Third Edition.
International Electrotechnical Commission. 1996.
44. IEC Std. 60826. Design criteria of overhead transmission lines. Ed. 3.0. International
Electrotechnical Commission. 2003.
45. IEEE Std. 1313.1-1996. IEEE Standard for Insulation Coordination – Definitions, Principles,
and Rules. The Institute of Electrical and Electronics Engineers, Inc. 1996.
46. IEEE Std. 1313.2-1999. IEEE Guide for the Application of Insulation Coordination. The
Institute of Electrical and Electronics Engineers, Inc. 1999.
47. IEEE Taskforce on Electrical Fields and Composite Insulators. Electric Fields on AC Composite
Transmission Line Insulators. IEEE Transactions on Power Delivery, Vol. 23, No. 2, April
2008. Pages 823-830.
48. Ilhan, S. and Ozdemir, A. Effects of corona ring design on the electric field intensity and
potential distribution along an insulator string. ELECO 2007, 5th International Conference on
Electrical and Electronics Engineering. 5-12 December 2007, Bursa, Turkey.2007.
49. Izgi, E., Inan, A., and Ay, S. The Analysis and Simulation of Voltage Distribution over String
Insulators using Matlab/Simulink. Electric Power Components and Systems, 36:2, 109 – 123.
86
2008.
50. Khalifa, M. High-Voltage Engineering. Theory and Practice. Marcel Dekker, Inc. 1990. ISBN:
978-0824-78128-6.
51. Kiessling, F., Nefzger, P., Nolasco, J.P., and Kaintzyk, U. Overhead Power Lines. Planning,
Design, Construction. Springer-Verlag. 2003. ISBN: 978-3540-00297-0.
52. Kontargyri V.T., Gonos I.F., Ilia N.C., Stathopulos I.A. Simulation of the electric field on
composite insulators using the finite elements method. WSEAS Transactions on Circuits and
Systems, Issue 5, Vol. 3, July 2004, pp. 1318-1322.2004a.
53. Kontargyri V.T., Gonos I.F., Ilia N.C., Stathopulos I.A. Electric Field and Voltage Distribution
along Insulators under Pollution Conditions. Proceedings of the 4th Mediterranean IEE
Conference and Exhibition on Power Generation, Transmission, Distribution and Energy
Conversion (Med Power 2004), Lemesos, Cyprus, November 15-17, 2004. 2004b.
54. Kontargyri V.T., Gonos I.F., Stathopulos I.A., Michaelides A.M. Measurement and
verification of the voltage distribution on high voltage insulators. Proceedings of the 12th
Biennial IEEE Conference on Electromagnetic Field Computation (CEFC 2006), Miami,
Florida, USA, April 30 – May 3, 2006. 2006a.
55. Kontargyri V.T., Gonos I.F., Stathopulos I.A., Michaelides A.M. Simulation of
the electric field on high voltage insulators using finite element method. Proceedings of the
12th Biennial IEEE Conference on Electromagnetic Field Computation (CEFC 2006), Miami,
Florida, USA, April 30 – May 3, 2006. 2006b.
56. Kontargyri V.T., Plati L.N., Gonos I.F., Stathopulos I.A., Michaelides A.M. Measurement and
simulation of the voltage distribution on an insulator string. Proceedings of the 15th
International Symposium on High-Voltage Engineering (ISH 2007), Ljubljana, Slovenia,
August 27-31, 2007, (T1-272).
57. Kontargyri V.T., Plati L.N., Gonos I.F., Stathopulos I.A. Measurement and simulation of the
voltage distribution and the electric field on a glass insulator string. Measurement, Vol. 41,
Issue 5, pp. 471-480. 2008.
59. Marungsri, B., and Onchantuek, W. Simulation of Electric Field and Potential Distributions
on Silicone Rubber Polymer Insulators under Contamination Conditions Using Finite Element
Method. WSEAS Transactions on Power Systems, Issue 9, Volume 3, September 2008. Pages
608-621.
61. Maxwell.http://www.ansys.com/Products/Simulation+Technology/Electromagnetics/Electr
omechanical+Design/ANSYS+Maxwell. Fecha de consulta: 10 de marzo de 2012.
62. Mehdi Ashouri, Mohammad Mirzaie & Ahman Gholami. Calculation of voltage distribution
87
along porcelain suspension insulators based on finite element method. Electric Power
Components and Systems. 38:7, 820-831. UK. 2010.
68. Onchantuek, W. and Marungsri, B. Comparison of Electric Field and Potential Distribution on
Silicone Rubber Polymer Insulators under clean and various Contamination Conditions using
Finite Element Method. WSEAS Transactions on Power Systems, Issue 2, Volume 4, February
2009.
71. Que, W. Electric Field and Voltage Distributions along Non-Ceramic Insulators. Ph. Degree
Theses. The Ohio State University. USA. 2002.
73. Rasolonjanahary, J.L. et al. Computation of Electric Fields and Potential on Polluted
Insulators using a Boundary Element Method. IEEE Transactions on Magnetics. Vol. 28, No.
2. 1992. pp. 1473-1476.
74. Schumann, U., Barcikowski, F., Schreiber, M., and Karner, H.C. FEM Calculation and
Measurement of the Electric Field Distribution of HV Composite Insulator Arrangements.
39th CIGRE Session, Paris, France, 25-30 August, 2002.
75. Still, Alfred. Electric Power Transmission. Third Edition, Seventh Impression. McGraw-Hill
Book Company, Inc. USA. 1927. ISBN: 978-1271-54229-1.
76. SubbaR., and Udaya Kumar. Enhancement of surface flashover performance of high voltage
ceramic disc insulators. Journal of Materials Engineering and Performance. Volume 20.
February 2011. JMEPEG (2011) 20:24-30. pp. 24-30. ASM International.
88
78. Zhao, T., and Comber, G. Calculation of Electric Field and Potential Distribution along
Nonceramic Insulators Considering the Effects of Conductors and Transmission Towers. IEEE
Transactions on Power Delivery, Vol. 15, No. 1, January 2000. pp. 313-318.
79. Zhou, P.B. Numerical Analysis of Electromagnetic Fields. Springer-Verlag. Berlin. 1993. ISBN:
0-387-54722-3.
80. Zienkiewicz, O. C., Taylor, R. L., and Zhu, J. Z. The Finite Element Method: Its Basis and
Fundamentals. Sixth edition. Elsevier Butterworth-Heinemann. UK. 2005. ISBN: 978-0750-
66320-5.
89
APENDICE A
Las características del equipo de cómputo empleado para el desarrollo de los casos de
estudio del presente trabajo son: Procesador Intel Core i3 a 2,13 GHz; Sistema de 64 bits; Memoria
4 GB RAM; Sistema operativo Windows 7 Home Premium. La Tabla A.1 muestra las características
de los casos de estudio realizados en cadenas de aisladores mediante el FEM.
Tabla A.1. Características de los casos de estudio realizados en cadenas de aisladores mediante el FEM.
Tiempo de solución
Caso de estudio Número de elementos Grados de libertad
[seg.]
0° 128794 171230 51,17
Poste
15° 202361 268966 97,313
1210SMP
30° 202000 268479 96,814
0° 208978 277816 99,216
Torre
15° 336134 448265 227,011
E71A21
30° 336632 448923 227,323
0° 112995 128856 84,084
Poste
15° 114064 150187 58,157
AS2
30° 113063 148901 71,261
0° 183292 243116 99,091
Torre
15° 183696 243464 99,872
2B2
30° 184179 244026 98,826
Torre EA4A22MA 33086 44064 7,862
Poste 1210SMP
127084 169029 49,312
(conductor 4 metros)
Poste 1210SMP
4749 6494 0,593
(frontera 20 metros)
Poste 1210SMP
4388 5993 0,5
(frontera 16 metros)
Poste 1210SMP
4424 6057 0,499
(frontera 12,8 metros)
Poste 1210SMP
4285 5897 0,437
(frontera 10,24 metros)
Poste 1210SMP
4146 5739 0,452
(frontera 8,192 metros)
90