Propagación del sonido. Impacto acústico

Propagación del sonido. Impacto acústico.

intro-impacto-sonoro

¿Qué es el sonido? Definición y unidades de medida.

  Pulsa en cada pestaña

  • El sonido es una onda de presión, que se transporta energía sin transporte de masa por medio de un medio elástico. 
  • Normalmente, nos interesará la propagación de estas ondas a través del aire.
  • Llamaremos presión sonora a la desviación local de la presión con respecto a la presión media o de equilibrio causada por la onda sonora $$p_{os} = p - p_a$$
Onda presión sonora
  • Esta desviación de la presión suele medirse mediante micrófono o un hidrófono
Underwater microphone (CLI & HannesGrobe, CC-BY-SA)
sonido-y-presion-sonora

  • Para medir las ondas sonoras se suele definir un nivel sonoro en escala logarítmica a partir de una presión de referencia $p_\text{ref} = 20\cdot 10^{-6} Pa$, que se considera el umbral de audición humana (en realidad a 1kHz), de la siguiente formato:

$$L_P = 10 \log_{10} \left( \frac{p_{os_\text{rms}}^2}{p_\text{ref}^2}\right) = 20 \log_{10} \left( \frac{p_{os_\text{rms}}}{p_\text{ref}}\right) dB$$

siendo $p_{os_\text{rms}} = \left(\frac{1}{T}\int_0^T p_{os}(t)^2 dt\right)^{1/2}$

  • A veces también se define a partir de la potencia sonora como $L_w = 10 \log_{10} \left( \frac{W_{os}}{W_\text{ref}}\right)$, siendo la potencia de referencia $W_\text{ref} = 10\cdot 10^{-12}{W}$
  • Se suele restringir el análisis a las frecuencias audibles, concretamente entre los 60Hz y los 8kHz
  • Promediado A, consiste en el filtrado espectral de la presión según la figura:
Figura con las funciones de promediado A (azul), B (amarillo), C (rojo) y D (negro)
nivel-sonoro

Ejemplo de procedimiento para la determinación del promediado A en un intervalo $T$ de tiempo

  1. Se elige los siguientes centros para las bandas de las octavas (frecuencias) a analizar: $$[62.5,\, 125,\, 250,\, 500,\, 1000,\, 2000,\, 4000,\, 8000]\,(Hz)$$
  2. Se fijan los límites de las bandas de frecuencia de cada octava $$[44.2,\, 88.4,\, 176.8,\, 353.6,\, 707.1,\, 1414,\, 2828,\, 5657,\, 11314]\,(Hz)$$
  3. Se aplica a la señal de $p_{os}(t)$ un filtro pasabanda centrado en la frecuencia $f$ de la octava correspondiente obteniendose la presión sonora en dicha banda de frecuencia $p_f(t)$
  4. Determinación del nivel sonoro medio para cada una de las octavas $L_{fT} = 10 \log_{10}\left(\frac{1}{T}\frac{\int_0^T p_f^2 (t) dt}{p_\text{ref}^2} \right) dB$
  5. La función de promediado $A_f$ en dB para estas frecuencias vendría dada por 
    Frec. ($Hz$) 62.5 125 250 500 1000 2000 4000 8000
    $A_f (dB)$ -26.2 -16.1 -8.6 -3.2 0 1.2 1 -1.1
  6. A cada nivel sonoro de cada banda le sumamos el valor en dB de la función de promediado A $A_f$:
    $$L_{Aft} = L_{fT} + A_f$$
  7. Agregamos todas las bandas para calcular el nivel sonoro promediado A en el intervalo T $$L_{AT} = 10 \log_{10} \left( \sum_{j=1}^{8} 10^{0.1\cdot L_{Af_jT}}\right)$$
  • Para más detalles o una descripción más precisa acudir a IEC 651 y IEC 1998-1
  • Este valor $L_{AT}$ es el que suele recoger las normativas y regulaciones
Procedimiento para el cálculo del promediado A del nivel sonoro

Se ha medido durante 10 minutos el valor de $\int_0^T p_{os}(t)^2 dt$ de un martillo neumático,  ascendiendo a $9.6\cdot 10^6 Pa^2 s$.

 Resuelve el ejercicio e indica las respuestas en el campo indicado

