Articulo de referencia

Tasa de estrés objetiva

Predicciones a partir de tres tasas de tensión objetivas bajo cizallamiento En mecánica de medios continuos , las tasas de tensión objetivas son derivadas temporales de la tensi...

Predicciones a partir de tres tasas de tensión objetivas bajo cizallamiento

En mecánica de medios continuos , las tasas de tensión objetivas son derivadas temporales de la tensión que no dependen del sistema de referencia . [ 1 ] Muchas ecuaciones constitutivas se diseñan en forma de una relación entre una tasa de tensión y una tasa de deformación (o el tensor de velocidad de deformación ). La respuesta mecánica de un material no debe depender del sistema de referencia. En otras palabras, las ecuaciones constitutivas de los materiales deben ser independientes del sistema de referencia (objetivas) . Si las medidas de tensión y deformación son cantidades materiales , la objetividad se satisface automáticamente. Sin embargo, si las cantidades son espaciales , la objetividad de la tasa de tensión no está garantizada, incluso si la tasa de deformación es objetiva.

En mecánica de medios continuos existen numerosas tasas de tensión objetivas, todas las cuales pueden demostrarse como formas especiales de derivadas de Lie . Algunas de las tasas de tensión objetivas más utilizadas son:

  1. la tasa de Truesdell del tensor de tensiones de Cauchy ,
  2. la tasa de Green-Naghdi de la tensión de Cauchy y
  3. la tasa de Zaremba-Jaumann de la tensión de Cauchy. [ 2 ]

La figura adjunta muestra el desempeño de varias tasas objetivas en un ensayo de corte simple donde el modelo de material es hipoelástico con módulos elásticos constantes . La relación entre el esfuerzo cortante y el desplazamiento se grafica en función del tiempo. Se utilizan los mismos módulos con las tres tasas de esfuerzo objetivas. Claramente se observan oscilaciones espurias para la tasa de esfuerzo de Zaremba-Jaumann. [ 3 ] Esto no se debe a que una tasa sea mejor que otra, sino a que es un uso incorrecto de los modelos de material usar las mismas constantes con diferentes tasas objetivas. [ 4 ] Por esta razón, una tendencia reciente ha sido evitar por completo las tasas de esfuerzo objetivas siempre que sea posible.

No objetividad de la derivada temporal de la tensión de Cauchy

Bajo rotaciones de cuerpo rígido (Q{\displaystyle {\boldsymbol {Q}}}), el tensor de tensiones de Cauchyσ{\displaystyle {\boldsymbol {\sigma }}}se transforma como σr=QσQT ;  QQT=1{\displaystyle {\boldsymbol {\sigma }}_{r}={\boldsymbol {Q}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {Q}}^{T}~;~~{\boldsymbol {Q}}\cdot {\boldsymbol {Q}}^{T}={\boldsymbol {\mathit {1}}}} Desdeσ{\displaystyle {\boldsymbol {\sigma }}}es una magnitud espacial y la transformación sigue las reglas de las transformaciones tensoriales ,σ{\displaystyle {\boldsymbol {\sigma }}}es objetivo. Sin embargo, ddt(σr)=σ˙r=Q˙σQT+Qσ˙QT+QσQ˙TQσ˙QT.{\displaystyle {\cfrac {d}{dt}}({\boldsymbol {\sigma }}_{r})={\dot {\boldsymbol {\sigma }}}_{r}={\dot {\boldsymbol {Q}}}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {Q}}^{T}+{\boldsymbol {Q}}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {Q}}^{T}+{\boldsymbol {Q}}\cdot {\boldsymbol {\sigma }}\cdot {\dot {\boldsymbol {Q}}}^{T}\neq {\boldsymbol {Q}}\cdot {\dot {\boldsymbol {\sigma }}}\cdot {\boldsymbol {Q}}^{T}\,.} Por lo tanto, la tasa de tensión no es objetiva a menos que la velocidad de rotación sea cero, es decirQ{\displaystyle {\boldsymbol {Q}}}es constante.

Figura 1. Elemento de material sin deformar y deformado, y un cubo elemental recortado del elemento deformado.

