stringtranslate.com

Tensor de tensión de Cauchy

En mecánica continua , el tensor de tensiones de Cauchy (símbolo , llamado así por Augustin-Louis Cauchy ), también llamado tensor de tensiones verdadero [1] o simplemente tensor de tensiones , define completamente el estado de tensiones en un punto dentro de un material en estado deformado , colocación o configuración. El tensor de segundo orden consta de nueve componentes y relaciona un vector de dirección de longitud unitaria e con el vector de tracción T ( e ) a través de una superficie imaginaria perpendicular a e :

[a]

Las unidades básicas del SI tanto del tensor de tensión como del vector de tracción son newton por metro cuadrado (N/m 2 ) o pascal (Pa), correspondiente al escalar de tensión. El vector unitario no tiene dimensiones .

El tensor de tensión de Cauchy obedece a la ley de transformación del tensor ante un cambio en el sistema de coordenadas. Una representación gráfica de esta ley de transformación es el círculo de tensión de Mohr .

El tensor de tensiones de Cauchy se utiliza para el análisis de tensiones de cuerpos materiales que experimentan pequeñas deformaciones : es un concepto central en la teoría lineal de la elasticidad . Para deformaciones grandes, también llamadas deformaciones finitas , se requieren otras medidas de tensión, como el tensor de tensiones de Piola-Kirchhoff , el tensor de tensiones de Biot y el tensor de tensiones de Kirchhoff .

Según el principio de conservación del momento lineal , si el cuerpo continuo está en equilibrio estático, se puede demostrar que los componentes del tensor de tensión de Cauchy en cada punto material del cuerpo satisfacen las ecuaciones de equilibrio ( ecuaciones de movimiento de Cauchy para aceleración cero). . Al mismo tiempo, según el principio de conservación del momento angular , el equilibrio requiere que la sumatoria de momentos respecto de un punto arbitrario sea cero, lo que lleva a la conclusión de que el tensor de tensiones es simétrico , teniendo así sólo seis componentes de tensión independientes. , en lugar de los nueve originales. Sin embargo, en presencia de tensiones de par, es decir, momentos por unidad de volumen, el tensor de tensiones no es simétrico. Este también es el caso cuando el número de Knudsen es cercano a uno, o el continuo es un fluido no newtoniano, lo que puede conducir a fluidos rotacionalmente no invariantes, como los polímeros .

Hay ciertas invariantes asociadas con el tensor de tensión, cuyos valores no dependen del sistema de coordenadas elegido ni del elemento de área sobre el que opera el tensor de tensión. Estos son los tres valores propios del tensor de tensiones, que se denominan tensiones principales.

Principio de tensión de Euler-Cauchy - vector de tensión

Figura 2.1a Distribución interna de fuerzas de contacto y tensiones de par en un diferencial de la superficie interna en un continuo, como resultado de la interacción entre las dos porciones del continuo separadas por la superficie
Figura 2.1b Distribución interna de fuerzas de contacto y esfuerzos de par en un diferencial de la superficie interna en un continuo, como resultado de la interacción entre las dos porciones del continuo separadas por la superficie
Figura 2.1c Vector de tensión en una superficie interna S con vector normal n. Dependiendo de la orientación del plano considerado, el vector de tensión puede no ser necesariamente perpendicular a ese plano, es decir , paralelo a , y puede resolverse en dos componentes: una componente normal al plano, llamada tensión normal , y otra componente paralela a este plano, llamado esfuerzo cortante .

El principio de tensión de Euler-Cauchy establece que sobre cualquier superficie (real o imaginaria) que divide el cuerpo, la acción de una parte del cuerpo sobre la otra es equivalente (equipolent) al sistema de fuerzas y pares distribuidos en la superficie que divide el cuerpo. cuerpo , [2] y está representado por un campo , llamado vector de tracción , definido en la superficie y que se supone depende continuamente del vector unitario de la superficie . [3] [4] : páginas 66–96 

Para formular el principio de tensión de Euler-Cauchy, considere una superficie imaginaria que pasa a través de un punto material interno que divide el cuerpo continuo en dos segmentos, como se ve en la Figura 2.1a o 2.1b (se puede usar el diagrama del plano de corte o el diagrama con el volumen arbitrario dentro del continuo encerrado por la superficie ).