Se ha medido el  valor rms de la presión sonora de un compresor en las 8 bandas de frecuencia obteniéndose los siguientes valores:

 Frec. (Hz) | 62.5 | 125 | 250 | 500 | 1000 | 2000 | 4000 | 8000
 $p_{\text{os}_{rms}} (\mbox{Pa})$  | 0.45 |  0.63 |   0.89 |   1.78 |   2.00 |  7.10 |   3.99 |  1.59

 

Determina:

a) El nivel sonoro de para cada banda u octave

b) El nivel sonoro promediado A ($L_{AT}$)

 Resuelve el ejercicio e indica las respuestas en el campo indicado

El nivel sonoro de una turbina viene dado por:

Frec. (Hz) | 62.5 | 125 | 250 | 500 | 1000 | 2000 | 4000 | 8000
$L_{fT}$ (dB) | 122 | 116 | 116 |115 |97 | 85 | 60 | 57

 

Determina:

a) El nivel sonoro promediado A ($L_{AT}$)

b) El valor RMS de la presión sonora para la octava centrada en 1kHz

 

 

Resolución: Enlace a pdf

Resolución ejercicio 3 (incluye al final resolución con GNU Octave)

Regulaciones sobre el nivel sonoro

  Pulsa en cada pestaña

  • Organización Mundial de la Salud (OMS - WHO)
    En su documento de recomendaciones para el nivel sonoro nocturno en Europa (2009) recomendó fijar como límite los 40 dB, estableciendo un objetivo intermedio/transitorio de 55 dB aplicable por ejemplo en aquellos lugares con niveles altos de ruido ambiental, o donde resulte difícil su aplicación rápida.

  • Región de Murcia
    Regula los niveles de ruidos en el Decreto 48/1998 de Protección del Medio Ambiente frente al ruido. 

  • Atención: Ambas regulaciones/propuestas están basadas en la afección del ruido sobre la salud humana. Dependiendo del proyecto considerado, podría ser necesario la evaluación del efecto del ruido en los ecosistemas.

El Decreto 48/1998 de Protección del Medio Ambiente frente al ruido regula los límites y condiciones que deben verificarse para limitar el impacto acústico.

Aunque se aconseja leerlo completamente, se muestran aquí algunos de los contenidos principales de la ley:

  • Artículo 3:
    Todo proyecto de actividades e infraestructuras sometidas a evaluación de impacto ambiental o al procedimiento de calificación ambiental, susceptibles de producir impacto por ruido, adaptará y diseñará las medidas correctoras, en su caso que garanticen que el nive de lruido recibido por los receptores y usos del suelo afectados no supere los límites especificados en el Anexo I, en cuanto a medio ambiente exterior, y los especificados en el Anexo II en cuanto al interior de los edificios. 

  • Artículo 5.

    Las medidas correctoras necesarias para cumplir con el Anexo I se limitarán al control del ruido en la fuente o en su propagación.

  • Artículo 10.

    No se permitirá el funcionamiento de actividades, máquinas o instalaciones que generen un nivel sonoro en el interior de las viviendas colindantes o receptores superior al señalado en el Anexo II.

  • Artículo 11.

    1. No se permitirá el funcionamiento de actividades, máquinas o instalaciones cuyo nivel sonoro exterior a las viviendas, en patios de manzana cerrados, sea superior a 45 dB(A) durante la noche y a 55 dB(A) durante el día.

    2. No se permitirá el funcionamiento de actividades, máquinas o instalaciones cuyo nivel sonoro exterior sea superior a los niveles establecidos en el Anexo I.

    3. Cuando el nivel de ruido de fondo en la zona de ubicación sea superior a los valores del Anexo I éste podrá ser considerado por la Administración autorizante con carácter excepcional como nuevo valor de referencia a no superar.

  • ANEXO 1.
    Valores límite de ruido en el medio ambiente exterior (LeqdB(A))

Actividad

Día

Noche

Sanitario, docente, cultural, espacios naturales protegidos, parques públicos y jardines locales

60

50

Viviendas, residencias temporales, áreas recreativas y deportivas no masivas

65

55

Oficinas, locales y centros comerciales, restaurantes, bares y similares, áreas deportivas de asistencia masiva

70

60

Industria, estaciones de viajeros

75

65

Definiéndose como día las horas entre las 7:00 y las 22:00 h, y como noche las restantes

 

  • ANEXO II.
    Valores límite de ruido en el interior de los edificios.
- Día (dB) Noche (dB)
Sanitario, Docente y Cultural   45     35  
Viviendas y hoteles   50     40  

 

Propagación del ruido