Para comprender físicamente lo anterior, considere la situación que se muestra en la Figura 1. En la figura, los componentes del tensor de tensiones de Cauchy (o verdadero) se denotan mediante los símbolosSij{\displaystyle S_{ij}}Este tensor, que describe las fuerzas sobre un pequeño elemento material imaginado como si se hubiera extraído del material deformado, no es objetivo para grandes deformaciones porque varía con las rotaciones de cuerpo rígido del material. Los puntos materiales deben caracterizarse por sus coordenadas lagrangianas iniciales.incógnitai{\displaystyle X_{i}}En consecuencia, es necesario introducir la denominada tasa de estrés objetiva.Sij{\displaystyle {\overset {\circ }{S}}_{ij}}o el incremento correspondienteΔSij=SijΔt{\displaystyle \Delta S_{ij}={\overset {\circ }{S}}_{ij}\Delta t}La objetividad es necesaria paraSij{\displaystyle {\overset {\circ }{S}}_{ij}}estar funcionalmente relacionado con la deformación del elemento. Significa queSij{\displaystyle {\overset {\circ }{S}}_{ij}}debe ser invariante con respecto a las transformaciones de coordenadas, en particular las rotaciones de cuerpo rígido, y debe caracterizar el estado del mismo elemento material a medida que se deforma.

La tasa de estrés objetiva se puede obtener de dos maneras:

  • mediante transformaciones de coordenadas tensoriales, [ 5 ] que es la forma estándar en los libros de texto de elementos finitos [ 6 ]
  • Variacionalmente, a partir de la densidad de energía de deformación en el material expresada en términos del tensor de deformación (que es objetivo por definición) [ 7 ] [ 8 ]

Si bien el primer método es instructivo y proporciona información geométrica útil, el segundo es matemáticamente más corto y tiene la ventaja adicional de garantizar automáticamente la conservación de la energía , es decir, garantizar que el trabajo de segundo orden del tensor de incremento de tensión sobre el tensor de incremento de deformación sea correcto (requisito de conjugación de trabajo).

Tasa de estrés de Truesdell de la tensión de Cauchy

La relación entre la tensión de Cauchy y la tensión de segundo orden de Piola se denomina transformación de Piola . Esta transformación puede escribirse en términos del retroceso deσ{\displaystyle {\boldsymbol {\sigma }}}o el impulso hacia adelante deS{\displaystyle {\boldsymbol {S}}}como S=J ϕ[σ] ;  σ=J1 ϕ[S]{\displaystyle {\boldsymbol {S}}=J~\phi ^{*}[{\boldsymbol {\sigma }}]~;~~{\boldsymbol {\sigma }}=J^{-1}~\phi _{*}[{\boldsymbol {S}}]}

La tasa de Truesdell de la tensión de Cauchy es la transformación de Piola de la derivada temporal material de la segunda tensión PK. Por lo tanto, definimos σ=J1 ϕ[S˙]{\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~\phi _{*}[{\dot {\boldsymbol {S}}}]}

Ampliado, esto significa que

σ=J1 FS˙FT=J1 F[ddt(J F1σFT)]FT=J1 Lφ[τ]{\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}\cdot {\boldsymbol {F}}^{T}=J^{-1}~{\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left(J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}=J^{-1}~{\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]} donde el estrés de Kirchhoffτ=J σ{\displaystyle {\boldsymbol {\tau }}=J~{\boldsymbol {\sigma }}}y la derivada de Lie de la tensión de Kirchhoff es Lφ[τ]=F[ddt(F1τFT)]FT .{\displaystyle {\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\tau }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}~.}

Esta expresión se puede simplificar a la conocida expresión para la tasa de Truesdell de la tensión de Cauchy.

Se puede demostrar que la tasa de Truesdell es objetiva.

Tasa de Truesdell de la tensión de Kirchhoff

La tasa de Truesdell de la tensión de Kirchhoff se puede obtener observando que S=ϕ[τ] ;  τ=ϕ[S]{\displaystyle {\boldsymbol {S}}=\phi ^{*}[{\boldsymbol {\tau }}]~;~~{\boldsymbol {\tau }}=\phi _{*}[{\boldsymbol {S}}]} y definiendo τ=ϕ[S˙]{\displaystyle {\overset {\circ }{\boldsymbol {\tau }}}=\phi _{*}[{\dot {\boldsymbol {S}}}]} Ampliado, esto significa que τ=FS˙FT=F[ddt(F1τFT)]FT=Lφ[τ]{\displaystyle {\overset {\circ }{\boldsymbol {\tau }}}={\boldsymbol {F}}\cdot {\dot {\boldsymbol {S}}}\cdot {\boldsymbol {F}}^{T}={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\tau }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}={\mathcal {L}}_{\varphi }[{\boldsymbol {\tau }}]} Por lo tanto, la derivada de Lie deτ{\displaystyle {\boldsymbol {\tau }}}es lo mismo que la tasa de Truesdell de la tensión de Kirchhoff .