Siguiendo la dinámica clásica de Newton y Euler , el movimiento de un cuerpo material se produce por la acción de fuerzas aplicadas externamente que se supone son de dos tipos: fuerzas superficiales y fuerzas corporales . [5] Así, la fuerza total aplicada a un cuerpo o a una porción del cuerpo se puede expresar como:

En este artículo solo se analizarán las fuerzas superficiales, ya que son relevantes para el tensor de tensión de Cauchy.

Cuando el cuerpo se somete a fuerzas superficiales externas o fuerzas de contacto , siguiendo las ecuaciones de movimiento de Euler , las fuerzas y momentos de contacto internos se transmiten de un punto a otro del cuerpo, y de un segmento a otro a través de la superficie divisoria , debido a la acción mecánica. contacto de una porción del continuo con la otra (Figura 2.1a y 2.1b). En un elemento de área que contiene , con vector normal , la distribución de fuerzas es equipotente a una fuerza de contacto ejercida en el punto P y al momento superficial . En particular, la fuerza de contacto está dada por

¿ Dónde está la tracción superficial media ?

El principio de tensión de Cauchy afirma [6] : p.47-102  que cuando se vuelve muy pequeño y tiende a cero, la relación se vuelve y el vector de tensión de pareja desaparece. En campos específicos de la mecánica del continuo se supone que la tensión de pareja no desaparece; sin embargo, las ramas clásicas de la mecánica continua abordan materiales no polares que no consideran tensiones de par ni momentos corporales.

El vector resultante se define como la tracción superficial , [7] también llamado vector de tensión , [8] tracción , [4] o vector de tracción . [6] dado por en el punto asociado a un plano con un vector normal :

Esta ecuación significa que el vector tensión depende de su ubicación en el cuerpo y de la orientación del plano sobre el que actúa.

Esto implica que la acción de equilibrio de las fuerzas de contacto internas genera una densidad de fuerzas de contacto o campo de tracción de Cauchy [5] que representa una distribución de las fuerzas de contacto internas en todo el volumen del cuerpo en una configuración particular del cuerpo en un momento dado . No es un campo vectorial porque depende no sólo de la posición de un punto material particular, sino también de la orientación local del elemento de superficie definida por su vector normal . [9]

Dependiendo de la orientación del plano considerado, el vector de tensión puede no ser necesariamente perpendicular a ese plano, es decir , paralelo a , y puede resolverse en dos componentes (Figura 2.1c):

¿Dónde está la componente normal de la fuerza al área diferencial?
donde es la componente tangencial de la fuerza a la superficie diferencial . El esfuerzo cortante se puede descomponer aún más en dos vectores mutuamente perpendiculares.

postulado de cauchy

Según el Postulado de Cauchy , el vector tensión permanece sin cambios para todas las superficies que pasan por el punto y tienen el mismo vector normal en , [7] [10] es decir, tienen una tangente común en . Esto significa que el vector de tensión es función únicamente del vector normal y no está influenciado por la curvatura de las superficies internas.

El lema fundamental de Cauchy

Una consecuencia del postulado de Cauchy es el Lema Fundamental de Cauchy , [1] [7] [11] también llamado teorema recíproco de Cauchy , [12] : p.103–130  que establece que los vectores de tensión que actúan en lados opuestos de la misma superficie son iguales en magnitud y opuestas en dirección. El lema fundamental de Cauchy es equivalente a la tercera ley del movimiento de acción y reacción de Newton , y se expresa como

Teorema de tensión de Cauchy: tensor de tensión

El estado de tensión en un punto del cuerpo se define entonces por todos los vectores de tensión T ( n ) asociados con todos los planos (en número infinito) que pasan por ese punto. [13] Sin embargo, de acuerdo con el teorema fundamental de Cauchy , [11] también llamado teorema de tensión de Cauchy , [1] simplemente conociendo los vectores de tensión en tres planos mutuamente perpendiculares, el vector de tensión en cualquier otro plano que pase por ese punto se puede encontrar a través de ecuaciones de transformación de coordenadas.

El teorema de tensión de Cauchy establece que existe un campo tensor de segundo orden σ ( x , t), llamado tensor de tensión de Cauchy, independiente de n , tal que T es una función lineal de n :

Esta ecuación implica que el vector de tensión T ( n ) en cualquier punto P en un continuo asociado con un plano con vector unitario normal n se puede expresar como una función de los vectores de tensión en los planos perpendiculares a los ejes de coordenadas, es decir, en términos de los componentes σ ij del tensor de tensión σ .