Principales métodos o estándares de cálculo

  • Norma ISO 9613-2:1996 "Acoustics --Attenuation of sound during propagation outdoors-- Part 2 General Method of Calculation".
    • Actualización ISO 9613-2:2024
  • Common noise assessment methods in Europe (CNOSSOS-EU), recogida en el anexo II de la Directiva 2002/49/EC
  • Otros métodos más avanzados.

 

Proceso de cálculo de acuerdo a la norma ISO 9613-2

Propagación Sonido Intro

  Pulsa en cada pestaña

La norma presupone unas condiciones meteorológicas moderadas pero favorables a la propagación del sonido entre la fuente y el receptor. A estas condiciones se les denomina dentro de la norma DW (de ``downwind''). 

Las condiciones son las siguientes:

  • Dirección de viento $\pm 45º$ respecto la dirección que une fuente y receptor. Por ello se suele apodar el ruido calculado como DW, al estar el receptor ``aguas abajo'' de la fuente.
  • Velocidad del viento entre 1 y 5 m/s a una altura entre 3 y 11m sobre el nivel del suelo.
    • Solo en el caso de aerogeneradores se recogen otras velocidades (Anexo D de la nueva edición de la norma 2024).

 

Si desea conocerse el comportamiento medio durante un período amplio donde hayan condiciones favorables y desfavorables puede corregirse el valor DW calculado mediante la estimación a largo plazo LT recogida en la norma. 

Propagación del sonido: Condiciones meteorológicas (DW, LT)

La norma asume que tenemos un receptor (r) y una o varias fuentes (s de "source").

  • Agregación de multiples fuentes en una fuente equivalente.
    Un grupo de fuentes puntuales puede ser descrita como una fuente equivalente puntual localizada en el centro del grupo si se cumplen las siguientes condiciones:
    • Las fuentes tienen una intensidad y altura sobre el suelo similares.
    • Las condiciones de propagación desde cada fuente al receptor son similares.
    • La distancia desde el punto de la fuente agregada equivalente hasta el receptor es más del doble de la máxima distancia en el conjunto de fuentes.

  • Agregación del nivel sonoro recibido desde $N_s$ fuentes:
    \[L_{AT} = L_{AT_1} \oplus L_{AT_2} \oplus \dots \oplus L_{AT_{N_s}} = 10 \log_{10} \left(\sum_{n=1}^{\text{N_s}} \sum_{j=1}^{8} 10^{0.1\cdot L_{{Af_jT}_n}}\right) \]

  • Fuentes líneales y superficiales pueden ser divididas en fuentes puntuales desde las que se puede calcular la propagación del sonido. Esas fuentes puntuales equivalentes deben estar a una distancia del punto de recepción de más del doble de la dimensión máxima que representan y tener las mismas condiciones de propagación. Si las condiciones de propagación difieren por ejemplo por la existencia de barreras, estas deberán tenerse en cuenta previamente para dividir el elemento.
Esquema de división de una fuente lineal en un conjunto de fuentes puntuales (figura de la norma ISO)
  Pulsa en la imagen
Propagación del sonido: Agregación de fuentes. Fuentes lineales y superficiales.
Propagación sonido: Procedimiento simplificado

Podemos calcular el nivel sonoro propagado ($L_{fT}(DW)$) hasta el receptor para cada banda de frecuencia como:
\[
    L_{fT}(DW)=L_{Wf} + D_c - A_{b_f}
\]
donde $L_{Wf}$ es el nivel de sonido en la banda emitido por la fuente, $D_c$ la corrección direccional y $A_{b_f}$ es la atenuación en la banda en la propagación de la fuente al receptor dada por la ecuación.


Esquema del proceso de cálculo de la atenuación para cada banda

 

 

  • Corrección direccional: $$D_c = \underbrace{D_I}_{\scriptsize\mbox{índice direccional}} + \underbrace{D_\Omega}_{\scriptsize\mbox{ índice de propagación para ángulos sólidos }<4\pi}$$

    Si la fuente es omnidireccional y la propagación es libre (áng. sólido $=4\pi$), entonces $D_c =  0 dB$. En general será $D_c = 10 log_{10}\left(\frac{4\pi}{\Omega}\right) $dB donde $\Omega$ es el ángulo solido en la que se propaga el sonido. 
    Por ejemplo, si una superficie cercana bloque la propagación, el angulo sólido restante sería la mitad $\Omega=\frac{4\pi}{2} \,\rightarrow D_c \approx 3.0 dB$. Nota: En este caso, la reflexión en esa superficie no debería tenerse en cuenta de forma independiente, al estar ya incluida.


  • La atenuación en la banda puede descomponerse como:

\[A_{b_f} = \underbrace{A_{div}}_{\scriptsize\mbox{divergencia geométrica}} +
    \underbrace{A_{atm}}_{\scriptsize\mbox{absorción atmosférica}} +
    \underbrace{A_{gr}}_{\scriptsize\mbox{efecto suelo}} + \underbrace{A_{bar}}_{\scriptsize\mbox{pantalla/barrera}} +
    \underbrace{A_{misc}}_{\scriptsize\mbox{Otros}}
  \]

 

Esquema del proceso de cálculo
Propagación del sonido: Atenuación

  • Tiene en cuenta la reducción de la potencia sonora en la onda debido a su expansión.
  • Supone propagación esférica desde la fuente.
  • Los efectos de reflexión o suelo son incluidos posteriormente.

Así, a una distancia $d$, la potencia sonora $W_{os}$ estaría repartida en un área de  $4\pi d^2$. $A_{div}$ sería la diferencia entre la potencia sonora emitida y la recibida a una distancia $d$:
$A_{div}=
  10\log_{10}\left(\frac{W_{os}}{W_{\text{{ref}}}}\right)-10\log_{10}\left(\frac{\frac{W_{os}}{4\pi
          d^2}}{W_{\text{{ref}}}}\right)=
\cancel{10\log_{10}\left(\frac{W_{os}}{W_{\text{{ref}}}}\right)}-\left(\cancel{10\log_{10}\left(\frac{W_{os}}{W_{\text{{ref}}}}\right)}
- 10\log_{10}\left(4\pi d^2\right)\right) =
\underbrace{10\log_{10}\left(4\pi\right)}_{10.99\approx 11}+ 20\log_{10}\left(d\right)
$

\[A_{div} = 20 \log_{10} \left(\frac{d}{d_0}\right) + 11\; dB)\]

donde $d$ es la distancia entre la fuente y el receptor y $d_0 = 1m$ es la distancia de referencia

 

Adiv

  • Atenuación por absorción atmosférica (transporte, rotación, vibraciones moleculares, etc.)
    \[ A_{atm} = \frac{\alpha \cdot d}{1000}\]
     donde $\alpha(f,\,H.R.,\,T,\,p_a)$ es el coeficiente de atenuación atmosférica en dB/km, tal y como se indica en la norma ISO 9613-1. Los parámetros con más influencia en $\alpha$ son la frecuencia y la humedad relativa. Algunos ejemplos de dicho coeficiente se muestran en la tabla siguiente que aparecía como ejemplo en la ISO9613-2:1996:

Tabla con valores de alpha Aatm en función de la temperatura y la humedad relativa que viene como ejemplo en ISO9613-2:1996 que viene en realidad de ISO9613-1
Aatm

  • Es el efecto de la interferencia entre el sonido propagado directamente de fuente a receptor y el reflejado en el suelo.
  • Debido a la dirección de reflexión y atenuación, solo el suelo cercano a la zona de recepción y de la fuente es el más relevante.
  • Se considera suelo plano o con pendiente aproximadamente constante.

¡Atención! $A_{gr}$ contiene tanto el efecto de atenuación como de reflexión por lo que su signo puede cambiar.

\[     A_{gr} = \underbrace{A_s}_{\scriptsize\mbox{  región fuente}}+ \underbrace{A_r}_{\scriptsize\mbox{  región receptora}}+ \underbrace{A_m}_{\scriptsize\mbox{  región intermedia}}\]

Figura de la norma que explica la distribuciones de las regiones en las que se divide el suelo bajo el canal de tranmisión.

$A_{s,r,m}(f,\,G_{s,r,m})$, donde $G$ es el factor del tipo de suelo

  • Suelo duro/compacto $\Rightarrow G=0$
        
        Pavimentos, agua, hielo, cemento, y en general, todos los suelos con porosidad baja.

  • Suelo poroso $\Rightarrow G=1$

    Suelo cubierto por hierba, árboles u otra vegetación, incluyendo todos los suelos adecuados para la producción de vegetales, pastos, etc.

  • Suelo mixto  $0 <G <1$

    La región incluye parte de suelo poroso y parte compactado, $G$ estará entre 0 y 1, empleándose como valor de $G$ la proporción de suelo poroso en la región.

