Revision as of 13:43, 12 July 2017

N. Frechilla - Universidad de Valladolid

M. Cacho-Pérez - Universidad de Valladolid

A. Lorenzana - Instituto de Tecnologías Avanzadas de la Producción

Resumen

Obras de ingeniería civil como las pasarelas peatonales son estructuras cada vez más esbeltas, ambiciosas y complejas en las que los efectos dinámicos producidos por las solicitaciones a las que están sometidas (peatones) pueden inducir problemas de servicio y/o amplificaciones importantes o repetitivas de los esfuerzos y deformaciones. Ésto junto con las actuales demandas sociales relativas a la percepción y confort suscitan la necesidad de integrar el análisis y diseño asistido por ordenador en la evaluación de las vibraciones generadas por las cargas a las que están sometidas. El empleo de los programas comerciales de elementos finitos no ha proliferado tanto como en otros campos de la ingeniería por la existencia de metodologías simplificativas en numerosas guías y normativas nacionales e internacionales y por la dificultad de hacer modelos realistas en el caso de estructuras complejas. Sin embargo, el conocimiento generado por el análisis de este tipo de modelos numéricos resultaría especialmente útil en las fases más tempranas del diseño de estructuras similares, tanto inicial como de las modificaciones necesarias para solucionar problemas dinámicos existentes, ya que es en estas etapas cuando se podría aplicar la mejor solución desde un punto de vista técnico y económico. En el presente artículo se describe el proceso de generación y ajuste del modelo mecánico de una pasarela actualmente en servicio y la simulación de la solicitación de dos tipos de carga: una excitación generada por un shaker y una carga peatonal, contrastando los resultados generados mediante el modelo numérico con los obtenidos experimentalmente.

Palabras clave: Pasarela peatonal, vibraciones, elementos finitos, ajuste modelo mecánico.

Abstract

Civil engineering constructions like pedestrian footbridges are structures increasingly slender, ambitious and complex in which the dynamic effects caused by the loads to which they are subjected (pedestrians) may induce service problems and / or significant or repetitive amplifications of the efforts and deformations. This fact together with the current social demands related to the perception and comfort raise the need to integrate the analysis and computer-aided design in the evaluation of the vibrations generated by the loads this type of structures supports. In the specific case of finite element commercial programs, their use has not proliferated as much as in other engineering fields due to the existence of simplified methodologies in numerous national and international guidelines and regulations and the difficulty of making realistic models in the case of complex structures. However, the knowledge generated by the analysis of this type of numerical models would be especially useful in the earlier phases of the design of similar structures, so much initial as of the necessary modifications to solve existing dynamic problems, since it is in these stages when the best solution could be applied from both a technical and an economic point of view. This article describes the process of generation and updating of the mechanical finite element model (FEM) of a footbridge in service and the simulation of two types of load: an excitation generated by a shaker and a pedestrian load, contrasting the results generated by the numerical model with those obtained experimentally.

1. Introducción

El diseño y construcción de un gran número de obras de ingeniería industrial y civil ya no está apenas limitada por los materiales –cada vez más resistentes- y los sistemas constructivos –cada vez con mayor componente tecnológico- lo que ha provocado que en la actualidad éstas sean más esbeltas, ambiciosas y complejas.

Este hecho, junto con las actuales demandas sociales relativas a la percepción y confort suscitan la necesidad de integrar el análisis y diseño asistido por ordenador en la evaluación de la respuesta dinámica de estructuras esbeltas que son las más susceptibles de sufrir movimientos o vibraciones significativas.

Un caso particularmente relevante es el de las pasarelas peatonales, ya que generalmente son propensas a ser excitadas por las personas que las ocupan, pudiendo los efectos dinámicos producidos inducir, además de problemas de servicio, amplificaciones importantes y repetitivas de los esfuerzos y deformaciones que pueden afectar a la seguridad resistente o la vida a fatiga de la estructura.

Aunque las normativas recogen diversos criterios de evaluación de estos fenómenos, también reflejan la dificultad de evaluar con precisión el comportamiento dinámico de las estructuras, exigiendo en algunos casos un estudio dinámico riguroso.

