Solución De La Ecuación De Estado

Solución de la ecuación de estado

Caso homogéneo

-Métodos de solución

Solución por fracciones parciales

Series

Cayley-Hamilton

Valores propios

Caso no-homogéneo

-Solución caso escalar

-Solución caso vectorial

Solución caso homogéneo

Considere x˙ = Ax ,x ∈ R^n ,A ∈ R^(n x n), asumiendo x(0) conocido

Obtenga la transformada de Laplace L

-Se sabe que sería infinita (de Taylor) F(x)= ∑_(n=0)^∞▒(f^((n) ) (a))/n!(〖x-a)〗^n

e^At se llama matriz de transición de estado, y se calcula obteniendo la anti transformada de Laplace

Ejemplo de caso homogéneo

Fracciones parciales

Tiene como solución

Solución caso no-homogéneo (caso escalar)

Suponga que hay una entrada u(t), la cual representa el caso homogéneo del problema de respuesta en el tiempo. Primero se considera un caso escalar , en donde x ∈ R^1,u ∈ R^1

x˙ = ax + bu, x(0) conocido

Resolviendo la ecuación de la siguiente manera podemos determinar x(t)

Solución caso no-homogéneo (caso matricial)

Teorema de Cayley-Hamilton

Ejemplo Teorema de Cayley-Hamilton

Inversión de matrices

Ejemplo Inversión de matrices

Solución de funciones analíticas de matrices

Ejemplo Solución de funciones analíticas de matrices

Solución por series

Ejemplo Solución de matriz transición por series

Descomposición modal

El objetivo del Análisis modal en la mecánica estructural es determinar las frecuencias naturales y modos de vibrar de un objeto o estructura durante vibración libre. Es común utilizar el Método de los elementos finitos (MEF, o FEM por sus siglas en inglés) para desarrollar el análisis porque, como en otros cálculos usando el MEF, el objeto que se analiza puede tener formas arbitrarias y los resultados de los cálculos son aceptables. Los tipos de ecuaciones que surgen del análisis modal son vistas en Sistemas propios. La interpretación física de los valores propios y vectores propios, los cuales vienen de resolver el sistema, representan las frecuencias y modos de vibrar correspondientes. A veces, los únicos modos deseados son los correspondientes a las menores frecuencias porque pueden ser los modos predominantes en la vibración del objeto.

Sistemas propios FEA

Para los problemas más básicos envolviendo un material linealmente elástico el cual obedece la Ley de elasticidad de Hooke, las ecuaciones de matriz tomando la forma de un sistema de masas de resorte tridimensional dinámico. La ecuación generalizada de movimiento es dada como:

Donde:

es la matriz de masa,

es la 2a derivada de tiempo del desplazamiento,

es el desplazamiento,

es la velocidad,

es la matriz de amortiguación,

es la matriz de rigidez, y

es el vector fuerza.

El problema general, con la amortiguación diferente de cero, es un problema de valor propio cuadratico. Sin embargo, para análisis modal vibracional, la amortiguación es generalmente ignorada, dejando solo el primer y tercer términos en el lado izquierdo:

Esta es la forma general de los sistemas propios encontrados en ingeniería estructural usando el Método de los elementos finitos. Adicionalmente, el movimiento armónico es típicamente asumido para la estructura que si es tomada a ser igual , donde es un valor propio (con unidades de cuadrado de tiempo reciproco, e.g., ), y la ecuación se reduce a:

en contraste, la ecuación para el problema estático es:

La cual es esperada cuando todos los términos teniendo un tiempo derivativo son fijados a cero.

Comparación al algebra lineal

En Álgebra lineal, es más común ver la forma estándar de un sistema propio el cual es expresado como:

Ambas ecuaciones pueden ser vistas como la misma porque si la ecuación general es multiplicada por el inverso de la masa, , tomara la forma de este último. debería notarse que los modos inferiores son deseados, resolviendo el sistema más probablemente envolvería el equivalente de multiplicando por el inverso de la rigidez, , un proceso llamado iteración inversa. Cuando esto es hecho, los valores propios resultantes, , relacionan a aquel del original por:

pero los vectores propios son los mismos.

Métodos de solución

Para problemas elásticos lineales las matrices de rigidez y masa y el sistema en general son definidas positivas. Esa propiedad permite aplicar algunos algoritmos numéricos comúnmente aplicados son usados para hallar una solución. Cuando todas las cualidades del sistema son consideradas:

Solo los pequeños valores propios y vectores propios de los bajos modos son deseados.

Las matrices de masa y rigidez son dispersas y altamente bandeadas.

La matriz del sistema es definida definitiva.

una descripción típica de la solución es primero para tridiagonalizar el sistema usando el algoritmo de Lanczos. Después, usa el algoritmo QR para encontrar los vectores propios y valores propios de este primer sistema tridiagonal. Si la iteración inversa es usada, los nuevos valores propios serán relacionados a los viejos por , mientras los vectores propios del original pueden ser calculados de aquellos de las matrices tridiagonalizadas por:

donde es un vector Ritz aproximadamente igual al vector propio del sistema original, es la matriz de vectores Lanczos, y es el vector propio de la matriz tridiagonal.

Bibliografías:

http://www.lawebdefisica.com/dicc/lagrange/

http://www.fime.uanl.mx/salinas/clase8CM.pdf

http://www.dea.icai.upco.es/ Solución de la ecuación de

...