Tarea computacional

Campos Electromagnéticos ELI-215

Guillermo Huerta Soto Rafael Möller Montero

201423061-3 201423012-5

Resumen—En el presente informe se aborda el problema de

incorporar en un arreglo de aisladores un electrodo de geometrı́a
toroidal, para lo cual se nos encargó diseñar la configuración
en base a una restricción matemática que llamamos ı́ndice de
desempeño y que se utiliza para buscar la mejor disposición
espacial del toroide respecto al arreglo de aisladores, con el fin
de homogeneizar las caı́das de tensión en estos.

I. I NTRODUCCI ÓN

El problema expuesto consta de un arreglo de 15 aisladores,

el cual utilizaremos en lı́neas de transmisión de 230 kVLLrms .
Figura 1. Base de diseño.
Los materiales que se incorporan para diseñar dichos acce-
sorios son los que se detallan en la siguiente tabla de datos:
Al resultado del trazado computacional, se le suavizaron los
bordes y además, se le redujeron las tolerancias para disminuir
el detalle del enmallado y por ende el coste computacional.
Material Permeavilidad relativaConductividad eléctrica [S/m]
Porcelana 6 2 · 10−13 El resultado del dibujo en Femm es el siguiente:
Cemento 14 1 · 10−13
Aleación de acero 1 4,03 · 106
Aire 1 0
Cuadro I
M ATERIALES DEL AISLADOR .

Sumado a esta configuración se incorpora un electrodo de

geometrı́a toroidal, que se sitúa de tal manera que encierre en
una posición dada a la lı́nea de 15 aisladores. La restricción
de diseño para este aparato es que tenga que ser construido
con 1 metro cuadrado de una lámina de aluminio. La manera
en la cual, en un primer momento, se aborda el problema es
con la utilización de programas de modelación matemática y
de diseño geométrico, es decir, con los programas Femm y
AUTOCAD, respectivamente.
Figura 2. Base de diseño computacional.

Mientras que para el diseño del aislador toroidal se incor-

II. G EOMETR ÍA DEL PROBLEMA pora una circunferencia con cierto radio, distante al eje de
simetrı́a, el cual se presentará de manera mas detallada en
Se nos presenta la imagen de un aislador tipo A, el cual la explicación posterior para los parámetros utilizados en la
con la geometrı́a mostrada es difı́cil dibujar con precisión en optimización. La figura con el toroide incorporado es de la
el programa Femm. Por ello es que se hizo relevante aprender siguiente forma:
a utilizar el programa de diseño geométrico AUTOCAD. Se
rediseña la imagen para obtener las dimensiones H=146[mm],
D=255[mm] y L=295 [mm], que se muestran en la siguiente
figura.
 III. M ODELO PARAM ÉTRICO PARA EL AN ÁLISIS
BIDIMENSIONAL

Se realiza un análisis bidimensional axisimétrico; los

parámetros que conectan los programas son geométricos, los
cuales se van variando y se relacionan con la condición del
área del toroide. El método iterativo utilizado, y mostrado
Figura 3. Toroide con respecto al eje se simetrı́a. como código 1 en el anexo, señala en un primer momento
valores constantes, que son puntos en los que se deben medir
el potencial. Además de ello se muestra que a la entrada del
Luego, haciendo uso de la geometrı́a de los aisladores, se primer ciclo for se entra con un valor de altura respecto al
hizo un modelo equivalente para el dibujo de los 15 aisladores. primer aislador, para luego entrar en el segundo ciclo for con
En este punto se presenta el problema de incorporar los el valor del radio mayor del toroide, y luego se calcula el
elementos que son parte de lineas de transmisión como cables valor del radio menor. la fórmula que relaciona las áreas es la
de fases, neutro y de guardia, además de la misma torre siguiente:
y el suelo, pero dada la opción de asignar condiciones de Atoroide = 4π 2 rR (1)
borde, a nuestro problema, podemos respresentar parte de esos
Donde r es el radio menor y R es el radio mayor.
elementos con esa opción del programa en nuestro análisis.
Siendo esa la parte de desarrollo en MatLab para entrar en
El enmallado se realiza sobre un cuarto de circunferencia, Femm, detallamos el modelo a desarrollar por el programa
donde en el arco se asigna la condición de borde en el infinito Femm en su proceso de análisis.
de potencial 0 [V]. Para incorporar la condición de borde del El método de elementos finitos o FEM consiste en discreti-
cable a tierra se asigna potencial 0 al tramo superior del cuarto zar el dominio de la solución a un número finito de elementos,
de circunferencia. Para el caso de la condición de borde que obteniéndose las ecuaciones que predominan sobre el elemento
impone el cable a potencial asignamos 230 √
3
[KVLL ] en el de estudio en cuestión, para luego ensamblar los diferentes
material inferior del aislador que tiene contacto con la fase, elementos del dominio de la solución para finalmente resolver
además de llevar a ese mismo potencial al toroide. De forma un sistema de ecuaciones. Además, se requiere la información
ilustrativa la geometrı́a del contorno donde se desarrolla el de las condiciones de borde las cuales serán mencionadas al
problema es la siguiente. realizar el modelo equivalente.

