Puede marcar fragmentos interesantes de texto que estarán disponibles a través de un enlace único en su navegador.

Método matricial de la rigidez

Método matricial de la rigidez

El método matricial de la rigidez es un método de cálculo aplicable a estructuras hiperestáticas de barras que se comportan de forma elástica y lineal.

Contenido

Introducción

El método consiste en asignar a la estructura de barras un objeto matemático, llamado matriz de rigidez, que relaciona los desplazamientos de un conjunto de puntos de la estructura, llamados nodos, con las fuerzas exteriores que es necesario aplicar para lograr esos desplazamientos (las componentes de esta matriz son fuerzas generalizadas asociadas a desplazamientos generalizados). La matriz de rigidez relaciona las fuerzas nodales equivalentes y desplazamientos sobre los nodos de la estructura, mediante la siguiente ecuación:

(1) \begin{Bmatrix} F_1 + R_1 \\ F_2 + R_2 \\ ... \\ F_n + R_n \end{Bmatrix}_G = 
\begin{bmatrix}
  k_{11} & k_{12} & \cdots & k_{1n} \\
  k_{21} & k_{22} & \cdots & k_{2n} \\
  \cdots & \cdots & \ddots & \cdots \\
  k_{n1} & k_{n2} & \cdots & k_{nn}
\end{bmatrix}_G
\begin{Bmatrix} \delta_1 \\ \delta_2 \\ ... \\ \delta_n \end{Bmatrix}_G

Donde: F_i \, son las fuerzas nodales equivalentes asociadas a las fuerzas exteriores aplicadas sobre la estructura; R_i \, son las reacciones hiperestáticas inicialmente desconocidas sobre la estructura; \delta_i \, los desplazamientos nodales incógnita de la estructura y n \, el número de grados de libertad de la estructura.

La energía de deformación elástica también puede expresarse en términos de la matriz de rigidez mediante la relación:

E_{def} = \frac{1}{2} \boldsymbol{\delta} \cdot [\mathbf{K}(\boldsymbol{\delta})] = 
\frac{1}{2} \sum_{i,j} k_{ij}\delta_i\delta_j

Del teorema de Maxwell-Betti se deduce que la matriz de rigidez debe ser simétrica y por tanto:

k_{ij} = k_{ji}\,

Fundamento teórico

En general, un sólido deformable real, como cualquier medio continuo es un sistema físico con un número infinito de grados de libertad. Así sucede que en general para describir la deformación de un sólido necesitándose explicitar un campo vectorial de desplazamientos sobre cada uno de sus puntos. Este campo de desplazamientos en general no es reductible a un número finito de parámetros, y por tanto un sólido deformable de forma totalmente general no tiene un número finito de grados de libertad.

Sin embargo, para barras largas elásticas o prismas mecánicos de longitud grande comparada con el área de su sección transversal, el campo de desplazamientos viene dado por la llamada curva elástica cuya deformación siempre es reductible a un conjunto finito de parámetros. En concreto, fijados los desplazamientos y giros de las secciones extremas de una barra elástica, queda completamente determinada su forma. Así, para una estructura formada por barras largas elásticas, fijados los desplazamientos de los nudos, queda completamente determinada la forma deformada de dicha estructura. Esto hace que las estructuras de barras largas puedan ser tratadas muy aproximadamente mediante un número finito de grados de libertad y que puedan ser calculadas resolviendo un número finito de ecuaciones algebráicas. El método matricial proporciona esas ecuaciones en forma de sistema matricial que relaciona los desplazamientos de los extremos de la barras con variables dependientes de las fuerzas exteriores.

Esto contrasta con la situación general de los sólidos elásticos, donde el cálculo de sus tensiones internas y deformaciones involucra la resolución de complejos sistemas de ecuaciones diferenciales en derivadas parciales.

Descripción del método

El método matricial requiere asignar a cada barra elástica de la estructura una matriz de rigidez, llamada matriz de rigidez elemental que dependerá de sus condiciones de enlace extremo (articulación, nudo rígido,...), la forma de la barra (recta, curvada, ...) y las constantes elásticas del material de la barra (módulo de elasticidad longitudinal y módulo de elasticidad transversal). A partir del conjunto de matrices elementales mediante un algoritmo conocido como acoplamiento que tiene en cuenta la conectividad de unas barras con otras se obtiene una matriz de rigidez global, que relaciona los desplazamientos de los nudos con las fuerzas equivalentes sobre los mismos.