El método CNOSSOS-EU incluye está tabla para ayudar en la asignación del coeficiente al suelo:

Tipos de suelo de acuerdo al CNOSSOS-EU
  •  $A_s(f,\,G_s)$ y $A_r(f,\,G_r)$ pueden calcularse según la frecuencia con las siguientes expresiones

     

    Freq   (Hz) $A_{s,r}\, (dB)$
    63  $-1.5$
    125 $-1.5 + G\cdot\left(1.5+3 e^{-0.12(h-5)^2}\left(1-e^{-d_p/50}          \right)+5.7e^{-0.09h^2}\left( 1-e^{-2.8\cdot 10^{-6} d_p^2}          \right) \right)$
    250  $-1.5 + G\cdot\left(1.5+8.6 e^{-0.09h^2}\left(1-e^{-d_p/50} \right) \right)$ 
    500 $-1.5 + G\cdot\left(1.5+ 14.0e^{-0.46h^2}\left(1-e^{-d_p/50} \right) \right)$
    1000 $-1.5 + G\cdot\left(1.5+ 5.0e^{-0.9h^2}\left(1-e^{-d_p/50} \right) \right)$
    2k-8k $-1.5(1-G)$

     




  • $A_m(f,\,G_m)$ sólo existirá (será $\neq 0$) si $30(h_s+h_r) < d_p$

    $A_{s,r}$(dB) Freq. (Hz)
    63 $-3 \left(1-\frac{30(h_s+h_r)}{d_p}\right)$
    125-8000  $-3 \left(1-\frac{30(h_s+h_r)}{d_p}\right)\left(1-G_m\right)$

     



Nota: ¡La nueva revisión ISO9613-2:2024 va a introducir algunos cambios en este término y el de la difracción!

 

Agr

  • Consideraremos que tenemos una barrera o pantalla acústica cuando una barrera cumpla con las siguientes condiciones:
    • Densidad "superficial" $\geq 10kg/m^2$
    • Superficie cerrada sin huecos de tamaño significativo
    • La longitud de su proyección horizontal sobre la dirección $\perp$ a la línea fuente-receptor es mayor que la londigud de onda de la octava corespondiente dada por $\lambda \approx \frac{340m/s}{f}$:
            $$ l_l + l_r > \lambda$$

Figura de la norma explicativa sobre la última condición:

Esquema de la longitud mínima de una pantalla acústica


Podrían quedar como vías de propagación por difracción uno o dos ejes verticales de la pantalla y el eje superior de la misma.

Si solo quedase una vía de propagación:

  • Difracción a través de un eje vertical:
    \[A_{bar} = \max\left\{D_z,\,0\right\}\]
    con $D_z$ dado por
    \[D_z =10\log_{10}\left(3+\frac{C_2}{\lambda}C_3 z k_\text{met}\right)\]
    donde $C_2 = 20$ e incluiría los efectos de reflexión en el suelo, si estos fuesen a introducirse como "fuentes imaginarias" adicionales, entonces debe emplearse $C_2 = 40$, $C_3 = 1$ para difracción simple, $z$ es un coeficiente que mide la diferencia entre las longitudes de los caminos de propagación de sonido directo (si no hubiese barrera) y con difracción, mientras que $k_\text{met}$ es un factor de corrección por efectos meteorológicos dada por
    • \( k_\text{met} = e^{-\frac{\sqrt{d_{ss}d_{sr}\frac{d}{2z}}}{2000}} \text{cuando } z>0\)
    • \( k_\text{met} = 1 \text{cuando } z\leq 0\)
    •  $k_\text{met} \approx 1$ para difracción lateral y para distancias menores que 100m.
  • Difracción a través del eje superior:
    \[A_{bar} = \max\left\{D_z-A_{gr},\,0\right\}\]
    donde el término $A_{gr}$ anula el efecto de $A_{gr}$, al ya estar incluido dentro de $D_z$



  • En el caso de difracción simple podemos calcular $z$ a partir de la siguiente ecuación:
    \[          z = \sqrt{(d_{ss}+d_{sr})^2+a^2}-d \]

  • Si tenemos doble difracción (por ejemplo en una barrera bastante ancha o dos barreras separadas)
    \[          z = \sqrt{(d_{ss}+d_{sr}+e)^2+a^2}-d \]
    y además, deberemos cambiar el valor de $C_3$ que será mayor que 1, concretamente tendremos $1<C_3\leq 3$:
    \[           C_3 =   \frac{1+\left(\frac{5\lambda}{e}\right)^2}{\frac{1}{3}+\left(\frac{5\lambda}{e}\right)^2}\]