IV. I NDICADORES DE DESEMPE ÑO

Para este punto se pensó en tener una referencia a la elécción
del mejor punto de operación, esa fue obtener los parámetros
concentrados previo a la utilización del toroide.

IV-A. Modelo equivalente

Representar nuestro problema de manera matricial nos per-
mite visualizar la relación lineal que tenemos del sistema, si
es que los vemos como electrodos, pese a que la naturaleza
no es por ningún motivo lineal. Por ello es que es importante
mencionar que al analizar el sistema con parámetros concen-
trados, la solución es del tipo exacta y por ende representa de
manera aproximada nuestro conjunto.
Lo dicho anteriormente se reduce en la siguiente relación
matricial:
[q] = [C] [v] (2)
Para la obtención de este modelo equivalente, la tarea de
calcular las capacitancias no es trivial y requiere, nuevamente,
la conexión entre MatLab y Femm para hacer la tarea lo más
Figura 4. Geometrı́a completa del problema. simplificada posible. Para ello se utilizo el código 2 mostrado
en el anexo, el cual hace el siguiente procedimiento:
En un primer momento se asigna a cada electrodo un
Se selecciona simetrı́a axisimétrico para el cálculo de las grupo, a los cuales se le van cambiando las condiciones
variables a analizar posteriormente, por ello es que se justifica de potencial. Estos se definen como conductores para
el trabajo en solo un cuarto de circunferencia lo cual permite lograr obtener, mediante comandos, la carga acumulada
ahorrar tiempo de iteración y esfuerzo del software. en ellos.
 Para una primera iteración, se lleva a potencial 1[V] a un a una altura h=337.46[mm], donde la coordenada (0,0) se
electrodo definido como conductor perteneciente a cierto define en la parte inferior del aislador 1.
grupo, y los demás electrodos a potencial 0 [V]. Las dimensiones se muestran en la siguiente figura:
Luego de ello se extrae el dato de la carga almacenada en
cada electrodo, posteriormente se relacionan las variables
tensión y carga para representar la interacción entre estas
en una capacitancia.
Teniendo esa primera columna de la matriz de capaci-
tancia, se procede a elevar el potencial de otro electrodo,
aterrizando el potencial del electrodo antes energizado
mediante un proceso iterativo.
Haciendo el procedimiento anterior en todos los electrodos se
obtiene la matriz de capacitancia sin el electrodo toroidal que
se entrega en el anexo.
Con el toroide tenemos la matriz de capacitancia que se
muestra también en el anexo. Ambos modelos se compararán
en el siguiente apartado, donde, de manera gráfica, se distingue
una diferencia en las curvas de potencial de la figura 6.
Al realizar una comparación entre los dos casos, se tiene que
al agregar el toroide en la configuración las capacitancias
aumentan, lo cual concuerda con la teorı́a donde se busca Figura 5. Posición óptima.
mejorar la distribución de potencial y reducción de la exigencia
de la estructura de aisladores, una mayor capacitancia implica Los resultados obtenido de las diferencias de potencial para
una mayor capacidad para almacenar carga a una misma los aisladores son las siguientes:
tensión, lo cual mejora el desempeño en aislación.
Aislador ∆Vi [Vrms ]
IV-B. Índice de desempeño de la configuración 1 10311
2 8699
Se nos pide evaluar el siguiente ı́ndice: 3 7803
15  2 4 9083
X Vnom 5 9739
I= ∆Vi − (3)
i∈
|N | 6 9676
7 9723
El cual mediante el código 1, antes mencionado, se evalúa de 8 9006
9 8892
la siguiente forma: 10 8414
En un primer momento se toma el potencial entre los 11 8008
12 8183
puntos extremos de un aislador.
13 7977
Se guarda en una matriz. 14 8947
Se evaluán las diferencias de tensión y se realiza la 15 8124
sumatoria antes descrita. Cuadro II
D IFERENCIAS DE POTENCIAL EN AISLADORES
Este procedimiento se hace de manera iterativa para cada
punto en el que vamos situando el toroide. El menor ı́ndice de
desempeño para nuestra configuración, es de 8210765 (orden Datos que comparamos con los obtenidos al trabajar sin
de magnitud 6). El cual se obtuvo en tres etapas de mapeo: toroide en la siguiente gráfica:
En la primera se hizo un proceso iterativo grueso de 2100
puntos en un área de 1.204 [m2 ]. 104
En la segunda se hizo un proceso iterativo fino de 1204 13.2791
12.2507
Con toroide

