stringtranslate.com

Tensor de tensión de Cauchy

En mecánica de medios continuos , 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 tensión en un punto dentro de un material en el 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), que corresponde al escalar de tensión. El vector unitario es adimensional .

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

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 grandes deformaciones, 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 tensiones 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 suma de momentos con respecto a un punto arbitrario sea cero, lo que lleva a la conclusión de que el tensor de tensiones es simétrico , por lo que tiene solo seis componentes de tensión independientes, en lugar de los nueve originales. Sin embargo, en presencia de tensiones de pareja, 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 .

Existen ciertas invariantes asociadas al tensor de tensiones, cuyos valores no dependen del sistema de coordenadas elegido ni del elemento de área sobre el que actúa el tensor de tensiones. Se trata de 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 sobre 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 tensiones de par sobre 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 sobre 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 descomponerse en dos componentes: un componente normal al plano, llamado tensión normal , y otro componente paralelo a este plano, llamado tensión cortante .

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

Para formular el principio de tensión de Euler-Cauchy, considérese una superficie imaginaria que pasa por 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 utilizar 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 que 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 puede expresarse 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 está sometido 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 en el cuerpo, y de un segmento al otro a través de la superficie divisoria , debido al contacto mecánico de una porción del continuo sobre la otra (Figura 2.1a y 2.1b). En un elemento de área que contiene , con vector normal , la distribución de fuerza es equipolenta a una fuerza de contacto ejercida en el punto P y 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 a medida que se vuelve muy pequeño y tiende a cero, la relación se convierte en y el vector de tensión de par se desvanece. En campos específicos de la mecánica de medios continuos se supone que la tensión de par no se desvanece; sin embargo, las ramas clásicas de la mecánica de medios continuos abordan materiales no polares que no consideran las tensiones de par y los momentos corporales.

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

Esta ecuación significa que el vector de 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 fuerza 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 solo de la posición de un punto material particular, sino también de la orientación local del elemento de superficie según lo definido 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 descomponerse en dos componentes (Figura 2.1c):

¿Dónde está el componente normal de la fuerza al área diferencial?
donde es el componente tangencial de la fuerza al área de superficie diferencial . La tensión cortante se puede descomponer además en dos vectores mutuamente perpendiculares.

Postulado de Cauchy

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

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 sobre lados opuestos de la misma superficie son iguales en magnitud y opuestos en dirección. El lema fundamental de Cauchy es equivalente a la tercera ley de Newton del movimiento de acción y reacción, 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 a todos los planos (infinitos en número) que pasan por ese punto. [13] Sin embargo, según 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 tensorial 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 de un continuo asociado a un plano con vector unitario normal n puede expresarse 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 demostrar esta expresión, considere un tetraedro con tres caras orientadas en los planos de coordenadas, 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 una normal unitaria 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 de movimiento de Euler (segunda ley de movimiento de Newton), da:

Figura 2.2. Vector de tensión que actúa sobre un plano con un 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 . Por lo tanto, la fuerza que actúa sobre el plano n es la reacción ejercida por la otra mitad del paralelepípedo y tiene un 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 perpendiculares a los ejes se puede hallar 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 a medida que el tetraedro se encoge hasta un punto, h debe tender 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 (ver figura en la parte superior de la página) con planos perpendiculares a los ejes de coordenadas de un sistema de coordenadas cartesianas, los vectores de tensión asociados a cada uno de los planos del elemento, es decir T ( e 1 ) , T ( e 2 ) y T ( e 3 ) se pueden descomponer en un 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 vector unitario normal orientado en la dirección del eje x 1 , denotemos la tensión normal por σ 11 , y las dos tensiones de corte 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 se puede utilizar 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 sobre 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 sobre el plano que es normal al 1 er eje es decir; X 1 y actúa a lo largo del 2 º eje es decir; X 2 ). Un componente de tensión es positivo si actúa en la dirección positiva de los ejes de coordenadas, y si el plano donde actúa tiene un vector normal hacia afuera que apunta en la dirección de coordenadas positiva.

Así, utilizando los componentes del tensor de tensión

