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 :
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.
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
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á representada 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):
una normal al plano, llamada tensión normal
¿Dónde está el componente normal de la fuerza al área diferencial?
donde es el componente tangencial de la fuerza al área superficial 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:
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 se acerca 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
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
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:
Consideremos un cuerpo continuo (véase la Figura 4) que ocupa un volumen , que tiene un área de superficie , con fuerzas de tracción o de superficie definidas por unidad de área que actúan sobre cada punto de la superficie del cuerpo, y fuerzas corporales por unidad de volumen sobre cada punto dentro del volumen . Por lo tanto, si el cuerpo está en equilibrio , la fuerza resultante que actúa sobre el volumen es cero, por lo tanto:
Para un volumen arbitrario la integral se desvanece y tenemos las ecuaciones de equilibrio
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:
Derivación de la simetría del tensor de tensión
Sumando los momentos respecto al punto O (Figura 4), el momento resultante es cero, ya que el cuerpo está en equilibrio. Por lo tanto,
donde es el vector de posición y se expresa como
Sabiendo eso y utilizando el teorema de divergencia de Gauss para pasar de una integral de superficie a una integral de volumen, tenemos
La segunda integral es cero, ya que contiene las ecuaciones de equilibrio. Esto deja la primera integral, donde , por lo tanto
Para un volumen arbitrario V, entonces tenemos
que se cumple en cada punto del cuerpo. Desarrollando esta ecuación tenemos
, , y
o en general
Esto demuestra que el tensor de tensión es simétrico.
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
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
Derivación de las tensiones cortantes máximas y mínimas [8] : p.45–78 [11] : p.1–46 [13] [15] : p.111–157 [16] : p.9–41 [17] : p.33–66 [18] : p.43–61
La tensión normal se puede escribir en términos de tensiones principales como
Sabiendo que , la tensión cortante en términos de componentes de tensiones principales se expresa como
El esfuerzo cortante máximo en un punto de un cuerpo continuo se determina maximizando el sujeto a la condición de que
Este es un problema de maximización con restricciones, que se puede resolver utilizando la técnica del multiplicador de Lagrange para convertir el problema en un problema de optimización sin restricciones. Por lo tanto, los valores estacionarios (valores máximo y mínimo) de ocurren donde el gradiente de es paralelo al gradiente de .
La función lagrangiana para este problema se puede escribir como
donde es el multiplicador lagrangiano (que es diferente del que se usa para denotar valores propios).
Los valores extremos de estas funciones son
de allí
Estas tres ecuaciones junto con la condición se pueden resolver para y
Al multiplicar las tres primeras ecuaciones por y , respectivamente, y sabiendo que obtenemos
Sumando estas tres ecuaciones obtenemos
Este resultado se puede sustituir en cada una de las tres primeras ecuaciones para obtener
Haciendo lo mismo para las otras dos ecuaciones que tenemos
Una primera aproximación para resolver estas tres últimas ecuaciones es considerar la solución trivial . Sin embargo, esta opción no cumple la restricción .
Considerando la solución donde y , se determina a partir de la condición de que , entonces de la ecuación original para se ve que . Los otros dos valores posibles para se pueden obtener de manera similar suponiendo
y
y
Por lo tanto, un conjunto de soluciones para estas cuatro ecuaciones es:
Estos corresponden a valores mínimos y verifican que no existen tensiones cortantes en los planos normales a las direcciones principales de tensión, como se mostró anteriormente.
Un segundo conjunto de soluciones se obtiene suponiendo que y . Por lo tanto, tenemos
Para encontrar los valores de y primero sumamos estas dos ecuaciones
Sabiendo que para
y
tenemos
y resolviendo tenemos
Luego, resolviendo tenemos
y
Los otros dos valores posibles para se pueden obtener de manera similar suponiendo
y
y
Por lo tanto, el segundo conjunto de soluciones para , que representa un máximo para es
Por lo tanto, suponiendo , el esfuerzo cortante máximo se expresa por
y puede afirmarse que es igual a la mitad de la diferencia entre las tensiones principales mayor y menor, actuando sobre el plano que biseca el ángulo entre las direcciones de las tensiones principales mayor y menor.
Tensor desviador de tensión
El tensor de tensión se puede expresar como la suma de otros dos tensores de tensión:
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
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,
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
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
^ abc Fridtjov Irgens (2008), "Mecánica continua". Saltador. ISBN 3-540-74297-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 , ISBN978-3-540-02547-4, MR 0118005, Zbl 0118.39702.
^ 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
^ ab Yuan-cheng Fung y Pin Tong (2001) "Mecánica de sólidos clásica y computacional". World Scientific. ISBN 981-02-4124-0
^ de Smith y Truesdell, pág. 97
^ 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
^ abc I-Shih Liu (2002), "Mecánica del medio continuo". Springer ISBN 3-540-43019-9
^ de Han-Chin Wu (2005), "Mecánica del medio continuo y plasticidad". CRC Press. ISBN 1-58488-363-4
^ Revestimiento de pintura
^ Basar
^ 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
^ Keith D. Hjelmstad (2005), "Fundamentos de mecánica estructural" (2.ª edición). Prentice-Hall. ISBN 0-387-23330-X
^ ab Wai-Fah Chen y Da-Jian Han (2007), "Plasticidad para ingenieros estructurales". J. Ross Publishing ISBN 1-932159-75-4
^ Bernard Hamrock (2005), "Fundamentos de los elementos de las máquinas". McGraw-Hill. ISBN 0-07-297682-9
^ Rabindranath Chatterjee (1999), "Teoría matemática de la mecánica del medio continuo". Alpha Science. ISBN 81-7319-244-8
^ 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
^ 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
^ William Prager (2004), "Introducción a la mecánica de los continuos". Publicaciones de Dover. ISBN 0-486-43809-0