puntos en un área de 0.301 [m2 ]. 11.3807

10.6004
Sin toroide

En la tercera se hizo un proceso iterativo fino de 1360 9.6922

Potencial [V]

8.7183
puntos en un área de 0.00136 [m2 ]. 7.7507
6.7784
5.8778
V. C ONFIGURACI ÓN ÓPTIMA 4.9886
4.1472
3.3464
La mejor configuración incorpora la geometrı́a antes 2.5281
1.7305
descrita para el conjunto de aisladores, además de los 0.8358
0
0 145 289 434 579 723 868 1013 1158 1302 1447 1592 1736 1881 2026 2170
parámetros geométricos para el toroide dado su respectivo Altura z [mm]
indice de desempeño que mencionamos anteriormente, que
consta de un valor de; R=483 [mm] y un valor de r=52.44[mm] Figura 6. Potencial en función de altura.
 Figura que muestra el efecto de incorporar el toroide a la Todo ello nos hace inferir que, con la consideración de
configuración, que se traduce en una mejora importante en la que el modelo es altamente no lineal, el análisis que
distribución de potencial a lo largo del arreglo. proponemos es de calidad; donde se tiene un error máximo
de 4,75 % en la zona de mayor dificultad de análisis, debido
V-A. Criterios de selección a que el electrodo 1 incluye al toroide, complicando la
Estos constan de condiciones pedidas por la tarea y otras in- geometrı́a y por ende afectando la resolución numérica. Sin
vestigadas para el diseño que incorpore una elección correcta. embargo en el resto del arreglo se tiene un error menor al
Estas se muestran a continuación: 2,5 %, lo cual está dentro de lo razonable.
Utilización de 1 [m2 ] para la construcción del toroide.
Tensiones nominales para cada aislador. Electrodo qmatriz [C] qf emm [C] ε[ %]
Menor ı́ndice de desempeño obtenido. 1 -5,69E-07 -5,43E-07 -4,75
Distribución lineal de potencial. 2 -9,08E-07 -8,93E-07 -1,64
3 -1,36E-06 -1,33E-06 -2,17
VI. M ODELO SIMPLIFICADO VERSUS MODELO EN FEMM 4 -1,89E-06 -1,85E-06 -2,15
5 -2,56E-06 -2,51E-06 -1,94
Los potenciales que vamos a considerar como parte de 6 -3,40E-06 -3,36E-06 -0.98
nuestra matriz de [v] de la ecuación (2) son los que permitirán 7 -4,55E-06 -4,49E-06 -1,14
validar el modelo de parámetros concentrados. Los valores de 8 -6,20E-06 -6,11E-06 -1,42
9 -8,58E-06 -8,5E-06 -0,93
potenciales escogido son los siguientes: 10 -1,23E-05 -1,21E-05 -1,63
  11 -1,76E-05 -1,71E-05 -2,56
0 12 -2,39E-05 -2,36E-06 -1,02
 0  13 -3,00E-05 -2,97E-06 -0,78
 
 0  14 -3,73E-05 -3,75E-05 0,73
15 -2,5E-03 -2,5E-03 -1,02
 
 0 

 0 
 16 2,7E-03 2,7E-03 -0,44
  Cuadro III
 0  C OMPARACI ÓN DE PAR ÁMETROS CONCENTRADOS .
 
 0 
 
 0  VII. F EN ÓMENOS Y EFECTOS DESPRECIADOS
[v] =  0 

 [V ] (4)
  En las lineas de transmisión se tienen fenómenos que vienen
 0 
  asociados a la presencia de los campos eléctricos y magnéticos
 0 
  provocados por el flujo de potencia.
 0 

 0 
 Uno de ellos es el efecto capacitivo que provoca la tierra

 0 
 con los elementos estudiados, que pese a considerarla en
  las condiciones de borde, el modelo capacitivo tambien
 0 
