UNIVERSIDAD NACIONAL DE INGENIERIA

Centro de Tecnologı́as de Información y

Comunicaciones

PROGRAMAS DE ALTA ESPECIALIZACION TECNOLÓGICA

Medición de radiaciones electromagnéticas no

ionizantes en CTIC

Alumnos:

Romero Catay Eder Angelo

Ticse Torres Roger

Asesor:

Dr. Carlos Javier Solano Salinas

1
 Resumen

Una de las formas de transmisión de energı́a es la que se realiza a través de

ondas electromagnéticas, estas se diferencian por la cantidad de energı́a que
puedan transmitir y esto a su vez por su frecuencia. En el presente trabajo
daremos una visión teórica rápida de las ondas electromagnéticas(OEM)
que producen radiación y pasaremos a explicar los métodos y procedimientos
con los cuales determinaremos teóricamente la cantidad de radiación no
ionizante presente en las instalaciones de CTIC.
Índice general

I 3

1. Introducción 4
1.1. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2. Fundamento Teórico 5
2.1. El espectro electromagnético . . . . . . . . . . . . . . . . . . 5
2.2. Ecuaciones de Maxwell. . . . . . . . . . . . . . . . . . . . . . 6
2.2.1. Ecuaciones Diferenciales . . . . . . . . . . . . . . . . . 6
2.2.2. Casos particulares de las ecuaciones de Maxwell . . . 6
2.2.3. Ecuaciones en forma integral . . . . . . . . . . . . . . 7
2.3. Ecuaciones de onda para los campos . . . . . . . . . . . . . . 7
2.3.1. Ecuaciones de Maxwell(Variacion armonica) . . . . . . 7
2.3.2. Ecuacion de continuidad . . . . . . . . . . . . . . . . . 7
2.3.3. Ecuación de onda para el campo eléctrico . . . . . . . 8
2.3.4. Ecuación de onda para el campo magnético . . . . . . 8
2.3.5. Definición de los potenciales . . . . . . . . . . . . . . . 8
2.3.6. Ecuaciones de onda de los potenciales . . . . . . . . . 9
2.3.7. Soluciones Integrales para los Potenciales . . . . . . . 9
2.3.8. Expresiones generales de los campos . . . . . . . . . . 9
2.3.9. Campos inducidos . . . . . . . . . . . . . . . . . . . . 10
2.3.10. Campos radiados . . . . . . . . . . . . . . . . . . . . . 10
2.3.11. Expresiones aproximadas para los campos radiados . . 11
2.3.12. Vector de radiacion . . . . . . . . . . . . . . . . . . . . 11
2.4. Antenas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.4.1. Diagrama de radiacion . . . . . . . . . . . . . . . . . . 12
2.4.2. Diagramas bidimencionales . . . . . . . . . . . . . . . 12
2.5. Parametros fundamentales de las antenas . . . . . . . . . . . 13
2.5.1. Densidad de potencia radiada . . . . . . . . . . . . . . 13
2.5.2. Directividad . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5.3. Ganancia . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.5.4. Polarización . . . . . . . . . . . . . . . . . . . . . . . . 15
2.5.5. Impedancia . . . . . . . . . . . . . . . . . . . . . . . . 16
2.5.6. Area y longitud efectivas . . . . . . . . . . . . . . . . 16

3. Esquemas de Regulación 17
3.1. Método de Predicción . . . . . . . . . . . . . . . . . . . . . . 18
3.2. Método de Medición . . . . . . . . . . . . . . . . . . . . . . . 19
3.3. Monitoreo Continuo . . . . . . . . . . . . . . . . . . . . . . . 19
3.4. Mediciones en CTIC . . . . . . . . . . . . . . . . . . . . . . . 20

1
 3.4.1. Antenas de internet inalámbrico(wireless) . . . . . . . 20
3.4.2. Monitores LCD . . . . . . . . . . . . . . . . . . . . . . 21
3.4.3. CPU (los circuitos que lo componen) . . . . . . . . . . 23
3.4.4. Aire acondicionado . . . . . . . . . . . . . . . . . . . . 23
3.5. Cálculo aproximado de radiación . . . . . . . . . . . . . . . . 24

II 26

4. Software para simulación de Antenas 27

4.1. Métodos computacionales . . . . . . . . . . . . . . . . . . . . 31
4.1.1. Método de los momentos (MoM) . . . . . . . . . . . . 31
4.1.2. Método de diferencias finitas en el dominio del tiempo
(FDTD) . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.1.3. El método de los elementos finitos . . . . . . . . . . . 35
4.2. Medición de Campos Magneticos en CTIC . . . . . . . . . . . 36

5. conclusiones 51

A. Medidores de campos electromagneticos 52

A.0.1. Medidor de campo eléctrico PCE-EM 29 . . . . . . . . 52
A.0.2. MEDIDOR DE RADIOFRECUENCIA . . . . . . . . 52
A.0.3. ANALIZADOR DE INTENSIDAD DE CAMPOS RF 53
A.0.4. MEDIDOR DE CAMPO TAOMA . . . . . . . . . . . 55
A.0.5. Analizador de Campos electromagnéticos EMR-300 de
Narda . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
A.0.6. Analizador de Espectros: Agilent HP 8594E . . . . . . 57

Bibliografı́a 62

2
Parte I

3
Capı́tulo 1

Introducción

El mundo en que vivimos tiene la tecnologı́a e información cada vez mas

al alcance de las personas, esto es muy bueno pero en cierta medida ya que
los sistemas electrónicos generan radiación, felizmente la mayorı́a de estos
aparatos emiten radiación no ionizante (RNI) que es menos perfudicial que
la radiación ionizante pero no por eso hay que dejar de lado los posibles
efectos que puedan ocasionar el estar inmersos en ese mar de radiaciones
que nos están llegando a cada instante.
Nosotros en este trabajo sólo nos dedicaremos a presentar algunas con-
sideraciones teóricas y dar un resultado de medida teórica de la potencia
de radiación no ionizante en las instalaciones del Centro de Tecnologı́as de
Investigación y Comunicaciones CTIC.
Este informe está basado en las consideraciones tanto teóricas como
prácticas de un estudio orientado a determinar la densidad de potencia esper-
ada y la medida en forma experimental en las inmediaciones de una antena
de una estación base de tecnologı́a GSM (Global System for Mobile Com-
munications) que fue desarrollado por el Dr. Cesar Javier Solano Salinas y
su tesista Loreleyn Flores

1.1. Objetivo
Lo que nos proponemos en este informe es estudiar y encontrar un
método para determinar teóricamente la potencia de radiación emitida por
fuentes que en nuestro caso serán aparatos electrónicos .
Lo cual servirá para luego determinar esperimentalmente la cantidad de
radiación real a la cual estamos expuestos los integrantes del CTIC.

4
Capı́tulo 2

Fundamento Teórico

2.1. El espectro electromagnético

Las ondas electromagnéticas se caracterizan por su frecuencia y longitud
de onda. El conjunto de todas las frecuencias se denomina espectro.

Figura 2.1: Espectro Electromagnetico.

Las ondas se clasifican por bandas. Las denominaciones de las bandas de

frecuencia se pueden realizar por décadas, como por ejemplo MF, HF, VHF,
UHF.

Figura 2.2: Frecuencias

5
2.2. Ecuaciones de Maxwell.
2.2.1. Ecuaciones Diferenciales
Ls fenomenos electromagneticos se pueden describir a partir de cuatro
ecuaciones de Maxwell.
Ley de Ampere
~
∇×B ~ = µo J~ + µo o ∂ E
∂t
Ley de Faraday
~
∇×E ~ = ∂B
∂t
Ley de Gauss
~ = ρ
∇.E
o
Ley de Gauss
~ =0
∇.B

2.2.2. Casos particulares de las ecuaciones de Maxwell

En el espacio libre las corrientes y las cargas son cero y las ecuaciones
de Maxwell se pueden simplificar eliminando los terminos correspondientes.
Asimismo si las fuentes varı́an armónicamente con el tiempo, las ecuaciones
electromagnéticas y sus soluciones se simplifican, utilizando para ello una
notación fasorial, de forma que las derivadas respecto al tiempo se transfor-
man en productos por el factor jω .

Finalmente para casos sin variación temporal, las ecuaciones toman las
formas de electrostática y magnetostática.

Figura 2.3: Unidades

6
2.2.3. Ecuaciones en forma integral
Las ecuaciones de Maxwell se pueden escribir en forma integral, aplican-
do para ello los teoremas de Stokes y de la divergencia

~
!
∂D
I Z Z
~ ~l =
H.d J~ + .d~s
∂t

~
!
∂B
I Z Z
~ ~l = −
E.d .d~s
∂t

Z Z Z Z Z
D. ~ =
~ ds ρdv

Z Z
~ =0
~ ds
B.

2.3. Ecuaciones de onda para los campos

2.3.1. Ecuaciones de Maxwell(Variacion armonica)
Considerando que las fuentes varian armonicamente con el tiempo, las
ecuaciones de Maxwell se pueden simplificar de la siguiente forma:

∇.E~ =ρ

~ =0
∇.H
~ = −jωµH
∇×E ~
~ = J~ + jωE
∇×H ~

2.3.2. Ecuacion de continuidad