o, equivalentemente,

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 Voigt se utiliza ampliamente para representar relaciones de 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, lo que es una declaración de cómo se transforma bajo un cambio del sistema de coordenadas. De un sistema x i a un sistema x i ' , las componentes σ ij en el sistema inicial se transforman en las componentes σ ij ' en el nuevo sistema de acuerdo con la regla de transformación de tensores (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

Desarrollando la operación matricial y simplificando los términos utilizando la simetría del tensor de tensión , obtenemos

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 del componente de tensión normal σ n de cualquier vector de tensión T ( n ) que actúa sobre un plano arbitrario con vector unitario normal n en un punto dado, en términos de los componentes σ ij del tensor de tensión σ , es el producto escalar del vector de tensión y el vector unitario normal:

La magnitud del componente de esfuerzo cortante τ n , que actúa ortogonalmente al vector n , se puede encontrar utilizando el teorema de Pitágoras :

dónde

Leyes de equilibrio: ecuaciones de movimiento de Cauchy

Figura 4. Cuerpo continuo en equilibrio

Primera ley del movimiento de Cauchy

De acuerdo con 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 tensiones 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 suma de momentos con respecto a 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 pares de tensiones, es decir, momentos por unidad de volumen, el tensor de tensiones no es simétrico. Esto también sucede cuando el número de Knudsen es cercano a uno , o el continuo es un fluido no newtoniano, lo que puede dar lugar a fluidos rotacionalmente no invariantes, como los polímeros .

Tensiones principales e invariantes de tensión

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

En cada punto de un cuerpo sometido a tensión 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 tensiones cortantes normales . Las tres tensiones normales a estos planos principales se denominan tensiones principales .

Los componentes del tensor de tensión dependen de la orientación del sistema de coordenadas en el punto considerado. Sin embargo, el tensor de tensión en sí mismo es una cantidad física y, como tal, es independiente del sistema de coordenadas elegido para representarlo. Hay ciertos invariantes asociados 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 cartesianas 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 con él. Un conjunto de tales invariantes son las tensiones principales del tensor de tensión, que son simplemente los valores propios del tensor de tensión. 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 tensión normal o tensiones principales.

Sabiendo que y , tenemos

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

La expansión del 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. Las , 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 dado. Por lo tanto, a partir de la ecuación característica, los coeficientes , y , llamados primero, segundo y tercer invariantes de tensiones , respectivamente, siempre tienen el mismo valor independientemente de la orientación del sistema de coordenadas.

Para cada valor propio, existe una solución no trivial para 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 los invariantes de tensión, , , y . El primer y tercer invariante son la traza y el determinante, respectivamente, del tensor de tensión. Por lo tanto,

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

Si se utiliza solo la parte de la ecuación que se encuentra debajo de la raíz cuadrada , se obtiene el esfuerzo cortante máximo y mínimo para los valores positivos y negativos. 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 el esfuerzo principal mayor y menor, y actúa sobre el plano que biseca el ángulo entre las direcciones del esfuerzo principal mayor y menor, es decir, el plano del esfuerzo cortante máximo está orientado a partir de los planos de esfuerzo principal. El esfuerzo cortante máximo se expresa como

Suponiendo entonces

Cuando el tensor de tensión no es cero, el componente de tensión normal que actúa sobre el plano para la tensión cortante máxima no es 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ática media o un tensor de tensión volumétrica o un tensor de tensión normal media , 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 ( ) se define generalmente como un tercio negativo de la traza del tensor de tensión menos cualquier tensión con la que contribuye 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 el k :ésimo componente cartesiano 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

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

donde , y son los invariantes de tensión desviatorios primero, segundo y tercero , respectivamente. Sus valores son los mismos (invariantes) independientemente de la orientación del sistema de coordenadas elegido. Estos invariantes de tensión desviatorios se pueden expresar como una función de los componentes de o sus valores principales , , y , o alternativamente, como una función de o sus valores principales , , y . Por lo tanto,

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

En mecánica de sólidos se utiliza habitualmente una cantidad denominada tensión equivalente o tensión de von Mises . La tensión equivalente se define como

Tensiones octaédricas

Figura 6. Planos de tensión octaédricos

Considerando las direcciones principales como los ejes de coordenadas, 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 denotada 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 de tensión en un plano octaédrico viene dado por:

El 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. La tensión cortante en el plano octaédrico es entonces

Véase 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, MR  0118005, Zbl  0118.39702.
  3. ^ Peter Chadwick (1999), "Mecánica del medio continuo: teoría concisa y problemas". Dover Publications, serie "Libros de física". ISBN 0-486-40180-4 . páginas 
  4. ^ ab Yuan-cheng Fung y Pin Tong (2001) "Mecánica de sólidos clásica y computacional". World Scientific. ISBN 981-02-4124-0 
  5. ^ de Smith y Truesdell, pág. 97
  6. ^ ab G. Thomas Mase y George E. Mase (1999), "Mecánica de medios continuos para ingenieros" (2.ª edición). CRC Press. ISBN 0-8493-1855-6 
  7. ^ abc I-Shih Liu (2002), "Mecánica del medio continuo". Springer ISBN 3-540-43019-9 
  8. ^ de Han-Chin Wu (2005), "Mecánica del medio continuo y plasticidad". CRC Press. ISBN 1-58488-363-4 
  9. ^ Revestimiento de pintura
  10. ^ Basar
  11. ^ abc Teodor M. Atanackovic y Ardéshir Guran (2000), "Teoría de la elasticidad para científicos e ingenieros". Springer. ISBN 0-8176-4072-X 
  12. ^ Keith D. Hjelmstad (2005), "Fundamentos de mecánica estructural" (2.ª edición). Prentice-Hall. ISBN 0-387-23330-X 
  13. ^ ab Wai-Fah Chen y Da-Jian Han (2007), "Plasticidad para ingenieros estructurales". J. Ross Publishing ISBN 1-932159-75-4 
  14. ^ Bernard Hamrock (2005), "Fundamentos de los elementos de las máquinas". McGraw-Hill. ISBN 0-07-297682-9 
  15. ^ Rabindranath Chatterjee (1999), "Teoría matemática de la mecánica del medio continuo". Alpha Science. ISBN 81-7319-244-8 
  16. ^ John Conrad Jaeger, NGW Cook y RW Zimmerman (2007), "Fundamentos de la mecánica de rocas" (4.ª 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 de contorno" (libro). Alpha Science, ISBN 1-84265-201-X 
  18. ^ William Prager (2004), "Introducción a la mecánica de los continuos". Publicaciones de Dover. ISBN 0-486-43809-0