La simulación mediante programas basados en el análisis de elementos finitos es una de las herramientas de cálculo con mayor profusión ya que los actuales softwares son muy competitivos respecto al ahorro de tiempo y costes. Sin embargo en el caso del análisis dinámico de estructuras no ha proliferado tanto debido a la existencia de normativas y guías nacionales e internacionales [1-3] que proporcionan metodologías simplificativas, que siendo útiles para la acotación del problema dinámico no ofrecen la enorme versatilidad del análisis numérico.

En general las predicciones logradas con el uso de modelos analíticos, especialmente en el caso de estructuras complejas como el que nos ocupa, difieren de los resultados obtenidos mediante ensayos experimentales por lo que es necesario efectuar un ajuste del modelo de elementos finitos (denominado en inglés FEM Updating), mejorando las hipótesis efectuadas y aproximando sus propiedades a las de la estructura real. Aunque los datos experimentales poseen errores que puedan provocar una baja correlación entre éstos y las predicciones, la teoría de ajuste de modelos supone que los principales causantes del desacuerdo teórico experimental son los errores en el modelo [4].

En el caso que nos ocupa se emplean los resultados experimentales procedentes de un análisis modal operacional, en inglés Operational Modal Analysis (OMA), realizado en la pasarela completa y un análisis modal clásico, del inglés Experimental Modal Analysis (EMA), realizado en el vano más sensible a excitaciones dinámicas producidas por humanos (vano 2), para conseguir que nuestro modelo numérico se ajuste lo mejor posible a la realidad. Los parámetros modales que identifican el comportamiento dinámico de la pasarela objeto de estudio (pasarela del Museo de la Ciencia de Valladolid) se han obtenido en estudios anteriores [5-7].

Una vez ajustado el modelo analítico (creado mediante el software de elementos finitos ANSYS) se han obtenido las respuestas dinámicas simuladas, en forma de aceleraciones, al aplicar la excitación generada durante el EMA y una carga peatonal según el modelo convencional de flujo de peatones [1]. Los resultados numéricos se han contrastado con los experimentales a fin de validar el comportamiento dinámico de la estructura modelizada.

El presente artículo se ha organizado del siguiente modo: en primer lugar, tras esta introducción, se resume la metodología empleada. A continuación en el tercer epígrafe se describe la pasarela estudiada, los parámetros modales que identifican su comportamiento dinámico y los procesos de modelización y ajuste de la estructura. En el cuarto apartado se presentan ejemplos de la respuesta dinámica simulada, y se comparan con la registrada experimentalmente para finalmente presentar las principales conclusiones que se deducen de este trabajo.

2. Metodología

2.1. Modelización y ajuste

La metodología seguida en el presente trabajo ha consistido en primer lugar en la modelización de la pasarela, partiendo de la geometría de proyecto, mediante el programa de elementos finitos ANSYS y procediendo posteriormente a su ajuste empleando para ello un esquema iterativo que permitió analizar y seleccionar aquellos parámetros físicos cuyo efecto era más relevante para tal fin. Respecto al ajuste del modelo mecánico (conocido como FEM Updating en la terminología inglesa), varios autores [8-15] han publicado y comparado varias técnicas sin alcanzarse ningún acuerdo sobre la existencia de un procedimiento general que pueda aplicarse a todos los casos, por lo que la elección ha dependido de los parámetros seleccionados.

Ésta se ha realizado en función de los parámetros modales extraídos de las campañas experimentales, empleando como referencia las frecuencias y las formas modales en el caso de la pasarela completa y complementándolo con las masas modales y/o modos normalizados respecto a la matriz de masa.

Dentro de los métodos iterativos de ajuste, se ha utilizado una aproximación para correlacionar el modelo experimental y el modelo matemático que consiste en la comparación entre las frecuencias naturales [16], pudiéndose la diferencia en porcentaje (${\textstyle {\epsilon }_{f}}$) expresarse tal y como se muestra a continuación:

siendo f_Xi el parámetro experimental y f_Ai el parámetro calculado.

Para medir la correlación entre los modos normalizados respecto a la matriz de masa experimentales y los obtenidos mediante el modelo en ANSYS se ha empleado el criterio de confianza modal (MAC según la abreviatura de su denominación en inglés) que proporciona una medida del grado de linealidad entre estimaciones de los modos de vibración originados desde diferentes fuentes [16]. Se define como un escalar y su expresión se muestra a continuación:

siendo {Ψ_X} el modo de vibración experimental y {Ψ_A} el modo de vibración calculado, donde la notación {}^H significa transpuesta compleja conjugada (matriz hermítica).