De las ecuaciones anteriores se deduce la ecuación de continuidad, toman-
do para ello la divergencia de la Ley de Ampère, y teniendo en cuenta que
la divergencia del rotacional es cero.
~ = J~ + jωE
∇×H ~
~ = ∇.J~ + jω∇.E
∇.∇ × H ~
ρ
0 = ∇.J~ + jω

∇.J~ + jωρ = 0

Las ecuaciones de Maxwell, desde un punto de vista matemático son un

sistema de ecuaciones diferenciales vectoriales de primer orden, apareciendo
entremezclados los campos eléctricos y magnéticos. A continuación se van
a obtener unas nuevas ecuaciones diferenciales, de segundo orden donde se
encuentren separados los campos

7
2.3.3. Ecuación de onda para el campo eléctrico
Tomando el rotacional de la Ley de Faraday se obtiene la ecuación de
onda para el campo eléctrico
~ = −jωµH
∇×E ~

~ = −jωµ∇ × H
∇×∇×E ~
 
∇∇.E ~ = −jωµ J~ + jωE
~ − ∇2 E ~

~ = jωµJ~ + ∇ ρ
 
2~ 2
∇ E + ω µE


2.3.4. Ecuación de onda para el campo magnético

Tomando el rotacional de la Ley de Ampère se obtiene la ecuación de
onda para el campo magnético.
~ = J~ + jωE
∇×H ~

~ = ∇ × J~ + jω∇ × E
∇×∇×H ~
~ − ∇2 H
∇∇.H ~ = ∇ × J~ + µω 2 H
~
~ + ω 2 µH
∇2 H ~ = −∇ × J~

2.3.5. Definición de los potenciales

Para simplificar el cálculo de los campos eléctricos y magnéticos se puede
recurrir a una funciones potenciales escalar y vector, que simplifiquen los
cálculos. El potencial vector se puede definir teniendo en cuenta la ley de
Gauss para el flujo magnético. Si se define B como el rotacional de un vector,
automáticamente se cumple que la divergencia es cero
~ =0
∇.B

~ =∇×A
B ~

El potencial escalar se puede definir a partir de la Ley de Faraday y de

la definición del potencial vector. Teniendo en cuenta que el rotacional del
gradiente de una función es cero, se puede definir el potencial escalar como:
~ = −jωµH
∇×E ~
 
~ + jω A
∇× E ~ =0

~ + jω A
E ~ = −∇Φ

8
2.3.6. Ecuaciones de onda de los potenciales
La ley de Ampère y las definiciones anteriores nos permiten obtener la
siguiente ecuación de onda para el potencial vector eléctrico
 
~ + ω 2 µA
∇2 A ~ = −µJ~ + ∇ ∇.A
~ + jωµΦ

Utilizando la Ley de Gauss para el campo eléctrico se obtiene la ecuación

de onda para el potencia escalar
ρ  ~ + ω 2 µΦ

∇2 Φ + ω 2 µΦ = − + −jω∇.A

Se ha definido el campo magnético a partir del rotacional del potencial
vector, pero es necesario definir también su divergencia. Esta relación se
denomina condición de Lorentz.
~ + jωµΦ = 0
∇.A

Es posible simplificar las expresiones de las ecuaciones de onda para los

potenciales
ρ
∇2 Φ + k 2 Φ = −

∇ A + k A = −µJ~
2~ 2~

Donde K = ω 2 µ.

2.3.7. Soluciones Integrales para los Potenciales

La solución a estas ecuaciones tiene la forma: Potencial escalar.
~0
1 e−jk|~r−r |
Z Z Z
0

Φ= ρ ~
r dv 0
 4π |~r − ~r0 |

Potencial vectorial.
~0
e−jk|~r−r |
Z Z Z
~=µ J ~r0 dv 0


A 0
4π |~r − ~r |

2.3.8. Expresiones generales de los campos

Los campos eléctricos y magnéticos se pueden obtener a partir de los
potenciales mediante las expresiones siguientes

~ = 1∇ × A
H ~
µ

~ = −jω A
E ~ − ∇Φ

Para calcular los campos es necesario calcular los gradientes y rotacionales

de los potenciales.
1
Z Z Z
∇G ~r, ~r0 ρ ~r0 dv 0



∇Φ (~r) =


9
 !
0
e−jkR dG(R) b
∇G ~r, ~r = ∇ = R
4πR dR
1 e−jkR b
 
∇G ~r, ~r0 = −jk −


R
R 4πR
Z Z Z Z Z Z
0
0
0
~ = µ∇ × ∇G ~r, ~r0 × J ~r0 dv 0



∇×A G ~r, ~r J ~r dv = µ ×

Finalmente, simplificando las anteriores expresiones se obtiene el campo

eléctrico, válido en todos los puntos del espacio.

e−jkR 0 e−jkR 0
~ = 1 jk b − ωµ
Z Z Z Z Z Z  
E Rρ
b ~r0 2
dv + ρ ~r R0

dv
4π 4πR 4π k 4πR
El campo magnético es

e−jkR 0
 e−jkR
~ =− 1 jk
Z Z Z Z Z Z 
H R r0
b×J ~
2
dv − r0
b×J ~
R dv 0
4π R 4π R
A distancias muy cercanas a las fuentes predominan los campos que son
proporcionales a 1/r2 , mientras que a grandes distancias predominan los
proporcionales a 1/r.

2.3.9. Campos inducidos

Son proporcionales a 1/r2 y corresponden a las leyes de Coulomb y Biot
y Savart, con un término adicional de fase. Si k=0 se obtienen las expresiones
conocidas de los campos de una distribución de cargas y de corrientes.

e−jkR
~i = 1
Z Z Z
E Rρ
b ~r0 dv 0
4π R2

e−jkR 0
~i = − 1
Z Z Z
H R r0
b × J~ ~ dv
4π R2

2.3.10. Campos radiados

Los campos radiados son proporcionales a 1/r y son los que contribuyen
a la radiación a grandes distancias de las fuentes.

~ r = jk b − ωµ e−jkR 0
Z Z Z  
0

E ρ ~r R dv
4π k R

e−jkR 0
~ r = − jk
Z Z Z
H Rb × J~ ~r0 dv
4π R
Los campos radiados se pueden calcular a partir de las expresiones generales
de los campos, realizando las correspondientes integrales. Si estamos a una
distancia suficientemente grande de la antena, podemos hacer una serie de
aproximaciones que nos simplificarán los cálculos. Se puede considerar que
todas las ondas originadas en la antena siguen trayectorias paralelas hasta
el punto de campo, es decir Al ser las trayectorias de los rayos paralelas, la

10
diferencia de caminos recorridos por las diferentes ondas se puede calcular
como
R = r − rb.~r0
Las ondas producidas en cada punto se pueden aproximar por una onda
centrada en el origen de coordenadas con un desfase adicional equivalente a
la diferencia de caminos.

Los campos magnéticos inducidos y radiados son iguales cuando se cumple

la condición
1 k
2
=
R R
λ
R=
2π

2.3.11. Expresiones aproximadas para los campos radiados

Los campos eléctricos y magnéticos son perpendiculares entre sı́ y sólo
tienen componentes tangenciales. La relación entre sus módulos es la impedan-
cia caracterı́stica del medio, que no debe confundirse con la eficiencia.

~ = − jk rb × A
H ~
µ
   
E ~ − rbA
~ = −jω A ~ = jω rb × r × A
~
 
~ =η H
E ~ × rb

Desarrollando los productos vectoriales, se pueden obtener las componentes

de los campos radiados
Er = 0, Hr = 0
Eφ
Eθ = −jωAθ , Hθ = −
η
Eθ
Eφ = −jωAφ , Hθ =
η

2.3.12. Vector de radiacion

El potencial vector es el producto de dos términos, por una parte tenemos
una onda esférica centrada en el origen de coordenadas, y por otra parte una
integral que vamos a denominar vector de radiación N.
jkr
~ = µe
Z Z Z
J~ ~r0 ejkr dv 0


A
4πr
jkr
~ = µe N
A ~
4πr
Z Z Z
0
~ = J~ ~r0 ejkr dv 0


N

El vector de radiación se puede considerar como la suma vectorial de

todas las corrientes multiplicadas por un término que representa la diferencia
de fase entre la onda que produce y la que tendrı́a si la onda estuviera situada
en el origen de coordenadas.

11
2.4. Antenas
Las antenas son las partes de los sistemas de telecomunicación especı́fi-
camente diseñadas para radiar o recibir ondas electromagnéticas.
También se pueden definir como los dispositivos que adaptan las ondas
guiadas, que se transmiten por conductores o guı́as, a las ondas que se pro-
pagan en el espacio libre.

2.4.1. Diagrama de radiacion

El diagrama de radiación de una antena se define como la representación
gráfica de las caracterı́sticas de radiación en función de la dirección angular.
Sistema de coordenadas
Se utilizará habitualmente un sistema de coordenadas esférico. Las tres vari-
ables de un sistema esférico son r, θ, φ .

Figura 2.4: Sistema de coordenadas

2.4.2. Diagramas bidimencionales

Se puede representar el campo eléctrico, magnético o la densidad de po-
tencia radiada. Dado que los campos son magnitudes vectoriales se pueden
representar el módulo o la fase de sus componentes. Las formas de repre-
sentación pueden ser tridimensionales o bidimensionales, en escalas lineal
o logarı́tmica. La siguiente figura es la representación tridimensional de los
campos radiados por una antena.

