ES2881831B2 - Sistema de sintesis y control de senales aleatorias multivariadas - Google Patents
Sistema de sintesis y control de senales aleatorias multivariadas Download PDFInfo
- Publication number
- ES2881831B2 ES2881831B2 ES202030500A ES202030500A ES2881831B2 ES 2881831 B2 ES2881831 B2 ES 2881831B2 ES 202030500 A ES202030500 A ES 202030500A ES 202030500 A ES202030500 A ES 202030500A ES 2881831 B2 ES2881831 B2 ES 2881831B2
- Authority
- ES
- Spain
- Prior art keywords
- control
- kurtosis
- signals
- matrix
- control system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000015572 biosynthetic process Effects 0.000 title claims description 48
- 238000003786 synthesis reaction Methods 0.000 title claims description 48
- 239000011159 matrix material Substances 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 31
- 239000013598 vector Substances 0.000 claims description 31
- 230000003595 spectral effect Effects 0.000 claims description 29
- 230000009466 transformation Effects 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 description 19
- 230000033001 locomotion Effects 0.000 description 13
- 230000008569 process Effects 0.000 description 13
- 238000009826 distribution Methods 0.000 description 11
- 238000000354 decomposition reaction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 230000005520 electrodynamics Effects 0.000 description 3
- 238000012804 iterative process Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 241000238631 Hexapoda Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/022—Vibration control arrangements, e.g. for generating random vibrations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/06—Multidirectional test stands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4472—Mathematical theories or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D19/00—Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase
- G05D19/02—Control of mechanical oscillations, e.g. of amplitude, of frequency, of phase characterised by the use of electric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Signal Processing (AREA)
- Probability & Statistics with Applications (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Acoustics & Sound (AREA)
- Feedback Control In General (AREA)
Description
DESCRIPCIÓN
SISTEMA DE SÍNTESIS Y CONTROL DE SEÑALES ALEATORIAS
MULTIVARIADAS
Campo de la invención
La presente invención tiene como objetivo principal la síntesis y control de señales aleatorias multivariadas Gaussianas y no Gaussianas y su aplicación en el campo de los ensayos acústicos y de vibraciones para componentes de vehículos automóviles. No obstante, su uso no queda limitado a este sector tecnológico, siendo aprovechable en otras áreas de la técnica en las que la síntesis y control de señales aleatorias multivariadas No Gaussianas sean un requerimiento para la realización de ensayos.
Antecedentes de la invención
La síntesis y control de señales aleatorias se utiliza ampliamente en el sector de automoción con el propósito de reproducir en laboratorio el entorno vibratorio del vehículo. En la actualidad está muy extendido el uso de sistemas univariadas o monoaxiales, que permiten sintetizar y controlar señales estacionarias unidimensionales en una única dirección del espacio. El sistema de síntesis y control genera una señal que es transmitida a un amplificador de potencia que, tras amplificar la señal, la transmite a un actuador electrodinámico o electrohidráulico que la convierte en un movimiento. Este movimiento se registra mediante un acelerómetro dando lugar a una señal de control que se transmite nuevamente al sistema de síntesis y control. Éste ajusta su señal de salida para compensar las distorsiones introducidas por el amplificador y el actuador de modo que la señal de control cumple con una densidad espectral de potencia (PSD) objetivo.
En lo tocante a la distribución estadística de sus amplitudes, las señales generadas por el sistema de síntesis y control pueden ser clasificadas como Gaussianas (patentes US3710082A y US5012428A) o no Gaussianas. En este último caso, se permite ajustar el parámetro de la distribución de amplitudes denominado kurtosis, de forma que se incrementa o reduce la probabilidad de amplitudes altas respecto de un proceso Gaussiano. Para lograr esta variación de la kurtosis, se pueden emplear diferentes métodos, como son la inclusión de ruido de disparo
filtrado (patente US7426426B2), la generación de un espectro auxiliar asociado a la parte no Gaussiana de la señal (patente US20100305886A1) o mediante manipulación de la fase de la señal (patentes US20130231878A1 y JP5421971B2). En el primer y segundo métodos, la alteración de la kurtosis produce una distorsión no deseada de la densidad espectral de potencia, que debe ser corregida. En el caso de manipulación de la fase de la señal, la densidad espectral de potencia no se ve afectada, lo cual representa una ventaja desde el punto de vista de la síntesis y control de la señal, cuyo objetivo es la síntesis de señales con una densidad de potencia determinada así como una kurtosis determinada.
Los sistemas de síntesis y control de señales aleatorias univariadas o monoaxiales tienen como principal desventaja el no representar completamente el entorno real del vehículo, el cual presenta movimientos en varias direcciones del espacio de forma simultánea. Esto implica que un sistema de síntesis y control de señales aleatorias debe ser multivariado o multiaxial para reproducir adecuadamente el entorno del vehículo. A diferencia de los sistemas univariados o monoaxiales, los sistemas multivariados o multiaxiales sintetizan y controlan múltiples señales que se transmiten a uno o varios amplificadores de potencia, los cuales están conectados a un conjunto de actuadores electrodinámicos o electrohidráulicos que transforman esas señales en movimientos en varias direcciones del espacio simultáneamente. Esos movimientos son captados por un conjunto de acelerómetros, dando a lugar a varias señales de control que se transmiten al sistema de síntesis y control, donde se utilizan para ajustar sus correspondientes señales de salida de modo que se compensen las distorsiones introducidas por los amplificadores y actuadores. A diferencia de los sistemas monoaxiales, en los sistemas multiaxiales se persigue ajustar una matriz de densidad espectral de potencia (SDM) en lugar de una densidad espectral de potencia (PSD). El elemento de la fila n y columna m de dicha matriz se corresponde con el espectro cruzado de potencia entre las componentes n y m de la señal multivariada.
En [1], M. Shinozuka presenta un método para la síntesis de señales multivariadas aleatorias, aunque las señales generales tienen necesariamente una distribución Gaussiana. En [2], DO. Smallwood y T.L. Paez proponen un método para la síntesis de señales aleatorias Gaussianas multivariadas que opera en el dominio de la frecuencia. En [3], DO. Smallwood propone una mejora al método anterior que permite la síntesis de señales aleatorias multivariadas y no Gaussianas, basado en
la aplicación de una transformación no lineal sin memoria (ZMNL). Este método presenta la desventaja de distorsionar el espectro de las señales de salida, lo que reduce el rango dinámico de un posible sistema de control basado en este método. En [4], P. Borgnat, P. Abry y P. Flandrin describen un método para la síntesis de señales aleatorias multiaxiales no Gaussianas que se ajustan a una función de densidad de probabilidad conjunta determinada. Se trata de un método iterativo conceptualmente similar a aquellos basados en la aplicación de transformaciones no lineales sin memoria, en el que se parte de una señal aleatoria Gaussiana a la que se aplica una transformación que la convierte en no Gaussiana, con la salvedad de que en este caso el método asegura que la SDM del proceso se corresponde al objetivo, por lo que no existe distorsión a nivel espectral. En [1]-[4] se describen métodos de síntesis de señales aleatorias multivariadas, pero en ningún momento se trata su control en lazo cerrado. En [5] y [6], R. Zheng, H. Chen y X. He proponen un método para la síntesis y control de señales aleatorias multivariadas y no Gaussianas, basadas en la aplicación del método de transformación no lineal sin memoria (ZMNL) de un proceso Gaussiano. El ajuste de la kurtosis del proceso no Gaussiano se logra mediante un proceso iterativo en el que se modifica sucesivamente la fase de la señal de salida mediante esta transformación no lineal. La distorsión resultante se evita forzando el módulo de los coeficientes de Fourier de la matriz de densidad espectral de potencia (SDM) al valor objetivo, aunque el ajuste de la kurtosis requiere necesariamente de un proceso iterativo en el que se modifica la transformación ZMNL.
Referencias:
[1] Masanobu Shinozuka. Simulation of multivariate and multidimensional random processes. The Journal of the Acoustical Society of America, 49(1B):357-368,1971.
[2] DO Smallwood and Thomas L Paez. A frequency domain method for the generation of partially coherent normal stationary time domain signals. Shock and Vibration, 1(1):45-53, 1993.
[3] DO Smallwood. Generation of partially coherent stationary time histories with non-gaussian distributions. Technical report, Sandia National Labs., Albuquerque, NM (United States), 1996.
[4] P. Borgnat, P. Abry, and P. Flandrin. Using surrogates and optimal transport for synthesis of stationary multivariate series with prescribed covariance function and non-gaussian joint-distribution. In Acoustics, Speech and Signal Processing (ICASSP), 2012 IEEE International Conference on, pages 3729-3732, March 2012.
[5] Ronghui Zheng, Huaihai Chen, and Xudong He. Control method for multiinput multi-output non-gaussian random vibration test with cross spectra consideration. Chinese Journal of Aeronautics, 30(6):1895-1906, 2017.
[6] Ronghui Zheng, Huaihai Chen, and Xudong He. Control method for multiple-input multiple-output non-gaussian random vibration test. Packaging Technology and Science, 30(7):331-345, 2017.
Descripción de la invención
La presente invención propone como novedad la síntesis y control de señales aleatorias multivariadas no Gaussianas mediante la manipulación de la fase de las señales sin necesidad de empleo de una transformación no lineal sin memoria ni de un método iterativo de ajuste de dicha transformación. Para ello las componentes de fase del proceso se generan de forma aleatoria y se ajustan variando los parámetros de un generador aleatorio de fase. Las señales sintetizadas pueden utilizarse para generar vibraciones en una máquina de vibraciones multiaxial o pueden ser registradas para su utilización posterior, como señal de entrada en simulaciones numéricas, por ejemplo. Mediante el adecuado control de las señales de mando generadas es posible ajustar la matriz de densidad espectral de potencia (SDM) y las kurtosis de las señales de control medidas sobre la mesa de vibraciones o sobre el espécimen de ensayo. Las señales generadas son de tipo aleatorio, es decir, no es posible predecir de forma determinista el valor de la señal en cada instante de tiempo, sino que ésta va variando en base a un conjunto de distribuciones estadísticas. Estas señales pueden representar diferentes tipos de procesos vibratorios comunes en la industria del transporte y automoción, como son por ejemplo la rodadura por carretera, el desplazamiento por raíles o el vuelo de cohetes entre otros.
Una máquina de vibración multiaxial es un conjunto de actuadores hidráulicos o electrodinámicos unidos a una plataforma móvil (mesa de vibraciones) con capacidad para generar movimiento en varias direcciones del espacio simultáneamente y en un amplio rango de frecuencias. Las direcciones del
movimiento que se genera dependerán de la disposición de los actuadores y del diseño de la plataforma móvil, y pueden ser dos (sistema biaxial), tres (sistema triaxial) o un número arbitrariamente elevado de direcciones (sistema multiaxial). Para una máquina de vibración multiaxial con una configuración de tipo hexápodo o plataforma de Stewart, se pueden alcanzar hasta seis grados de libertad, correspondientes a movimiento longitudinal, lateral, vertical, alabeo, cabeceo y guiñada. Otros sistemas con más grados de libertad son posibles.
Este tipo de sistemas de vibración pueden reproducir múltiples tipos de señales, como son señales sinusoidales, señales temporales registradas en campo o señales de tipo aleatorio, entre otros. La presente invención aborda la síntesis y control de señales aleatorias multivariadas.
Para controlar el movimiento de este tipo de sistemas, es necesario un sistema de síntesis y control de señales multiaxiales. Este sistema tiene por finalidad sintetizar y ajustar un conjunto de señales de mando de modo que las señales de control resultantes posean los valores de kurtosis y una SDM definida por el usuario, compensando para ello las posibles distorsiones causadas por los diferentes elementos del sistema de vibración.
El sistema de síntesis y control de señales aleatorias multivariadas genera tantas señales de mando como grados de libertad tenga el sistema de vibración multiaxial. Las señales de mando son un conjunto de señales aleatorias con una matriz de densidad espectral de potencia (SDM) y unas kurtosis determinadas. Estas señales se envían a los amplificadores de potencia o a un sistema de control secundario específico del sistema de vibración multiaxial, que generan a su vez las señales de excitación utilizadas para provocar el movimiento de los actuadores. El movimiento producido se mide con un conjunto de sensores ubicados sobre la mesa de vibraciones o sobre el espécimen de ensayo montado sobre ella, dando lugar a unas señales de control que se transmiten al sistema de síntesis y control de señales aleatorias multivariadas. Todos los elementos ubicados entre las salidas del sistema de síntesis y control y sus entradas integran lo que comúnmente se denomina como ‘planta’. Las señales de control se utilizan para modificar las señales de mando hasta que se alcanza la consigna de SDM y las kurtosis deseadas, compensando la distorsión introducida por la planta.
La presente invención propone como solución novedosa el ajuste de las kurtosis individuales de las señales de control mediante la variación de los
parámetros de síntesis de un conjunto de generadores de números aleatorios, sin que dicho ajuste afecte a la SDM. El sistema dispone de tantos generadores de números aleatorios como señales de mando. Las secuencias de números aleatorios se utilizan para calcular la fase de cada una de las señales de control a cada frecuencia, de modo que las kurtosis de las señales de control se ajusten al objetivo definido por el usuario. La SDM de las señales de control se ajusta mediante un proceso de ecualización que se calcula considerando la matriz objetivo definida por el usuario y la distorsión introducida por la planta, integrada entre otros por el amplificador de potencia, la máquina de vibración multiaxial, el espécimen de ensayo y los sensores utilizados para medir las señales de control. Este proceso de ecualización de la SDM es independiente del proceso de ajuste de las kurtosis.
La solución propuesta evita la necesidad de aplicar una función de transformación no lineal sin memoria y el uso de un proceso iterativo para ajustar dicha función de transformación, ya que el método de ajuste de la kurtosis no afecta a la SDM del proceso. Esto conlleva una simplificación del sistema y reducción del coste computacional a la hora de sintetizar y controlar las señales.
El sistema de síntesis y control de señales aleatorias multivariadas no está limitado en el número de grados de libertad que puede controlar, aunque es necesario que el número de señales de control sea mayor o igual que el número de grados de libertad a controlar.
Breve descripción de las figuras
A continuación se describe la invención en una serie de figuras que ayudan a comprender la misma y que se relacionan expresamente con una implementación de dicha invención, que se presenta como un ejemplo no limitativo de ésta, y en las que se representa:
Figura 1 muestra un diagrama de bloques general incluyendo la máquina de vibraciones multiaxial y los sensores que miden las señales de control.
Figura 2 muestra un diagrama de bloques del sistema de síntesis y control de señales multivariadas.
Figura 3 muestra un diagrama de bloques del sistema de síntesis de señales aleatorias no-Gaussianas.
Figura 4 muestra un diagrama de bloques del generador de fases aleatorias. Figura 5 muestra un diagrama de bloques del sistema de control de las
kurtosis.
Figura 6 muestra una posible implementación del dispositivo de síntesis y control de señales aleatorias multivariadas.
Descripción de una realización preferida de la invención
La figura 1 muestra el sistema de síntesis y control de señales aleatorias multivariadas (1.a) y su conexión a la planta (1.b). El sistema de síntesis y control envía a través de sus canales de salida tantas señales de mando (1.c) como grados de libertad se deseen controlar en la máquina de vibraciones (1.d). Las señales de mando son señales de voltaje que se transmiten mediante un conjunto de cables a los amplificadores de potencia o a un sistema de control específico de la máquina de vibración multiaxial (1.e), la cual genera unas señales de excitación (1.f) de elevada potencia que se utilizan para mover los actuadores de la máquina de vibración multiaxial (1.d). Los sensores de aceleración, velocidad o desplazamiento (1.g) ubicados sobre la máquina de vibración multiaxial o sobre el espécimen de ensayo miden el movimiento generado en sus diferentes grados de libertad, y transmiten las señales de control (1.h) a través de un conjunto de cables a los canales de entrada del sistema de síntesis y control de señales aleatorias multivariables (1.a). Este sistema analiza las señales de control y corrige en consecuencia las señales de mando sintetizadas (1.c) hasta que las señales de control recibidas se ajusten a la consigna de matriz de densidad espectral de potencia (SDM) y kurtosis, cerrándose de este modo el lazo de control.
En la figura 2 se muestra un diagrama de operación del sistema de síntesis y control de señales aleatorias multivariadas. El usuario define como parámetros de entrada la matriz de densidad espectral de potencia Sk objetivo, para k = 0,1,..., N - 1 siendo N el número de líneas espectrales definidas por el usuario; y el vector de kurtosis k deseadas para los grados de libertad del sistema. Generalmente, para un sistema con P grados de libertad el vector de kurtosis k será un vector de tamaño P, mientras que la matriz de densidad espectral de potencia Sk será una matriz de tamaño P x P por cada línea espectral definida por el usuario y correspondiente al índice frecuencial k. Para cada línea espectral, la matriz de densidad espectral de potencia debe poseer simetría hermítica, de forma que Sk es igual a su traspuesta conjugada, y además Sk debe ser semidefinida positiva para constituir una SDM válida.
Las señales de control r[n ] pasan inicialmente al módulo de transformación inversa ITM (2.a), que convierte el conjunto de Q señales de control r[n] a un máximo de P señales y[n], correspondientes a los diferentes grados de libertad que se desean controlar, de modo que Q > P. Para ello multiplica al conjunto de Q señales de control por una matriz de transformación, definida por el usuario en base a la disposición y orientación de los sensores utilizados para medir el movimiento generado por la máquina de vibración multiaxial [7]. Las señales resultantes y[n ] pasan al estimador de canal (2.b), que utiliza tanto las señales de control transformadas y[n] como las señales de control x[n] enviadas a la planta para estimar la matriz función de transferencia de la planta H k, la cual es una matriz de tamaño Q x P, para cada línea espectral definida por el usuario y correspondiente al índice frecuencial k. Existen múltiples algoritmos disponibles para la realización de dicha estimación, entre los cuales se encuentra el obtenido a partir del método de estimación espectral de Welch. La matriz función de transferencia de la planta estimada se utiliza en el módulo de control espectral (2.c) para modificar la matriz de densidad espectral de potencia y compensar la distorsión introducida por la planta. La matriz SDM ecualizada Sk se obtiene aplicando por cada índice frecuencial k la expresión
~Sk = [Hk]*Sk ([H k]* )H, (1)
donde [Hk]# denota la pseudoinversa de la matriz de función de transferencia de la planta correspondiente al índice frecuencial k, ([ Hk]#)H denota la matriz traspuesta y conjugada de [H k]#, y Sk la SDM objetivo introducida por el usuario para el índice frecuencial k .
Para cada índice frecuencial k = 0,1,...,N-1, a Sk se le aplica una factorización dando lugar a otra matriz Ak de tamaño P x P, la cual se utiliza como matriz objetivo para la síntesis de la señal de control ecualizada. En general, será válida cualquier descomposición del tipo
s k = A kA « . (2)
Las descomposiciones de Cholesky o en autovalores y autovectores (EVD) son ejemplos que cumplen la ecuación (4) y son por lo tanto descomposiciones
válidas. En el primer caso la matriz Ak es directamente el factor de Cholesky, el cual es triangular inferior. En caso de emplear la EVD, la factorización resultante es de la forma dada en la expresión (3):
Sk = WkDkW», (3)
de modo que las columnas de la matriz Wk de tamaño P x P son los autovectores correspondientes a los autovalores recogidos en la matriz diagonal no negativa Dk de tamaño P x P. La matriz Ak está dada en tal caso por
Ak = WkD1k/2 , (4)
siendo Dk/2 una matriz diagonal cuyos elementos están dados por la raíz cuadrada de los de Dk.
Las matrices A k, k = 0,1,...,N-1, se utilizan tanto para la síntesis de las señales de mando en el sintetizador de ruido no Gaussiano (2.f) como para el ajuste de las kurtosis en lazo abierto en el módulo de control de kurtosis (2.e).
El módulo de estimación de las kurtosis (2.d) estima el valor de este estadístico para cada uno de los P grados de libertad de la señal de control transformada y[n], agrupándolos en el vector k . El módulo de control de kurtosis calcula el parámetro de control de la kurtosis X que se utiliza en el sintetizador de ruido no Gaussiano. Este parámetro es un vector con tantos elementos como grados de libertad P tenga la señal de mando. Para su cálculo puede utilizar el valor de kurtosis objetivo k indicado por el usuario y el valor de kurtosis estimado k cuando se opera en lazo cerrado, o el valor de la kurtosis objetivo k y las matrices A k, k = 0, 1 , . , N-1 calculadas a partir de la SDM objetivo y la función de transferencia de la planta cuando sólo se desea realizar un registro de la señal sintetizada. El uso de unos parámetros u otros dependerá de la elección del usuario.
Con el parámetro de control de la kurtosis X y las matrices A k, en el sintetizador de ruido no Gaussiano (2.f) se sintetizan las señales de mando x[n]. En general, se generarán tantas señales de mando como grados de libertad P tenga la máquina de vibración multiaxial. En ciertos casos puede ser necesario transformar las señales de mando antes de su envío a la planta. Esto ocurre, por ejemplo, en el caso en el que varios actuadores trabajan conjuntamente en un mismo grado de
libertad, o la disposición de los actuadores genera acoplamientos entre los diferentes grados de libertad. Estas señales se trasladan al módulo de transformación matricial TM (2.g) que permite transformar el conjunto de señales de mando a partir de una matriz definida por el usuario, que considera la ubicación y orientación de los diferentes actuadores en la máquina de vibración multiaxial. De este módulo salen las señales de mando transformadas d[n], que o bien se transmiten a la planta con el propósito de generar el movimiento; o bien son registradas para su uso posterior.
La figura 3 ilustra el diagrama de bloques del sintetizador de ruido no Gaussiano. La síntesis se basa en la generación de un conjunto de L componentes independientes entre sí. La generación de L componentes asegura que cada bloque de señal sintetizado puede poseer hasta L picos de señal, lo que permite un ajuste adicional de la distribución estadística de la señal sintetizada más allá de su valor de kurtosis. A partir del parámetro de control de las kurtosis X, y para cada índice frecuencial k = 0,1, ...,N - 1 , se generan los L vectores P-dimensionales de fase aleatoria v (í), i = 1 , ... , ! . Para cada índice frecuencial k y cada componente i, el vector de fase aleatoria v (í) posee P elementos (tantos como grados de libertad tenga la señal de mando). El número de componentes L puede ser especificado manualmente por el usuario o puede ser determinado automáticamente por el sistema. En este último caso, el número de componentes L será el mayor número que permita alcanzar el valor de kurtosis deseado.
El vector de fase aleatoria v (í) correspondiente a la componente i y a la línea espectral k es multiplicado por la matriz 4 (í). Para cada línea espectral k, las matrices 4 (í), i = 1,...,L, se obtienen distribuyendo la energía de la SDM entre las L componentes. La distribución de la energía entre las diferentes componentes puede realizarse en base a múltiples criterios, siempre y cuando se cumpla la condición
El reparto de la energía puede ser determinista o realizarse de forma aleatoria. Como ejemplo y de forma no limitativa, una posibilidad es realizar un reparto igualitario entre las L componentes, de modo que para cada i = 1,...,L, se tiene
Cada vector de fase v (í) se multiplica por su matriz correspondiente, obteniendo así los L vectores de coeficientes complejos F(í) según la expresión (7):
Para cada línea espectral k, los L vectores de coeficientes complejos F( í ) , i = 1,..,,L se suman según la expresión (8):
como paso previo al cálculo de la transformada rápida de Fourier inversa (IFFT, 3.c) según la expresión (9):
Los vectores y[n ] generados mediante la expresión (9), y que toman valores complejos en general, tienen por tanto P componentes y constituyen un bloque de N muestras de señal temporal multivariada. Tanto la parte real como la parte imaginaria de dicho bloque poseen una SDM correspondiente a la SDM ecualizada dada por la expresión (1), así como una kurtosis que se ajusta al valor objetivo k . Cada nuevo bloque sintetizado se pasa al módulo de solapamiento y suma (OLA, 3.d) que se ocupa de aplicar esta técnica, la cual es habitual en sistemas de síntesis, para concatenar los sucesivos bloques y generar una señal de longitud arbitrariamente larga. Como resultado, se sintetizan las señales de mando x[n].
El control de las kurtosis se efectúa mediante la adecuada generación de los vectores de fase v (í), en los generadores de fase. El objetivo de estos generadores es reducir la diferencia de fase entre las componentes frecuenciales para que se produzca una interferencia constructiva entre ellas y se produzca un incremento del nivel de pico de la señal. Esto resulta en última instancia en un incremento de la kurtosis. Alterando los parámetros de los diferentes generadores de números
aleatorios, se puede controlar el grado de alineamiento entre las diferentes componentes frecuenciales y, por ende, el nivel de kurtosis. Como ejemplo y de forma no limitativa, la figura 4 muestra una posible implementación de estos generadores de fase. El vector X de control de las kurtosis consta de un total de P elementos, es decir, tantos como grados de libertad se desean controlar. Para cada uno de esos grados de libertad, el parámetro A¿, i = 1,...,P, controla un generador de números aleatorios con distribución exponencial (4.a), que produce un número aleatorio por cada bloque de señal. El valor aleatorio así generado se utiliza como parámetro de varianza de un segundo generador de números aleatorios con distribución normal (4.b), mientras que la media de este segundo generador se mantiene constante e igual a cero. El generador de números aleatorios normal (4.b) genera los números ^ (í), k = 1, siendo N el número de líneas espectrales definidas por el usuario. Un tercer generador de números aleatorios (4.c) genera un valore® por cada bloque, con una distribución uniforme entre los valores 0 y 1, siendo multiplicado por un coeficiente - 2 nk dependiente del término frecuencial k. Este término - 2 nke(V) tiene por efecto un desplazamiento aleatorio del pico de mayor amplitud del bloque de señal generado en el sintetizador. Sin su uso, el pico de mayor amplitud tendería a situarse siempre en las primeras muestras del bloque dando lugar a un proceso periódico, lo cual no es deseable en la práctica. El generador (4.d) genera un número aleatorio ju® por cada bloque, con una distribución uniforme entre 0 y 2n. Las salidas de los tres generadores aleatorios se combinan en (4.e) dando lugar al término de fase = ju® - 2n±e® ^ (í). De este modo, los vectores de fase v (í) se obtienen según la expresión (10):
,(0 = e J ,- 9 „( k Q (10)
y se integran finalmente (4.g) en el vector v k con P valores por cada índice frecuencial k .
Tal y como se evidencia, el ajuste de las kurtosis no requiere de ningún tipo de transformación no lineal, sino que se logra modificando los parámetros de los generadores de números aleatorios con distribución exponencial. Tal y como se demuestra en esta descripción, el ajuste de las kurtosis no afecta a las matrices Ak que definen la característica espectral del proceso. Por lo tanto, cualquier cambio en
las kurtosis no afectará a la SDM del proceso generado, lo que evita la necesidad de un lazo interno iterativo.
El vector X de control de las kurtosis se calcula en el módulo de control de las kurtosis. Este módulo puede operar de dos modos diferentes, a elección del usuario. En modo de lazo cerrado, el sistema ajusta el vector X de control de las kurtosis únicamente en base a los valores de kurtosis objetivo k y los valores de kurtosis k estimados a partir de la señal de control transformada y[n]. Como ejemplo y de forma no limitativa, el valor del vector de control se puede calcular como se muestra en la figura 5: mediante el sumador (5.d), al valor anterior del vector de control, almacenado en el registro (5.e) se le añade un término proporcional a la diferencia entre los vectores de kurtosis objetivo k y estimadas k calculada en el bloque (5.b), tal y como se muestra en la ecuación (11):
donde es un factor de convergencia ajustable.
En el otro modo de operación posible, consistente en la operación en lazo abierto, un calculador de lazo abierto (5.a) se ocupa del cálculo del vector de control de la kurtosis en base a la kurtosis objetivo k y a las matrices Ak dependientes de la SDM objetivo. La relación entre la kurtosis objetivo k , las matrices Ak y el vector de control de la kurtosis viene dado por la expresión (12):
K — Col C ^ C ^ (12)
en la cual el símbolo 1 representa un vector de tamaño P en el que todas sus componentes son igual a 1; los vectores 72 y 73 están dados por la expresión (13):
en la que el superíndice T indica transposición, y el parámetro A¿, i = 1,...,P, es la componente /-ésima del vector X de control de las kurtosis; las matrices C0,C2 ,C3 de tamaño P x P están dadas a través de sus elementos en la fila i y columna s mediante la expresión (14):
donde
siendo Ak(i, s) el elemento situado en la fila i y columna s de la matriz Ak.
La ecuación (12) constituye una relación no lineal que expresa la kurtosis objetivo k en función de los parámetros de control X y las matrices Ak correspondientes a la descomposición de la SDM dada en la expresión (2). De forma general, la relación (12) puede resolverse numéricamente para la obtención de los valores de X, dados k y las matrices Ak, siendo posible para algunos tipos de descomposición la obtención de una solución cerrada. Como ejemplo y de forma no limitativa, las siguientes expresiones muestran una solución cerrada para el caso de una descomposición de Cholesky de la SDM, aunque otras descomposiciones son aplicables, resultando en expresiones diferentes: cada parámetro de control X es la solución positiva de una ecuación de segundo grado, dada en la expresión (18):
atX\ bíXí + = 0, í = 1, ...,P (18) cuyos coeficientes están dados por las expresiones (19), (20) y (21):
ai = C2( í , í ) C3(U ) + Q i - Ki, (19) bt = 3C2(i, 0 2C3( í , 0 5(Qi - Ki), (20) Ci = 6(Qi - Ki), (21)
siendo
y con los términos Cn(i,s ) dados por la expresión (14). Nótese que estos términos, y por tanto los coeficientes a¡, b¡, c¡ en la expresión (18), dependen únicamente de las matrices Ak y, por lo tanto, de la SDM. A partir las ecuaciones (18) es posible calcular el valor de los parámetros de control de la kurtosis para cada grado de libertad i = 1, ...,P según la expresión (23):
que, bajo la condición A i > 0 ya que los parámetros de control han de ser positivos, da lugar a una única solución válida.
Mediante la operación en lazo abierto, es posible el cálculo preciso del vector de control de las kurtosis, proporcionando una solución exacta.
En ambos modos de operación (lazo cerrado o lazo abierto), es posible corregir el valor de las kurtosis objetivo para tener en cuenta el posible efecto de los parámetros utilizados en el módulo de procesado OLA del sintetizador de ruido no Gaussiano.
La figura 6 muestra una posible implementación del dispositivo, que incluye un módulo de salidas (6.a) al que se conectan las entradas a la planta para la transmisión de las señales de mando; un módulo de entradas (6.b) al que se conectan los cables de los sensores utilizadas para la medida de las señales de control; un módulo de comunicaciones digitales (6.c) para el intercambio de datos a un ordenador de control; un botón de encendido y apagado (6.d) y un pulsador de parada de emergencia (6.e) que provoca una parada controlada del sistema de síntesis y control de señales multivariadas.
Referencias:
[7] Norman Fitz-Coy, Vivek Nagabhushan, and Michael T Hale. Benefits and challenges of over-actuated excitation systems. Shock and Vibration, 17(3):285-303, 2010.
Claims (9)
1. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) caracterizado por que comprende:
- un estimador de canal (2.b) que calcula una estimación de la matriz de funciones de transferencia de la planta;
- un módulo de control espectral (2.c) que calcula una matriz espectral A a partir de una matriz de densidades espectrales de potencia (SDM) objetivo y una estimación de la matriz de funciones de transferencia de la planta;
- un estimador de kurtosis (2.d) que calcula un valor estimado de las kurtosis de las señales de control transformadas;
- un módulo de control de kurtosis (2.e) que calcula un vector de control de las kurtosis A;
- un sintetizador de ruido no Gaussiano (2.f) que genera tantas señales de mando como grados de libertad se desean controlar de acuerdo a una matriz espectral y a un vector de control de las kurtosis de los diferentes grados de libertad.
2. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según la reivindicación 1, en la que el ajuste de las kurtosis durante la síntesis de las señales de mando (2.f) se realiza modificando los parámetros de los generadores de fase aleatoria (3.a).
3. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según cualquiera de las anteriores reivindicaciones, en la que la energía de la señal se distribuye entre una o varias componentes independientes (3.b).
4. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según cualquiera de las anteriores reivindicaciones, en la que los bloques de señal se concatenan mediante una técnica de solapamiento y suma para generar secuencias arbitrariamente largas (3.d).
5. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según cualquiera de las anteriores reivindicaciones, que comprende además un módulo de matriz de transformación inversa (2.a) que convierte las señales de control en unas
señales de control transformadas a partir de una matriz definida por el usuario.
6. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según cualquiera de las anteriores reivindicaciones, que comprende además un módulo de matriz de transformación matricial (2.g) que transforma las señales de mando generadas por el sintetizador en otras señales adaptadas a la configuración de la máquina de vibración multiaxial.
7. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según cualquiera de las anteriores reivindicaciones, en la que el módulo de control de kurtosis (2.e) posee dos modos de operación, en lazo abierto y en lazo cerrado.
8. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según la reivindicación 7, en la que en el modo de control en lazo cerrado el vector de control de la kurtosis se calcula a partir de la kurtosis objetivo k y la kurtosis estimada k .
9. - Sistema de síntesis y control de señales aleatorias multivariadas (1.a) según la reivindicación 7, en la que en el modo de control en lazo abierto el vector de control de las kurtosis se calcula a partir de las kurtosis objetivo k y las matrices frecuenciales Ak .
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ES202030500A ES2881831B2 (es) | 2020-05-29 | 2020-05-29 | Sistema de sintesis y control de senales aleatorias multivariadas |
PCT/ES2021/070381 WO2021240040A1 (es) | 2020-05-29 | 2021-05-27 | Sistema de síntesis y control de señales aleatorias multivariadas |
EP21811749.7A EP4160176A4 (en) | 2020-05-29 | 2021-05-27 | SYSTEM FOR SYNTHESIS AND CONTROL OF MULTIVARIATE RANDOM SIGNALS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ES202030500A ES2881831B2 (es) | 2020-05-29 | 2020-05-29 | Sistema de sintesis y control de senales aleatorias multivariadas |
Publications (2)
Publication Number | Publication Date |
---|---|
ES2881831A1 ES2881831A1 (es) | 2021-11-30 |
ES2881831B2 true ES2881831B2 (es) | 2023-01-09 |
Family
ID=78703349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
ES202030500A Active ES2881831B2 (es) | 2020-05-29 | 2020-05-29 | Sistema de sintesis y control de senales aleatorias multivariadas |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP4160176A4 (es) |
ES (1) | ES2881831B2 (es) |
WO (1) | WO2021240040A1 (es) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230093741A1 (en) * | 2021-09-17 | 2023-03-23 | Triad National Security, Llc | Signal processing methods and apparatus |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3710082A (en) * | 1970-03-03 | 1973-01-09 | Time Data Corp | System for digitally controlling a vibration testing environment or apparatus |
JPH065192B2 (ja) * | 1988-11-25 | 1994-01-19 | アイエムブイ株式会社 | 振動制御装置 |
JP4429358B2 (ja) * | 2004-07-02 | 2010-03-10 | ヴァイブレイション・リサーチ・コーポレーション | ランダム振動のスペクトルと尖度を同時に制御するためのシステム及び方法 |
CN102449455A (zh) * | 2009-06-01 | 2012-05-09 | 布鲁尔及凯尔声音及振动测量公司 | 峰度调节振动控制仪的装置和方法 |
JP5421971B2 (ja) * | 2011-10-13 | 2014-02-19 | 地方独立行政法人大阪府立産業技術総合研究所 | 非ガウス特性振動制御装置 |
US8942930B2 (en) * | 2012-03-05 | 2015-01-27 | Crystal Instruments Corporation | System and method for generation and control of mechanical vibration |
US9317044B2 (en) * | 2012-11-07 | 2016-04-19 | Crystal Instruments Corporation | Mechanical vibration system and control method with limited displacement |
JP2018004504A (ja) * | 2016-07-05 | 2018-01-11 | 地方独立行政法人大阪府立産業技術総合研究所 | 多軸振動制御装置 |
JP6518631B2 (ja) * | 2016-08-02 | 2019-05-22 | 地方独立行政法人大阪産業技術研究所 | 非ガウス性振動制御装置 |
CN106546400B (zh) * | 2016-09-30 | 2018-11-09 | 南京航空航天大学 | 一种多输入多输出非高斯随机振动试验系统及试验方法 |
-
2020
- 2020-05-29 ES ES202030500A patent/ES2881831B2/es active Active
-
2021
- 2021-05-27 WO PCT/ES2021/070381 patent/WO2021240040A1/es unknown
- 2021-05-27 EP EP21811749.7A patent/EP4160176A4/en active Pending
Also Published As
Publication number | Publication date |
---|---|
EP4160176A1 (en) | 2023-04-05 |
EP4160176A4 (en) | 2024-06-26 |
ES2881831A1 (es) | 2021-11-30 |
WO2021240040A1 (es) | 2021-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Elliott | Signal processing for active control | |
Yang et al. | Dynamic modeling and control of a 6-DOF micro-vibration simulator | |
ES2881831B2 (es) | Sistema de sintesis y control de senales aleatorias multivariadas | |
Šika et al. | Two-dimensional delayed resonator for entire vibration absorption | |
Halim et al. | Experimental implementation of spatial H/sub/spl infin//control on a piezoelectric-laminate beam | |
Muthalif et al. | Active control of high-frequency vibration: Optimisation using the hybrid modelling method | |
Zilletti et al. | Optimisation of a velocity feedback controller to minimise kinetic energy and maximise power dissipation | |
Camperi et al. | Parametric study on the optimal tuning of an inertial actuator for vibration control of a plate: theory and experiments | |
Temimi et al. | Time-delay effects on controlled seismically excited linear and nonlinear structures | |
Mousavi | Modeling and controlling a semi-active nonlinear single-stage vibration isolator using intelligent inverse model of an MR damper | |
Kim et al. | Wavefront correction of optical beam for large space mirrors using robust control techniques | |
Naseri et al. | A robust active control scheme for automotive engine vibration based on disturbance observer | |
Roberts et al. | Integrated modeling of the Canadian Very Large Optical Telescope | |
Deshmukh et al. | Primary mirror active control system simulation of prototype segmented mirror telescope | |
Scholte et al. | Active vibro-acoustic control of a flexible beam using distributed control | |
Xu et al. | Frequency shaping-based H∞ control for active pneumatic vibration isolation with input voltage saturation | |
Sedghi et al. | Dynamical simulation of E-ELT segmented primary mirror | |
Ungnad et al. | Active control of target sound fields using structural-acoustic brightness applied to a ship model's acoustic signature | |
Hu et al. | Vibration control of piezoelectric flexible structure using robust control methodology | |
Navarra et al. | Design of a NES device via efficient stochastic linearisation | |
Yan et al. | Neural networks-based model predictive control for precision motion control of a micropositioning system | |
Quirós-Pacheco et al. | The Giant Magellan Telescope phasing system: Algorithms and performance simulations | |
Dayou et al. | Control of kinetic energy of a one-dimensional structure using multiple vibration neutralizers | |
Aouragh et al. | Boundary stabilization of an overhead crane with beam model | |
Engstrom et al. | Misregistration in adaptive optics systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
BA2A | Patent application published |
Ref document number: 2881831 Country of ref document: ES Kind code of ref document: A1 Effective date: 20211130 |
|
FG2A | Definitive protection |
Ref document number: 2881831 Country of ref document: ES Kind code of ref document: B2 Effective date: 20230109 |