Para probar esta expresión, considere un tetraedro con tres caras orientadas en los planos coordenados y con un área infinitesimal d A orientada en una dirección arbitraria especificada por un vector unitario normal n (Figura 2.2). El tetraedro se forma cortando el elemento infinitesimal a lo largo de un plano arbitrario con unidad normal n . El vector de tensión en este plano se denota por T ( n ) . Los vectores de tensión que actúan sobre las caras del tetraedro se denotan como T ( e 1 ) , T ( e 2 ) y T ( e 3 ) , y son, por definición, los componentes σ ij del tensor de tensión σ . Este tetraedro a veces se llama tetraedro de Cauchy . El equilibrio de fuerzas, es decir , la primera ley del movimiento de Euler (segunda ley del movimiento de Newton), da:

Figura 2.2. Vector de tensión que actúa sobre un plano con vector unitario normal n .
Una nota sobre la convención de signos: el tetraedro se forma cortando un paralelepípedo a lo largo de un plano arbitrario n . Entonces, la fuerza que actúa sobre el plano n es la reacción ejercida por la otra mitad del paralelepípedo y tiene signo opuesto.

donde el lado derecho representa el producto de la masa encerrada por el tetraedro y su aceleración: ρ es la densidad, a es la aceleración y h es la altura del tetraedro, considerando el plano n como base. El área de las caras del tetraedro perpendicular a los ejes se puede encontrar proyectando d A en cada cara (usando el producto escalar):

y luego sustituyendo en la ecuación para cancelar d A :

Para considerar el caso límite cuando el tetraedro se reduce a un punto, h debe llegar a 0 (intuitivamente, el plano n se traslada a lo largo de n hacia O ). Como resultado, el lado derecho de la ecuación tiende a 0, por lo que

Suponiendo un elemento material (Figura 2.3) con planos perpendiculares a los ejes de coordenadas de un sistema de coordenadas cartesiano, los vectores de tensión asociados con cada uno de los planos del elemento, es decir, T ( e 1 ) , T ( e 2 ) y T ( e 3) . ) se puede descomponer en una componente normal y dos componentes de corte, es decir, componentes en la dirección de los tres ejes de coordenadas. Para el caso particular de una superficie con un vector unitario normal orientado en la dirección del eje x 1 , denotemos la tensión normal como σ 11 , y las dos tensiones cortantes como σ 12 y σ 13 :

En notación de índice esto es

Los nueve componentes σ ij de los vectores de tensión son los componentes de un tensor cartesiano de segundo orden llamado tensor de tensión de Cauchy , que puede usarse para definir completamente el estado de tensión en un punto y está dado por

donde σ 11 , σ 22 y σ 33 son tensiones normales, y σ 12 , σ 13 , σ 21 , σ 23 , σ 31 y σ 32 son tensiones cortantes. El primer índice i indica que la tensión actúa en un plano normal al eje X i , y el segundo índice j denota la dirección en la que actúa la tensión (por ejemplo, σ 12 implica que la tensión actúa en el plano que está normal al primer eje, es decir, X 1 y actúa a lo largo del segundo eje , es decir, X 2 ). Una componente de tensión es positiva si actúa en la dirección positiva de los ejes coordenados y si el plano donde actúa tiene un vector normal hacia afuera que apunta en la dirección coordenada positiva.

Por lo tanto, usando los componentes del tensor de tensión.

o equivalente,

Alternativamente, en forma matricial tenemos

La representación en notación Voigt del tensor de tensión de Cauchy aprovecha la simetría del tensor de tensión para expresar la tensión como un vector de seis dimensiones de la forma:

La notación de Voigt se utiliza ampliamente para representar las relaciones tensión-deformación en mecánica de sólidos y para la eficiencia computacional en software de mecánica estructural numérica.

Regla de transformación del tensor de tensión.

Se puede demostrar que el tensor de tensión es un tensor contravariante de segundo orden, que es una declaración de cómo se transforma ante un cambio del sistema de coordenadas. De un sistema x i a un sistema x i ' , los componentes σ ij en el sistema inicial se transforman en los componentes σ ij ' en el nuevo sistema de acuerdo con la regla de transformación tensorial (Figura 2.4):