Figura 2.5: Diagrama tridimencional

Dada la dificultad de representar gráficamente el diagrama tridimen-

sional se opta por representar cortes del diagrama en coordenadas polares

12
o cartesianas. Los cortes corresponden a la intersección del diagrama 3D
con planos. Un corte bidimensional en coordenadas polares se representarı́a
como

Figura 2.6: Diagrama bidimencional

En coordenadas cartesianas y escala logarı́tmica

Figura 2.7: Diagrama logaritmico

2.5. Parametros fundamentales de las antenas

2.5.1. Densidad de potencia radiada
La densidad de potencia radiada se define como la potencia por unidad
de superficie en una determinada dirección. Las unidades son watios por
metro cuadrado. Se puede calcular a partir de los valores eficaces de los
campos como:  
P~ (θ, φ) = Re E~ ×H~

La relación entre el módulo del campo eléctrico y el módulo del campo

magnético es la impedancia caracterı́stica del medio
 
~
E 
 =η
~
H

Por lo tanto, la densidad de potencia radiada también se puede calcular a

partir de las dos componentes del campo eléctrico.

Eθ2 + Eφ2
P~ (θ, φ) =
η

13
La potencia total radiada se puede obtener como la integral de la densidad
de potencia en una esfera que encierre a la antena.
Z Z
Wr = P~ (θ, φ) d~s

La intensidad de radiación es la potencia radiada por unidad de ángulo

sólido en una determinada dirección. Las unidades son watios por estereora-
dián. Dicho parámetro es independiente de la distancia a la que se encuentre
la antena emisora.
La relación entre la intensidad de radiación y la densidad de potencia radiada
es:
K (θ, φ) = P (θ, φ) r2
La potencia total radiada se puede calcular integrando la intensidad de
radiación en todas las direcciones del espacio.
Z Z Z Z
Wr = K (θ, φ) dΩ = K (θ, φ) sin θdΩdφ

2.5.2. Directividad
La Directividad de una antena se define como la relación entre la densi-
dad de potencia radiada en una dirección, a una distancia, y la densidad de
potencia que radiarı́a a la misma distancia una antena isotrópica,, a igualdad
de potencia total radiada.

P (θ, φ)
D (θ, φ) = Wr
4πr 2

Si no se especifica la dirección angular, se sobreentiende que la Directividad

se refiere a la dirección de máxima radiación.
La directividad se puede obtener en general a partir del diagrama de
radiación de la antena.
Pmax
D (θ, φ) = R R P (θ,φ) sin θdΩdφ
4πr 2

Simplificando términos, resulta

4π 4π
D (θ, φ) = R R P (θ,φ)
=
sin θdΩdφ Ωe
Pmax

Ωe se define como el ángulo sólido equivalente.

Para antenas directivas, con un solo lóbulo principal y lóbulos secun-

darios de nivel despreciable, se puede obtener una directividad aproximada
considerando que se produce radiación uniforme en el ángulo sólido definido
a partir de los anchos del haz a −3dB en los dos planos principales del
diagrama de radiación.
4π 4π
D= =
Ωe θ1 θ2

14
2.5.3. Ganancia
La ganancia de una antena se define como la relación entre la densidad
de potencia radiada en una dirección y la densidad de potencia que radiarı́a
una antena isotrópicaideal, a igualdad de distancias y potencias entregadas
a la antena.
P (θ, φ)
D (θ, φ) = We
4πr 2
Si no se especifica la dirección angular, se sobreentiende que la Ganancia se
refiere a la dirección de máxima radiación.
Pmax
D (θ, φ) = We
4πr2

En la definición de Directividad se habla de potencia radiada por la antena,

mientras que en la definición de ganancia se habla de potencia entregada a
la antena. La diferencia entre ambas potencias es la potencia disipada por
la antena, debida a pérdidas óhmicas. La eficiencia se puede definir como la
relación entre la potencia radiada por una antena y la potencia entregada a
la misma. La eficiencia es un número comprendido entre 0 y 1. La relación
entre la ganancia y la directividad es la eficiencia

G (θ, φ) = D (θ, φ) η

Si una antena no tiene pérdidas óhmicas, la Directividad y la Ganancia son

iguales.

2.5.4. Polarización
La polarización de una antena es la polarización de la onda radiada por
dicha antena en una dirección dada.
La polarización de una onda es la figura geométrica determinada por el ex-
tremo del vector que representa al campo eléctrico en función del tiempo,
en una posición dada. Para ondas con variación sinusoidal dicha figura es en
general una elipse. Hay una serie de casos particulares
Si la figura trazada es una recta, la onda se denomina linealmente polariza-
da, si es un cı́rculo circularmente polarizada. El sentido de giro del campo
eléctrico, para una onda que se aleja del observador, determina si la onda
está polarizada circularmente a derechas o a izquierda. Si el sentido de giro
coincide con las agujas del reloj, la polarización es circular a derechas. Si el
sentido de giro es contrario a las agujas del reloj, la polarización es circular
a izquierdas. El mismo convenio aplica a las ondas con polarización elı́ptica.
Se define la relación axial de una onda polarizada elı́pticamente, como la
relación entre los ejes mayor y menor de la elipse de polarización. La relación
axial toma valores comprendidos entre 1 e infinito.
Los campos se pueden representar en notación fasorial. Para determinar la
variación temporal es suficiente con determinar el valor real de cada una de
las componentes.

15
2.5.5. Impedancia
La impedancia de una antena se define como la relación entre la tensión
y la corriente en sus terminales de entrada. Dicha impedancia es en general
compleja. La parte real se denomina resistencia de antena y la parte imagi-
naria, reactancia de antena.

Vi
Zi = = Ra + jXa
Ii
Se define la resistencia de radiación como la relación entre la potencia total
radiada por una antena y el valor eficaz de la corriente en sus terminales de
entrada, elevada al cuadrado.
Se refine la resistencia óhmica de una antena como la relación entre la poten-
cia disipada por efecto de pérdidas resistivas y la corriente en sus terminales
al cuadrado.
Por lo tanto la resistencia de antena la podemos considerar como la suma
de la resistencia de radiación y la resistencia óhmica.

Vi
Zi = = Ra + jXa = (Rr + RΩ ) + jXa
Ii
La eficiencia de una antena se puede obtener a partir de las resistencias de
radiación y óhmicas, teniendo en cuenta que es la relación entre la potencia
total radiada y la potencia entregada a la antena.

Wt Wt I 2 Rr Rr
η= = = 2 =
We Wt + WΩ I (Rr + RΩ ) (Rr + RΩ )

2.5.6. Area y longitud efectivas

El área efectiva se define como la relación entre la potencia recibida
y la densidad de potencia incidente en una antena. La antena debe estar
adaptada a la carga, de forma que la potencia transferida sea la máxima.
La onda recibida debe estar adaptada en polarización a la antena
Wr
Aef =
Pi
La longitud efectiva de una antena linealmente polarizada se define como
la relación entre la tensión inducida en una antena en circuito abierto y el
campo incidente en la misma.
Va
lef =
Ei

16
Capı́tulo 3

Esquemas de Regulación

Se entiende por RNI a aquellas radiaciones del espectro electromagnético

que no tienen energı́a suficiente para ionizar la materia. En tal caso, la
Emisión es la radiación producida por una única fuente de radiofrecuencia,
mientras que la Inmisión es la radiación resultante del aporte de todas las
fuentes de radiofrecuencias, cuyos campos están presentes en el lugar.

La Exposición Poblacional o No Controlada, corresponde a situa-

ciones en donde el público en general puede estar expuesto a campos eléctri-
cos, magnéticos o electromagnéticos o en donde las personas expuestas, de-
bido a su trabajo, pueden no saber de la potencial exposición y no pueden
ejercer control sobre ella.

Los esquemas de regulación con que se cuentan para hacer las mediciones
de las radiaciones no ionizantes, en este caso OEM, se clasifican en dos
partes:

Método predictivo.

Método de medición.

Algunas compañias proponen emplear ambos métodos, para evitar medi-

ciones y gastos innecesarios en aquellos casos simples y que, con ayuda de
cálculos teóricos, se puede determinar que las RNI están por debajo de nive-
les predefinidos como seguros. Por otra parte, hay otras que aseguran que,
dada la complejidad de las emisiones y el crecimiento del número de fuentes
de radiofrecuencias, el único método posible para garantizar la no superación
de lı́mites predefinidos, es realizar mediciones según un protocolo especı́fico.
Segun nuestro punto de vista, la aplicación de uno u otro método debe
ser objeto de estudios particulares y dependen de:

Los tipos de fuentes de radiofrecuencias.

La densidad de las fuentes intervinientes y su emplazamiento, en el

ámbito de medición.

La potencia y las bandas de frecuencia de las fuentes intervinientes.

La definición polı́tica en el esquema de regulación.

17
 Viendo más en detalle el problema, podemos establecer la siguiente di-
visión conceptual, basado en la definición de Campo Cercano y Campo Le-
jano indicada en la siguiente figura:

Si se considera Campo Lejano −→ evaluación de valores de RNI por

cálculo, basado en Método de Predicción.

Si se superan los lı́mites de Máxima Exposición Permitida Poblacional