230000 podrı́a incorporar la interaccion entre los accesorios
√
3 estudiados y la misma tierra.
Los cuales son los potenciales del electrodo superior hasta el Se despreció el efecto Ferranti que es la sobretensión que
inferior. se provoca en las lineas de transmisión producto del nivel
Para contrastar el resultado anterior, obtenemos las cargas de demanda que se tenga en la carga. En nuestro análisis
para dicho potenciales mediante la utilización del sistema que no se trabajo en rangos de tensión, que serı́a algo ideal
anteriormente citamos (2). Se usa la matriz de capacitancia del y relevante a a la hora de diseñar los aisladores de tal
cuadro (V) para resolver el sistema, el cual es ejecutado por forma que puedan soportar todos los cambios de tensión
el programa MatLab. Los resultados arrojados, en contraste en los que trabajan las lineas de transmisión.
con las cargas reales dadas por Femm, además del cálculo del Otro fenómeno que no fue considerado es la suciedad
error hecho con la siguiente fórmula: que se impregna en los aisladores a medida que pasa el
tiempo, esto provoca que se tenga otro material sobre
qf emm − qmatriz la superficie de los accesorios diseñados, haciendo que
ε= · 100 % (5)
qf emm varı́en las propiedades del aislador y su desempeño. [1]
Se muestran en el cuadro III. Además, los cambios de clima, según la estación, tam-
Los errores se justifican debido a que estamos frente a un bién cambian las propiedades del aire, como es la hu-
modelo lineal, siendo entonces una aproximación de nuestro medad; parámetro a considerar si queremos diseñar el
sistema, ya que entrega soluciones del tipo exacta. Pese a ello arreglo para todas las estaciones del año, donde los
los resultados obtenidos son bastante cercanos a lo que el aisladores se verı́an sometidos a un estrés adicional al
programa computacional entrega, debido a que el análisis se considerar esta idea, siendo por ejemplo la temperatura
sustenta en las leyes de Maxwell, que nos permiten relacionar en verano o la humedad en un dı́a lluvioso un factor
la tensión con la carga en los electrodos. influyente en el desempeño de los aisladores que no fue
considerado.
 R EFERENCIAS VIII. A NEXOS
