Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
Compartir
Información tecnológica
versión On-line ISSN 0718-0764
Inf. tecnol. v.16 n.2 La Serena 2005
http://dx.doi.org/10.4067/S0718-07642005000200009
Información Tecnológica-Vol. 16 N°2-2005, págs.: 53-60 ELECTROMECÁNICA Utilización de Onditas (Wavelets) en Sistemas Dinámicos Estructurales. Un Enfoque Didáctico Wavelets in Structural Dynamic Systems. A Didactic Perspective F.O. Cappellari, L.A. Lifschitz, H.H. Brito y L.S. Maglione Resumen Este trabajo tiene como objetivo mostrar la utilización de onditas (wavelets) como herramienta de aprendizaje en sistemas dinámicos estructurales y su contribución a la enseñanza de ingeniería mecánica. Se estableció un procedimiento que, mediante onditas, permite la restitución de las fuerzas de entrada de un sistema dinámico estructural, en función de datos correspondientes a las salidas del mismo, siendo conocidas las funciones de respuesta al impulso del sistema. El caso estudiado tiene vinculación con la determinación del efecto propulsivo de un motor de plasma, aplicable a micro- y nano-satélites. Los resultados obtenidos en la restitución de fuerzas de entrada mostraron la conveniencia de la aplicación de wavelets, no detectándose las limitaciones que suelen presentarse con otros métodos descritos en la literatura. Abstract The objective of the present study is to show the use of wavelets as a learning tool in structural dynamic systems and their contribution to teaching in mechanical engineering. A procedure was established using wavelets for the restitution of input forces to the systems' structural dynamics as a function of data representing outputs of the system, when the impulse response functions are known. The case studied was related to the determination of the propulsive effect of a plasma motor applicable to micro- and nano-satellites. The results obtained in the input force restitution showed the advantages of the use of wavelets; and the limitations usually present in other methods described in the literature were not detected. Keywords: wavelets, dynamic systems, microgravity, mechanical devices, nano-satellites.
INTRODUCCIÓN La utilización de Onditas en diversos campos tales como medicina, geología, biología, ingeniería aeroespacial, compresión de datos e imágenes, etc. (Meyer, 1993; Williams y Amaratunga, 1994), ha tenido un desarrollo muy intenso en los últimos años, tomando en consideración que la teoría de Onditas y sus aplicaciones son de relativamente reciente desarrollo (Daubechies, 1992). La utilización de Onditas en ingeniería estructural ha tenido un importante desarrollo en los últimos tiempos (Douka et al., 2003), principalmente en la identificación de sistemas dinámicos (Robertson et al., 1995), detección de fallas en dispositivos mecánicos y como herramienta alternativa en mantenimiento predictivo (Peng y Chu, 2004). Particularmente en el tema de vibraciones se pueden encontrar aplicaciones de Onditas (Newland, 1993). Entre los temas, en el área de vibraciones, que merecen una atención considerable se encuentra la recuperación de fuerzas de entrada. Es propósito del presente trabajo presentar una aplicación de Onditas (wavelets, ondelettes), como herramienta didáctica, teniendo en cuenta que aparece cada vez con mayor énfasis la necesidad de poner en contacto a los estudiantes de carreras de Ingeniería con temas que en el presente son objeto de activa investigación y demuestran tener importantes y variadas aplicaciones (D´Attellis et al.,1995) considerando, asimismo, que las Onditas han ido cobrando mayor importancia en los últimos años representando una alternativa al análisis de Fourier con ventanas. La metodología utilizada consistió en la de-convolución, para sistemas multi-entrada multi-salida, en el dominio wavelet, a partir de las relaciones entre la matriz de entrada, la matriz de los coeficientes de Markov discretos y la matriz de datos de salida, en dicho dominio. La función ondita aplicada es Daubechies 4. El problema de prueba consistió en un modelo, ampliamente difundido en la enseñanza de la ingeniería mecánica, de un sistema dinámico con dos entradas y dos salidas. Los resultados obtenidos consistieron en la restitución de las entradas, para diversas salidas conocidas, consignándose algunas de ellas en el presente trabajo.
MARCO TEÓRICO Transformada Ondita En tanto que los elementos básicos de la Transformada de Fourier con ventanas son las funciones seno y coseno multiplicadas por una ventana deslizante, el elemento básico de la Transformada Ondita es una función de soporte acotado o de rápido decaimiento, tal que:
lo que significa, en cierto sentido, que la misma función y(t) es oscilante. Esto establece una diferencia respecto a la transformación de Fourier, ya que en ella las funciones oscilante son el seno y el coseno, y la ventana es la función que las acota en el tiempo, en tanto que la función Ondita lleva en si misma las oscilaciones y los intervalos de tiempo en los que actúa. La familia de funciones necesaria para definir la Transformada Ondita de una señal f(t), surge de la función y(t) anterior a través de un cambio de variable que depende de dos parámetros:
donde el parámetro a (parámetro de dilatación) se relaciona con la frecuencia de oscilación, y b (parámetro de corrimiento) con la posición en el tiempo de la Ondita. Considerando señales de energía finita, es decir funciones f(t) (donde t corresponde a la variable tiempo) de cuadrado integrable -en el sentido de Lebesgue-, y al espacio -de Hilbert- que ellas forman (L2(Â)), se define como Transformada Ondita:
Análisis multirresolución Para la implementación computacional de la Transformada Ondita se efectúa una discretización, tomando a=2-j, En consecuencia la familia generada por
Se denomina análisis multirresolución de L2(Â) a una sucesión decreciente de subespacios cerrados
que satisface las siguientes propiedades:
existe Es posible demostrar que cada vector
con La función Llamando Wj al complemento ortogonal de Vj relativo a Vj+1, se puede obtener la siguiente descomposición del espacio
resultando:
Si existe una función
A su vez la componente (proyección) de f(t) en Vj estará dada por:
El cambio del índice j significa cambiar la octava de frecuencias analizadas, y desde el punto de vista de los subespacios construidos, es pasar de la aproximación de la señal en Vj a la que puede obtenerse en Vj+1. La función de Escala La ecuación de dilatación: Dado que
La expresión anterior se denomina ecuación de dilatación, o ecuación de refinamiento. La ecuación Ondita: Vincula la función Escala con la función Ondita. Según (7):
y en consecuencia
La expresión anterior se conoce como la Ecuación Ondita.
Onditas Ortogonales Un ejemplo de Ondita ortogonal, utilizada en este trabajo, es la Ondita Daubechies 4,cuyos coeficientes c(k) son:
Las funciones Escala y Ondita se consignan en las Figuras 1 y 2.
Algoritmo de Mallat Si se considera la proyección de una función f(t) en el subespacio Vj, es decir la aproximación de f(t) al nivel j = 1, la misma será, según la expresión (10):
a su vez, como consecuencia de la expresión (7) y considerando asimismo la expresión (9), resulta:
Haciendo n = l 2k, en la ecuación de dilatación (11), resulta:
y de la ecuación Ondita, (13), resulta:
Siendo las bases involucradas ortonormales:
Lo anterior puede generalizarse entre dos niveles cualesquiera consecutivos:
Las expresiones (22) y (23) implican un cambio de base, de la base Operando en sentido inverso se tiene el proceso de síntesis de la señal, en el cual se pasa de la base ortonormal
METODOLOGÍA Restitución de fuerzas de entrada La salida y(t) de un sistema LTI, causal, bajo la acción de una entrada u(t) dada (con condiciones iniciales nulas), puede ser expresada en forma genérica por la integral de convolución [13]:
siendo h(t) la función respuesta al impulso del sistema (h(t) = 0, para t < 0) Para un valor de t determinado, por ejemplo t = tn, el valor de la salida y(tn) resulta:
Considerando la expansión de una función en serie de Onditas, según (9), la entrada puede ser expresada en la forma:
en la expansión de
Sustituyendo (27) y (28) en (26), resulta:
que puede expresarse como:
a su vez, teniendo en cuenta las condiciones de ortonormalidad, resulta:
la expresión (30) se puede reescribir como:
donde Para todos las datos de la respuesta [y(t0 ), y(t1), y(t2)], la relación entrada-salida puede expresarse como:
donde: y(0) es el valor de la salida para t = 0, etc.; h(0)DWT es el vector que contiene los coeficientes del desarrollo en serie Ondita de De lo anterior, resolviendo, se puede expresar u DWT de la forma:
luego se obtiene u utilizando la Transformada Inversa Ondita, es decir:
Generalizando para sistemas multientrada-multisalida:
Por ejemplo, para el caso de dos entradas, m y n, y dos salidas, p y q:
donde, por ejemplo:
Modelo analizado El modelo físico-matemático adoptado contempla las siguientes hipótesis: H1):Se utiliza, en primera instancia, un modelo discreto del sistema. La estructura formada por la lámina flexible que sostiene al motor de plasma, vg. en ambiente de microgravedad, considerado como un bloque cilíndrico de masa M y un momento de inercia de masa Jc (respecto a un eje normal al plano de movimiento que pasa por el centro de masa), se discretiza como un sistema de dos grados de libertad dinámicos, uno de traslación y otro de rotación. La masa inercial -compuesta por la masa del módulo propulsivo y de las piezas que sirven de interface entre el mismo y la lámina resonante- se modela como posicionada en el extremo libre de la lámina resonante. H2):Se considera un sistema lineal invariante en el tiempo. H3):Se supone que el sistema posee un amortiguamiento proporcional de Rayleigh. Un esquema simplificado del modelo adoptado se muestra en la Figura 3. Las ecuaciones de movimiento pueden deducirse utilizando las ecuaciones dinámicas de Lagrange en la forma:
donde qj son las coordenadas generalizadas, T es la energía cinética del sistema, V es la energía potencial (asociada a las fuerzas conservativas), F es la función de disipación de Rayleigh y Qi son las fuerzas generalizadas, asociadas a las acciones no conservativas que no sean proporcionales a la velocidad. Se adopta q1 = x, y q2 = q, desplazamiento y rotación del extremo superior de la lámina resonante, respectivamente.
RESULTADOS Y DISCUSIÓN A los fines de comprobar la eficacia del procedimiento para la restitución de fuerzas de entrada, consignado en el presente trabajo, en el caso de que una de las entradas corresponda a un tren de pulsos, caso semejante a la excitación a la que estará sometido el banco de ensayos, siendo la otra entrada idénticamente nula (siendo ambas salidas no nulas), se efectuó una simulación a tal efecto a) sobre el grado de libertad traslacional se aplicó una señal consistente en un tren de pulsos, según se muestra en Figura 4, y una entrada nula sobre el grado de libertad rotacional. b) se obtuvieron ambas respuestas (salidas 1 y 2) del sistema, correspondientes a las señales de entrada dadas, mediante Simulink de Matlab. Las salidas correspondientes se muestran en las Figuras 5 y 7. c) se aplicó el procedimiento de restitución en base a Onditas, y con las respuestas de las Figuras 5 y 7, y los parámetros de Markov del sistema que se analiza, se obtuvo la restitución de la entrada según se muestra en la Figura 6.
De la comparación de la entrada aplicada, en primera instancia, y la entrada restituida correspondiente, se observa una adecuada concordancia entre ambas. Cabe señalar que no se encontró en la bibliografía consultada trabajos referidos a restitución de fuerzas de entrada en el dominio wavelet.
CONCLUSIONESEl presente trabajo, desarrollado desde una perspectiva educativa, permitió introducir a los estudiantes de carreras de Ingeniería en áreas de conocimiento que son objeto de activa investigación en el presente, despertando el interés de los alumnos hacia temas de actualidad que tienen variadas aplicaciones. El tema tratado permitió relacionar diversos tópicos y contenidos curriculares. A pesar de la sencillez del modelo de dos grados de libertad adoptado para la restitución de fuerzas de entrada, se logró abordar un problema específico de interés en ingeniería como asimismo introducir a los estudiantes en Onditas, a partir de la integral de convolución que relaciona la entrada del sistema, los parámetros de Markov y la salida. De los resultados obtenidos con el procedimiento expuesto en este trabajo se concluye que el mismo ha demostrado un comportamiento muy satisfactorio en los casos analizados, habiéndose logrado una reconstitución total de la entrada. Asimismo, cabe destacar que mediante la utilización de este procedimiento, no se detectaron las limitaciones que suelen presentarse con los procedimientos en el dominio de la frecuencia (Kammer, 1996).
REFERENCIAS Daubechies I., Ten Lectures on Wavelets, SIAM(1992). D´Attellis C.E., M.T. Anaya, M.I. Cavallaro, y F.F. Villaverde, Introducción a las Onditas- Una presentación para Cursos de Grado de Ingeniería con MatLab. Nueva Librería (1995). Douka, E. S. Loutridis, y A. Trochidis, Crack identification in beams using wavelet analysis, Journal of sound and vibration (2003). Kammer, D., Input force reconstruction using a time domain technique, AIAA-96-1201-CP (1996). Meyer Y., Wavelets, Algorithms & Applications, SIAM(1993). Newland D.E., An introduction to random vibrations, spectral & Wavelet Analysis. Longman Group Limited (1993). Peng, Z.K., F.L. Chu, Application of the wavelet transform in a machine condition monitoring and faults diagnostics: a review with bibliography, Mechanical Systems and Signal Processing (2004). Robertson A., K. Park y K. Alvin, Extraction of impulse response data via wavelet transform for structural system identification, ASME DE: Vol. 84-1 (1995). Williams J. y K. Amaratunga, Introduction to wavelets in engineering, International Journal for Numerical Methods in Engineering: Vol. 37, 2365-2388 (1994).
|