Teoría de la deformación finita

De testwiki
Ir a la navegación Ir a la búsqueda
Figura 1. Movimiento de un cuerpo continuo
Figura 2. Deformación de un cuerpo continuo

En mecánica de medios continuos, la teoría de la deformación finita, también llamada teoría de grandes deformaciones analiza el comportamiento de los sólidos frente a deformaciones y/o rotaciones lo suficientemente grandes como para invalidar los supuestos inherentes a la teoría de la deformación infinitesimal. En este caso, las configuraciones no deformadas y deformadas del medio continuo son significativamente diferentes, lo que requiere una distinción clara entre ellas. Este tipo de deformaciones son propias de materiales como elastómeros, plásticos, fluidos viscosos y, en el campo de la biología, los tejidos blandos.

Campo de desplazamiento

Plantilla:AP El desplazamiento de un cuerpo puede tener dos componentes: un desplazamiento rígido y una deformación.

  • El desplazamiento de un cuerpo rígido consiste en una traslación y/o una rotación simultáneas del cuerpo sin cambiar su forma o tamaño.
  • La deformación implica el cambio de forma y/o tamaño del cuerpo desde una configuración inicial o no deformada κ0() a una configuración actual o deformada κt() (Figura 1).

Un cambio en la configuración de un cuerpo continuo puede describirse mediante un campo de desplazamiento. Un campo de desplazamiento es un campo vectorial que describe todos los vectores de desplazamiento de todas las partículas del cuerpo, y que relaciona la configuración deformada con la configuración no deformada. La distancia entre dos partículas cualesquiera cambia si y solo si se ha producido una deformación. Si el desplazamiento se produce sin deformación, entonces se trata del desplazamiento de un cuerpo rígido.

Tensor del gradiente de deformación

El tensor del gradiente de deformación 𝐅(𝐗,t)=FjK𝐞j𝐈K está relacionado tanto con la configuración de referencia como con la actual, expresada según los vectores unitarios 𝐞j y 𝐈K. Por lo tanto, es un tensor de punto a punto.

Se pueden definir dos tipos de tensor de gradiente de deformación.

Debido al supuesto de continuidad de χ(𝐗,t), 𝐅 tiene el 𝐇=𝐅1 inverso, donde 𝐇 es el tensor de gradiente de deformación espacial. Entonces, por el teorema de la función implícita,[1] el determinante jacobiano J(𝐗,t) debe ser no singular, es decir J(𝐗,t)=det𝐅(𝐗,t)0

El tensor del gradiente de deformación del material 𝐅(𝐗,t)=FjK𝐞j𝐈K es un tensor de segundo orden que representa el gradiente de la función de aplicación o relación funcional χ(𝐗,t), que describe el desplazamiento en un medio continuo. El tensor de gradiente de deformación del material caracteriza la deformación local en un punto material con el vector de posición 𝐗, es decir, la deformación en puntos vecinos, transformando (mediante una aplicación lineal) un elemento lineal del material que conecta la configuración de referencia de cada punto con su configuración actual o deformada, asumiendo la continuidad de la función χ(𝐗,t), es decir, considerando que es diferenciable respecto a la posición 𝐗 y al tiempo t, lo que implica que no se producen fracturas y que no se abren ni cierran huecos durante la deformación. Así, se tiene que

d𝐱=𝐱𝐗d𝐗dxj=xjXKdXK=χ(𝐗,t)d𝐗dxj=FjKdXK.=𝐅(𝐗,t)d𝐗

Vector de desplazamiento relativo

Considérese una partícula o punto material P con el vector de posición 𝐗=XI𝐈I en la configuración no deformada (Figura 2). Después de un desplazamiento del cuerpo, la nueva posición de la partícula indicada por p en la nueva configuración viene dada por la posición del vector 𝐱=xi𝐞i. Los sistemas de coordenadas para la configuración deformada y no deformada se pueden superponer si así conviene.