Fig 6 y 7 calculo de z pantallas acusticas ISO9613-2:1996

 

  • En todo caso, el efecto máximo vendría limitado por $D_z \leq 20 dB $ para difracción simple o $ D_z \leq 25dB$ para difracción doble.

  • Si hay más de dos barreras, se puede calcular tomando las difracción doble entre aquellas dos que sean más efectivas y despreciar el resto.

  • Si hay elementos laterales que reflejen, la atenuación puede ser menor (ver reflexión).

  • En casos, como una autovía por debajo del nivel del suelo la atenuación real sería mayor que la aquí calculada.

  • En el caso de que la línea una fuente y receptor pase por encima del borde superior de la barrera, deberemos cambiarle el signo al coeficiente $z$ definido anteriormente.
Fig 5 ISO 9613-2:1996 varios ejes de difracción pantalla

 

  • Si hubiese varios caminos de difracción, se debe calcular los niveles de sonido por separado por cada eje de difracción agregándolas posteriormente (como si fuesen fuentes distintas).
Propagación sonora: Abar

Nota: el vídeo contiene una errata, cuando habla de z negativo. Z no puede salir negativo calculando, pero si la pantalla no interrumpe la vía de propagación directa, una vez calculado z se le cambia el signo. 

El vídeo será revisado en una próxima edición. 

En cualquier caso, mientras tanto si queréis escuchar la explicación también aparece mientras hago el apartado c del ejercicio 1.

  • En este apartado se quiere tratar el efecto de reflexión que pueden provocar paredes más o menos verticales, techos, cubiertas, etc. que pueden reflejar el sonido hacia el receptor

  • No se incluye la reflexión en el suelo pues ya esta considerada

  • Para que se necesario incluir sus efectos se deben dar estas condiciones:
    • El sonido reflejado puede alcanzar el receptor con un ángulo de incidencia ($\beta$) idéntica al ángulo reflejado
    • Coeficiente de reflexión de la superficie $> 0.2$
    • Longitud de la superficie suficientemente grande ($l>l_\text{min}$) para reflejar en la octava o banda de frecuencia considerada:
      \[      \frac{1}{\lambda} > \frac{2}{\left(l_\text{min} \cos\beta\right)^2}\frac{d_{so}d_{ro}}{d_{so}+d_{ro}}\]

 

Figura 8 norma ISO 9613-2:1996 muestra el esquema de una reflexión

 

  • En ese caso se debe simular los efectos de reflexión mediante una ``fuente imaginaria'' situada en el punto $S_i$, calculando la propagación hasta el receptor, y agregando sus efectos posteriormente.

  • La fuente imaginaria tendrá un nivel de potencia sonora dada por

\[
    L_{wf,im} = L_{wf} + 10 \log_{10} \rho_{rc} + D_{Ir}
  \]
donde $D_{Ir}$ es el índice direccional de la fuente, teniendo en cuenta el camino real (hacia la superficie de reflexión), y $\rho_{rc}$ sería el coeficiente de reflexión del sonido para el ángulo de incidencia $\beta$.

Si no está disponible (ensayos, etc.), se puede estimar de esta tabla:
  

 

$\rho_{rc}$ Tipo de superficie
Muro planos y compactos 1
Muros de edificios con ventanas 0.8
Muro de fabrica con el 50% de su superficie formada por espacios abiertos, instalaciones, tuberías, … 0.4
Instalaciones totalmente abiertas (tuberías, torres, etc.) 0

 

  • La atenuación desde $S_i$ hasta el receptor deberá calcularse según el procedimiento normal.
Propagación sonora: Reflexión

  • La norma cubre tres tipos de efectos:
    1. Atenuación a través de masas arbóreas $A_{fol}$
    2. Atenuación a través de un emplazamiento industrial $A_{site}$
    3. Atenuación a través de zonas residenciales $A_{hous}$

En esta presentación cubriremos la forma simplificada de $A_{fol}$. La nueva versión de la norma ha añade una forma más detallada de estimación de la atenuación a través de masas arbóreas cuando se dispone de un mayor número de datos sobre la misma.

 

 