2.2. Simulación de cargas

Una vez ajustado el modelo se ha estudiado el comportamiento del mismo ante dos escenarios reales diferentes. El primero de ellos ha consistido en la simulación de las fuerzas generadas durante el EMA realizado en el vano más susceptible a las vibraciones (vano 2). En el segundo supuesto se ha simulado el tránsito peatonal a lo largo de toda la pasarela. En ambos casos el objetivo es obtener la respuesta dinámica en forma de aceleraciones del modelo numérico y poder comparar los resultados con los registrados experimentalmente con el fin de validarlo.

Para simular el efecto de varios peatones (solicitación del segundo escenario estudiado) se ha empleado el modelo convencional que considera que la respuesta de un flujo de peatones (sin ninguna sincronización particular) se obtiene multiplicando el efecto de un solo peatón por un factor igual a ${\textstyle {\sqrt {N}}}$, siendo N el número de peatones presentes en la pasarela en un momento determinado [1]. La velocidad del peatón se ha considerado igual a 0.9 veces la frecuencia de paso en Hz.

La fuerza de excitación generada por un solo peatón (movimiento rítmico al caminar) se ha expresado matemáticamente mediante una serie de Fourier de la forma:

donde

G₀ = fuerza estática (componente vertical del peso de un peatón)

G₁ = amplitud del primer armónico

G_i = amplitud del i-ésimo armónico

f_m = Frecuencia o cadencia de la actividad expresada en Herzios (Hz).

ϕ_i = ángulo de desfase del i-ésimo armónico respecto al primero

n = número de armónicos

En el caso que estamos analizando, consideramos la actividad de caminar/marcha (contacto continuo con el suelo), estableciendo una frecuencia de paso en torno a 1.8 Hz.

Se ha estimado que con n=3 se obtiene una precisión adecuada, siendo los valores de los coeficientes los siguientes: G₀= 700 N, G₁= 0,4 G₀, G₂=G₃≈0,1 G₀, ϕ₂=ϕ₃≈π/2.

3. Resultados numéricos y discusión

3.1. Resultados del análisis modal experimental

La Pasarela del Museo de la Ciencia es una estructura en servicio situada en Valladolid sobre el río Pisuerga (figura 1). Fue inaugurada en 2004 y tiene una longitud total de 234 metros repartidos en cuatro vanos: tres formados por barras metálicas de 51, 111 y 21 metros de luz (vano 2, 3 y 4 respectivamente), y uno de hormigón de 45 metros (vano 1). Éste último es muy rígido por lo que no se considera en el análisis objeto del presente trabajo. La tipología resistente de los vanos metálicos es de entramado triangulado de barras tubulares con sección pretensado formado por un conjunto de cables longitudinales y transversales en el vano 3 (figura 2 y 3).

Figura 1. Pasarela del Museo de la Ciencia. Vista general.

Figura 2. Pasarela del Museo de la Ciencia. Vista del vano 3.

Figura 3. Pasarela del Museo de la Ciencia. Detalle del pretensado perimetral en el vano 3.

En el trabajo de investigación llevado a cabo por J. de Sebastián Sanz [5] se monitorizaron durante varios meses las aceleraciones verticales y laterales en el centro de cada uno de los tres vanos metálicos de la pasarela observándose que los niveles de vibración en la dirección lateral eran considerablemente menores que en la dirección vertical y que el nivel de vibraciones del vano 2 era notablemente mayor que el del vano 3. Más adelante se identificaron mediante un OMA los parámetros dinámicos de la estructura a partir de los registros de aceleración en dirección vertical monitorizados a lo largo de los tres vanos metálicos de la pasarela. Dichos parámetros modales se obtuvieron mediante el programa comercial ARTeMIS empleando tres algoritmos en el dominio de la frecuencia (FDD (Frequency Domain Decomposition), EFDD (Enhanced Frequency Domain Decomposition) y CFDD (Curve-Fit Frequency Domain Decomposition)) [17-19]. Éstos se recogen en la tabla 1 para los primeros cuatro modos de vibración:

Los cuatro primeros modos de flexión de la pasarela completa se presentan a continuación.

|

(b) Segundo modo

Figura : Evolución de los errores: (a) En simulación, (b) En experimento.

a) Primer modo b) Segundo modo