Considérese ahora un punto material Q vecino de P, con vector de posición 𝐗+Δ𝐗=(XI+ΔXI)𝐈I. En la configuración deformada, esta partícula tiene una nueva posición q dada por el vector de posición 𝐱+Δ𝐱. Suponiendo que los segmentos de línea ΔX y Δ𝐱 que unen las partículas P y Q en la configuración no deformada y deformada, respectivamente, son muy pequeños, entonces se puede expresarlos como d𝐗 y d𝐱. Así, de la Figura 2 se deduce que

𝐱+d𝐱=𝐗+d𝐗+𝐮(𝐗+d𝐗)d𝐱=𝐗𝐱+d𝐗+𝐮(𝐗+d𝐗)=d𝐗+𝐮(𝐗+d𝐗)𝐮(𝐗)=d𝐗+d𝐮

donde 𝐝𝐮 es el vector de desplazamiento relativo, que representa el movimiento del punto Q con respecto a P en la configuración deformada.

Aproximación de Taylor

Para un elemento infinitesimal d𝐗, y suponiendo la continuidad en el campo de desplazamiento, es posible utilizar una serie de Taylor alrededor del punto P, despreciando términos de orden superior, para aproximar las componentes del vector de desplazamiento relativo para la partícula vecina Q como

𝐮(𝐗+d𝐗)=𝐮(𝐗)+d𝐮ui*=ui+dui𝐮(𝐗)+𝐗𝐮d𝐗ui*ui+uiXJdXJ.

Así, la ecuación anterior d𝐱=d𝐗+d𝐮 se puede escribir como

d𝐱=d𝐗+d𝐮=d𝐗+𝐗𝐮d𝐗=(𝐈+𝐗𝐮)d𝐗=𝐅d𝐗

Derivada respecto al tiempo del gradiente de deformación

Los cálculos que implican la deformación de un cuerpo dependiente del tiempo a menudo requieren que se calcule una derivada respecto al tiempo del gradiente de deformación. Una definición geométricamente consistente de dicha derivada requiere una incursión en la geometría diferencial,[2] aunque en el presente artículo no se considera este tipo de problemas.

La derivada respecto al tiempo de 𝐅 es

𝐅˙=𝐅t=t[𝐱(𝐗,t)𝐗]=𝐗[𝐱(𝐗,t)t]=𝐗[𝐕(𝐗,t)]

donde 𝐕 es la velocidad de deformación del material, por lo que la derivada del lado derecho representa el gradiente de la velocidad del material. Es común convertirlo en un gradiente espacial aplicando la regla de la cadena de las derivadas, es decir,

𝐅˙=𝐗[𝐕(𝐗,t)]=𝐗[𝐯(𝐱(𝐗,t),t)]=𝐱[𝐯(𝐱,t)]|𝐱=𝐱(𝐗,t)𝐱(𝐗,t)𝐗=𝒍𝐅

donde 𝒍=(𝐱𝐯)T es el gradiente de la velocidad espacial y donde 𝐯(𝐱,t)=𝐕(𝐗,t) es la velocidad espacial (euleriana) en 𝐱=𝐱(𝐗,t). Si el gradiente de la velocidad espacial es constante en el tiempo, la ecuación anterior se puede resolver exactamente para obtener

𝐅=e𝒍t

suponiendo que 𝐅=𝟏 en t=0. Existen varios métodos para calcular el valor exponencial anterior.

Cantidades relacionadas que se utilizan a menudo en mecánica continua son el tensor de velocidad de deformación y el tensor de giro definidos, respectivamente, como:

𝒅=12(𝒍+𝒍T),𝒘=12(𝒍𝒍T).

El tensor de velocidad de deformación proporciona la velocidad de estiramiento de los elementos lineales, mientras que el tensor de espín indica la velocidad de rotación o vorticidad del movimiento.