Figura A1 de la norma ISO9613-2:1996  (distancias masa arbórea)
  • La masa arbórea debe ser suficientemente denso para bloquear totalmente el camino directo.
        
    Orientativamente, se puede considerar que se da esta circunstancia cuando es imposible ver una distancia corta a través del mismo

  • La vegetación puede atenuar por estar
    • Cercana a la fuente
    • Cercana al receptor
    • En ambas situaciones

  • Para calcular $d_1$ y $d_2$ puede trazarse una curva con un radio de curvatura de 5km como en la figura anterior, o bien aproximarse mediante líneas a 15º del suelo

  • La distancia de propagación a través de masa arbórea queda definida como
    \[d_f = d_1 + d_2\]

  • Si $10\leq d_f \leq 20m$

    Frec. (Hz) 63 125 250 500 1000 2000 4000 8000
    $A_{fol}(dB)$ 0 0 1 1 1 1 2 3

     

  •  Si $20\leq d_f \leq 200m$

    Frec. (Hz) 63 125 250 500 1000 2000 4000 8000
    $\frac{A_{fol}}{d_f} (dB/m)$ 0.02 0.03 0.04 0.05 0.06 0.08 0.09 0.12

 

Propagación sonora: Afol

  • Cuando se quiere calcular $L_{AT}(LT)$ a largo plazo las condiciones meteorológicas variarán en el tiempo, por ejemplo, el viento será siempre favorable a la transmisión.
  • Para obtener un valor más cercano a la media en un período largo podemos emplear esta expresión (suponiendo que el sonido es cte. en el tiempo)
    \[    L_{AT}(LT) =  \underbrace{L_{AT}(DW)}_{\mbox{calculado anteriormente}} - C_\text{met}
    \]
    donde \(C_\text{met} \) será 
    • \( C_\text{met} = 0  \text{ si } d_p\leq 10(h_s+h_r)\)
    • \( C_\text{met} = C_0\left[1-10\frac{h_s+h_r}{d_p}\right] \text{ si } d_p > 10(h_s+h_r)
      \)
      siendo $C_0$ un factor que depende de las variaciones del viento y las temperatura.
      • $0\leq C_0\leq 5 dB$, normalmente será $C_0 \leq 2dB$
      • $C_0$ podría ser aportado por autoridades locales o estimarse a partir de la climatología.

 

Propagación sonora: Estimación a largo plazo

Aunque hay similitudes entre el método descrito y la guía de Common noise assessment methods in Europe (CNOSSOS-EU), por ejemplo en cuanto a la agregación de fuentes, cálculo de atenuación por absorción atmosférica o divergencia geométrica, y en muchas definiciones adicionales, por ej. relativas al tipo de suelo, existen ciertas diferencias. Por ejemplo:

  • Propone un método diferente para el cálculo de la atenuación por efecto suelo, difracción o reflexión.
  • CNOSSOS-EU contiene distintas propuestas para estimar fuentes de ruido habituales como las de tráfico rodado, ferroviario o aereo.

Ejercicios y herramientas de apoyo

  Pulsa en cada pestaña

Enunciado del ejercicio junto a resolución.

A continuación se muestra la resolución de cada apartado en un vídeo.

a) Cálcula el nivel sonoro que llegaría para cada banda y el ponderado A despreciando cualquier 
efecto de reflexión en el edificio de oficinas.

Ejercicio 1. Apartado a.

b) Se producirá reflexión en el edificio de oficinas? Recalcula el nivel sonoro A en el 
restaurante teniendo en cuenta dicho efecto.

Ejercicio 1. Apartado b.

b) Se producirá reflexión en el edificio de oficinas? Recalcula el nivel sonoro A en el 
restaurante teniendo en cuenta dicho efecto.

Ejercicio 1. Apartado c

"Script" de GNU Octave con la resolución del ejercicio:

###############################################
cargafuncionesruido
###############################################
% Preproblema obtenemos ruido emitido 1º a partir de datos a 10m
Lf10m_excavadora40t = [85 78 77 77 73 71 68 63];
hs = 0.15;
hr_10m = 0.15;
%we negglect Aatm 10m
Lw_excavadora = Lf10m_excavadora40t + Adivf(10)+Asorf(10,hs,1)+Asorf(10,hr_10m,1)+Amf(10,hs,hr_10m,1)*(30*(hs+hr_10m)<10);

###############################################
### Aqui empieza la resolución del ejercicio
###############################################
## Lw_excavadora =

##    112.700   109.068   111.112   113.023   105.776   102.000    99.000    94.000