[1] S. M. A. Dhalaan, “Investigation on the characteristics of
a string of insulator due to the effect of dirt,” IEEE, and
M. A. Elhirbawy, 2003.
[2] S. M. A. Dhalaan, “Simulation of voltage distribution cal-
culation methods over a string of suspension insulators,”
IEEE, and M. A. Elhirbawy, 2003.
[3] A. J. Phillips, “Electric fields on ac composite transmission
line insulators,” IEEE Taskforce on Electric Fields and
Composite Insulators, 2007.
[4] M. M. A.Rahimnejad, “Optimal corona ring selection for
230 kv ceramic i-string insulator using 3d simulation,”
International Journal of Scientific Engineering Research,
2012.
[5] M. M. Ebrahim Akbari, “Investigating the effects of disc
insulator type and corona ring on voltage distribution over
230- kv insulator string using 3-d fem,” Dept. of Elec. and
Comp. Eng., Babol Univ. of Tech., Babol, Iran, 2012.
[6] S. Karunaratne, “Linear analysis of texture-property rela-
tionships using process-based representations of rodrigues
space,” A New Approach to the Study of Voltage Distribu-
tion along a Suspension Insulator String, 2017.
VIII-A. Códigos %R e s u l t a d o s
VIII-A1. Código 1 principal para optimizar: I m i n =min ( min ( I ) )
[ c o l u m n a o p t i m a , f i l a o p t i m a ] = f i n d ( I ==min ( min ( I ) ) )
% %P r o c e s o Roptimo =RR( c o l u m n a o p t i m a , f i l a o p t i m a )
Zoptimo =ZZ ( c o l u m n a o p t i m a , f i l a o p t i m a )
I = [ ] ; %M a t r i z que g u a r d a l o s ndices
de l a s i t e r a c i o n e s .
c = 0 ; %C o n t a d o r .
r o = 1 4 2 0 2 . 5 ; %p u n t o en e l e j e de a r r e g l o de a i s l a d o r e s .
z o i = 3 5 3 8 . 5 ; %P u n t o i n f e r i o r d e l p r i m e r a i s l a d o r .
z o s = 3 7 0 2 . 5 ; %P u n t o s u p e r i o r d e l p r i m e r a i s l a d o r .
V = [ ] ; %D i f e r e n c i a s de p o t e n c i a l .
h = 3 8 4 7 . 2 − 3 7 0 2 . 5 ; %A l t u r a de un a i s l a d o r .
Vnom = 2 3 0 0 0 0 / 3 ˆ 0 . 5 ; %T e n s i n n o m i n a l .
RR = [ ] ; %M a t r i z que g u a r d a r a d i o s .
ZZ = [ ] ; %M a t r i z que g u a r d a a l t u r a s .
cuenta =0;
f o r z = 3 8 5 4 : 1 : 3 8 8 8 ; %V a l o r v e r t i c a l de n o d o s
c=c + 1 ;
f =0;
f o r R= 4 6 0 : 1 : 5 0 0 ;
c u e n t a = c u e n t a +1
openfemm
o p e n d o c u m e n t ( ’ modelo . FEE ’ ) ;
e i s a v e a s ( ’ temp . f e e ’ ) ;
f = f +1;
RR( c , f ) =R ;
ZZ ( c , f ) = z ;
r = ( 1 * 1 0 ˆ 6 ) / ( 4 * p i ˆ 2 *R ) ;
r 2 = ro−R+ r ; %V a l o r h o r i z o n t a l de nodo d e r e c h o .
r 1 = ro−R−r ; %V a l o r h o r i z o n t a l de nodo i z q u i e r d o .
e i a d d n o d e ( r1 , z ) ; %Nodo i z q u i e r d o
e i a d d n o d e ( r2 , z ) ; %Nodo d e r e c h o
e i a d d b l o c k l a b e l ( ( r 1 + r 2 ) / 2 , z ) ; %C a r a c t e r s t i c a de m a t e r i a l
e i a d d a r c ( r1 , z , r2 , z , 1 8 0 , 1 ) ; %Arco g e n e r a d o s u p e r i o r
e i a d d a r c ( r2 , z , r1 , z , 1 8 0 , 1 ) ; %Arco g e n e r a d o i n f e r i o r
e i s e l e c t a r c s e g m e n t ( ( r 1 + r 2 ) / 2 , z+ r ) ;
e i s e l e c t a r c s e g m e n t ( ( r 1 + r 2 ) / 2 , z−r ) ;
e i s e t a r c s e g m e n t p r o p ( 1 , ’ None ’ , 0 , 4 , 2 3 0 ) ;
ei clearselected ;
e i s e l e c t l a b e l ( ( r 1 + r 2 ) / 2 , z ) %s e l e c c i o n de l a b e l ;
e i s e t b l o c k p r o p ( ’ <No Mesh > ’ , 0 , 0 , 1 ) ;
ei clearselected ;
ei analyze ;
ei loadsolution ;
% Postproceso
I I = [ ] ; %M a t r i z que g u a r d a l a s d i f e r e n c a s c u a d r t i c a s .
for i =1:1:15;
V( i ) = e o g e t v ( ro , ( zos−h ) + ( i −1) * h)− e o g e t v ( ro , z o s + ( i −1) * h ) ; %D i f e r e n c i a de t e n s i
end
for i =1:1:15;
I I ( i ) = (V( i )−Vnom / 1 5 ) ˆ 2 ; %E l e m e n t o i de l a s u m a t o r i a .
end
I ( c , f ) = sum ( I I ) ; %S u m a t o r i a de c a s o n .
closefemm
end
end
 VIII-A2. Código 2 para cálculo de capacitancias: ei clearselected

C=[]; ei selectgroup (18)