La derivada material-tiempo de la inversa del gradiente de deformación (manteniendo fija la configuración de referencia) a menudo se requiere en análisis que involucran deformaciones finitas. Esta derivada es

t(𝐅1)=𝐅1𝐅˙𝐅1.

La relación anterior se puede verificar tomando la derivada respecto al tiempo del material de 𝐅1d𝐱=d𝐗 y observando que 𝐗˙=0.

Descomposición polar del tensor de gradiente de deformación

Figura 3. Representación de la descomposición polar del gradiente de deformación

El gradiente de deformación 𝐅, como cualquier tensor invertible de segundo orden, se puede descomponer (utilizando el teorema de descomposición polar) en un producto de dos tensores de segundo orden (Truesdell y Noll, 1965): un tensor ortogonal y un tensor simétrico definido positivo, es decir,

𝐅=𝐑𝐔=𝐕𝐑

donde el tensor 𝐑 es un tensor ortogonal propio, es decir, 𝐑1=𝐑T y det𝐑=+1, que representan una rotación; el tensor 𝐔 es el tensor de estiramiento derecho; y 𝐕 el tensor de estiramiento izquierdo. Los términos derecha e izquierda significan que están a la derecha e izquierda del tensor de rotación 𝐑, respectivamente. 𝐔 y 𝐕 son ambos matrices definidas positivas, es decir, 𝐱𝐔𝐱>0 y 𝐱𝐕𝐱>0 para todos los 𝐱3 distintos de cero, y tensores simétricos, es decir, 𝐔=𝐔T y 𝐕=𝐕T, de segundo orden.

Esta descomposición implica que la deformación de un elemento lineal d𝐗 en la configuración no deformada sobre d𝐱 en la configuración deformada, es decir, d𝐱=𝐅d𝐗, se puede obtener estirando primero el elemento en 𝐔, es decir, d𝐱=𝐔d𝐗, seguido de una rotación 𝐑, es decir, d𝐱=𝐑d𝐱; o de manera equivalente, aplicando primero una rotación rígida 𝐑, es decir, d𝐱=𝐑d𝐗, seguida a continuación de un estiramiento 𝐕, es decir, d𝐱=𝐕d𝐱 (véase la Figura 3).

Debido a la ortogonalidad de 𝐑

𝐕=𝐑𝐔𝐑T

de modo que 𝐔 y 𝐕 tienen los mismos autovalores o deformaciones principales, pero diferentes vectores propios o direcciones principales 𝐍i y 𝐧i, respectivamente. Las direcciones principales están relacionadas por

𝐧i=𝐑𝐍i.

Esta descomposición polar, que es única ya que 𝐅 es invertible con un determinante positivo, es un corolario de la descomposición en valores singulares.

Transformación de un elemento de superficie y volumen

Para transformar cantidades definidas con respecto a áreas en una configuración deformada a aquellas relativas a áreas en una configuración de referencia, y viceversa, se usa la relación de Nanson, expresada como

da𝐧=JdA𝐅T𝐍

donde da es un área de una región en la configuración deformada, dA es la misma área en la configuración de referencia, y 𝐧 es la normal exterior al elemento de área en la configuración actual, mientras que 𝐍 es la normal exterior en la configuración de referencia, 𝐅 es el gradiente de deformación y J=det𝐅.

La fórmula correspondiente para la transformación del elemento de volumen es:dv=JdV

Plantilla:Demostración

Tensores de deformación fundamentales

Un tensor de deformación está definido según la Unión Internacional de Química Pura y Aplicada como:[3]

Plantilla:Cita

Dado que una rotación pura no debería inducir ninguna deformación en un cuerpo deformable, a menudo es conveniente utilizar medidas de deformación independientes de la rotación en mecánica de medios continuos. Como una rotación seguida de su rotación inversa no produce ningún cambio (𝐑𝐑T=𝐑T𝐑=𝐈), se puede excluir la rotación multiplicando el tensor de gradiente de deformación 𝐅 por su matriz transpuesta.