Siguiendo el mismo proceso que para la tensión de Cauchy anterior, podemos demostrar que

Tasa de estrés de Cauchy de Green-Naghdi

Esta es una forma especial de la derivada de Lie (o la tasa de Truesdell de la tensión de Cauchy). Recordemos que la tasa de Truesdell de la tensión de Cauchy viene dada por σ=J1 F[ddt(J F1σFT)]FT .{\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}=J^{-1}~{\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left(J~{\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}~.} Del teorema de descomposición polar tenemos F=RU{\displaystyle {\boldsymbol {F}}={\boldsymbol {R}}\cdot {\boldsymbol {U}}} dóndeR{\displaystyle {\boldsymbol {R}}}es el tensor de rotación ortogonal (R1=RT{\displaystyle {\boldsymbol {R}}^{-1}={\boldsymbol {R}}^{T}}) yU{\displaystyle {\boldsymbol {U}}}es el estiramiento derecho simétrico, definido positivo.

Si asumimos queU=1{\displaystyle {\boldsymbol {U}}={\boldsymbol {\mathit {1}}}}obtenemosF=R{\displaystyle {\boldsymbol {F}}={\boldsymbol {R}}}Además, como no hay estiramientoJ=1{\displaystyle J=1}y tenemosτ=σ{\displaystyle {\boldsymbol {\tau }}={\boldsymbol {\sigma }}}. Tenga en cuenta que esto no significa que no haya estiramiento en el cuerpo real; esta simplificación es solo para definir una tasa de estrés objetiva. Por lo tanto, σ=R[ddt(R1σRT)]RT=R[ddt(RTσR)]RT{\displaystyle {\overset {\circ }{\boldsymbol {\sigma }}}={\boldsymbol {R}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {R}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {R}}^{-T}\right)\right]\cdot {\boldsymbol {R}}^{T}={\boldsymbol {R}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {R}}^{T}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {R}}\right)\right]\cdot {\boldsymbol {R}}^{T}} Podemos demostrar que esta expresión se puede simplificar a la forma comúnmente utilizada de la tasa de Green-Naghdi .

La tasa de Green-Naghdi de la tensión de Kirchhoff también tiene la forma ya que no se tiene en cuenta el estiramiento, es decir, τ=τ˙+τΩΩτ{\displaystyle {\overset {\square }{\boldsymbol {\tau }}}={\dot {\boldsymbol {\tau }}}+{\boldsymbol {\tau }}\cdot {\boldsymbol {\Omega }}-{\boldsymbol {\Omega }}\cdot {\boldsymbol {\tau }}}

Tasa de Zaremba-Jaumann de la tensión de Cauchy

La tasa de Zaremba-Jaumann de la tensión de Cauchy es una especialización adicional de la derivada de Lie (tasa de Truesdell). Esta tasa tiene la forma

La tasa de Zaremba-Jaumann se utiliza ampliamente en cálculos principalmente por dos razones.

  1. Su implementación es relativamente sencilla.
  2. Esto da lugar a módulos tangentes simétricos.