(MEP) −→ se emplea Método de Medición.

3.1. Método de Predicción

Planteamos aquı́ un Sitio Mono-Antena, en donde las ecuaciones sola-
mente son válidas para cálculos en el campo lejano, pero pueden utilizarse
para predecir el peor de los casos:
s
(P RA)(1,64)(2,56)(F 2 )
r=
4πS
s
(P IRE)(2,56)(F 2 )
r=
4πS
Donde:
W
S: Densidad de Potencia Máxima (MEP) [ m2 ].

PRA y PIRE: Potencia en antena [W ].

F: Atenuación en veces de la radiación para un ángulo dado de incidencia

en el plano vertical, (F=1, peor caso).

2,56: Factor de reflexión empı́rico (campos reflejados adicionados en fase

con campo incidente directo −60 %).

r: Distancia desde la antena [m].

Si la distancia desde la antena a todo punto accesible por el público en

general es ≥ r, no se requiere verificar el sitio mediante mediciones.

18
3.2. Método de Medición
Para su aplicación se deben realizar los siguientes pasos:

Definición e inspección del lugar escogido para la medición y determi-

nación de las fuentes de radiofrecuencia, tipos de emisión, caracterı́sti-
cas de irradiación y entorno circundante;
Determinación del esquema de medición (campo cercano, campo le-
jano);
Determinación del protocolo de medición a aplicar y el instrumental y
sondas a emplear en el proceso de medición;
Definición de puntos a medir, según el protocolo escogido y los puntos
de mayor riesgo (externos e internos);
Desarrollo de las mediciones y confección de informes.

El procedimiento particular de medición que se sigue, por ejemplo, en la

regulación peruana, es:
Campo Cercano se mide E, H o ambos .
Campo Lejano se mide E o H y se obtiene S
Secuencia: Se mide inmisión.
Inmisión: empleo de instrumentos de banda ancha (detectores de ra-
diación electromagnética no sintonizables), con sondas de medición E
y H isotrópicas.
Emisión: empleo de instrumentos de banda angosta (medidores de in-
tensidad de campo, analizadores de espectros, etc., sintonizables), con
antenas aptas para rangos de frecuencia de medición.
Todos los instrumentos, antenas y sondas deben tener certificado de
calibración (fabricante o laboratorio acreditado en paı́s de origen).
Registro del valor de la medición realizada, más las incertidumbres
especificadas (fabricante), más el error del método empleado.

3.3. Monitoreo Continuo

Los métodos evaluados son de aplicaciones puntuales o repetitivas, pero
con periodos largos de tiempo sin mediciones. Hoy hay nuevos modelos de
mediciones, basados en detecciones continuas (7 x 24) de zonas potencial-
mente riesgosas, con datos publicados en Internet, de acceso libre para los
ciudadanos, que deberı́an ser objeto de análisis y evaluación para su imple-
mentación, por parte de las autoridades con competencia en las RNI.
Aplicar este tipo de monitoreo, trae una serie de ventajas, a saber:
Población: dispone de datos objetivos, las 24 horas del dı́a, para asegu-
rarse que los valores de radiaciones están por debajo de los lı́mites
normativos.

19
Municipios: ofrecen a sus vecinos la tranquilidad que supone disponer de
una red exhaustiva de monitorización de la radiación electromagnética.

Prestadores de Servicios: se reduce la percepción de peligro o alarma

social que generan sus estaciones base.

3.4. Mediciones en CTIC

Para empezar a hacer mediciones en las instalaciones del CTIC primero
debemos de hacer un calculo teórico aproximado de la radiación existente
en un laboratorio, para esto vemos los posibles generadores (fuentes) de
radiacion, los cuales son:
Antenas del internet inalámbrico(wireless).

Monitores LCD.

CPU (los circuitos que lo componen).

Aire acondicionado.
De los dispositivos electrónicos mencionados, cada uno de estos inter-
viene en mayor o menor medida a la radiacion no ionizante existente en un
laboratorio, pasaremos a analizar cada uno de estos dispositivos:

3.4.1. Antenas de internet inalámbrico(wireless)

El uso de la tecnologı́a basada en enlaces inalámbricos entre las redes

de computadoras de los laboratorios de CTIC se encuentran en la banda
de frecuencia de 2,45 GHz (2.450.000 Hertz o ciclos por segundo) operando
bajo el estándar IEEE 802.11b.(véase la siguiente figura)

20
 Figura 3.1: Standart IEEE

Y acontinuación se muestra como está distribuido las frecuencias de las

ondas ;que comunmente se utilizan para comunicación por internet, en el
espectro electromagnético.

3.4.2. Monitores LCD

Los monitores de pantalla plana (también llamados FPD que significa
pantallas de panel plano) se popularizan cada vez más, ya que ocupan menos
espacio y son menos pesados que las tradicionales pantallas CRT.(tubo de
rayos catodicos)

La LCD (pantalla de cristal lı́quido) se basa en una pantalla hecha de

dos placas paralelas transparentes ranuradas y orientadas a 90o una de otra.
El espacio entre ellas alberga una fina capa de lı́quido que contiene ciertas
moléculas (cristales lı́quidos) que poseen la propiedad de orientarse cuando

21
se ven expuestas a la corriente eléctrica.Combinada con una fuente de luz,
la primera placa actúa como un filtro de polarización, permitiendo el paso
sólo de aquellos componentes de luz cuya oscilación es paralela a las ranuras.

Figura 3.2: Monitor LCD

Durante la ausencia de corriente eléctrica, la segunda placa bloquea la

luz, actuando como un filtro de polarización perpendicular.

Figura 3.3: Monitor LCD

Cuando se encuentra encendida, los cristales se alinean de manera pro-

gresiva en la dirección del campo eléctrico y de esta manera pueden cruzar
la segunda placa. Al controlar localmente la orientación de los cristales, es
posible crear pı́xeles.

Además, la tecnologı́a utilizada por los monitores de pantalla plana suele