Debido a la mayor susceptibilidad del vano 2 de la pasarela, se realizó un EMA en el mismo identificando, entre otros, la masa modal vertical mediante el programa comercial OROS Modal 2. Se realizaron asimismo en este mismo estudio varias estimaciones de dicha masa modal empleando diferentes algoritmos, presentándose en la tabla 2 los resultados obtenidos [5].

Tabla2.

Estimaciones de la masa modal del modo vertical 3.

Posteriormente empleando los registros del EMA realizado en el vano 2 se calcularon las masas modales y los modos normalizados asociados a cada modo de vibración de dicho vano mediante dos técnicas de escalado. Ambas parten de la información obtenida al aplicar el método de identificación de subespacios estocásticos, en inglés Stochastic Subspace Identification (SSI), consistiendo la primera de ellas en la resolución de un problema simétrico de valores propios y la segunda en el escalado de los modos de vibración a partir de la estimación de los residuos de la descomposición en fracciones parciales de la función de respuesta en frecuencia (FRF) de las mediciones [6] y [7]. Aplicando ambas técnicas se calcularon las masas modales y los valores de los modos normalizados respecto a la matriz de masa en las coordenadas L/3, L/2 y 2L/3 del vano 2 (L=51 m).

En la tabla 3 se indican los modos de vibración normalizados y la masa modal del modo vertical 3 de la pasarela estimado mediante el problema simétrico de valores propios expresado en el espacio de estados (método 1) [6].

Tabla 3.

Modos de vibración normalizados y masa modal del modo vertical 3 de la pasarela completa (según método 1).

Tabla 4.

Modos de vibración normalizados y masa modal del modo vertical 3 de la pasarela (según el método 2).

En la tabla 4 se indican los modos de vibración normalizados y la masa modal del modo vertical 3 de la pasarela mediante escalado a partir de la estimación de los residuos de la descomposición en fracciones parciales de la Función de Respuesta en Frecuencia (FRF) de las mediciones (método 2) [7].

3.2. Generación del modelo de la pasarela

La estructura del tramo metálico de la pasarela (vanos 2, 3 y 4) consiste en una celosía tubular tridimensional de sección transversal hexagonal (figura 5), que en el vano 3 se sitúa en el interior de otra celosía formada por cables pretensados longitudinales y transversales. Ambas celosías se enlazan mediante radios tubulares que permiten mantener la forma de los cables y aumentar su excentricidad en el centro del vano. La sección transversal es un hexágono irregular cuya inclinación y dimensiones cambian a lo largo del puente [20].

En cuanto al pretensado transversal, cada sección hexagonal del vano 3 (son 15) consta de un tendón periférico en forma de hexágono, cuyos vértices son los extremos de los seis radios. El pretensado longitudinal consta de seis tendones que pasan por los extremos de los radios de las secciones transversales [21]. Se ha modelizado la estructura de forma que los nudos que conforman tanto el hexágono exterior como el interior, sean coplanarios (figura 6).

La estructura principal está formada por perfiles tubulares de sección circular (figura 7). El diámetro de los perfiles longitudinales es de 273 mm y las diagonales son perfiles de diámetro 135 mm. Los montantes y radios son perfiles tubulares de sección cuadrada de dimensiones 200 mm x 200 mm.

Las secciones de la estructura del tablero (figura 8), constituida por travesaños y arriostramientos, son cuadradas de dimensiones 120 mm x 100 mm y 120 mm x 120 mm respectivamente.

Los apoyos centrales son dos torres constituidas por celosías resueltas también con perfiles tubulares de 355 mm de diámetro (figura 9).

En cuanto a las propiedades del acero se ha considerado un módulo de elasticidad de 2,1·10¹¹ N/m², un coeficiente de Poisson igual a 0,3 y una densidad de 7.700 Kg/m³.

La modelización de la estructura se ha realizado mediante el programa ANSYS empleando dos tipos de elementos: BEAM188 para simular los perfiles tubulares, definiendo para cada tipo las secciones indicadas anteriormente y LINK180 para los tendones de pretensado (figura 10). En el caso del pretensado longitudinal se ha impuesto que los movimientos en la dirección lateral y vertical de los nudos que lo modelizan sean iguales a los de los nudos que representan los extremos de los radios de las secciones transversales. La fuerza de pretensado se ha modelizado mediante el comando INISTATE.