Recordemos que el tensor de espínw{\displaystyle {\boldsymbol {w}}}(la parte sesgada del gradiente de velocidad) se puede expresar como w=R˙RT+12 R(U˙U1U1U˙)RT{\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}+{\frac {1}{2}}~{\boldsymbol {R}}\cdot ({\dot {\boldsymbol {U}}}\cdot {\boldsymbol {U}}^{-1}-{\boldsymbol {U}}^{-1}\cdot {\dot {\boldsymbol {U}}})\cdot {\boldsymbol {R}}^{T}} Por lo tanto, para el movimiento de cuerpo rígido puro w=R˙RT=Ω{\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}={\boldsymbol {\Omega }}} Alternativamente, podemos considerar el caso de carga proporcional cuando las direcciones principales de deformación permanecen constantes. Un ejemplo de esta situación es la carga axial de una barra cilíndrica. En esa situación, dado que U=[λincógnitaλYλZ]{\displaystyle {\boldsymbol {U}}={\begin{bmatrix}\lambda _{X}\\&\lambda _{Y}\\&&\lambda _{Z}\end{bmatrix}}} tenemos U˙=[λ˙incógnitaλ˙Yλ˙Z]{\displaystyle {\dot {\boldsymbol {U}}}={\begin{bmatrix}{\dot {\lambda }}_{X}\\&{\dot {\lambda }}_{Y}\\&&{\dot {\lambda }}_{Z}\end{bmatrix}}} También, U1=[1/λincógnita1/λY1/λZ]{\displaystyle {\boldsymbol {U}}^{-1}={\begin{bmatrix}1/\lambda _{X}\\&1/\lambda _{Y}\\&&1/\lambda _{Z}\end{bmatrix}}} del estrés de Cauchy. Por lo tanto, U˙U1=[λ˙incógnita/λincógnitaλ˙Y/λYλ˙Z/λZ]=U1U˙{\displaystyle {\dot {\boldsymbol {U}}}\cdot {\boldsymbol {U}}^{-1}={\begin{bmatrix}{\dot {\lambda }}_{X}/\lambda _{X}\\&{\dot {\lambda }}_{Y}/\lambda _{Y}\\&&{\dot {\lambda }}_{Z}/\lambda _{Z}\end{bmatrix}}=U^{-1}{\dot {U}}} Esto una vez más da w=R˙RT=Ω{\displaystyle {\boldsymbol {w}}={\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}={\boldsymbol {\Omega }}} En general, si aproximamos wR˙RT{\displaystyle {\boldsymbol {w}}\approx {\dot {\boldsymbol {R}}}\cdot {\boldsymbol {R}}^{T}} La tasa de Green-Naghdi se convierte en la tasa de Zaremba-Jaumann de la tensión de Cauchy. σ=σ˙+σwwσ{\displaystyle {\overset {\triangle }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {w}}-{\boldsymbol {w}}\cdot {\boldsymbol {\sigma }}}

Otras tasas de estrés objetivas

Puede existir una variedad infinita de tasas de estrés objetivas. Una de ellas es la tasa de estrés de Oldroyd.σ=Lφ[σ]=F[ddt(F1σFT)]FT{\displaystyle {\overset {\triangledown }{\boldsymbol {\sigma }}}={\mathcal {L}}_{\varphi }[{\boldsymbol {\sigma }}]={\boldsymbol {F}}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{-1}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}^{-T}\right)\right]\cdot {\boldsymbol {F}}^{T}} En forma más simple, la tasa de Oldroyd viene dada por σ=σ˙lσσlT{\displaystyle {\overset {\triangledown }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}-{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}}

Si se asume que la configuración actual es la configuración de referencia, entonces las operaciones de retroceso y avance se pueden realizar utilizandoFT{\displaystyle {\boldsymbol {F}}^{T}}yFT{\displaystyle {\boldsymbol {F}}^{-T}}respectivamente. La derivada de Lie de la tensión de Cauchy se denomina entonces tasa de tensión convectiva.σ=FT[ddt(FTσF)]F1{\displaystyle {\overset {\diamond }{\boldsymbol {\sigma }}}={\boldsymbol {F}}^{-T}\cdot \left[{\cfrac {d}{dt}}\left({\boldsymbol {F}}^{T}\cdot {\boldsymbol {\sigma }}\cdot {\boldsymbol {F}}\right)\right]\cdot {\boldsymbol {F}}^{-1}} En forma más simple, la tasa de convección viene dada por σ=σ˙+lσ+σlT{\displaystyle {\overset {\diamond }{\boldsymbol {\sigma }}}={\dot {\boldsymbol {\sigma }}}+{\boldsymbol {l}}\cdot {\boldsymbol {\sigma }}+{\boldsymbol {\sigma }}\cdot {\boldsymbol {l}}^{T}}

Tasas de tensión objetivas en la inelasticidad de deformación finita

Muchos materiales experimentan deformaciones inelásticas causadas por plasticidad y daños. Estos comportamientos no pueden describirse mediante un potencial. Además, a menudo no existe memoria del estado inicial, especialmente cuando se producen grandes deformaciones. En estos casos, la relación constitutiva se define generalmente de forma incremental para facilitar el cálculo de tensiones y deformaciones.

El procedimiento de carga incremental