En mecánica se utilizan varios tensores de gradiente de deformación independientes de la rotación (o tensores de deformación, para abreviar). En mecánica de sólidos, los más habituales son los tensores de deformación de Cauchy-Green derecho e izquierdo.

Tensor de deformación de Cauchy (tensor de deformación de Cauchy-Green derecho)

En 1839, George Green introdujo un tensor de deformación conocido como tensor de deformación derecho de Cauchy-Green o tensor de deformación de Green (la Unión Internacional de Química Pura y Aplicada recomienda que este tensor se llame "tensor de deformación de Cauchy"),[3] definido como:

𝐂=𝐅T𝐅=𝐔2CIJ=FkIFkJ=xkXIxkXJ.

Físicamente, el tensor de Cauchy-Green proporciona el cuadrado del cambio local en distancias debido a la deformación, es decir, d𝐱2=d𝐗𝐂d𝐗

Las invariantes de 𝐂 se utilizan a menudo en las expresiones para caracterizar la función densidad de energía de deformación. Los invariantes más utilizados son

I1C:=tr(𝐂)=CII=λ12+λ22+λ32I2C:=12[(tr𝐂)2tr(𝐂2)]=12[(CJJ)2CIKCKI]=λ12λ22+λ22λ32+λ32λ12I3C:=det(𝐂)=J2=λ12λ22λ32.

donde J:=det𝐅 es el determinante del gradiente de deformación 𝐅 y λi son relaciones de estiramiento para las fibras unitarias que inicialmente están orientadas en las direcciones del vector propio del tensor de estiramiento derecho (de referencia, que generalmente no están alineados con los tres ejes de los sistemas de coordenadas).

Tensor de tensión de Finger

La Unión Internacional de Química Pura y Aplicada recomienda[3] que el inverso del tensor de deformación de Cauchy-Green derecho (llamado tensor de deformación de Cauchy en ese documento), es decir, 𝐂1, se denomine tensor de tensión de Finger. Sin embargo, esa nomenclatura no es universalmente aceptada en la mecánica aplicada.

𝐟=𝐂1=𝐅1𝐅TfIJ=XIxkXJxk

Tensor de deformación de Green (tensor de deformación de Cauchy-Green izquierdo)

Invertir el orden de multiplicación en la fórmula del tensor de deformación de Green-Cauchy derecho conduce al tensor de deformación de Cauchy-Green izquierdo, que se define como:

𝐁=𝐅𝐅T=𝐕2Bij=xiXKxjXK

El tensor de deformación de Cauchy-Green izquierdo a menudo se denomina tensor de deformación de Finger, en memoria de Josef Finger (1894).[4]

La Unión Internacional de Química Pura y Aplicada recomienda que este tensor se llame tensor de deformación de Green.[3]

Los invariantes de 𝐁 también se utilizan en las expresiones para caracterizar la función densidad de energía de deformación. Los invariantes usuales se definen como

I1:=tr(𝐁)=Bii=λ12+λ22+λ32I2:=12[(tr𝐁)2tr(𝐁2)]=12(Bii2BjkBkj)=λ12λ22+λ22λ32+λ32λ12I3:=det𝐁=J2=λ12λ22λ32

donde J:=det𝐅 es el determinante del gradiente de deformación.

Para materiales comprimibles, se utiliza un conjunto de invariantes ligeramente diferente:

(I¯1:=J2/3I1;I¯2:=J4/3I2;J1).

Tensor de deformación de Piola (tensor de deformación de Cauchy)

A principios de 1828,[5] Augustin Louis Cauchy introdujo un tensor de deformación definido como el inverso del tensor de deformación izquierdo de Cauchy-Green, 𝐁1. Este tensor también ha sido llamado tensor de deformación de Piola por la IUPAC[3] y tensor de Finger[6] en la literatura de reología y dinámica de fluidos.