donde A es una matriz de rotación con componentes a ij . En forma matricial esto es

Figura 2.4 Transformación del tensor de tensión.

Expandiendo la operación matricial y simplificando términos usando la simetría del tensor de tensión , se obtiene

El círculo de Mohr para tensiones es una representación gráfica de esta transformación de tensiones.

Esfuerzos normales y cortantes.

La magnitud de la componente de tensión normal σ n de cualquier vector de tensión T ( n ) que actúa sobre un plano arbitrario con el vector unitario normal n en un punto dado, en términos de las componentes σ ij del tensor de tensión σ , es el producto escalar de el vector de tensión y el vector unitario normal:

La magnitud del componente de tensión cortante τ n , que actúa ortogonal al vector n , se puede encontrar usando el teorema de Pitágoras :

dónde

Leyes del equilibrio: ecuaciones de movimiento de Cauchy

Figura 4. Cuerpo continuo en equilibrio.

Primera ley del movimiento de Cauchy

Según el principio de conservación del momento lineal , si el cuerpo continuo está en equilibrio estático se puede demostrar que las componentes del tensor de tensión de Cauchy en cada punto material del cuerpo satisfacen las ecuaciones de equilibrio:

,

dónde

Por ejemplo, para un fluido hidrostático en condiciones de equilibrio, el tensor de tensión toma la forma:

donde es la presión hidrostática y es el delta de kronecker .

Segunda ley del movimiento de Cauchy

Según el principio de conservación del momento angular , el equilibrio requiere que la sumatoria de momentos respecto de un punto arbitrario sea cero, lo que lleva a la conclusión de que el tensor de tensiones es simétrico , teniendo así sólo seis componentes de tensión independientes, en lugar del original. nueve:

Sin embargo, en presencia de tensiones de par, es decir, momentos por unidad de volumen, el tensor de tensiones no es simétrico. Este también es el caso cuando el número de Knudsen es cercano a uno, o el continuo es un fluido no newtoniano, lo que puede conducir a fluidos rotacionalmente no invariantes, como los polímeros .

Estreses principales e invariantes de estrés.

Componentes de tensión en un elemento giratorio 2D . Ejemplo de cómo varían los componentes de tensión en las caras (bordes) de un elemento rectangular a medida que varía el ángulo de su orientación. Las tensiones principales ocurren cuando las tensiones cortantes desaparecen simultáneamente de todas las caras. La orientación en la que esto ocurre da las direcciones principales . En este ejemplo, cuando el rectángulo es horizontal, las tensiones vienen dadas por

En cada punto de un cuerpo estresado existen al menos tres planos, llamados planos principales , con vectores normales , llamados direcciones principales , donde el vector de tensión correspondiente es perpendicular al plano, es decir, paralelo o en la misma dirección que el vector normal . y donde no existen esfuerzos cortantes normales . Las tres tensiones normales a estos planos principales se denominan tensiones principales .

Las componentes del tensor de tensiones dependen de la orientación del sistema de coordenadas en el punto considerado. Sin embargo, el tensor de tensión en sí es una cantidad física y, como tal, es independiente del sistema de coordenadas elegido para representarlo. Hay ciertas invariantes asociadas con cada tensor que también son independientes del sistema de coordenadas. Por ejemplo, un vector es un tensor simple de rango uno. En tres dimensiones, tiene tres componentes. El valor de estos componentes dependerá del sistema de coordenadas elegido para representar el vector, pero la magnitud del vector es una cantidad física (un escalar) y es independiente del sistema de coordenadas cartesiano elegido para representar el vector (siempre que sea normal ). De manera similar, cada tensor de segundo rango (como los tensores de tensión y deformación) tiene tres cantidades invariantes independientes asociadas. Un conjunto de tales invariantes son las tensiones principales del tensor de tensiones, que son simplemente los valores propios del tensor de tensiones. Sus vectores de dirección son las direcciones principales o vectores propios .

Un vector de tensión paralelo al vector unitario normal viene dado por:

donde es una constante de proporcionalidad, y en este caso particular corresponde a las magnitudes de los vectores de tensiones normales o tensiones principales.

Sabiendo que y tenemos

Se trata de un sistema homogéneo , es decir igual a cero, de tres ecuaciones lineales donde están las incógnitas. Para obtener una solución no trivial (distinta de cero) para , la matriz determinante de los coeficientes debe ser igual a cero, es decir, el sistema es singular. De este modo,