Igualmente a partir de las fuerzas aplicadas sobre cada barra se construye el llamado vector de fuerzas nodales equivalentes que dependen de las acciones exteriores sobre la estructura. Junto con estas fuerzas anteriores deben considerarse las posibles reacciones sobre la estructura en sus apoyos o enlaces exteriores (cuyos valores son incógnitas).

Finalmente se construye un sistema lineal de ecuaciones, para los desplazamientos y las incógnitas. El número de reacciones incógnita y desplazamientos incógnita depende del número de nodos: es igual a 3N para problemas bidimensionales, e igual a 6N para un problema tridimensional. Este sistema siempre puede ser dividido en dos subsistemas de ecuaciones desacoplados que cumplen:

  • Subsistema 1. Que agrupa todas las ecuaciones lineales del sistema original que sólo contienen desplazamientos incógnita.
  • Subsistema 2. Que agrupa al resto de ecuaciones, y que una vez resuelto el subsistema 1 y substituido sus valores en el subsistema 2 permite encontrar los valores de las reacciones incógnita.

Una vez resuelto el subsistema 1 que da los desplazamientos, se substituye el valor de estos en el subsistema 2 que es trivial de resolver. Finalmente a partir de las reacciones, fuerzas nodales equivalentes y desplazamientos se encuentran los esfuerzos en los nudos o uniones de las barras a partir de los cuales pueden conocerse los esfuerzos en cualquier punto de la estructura y por tanto sus tensiones máximas, que permiten dimensionar adecuadamente todas las secciones de la estructura.


Matrices de rigidez elementales

Para construir la matriz de rigidez de la estructura es necesario asignar previamente a cada barra individual (elemento) una matriz de rigidez elemental. Esta matriz depende exclusivamente de:

  1. Las condiciones de enlace en sus dos extremos (barra bi-empotrada, barra empotrada-articulada, barra biarticulada).
  2. Las características de la sección transversal de la barra: área, momentos de área (momentos de inercia de la sección) y las características geométricas generales como la longitud de la barra, curvatura, etc.
  3. El número de grados de libertad por nodo, que depende de si se trata de problemas bidimensionales (planos) o tridimensionales.

La matriz elemental relaciona las fuerzas nodales equivalentes a las fuerzas aplicadas sobre la barra con los desplazamientos y giros sufridos por los extremos de la barra (lo cual a su vez determina la deformada de la barra).

Barra recta bidimensional de nudos rígidos

Un nudo donde se unen dos barras se llama rígido o empotrado si el ángulo formado por las dos barras después de la deformación no cambia respecto al ángulo que formaban antes de la deformación. Aún estando imposibilitado para cambiar el ángulo entre barras las dos barras en conjunto, pueden girar respecto al nodo, pero manteniendo el ángulo que forman en su extremo. En la realidad las uniones rígidas soldadas o atornilladas rígidamente se pueden tratar como nudos rígidos. Para barra unida rígidamente en sus dos extremos la matriz de rigidez elemental que representa adecuadamente su comportamiento viene dada por:


\left [ K^{(e)} \right ] =
\begin{bmatrix}
  \frac{EA}{L} & 0 & 0 & -\frac{EA}{L} & 0 & 0 \\
  0 & \frac{12EI}{L^3} & \frac{6EI}{L^2} & 0 & -\frac{12EI}{L^3} & \frac{6EI}{L^2} \\
  0 & \frac{6EI}{L^2} & \frac{4EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{2EI}{L} \\
  -\frac{EA}{L} & 0 & 0 & \frac{EA}{L} & 0 & 0 \\
  0 & -\frac{12EI}{L^3} & -\frac{6EI}{L^2} & 0 & \frac{12EI}{L^3} & -\frac{6EI}{L^2} \\
  0 & \frac{6EI}{L^2} & \frac{2EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{4EI}{L}
\end{bmatrix}

Donde:
L, A, I, \, son las magnitudes geométricas (longitud, área y momento de inercia).
E \, la constante de elasticidad longitudinal (módulo de Young).

Alternativamente la matriz de rigidez de una barra biempotrada recta puede escribirse más abreviadamente, introduciendo la esbeltez mecánica característica:


\left [ K^{(e)} \right ] =
\frac{EI}{L^3}
\begin{bmatrix}
  \lambda_k^2 & 0 & 0 & -\lambda_k^2 & 0 & 0 \\
  0 & 12 & 6L & 0 & -12 & 6L \\
  0 & 6L & 4L^2 & 0 & -6L & 2L^2 \\
  -\lambda_k^2 & 0 & 0 & \lambda_k^2 & 0 & 0 \\
  0 & -12 & -6L & 0 & 12 & -6L \\
  0 & 6L & 2L^2 & 0 & -6L & 4L^2
\end{bmatrix}


Donde: \lambda_k := \sqrt{\frac{AL^2}{I}} = \frac{L}{i_{giro}^2} es la es esbeltez mecánica característica.

Barra recta bidimensional con un nudo articulado y otro rígido

En este caso cuando se imponen giros en el nudo articulado no se transmiten esfuerzos hacia el nudo no articulado. En ese caso la matriz de rigidez, usando la misma notación que en la sección anterior, viene dada por:

 \left [ K^{(e)} \right ] = \frac{EI}{L^3} \begin{bmatrix}
  \lambda_k^2 & 0 & 0 & -\lambda_k^2 & 0 & 0 \\
  0 & 3 & 3L & 0 & -3 & 0 \\
  0 & 3L & 3L^2 & 0 & -3L & 0 \\
  -\lambda_k^2 & 0 & 0 & \lambda_k^2 & 0 & 0 \\
  0 & -3 & -3L & 0 & 3 & 0 \\
  0 & 0 & 0 & 0 & 0 & 0 \end{bmatrix}

Donde se ha supuesto que el nudo articulado es el segundo. Si fuera el primero, habría que permutar los elmentos de la matriz anterior para obtener:

 \left [ K^{(e)} \right ] = \frac{EI}{L^3} \begin{bmatrix}
  \lambda_k^2 & 0 & 0 & -\lambda_k^2 & 0 & 0 \\
  0 & 3 & 0 & 0 & -3 & -3L \\
  0 & 0 & 0 & 0 & 0 & 0 \\
  -\lambda_k^2 & 0 & 0 & \lambda_k^2 & 0 & 0 \\
  0 & -3 & 0 & 0 & 3 & 3L \\
  0 & -3L & 0 & 0 & 3L & 3L^2 \end{bmatrix}

Barra recta bidimensional con dos nudos articulados

Puesto que una barra recta de nudos articulados sólo puede transmitir esfuerzos a lo largo de su eje, la correpondiente matriz de rigidez de esa barra sólo tiene componentes diferentes para los grados de libertad longitudinales. En ese caso la matriz de rigidez, usando la misma notación que en la sección anterior, viene dada por:

 \left [ K^{(e)} \right ] = \frac{EA}{L} \begin{bmatrix}
  +1 & 0 & 0 & -1 & 0 & 0 \\
  0 & 0 & 0 & 0 & 0 & 0 \\
  0 & 0 & 0 & 0 & 0 & 0 \\
  -1 & 0 & 0 & +1 & 0 & 0 \\
  0 & 0 & 0 & 0 & 0 & 0 \\
  0 & 0 & 0 & 0 & 0 & 0   \end{bmatrix}

Arco circular bidimensional de nudos rígidos

Barra recta tridimensional de nudos rígidos

Una barra recta tridimensional tiene 6 grados de libertad por nudo (3 de traslación y 3 de orientación), como la barra tiene dos nudos la matriz de rigidez es una matriz de 12 x 12. Además una barra tridimensional puede transmitir torsiones, y también flexión y esfuerzo cortante en dos direcciones diferentes, esa mayor complejida de comportamiento estructural es lo que hace que una barra tridimensional requiera más grados de libertad y un matriz de rigidez más compleja para describir su comportamiento, esta matriz está compuesta de 3 submatrices:

\left [ K^{(e)} \right ] = \begin{bmatrix} \mathbf{A}_1 & \mathbf{B}^T \\
\mathbf{B} & \mathbf{A}_2 \end{bmatrix}

Donde las submatrices son:

\left [ \mathbf{A}_i \right ] = \begin{bmatrix}
   \frac{EA}{L} & 0 & 0 & 0 & 0 & 0 \\
   0 & \frac{12EI_y}{L^3} & 0 & 0 & 0 & \epsilon_i\frac{6EI_y}{L^2} \\
   0 & 0 & \frac{12EI_z}{L^3} & 0 & -\epsilon_i\frac{6EI_z}{L^2} & 0 \\
   0 & 0 &  0 &\frac{GJ}{L} & 0 & 0 \\
   0 & 0 & -\epsilon_i\frac{6EI_z}{L^2} & 0 &\frac{4EI_y}{L} & 0 \\
   0 & \epsilon_i\frac{6EI_y}{L^2} & 0 & 0 & 0 &\frac{4EI_z}{L} \\
\end{bmatrix}
\left [ \mathbf{B} \right ] = \begin{bmatrix}
   -\frac{EA}{L} & 0 & 0 & 0 & 0 & 0 \\
   0 & -\frac{12EI_y}{L^3} & 0 & 0 & 0 & -\frac{6EI_y}{L^2} \\
   0 & 0 & -\frac{12EI_z}{L^3} & 0 & +\frac{6EI_z}{L^2} & 0 \\
   0 & 0 &  0 &-\frac{GJ}{L} & 0 & 0 \\
   0 & 0 & -\frac{6EI_z}{L^2} & 0 &\frac{2EI_y}{L} & 0 \\
   0 & \frac{6EI_y}{L^2} & 0 & 0 & 0 &\frac{2EI_z}{L} \\
\end{bmatrix}

Y las magntiudes geométricas y mecánicas asociadas a la barra son:

L, A; I_y, I_z; J \, son las magnitudes geométricas: longitud de la barra y su área transversal, momentos de área en las direcciones y y z y módulo de torsión, respectivamente.
E, G \, la el módulo de elasticidad longitudinal y el módulo de elasticidad transversal.
\epsilon_1 =+1, \epsilon_2=-1 \, son signos relativos.

Fuerzas nodales

Para cada barra se define un vector elemental de fuerzas nodales generalizadas, que sea estáticamente equivalente, a las fuerzas aplicadas sobre la barra. El tamaño del vector de fuerzas nodales depende de la dimensionalidad de la barra:

\{\mathbf{F}^{(e)}\}\in\begin{cases}\R^6 & \mbox{bidimensional} \\
\R^{12} & \mbox{tridimensional} \end{cases}


Las componentes de este vector conforman un sistema de fuerzas y momentos de fuerza, tal que la fuerza resultante y el momento resultante de las mismas coinciden con la fuerza y momento del sistema de fuerzas original sobre la barra.

Ejemplo

Ejemplo de carga sobre una viga, P es una carga puntual, y q representa una carga por unidad de longitud.

Para las cargas mostradas en la figura adjunta sobre una barra o viga bidimensional el vector de fuerzas nodales consiste en dos fuerzas verticales (FVd, FVi) aplicadas en cada uno de los dos extremos, dos fuerzas horizontales (FHd, FHi) aplicadas en cada uno de los extremos y dos momentos de fuerza (Md, Mi) aplicados en cada uno de los extremos. Esas seis componentes forman el vector de fuerzas nodales. Es sencillo comprobar que la fuerza y el momento resultantes de estas seis componentes son estáticamente equivalentes al sistema de fuerzas original formado por P y q si se toman los siguientes valores:


\{\mathbf{F}^{(e)}\} =
\begin{Bmatrix}F_{Hd}\\ F_{Vd}\\ M_{d}\\ F_{Hi}\\ F_{Vi}\\ M_{i} \end{Bmatrix} =
\frac{1}{\sqrt{2}}\begin{Bmatrix}\frac{2}{3}P\\ -\frac{20}{27}P \\ -\frac{4}{27}PL \\ \frac{1}{3}P \\ -\frac{7}{27}P\\ +\frac{2}{27}PL \end{Bmatrix} + \begin{Bmatrix}0 \\ +\frac{13}{54}qL \\ -\frac{1}{108}qL^2 \\ 0 \\ -\frac{31}{54}qL \\ \frac{1}{324}qL^2 \end{Bmatrix}

Cálculo de desplazamientos

Una vez encontrada la matriz de rigidez global y el vector de fuerzas nodales global se construye un sistema de ecuaciones como (1). Este sistema tiene la propiedad de que puede descomponerse en dos subsistemas de ecuaciones:

  1. El primero de estos sistemas relaciona únicamente los desplazamientos incógnita con algunas de las componentes del vector de fuerzas nodales global y constituye siempre un sistema compatible determinado
  2. El segundo subsistema contiene también las reacciones incógnita y una vez resuelto el primer subsistema es de resolución trivial.

Resolviendo el primer subsistema compatible determinado, se conocen los desplazamientos incógnita de todos los nudos de la estructura. Insertando la solución del primer subsistema en el segundo resultan las reacciones.

Podemos ilustrar el cálculo de desplazamientos con un ejemplo. Por ejemplo si consideramos la flexión en el plano XY de la viga recta de la sección anterior considerando que se trata de una viga biarticulada unida en sus extremos a dos rótulas fijas tendríamos que el sistema general (1) tendría la forma para este caso particular:


\frac{1}{\sqrt{2}}\begin{Bmatrix}\frac{2}{3}P\\ -\frac{20}{27}P \\ -\frac{4}{27}PL \\ \frac{1}{3}P \\ -\frac{7}{27}P\\ +\frac{2}{27}PL \end{Bmatrix} + \begin{Bmatrix}0 \\ +\frac{13}{54}qL \\ -\frac{1}{108}qL^2 \\ 0 \\ -\frac{31}{54}qL \\ \frac{1}{324}qL^2 \end{Bmatrix} +
\begin{Bmatrix} R_{H1}\\ R_{V1}\\ 0\\ R_{H2}\\ R_{V2}\\ 0 \end{Bmatrix} = \begin{bmatrix}
    \frac{EA}{L} & 0 & 0 & -\frac{EA}{L} & 0 & 0 \\
    0 & \frac{12EI}{L^3} & \frac{6EI}{L^2} & 0 & -\frac{12EI}{L^3} & \frac{6EI}{L^2} \\
    0 & \frac{6EI}{L^2} & \frac{4EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{2EI}{L} \\
    -\frac{EA}{L} & 0 & 0 & \frac{EA}{L} & 0 & 0 \\
    0 & -\frac{12EI}{L^3} & -\frac{6EI}{L^2} & 0 & \frac{12EI}{L^3} & -\frac{6EI}{L^2} \\
    0 & \frac{6EI}{L^2} & \frac{2EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{4EI}{L} \end{bmatrix}
\begin{Bmatrix} 0\\ 0\\ \theta_{1}\\ 0\\ 0\\ \theta_{2} \end{Bmatrix}

Las filas 3 y 6 contienen los giros (desplazamientos) incógnita de los extremos de la viga y tomadas en conjunto conforman el primer subsistema para los desplazamientos. Ignorando los términos nulos y reescrito en forma matricial el subsistema de ecuaciones para los desplazamientos es simplemente:


\frac{1}{\sqrt{2}}\begin{Bmatrix}-\frac{4}{27}PL \\ +\frac{2}{27}PL \end{Bmatrix} +
\begin{Bmatrix} -\frac{1}{108}qL^2 \\ \frac{1}{324}qL^2 \end{Bmatrix}
= \begin{bmatrix}
    \frac{4EI}{L} & \frac{2EI}{L} \\
    \frac{2EI}{L} & \frac{4EI}{L} \end{bmatrix}
\begin{Bmatrix} \theta_{1}\\ \theta_{2} \end{Bmatrix}

Cuya solución nos da el valor del ángulo girado por el extremo derecho e izquierdo de la viga bajo esas cargas:

\theta_1 = -\frac{5}{81\sqrt{2}}\frac{PL^2}{EI}-\frac{7}{1944}\frac{qL^3}{EI} \qquad
\theta_2 = +\frac{4}{81\sqrt{2}}\frac{PL^2}{EI}+\frac{5}{1944}\frac{qL^3}{EI}

Una vez conocidos estos valores e insertados en la matriz las filas 1, 2, 4 y 5 nos proporcionan en valor de las cuatro reacciones hiperestáticas desconocidas previamente.

Cálculo de reacciones

Una vez calculados los desplazamientos resolviendo un sistema de ecuaciones, el cálculo de las reacciones es sencillo. A partir de la ecuación (1) tenemos simplemente:

\begin{Bmatrix} R_1\\ R_2\\ ... \\ R_n \end{Bmatrix}_G = 
\begin{bmatrix}
  k_{11} & k_{12} & \cdots & k_{1n} \\
  k_{21} & k_{22} & \cdots & k_{2n} \\
  \cdots & \cdots & \ddots & \cdots \\
  k_{n1} & k_{n2} & \cdots & k_{nn}
\end{bmatrix}_G
\begin{Bmatrix} \delta_1 \\ \delta_2 \\ ... \\ \delta_n \end{Bmatrix}_G -
\begin{Bmatrix} F_1\\ F_2\\ ... \\ F_n\end{Bmatrix}_G

Tomando el mismo ejemplo que en la última sección el cálculo de reacciones sobre la viga biarticulada con carga P y q sería:


\begin{Bmatrix} R_{H1}\\ R_{V1}\\ 0\\ R_{H2}\\ R_{V2}\\ 0 \end{Bmatrix} = \begin{bmatrix}
    \frac{EA}{L} & 0 & 0 & -\frac{EA}{L} & 0 & 0 \\
    0 & \frac{12EI}{L^3} & \frac{6EI}{L^2} & 0 & -\frac{12EI}{L^3} & \frac{6EI}{L^2} \\
    0 & \frac{6EI}{L^2} & \frac{4EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{2EI}{L} \\
    -\frac{EA}{L} & 0 & 0 & \frac{EA}{L} & 0 & 0 \\
    0 & -\frac{12EI}{L^3} & -\frac{6EI}{L^2} & 0 & \frac{12EI}{L^3} & -\frac{6EI}{L^2} \\
    0 & \frac{6EI}{L^2} & \frac{2EI}{L} & 0 & -\frac{6EI}{L^2} & \frac{4EI}{L} \end{bmatrix}
\begin{Bmatrix} 0\\ 0\\ \theta_{1}\\ 0\\ 0\\ \theta_{2} \end{Bmatrix} -
\frac{1}{\sqrt{2}}\begin{Bmatrix}\frac{2}{3}P\\ -\frac{20}{27}P \\ -\frac{4}{27}PL \\ \frac{1}{3}P \\ -\frac{7}{27}P\\ +\frac{2}{27}PL \end{Bmatrix} - \begin{Bmatrix}0 \\ +\frac{13}{54}qL \\ -\frac{1}{108}qL^2 \\ 0 \\ -\frac{31}{54}qL \\ \frac{1}{324}qL^2 \end{Bmatrix}

Introduciendo los valores de los giros en los extremos y multiplicando la matriz de rigidez por el vector de desplazamientos se tiene finalmente que:


\begin{Bmatrix} R_{H1}\\ R_{V1}\\ R_{H2}\\ R_{V2} \end{Bmatrix} = 
\begin{Bmatrix} 0\\ \frac{6EI}{L^2}(\theta_1+\theta_2)\\ 0\\ \frac{6EI}{L^2}(\theta_1+\theta_2) \end{Bmatrix} -
\frac{1}{\sqrt{2}}\begin{Bmatrix}\frac{2}{3}P\\ -\frac{20}{27}P\\ \frac{1}{3}P \\ -\frac{7}{27}P \end{Bmatrix} - \begin{Bmatrix}0 \\ +\frac{13}{54}qL \\ 0 \\ -\frac{31}{54}qL \end{Bmatrix} = 
\begin{Bmatrix}\frac{2}{3\sqrt{2}}P \\ \frac{18}{27\sqrt{2}}P - \frac{20}{81}qL \\
 \frac{1}{3\sqrt{2}}P\\ \frac{5}{27\sqrt{2}}P + \frac{46}{81}qL \end{Bmatrix}

Esto completa el cálculo de reacciones.

Cálculo de esfuerzos

El cálculo de esfuerzos se realiza examinando en coordenadas locales de las barras el esfuerzo axial, los esfuerzos cortantes, los momentos flectores y el momento torsor generados en cada una de las barras, conocidos los desplazamientos de todos los nudos de la estructura. Esto puede hacerse usando las matrices de rigidez expresadas en coordenadas locales y los desplazamientos nodales expresados también en coordenadas locales.

Análisis dinámico

Artículo principal: Análisis dinámico

El análisis estático discutido anteriormente puede generalizarse para encontrar la respuesta dinámica de una estructura. Para ello se require representar el comportamiento inercial de la estructura mediante una matriz de masa \scriptstyle \mathbf{M}, modelizar las fuerzas disipativas mediante una matriz de amortiguamiento \scriptstyle \mathbf{C}, que junto con la matriz de rigidez \scriptstyle \mathbf{K} permiten plantear un sistema de ecuaciones de segundo orden del tipo:

\mathbf{M\ddot{\delta}}(t) + \mathbf{C\dot{\delta}}(t) + \mathbf{K{\delta}}(t) = \mathbf{f}(t)

La solución del sistema anterior pasa por un cálculo de las frecuencias propias y los modos propios. Admitiendo que las fuerzas disipativas son poco importantes las frecuencias propias se pueden determinar resolviendo la siguiente ecuación polinómica en \scriptstyle \omega^2:

\det (\mathbf{K} - \omega^2\mathbf{M}) = 0

Esas magnitudes permiten realizar un análisis modal que reproduce el comportamiento de la estructura bajo diferentes tipos de situaciones.

Referencia

Bibliografía

  • Timoshenko, Stephen; Godier J.N.. McGraw-Hill. ed. Theory of elasticity. 
  • Ortiz Berrocal, Luis. McGraw-Hill. ed. Resistencia de Materiales. Aravaca (Madrid). ISBN 84-7651-512-3. 
  • Monleón Cremades, S., Análisis de vigas, arcos, placas y láminas, Ed. UPV, 1999, ISBN 84-7721-769-6.

Enlaces externos

Programas

  • EBEs (Estructuras de Barras Espaciales - Diagramas de elástica y solicitaciones -Gratuito)
  • W-Trite (programa gratuito de análisis de estructuras de barras)
  • Cype
  • Sap2000
  • ED-Tridim

Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Rigidez — Saltar a navegación, búsqueda En ingeniería, la rigidez es la capacidad de un objeto sólido o elemento estructural para soportar esfuerzos sin adquirir grandes deformaciones o desplazamientos. Los coeficientes de rigidez son magnitudes físicas… …   Wikipedia Español

  • Método de flexibilidad — En ingeniería estructural, el Método de flexibilidad es el clásico método consistente en deformación para calcular fuerzas en miembros y desplazamientos en sistemas estructurales. Su versión moderna formulada en términos de la matriz de… …   Wikipedia Español

  • Método de los elementos finitos en la mecánica estructural — Este artículo o sección necesita referencias que aparezcan en una publicación acreditada, como revistas especializadas, monografías, prensa diaria o páginas de Internet fidedignas. Puedes añadirlas así o avisar …   Wikipedia Español

  • Método de matriz — El método de matriz es un método de análisis estructural usado como un principio fundamental en muchas aplicaciones en la ingeniería civil. El metodo es llevado a cabo mediante el uso de cualquier matriz de rigidez o de flexibilidad. Véase… …   Wikipedia Español

  • Método de los elementos finitos — Solución de MEF en 2D para una configuración de un magnetostato, (las líneas muestran la dirección de la densidad de flujo calculada, y el color, su magnitud) …   Wikipedia Español

  • Método de distribución de momentos — Este artículo o sección necesita una revisión de ortografía y gramática. Puedes colaborar editándolo (lee aquí sugerencias para mejorar tu ortografía). Cuando se haya corregido, borra este aviso por favor. NOTA: este artículo está siendo… …   Wikipedia Español

  • Matriz de rigidez — La matriz de rigidez aparece en varios contextos: En el método matricial de la rigidez, la matriz de rigidez conceptualmente relaciona los desplazamientos de una serie de puntos o nodos, con los esfuerzos puntuales efectivos en dichos puntos.… …   Wikipedia Español

  • Resistencia de materiales — Saltar a navegación, búsqueda Para otros usos de este término, véase Resistencia. La resistencia de materiales clásica es una disciplina de la ingeniería mecánica y la ingeniería estructural que estudia los sólidos deformables mediante modelos… …   Wikipedia Español

  • Momentos de empotramiento — Los momentos de empotramiento son momentos de reacción sobre una viga cuyos extremos están fijos al ser coaccionados para no moverse. Varios métodos de análisis estructural incluyendo el método de distribución de momentos, el método de las… …   Wikipedia Español

  • Celosía (ingeniería) — Saltar a navegación, búsqueda Puente a base celosías planas en sus caras construido para un antiguo ferrocarril (ahora convertido en puente peatonal) En ingeniería estructural, una celosía es una estructura reticular de barras rectas… …   Wikipedia Español