utilizar menos energı́a (inferior a 10 W, a diferencia de los 100 W de las
pantallas CRT(Tubo de rayos catodicos) y emite menos radiación electro-
magnética que generan campos de frecuencia intermedia que van desde 300
Hz a 100 kHz, y que son más bajas que las radiofrecuencias y más altas
que las frecuencias extremadamente bajas. La intensidad del campo electro-
magnético disminuye rápidamente el alejarse de la fuente de emisión. En las
partes laterales y traseras de algunos monitores la intensidad de campo es
menor por el blindaje metalico que lo protege.

Las radiaciones que inciden en un objeto pueden atravesarlo, ser ab-

sorbidas o reflejadas por él. La capacidad de penetración en un objeto de

22
una radiación dada depende de su longitud de onda y de las carácterı́sticas
del material.
La radiación ionizante producida por los monitores representa una fracción
pequeña de la radiación de fondo.
Asi, estos monitores solamente emiten una pequeña porcion del espectro
electromagnetico, los posibles efectos en la salud debido a campos electricos
o magneticos solamente se han demostrado en el trabajo con CRT sin em-
bargo. Los efectos biológicos de esta radiación de bajo nivel siguen siendo
tema de controversia y su análisis está fuera del objetivo de este trabajo,
por tanto es conveniente apagar el monitor cuando no se use, de manera de
minimizar el tiempo de exposición.

3.4.3. CPU (los circuitos que lo componen)

La inmensa mayorı́a de estos, funcionan mediante mecanismos electrónicos
que emiten radiaciones electromagnéticas de todo tipo de frecuencia y po-
tencia. Una PC emite radiación electromagnetica VLF (frecuencias interme-
dias)y radiación electromagnetica ELF(frecuencias extremadamente bajas)
muchos de ellos generan niveles de radiación que están por debajo de los
niveles máximos permitidos debido al blindaje del case metalico que lo
cubre aun asi estas radiaciones suelen propagarse en los orificios de venti-
lación y mal garantizados paneles.

Figura 3.4: Campo electrico y magnetico de una PC

3.4.4. Aire acondicionado

Existe una gran cantidad de elementos y situaciones que vician el aire y
que, sobre todo, reducen de forma alarmante la concentración de iones neg-
ativos, altamente favorables para la salud. El empleo del aire acondicionado
emiten partı́culas cargadas positivamente (+) o ciertas dosis de radiactivi-
dad ligeramente ionizante, especialmente en ambientes cerrados, se ha de-
mostrado que todos estos dispositivos eléctricos y situaciones varias generan
grandes concentraciones de iones positivos (+), altamente perjudiciales para
la salud. su enorme carga eléctrica positiva o el predominio de iones posi-
tivos (+) con respecto a los iones negativos (-); pero este desequilibrio iónico

23
también lo hallamos permanentemente en el interior de industrias, viviendas
y oficinas.

3.5. Cálculo aproximado de radiación

Usando las ecuaciones válidas para calcular la potencia de radiación de
campos lejanos producidos por una antena se obtienen los resultados que se
presentará en la siguiente tabla, estos son calculados en una laboratorio de
CTIC.
se uso la siguiente relación:
s
(P IRE)(2,56)(F 2 )
r=
4πS

y se obtuvo los siguientes resultados de los cálculos de predicción:

Apartir de estos datos se pudo obtener la siguiente figura:

S (mW/cm2 ) Vs R(m)

24
 Y de la tabla anterior se puede inferir un gráfico de densidad de radiación
en un ambiente aproximado a un laboratorio de CTIC de 5x15m.

Acontinuación se muestra la tabla de valores limites permitidos para el

Perú dado por el D.S. 038-2003-MTC que establece los “Lı́mites Máximos
Permisibles de Radiaciones No Ionizantes en Telecomunicaciones”.

25
Parte II

26
Capı́tulo 4

Software para simulación de

Antenas

Para el trabajo a realizar necesitamos software que realice simulaciones

de antenas, estos deben ser capaces de calcluar las potencias y graficar la
radiación electromagnética que emite dicha antena. En la actualidad exis-
ten muchos programas con las virtudes antes mencionadas, la mayorı́a de
estos estan basados o son versiones mejoradas del paquete NEC(Numerical
Electromagnetics Code) que es un programa para modelación de antenas, es
de código libre y esta hecho en un inicio en el lenguaje Fortran y luego fue
trasncrito al lenguaje C.

MMANA-GAL
MMANA es un analizador de antenas basado en el método de momentos.
MMANA-GAL funciona bajo Windows, es software libre, no consume mu-
chos recursos de máquina. Es especializado en antenas hechas de conductores
rectilı́neos sin necesidad de plano de masa.
Algunas de las limitaciones de este programa, es el hecho que sólo puede
modelizar antenas compuestas de conductores rectilı́neos, y desnudos. Y
aunque a mejorado en el tratamiento de la tierra o masa, no funciona bien
cuando un conductor está a pocos centı́metros del suelo.
La interfaz del programa es bastante simple, lo que si puede resultar un
tanto complejo es la introducción de la estructura de la antena dependiendo
de la misma. Permite obtener los parámetros básicos de una antena, como
los diagramas de radiación, ganancia, adaptación de impedancia, etc.

27
GAL-ANA
GAL-ANA no es software libre. Podrı́a decirse que es una versión mejo-
rada de MMANA-GAL, puesto que es realizado por los radioaficionados ale-
manes que mejoraron el MMANA, por cuanto posee caracterı́sticas bastante
parecidas a este.

4NEC2
4NEC2 es de libre distribución, usa Nec2, Nec4 y es una herramien-
ta basada en ventanas para crear, ver, optimizar y chequear estilos de ge-
ometrı́as de antenas en 2D y 3D y generar, presentar y/o comparar diagra-
mas de radiación de campo cercano y lejano. Permite al igual que la mayorı́a
de programas para esto, hallar diagramas de radiación, acople, eficiencia, etc.
Presenta ventanas independientes para la geometrı́a, vista en 3D, plots de
eficiencia, etc. Presenta también una interfaz especial Built para facilitar la
construcción de los modelos de antenas.

28
EZNEC
EZNEC, no es libre, se encuentra en su versión 5.0, se presentan tres
distribuciones, uná básica EZNEC, EZNEC+ y EZNEC pro. Parecido en
su interfaz y proceso a seguir para simular a MMANA-GAL. El EZNEC+
está diseñado para una experiencia avanzada. La interfaz es identica a EZNEC,
pero presenta mejores caracterı́sticas, como: permite analizar antenas com-
plejas (1500 segmentos versus 500 del programa estándar), análisis de campo
lejano para polarización circular, double-precision, etc.

SuperNEC
SuperNEC utiliza el método de momentos electromagnéticos (EM) para
la simulación, puede correr bajo plataformas Windows o Linux. Presenta una
entrada GUI 3D fácil de usar basada en MatLab, hace uso de ensamblaje
multi-nivel, provee una herramienta para crear la estructura inicial bastante
sencilla. La salida provee toda la información necesaria para analizar an-
tenas incluido diagramas de radiación en 3D&2D, diagrama de la carta de
smith, acoplación, eficiencia, etc. SuperNEC incluye optimización de código
para procesadores Intel, varios métodos de resolución rápida y ejecución en
paralelo para optimizar el proceso en caso de problemas largos. SuperNEC
puede ser utilizado para determinar el lugar óptimo para las antenas en
estructuras como por ejemplo en aviones, etc.
No es software libre, y el costo depende del producto que se desee y
de la utilización que se de al mismo, sea esta comercial o académica. El
superNEC básico en versión comercial cuesta $6080 y el académico $1330.
Se puede conseguir una versión Lite para utilización académica libre.

29
30
4.1. Métodos computacionales
Los principales métodos numéricos usadas en electormagnetismo son:

1. El método de los momentos(MoM)

2. El método de las diferencias finitas en el dominio del tiempo(FDTD)

3. El método de los elementos finitos(FEM)

Problemas electromagneticos que tengan geometrias simples pueden ser de-

scritas usando sistemas de coordenadas básicas(rectangular, cilindrico, esféri-
cas)pueden algunas veces pueden ser solucionadas analı́ticamente.Dependiendo
de la complejidad del problema, la solución puede ser obtenida por inte-
gración directa o usando técnicas como la separacion de variables. Para la
solución numérica mostramos los siguientes métodos:

4.1.1. Método de los momentos (MoM)

Este método de cálculo es uno de los más usados en la actualidad para
determinar los campos emitidos o recibidos por estructuras radiantes. En
el primer caso se supone que una fuente externa establece corrientes vari-
ables en el tiempo en los conductores de la antena, y se calcula los cam-
pos generados por ellas(antena emisora), mientras que el segundo supone
que una onda (habitualmente plana) incide sobre una estructura inducien-
do corrientes variables en el tiempo sobre ella. Este ultimo caso se puede
analizarse desde el punto de vista de una fuerza electromotriz (fem) util in-
ducida sobre los conductores(antena receptora) o con el objetivo de calcular
los campos reirradiados por la estructura para analizar sus propiedades de
dispersion(scattering) de ondas electromagneticas.
El corazón del metodo recide en la linealidad de las ecuaciones de Maxwell,
que desde el punto de vista matemático se puede representar mediante op-
eradores lineales aplicados a funciones de las coordenadas espaciales y el
tiempo.Por ejemplo las ecuaciones de onda para el campo magnetico y elec-
trico toman la forma:
!
2 1 ∂2f 2 1 ∂2
∇ f− 2 2 =g ⇒ ∇ − 2 2 f = g ⇒ Lf = g
c ∂t c ∂t
en esta ecuación g es la exitación o fuente(funcion conocida) y f es el campo
o respuesta(incognita a ser determinada) Este es un caso particular de la
ecuacion lineal: L(f ) = g
donde L es un operador lineal.
Existen dos problemas básicos:

El problema del análisis involucra hallar f (r, t) dados L y g(r, t), en

este caso se conoce el sistema y la exitación, y se desea hallar la re-
spuesta.

El problema de sintesis involucra determinar L cuando se dan f (r, t)

y g(r, t), en este caso dada la exitación y la respuesta deseada, se pide
hallar el sistema que la produce.

31
Las funciones f (r, t) satisfacen determinadas condiciones.Por una parte, es
necesario poder definir un producto interno:
< f 1, f 2 >=< f 2, f 1 >
< αf 1 + βf 2, f 3 >= α < f 1, f 3 > +β < f 2, f 3 >
< f, f ∗ >> 0 si f 6= 0
< f, f ∗ >= 0 si f = 0, donde α y β son escalares. Si la solución de la
ecuacion operacional L(f ) = g existe, entonces existe el operador inverso
L−1 definido por: f = L−1 (g).
El método de los momentos es un procedimiento general de resolucion
de ecuaciones operacionales lineales. Las funciones solución se representan
como un desarrollo en serie de funciones base f1 , f2 , ...
X
f= αn fn
n

donde los factores αn son coeficientes constantes. Generalmente las series

son infinitas, y las funciones base forman un conjunto completo, aunque en
los casos prácticos se obtiene series finitas que dan una solución aproximada
al problema.
Definido el conjunto de funciones base, se tiene:
X
αn L(fn ) = g
n

Definimos ahora un conjunto de funciones peso o funciones de testeoW1 , W2 , ...

con lo que nos queda:
X
< Wm , L(fn ) >=< Wn , g > m = 1, 2, ...
n

Esto representa un conjunto de n ecuaciones lineales en los coeficientes αn

que puede escribirse en forma matricial como :

[Lmn ][αn ] = [gm ]

Si la matriz [Lmn ] es no singular, existe su inversa y podemos calcular los

coeficientes αn :
[αn ] = [Lmn ]−1 [gm ]
Esta solución se conoce como el método de Galerkin cuando se toma wn = fn
La eleccion del conjunto de funciones base es la clave para el adecuado
comportamineto del método. Existen al menos dos formas de elegir:
Ajuste por puntos. La integración involucrada en el calculo de los el-
emtos de la matriz [Lmn ] puede ser de dificil evaluacion en terminos
practicos. Una forma de obtener soluciones aproximadas es requerir
que la condición L(f ) = g se cumpla solamente en un conjunto discre-
tos de puntos de la region de interés.

Bases subseccionales. En muchos casos se define funciones base cada

una de las cuales existe solo en subsecciones de la region de interés.
Entonces cada coeficiente αn afecta la aproximacion a la solucion f
sólo sobre una subsección, lo que lleva a matrices más despobladas y
de más sencilla inversión numérica.

32
Existen actualmente muchas implementaciones del método de momentos,
algunas de ellas comerciales bajo diversas plataformas y otras de uso públi-
co. De estas últimas son relativamente accesibles y están razonablemente
documentadas tres versiones:NEC-2, XLNEC y WIPL-DEMO.

4.1.2. Método de diferencias finitas en el dominio del tiempo

(FDTD)
De entre todas las técnicas numéricas que resuelven las ecuaciones de
Maxwell con ayuda de un ordenador, el método de diferencias fnitas (FD,
del inglés Finite Difference) es una de las más antiguas y posiblemente la
que menor contenido analı́tico tiene. En esencia, la aplicación de este método
consiste en calcular la solución del problema diferencial original en un con-
junto discreto de puntos, que suelen asociarse con los nudos de una red.En
el método FD, este proceso de discretización se realiza mediante la aproxi-
mación de las derivadas parciales que aparecen en las ecuaciones de Maxwell
-en cualquiera de sus versiones temporal o frecuencial- por operadores en
diferencias fnitas. Como resultado se obtiene un sistema de ecuaciones lin-
eales que relacionan los campos en los distintos nudos de la red. El método
FDTD, introducido por K. S. Yee en el año 1966, consiste en la sustitu-
ción de las derivadas parciales que aparecen en las ecuaciones de Maxwell
del rotacional por cocientes de diferencias fnitas centradas. Para poder ll-
evar a cabo esta sustitución de diferenciales por incrementos, se discretiza
el espacio-tiempo mediante una malla discreta. En consecuencia, el prob-
lema diferencial se transforma en un problema en diferencias (algebraico),
es decir, en un sistema de ecuaciones lineales. El esquema de Yee para la
discretización de las ecuaciones de Maxwell es el más famoso y de uso más
frecuente dentro de la gama de métodos FDTD.El cálculo del campo eléctri-
co (magnético) en un nudo de la malla y en un instante de tiempo, sólo
precisa de los valores del campo magnético (eléctrico) en los nudos adya-
centes y en el instante de tiempo anterior. De esta forma, partiendo de unas
condiciones iniciales, el método permite calcular la evolución temporal del
campo electromagnético en la región de interés.
Comenzemos con las ecuaciones de Maxwell en un recinto sin fuentes de
campo. Las ecuaciones de la divergencia:∇ • E = 0 ; ∇ • H = 0 de-
scriben la transversabilidad de la propagacion, mientras que las ecuaciones
del rotor:∇ × E + ∂H ∂E
∂t = 0 ; ∇ × H − σE −  ∂t = 0 establecen las relaciones
de autosostenimiento de los campos en la propagación.

Discretizacion de las ecuaciones de Maxwell según el esquema de Yee.

Las ecuaciones del rotor se pueden escribir por componentes en coordenadas

cartesianas:
     
∂Hx 1 ∂Ey ∂Hy ∂E
∂t = − ∂E z
= 1 ∂Ez
− ∂E x ∂Hz
= 1 ∂Ex
− ∂xy
 µ ∂z ∂y  ∂t  µ ∂x ∂z  ∂t  µ ∂y 
∂Ex 1 ∂Hz ∂Hy ∂Ey 1 ∂Hx ∂Hz ∂Ez 1 ∂Hy ∂Hx
∂t =  ∂y − ∂z − σE x ∂t =  ∂z − ∂x − σE y ∂t =  ∂x − ∂y − σE z

Definamos una grilla tridimencional cubica tal que el punto (i, j, k) está definido
por:xi = ihyj = jhzk = kh y sea f n (i, j, k) el valor de cualquiera de las
componentes del campo en ese punto y en el instante n∆t. Las derivadas

33
espaciales y temporales se pueden escribir por diferencias en la forma:
∂f n (i, j, k) f n (i + 1/2, j, k) − f n (i − 1/2, j, k)
=
∂x h
∂f n (i, j, k) f n+1/2 (i, j, k) − f n−1/2 (i, j, k)
=
∂t δt
Yee define una celda de calculo de modo que las componentes de los campos
E y H son calculadas en posiciones interlaceadas y en pasos alternados de
tiempo.En tres dimenciones se tiene la celda de Yee de la figura se observa
que todas las componentes están presentes en cada cuarto de celda. Las seis
ecuaciones diferenciales en las componentes, que surgen de las ecuaciones de
rotor, llevan a las ecuaciones de diferencias de la forma:
 
∂Hx 1 ∂Ey ∂Ez n+1/2 n−1/2
∂t = µ ∂z − ∂y ⇒ Hx (i, j + 1/2, k + 1/2) = Hx (i, j + 1/2, k + 1/2)+
∆t n
µh [Ey (i, j + 1/2, k + 1) − Ey (i, j + 1/2, k) + Ez (i, j, k + 1/2) − Ezn (i, j + 1, k + 1/2)]
n n

,de forma análoga para las demás componentes y para las componentes del
campo eléctrico. E.

Figura 4.1: Celda de Yee 3D

El método FDTD se aplica habitualmente a problemas de propagación

y dispersión (scattering) de ondas. Los dos problemas básicos son:

determinar las condiciones de frontera de los campos en los bordes del

recinto de integración,

determinar la inclusión de las fuentes que dan origen al campo dentro

del recinto de integración.

Dado que en la práctica no es posible trabajar con recintos infinitos, aún

cuando el problema fı́sico requiera una solución en un recinto ilimitado,
se requiere introducir algún tipo de frontera artificial. Las fronteras deben
colocarse lejos de las estructuras cuyo comportamiento se desea analizar.
Las condiciones de frontera usadas para simular un recinto indefinido son

34
las condiciones de no reflexión o absorción. Estas condiciones simulan ondas
que se propagan más allá de la frontera sin reflexión.
Hay numerosas implementaciones de FDTD en uso actualmente, ya que se
trata del método más popular para la simulación de ondas que se propagan
en guı́as, microstrips, y para cálculo de efectos de compatibilidad electro-
magnética (EMC), especialmente en circuitos impresos (PCB). Sin embargo,
la mayorı́a son programas comerciales de alto costo. Hay una implementación
gratuita, llamada WINFDTD

4.1.3. El método de los elementos finitos

Este método constituye un método numérico destinado a resolver me-
diante ecuaciones matriciales las ecuaciones diferenciales que se plantean
en sistemas discretos (estructuras) o continuos (campos).En el caso de sis-
temas contı́nuos, el método consiste en discretizar el dominio de interés en
Elementos Finitos y resolver, mediante una función de prueba o de aproxi-
mación, la ecuación que rige el sistema en cada EF para luego sumar todas
las soluciones.
Procedimiento de solucion usando el método de elementos finitos.
El método de los elementos finitos es un método de aproximación de
problemas continuos y la solución de cualquier problema utilizando este
método contempla los siguientes pasos:

El continuo se divide en un numero finito de partes, “elementos”, cuyo

comportamiento se especifica mediante un numero finito de parámet-
ros asociados a ciertos puntos caracterı́ticos llamados “nodos”. Estos
nodos son los puntos de union de cada elemento con sus adyacentes.

Aplicar las condiciones de frontera (restricciones) y las cargas externas.

Las incognitas del problema dejan de ser funciones matématicas y

pasan a ser el valor de estas funciones en los nodos.

El comportamiento en el interior de cada elemento queda definido a

partir del comportamiento de los nodos mediante las adecuadas fun-
ciones de interpolacion o funciones de forma.

El método de los elementos finitos , por tanto, se basa en transformar

un cuerpo de naturaleza continua en un método discreto aproximado, es-
ta trasformación se denomina discretización del modelo. El conocimiento de
lo que sucede en le interior de este modelo del cuerpo aproximado, se obtiene
mediante la interpolacı́ón de los valores conocidos en los nodos. Es por tanto
una aproximación de los valores de una función a partir del conocimiento de
un numero determinado y finito de puntos.
aunque puede realizarse una implementacion del método adecuada a necesi-
dades propias, ya existen comercialmente paquetes que implementan el méto-
do y que permiten acceder rapidamente a la solución de un análisis especı́fi-
co. Entre los numerosos paquetes comerciales disponibles destacan: ANSYS,
COMSOL, ALGOR, NASTRAN, ABAQUS, DYNA-3D.

35
4.2. Medición de Campos Magneticos en CTIC
Para tener una idea de la existencia de los campos magnéticos en las
instalaciones de CTIC, nosotros usamos el medidor de campos magneticos
(a falta de un medidor de campos electromagneticos) modelo Magnetic
Field Logger ML-1 de Enviro Mentor

Figura 4.2: Magnetic Filed Logger ML-1

Hicimos las medidas tomando en cuenta las dimensiones de cada am-

biente sean estos los laboratorios, oficinas administrativas, y el auditorio.
Como mostramos en las siguientes figuras la toma de datos se hizo a 3 al-
turas las cuales consideramos son representativas para mostrar una gráfica
de radiación. Las medidas fueron hechas a 30 , 90 y 180 cm .

Puntos a medir primer piso

Figura 4.3: Primer Piso A

36
 Figura 4.4: Primer Piso B

Puntos a medir segundo piso

Figura 4.5: Segundo Piso A

37
 Figura 4.6: Segundo Piso B

Medidas para el data center

Equipo de medida: Campo Magnetico, 30Hz-2KHz
Magnetic Field Logger ML-1 Medidas en µT
Lugar: Data Center Fecha: 05-03-09
Medidas en el punto: Altura sobre el piso
X(m) Y(m) 30cm 90cm 180cm
0.8 1.5 0.07 0.04 0.03
0.8 3 0.17 0.15 0.16
0.8 4.5 0.02 0.02 0.02
1.6 1.5 0.6 0.06 0.03
1.6 3 0.42 0.27 0.07
1.6 4.5 0.09 0.08 0.05
2.4 1.5 0.02 0.06 0.03
2.4 3 0.14 0.09 0.03
2.4 4.5 0.11 0.08 0.19
3.2 1.5 0.01 0.01 0.02
3.2 3 0.04 0.03 0.01
3.2 4.5 0.06 0.1 0.52
*Punto de origen extremo superior derecha

38
Gráficos de Radiación para el Primer Piso

Figura 4.7: Primer Piso A-30 cm

39
Figura 4.8: Primer Piso A-90 cm

40
Figura 4.9: Primer Piso A-180 cm

41
Figura 4.10: Primer Piso B-30 cm

42
Figura 4.11: Primer Piso B-90 cm

43
Figura 4.12: Primer Piso B-180 cm

44
Gráficos de Radiación para el Segundo Piso

Figura 4.13: Segundo Piso A-30 cm

45
Figura 4.14: Segundo Piso A-90 cm

46
Figura 4.15: Segundo Piso A-180 cm

47
Figura 4.16: Segundo Piso B-30 cm

48
Figura 4.17: Segundo Piso B-90 cm

49
Figura 4.18: Segundo Piso B-180 cm

50
Capı́tulo 5

conclusiones

Los equipos que contribuyen a incrementar la cantidad de radiación no

ionizante en un punto dado son:

Monitores LCD: Utilizan cristales liquidos que practicamente no producen

radiación.

CPU (los circuitos que lo componen): Está recubierto por una carcaza
metálica (case) que hace las veces de una Jaula de Faraday y esto impi-
de la salida y entrada masiva de ondas electromagnéticas. Salvo por
los orificios por los cuales podrian irradiar.

Aire acondicionado: Solo produce radiación no ionizante en niveles bajos.

Asi mismo se hizo la predicción teórica de cuales serı́an los niveles de

radiación en el caso de una antena utilizada para proporcionar intenet in-
alámbrico, y comparando estos datos con los niveles permitidos en el Perú se
deduce que la densidad de potencia a la que podrı́a estar sometida una per-
sona en las cercanı́as de una de estas antenas está por debajo de la densidad
de potencia establecida por la norma peruana.

Luego de realizar los gráficos con las medidas tomadas por el medidor
Magnetic Field Logger ML-1 para campos magneticos variables, podemos
afirmar que los campos magnéticos en las instalaciones del ctic son menores
a 0.6 µ T , teniendo mayor presencia en el data center. En los laboratorios
apenas llegan a 0.2 µ T al igual que en Telecom o el servicio técnico y en
las demás instalaciones las medidas indican valores próximos a cero; valores
que estan muy por debajo de los lı́mites de exposición permitidas.
Además debemos mencionar que el medidor de campo magnetico que
usamos trabaja en frecuencias de 30Hz-2KHz que estan por debajo de las
frecuencias producidas por la red inalambrica(2.4GHz), para futuras medi-
ciones se recomienda utilizar un equipo que pueda registrar en estas frecuen-
cias (presentamos en el anexo algunas proformas de equipos a utilizar).

51
Apéndice A

Medidores de campos
electromagneticos

Aquı́ varios instrumentos de medida para realizar mediciones de la in-

tensidad de los campos electromagnéticos de frecuencia baja.

A.0.1. Medidor de campo eléctrico PCE-EM 29

El medidor de campo eléctrico dispone de una sonda triaxial esférica
para la detección de la radiación electromagnética. El medidor de campo
eléctrico es igual de idóneo para medir en transformadores, como también
para una evaluación de campos magnéticos que son producidos por pantallas
de ordenador, televisores, instalaciones industriales (separador magnético,
electromotores, soldadores bajo gas protector, etc.).
Además, también está preparado para detectar radiaciones en el ámbito de
sistemas inalámbricos (Wireless LAN), GSM o para determinar la radiación
de microondas.

Figura A.1: Precio:199 euros(sin IGV) [5]

A.0.2. MEDIDOR DE RADIOFRECUENCIA

Para frecuencias de 800 MHz - 2.5 GHz. El medidor de radiofrecuencia
permite tomar medidas en ambas lecturas, analogico y digital. Dos rangos 1−

52
 Figura A.2: Especificaciones tecnicas

1999µW/m2 0,1 − 199,9nW/cm2 y 0,1−199,9µW/m2 0,01 − 19,99nW/cm2 .

Adicionalmente incluye tonos de audio proporcionales a la intensidad de

campo radiadado o proporcional a la frecuencia, para elanalisis de signos de
pulso.
Este es un medidor diseñado especialmente para la medida del campo elec-
tromagnético o su monitorización, por ejemplo: estaciones de telefonı́a móvil,
equipamiento de hospitales, radar, hornos microondas, actividades laborales
bajo radiación, antenas de TV, estaciones de radio, equipos de soldadura,
equipos de cocción, televisión, ordenadores, fábricas, laboratorios y otros
ambientes. Modelo HF35C.

Figura A.3: Precio 499 dolares [5]

A.0.3. ANALIZADOR DE INTENSIDAD DE CAMPOS RF

Esta poderosa herramienta lo tiene todo: contador de frecuencia, exposi-
ción del espectro, e intensidad RF. Escanee el rango completo de frecuencias
¡desde los 100 kHz hasta arriba de los 2060 MHz!
¿Quiere saber si los teléfonos celulares están operando en su área? Observe
el espectro de frecuencia, luego sintonı́cese la banda de los 800 MHz o 1900

53
MHz para rastrear la fuerza de la señal. ¿Quiere saber cuánta radiación
recibe de una torre de radiodifusión? ¿Qué tan bien funciona su blindaje
para las diferentes frecuencias? ¿Qué tipo de interferencia RF le llega a ust-
ed de su computadora? Si alguna vez quiso analizar la mezcla de señales RF
en su ambiente, ésta es la unidad para usted. Con solo un tamaño de 9”x
4”, esta estupenda unidad a baterı́a está contenida en sólo 1.4 lbs.
Tiene un parlante incorporado para el monitoreo de las señales por audio.
Simple selección de las funciones sobre el menú en pantalla. Pantalla con
retro iluminación y capaz de proyectar el espectro, gráfica de barras y ex-
hibición digital. ¡Incluye antena desmontable para señales de radiodifusión,
(utilice un adaptador optativo F a BNC para medir señales en cable coaxial)
e incluso tiene salida RS-232 con software para PC y cable, de manera que
usted puede cargar y manipular datos en su computadora! Garantı́a de un
año del fabricante.

Figura A.4: Precio:1895 dolares [5]

Figura A.5: Especificaciones técnicas [5]

54
A.0.4. MEDIDOR DE CAMPO TAOMA
Medidor de campos electromagnéticos de 5Hz a 18GHz según la sonda
utilizada, FFT de 5Hz a 100kHz, integra un GPS y datos T◦ C/RH, software
interno para la relaización de informes.La avanzada tecnologı́a del sistema
de medida de campo electromagnético TAOMA, es el “estado del arte” de
la medición ambiental, resultado de aplicar las últimas tecnologı́as que, re-
duciendo el tiempo de trabajo, permiten obtener los más fiables resultados
a las regulaciones relativas a la exposición de la población y los trabajadores
a los campos electromagnéticos.
El sistema TAOMA posibilita medidas EM en el rango 5Hz - 40GHz y
dispone de toda una serie de accesorios opcionales pensados y diseñados
para facilitar la labor del usuario en distintas condiciones. Está construido
con materiales de la más alta calidad para conferirle una estructura robusta
y ligera, y su diseño ergonómico dispone además de protección ambiental
IP-54.Precio dado para sonda E 100kHz a 6GHz, 0.2 a 400V/m

Figura A.6: Precio:8060euros [5]

A.0.5. Analizador de Campos electromagnéticos EMR-300

de Narda
Sistema versátil para las medidas de campos magnéticos

Unidad central del instrumento con una gran cantidad de accesorios

Sondas intercambiables para multitud de aplicaciones y rangos de fre-

cuencia

Medidas no direccionales (isotrópicas) con sonda de tres canales

Gran margen dina mico debido al procesamiento digital de los resul-

tados de los tres canales

Interfaz óptico para la calibración y transferencia de los resultados de

medida

Excelente precisión, con ajuste de cero automático, incluso durante la

exposición al campo

55
 Sencillo manejo

Hermético y resistente a los golpes

Calibrado

Especificaciones del Detector de Radiación EMR-300

56
A.0.6. Analizador de Espectros: Agilent HP 8594E
Analizador de Espectros, marca Hewlett Packard, modelo HP-8594E, el
cual permite realizar mediciones en distintas unidades y modalidades, tiene
las siguientes caracterı́sticas técnicas principales:

Su costo es desde $4500.00 en NATIONAL TEST EQUIPMENT[6]

$6995.00 dolares americanos, y el proveedor es TEST EQUIPMENT[7]

57
PROPUESTA ECONÓMICA

Presupuesto 300-02-Extech Instruments Ciudad: LIMA

Señor ANGELO ROMERO Fecha: 12/03/2009

Encargado: YBE
Atención Angelo Romero Telefono: 2743414
Referencia Solicitud de Cotización E-mail Fax: 2743394
e-mail: angelosb2000@gmail.com

Estimados Señores:
Atendiendo a su solicitud, nos complace presentar a ustedes nuestra oferta por:

It CÓDIGO CANT. DESCRIPCIÓN PRECIO UNIT. TOTAL S/.

1 480836 1 RF EMF Strength Meter 1.020.00 1.020.00

http://www.extech.com/instruments/resources/manuals/480836_UMsp.pdf

CONDICIONES
Precio: En Nuevos Soles, incluye IGV
Forma de pago: 10% de adelanto, 90% contraentrega
Validez oferta: 7 días
Plazo entrega: 4 semanas a partir de la recepión de la orden de compra, consultar para equipos con certificado NIST
Lugar entrega: Nuestras Oficinas
Garantía: Un año contra defectos de fabricación
RUC: 20296021408
Nro Cta. Cte. Soles ARMOTEC S.A.C. Banco Continental 0011-0148-01000-27939-46

Atentamente,
Yefer Bazán E.
ARMOTEC sac
 PROVIENTO S.A.C.
RUC: 20521134586

Proforma #: 20-0309 Fecha: 7-Mar-09

Estimado Royer,
le podemos ofrecer un equipo todavia más económico para la medición de las ondas de baja y alta frecuencia.

Datos del cliente:

Centro de Tecnologias de Informacion y Comunicaciones Universidad Nacional de
Nombre/Razon Social:
Ingenieria (CTIC-UNI)
Dirección: Direccion Av. Túpac Amaru 210 - Rimac / Lima 25 - Perú
Responsable / Attn.: Royer Ticse Torres

Proforma por:
Equipos Medición: x Equipos Generación: Servicios internos: Servicios externos:

Cantidad: Valor (USD)

Descripción Item:
# En Stock Unidad Total
Maletín MK20 con medidores de contaminación por ondas
1 Pcs electro-magneticas de baja y alta frecuencia tipo GIGAHERTZ 620.00 620.00
SOLUTIONS
Este juego de equipos mide la intensidad de ondas
electromagnetcias causados por electrodomesticos
(baja frecuencia) con el equipo ME3830B y las
frecuencias altas (antenas celulares, WLAN etc).
con el equipo HF35C. El display de 3,5 digitos está
respaldado por una indicación acústica que es
proporcional a la intensidad de las ondas y sirve
detectar tambien ondas pulsadas. La antena es
logarítmica periodica y tiene caracteristica de
dirección.

Equipo de frecuencia baja (GIGAHERTZ

SOLUTIONS ME3830B):
Rango de frecuencia: 5Hz - 100 kHz,
Sensibilidad: 1 - 2000 nT (componente magnético)
Sensibilidad: 1 - 2000 V/m (componente eléctrico)
Precisión: +/- 2% (+/- 20 dígitos)

Equipo de frecuencia alta (GIGAHERTZ

SOLUTIONS HF35C):
Rango de frecuencia: 800MHz - 2500 MHz (GSM,
DECT, WLAN, BLUETOOTH, CDMA,...),
Sensibilidad:0,01 - 19990 µW /m2
Precisión: +/- 6dB (+/- 9 digitos)

SUBTOTAL: 620.00
IGV 19% 117.80
VALOR TOTAL (USD): 737.80

Información adicional:
Condiciones de Pago: 20% Anticipo, 80% contra entrega en Lima
Lugar de entrega: Lima
Tiempo de entrega: 1 mes
Validez de la oferta: 14 dias
Referencia bancaria: Banco del Credito del Peru BCP, Cuenta Dólares: 194-1791173165

Atentamente:

Olaf Schwetje, Apoderado ProViento S.A.C.

ProViento S.A.C Energías Renovables Perú

Av. Jose Pardo 540 y Ca. Gral. Bogoño, Miraflores, Lima - Perú
Tel#1: +51-1-4454386, Tel#2: +51-1-3029300, peru@proviento.com www.proviento.com.pe
 PROVIENTO S.A.C.
RUC: 20521134586

Proforma #: 12-0309 Fecha: 7-Mar-09

Datos del cliente:

Centro de Tecnologias de Informacion y Comunicaciones Universidad Nacional de
Nombre/Razon Social:
Ingenieria (CTIC-UNI)
Dirección: Direccion Av. Túpac Amaru 210 - Rimac / Lima 25 - Perú
Responsable / Attn.: Royer Ticse Torres

Proforma por:
Equipos Medición: x Equipos Generación: Servicios internos: Servicios externos:

Cantidad: Valor (USD)

Descripción Item:
# En Stock Unidad Total
Maletín MK30 con medidores de contaminación por ondas
1 Pcs electro-magneticas de baja y alta frecuencia (equipo semi- 1,299.00 1,299.00
profesional)
Este juego de equipos semi-profesionales mide la
intensidad de ondas electromagnetcias causados por
electrodomesticos (baja frecuencia) antenas
celulares, WLAN etc. en nT/Vm o respectivamente
µW /m2 . El display de 3,5 digitos está respaldado
por una indicación acústica que es proporcional a la
intensidad de las ondas y sirve detectar tambien
ondas pulsadas. La antena es logarítmica periodica y
tiene caracteristica de dirección.

Equipo de frecuencia baja:

Rango de frecuencia: 5Hz - 100 kHz,
Sensibilidad: 1 - 2000 nT /Vm
Precisión: +/- 2% (+/- 14 dígitos)

Equipo de frecuencia alta:

Rango de frecuencia: 800MHz - 2500 MHz (GSM,
DECT, WLAN, BLUETOOTH, CDMA,...),
Sensibilidad:0,01 - 19990 µW /m2
Precisión: +/- 6dB (+/- 7 digitos)

SUBTOTAL: 1,299.00
IGV 19% 246.81
VALOR TOTAL (USD): 1,545.81

Información adicional:
Condiciones de Pago: 20% Anticipo, 80% contra entrega en Lima
Lugar de entrega: Lima
Tiempo de entrega: 1 mes
Validez de la oferta: 20 dias
Referencia bancaria: Banco del Credito del Peru BCP, Cuenta Dólares: 194-1791173165

Atentamente:

Olaf Schwetje, Apoderado ProViento S.A.C.

ProViento S.A.C Energías Renovables Perú
Av. Jose Pardo 540 y Ca. Gral. Bogoño, Miraflores, Lima - Perú
Tel#1: +51-1-4454386, Tel#2: +51-1-3029300, peru@proviento.com www.proviento.com.pe
 www.pce-iberica.es

Calle Mayor, 53 Bajo

Angelo Romero Catay 02500 Tobarra
Av. Tupac Amaru 210-Rimac Lima 25 Albacete - España
Tel. : +34 967 543 548
Fax : +34 967 543 542
info@pce-iberica.es
PERU www.pce-iberica.es

Confirmación de Pedido
Fecha : 03/03/2009
Forma de pago : con anticipo
Fecha de entrega : 1-2 Dias
C.I.F. / N.I.F. N.° 42471455
Comprador : Angelo Romero
Su N.° de Pedido: Por Email
Fecha de Pedido: 03/03/2009 Para alguna pregunta estamos a su disposición.
Su N.° de Cliente: 23507 Tel. +34 967 543 548
Fax +34 967 543 542
E- Mail info@pce-iberica.es

Pos. Descripción N.° de Articulo Unidades Precio Precio Total

1 Medidor de radiación PCE-EM 29 PCE-EM 29 1 199,00 € 199,00 €

Importe Neto 199,00 €

Gastos de Envio 0,00 €
Importe Neto Total 199,00 €
Impuesto I.V.A. (0,00 %) 0,00 €

Bruto Importe Total 199,00 €

Bruto Importe En Pesetas 33.110,81 Pst.

GRACIAS POR SU COMPRA

C.I.F. B-02363497, R.M. de Albacete Tomo 771, Libro 535, Folio 107 Secc. 8, Hoja AB-13223, Inscrip.1
Dirección Postal PCE Group Ibérica S.L. Datos Bancarios
Calle Mayor, 53 Bajo Gerente BANCAJA (Caja de Ahorros de Valencia)
02500 Tobarra (Albacete) José Garcia Guerrero ENT.: 2077 SUC.: 0326 03
N.° Cuenta.: 310 081 0763
CODIGO BIG/SWIFT: CVALESVVXXX
CODIGO IBAN: ES26 2077 0326 0331 0081 0763
Bibliografı́a

[1] Brueckmann H Antenna Pattern Measurement by Satellite, IEEE

Transaction on Antennas Propagation,AP-11,pp. 143-147 (1963)

[2] Dr. Carlos Javier Solano Salinas Proyecto para la medición de radiación
electromagnética no ionizante de antenas de estaciones base.

[3] 2006, Limehouse Book Sprint Team Redes Inalámbricas en los Paı́ses
en Desarrollo

[4] Free antenna designs

http://www.freeantennas.com/

[5] http://www.lessemf.com/

[6] http://www.ntecusa.com/sales/dsp_model.cfm?modelID=2981

[7] http://spanish.testequipmentconnection.com/

62