Desarrollar el determinante conduce a la ecuación característica.

dónde

La ecuación característica tiene tres raíces reales , es decir, no imaginarias debido a la simetría del tensor de tensiones. Los , y , son las tensiones principales, funciones de los valores propios . Los valores propios son las raíces del polinomio característico . Las tensiones principales son únicas para un tensor de tensiones determinado. Por lo tanto, a partir de la ecuación característica, los coeficientes , y , llamados primera, segunda y tercera invariantes de tensión , respectivamente, siempre tienen el mismo valor independientemente de la orientación del sistema de coordenadas.

Para cada valor propio, hay una solución no trivial en la ecuación . Estas soluciones son las direcciones principales o vectores propios que definen el plano donde actúan las tensiones principales. Las tensiones principales y las direcciones principales caracterizan la tensión en un punto y son independientes de la orientación.

Un sistema de coordenadas con ejes orientados a las direcciones principales implica que las tensiones normales son las tensiones principales y el tensor de tensiones está representado por una matriz diagonal:

Las tensiones principales se pueden combinar para formar las invariantes de tensión, , y . El primer y tercer invariante son la traza y el determinante, respectivamente, del tensor de tensión. De este modo,

Debido a su simplicidad, el sistema de coordenadas principales suele ser útil cuando se considera el estado del medio elástico en un punto particular. Las tensiones principales a menudo se expresan en la siguiente ecuación para evaluar tensiones en las direcciones x e y o tensiones axiales y de flexión en una pieza. [14] : p.58–59  Las tensiones normales principales se pueden utilizar para calcular la tensión de von Mises y, en última instancia, el factor de seguridad y el margen de seguridad.

Usar solo la parte de la ecuación debajo de la raíz cuadrada es igual al esfuerzo cortante máximo y mínimo para más y menos. Esto se muestra como:

Esfuerzos cortantes máximos y mínimos.

El esfuerzo cortante máximo o esfuerzo cortante principal máximo es igual a la mitad de la diferencia entre los esfuerzos principales mayor y menor, y actúa sobre el plano que biseca el ángulo entre las direcciones de los esfuerzos principales mayor y menor, es decir, el plano de la El esfuerzo cortante máximo está orientado desde los planos de esfuerzos principales. El esfuerzo cortante máximo se expresa como

Suponiendo entonces

Cuando el tensor de tensión es distinto de cero, la componente de tensión normal que actúa en el plano para el esfuerzo cortante máximo es distinta de cero y es igual a

Tensor desviador de tensión

El tensor de tensión se puede expresar como la suma de otros dos tensores de tensión:

  1. un tensor de tensión hidrostático medio o un tensor de tensión volumétrico o un tensor de tensión normal medio , que tiende a cambiar el volumen del cuerpo estresado; y
  2. un componente desviador llamado tensor desviador de tensión , que tiende a distorsionarlo.

Entonces

¿Dónde está la tensión media dada por

La presión ( ) generalmente se define como un tercio negativo de la traza del tensor de tensión menos cualquier tensión con la que contribuya la divergencia de la velocidad, es decir

donde es una constante de proporcionalidad (es decir, el primero de los parámetros de Lamé ), es el operador de divergencia , es la k :ésima coordenada cartesiana , es la velocidad del flujo y es la k :ésima componente cartesiana de .

El tensor de tensión desviador se puede obtener restando el tensor de tensión hidrostática del tensor de tensión de Cauchy:

Invariantes del tensor desviador de tensiones

Al ser un tensor de segundo orden, el tensor desviador de tensiones también tiene un conjunto de invariantes , que se pueden obtener mediante el mismo procedimiento que se utiliza para calcular las invariantes del tensor de tensiones. Se puede demostrar que las direcciones principales del tensor desviador de tensión son las mismas que las direcciones principales del tensor de tensión . Por tanto, la ecuación característica es

donde , y son los invariantes de tensión desviadores primero, segundo y tercero , respectivamente. Sus valores son los mismos (invariantes) independientemente de la orientación del sistema de coordenadas elegido. Estas invariantes de tensión desviatorias se pueden expresar en función de los componentes de o sus valores principales , y , o alternativamente, en función de o sus valores principales , y . De este modo,

Porque , el tensor desviador de tensión está en un estado de corte puro.