%% Primero solo nivel sonoro directo!
hr = 1.5;
dp = 160;
d = sqrt(dp^2+(hs-hr)^2)
Gs = 0;
Gr = 0;
Gm = 0;
%Adiv
Adiv = Adivf(d)
%Aatm 15 deg, 70% % justo media entre 10 y 20 deg 70%
alpha = ([0.1, 0.4, 1.0, 1.9, 3.7, 9.7, 32.8, 117] + [0.1, 0.3, 1.1, 2.8, 5.0, 9.0, 22.9, 76.6])/2
Aatm = Aatmf(alpha,d)
Asorf(dp,hs,Gs)
Asorf(dp,hr,Gr)
Amf(dp,hs,hr,Gm)*(30*(hs+hr)<dp)
Agr = Asorf(dp,hs,Gs) + Asorf(dp,hr,Gr) + Amf(dp,hs,hr,Gm)*(30*(hs+hr)<dp)
Abf = Adiv + Aatm + Agr

% Nivel sonoro que llega (prop. directa)
LfT = Lw_excavadora - Abf

% Ponderación A
LAfT = LfT + Af
LAT = 10*log10(sum(10.^(LAfT/10)))


% reflexión
lambda = 340./f
dso = sqrt(25^2+80^2)
dro =  sqrt(25^2+80^2)
beta = acos(25/dso);
beta*180/pi
lmin = sqrt(lambda*2/cos(beta)^2*dso*dro/(dso+dro))
rhorc = 0.8
% solo banda 8
Lwf8im = Lw_excavadora(8) + 10*log10(rhorc)
dref = dso+dro
Adivref = Adivf(dref)
Aatmref = Aatmf(alpha,dref)
Agrref = Asorf(dref,hs,Gs) + Asorf(dref,hr,Gr) + Amf(dref,hs,hr,Gm)*(30*(hs+hr)<dref)
Abfref = Adivref + Aatmref(8) + Agrref(8)
Lf8ref = Lwf8im - Abfref
LAf8Tref = Lf8ref + Af(8)
LAT = 10*log10(sum(10.^(LAfT/10)) + 10.^(LAf8Tref/10)) 


% apartado c) Pantalla acústica
C2 = 20
C3 = 1
dss = sqrt(16^2 + (2-0.15)^2)
dsr = sqrt((160-16)^2 + (2-1.5)^2)
a = 0
z = sqrt((dss+dsr)^2+a^2)-d
kmet = e^(-sqrt(dss*dsr*d/2/z)/2000)
Dz = 10*log10(3+C2./lambda*C3*z*kmet)
Abar = Dz-Agr
Abar(Abar<0) = 0

Abf_pantalla = Adiv + Aatm + Agr + Abar

% Nivel sonoro que llega (prop. directa)
LfT_pantalla = Lw_excavadora - Abf_pantalla

% Ponderación A
LAfT_pantalla = LfT_pantalla + Af
LAT_pantalla = 10*log10(sum(10.^(LAfT_pantalla/10)))

Base de datos de máquinaria de obra civil. 

¡¡Atención son niveles sonoros medidos a 10m no los emitidos (requiere estimar atenuación a 10m para calcular los originales antes de usarlos)!!

Enlace a pdf.

Este script precarga algunas funciones a GNU Octave para agilizar cálculos o comprobar lo hecho a mano:

cargafuncionesruido.m

 

f=[62.5 125 250 500 1e3 2e3 4e3 8e3];
Af = [-26.2  -16.1  -8.6  -3.2  0  1.2  1  -1.1];

 

function adiv = Adivf(d)
  adiv = 20*log10(d)+11;
endfunction

 

function aatm = Aatmf(alfa,d)
  aatm = alfa*d/1000;
endfunction

 

function asr = Asorf(dp,h,G)
  asr(1) = -1.5;
  asr(2) = -1.5 + G*(1.5+3*e^(-0.12*(h-5)^2)*(1-e^(-dp/50))+5.7*e^(-0.09*h^2)*(1-e^(-2.8e-6*dp^2)));
  asr(3:5) = -1.5 + G.*(1.5 + [8.6 14 5].*e.^([-0.09 -0.46 -0.9]*h^2)*(1-e^(-dp/50)));
  asr(6:8)= -1.5*(1-G);
endfunction

 

function am = Amf(dp,hs,hr,G)
  if 30*(hs+hr)<dp,
    am(1) = -3*(1-30*(hs+hr)/dp);
    am(2:8) = -3*(1-30*(hs+hr)/dp)*(1-G);
  else
    am(1:8)=0;
  end
endfunction