Para un paso de carga suficientemente pequeño, la deformación del material se puede caracterizar mediante el tensor de incremento de deformación pequeño (o linealizado).mi=12[+()T]miij=12(i,j+j,i){\displaystyle {\boldsymbol {e}}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}\right]\quad \equiv \quad e_{ij}={\tfrac {1}{2}}(u_{i,j}+u_{j,i})} dónde{\displaystyle \mathbf {u} }es el incremento de desplazamiento de los puntos del continuo. La derivada temporal mit=mi˙=12[v+(v)T]mi˙ij=12(vi,j+vj,i){\displaystyle {\frac {\partial {\boldsymbol {e}}}{\partial t}}={\dot {\boldsymbol {e}}}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {v} +({\boldsymbol {\nabla }}\mathbf {v} )^{T}\right]\quad \equiv \quad {\dot {e}}_{ij}={\tfrac {1}{2}}(v_{i,j}+v_{j,i})} es el tensor de velocidad de deformación (también llamado deformación de velocidad) yv=˙{\displaystyle \mathbf {v} ={\dot {\mathbf {u} }}}es la velocidad o tasa de desplazamiento del punto material. Para deformaciones finitas, se pueden utilizar medidas de la familia de Seth-Hill (también llamadas tensores de Doyle-Ericksen):mi(metro)=12metro(U2metroI){\displaystyle \mathbf {E} _{(m)}={\frac {1}{2m}}(\mathbf {U} ^{2m}-\mathbf {I} )} dóndeU{\displaystyle \mathbf {U} }es el estiramiento correcto. Una aproximación de segundo orden de estos tensores es mi(metro)mi+12()T(1metro)mimi{\displaystyle \mathbf {E} _{(m)}\approx {\boldsymbol {e}}+{\tfrac {1}{2}}(\nabla \mathbf {u} )^{T}\cdot \nabla \mathbf {u} -(1-m){\boldsymbol {e}}\cdot {\boldsymbol {e}}}

Tasas de estrés objetivas consistentes con la energía

Consideremos un elemento material de volumen inicial unitario, que parte de un estado inicial bajo tensión de Cauchy (o tensión real) inicial.σ0{\displaystyle {\boldsymbol {\sigma }}_{0}}y dejarσ{\displaystyle {\boldsymbol {\sigma }}}Sea la tensión de Cauchy en la configuración final.W{\displaystyle W}Sea el trabajo realizado (por unidad de volumen inicial) por las fuerzas internas durante una deformación incremental desde este estado inicial. Entonces, la variaciónδW{\displaystyle \delta W}corresponde a la variación en el trabajo realizado debido a una variación en el desplazamientoδ{\displaystyle \delta \mathbf {u} }La variación del desplazamiento debe satisfacer las condiciones de contorno del desplazamiento.

DejarS(metro){\displaystyle {\boldsymbol {S}}_{(m)}}Sea un tensor de tensión objetivo en la configuración inicial. Defina el incremento de tensión con respecto a la configuración inicial comoS=S(metro)σ0{\displaystyle {\boldsymbol {S}}={\boldsymbol {S}}_{(m)}-{\boldsymbol {\sigma }}_{0}}Alternativamente, siPAG{\displaystyle {\boldsymbol {P}}}es la primera tensión de Piola-Kirchhoff asimétrica referida a la configuración inicial, el incremento de tensión se puede expresar comoT=PAGσ0{\displaystyle {\boldsymbol {T}}={\boldsymbol {P}}-{\boldsymbol {\sigma }}_{0}}.

Variación del trabajo realizado

Entonces, la variación en el trabajo realizado se puede expresar como δW=S(metro):δmi(metro)=PAG:δ{\displaystyle \delta W={\boldsymbol {S}}_{(m)}:\delta {\boldsymbol {E}}_{(m)}={\boldsymbol {P}}:\delta \nabla \mathbf {u} } donde la medida de deformación finitami(metro){\displaystyle {\boldsymbol {E}}_{(m)}}es energía conjugada a la medida de estrésσ(metro){\displaystyle {\boldsymbol {\sigma }}^{(m)}}. Expandido, δW=(S+σ0):δmi(metro)=(T+σ0):δ.{\displaystyle \delta W=\left({\boldsymbol {S}}+{\boldsymbol {\sigma }}_{0}\right):\delta {\boldsymbol {E}}_{(m)}=\left({\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}\right):\delta \nabla \mathbf {u} \,.} La objetividad del tensor de tensionesS(metro){\displaystyle {\boldsymbol {S}}_{(m)}}está garantizado por su transformación en un tensor de segundo orden bajo rotaciones de coordenadas (lo que hace que las tensiones principales sean independientes de las rotaciones de coordenadas) y por la corrección deS(metro):δmi(metro){\displaystyle {\boldsymbol {S}}_{(m)}:\delta {\boldsymbol {E}}_{(m)}}como una expresión de energía de segundo orden.

