ES2399624T3 - Procedimiento de evaluación del envejecimiento de un conjunto electrónico - Google Patents
Procedimiento de evaluación del envejecimiento de un conjunto electrónico Download PDFInfo
- Publication number
- ES2399624T3 ES2399624T3 ES09768196T ES09768196T ES2399624T3 ES 2399624 T3 ES2399624 T3 ES 2399624T3 ES 09768196 T ES09768196 T ES 09768196T ES 09768196 T ES09768196 T ES 09768196T ES 2399624 T3 ES2399624 T3 ES 2399624T3
- Authority
- ES
- Spain
- Prior art keywords
- kinematic
- points
- damage
- electronic assembly
- electronic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000032683 aging Effects 0.000 title claims abstract description 12
- 238000011156 evaluation Methods 0.000 title claims description 7
- 230000035882 stress Effects 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 230000005284 excitation Effects 0.000 claims description 30
- 230000001133 acceleration Effects 0.000 claims description 14
- 230000004044 response Effects 0.000 description 14
- 238000005070 sampling Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- 239000000956 alloy Substances 0.000 description 5
- 229910045601 alloy Inorganic materials 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 229910000679 solder Inorganic materials 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000003466 welding Methods 0.000 description 2
- 238000000429 assembly Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- 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/025—Measuring arrangements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0676—Force, weight, load, energy, speed or acceleration
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05K—PRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
- H05K3/00—Apparatus or processes for manufacturing printed circuits
- H05K3/30—Assembling printed circuits with electric components, e.g. with resistor
- H05K3/32—Assembling printed circuits with electric components, e.g. with resistor electrically connecting electric components or wires to printed circuits
- H05K3/34—Assembling printed circuits with electric components, e.g. with resistor electrically connecting electric components or wires to printed circuits by soldering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Procedimiento de evaluación del envejecimiento por fatiga a vibración de un conjunto electrónico, en el que- se miden una o varias magnitudes cinemáticas en uno o varios puntos particulares que experimentan elmovimiento vibratorio, caracterizado porque: - se establece un modelo dinámico que relaciona las magnitudes cinemáticas medidas y unas tensionesmecánicas experimentadas en unos puntos críticos con relación a la fatiga a vibración, denominándose estemodelo observador de estado, - se infiere de las tensiones mecánicas calculadas según el observador de estado un estado de daño decada uno de los puntos críticos, estados de daño que caracterizan el envejecimiento que se persigue, - se integra la cadena medidas cinemáticas - observador de estado - cálculo de los estados de daño en undispositivo embarcado que dispone de una capacidad de cálculo autónoma.
Description
Procedimiento de evaluación del envejecimiento de un conjunto electrónico.
La invención concierne a un procedimiento para evaluar el estado de envejecimiento denominado daño por fatiga a vibración en un sistema electrónico cualquiera. Puede tratarse de placas electrónicas individuales o de un conjunto de placas electrónicas. La evaluación estriba en un dispositivo de medida y de procesado digital de magnitudes físicas relativas a la solicitación que es causante del daño.
Actualmente existen dispositivos similares para piezas y estructuras mecánicas, pero no disponen de determinadas características necesarias para su aplicación en las placas electrónicas. Estas características son la capacidad de procesar, en tiempo real y de manera autónoma con toda la capacidad de cálculo embarcada, la estimación del daño sin recurrir a telemedida ni a post-procesado.
La telemedida y el post-procesado de datos son de uso corriente para supervisar estructuras aéreas y obras de fábrica, pero los correspondientes dispositivos no dan respuesta simultáneamente a los criterios de aptitud para ir embarcados y de autonomía de toda la potencia de cálculo.
La invención se basa en la asociación de uno o varios sensores cinemáticos de la vibración ubicados en puntos particulares, de un observador de estado para evaluar dinámicamente unas tensiones mecánicas localizadas en puntos críticos y de un cálculo de daño en esos puntos críticos, integrables los tres en un dispositivo embarcado y autónomo en cuanto a su función. Los puntos particulares de medida y los puntos críticos de cálculo de daño están diferenciados o coinciden, según las aplicaciones.
La noción de observador de estado halla su definición en las disciplinas de la automática y de la teoría de la información. Se aplica a sistemas dinámicos representados mediante magnitudes de estado. Cuando un estado de un sistema no es accesible a la medida, se construye un modelo dinámico que permite reconstruir el estado no mensurable a partir de medidas de otras magnitudes accesibles. Este modelo dinámico recibe el nombre de observador de estado de la magnitud. Sin embargo, la construcción de un observador de estado no siempre es posible y la posibilidad de construir uno recibe el nombre de condición de observabilidad de esa magnitud.
En la presente invención, con el concurso de sensores, por ejemplo de acelerómetros, son accesibles a la medida magnitudes cinemáticas (una aceleración o una velocidad o una posición o una deformación) de puntos particulares del conjunto electrónico que ha de supervisarse. Por el contrario, en puntos críticos del ensamble electrónico, casi siempre hay tensiones mecánicas inaccesibles a la medida, y más aún estados de daño, salvo en laboratorio. Se puede establecer un observador de estado de las tensiones mecánicas en los puntos críticos que utiliza la medida de una o varias magnitudes cinemáticas en puntos particulares, a condición de observabilidad. El carácter observable de las tensiones mecánicas depende directamente del o de los puntos particulares en los que se ubican el o los sensores cinemáticos. En general, basta con ubicar el o los sensores cinemáticos fuera de los nodos de vibración, esto es, ubicarlo(s) por ejemplo en la fijación del conjunto electrónico.
Cada vez se impone más la necesidad de saber evaluar correctamente el daño a vibración en placas electrónicas embarcadas. Esta necesidad se relaciona con un cierto número de puntos críticos de estas placas electrónicas, cuyo sometimiento a vibraciones es susceptible de acarrear el fallo del conjunto electrónico. Por ejemplo, los puntos críticos que han de supervisarse son juntas de soldadura de aleación de ciertos componentes electrónicos. En el caso de tales placas electrónicas, existen sistemas de supervisión embarcados. Estos pueden estar dotados de acelerómetros o de calibres de deformación para medir aceleraciones o deformaciones experimentadas por esa placa. No obstante, estos sensores casi nunca pueden ser ubicados en los puntos críticos de la placa a causa de las pequeñísimas dimensiones de los componentes y, más aún, de las juntas de soldadura de aleación. Por lo tanto, los sistemas existentes no pueden acceder al dato de las tensiones mecánicas en esos puntos críticos. Se sabe medir un estado de envejecimiento general de la placa, pero no se sabe medir un estado de fatiga local de un componente, al no poder emplazar en él un sensor que dé indicación de las tensiones experimentadas por ese componente.
La técnica anterior conoce, por la solicitud de patente PCT N.° WO2007/102155, un procedimiento y un sistema de simulación de fiabilidad.
La técnica anterior conoce asimismo, por la publicación científica LALL P ET AL «Failure-Envelope Approach to Modeling Shock and Vibration Survivability of Electronic and MEMS Packaging», ELECTRONIC COMPONENTS AND TECHNOLOGY, 2005. ECTC ’05, PROCEEDINGS LAKE BUENA VISTA, FL, USA MAY 31 - JUNE 3, 2005, PISCATAWAY, NJ, USA, IEEE, 31 de mayo de 2005, páginas 480-490, un procedimiento de evaluación de fiabilidad a los choques de componentes electrónicos.
Han sido propuestos procedimientos por el laboratorio CALCE, Center for Advanced Life Cycle Engineering, para establecer el nexo entre medida y tensiones en los puntos críticos. Pero la gran laguna de estos procedimientos es la de no tomar en cuenta los efectos dinámicos de las vibraciones. En efecto, el procedimiento que el CALCE
propone omite la mayoría de los modos propios de vibración y estriba únicamente en cálculos estáticos de curvatura, vibración y estriba únicamente en cálculos estáticos de curvatura.
Por el contrario, la presente invención estriba en observadores de estado que dan razón de los efectos dinámicos vinculados a las fuerzas de inercia y a los modos propios de vibración del conjunto electrónico. El procedimiento que la invención propone hace un uso mucho más amplio de las medidas cinemáticas para evaluar con la mayor fidelidad posible las tensiones mecánicas en los puntos críticos. Se respeta el carácter multiaxial de las tensiones mecánicas y la dinámica del conjunto electrónico. La ventaja de evaluar de manera precisa estas tensiones mecánicas está en poder calcular de manera mucho más precisa el daño sufrido por el conjunto electrónico, en el lugar en el que se ubican el o los componentes críticos. La utilización de tal procedimiento permite mejorar las previsiones de mantenimiento, ya que la indicación de daño que éste suministra sirve para anticipar los potenciales fallos. El planteamiento propuesto permite aplicar los principios de la supervisión de las estructuras mecánicas y el procesado digital de las medidas para conocer el estado de daño en varios lugares, merced a un dispositivo embarcado, con uno o varios sensores.
El dispositivo que integra el procedimiento completo es un pequeño sistema electrónico embarcado de supervisión (PHM por Prognostics and Health Management) provisto de uno o varios sensores cinemáticos (aceleración o velocidad o posición o deformación) y de una unidad de cálculo. El sistema de supervisión está destinado a estimar el daño vinculado a un potencial fallo de las placas electrónicas constitutivas del conjunto electrónico que ha de supervisarse. El o los sensores miden una magnitud cinemática en uno o respectivamente varios puntos del conjunto electrónico que ha de supervisarse. Una de las funciones del sistema de supervisión es la evaluación, a partir de la o de las medidas del o de los sensores, de la tensión mecánica a vibración que es causante de un potencial fallo. Esta función la realiza un observador de estado, que permite reconstruir la dinámica de la tensión mecánica a partir de un modelo dinámico y de las medidas cinemáticas. La función se opera en tiempo real. El o los puntos de medida y los puntos de evaluación de las tensiones no son forzosamente los mismos. La utilización de un observador de estado tan sólo requiere de la unidad una potencia de cálculo bastante limitada y, así, el dispositivo completo se puede estructurar por ejemplo en torno a un mero microcontrolador. Por lo tanto, el propio procedimiento de evaluación de envejecimiento por fatiga es el que permite la integración en un pequeño dispositivo embarcado que dispone de una capacidad de cálculo autónoma.
En la práctica, el observador de estado puede realizarse por ejemplo mediante una operación de convolución de la aceleración medida con una respuesta al impulso en tensión mecánica. De esta manera, se toma en cuenta el conjunto de la banda de frecuencias solicitada, devolviendo así todos los efectos dinámicos. Además, se respeta el carácter multiaxial de la tensión mecánica, toda vez que se utiliza una respuesta al impulso para cada una de las componentes del tensor de tensiones. El observador de estado se podría realizar de otro modo que no sea mediante una operación de convolución: por ejemplo, con un esquema de resolución de ecuación diferencial o incluso utilizando un modelo de derivada no entera (tesis de doctorado de Xavier Moreau, n.º 95BOR10512, 1995).
La estimación del daño estriba en la tensión mecánica dada por el observador de estado. Por ejemplo, la estimación del daño puede realizarse mediante cómputo rainflow (descrito en la norma NF-A03-406, «Fatigue sous sollicitations d’amplitude variable, méthode rainflow de comptage des cycles», noviembre de 1993, ISSN 0335-3931) y mediante aplicación de la regla de Palmgren-Miner. La potencia de cálculo requerida para estas operaciones es muy escasa, lo cual también permite la integración en un dispositivo embarcado que dispone de una capacidad de cálculo autónoma.
La invención tiene pues por objeto un procedimiento de evaluación del envejecimiento de un conjunto electrónico, típicamente una placa electrónica, sometido a un movimiento vibratorio, procedimiento en el cual:
- -
- se miden una o varias magnitudes cinemáticas, por ejemplo la aceleración o la deformación, en uno o varios puntos particulares que experimentan el movimiento vibratorio, por ejemplo en el lugar o en los lugares de su fijación a uno o unos soportes, -se establece un modelo dinámico que relaciona las magnitudes cinemáticas medidas y unas tensiones mecánicas experimentadas en unos puntos críticos con relación a la fatiga a vibración, denominándose este modelo observador de estado, -se infiere de las tensiones mecánicas calculadas según el observador de estado un estado de daño de cada uno de esos puntos críticos, estados de daño que caracterizan el envejecimiento que se persigue, -se integra la cadena medidas cinemáticas - observador de estado - cálculo de los estados de daño en un dispositivo embarcado que dispone de una capacidad de cálculo autónoma.
Se comprenderá mejor la invención con la lectura de la descripción siguiente y con la revisión de las figuras que la acompañan. Estas tan sólo se aportan a título indicativo y sin carácter limitativo alguno de la invención. Las figuras muestran:
figura 1: un esquema sinóptico del procedimiento de la invención,
figura 2: un ejemplo práctico de utilización de la invención,
figura 3: un ejemplo de observador de estado de una tensión mecánica de una junta de soldadura de aleación de un componente electrónico crítico, tomado de un caso real.
La figura 1 muestra cómo un conjunto electrónico, por ejemplo una placa electrónica 7, es sometido a vibraciones mecánicas por un entorno 6. De acuerdo con la invención, el procedimiento de supervisión y, por tanto, de evaluación del envejecimiento, se lleva a la práctica con ayuda de un sistema 3 que incluye un sensor cinemático 1, una función llamada observador de estado 4 que procesa las señales γ producidas por el sensor 1 para devolver una
estimación !ˆ de la tensión mecánica verdadera σ en un punto crítico del conjunto 7 y una función de cálculo de
daño a fatiga 5 acumulado en ese punto crítico del conjunto 7 que toma como variable los resultados !ˆ producidos por el observador de estado 4. El observador de estado 4 es capaz de simular las consecuencias de la excitación del entorno 6 en un lugar seleccionado, mientras que el sensor 1 se halla situado a priori en otro lugar.
Por ejemplo, figura 2, el conjunto 7 es una placa electrónica dotada de componentes útiles 8 a 10, en concreto de componentes críticos en el sentido de que su funcionamiento debe estar garantizado y, por lo tanto, supervisado, y sensibles en el sentido de que unas vibraciones son susceptibles de acarrear un fallo. La placa 7 va fijada en su soporte de utilización mediante unas fijaciones tales como 11 a 12. Unos sensores 1 a 2, en el presente caso acelerómetros, miden los efectos de excitaciones e1 a e2 aplicadas en las fijaciones 11 a 12. De acuerdo con la invención, con el observador de estado 4 somos capaces de evaluar las tensiones mecánicas en los lugares en los que se ubican los componentes 8 a 10 sobre la placa 7.
Para simplificar la explicación, consideraremos que las excitaciones e1 a e2 están orientadas en una misma dirección, perpendicularmente a la placa 7, pero no necesariamente son síncronas. No obstante, cabría la posibilidad de construir un observador de estado aun si no se cumpliera esta obligación, pero las vibraciones según direcciones en el plano de la placa 7 ocasionan deformaciones de amplitudes mucho menores que ocasionan menos daños a los componentes electrónicos. En principio, es el soporte el que vibra y el que transmite sus vibraciones a la placa mediante las fijaciones 11 a 12. En un ejemplo, la vibración perpendicularmente a la placa 7 acarrea mayoritariamente unas tensiones de tracción/compresión, o tensiones normales, sobre unas juntas de soldadura de aleación 0 de un componente electrónico crítico 8 (tipo BGA, Ball Grid Array) que potencialmente conducen a la ruptura por tracción/compresión de esa junta y al fallo de la placa electrónica 7. En la invención, especialmente si el acelerómetro es un acelerómetro de tres dimensiones, se puede calcular los daños resultantes de las tres direcciones simultáneamente.
Para simplificar la explicación, también consideraremos que el conjunto electrónico 7 reacciona de manera lineal a las excitaciones vibratorias. Esto significa que cualquier magnitud física relativa a la vibración se comporta de manera aditiva cuando superponemos varias excitaciones vibratorias aplicadas en lugares diferentes o concurrentes. Esta restricción no es limitativa, ya que, en la práctica, casi todos los conjuntos electrónicos se comportan así.
Un observador de estado de la tensión de tracción/compresión σ en el punto crítico 0 se realiza por ejemplo mediante una operación de convolución entre las medidas de aceleración γ1 y γ2 de los puntos particulares 1 y 2, y unas respuestas al impulso s1 y s2 cuya determinación se especifica más adelante. La operación de convolución que realiza el observador de estado se escribe:
n#1 n#1
!ˆ(t) ∃∀s1(k )∋1(t # k )%∀s2(k )∋ 2(t # k )
k ∃0 k ∃0
siendo n un número de puntos necesarios y suficientes para representar las respuestas al impulso s1 y s2, siendo k un número entero que toma sucesivamente todos los valores de 0 a n−1, siendo τ un paso de tiempo de muestreo
de las medidas cinemáticas, y denotando por !ˆ(t) la tensión normal evaluada en el punto crítico 0 en el instante presente t, por γ4(t−kτ) el k-ésimo último valor medido de γ1 según el muestreo al paso de tiempo τ, por γ2(t−kτ) el késimo último valor medido de γ2 según el muestreo al paso de tiempo τ, por s1(kτ) el k-ésimo valor tabulado de la respuesta al impulso s1 según el muestreo al paso de tiempo τ, por s2(kτ) el k-ésimo valor tabulado de la respuesta al impulso s2 según el muestreo al paso de tiempo τ, y utilizando la notación bien conocida E para designar un sumatorio discreto según una variable discreta y unos extremos indicados como índice y exponente de esta letra. La variable que indica el tiempo presente es muestreada a su vez al paso de tiempo τ y toma por tanto valores discretos múltiplos de τ. Por lo tanto, el observador presentado en el presente caso suministra en cada instante t múltiplo de τ una estimación !ˆ de la tensión mecánica verdadera σ en el punto crítico 0 cuya medida directa es casi siempre
imposible.
La figura 3 muestra un ejemplo de observador de estado en forma de convolución tomado del estudio de una placa electrónica real. En este ejemplo, no hay más que una sola excitación vibratoria e1 que se aplica de manera síncrona sobre todos los puntos de fijación de la placa en dirección perpendicular al plano de la placa. Por lo tanto, esta excitación se caracteriza por un número de grados de libertad igual a 1. Para construir un observador de estado
para cada componente del tensor de tensiones en un punto crítico, es necesario y suficiente un sólo sensor, a condición de observabilidad relacionada con las respectivas posiciones del sensor y del punto crítico y con las características del sensor. En el ejemplo de la figura 3, el sensor cinemático es un acelerómetro que mide la
aceleración γ perpendicularmente al plano de la placa, y la tensión mecánica evaluada !ˆ es la tensión normal en una junta con soldadura de aleación de un componente electrónico de tipo BGA. En este ejemplo, el número n de puntos que representan la respuesta al impulso es 512, el paso de tiempo es τ = 0,5 ms y, por lo tanto, la convolución recae sobre un intervalo de tiempo de 0,256 s. La dinámica de vibración tomada en cuenta de la placa electrónica está en un margen de frecuencias de 0 Hz a 1000 Hz, habida cuenta del paso de tiempo de muestreo τ = 0,5 ms. En este ejemplo tomado de un caso real, el número de operaciones aritméticas por segundo que ha de realizar la unidad de cálculo del dispositivo de supervisión (PHM) para la función de observador de estado del punto crítico 0 es de 1 024 000 multiplicaciones por segundo y 1 024 000 sumas por segundo. Por lo tanto, la capacidad de cálculo requerida para la función de observador de estado es bastante modesta, atendiendo a las arquitecturas electrónicas actuales. Por tanto, el dispositivo prescinde completamente de cualquier post-procesado o telemedida, ya que la técnica de cálculo utilizada proporciona el resultado a costa de un número de cálculos limitado y bastante pequeño.
Cabe destacar que el número de sensores necesarios generalmente debe ser mayor o igual que el número de grados de libertad de excitación. Si el número de sensores es inferior, generalmente no se verifica la condición de observabilidad. Esto clarifica el ejemplo dado en la figura 2 con dos sensores cinemáticos 1 y 2 para dos excitaciones vibratorias independientes e1 y e2.
Tal como se anunció anteriormente, a continuación se da un ejemplo de determinación de las respuestas al impulso s1 y s2 para el ejemplo de la figura 2. En lo que sigue, cada magnitud que lleve una tilde (~) designa la transformada de Fourier de una magnitud temporal. La pulsación, que es la variable de una transformada de Fourier, se denota por ω y es homogénea a la inversa de un tiempo. La unidad de ω es usualmente el radián por segundo.
~
Así, la función ∋ (() es resultado de la transformación de Fourier de la función γ1(t), aceleración suministrada por
~
el sensor 1. La función ∋(() es resultado de la transformación de Fourier de la función γ2(t), aceleración
~
suministrada por el sensor 2. La función e(() es resultado de la transformación de Fourier de la función e1(t),
~
excitación que se aplica en la fijación 11. La función e(() es resultado de la transformación de Fourier de la
función e2(t), excitación que se aplica en la fijación 12. La función !~(() es resultado de la transformación de
Fourier de la función σ0(t), tensión mecánica en el punto crítico 0. Se hace constar que las transformadas de Fourier son magnitudes escalares complejas.
Cuatro funciones H11(ω), H12(ω), H21(ω), H22(ω) de la pulsación ω quedan definidas cada una de ellas como la respuesta de uno de los puntos de medida, en términos de aceleración, a una excitación unitaria de pulsación ω en uno de los puntos de excitación, manteniéndose idénticamente nulas todas las demás excitaciones. Mediante la linealidad anteriormente especificada, estas funciones se escriben en términos generales como la razón de una aceleración a una excitación no nula que la origina, manteniéndose idénticamente nulas las demás excitaciones. Esta razón es independiente de la excitación aplicada a causa de la linealidad especificada antes. De manera explícita:
~
∋1(()
H11(() ∃
~
e1(()
e ∃0
~
∋ 2(()
H21(() ∃
~
e1(()
e ∃0
La notación en índice e2 = 0 significa que la razón es calculada para la excitación e1 aplicada sola.
~
∋1(()
H12(() ∃ ~ e2(()
e ∃0
~
∋ 2(()
H22(() ∃
~
e2(()
e ∃0
La notación en índice e1 = 0 significa que la razón es calculada para la excitación e2 aplicada sola.
Estas cuatro funciones son independientes de las excitaciones aplicadas. Pueden ser determinadas experimentalmente, por ejemplo procediendo sucesivamente a las excitaciones unitarias de que se ha tratado anteriormente, barriendo sucesivamente valores de pulsación ω. El método de los elementos finitos, técnica de cálculo muy usual y acreditada en la industria, la ingeniería civil y las ciencias, ofrece asimismo múltiples medios para determinar esas cuatro funciones, por ejemplo calculando numéricamente el procedimiento experimental que permite determinar esas cuatro funciones.
Análogamente, dos funciones G1(ω) y G2(ω) de la pulsación ω quedan definidas cada una de ellas como la respuesta del punto crítico 0, en términos de tensión mecánica, a una excitación unitaria de pulsación ω en uno de los puntos de excitación, manteniéndose idénticamente nulas todas las demás excitaciones. Mediante linealidad, al igual que anteriormente:
!~(()G (() ∃ 0
1 e ~1(()
e ∃0
!~(()G (() ∃ 0
2 e ~2(()
e ∃0
La notación en índice e1 = 0 o e2 = 0 significa que la razón es calculada para una excitación aplicada sola. Las dos funciones G1 y G2 pueden ser determinadas experimentalmente o mediante el método de los elementos finitos, al igual que para las cuatro funciones anteriores.
~~
Las transformadas de Fourier ∋(() y ∋(() de las dos aceleraciones γ1(t) y γ2(t) medidas por los sensores 1 y 2
se escriben superponiendo linealmente las respuestas a las dos excitaciones e1 y e2:
~
+∋ (() ∃ H (()~ e (()% H (()~ e (()
∗~
)∋ 2(() ∃ H21(()e ~1(()% H22(()e ~2(()
1 111 122
- Este sistema escalar se escribe en forma matricial:
- /0, ∃ 1 − ./0 1 ( ) ( ) ( )~( )~ 21 11 2 1 ( ( (∋ (∋ H H
- ,− . ( ) ( ) 22 12 ( ( H H ,− ./0 1 ( )~( )~ 2 1 ( ( e e
- El sistema matricial inverso se escribe:
- /0, ∃ 1 − ./0 1 ( ) ( ) ( )~( )~ 21 11 2 1 ( ( ( ( H H e e
- ,− . # ( ) ( ) 1 22 12 ( ( H H ,− ./0 1 ( )~( )~ 2 1 (∋ (∋
donde el exponente −1 indica la inversa de la matriz. Esto es, en forma de sistema escalar:
+ H (()∋~(()# H (()∋~(()
~ 221 122
e (() ∃
2 1
H (()H (()# H (()H (()
2 1122 1221
∗
2# H (()∋~(()% H (()∋~(()
~ 211 112
e2(() ∃
) H11(()H22(()# H12(()H21(()
La transformada de Fourier !~0(() de la tensión mecánica σ0 en el punto crítico 0 se escribe superponiendo
linealmente las respuestas a las dos excitaciones e1 y e2:
!~(() ∃ G (()~ e (()% G (()~ e (()
0 11 22
Esta expresión se escribe entonces según las transformadas de Fourier de las aceleraciones:
G (()H (()# G (()H (()
1 22 2 21~
!~0(() ∃∋1(()
H (()H (()# H (()H (()
1122 1221
# G (()H (()% G (()H (()
1 12 2 11~
%∋ 2(()
H (()H (()# H (()H (()
1122 1221
~~
Se definen entonces los dos factores delante de ∋ ((), respectivamente ∋(() en la expresión que antecede
~~
como dos funciones s(() y respectivamente s (():
G (()H (()# G (()H (()
~ 122 221
s (() ∃
1 H (()H (()# H (()H (()
1122 1221
# G (()H (()% G (()H (()
~ 112 211
s (() ∃
2 H (()H (()# H (()H (()
1122 1221
~~
Las respuestas al impulso s1(t) y s2(t) son simplemente las transformadas de Fourier inversas de s(() y s ((),
respectivamente. Se hace constar que s1(t) y s2(t) son magnitudes escalares reales, por su naturaleza física. Hay para esto una razón algebraica que obedece a las propiedades de la transformación de Fourier. s1(t) y s2(t) son las dos funciones que, convolucionadas respectivamente con las medidas γ1(t) y γ2(t), devuelven la tensión mecánica σ0(t) en el punto crítico 0 del conjunto electrónico 7:
! (t) ∃3ts (u)∋ (t # u)du %3ts (u)∋ (t # u)du
0 11 22
La variable u bajo el signo de integral es la variable muda de integración.
Una versión discretizada de esta ecuación integral es la expresión antes dada del observador de estado !ˆ(t). Para obtener tal versión discretizada, basta por ejemplo con sustituir las integrales por sumatorios discretos y con no utilizar más que los valores de s1 y s2 tomados en los puntos múltiplos del paso de tiempo τ:
n#1 n#1
!ˆ(t) ∃∀s (k )∋ (t # k )%∀s (k )∋ (t # k )
11 22 k ∃0 k ∃0
La discretización tan sólo es válida con un paso de tiempo suficientemente pequeño para representar correctamente las integrales y las respuestas s1 y s2, tal y como ilustra la figura 3.
Por supuesto, en todo lo expuesto, el denominador H11(ω)H22(ω)−H12(ω)H21(ω) no debe anularse para ninguno de los valores de pulsación ω en el margen de frecuencias de las vibraciones. Si así fuera, no se verificaría la condición de observabilidad: no sería posible deducir las excitaciones a partir de las medidas cinemáticas y habría que posicionar de otro modo los sensores cinemáticos.
El anterior cálculo no es más que una ilustración de una manera de desarrollar un observador de estado para una tensión mecánica de un punto crítico. Ésta es automáticamente generalizable a un sólo sensor o a un número de sensores superior a dos, siempre a condición de observabilidad, por tanto a condición de que el número de canales de medida sea al menos igual al número de grados de libertad de las excitaciones.
El cálculo mediante elementos finitos estriba en el mallado de la estructura geométrica del conjunto electrónico en elementos poliédricos y en la escritura de las ecuaciones de la mecánica en vibración según el método de Galerkin, también denominado formulación débil. El método de los elementos finitos está profusamente documentado en las revistas de matemáticas aplicadas y de mecánica. El cálculo se puede realizar con uno de los abundantes equipos
5 lógicos del mercado. Está disponible una presentación de esta técnica en un artículo titulado «Durability Modelling of a BGA Component under Random Vibration» y presentado en el congreso EuroSimE en Freiburg-im-Breisgau los días 21-23 de abril de 2008 y disponible en el sitio Internet http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4525047&isnumber=4 525005.
10 El cálculo 5 del daño es realizable, por ejemplo, según el método rainflow, bien conocido (norma NF-A03-406), cuando la tensión mecánica evaluada por el observador de estado es escalar. En este método, se identifican todos los extremos locales de la tensión mecánica y un algoritmo permite identificar ciclos de tensión mecánica, definiéndose cada ciclo mediante uno de los mínimos locales y mediante uno de los máximos locales de la tensión mecánica. Cada ciclo contribuye al daño del punto crítico 0. Un ejemplo de cálculo del daño acumulado consiste en
15 utilizar la hipótesis de Palmgren-Miner, según la cual el daño acumulado es una magnitud escalar que se incrementa con cada ciclo en un valor que tan sólo depende de las características del ciclo en cuestión, independientemente de los ciclos que lo han precedido. Cuando el daño alcanza un valor umbral, se puede considerar que hay una notable probabilidad de que se produzca un fallo del punto crítico 0 y el dispositivo de supervisión está en disposición entonces de generar una alarma. Entonces se pueden tomar medidas preventivas antes de que se produzca el fallo
20 en cuestión.
Cuando la tensión mecánica es más generalmente tensorial, existen varios procedimientos. Uno de ellos está presentado en el citado artículo «Durability Modelling of a BGA Component under Random Vibration». De lo contrario, también existen métodos rainflow generalizados para las magnitudes tensoriales. Finalmente, cabe citar
25 unos criterios multiaxiales de Dang Van y de Sines: véase «Mécanique des matériaux solides», J. Lemaître, J.-L. Chaboche, ISBN 97-1.
Claims (7)
- REIVINDICACIONES1. Procedimiento de evaluación del envejecimiento por fatiga a vibración de un conjunto electrónico, en el que5-se miden una o varias magnitudes cinemáticas en uno o varios puntos particulares que experimentan el movimiento vibratorio,caracterizado porque:10 -se establece un modelo dinámico que relaciona las magnitudes cinemáticas medidas y unas tensiones mecánicas experimentadas en unos puntos críticos con relación a la fatiga a vibración, denominándose este modelo observador de estado, -se infiere de las tensiones mecánicas calculadas según el observador de estado un estado de daño de cada uno de los puntos críticos, estados de daño que caracterizan el envejecimiento que se persigue,15 -se integra la cadena medidas cinemáticas - observador de estado - cálculo de los estados de daño en un dispositivo embarcado que dispone de una capacidad de cálculo autónoma.
- 2. Procedimiento según la reivindicación 1, caracterizado porque20 -se ubica un sensor cinemático en un lugar del conjunto electrónico que no es ni una fijación del conjunto electrónico al soporte ni un punto crítico que ha de supervisarse.
- 3. Procedimiento según la reivindicación 1, caracterizado porque el observador de estado se determina mediante cálculo según el método de los elementos finitos.
- 4. Procedimiento según una de las reivindicaciones 1 a 3, caracterizado porque se ubica el o los sensores en uno o unos lugares que permiten la observabilidad de las tensiones mecánicas en los puntos críticos basándose en el análisis vibratorio del conjunto electrónico por medio del método de los elementos finitos.30 5. Procedimiento según una de las reivindicaciones 1 a 4, caracterizado porque
- -
- el número de puntos de medida es mayor o igual que el número de grados de libertad de las excitaciones.
- 6. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque dicho conjunto electrónico 35 es una placa electrónica.
- 7. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque una o varias magnitudes cinemáticas comprende la aceleración.40 8. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque una o varias magnitudes cinemáticas comprende la deformación.
- 9. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque uno o varios puntosparticulares que experimentan el movimiento vibratorio corresponde al lugar o a los lugares de su fijación a uno o 45 unos soportes.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0857792A FR2938648B1 (fr) | 2008-11-17 | 2008-11-17 | Procede d'evaluation du vieillissement d'un ensemble electronique |
FR0857792 | 2008-11-17 | ||
PCT/FR2009/052189 WO2010055272A1 (fr) | 2008-11-17 | 2009-11-16 | Procédé d'évaluation du vieillissement d'un ensemble électronique |
Publications (1)
Publication Number | Publication Date |
---|---|
ES2399624T3 true ES2399624T3 (es) | 2013-04-02 |
Family
ID=41112494
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
ES09768196T Active ES2399624T3 (es) | 2008-11-17 | 2009-11-16 | Procedimiento de evaluación del envejecimiento de un conjunto electrónico |
Country Status (5)
Country | Link |
---|---|
US (1) | US8650957B2 (es) |
EP (1) | EP2356424B1 (es) |
ES (1) | ES2399624T3 (es) |
FR (1) | FR2938648B1 (es) |
WO (1) | WO2010055272A1 (es) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10691847B2 (en) * | 2017-01-13 | 2020-06-23 | Sap Se | Real-time damage determination of an asset |
US10880690B2 (en) | 2018-09-10 | 2020-12-29 | Bby Solutions, Inc. | Multi-axis temporal motion acceleration detection for device performance monitoring |
CN109991549B (zh) * | 2019-04-24 | 2021-02-09 | 东南大学 | 锂离子电池荷电状态及内阻的联合预测方法 |
CN110927249B (zh) * | 2019-11-28 | 2022-09-27 | 北京机电工程研究所 | 轻量化高强度复合材料应力共振疲劳试验方法及试验*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1738378A1 (ru) * | 1989-07-19 | 1992-06-07 | Научно-Исследовательский Институт Строительных Конструкций Госстроя Ссср | Устройство дл возбуждени колебаний |
US5958556A (en) * | 1996-12-19 | 1999-09-28 | Minnesota Mining And Manufacturing Company | Vibration damped and stiffened circuit articles |
US7260509B1 (en) * | 2001-07-06 | 2007-08-21 | Cingular Wireless Ii, Llc | Method for estimating changes in product life resulting from HALT using quadratic acceleration model |
US6704664B2 (en) * | 2001-12-18 | 2004-03-09 | Visteon Global Technologies, Inc. | Fatigue sensitivity determination procedure |
WO2006086021A2 (en) * | 2004-10-25 | 2006-08-17 | University Of Dayton | Method and system to provide improved accuracies in multi-jointed robots through kinematic robot model parameters determination |
IL174186A0 (en) | 2006-03-08 | 2006-08-01 | Expert Dynamics Ltd | Reliability simulation method and system |
DE102008008609A1 (de) * | 2008-02-12 | 2009-08-13 | Siemens Aktiengesellschaft | Verfahren zur Optimierung der Spannungsverteilung bei Anwendungen akustischer Thermographie |
-
2008
- 2008-11-17 FR FR0857792A patent/FR2938648B1/fr not_active Expired - Fee Related
-
2009
- 2009-11-16 EP EP09768196A patent/EP2356424B1/fr not_active Not-in-force
- 2009-11-16 US US13/129,416 patent/US8650957B2/en not_active Expired - Fee Related
- 2009-11-16 WO PCT/FR2009/052189 patent/WO2010055272A1/fr active Application Filing
- 2009-11-16 ES ES09768196T patent/ES2399624T3/es active Active
Also Published As
Publication number | Publication date |
---|---|
EP2356424A1 (fr) | 2011-08-17 |
WO2010055272A1 (fr) | 2010-05-20 |
US20110271761A1 (en) | 2011-11-10 |
EP2356424B1 (fr) | 2012-11-14 |
US8650957B2 (en) | 2014-02-18 |
FR2938648A1 (fr) | 2010-05-21 |
FR2938648B1 (fr) | 2012-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
ES2399624T3 (es) | Procedimiento de evaluación del envejecimiento de un conjunto electrónico | |
Wang et al. | A frequency and spatial domain decomposition method for operational strain modal analysis and its application | |
Humar et al. | Performance of vibration-based techniques for the identification of structural damage | |
US7680620B2 (en) | Dynamic matrix sensitivity measuring instrument for inertial sensors, and measuring method therefor | |
Chuah et al. | Enabling force sensing during ground locomotion: A bio-inspired, multi-axis, composite force sensor using discrete pressure mapping | |
Wang et al. | Damage identification in beam type structures based on statistical moment using a two step method | |
Ghemari et al. | Simulation and experimental validation of new model for the piezoresistive accelerometer displacement | |
Moradipour et al. | Benchmark studies for bridge health monitoring using an improved modal strain energy method | |
Zhang et al. | Substructure identification for plane frame building structures | |
Lall et al. | Prognostication of accrued damage in board assemblies under thermal and mechanical stresses | |
Gu et al. | Soft manipulator fault detection and identification using anc-based lstm | |
Zhu et al. | Sensitivity analysis of transmissibility functions for structural damage detection | |
US10890499B2 (en) | System and method for predicting strain power spectral densities of light machine structure | |
Fu et al. | Six-axis load head with application to electrical conductor nonlinear dynamics | |
Li et al. | Damage detection for structures under ambient vibration via covariance of covariance matrix and consistent regularization | |
Hummel et al. | Structural assessment and spatial damping identification using low-cost acceleration sensors | |
Giraldo et al. | Control‐oriented system identification using ERA | |
Gelber et al. | Testbed to Characterize MEMS Sensor Degradation During Extended Operation | |
Kunath et al. | Predictive reliability with signal based meta-models | |
RU2799168C1 (ru) | Способ и устройство для испытания стрелкового оружия и/или стендовых испытаний ракетных двигателей | |
Bogdanov | Inertial measurement unit AIST-350T: Results of mechanical investigation tests | |
Rohe et al. | Investigation of the effectiveness of using an experiment to validate experimental substructure models | |
Moron et al. | Automatic system for detection and positioning of impacts in metals based on Arduino | |
Menyhardt et al. | Design of modular vibration testing equipment | |
KR102579510B1 (ko) | 임팩트해머 시험에 의한 콘크리트의 강도 추정 시스템 |