𝐜=𝐁1=𝐅T𝐅1cij=XKxiXKxj

Representación espectral

Si hay tres deformaciones principales distintas λi, las descomposiciones espectrales de 𝐂 y 𝐁 vienen dadas por

𝐂=i=13λi2𝐍i𝐍i𝐁=i=13λi2𝐧i𝐧i

Además,

𝐔=i=13λi𝐍i𝐍i;𝐕=i=13λi𝐧i𝐧i
𝐑=i=13𝐧i𝐍i;𝐅=i=13λi𝐧i𝐍i

y se observa que

𝐕=𝐑𝐔𝐑T=i=13λi𝐑(𝐍i𝐍i)𝐑T=i=13λi(𝐑𝐍i)(𝐑𝐍i)

Por lo tanto, la unicidad de la descomposición espectral también implica que 𝐧i=𝐑𝐍i. El estiramiento izquierdo (𝐕) también se llama tensor de estiramiento espacial, mientras que el estiramiento derecho (𝐔) se llama tensor de estiramiento material.

El efecto de 𝐅 actuando sobre 𝐍i es estirar el vector λi y rotarlo a la nueva orientación 𝐧i, es decir,

𝐅𝐍i=λi(𝐑𝐍i)=λi𝐧i

En un sentido similar,

𝐅T𝐍i=1λi𝐧i;𝐅T𝐧i=λi𝐍i;𝐅1𝐧i=1λi𝐍i.

Ejemplos

Extensión uniaxial de un material incompresible.
Este es el caso en el que una muestra se estira en una dirección con un relación de estiramiento de α=α𝟏. Si el volumen permanece constante, la contracción en las otras dos direcciones es tal que α𝟏α𝟐α𝟑=𝟏 o α𝟐=α𝟑=α𝟎.𝟓. Entonces:
𝐅=[α000α0.5000α0.5] :𝐁=𝐂=[α2000α1000α1]
Cortante simple
𝐅=[1γ0010001] :𝐁=[1+γ2γ0γ10001] :𝐂=[1γ0γ1+γ20001]
Rotación del cuerpo rígido
𝐅=[cosθsinθ0sinθcosθ0001] :𝐁=𝐂=[100010001]=𝟏

Derivadas de deformación

Las derivadas del estiramiento con respecto al tensor de deformación de Cauchy-Green derecho se utilizan para caracterizar las relaciones tensión-deformación de muchos sólidos, particularmente los materiales hiperelásticos, según las expresiones siguientes:

λi𝐂=12λi𝐍i𝐍i=12λi𝐑T(𝐧i𝐧i)𝐑;i=1,2,3

de donde se deduce que

𝐂:(𝐍i𝐍i)=λi2;𝐂𝐂=𝖨(s);𝖨(s):(𝐍i𝐍i)=𝐍i𝐍i.

Interpretación física de los tensores de deformación

Sea 𝐗=Xi𝑬i un sistema de coordenadas cartesiano definido en el cuerpo no deformado y sea 𝐱=xi𝑬i otro sistema definido en el cuerpo deformado. A continuación, se parametriza una curva 𝐗(s) en el cuerpo no deformado usando s[0,1]. Su imagen en el cuerpo deformado es 𝐱(𝐗(s)).

La longitud no deformada de la curva viene dada por

lX=01|d𝐗ds|ds=01d𝐗dsd𝐗dsds=01d𝐗ds𝑰d𝐗dsds

Después de la deformación, la longitud pasa a ser

lx=01|d𝐱ds|ds=01d𝐱dsd𝐱dsds=01(d𝐱d𝐗d𝐗ds)(d𝐱d𝐗d𝐗ds)ds=01d𝐗ds[(d𝐱d𝐗)Td𝐱d𝐗]d𝐗dsds

Teniendo en cuenta que el tensor de deformación de Cauchy-Green derecho se define como

𝑪:=𝑭T𝑭=(d𝐱d𝐗)Td𝐱d𝐗