A partir de la simetría de la tensión de Cauchy, tenemos σ0:δ=σ0:δmi.{\displaystyle {\boldsymbol {\sigma }}_{0}:\delta \nabla \mathbf {u} ={\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {e}}\,.} Para pequeñas variaciones en la deformación, utilizando la aproximación S:δmi(metro)S:δ{\displaystyle {\boldsymbol {S}}:\delta {\boldsymbol {E}}_{(m)}\approx {\boldsymbol {S}}:\delta \nabla \mathbf {u} } y las expansiones σ0:δmi(metro)=σ0:[mi(metro):δ] ,  σ0:δmi=σ0:[mi:δ]{\displaystyle {\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {E}}_{(m)}={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]~,~~{\boldsymbol {\sigma }}_{0}:\delta {\boldsymbol {e}}={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]} obtenemos la ecuación σ0:[mi(metro):δ]+S:δ=σ0:[mi:δ]+T:δ.{\displaystyle {\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {E}}_{(m)}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]+{\boldsymbol {S}}:\delta \nabla \mathbf {u} ={\boldsymbol {\sigma }}_{0}:\left[{\frac {\partial {\boldsymbol {e}}}{\partial \nabla \mathbf {u} }}:\delta \nabla \mathbf {u} \right]+{\boldsymbol {T}}:\delta \nabla \mathbf {u} \,.} Imponiendo la condición variacional de que la ecuación resultante debe ser válida para cualquier gradiente de deformación.δ{\displaystyle \delta \nabla \mathbf {u} }, tenemos [ 7 ]

También podemos escribir la ecuación anterior como

derivadas temporales

La tensión de Cauchy y la primera tensión de Piola-Kirchhoff están relacionadas por (ver Medidas de tensión ) σ=PAGFTJ1=(PAG+PAGT)J1.{\displaystyle {\boldsymbol {\sigma }}={\boldsymbol {P}}\cdot {\boldsymbol {F}}^{T}J^{-1}=({\boldsymbol {P}}+{\boldsymbol {P}}\cdot \nabla \mathbf {u} ^{T})J^{-1}\,.} Para pequeñas deformaciones incrementales, J11.{\displaystyle J^{-1}\approx 1-\nabla \cdot \mathbf {u} \,.} Por lo tanto, Δσ=σσ0(PAG+PAGT)(1)σ0.{\displaystyle \Delta {\boldsymbol {\sigma }}={\boldsymbol {\sigma }}-{\boldsymbol {\sigma }}_{0}\approx ({\boldsymbol {P}}+{\boldsymbol {P}}\cdot \nabla \mathbf {u} ^{T})(1-\nabla \cdot \mathbf {u} )-{\boldsymbol {\sigma }}_{0}\,.} SustituyendoT+σ0=PAG{\displaystyle {\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}={\boldsymbol {P}}}, Δσ[T+σ0+(T+σ0)T](1)σ0.{\displaystyle \Delta {\boldsymbol {\sigma }}\approx [{\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0}+({\boldsymbol {T}}+{\boldsymbol {\sigma }}_{0})\cdot \nabla \mathbf {u} ^{T}](1-\nabla \cdot \mathbf {u} )-{\boldsymbol {\sigma }}_{0}\,.} Para pequeños incrementos de tensiónT{\displaystyle {\boldsymbol {T}}}en relación con la tensión inicialσ0{\displaystyle {\boldsymbol {\sigma }}_{0}}, lo anterior se reduce a

De las ecuaciones (1) y (3) tenemos

Recuerda queS{\displaystyle {\boldsymbol {S}}}es un incremento de la medida del tensor de tensiónS(metro){\displaystyle {\boldsymbol {S}}_{(m)}}. Definir la tasa de estrés S=:S(metro)Δt{\displaystyle {\boldsymbol {S}}=:{\overset {\circ }{\boldsymbol {S}}}_{(m)}\Delta t} y observando que Δσ=σ˙Δt{\displaystyle \Delta {\boldsymbol {\sigma }}={\dot {\boldsymbol {\sigma }}}\Delta t} podemos escribir la ecuación (4) como

