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 PDF

Info

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
Application number
ES09768196T
Other languages
English (en)
Inventor
Marc Grieu
Gregor Massiot
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Airbus Group SAS
Original Assignee
European Aeronautic Defence and Space Company EADS France
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by European Aeronautic Defence and Space Company EADS France filed Critical European Aeronautic Defence and Space Company EADS France
Application granted granted Critical
Publication of ES2399624T3 publication Critical patent/ES2399624T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/025Measuring arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0676Force, weight, load, energy, speed or acceleration
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05KPRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
    • H05K3/00Apparatus or processes for manufacturing printed circuits
    • H05K3/30Assembling printed circuits with electric components, e.g. with resistor
    • H05K3/32Assembling printed circuits with electric components, e.g. with resistor electrically connecting electric components or wires to printed circuits
    • H05K3/34Assembling 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)

  1. REIVINDICACIONES
    1. Procedimiento de evaluación del envejecimiento por fatiga a vibración de un conjunto electrónico, en el que
    5-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. 2. Procedimiento según la reivindicación 1, caracterizado porque
    20 -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. 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. 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.
  5. 6. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque dicho conjunto electrónico 35 es una placa electrónica.
  6. 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.
  7. 9. Procedimiento según una de las anteriores reivindicaciones, caracterizado porque uno o varios puntos
    particulares que experimentan el movimiento vibratorio corresponde al lugar o a los lugares de su fijación a uno o 45 unos soportes.
ES09768196T 2008-11-17 2009-11-16 Procedimiento de evaluación del envejecimiento de un conjunto electrónico Active ES2399624T3 (es)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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) 임팩트해머 시험에 의한 콘크리트의 강도 추정 시스템