La principal dificultad en la definición geométrica de la pasarela ha sido la modelización del pretensado, principalmente por la elección del tipo de elemento y de las conexiones de éste con el resto de la estructura con el fin de que la transmisión de esfuerzos se realizara de forma realista.

En cuanto a los apoyos, la estructura se considera apoyada en las torres centrales y empotrada en los extremos del tablero. Las bases de las torres se suponen empotradas.

3.3. Ajuste del modelo de la pasarela

Para el ajuste del modelo de la pasarela se emplearon todos los datos experimentales expuestos en los estudios anteriores [5] que afectaban a los cuatro primeros modos verticales de vibración de la pasarela completa, esto es, las frecuencias y las formas modales de dichos modos y la masa modal y los modos normalizados del tercer modo de vibración.

Se empleó un esquema iterativo que nos permitió analizar y seleccionar aquellos parámetros físicos cuyo efecto era más relevante para que las frecuencias y las formas modales de nuestro modelo se aproximaran lo más posible a las obtenidas experimentalmente. En nuestro caso se comprobó que la masa era la variable que mayor influencia tenían en el ajuste una vez modelizada la geometría de la estructura.

En el caso de la masa, ha sido la principal variable que ha afectado al ajuste del tercer modo vertical, no sólo de la frecuencia y la forma modal sino también de su masa modal.

Se ha seleccionado para el análisis modal realizado el método de Block Lanczos y se han obtenido los primeros cuatro modos verticales, que a continuación se reflejan:

a) Primer modo: f1 = 1.224 Hz

Tal y como puede observarse la correspondencia entre las formas modales del modelo en ANSYS (figura 11) y las del OMA (figura 4) es buena.

En la tabla 5 se comparan los resultados del modelo numérico con los obtenidos experimentalmente (tabla 1):

Tabla 5.

Comparativa entre las frecuencias obtenidas experimentalmente y las obtenidas mediante la modelización de la pasarela en ANSYS.

En todos los casos se obtienen diferencias menores del 5% lo que indica una buena correspondencia del modelo.

En el caso de la masa modal del tercer modo se ajustó a un valor de 18057 kg que corresponde aproximadamente al valor máximo que se ha obtenido en las estimaciones experimentales (tablas 2, 3 y 4). La diferencia en porcentaje con respecto al resto de las estimaciones no supera el 6%.

En la tabla 6 se indican los valores resultantes de la comparación entre los modos normalizados respecto a la matriz de masa del modelo en ANSYS y los experimentales obtenidos en [6] y [7]:

Tabla 6.

Valores de MAC comparando los modos normalizados obtenidos numéricamente y los experimentales empleando técnicas de escalado.

Los resultados son cercanos a la unidad lo que nos indica una buena correspondencia entre modos.

4. Análisis dinámico y validación del modelo

Una vez ajustado el modelo de elementos finitos de la pasarela, se obtienen los modos normalizados de los cuatro primeros modos de vibración y, con el resto de los parámetros modales, se procede a simular diferentes fuerzas de excitación para obtener su respuesta, en forma de aceleraciones. El objetivo es poder comparar los resultados con los obtenidos experimentalmente con el objetivo de validar la modelización. Para obtener la respuesta dinámica del modelo se emplearon los valores medios de los coeficientes de amortiguamiento obtenidos experimentalmente en [5] correspondientes a dichos modos. Estos coeficientes se muestran en la tabla 7.

Tabla 7.

Coeficientes de amortiguamiento identificados experimentalmente.

Para validar el comportamiento dinámico observado en el vano 2 del modelo se simularon las fuerzas generadas durante el EMA y se calcularon las aceleraciones en los mismos puntos en los que se obtuvieron mediciones experimentales (L/3, L/2 y 2L/3 del vano 2 siendo L=51). A continuación, se muestran algunos ejemplos de fuerzas simuladas (figura 12) y de las comparativas entre las respuestas experimentales y las obtenidas numéricamente (figuras 13, 14 y 15).

a) Ejemplo 1       b) Ejemplo 2 c) Ejemplo 3

a) Ejemplo 1 b) Ejemplo 2 c) Ejemplo 3

a) Ejemplo 1       b) Ejemplo 2 c) Ejemplo 3

a) Ejemplo 1       b) Ejemplo 2 c) Ejemplo 3

Tabla 8.

Valores de aceleraciones máximas (en m/s²).