Tomando el límite enΔt0{\displaystyle \Delta t\rightarrow 0}y señalando queσ0=σ{\displaystyle {\boldsymbol {\sigma }}_{0}={\boldsymbol {\sigma }}}En este límite, se obtiene la siguiente expresión para la tasa de tensión objetiva asociada con la medida de deformación.mi(metro){\displaystyle {\boldsymbol {E}}_{(m)}}:

Aquíσ˙ij=σij/t{\displaystyle {\dot {\sigma }}_{ij}=\partial \sigma _{ij}/\partial t}= tasa material de tensión de Cauchy (es decir, la tasa en coordenadas lagrangianas del estado de tensión inicial).

tasas de tensión conjugada al trabajo

Una tasa para la cual no existe un tensor de deformación finito legítimo.mi(metro){\displaystyle {\boldsymbol {E}}_{(m)}}La asociación según la ecuación (6) es energéticamente inconsistente, es decir, su uso viola el balance energético (es decir, la primera ley de la termodinámica ).

Evaluar la ecuación (6) para generalmetro{\displaystyle m}y parametro=2{\displaystyle m=2}, se obtiene una expresión general para la tasa de estrés objetiva: [ 7 ] [ 8 ]

dóndeS(2){\displaystyle {\overset {\circ }{\boldsymbol {S}}}_{(2)}}es la tasa de estrés objetiva asociada con la deformación de Green-Lagrangiano (metro=2{\displaystyle m=2}).

En particular,

  • metro=2{\displaystyle m=2}proporciona la tasa de estrés de Truesdell
  • metro=0{\displaystyle m=0}proporciona la tasa de estrés de Kirchhoff de Zaremba-Jaumann
  • metro=1{\displaystyle m=1}proporciona la tasa de estrés de Biot

(Tenga en cuenta que m = 2 conduce a la fórmula de Engesser para la carga crítica en el pandeo por cortante, mientras que m = -2 conduce a la fórmula de Haringx , que puede dar cargas críticas que difieren en más del 100%).

tasas de tensión no conjugadas con el trabajo

Otras tasas, utilizadas en la mayoría de los códigos comerciales, que no son conjugadas por trabajo a ningún tensor de deformación finito son: [ 8 ]

  • La tasa de tensión de Cauchy de Zaremba-Jaumann, o corrotacional : difiere de la tasa de tensión de Kirchhoff de Zaremba-Jaumann al no considerar la tasa de cambio de volumen relativo del material. La falta de conjugación del trabajo no suele ser un problema grave, ya que este término es insignificante para muchos materiales y cero para materiales incompresibles (pero en la indentación de una placa sándwich con núcleo de espuma, esta tasa puede generar un error superior al 30 % en la fuerza de indentación).
  • La tasa de Cotter-Rivlin corresponde ametro=2{\displaystyle m=-2}pero de nuevo omite el término volumétrico.
  • La tasa de Green-Naghdi : Esta tasa de tensión objetiva no es conjugada con el trabajo de ningún tensor de deformación finito, no solo por la ausencia del término volumétrico, sino también porque la velocidad de rotación del material no es exactamente igual al tensor de giro . En la gran mayoría de las aplicaciones, los errores en el cálculo de la energía, causados ​​por estas diferencias, son despreciables. Sin embargo, cabe señalar que ya se ha demostrado un gran error de energía para un caso con deformaciones de corte y rotaciones superiores a aproximadamente 0,25. [ 9 ]
  • la tasa Oldroyd .

Tasas objetivas y derivados de Lie

Las tasas de estrés objetivas también podrían considerarse como las derivadas de Lie de varios tipos de tensor de estrés (es decir, los componentes covariantes, contravariantes y mixtos asociados del estrés de Cauchy) y sus combinaciones lineales. [ 10 ] La derivada de Lie no incluye el concepto de conjugación de trabajo.

Módulos de rigidez tangencial y sus transformaciones para lograr consistencia energética.

La relación tensión-deformación tangencial generalmente tiene la forma

dóndedoijkl(metro){\displaystyle C_{ijkl}^{(m)}}son los módulos tangenciales (componentes de un tensor de cuarto orden) asociados con el tensor de deformaciónϵij(metro){\displaystyle \epsilon _{ij}^{(m)}}Son diferentes para diferentes elecciones demetro{\displaystyle m}y se relacionan de la siguiente manera:

Del hecho de que la ecuación (7) debe cumplirse para cualquier gradiente de velocidadvk,l{\displaystyle v_{k,l}}, de ello se deduce que: [ 7 ]

dóndedoijkl(2){\displaystyle C_{ijkl}^{(2)}}son los módulos tangenciales asociados con la deformación Green-Lagrangiana (metro=2{\displaystyle m=2}), tomado como referencia,Sij{\displaystyle S_{ij}}= tensión de Cauchy actual, yδij{\displaystyle \delta _{ij}}= Delta de Kronecker (o tensor unitario).

La ecuación (8) se puede utilizar para convertir una tasa de tensión objetiva en otra. Dado queSijmi˙kk=(Sijδkl)δmikl{\displaystyle S_{ij}{\dot {e}}_{kk}=(S_{ij}\delta _{kl})\delta e_{kl}}, la transformación [ 7 ] [ 8 ]

puede corregir aún más la ausencia del términoSijvk,k{\displaystyle S_{ij}v_{k,k}}(tenga en cuenta que el términoSijδkmetro{\displaystyle S_{ij}\delta _{km}}No permite intercambiar subíndices.ij{\displaystyle ij}conkl{\displaystyle kl}, lo que significa que su ausencia rompe la simetría principal del tensor de módulos tangencialesdoijklnorteonortedoonortej{\displaystyle C_{ijkl}^{\mathrm {nonconj} }}).

A menudo se desarrollan grandes deformaciones cuando el comportamiento del material se vuelve no lineal, debido a la plasticidad o al daño. Entonces, la causa principal de la dependencia de la tensión de los módulos tangenciales es el comportamiento físico del material. Lo que significa la ecuación (8) es que la dependencia no lineal dedoijkl{\displaystyle C_{ijkl}}sobre el estrés debe ser diferente para diferentes tasas de estrés objetivas. Sin embargo, ninguna de ellas es fundamentalmente preferible, excepto si existe una tasa de estrés, unametro{\displaystyle m}, para los cuales los módulos pueden considerarse constantes.

Véase también

Referencias

  1. ME Gurtin, E. Fried y L. Anand (2010). "La mecánica y la termodinámica de los continuos". Cambridge University Press , (ver págs. 151, 242).
  2. Zaremba, "Sur une forme perfectée de la théorie de la relax" , Toro. Int. Acad. Ciencia. Cracovia , 1903.
  3. Dienes, J. (1979). "Sobre el análisis de la rotación y la tasa de tensión en cuerpos deformables". Acta Mechanica . Vol.  32. pág.  217.
  4. Brannon, RM (1998). "Advertencias sobre las medidas conjugadas de tensión y deformación para la elasticidad anisotrópica indiferenciada al marco". Acta Mechanica . Vol. 129. pp. 107–116 .  
  5. HD Hibbitt, PV Marçal y JR Rice (1970). "Una formulación de elementos finitos para problemas de grandes deformaciones y grandes desplazamientos". Intern. J. of Solids Structures , 6, 1069–1086.
  6. T. Belytschko, W. K. Liu y B. Moran (2000). Elementos finitos no lineales para continuos y estructuras. J. Wiley & Sons, Chichester, Reino Unido.
  7. 1 2 3 4 5 Z.P. Bažant (1971). "Un estudio de correlación de formulaciones de deformación incremental y estabilidad de cuerpos continuos". J. of Applied Mechanics ASME , 38(4), 919–928.
  8. 1 2 3 4 Z.P. Bažant y L. Cedolin (1991). Estabilidad de estructuras. Teorías elásticas, inelásticas, de fractura y de daño. Oxford Univ. Press, Nueva York (2.ª ed. Dover Publ., Nueva York 2003; 3.ª ed., World Scientific 2010).
  9. ZP Bažant y J. Vorel (2013). Error de conservación de energía debido al uso de la tasa de tensión objetiva de Green-Naghdi en códigos comerciales de elementos finitos y su compensación." ASME Journal of Applied Mechanics , 80(4).
  10. JE Marsden y TJR Hughes (1983). Fundamentos matemáticos de la elasticidad. Prentice Hall, Englewood Cliffs, NJ (pág. 100).
  • Objetividad en la mecánica clásica
  • Wikiversidad: Elementos finitos no lineales/Enfoque lagrangiano actualizado