for n =1:1:16 ei setnodeprop ( ’18 ’ ,0 , ’18 ’)
m=n + 9 ; ei setsegmentprop ( ’18 ’ ,1 ,0 ,0 ,18 , ’18 ’)
openfemm ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 8 , ’ 1 8 ’ )
opendocument ( ’ e l e c t r o d o s s i n t o r o i d e . fee ’ ) ; ei clearselected
e i s a v e a s ( ’ temp . f e e ’ ) ;
e i m o d i f y c o n d u c t o r p r o p ( n u m 2 s t r (m) , 1 , 1 ) ei selectgroup (19)
ei setnodeprop ( ’19 ’ ,0 , ’19 ’)
ei selectgroup (10) ei setsegmentprop ( ’19 ’ ,1 ,0 ,0 ,19 , ’19 ’)
ei setnodeprop ( ’10 ’ ,0 , ’10 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 9 , ’ 1 9 ’ )
ei setsegmentprop ( ’10 ’ ,1 ,0 ,0 ,10 , ’10 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 0 , ’ 1 0 ’ )
ei clearselected ei selectgroup (20)
ei setnodeprop ( ’20 ’ ,0 , ’20 ’)
ei selectgroup (11) ei setsegmentprop ( ’20 ’ ,1 ,0 ,0 ,20 , ’20 ’)
ei setnodeprop ( ’11 ’ ,0 , ’11 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 0 , ’ 2 0 ’ )
ei setsegmentprop ( ’11 ’ ,1 ,0 ,0 ,11 , ’11 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 1 , ’ 1 1 ’ )
ei clearselected ei selectgroup (21)
ei setnodeprop ( ’21 ’ ,0 , ’21 ’)
ei selectgroup (12) ei setsegmentprop ( ’21 ’ ,1 ,0 ,0 ,21 , ’21 ’)
ei setnodeprop ( ’12 ’ ,0 , ’12 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 1 , ’ 2 1 ’ )
ei setsegmentprop ( ’12 ’ ,1 ,0 ,0 ,12 , ’12 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 2 , ’ 1 2 ’ )
ei clearselected ei selectgroup (22)
ei setnodeprop ( ’22 ’ ,0 , ’22 ’)
ei selectgroup (13) ei setsegmentprop ( ’22 ’ ,1 ,0 ,0 ,22 , ’22 ’)
ei setnodeprop ( ’13 ’ ,0 , ’13 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 2 , ’ 2 2 ’ )
ei setsegmentprop ( ’13 ’ ,1 ,0 ,0 ,13 , ’13 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 3 , ’ 1 3 ’ )
ei clearselected ei selectgroup (23)
ei setnodeprop ( ’23 ’ ,0 , ’23 ’)
ei selectgroup (14) ei setsegmentprop ( ’23 ’ ,1 ,0 ,0 ,23 , ’23 ’)
ei setnodeprop ( ’14 ’ ,0 , ’14 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 3 , ’ 2 3 ’ )
ei setsegmentprop ( ’14 ’ ,1 ,0 ,0 ,14 , ’14 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 4 , ’ 1 4 ’ )
ei clearselected ei selectgroup (24)
ei setnodeprop ( ’24 ’ ,0 , ’24 ’)
ei selectgroup (15) ei setsegmentprop ( ’24 ’ ,1 ,0 ,0 ,24 , ’24 ’)
ei setnodeprop ( ’15 ’ ,0 , ’15 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 4 , ’ 2 4 ’ )
ei setsegmentprop ( ’15 ’ ,1 ,0 ,0 ,15 , ’15 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 5 , ’ 1 5 ’ )
ei clearselected ei selectgroup (25)
ei setnodeprop ( ’25 ’ ,0 , ’25 ’)
ei selectgroup (16) ei setsegmentprop ( ’25 ’ ,1 ,0 ,0 ,25 , ’25 ’)
ei setnodeprop ( ’16 ’ ,0 , ’16 ’) ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 2 5 , ’ 2 5 ’ )
ei setsegmentprop ( ’16 ’ ,1 ,0 ,0 ,16 , ’16 ’) ei clearselected
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 6 , ’ 1 6 ’ )
ei clearselected
ei analyze ;
ei selectgroup (17) ei loadsolution ;
ei setnodeprop ( ’17 ’ ,0 , ’17 ’)
ei setsegmentprop ( ’17 ’ ,1 ,0 ,0 ,17 , ’17 ’) for x =1:1:16
ei s e t a r c s e g m e n t p r o p (1 , ’ < None > ’ , 0 , 1 7 , ’ 1 7 ’ ) Datos= e o g e t c o n d u c t o r p r o p e r t i e s ( x +9)
 C( n , x )= Datos ( 2 ) ; VIII-B. Capacitancias
end
closefemm
end
 1,87E-08 -1,84E-08 -8,43E-11 -2,86E-11 -1,50E-11 -9,73E-12 -6,84E-12 -5,08E-12 -3,87E-12 -3,06E-12 -2,36E-12 -1,78E-12 -1,35E-12 -9,70E-13 -6,59E-13 -4,23E-13