En mecánica de sólidos se utiliza comúnmente una cantidad llamada tensión equivalente o tensión de von Mises . La tensión equivalente se define como

Tensiones octaédricas

Figura 6. Planos de tensiones octaédricos

Considerando las direcciones principales como ejes coordenados, un plano cuyo vector normal forma ángulos iguales con cada uno de los ejes principales (es decir, que tiene cosenos directores iguales a ) se llama plano octaédrico . Hay un total de ocho planos octaédricos (Figura 6). Los componentes normal y cortante del tensor de tensión en estos planos se denominan tensión normal octaédrica y tensión cortante octaédrica , respectivamente. El plano octaédrico que pasa por el origen se conoce como plano π ( π no debe confundirse con la tensión media indicada por π en la sección anterior) . En el plano π , .

Sabiendo que el tensor de tensiones del punto O (Figura 6) en los ejes principales es

el vector tensión en un plano octaédrico viene dado por:

La componente normal del vector de tensión en el punto O asociado con el plano octaédrico es

que es la tensión normal media o tensión hidrostática. Este valor es el mismo en los ocho planos octaédricos. El esfuerzo cortante en el plano octaédrico es entonces

Ver también

Notas

  1. ^ En detalle:

Referencias

  1. ^ abc Fridtjov Irgens (2008), "Mecánica continua". Saltador. ISBN  3-540-74297-2
  2. ^ Truesdell, C .; Toupin, RA (1960), "The Classical Field Theories", en Flügge, Siegfried (ed.), Principios de la mecánica clásica y la teoría de campos/Prinzipien der Klassischen Mechanik und Feldtheorie , Handbuch der Physik (Enciclopedia de física), vol. III/1, Berlín– Heidelberg –Nueva York: Springer-Verlag , págs. 226–793, Bibcode :1960HDP.....2.....F, doi :10.1007/978-3-642-45943-6 , ISBN 978-3-540-02547-4, SEÑOR  0118005, Zbl  0118.39702.
  3. ^ Peter Chadwick (1999), "Mecánica continua: teoría y problemas concisos". Publicaciones de Dover, serie "Libros de física". ISBN 0-486-40180-4 . paginas 
  4. ^ ab Yuan-cheng Fung y Pin Tong (2001) "Mecánica de sólidos clásica y computacional". Científico mundial. ISBN 981-02-4124-0 
  5. ^ ab Smith y Truesdell p.97
  6. ^ ab G. Thomas Mase y George E. Mase (1999), "Mecánica continua para ingenieros" (segunda edición). Prensa CRC. ISBN 0-8493-1855-6 
  7. ^ abc I-Shih Liu (2002), "Mecánica continua". Saltador ISBN 3-540-43019-9 
  8. ^ ab Han-Chin Wu (2005), "Mecánica continua y plasticidad". Prensa CRC. ISBN 1-58488-363-4 
  9. ^ Lubliner
  10. ^ basar
  11. ^ abc Teodor M. Atanackovic y Ardéshir Guran (2000), "Teoría de la elasticidad para científicos e ingenieros". Saltador. ISBN 0-8176-4072-X 
  12. ^ Keith D. Hjelmstad (2005), "Fundamentos de la mecánica estructural" (segunda edición). Prentice Hall. ISBN 0-387-23330-X 
  13. ^ ab Wai-Fah Chen y Da-Jian Han (2007), "Plasticidad para ingenieros estructurales". Publicación J. Ross ISBN 1-932159-75-4 
  14. ^ Bernard Hamrock (2005), "Fundamentos de los elementos de la máquina". McGraw-Hill. ISBN 0-07-297682-9 
  15. ^ Rabindranath Chatterjee (1999), "Teoría matemática de la mecánica del continuo". Ciencia alfa. ISBN 81-7319-244-8 
  16. ^ John Conrad Jaeger, NGW Cook y RW Zimmerman (2007), "Fundamentos de la mecánica de rocas" (cuarta edición). Wiley-Blackwell. ISBN 0-632-05759-9 
  17. ^ Mohammed Ameen (2005), "Elasticidad computacional: teoría de la elasticidad y métodos de elementos finitos y límite" (libro). Ciencia alfa, ISBN 1-84265-201-X 
  18. ^ William Prager (2004), "Introducción a la mecánica de Continua". Publicaciones de Dover. ISBN 0-486-43809-0