entonces

lx=01d𝐗ds𝑪d𝐗dsds

lo que indica que los cambios de longitud se caracterizan mediante 𝑪.

Tensores de deformación finita

El concepto de deformación se utiliza para evaluar en qué medida un desplazamiento dado difiere localmente de un desplazamiento de un cuerpo rígido.[1][7][8] Una de esas deformaciones para el caso de grandes deformaciones es el tensor de deformaciones finitas lagrangiano, también llamado tensor de deformaciones de Green lagrangiano o tensor de deformaciones de Green-St-Venant, definido como

𝐄=12(𝐂𝐈)EKL=12(xjXKxjXLδKL)

o en función del tensor de gradiente de desplazamiento

𝐄=12[(𝐗𝐮)T+𝐗𝐮+(𝐗𝐮)T𝐗𝐮]

o

EKL=12(uKXL+uLXK+uMXKuMXL)

El tensor de deformación de Green lagrangiano es una medida de cuánto difiere 𝐂 de 𝐈.

El tensor de deformación finita de Euler, o el tensor de deformación finita de Almansi euleriano, con referencia a la configuración deformada (es decir, según su descripción euleriana) se define como

𝐞=12(𝐈𝐜)=12(𝐈𝐁1)ers=12(δrsXMxrXMxs)

o en función de los gradientes de desplazamiento se tiene que

eij=12(uixj+ujxiukxiukxj)

Plantilla:Demostración

Familia de Seth-Hill de tensores de deformación generalizados

B. R. Seth, del Instituto Indio de Tecnología de Kharagpur, fue el primero en demostrar que los tensores de deformación de Green y de Almansi son casos especiales de una medida de deformación más general.[9][10] La idea fue ampliada aún más por Rodney Hill en 1968.[11] La familia de medidas de deformación de Seth-Hill (también llamadas tensores de Doyle-Ericksen)[12] se puede expresar como

𝐄(m)=12m(𝐔2m𝐈)=12m[𝐂m𝐈]

Para diferentes valores de m se tiene que:

  • El tensor de deformación de Green lagrangiano 𝐄(1)=12(𝐔2𝐈)=12(𝐂𝐈)
  • El tensor de deformación de Biot 𝐄(1/2)=(𝐔𝐈)=𝐂1/2𝐈
  • La deformación logarítmica, deformación natural, deformación verdadera o deformación de Hencky 𝐄(0)=ln𝐔=12ln𝐂
  • La deformación de Almansi 𝐄(1)=12[𝐈𝐔2]

La aproximación de segundo orden de estos tensores es

𝐄(m)=ε+12(𝐮)T𝐮(1m)εTε

donde ε es el tensor de deformación infinitesimal.

Son admisibles muchas otras definiciones diferentes de tensores 𝐄, siempre que todas cumplan las condiciones siguientes:[13]

  • 𝐄 desaparece para todos los movimientos de cuerpos rígidos.
  • La dependencia de 𝐄 del tensor de gradiente de desplazamiento 𝐮 es continua, continuamente diferenciable y monótona
  • También se impone que 𝐄 se reduzca al tensor de deformación infinitesimal ε como la norma |𝐮|0

Un ejemplo es el conjunto de tensores

𝐄(n)=(𝐔n𝐔n)/2n

que no pertenecen a la clase de Seth-Hill, pero tienen la misma aproximación de segundo orden que las medidas de Seth-Hill en m=0 para cualquier valor de n.[14]

Interpretación física del tensor de deformación finita

Las componentes diagonales EKL del tensor de deformación finita lagrangiano están relacionados con la deformación normal, como por ejemplo en

E11=e(𝐈1)+12e(𝐈1)2

donde e(𝐈1) es la deformación normal o de ingeniería en la dirección 𝐈1.

Los componentes fuera de la diagonal EKL del tensor de deformación finita lagrangiano están relacionados con la deformación cortante, como por ejemplo en