-1,84E-08 3,75E-08 -1,87E-08 -1,24E-10 -4,01E-11 -2,14E-11 -1,37E-11 -9,61E-12 -7,06E-12 -5,52E-12 -4,17E-12 -3,10E-12 -2,30E-12 -1,63E-12 -1,10E-12 -7,05E-13
-8,79E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,01E-10 -2,89E-11 -1,42E-11 -8,78E-12 -5,93E-12 -4,43E-12 -3,20E-12 -2,28E-12 -1,65E-12 -1,14E-12 -7,67E-13 -5,01E-13
-3,02E-11 -1,23E-10 -1,87E-08 3,78E-08 -1,87E-08 -1,04E-10 -3,03E-11 -1,49E-11 -8,97E-12 -6,19E-12 -4,28E-12 -2,99E-12 -2,13E-12 -1,47E-12 -9,72E-13 -6,12E-13
-1,48E-11 -3,72E-11 -1,04E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,80E-11 -2,77E-11 -1,33E-11 -8,21E-12 -5,32E-12 -3,54E-12 -2,43E-12 -1,63E-12 -1,05E-12 -6,77E-13
-1,00E-11 -2,09E-11 -3,08E-11 -1,03E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,94E-11 -2,83E-11 -1,43E-11 -8,34E-12 -5,20E-12 -3,43E-12 -2,25E-12 -1,44E-12 -9,08E-13
-6,99E-12 -1,32E-11 -1,48E-11 -2,85E-11 -9,92E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,99E-11 -2,95E-11 -1,39E-11 -7,69E-12 -4,79E-12 -3,00E-12 -1,89E-12 -1,17E-12
-5,14E-12 -9,14E-12 -9,00E-12 -1,38E-11 -2,74E-11 -9,70E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,01E-10 -2,87E-11 -1,27E-11 -7,15E-12 -4,27E-12 -2,58E-12 -1,58E-12
-4,00E-12 -6,96E-12 -6,33E-12 -8,69E-12 -1,38E-11 -2,86E-11 -1,00E-10 -1,87E-08 3,78E-08 -1,87E-08 -1,08E-10 -2,90E-11 -1,29E-11 -6,98E-12 -4,02E-12 -2,35E-12
-3,01E-12 -5,01E-12 -4,27E-12 -5,50E-12 -7,81E-12 -1,31E-11 -2,74E-11 -9,65E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,01E-10 -2,58E-11 -1,12E-11 -5,89E-12 -3,37E-12

Cuadro IV
-2,38E-12 -3,92E-12 -3,24E-12 -3,98E-12 -5,34E-12 -8,04E-12 -1,36E-11 -2,85E-11 -1,03E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,75E-11 -2,56E-11 -1,09E-11 -5,66E-12
-1,85E-12 -3,04E-12 -2,44E-12 -2,90E-12 -3,69E-12 -5,26E-12 -7,95E-12 -1,33E-11 -2,91E-11 -1,06E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,60E-11 -2,50E-11 -1,06E-11

C APACITANCIA SIN T OROIDE

-1,42E-12 -2,33E-12 -1,83E-12 -2,15E-12 -2,69E-12 -3,59E-12 -5,10E-12 -7,73E-12 -1,34E-11 -2,95E-11 -1,03E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,39E-11 -2,47E-11
-1,02E-12 -1,60E-12 -1,23E-12 -1,41E-12 -1,71E-12 -2,26E-12 -3,09E-12 -4,42E-12 -6,95E-12 -1,25E-11 -2,67E-11 -9,76E-11 -1,87E-08 3,78E-08 -1,86E-08 -9,60E-11
-6,87E-13 -1,08E-12 -8,14E-13 -9,24E-13 -1,10E-12 -1,43E-12 -1,92E-12 -2,65E-12 -3,96E-12 -6,47E-12 -1,13E-11 -2,44E-11 -9,14E-11 -1,86E-08 3,78E-08 -1,87E-08
-4,12E-13 -6,76E-13 -5,13E-13 -5,80E-13 -6,93E-13 -8,87E-13 -1,17E-12 -1,58E-12 -2,28E-12 -3,54E-12 -5,70E-12 -1,02E-11 -2,29E-11 -9,37E-11 -1,87E-08 2,00E-08