Como puede observarse en la tabla 8, donde se muestran las aceleraciones máximas obtenidas tanto numérica como experimentalmente, las respuestas del modelo son ligeramente inferiores salvo a L/3 del vano 2, que son algo superiores debido a que el modo normalizado del modelo en esa coordenada es también superior al obtenido experimentalmente. En el resto de los casos las diferencias son probablemente consecuencia de que el modelo sólo tiene en cuenta la fuerza de excitación aplicada y no el resto de las cargas ambientales que pudieran existir. Las diferencias son más notables y variables en el caso de que la fuerza artificial aplicada es menor y por tanto el peso en la respuesta de las cargas ambientales es mayor, lo que corroboraría la hipótesis anterior. Tal y como queda reflejado en la tabla anterior (tabla 8) las diferencias en el valor máximo de las aceleraciones se encuentran en torno al 25% lo que indica que la correspondencia entre la respuesta del modelo y la experimental es buena.

Para estudiar el comportamiento del modelo completo se ha procedido a simular la principal solicitación a la que están sometidas este tipo de estructuras, esto es, el tránsito peatonal, para comparar la respuesta obtenida con la registrada experimentalmente en varios puntos de la pasarela y así verificar el comportamiento de la estructura completa.

Las aceleraciones que se emplean de referencia se extraen de registros de aceleraciones verticales que corresponden al tránsito peatonal de hasta 10 individuos¹.

La simulación de dicha solicitación sobre la pasarela se realiza empleando el programa Mathematica y se obtienen las respuestas dinámicas en aquellos puntos donde tenemos registros experimentales para poder comparar posteriormente ambos resultados.

A continuación se muestran algunos ejemplos de las comparativas obtenidas entre respuestas dinámicas en forma de aceleraciones del modelo numérico y las registradas experimentalmente (figuras 16, 17 y 18):

Siendo conscientes de las limitaciones del modelo de flujo peatonal y de la variabilidad de los registros experimentales de este tipo de aceleraciones, la finalidad del análisis de estos resultados es comprobar que los obtenidos numéricamente son del mismo orden de magnitud que los reales y que la estructura modelizada completa simula con adecuada precisión este tipo de solicitación. Tal y como puede desprenderse de las figuras anteriores (figuras 16, 17 y 18) los valores máximos alcanzados analíticamente son similares a los experimentales, situándose las mayores diferencias en torno al 20% en la mayoría de los casos.

5. Conclusiones

A la vista de los resultados obtenidos se puede llegar a la conclusión de que el modelo de elementos finitos representa de forma bastante precisa los parámetros modales de la estructura a pesar de la complejidad geométrica de la misma.

En el caso de la respuesta dinámica se obtienen buenos resultados considerando la multitud de factores –y errores inherentes a ellos- que intervienen tanto en la obtención de los resultados numéricos como los experimentales (en especial en el caso del tránsito peatonal). Se ha podido comprobar que si la geometría de partida es realista y completa el ajuste únicamente de las frecuencias y formas modales proporciona una precisión adecuada del modelo.

Una vez validado, este modelo nos permitirá estudiar los diferentes factores que provocan la susceptibilidad a las vibraciones de la pasarela para, por ejemplo, contemplar diferentes alternativas de diseño y/o predecir su comportamiento ante diferentes solicitaciones durante las diversas fases del Proyecto. Asimismo se puede emplear el conocimiento generado en el presente estudio en las fases de diseño de otros casos similares, reduciendo de este modo las inversiones que serían necesarias implantar para solucionar efectos dinámicos no deseados en etapas más tardías que habitualmente son mucho más elevadas. Este análisis puede abrir nuevas líneas de investigación que pueden ser recogidas en futuras publicaciones.

Todo ello contribuirá a una mejor compresión del comportamiento dinámico de este tipo de estructuras favoreciendo asimismo la proliferación del uso de este tipo de herramientas y así ayudar a los proyectistas y a las administraciones responsables de su diseño y mantenimiento.

Agradecimientos

Este trabajo ha sido realizado en parte con la financiación aportada al proyecto de investigación BIA2011-28493-C02-02 por el programa de apoyo a Proyectos de Investigación Fundamental No Orientada (PIFNO) del Ministerio de Ciencia.

Los autores agradecen al Ayuntamiento de Valladolid y al Museo de la Ciencia su apoyo en las pruebas de campo.

Bibliografía