E12=122E11+12E22+1sinϕ12

donde ϕ12 es el cambio en el ángulo entre dos elementos lineales que originalmente eran perpendiculares a las direcciones 𝐈1 y 𝐈2, respectivamente.

En determinadas circunstancias, es decir, desplazamientos pequeños y tasas de desplazamiento pequeñas, las componentes del tensor de deformación finito lagrangiano pueden aproximarse mediante las componentes del tensor de deformación infinitesimal.

Plantilla:Demostración

Condiciones de compatibilidad

Plantilla:AP

El problema de la compatibilidad en la mecánica continua implica la determinación de campos continuos de un solo valor permisibles en los cuerpos. Estas condiciones permitidas dejan el cuerpo sin espacios ni superposiciones no físicas después de una deformación. La mayoría de estas condiciones se aplican a cuerpos simplemente conexos. Se requieren condiciones adicionales para los límites internos en el caso de múltiples cuerpos conexos.

Compatibilidad del gradiente de deformación

La condición necesaria y suficiente para la existencia de un campo 𝑭 compatible sobre un cuerpo simplemente conexo es que

×𝑭=0

Compatibilidad del tensor de deformación de Cauchy-Green derecho

Las condiciones necesarias y suficientes para la existencia de un campo 𝑪 compatible sobre un cuerpo simplemente conexo son

Rαβργ:=Xρ[(X)Γαβγ]Xβ[(X)Γαργ]+(X)Γμργ(X)Γαβμ(X)Γμβγ(X)Γαρμ=0

Se puede demostrar que estos son los componentes mixtos del tensor de curvatura. Por lo tanto, las condiciones necesarias para la compatibilidad 𝑪 son que la curvatura de Riemann-Christoffel de la deformación sea cero.

Compatibilidad del tensor de deformación de Cauchy-Green izquierdo

Amit Acharya dedujo las condiciones generales de suficiencia para el tensor de deformación izquierdo de Cauchy-Green en tres dimensiones.[15] Janet Blume encontró las condiciones de compatibilidad para campos 𝑩 bidimensionales.[16]

Véase también

Referencias

Plantilla:Listaref

Lecturas adicionales

Enlaces externos

Plantilla:Control de autoridades

  1. 1,0 1,1 Plantilla:Cite book
  2. A. Yavari, J.E. Marsden, and M. Ortiz, On spatial and material covariant balance laws in elasticity, Journal of Mathematical Physics, 47, 2006, 042903; pp. 1–53.
  3. 3,0 3,1 3,2 3,3 3,4 Plantilla:Cite journal
  4. Eduardo N. Dvorkin, Marcela B. Goldschmit, 2006 Nonlinear Continua, p. 25, Springer Plantilla:ISBN.
  5. Jirásek,Milan; Bažant, Z. P. (2002) Inelastic analysis of structures, Wiley, p. 463 Plantilla:ISBN
  6. J. N. Reddy, David K. Gartling (2000) The finite element method in heat transfer and fluid dynamics, p. 317, CRC Press Plantilla:ISBN.
  7. Plantilla:Cite book
  8. Plantilla:Cite journal
  9. Plantilla:Citation
  10. Plantilla:Citation
  11. Plantilla:Citation
  12. T.C. Doyle and J.L. Eriksen (1956). "Non-linear elasticity." Advances in Applied Mechanics 4, 53–115.
  13. Z.P. Bažant and L. Cedolin (1991). Stability of Structures. Elastic, Inelastic, Fracture and Damage Theories. Oxford Univ. Press, New York (2nd ed. Dover Publ., New York 2003; 3rd ed., World Scientific 2010).
  14. Z.P. Bažant (1998). "Easy-to-compute tensors with symmetric inverse approximating Hencky finite strain and its rate." Journal of Materials of Technology ASME, 120 (April), 131–136.
  15. Plantilla:Cite journal
  16. Plantilla:Cite journal