2,02E-08 -1,86E-08 -2,82E-10 -2,22E-10 -1,76E-10 -1,28E-10 -9,01E-11 -6,33E-11 -4,53E-11 -3,33E-11 -2,50E-11 -1,86E-11 -1,37E-11 -9,75E-12 -6,59E-12 -4,02E-12
-1,86E-08 3,76E-08 -1,87E-08 -9,54E-11 -1,83E-11 -6,63E-12 -3,79E-12 -2,77E-12 -2,07E-12 -1,66E-12 -1,38E-12 -1,09E-12 -8,19E-13 -6,08E-13 -4,14E-13 -2,73E-13
-2,81E-10 -1,87E-08 3,78E-08 -1,87E-08 -8,54E-11 -1,73E-11 -6,11E-12 -3,26E-12 -1,92E-12 -1,30E-12 -9,67E-13 -7,07E-13 -4,99E-13 -3,58E-13 -2,37E-13 -1,57E-13
-2,21E-10 -9,44E-11 -1,87E-08 3,78E-08 -1,87E-08 -8,89E-11 -1,98E-11 -8,20E-12 -4,19E-12 -2,57E-12 -1,75E-12 -1,19E-12 -8,02E-13 -5,49E-13 -3,56E-13 -2,32E-13
-1,76E-10 -1,89E-11 -8,64E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,25E-11 -2,38E-11 -9,91E-12 -5,50E-12 -3,52E-12 -2,31E-12 -1,51E-12 -1,02E-12 -6,46E-13 -4,07E-13
-1,28E-10 -7,17E-12 -1,89E-11 -9,15E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,00E-10 -2,61E-11 -1,19E-11 -6,95E-12 -4,35E-12 -2,76E-12 -1,82E-12 -1,14E-12 -7,14E-13
-9,04E-11 -4,05E-12 -6,70E-12 -2,07E-11 -9,46E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,02E-10 -2,79E-11 -1,32E-11 -7,51E-12 -4,51E-12 -2,90E-12 -1,78E-12 -1,08E-12
-6,33E-11 -2,89E-12 -3,51E-12 -8,51E-12 -2,46E-11 -9,90E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,73E-11 -2,86E-11 -1,32E-11 -7,22E-12 -4,41E-12 -2,63E-12 -1,58E-12
-4,53E-11 -2,06E-12 -1,94E-12 -4,07E-12 -9,58E-12 -2,43E-11 -9,92E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,78E-11 -2,65E-11 -1,15E-11 -6,47E-12 -3,64E-12 -2,16E-12
-3,33E-11 -1,71E-12 -1,37E-12 -2,57E-12 -5,48E-12 -1,14E-11 -2,68E-11 -9,73E-11 -1,87E-08 3,78E-08 -1,87E-08 -1,04E-10 -2,62E-11 -1,18E-11 -6,08E-12 -3,46E-12

Cuadro V
-2,50E-11 -1,37E-12 -9,76E-13 -1,68E-12 -3,44E-12 -6,49E-12 -1,24E-11 -2,71E-11 -9,81E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,94E-11 -2,59E-11 -1,09E-11 -5,68E-12
-1,86E-11 -1,13E-12 -7,65E-13 -1,22E-12 -2,39E-12 -4,30E-12 -7,44E-12 -1,32E-11 -2,83E-11 -1,02E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,89E-11 -2,59E-11 -1,08E-11

C APACITANCIA CON T OROIDE

-1,37E-11 -8,69E-13 -5,50E-13 -8,36E-13 -1,59E-12 -2,77E-12 -4,53E-12 -7,33E-12 -1,25E-11 -2,75E-11 -1,03E-10 -1,87E-08 3,78E-08 -1,87E-08 -9,86E-11 -2,49E-11
-9,89E-12 -6,12E-13 -3,64E-13 -5,36E-13 -1,00E-12 -1,72E-12 -2,73E-12 -4,22E-12 -6,54E-12 -1,16E-11 -2,62E-11 -9,79E-11 -1,87E-08 3,78E-08 -1,87E-08 -9,63E-11
-6,61E-12 -4,18E-13 -2,43E-13 -3,45E-13 -6,41E-13 -1,07E-12 -1,67E-12 -2,50E-12 -3,71E-12 -5,98E-12 -1,10E-11 -2,45E-11 -9,57E-11 -1,87E-08 3,78E-08 -1,87E-08
-4,14E-12 -2,70E-13 -1,56E-13 -2,19E-13 -4,00E-13 -6,65E-13 -1,02E-12 -1,51E-12 -2,16E-12 -3,33E-12 -5,61E-12 -1,02E-11 -2,29E-11 -9,45E-11 -1,87E-08 2,00E-08