Suavizado de esfuerzos internos en superficies en RFEM

Puede seleccionar diferentes opciones en el programa para mostrar resultados sin suavizar o suavizados. Estos se describen en el manual en línea de RFEM (en inglés).

• Manual en línea de RFEM 5

Las opciones de suavizado para evaluar los esfuerzos internos y tensiones de superficies también se explican en este artículo técnico.

• KB 001571 | Diagrama de esfuerzos internos/tensiones superficiales - Opciones de suavizado

El siguiente artículo utiliza un ejemplo para ilustrar cómo se usan los esfuerzos internos en los elementos finitos con el suavizado para determinar los valores. Además de la posición de un elemento finito, la teoría de placas (Mindlin o Kirchhoff) juega un papel en este proceso.

Ejemplo: Placa según Mindlin y Kirchhoff

1 Ejemplo

Se crea una placa con las dimensiones 4 m ⋅ 3 m como modelo bidimensional. La placa está articulada en ambos lados a las líneas de límite más cortas. El giro de las líneas de contorno alargadas está restringido por una coacción. Para evitar los efectos de la deformación transversal, el coeficiente de Poisson del material se establece en 0. Para comparar la determinación de los esfuerzos internos según las teorías de flexión de la placa de Kirchhoff y Mindlin, se selecciona un espesor de la placa de 80 cm. Por lo tanto, la relación d/L = 0.2 se da en la zona límite entre las dos teorías de flexión (ver FAQ 003158).

• FAQ | ¿Es suficiente con realizar el cálculo con los elementos de superficie o es mejor utilizar elemento sólidos?

Para simplificarlo, el tamaño de la malla de EF se establece en 1 m. Por lo tanto, hay 4 · 3 elementos finitos.

Una carga de p = 5 kN/m² actúa sobre la placa. El peso propio no se aplica en el caso de carga.

Fuerzas internas sin suavizar

Los resultados del núcleo de cálculo para los esfuerzos internos v_x y m_x se consideran para una sección longitudinal que se da en el centro de la placa. La distribución de esfuerzos internos sin suavizado elemento por elemento se puede representar con la opción de visualización "No continuo". Según las teorías de flexión de Mindlin y Kirchhoff, están disponibles los diagramas mostrados en la figura 02.

2 Resultados sin suavizar según Mindlin (izquierda) y Kirchhoff (derecha)

Suavizado estándar para los nudos internos

Para los nudos de EF que se encuentren dentro de una superficie, la media aritmética se forma primero a partir de los valores en nudo de los elementos finitos adyacentes. Para este enfoque, los elementos se deben encontrar en una superficie y en el mismo lado de una posible línea interna. El nudo no debe ser un nudo definido por el usuario en la superficie.

Para la ubicación de la discontinuidad en x = 1,00 m, estos valores son los siguientes internamente (el resultado de este paso no está disponible para la salida).

• Mindlin:
  $${\mathrm{v}}_{\mathrm{x}}\left(1\right) = \frac{7.5 · 2.5}{2} = 5$$

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{8.772 · 5.351}{2} = 7.061$$

  Suavizado estándar para nudos internos según Mindlin
• Kirchhoff:
  $${\mathrm{v}}_{\mathrm{x}}\left(1\right) = \frac{7.785 · 2.595}{2} = 5.19$$

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{7.954 · 7.904}{2} = 7.929$$

  Suavizado estándar para nudos internos según Kirchhoff

Suavizado para nudos de exteriores

En los bordes de una superficie (opción de suavizados "Continua dentro de las superficies") o un modelo ("Continua total"), no hay elementos finitos adyacentes que se puedan utilizar para promediar los valores en los nudos. Por lo tanto, se utiliza un enfoque diferente, que se realiza en dos pasos.

En el primer caso, se calculan los valores medios para aquellos nudos que no están ubicados en el borde de la superficie o del modelo. Los valores de los nudos en los bordes se calculan de tal forma que los valores originales permanecen en los centros de los elementos finitos. En el segundo paso, se determinan los valores medios de los nudos exteriores.

Se obtienen los siguientes valores para el nudo de borde en x = 0,00 m.

• Mindlin:
  $${\mathrm{v}}_{\mathrm{x}}\left(0\right) = 7.5 · (7.5 - 5) = 10$$

  $${\mathrm{m}}_{\mathrm{x}}(0.5) = \frac{5.351 · 2.982}{2} = 4.167$$

  $${\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 · (4.167 - 7.061) = 1.272$$

  Suavizado para nudos de borde según Mindlin
• Kirchhoff:
  $${\mathrm{v}}_{\mathrm{x}}\left(0\right) = 7.785 · (7.785 - 5.19) = 10.381$$

  $${\mathrm{m}}_{\mathrm{x}}(0.5) = \frac{7.954 · 0.379}{2} = 4.167$$

  $${\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 · (4.167 - 7.929) = 0.404$$

  Suavizado para nudos de borde según Kirchhoff

Esfuerzos cortantes

Según Mindlin, los esfuerzos cortantes se calculan como la primera derivada de la deformación.

Estos valores están determinados por el cálculo del núcleo y se usan directamente. El suavizado de los esfuerzos cortantes se lleva a cabo como se ha descrito anteriormente, dependiendo de la posición del nudo del FE en la superficie.

En la teoría de flexión según Kirchhoff, los esfuerzos cortantes se calculan como una tercera derivada de la deformación.

Esta determinación del valor utilizando la tercera derivada da como resultado una pérdida considerable de precisión. Por esta razón, los esfuerzos cortantes del cálculo del núcleo no se utilizan, sino que se determinan de las derivaciones de los momentos mediante un planteamiento mejorado.

En las ecuaciones, los momentos m_x, m_y y m_xy representan los valores suavizados que se han determinado según los métodos descritos anteriormente. Esto da como resultado unos valores más precisos para los esfuerzos cortantes que realizar el cálculo del núcleo.

Momentos

Mientras que los valores de los momentos en los puntos de integración se corresponden con los valores teóricos, la extrapolación conlleva a una pérdida de precisión cuando se utiliza la función paraboloide hiperbólica para el suavizado. Un paraboloide hiperbólico sólo aproxima la distribución de momentos en el elemento. Por esta razón, se utiliza un algoritmo mejorado que sustituye la extrapolación con un método integral avanzado de esfuerzos cortantes. Dado que se supone que la distribución de esfuerzos internos en el elemento de superficie tiene la forma de un paraboloide hiperbólico (superficie de segundo orden), la integral de esta superficie representa una superficie del tercer orden, que muestra la distribución de momentos con mayor precisión.

Esta aproximación se corresponde con las ecuaciones descritas anteriormente para determinar los esfuerzos cortantes mediante las derivadas de los momentos. Los momentos se determinan así con las siguientes ecuaciones.

Las constantes de integración m_x,0 y m_y,0 se calculan desde las condiciones del valor en el centro de la barra.

Los valores de los esfuerzos internos se determinan mediante los métodos descritos anteriormente, las integrales por integración numérica. De esta forma, los momentos determinados de este modos se suavizan de esta forma según el método de los nudos internos o exteriores.

En la teoría de flexión según Mindlin, los siguientes valores resultan en este ejemplo.

• Primer elemento:
  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 10 - 5 \mathrm{x}$$

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} {\mathrm{m}}_{\mathrm{x},0} = 10 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(0.5) = 4.167 = {\int }_0^{1}10 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{10}{2} {\mathrm{x}}^2 - \frac{5}{6} {\mathrm{x}}^3 {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 4.167 {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 - 4.167 = 0$$

  $${\mathrm{m}}_{\mathrm{x}}{\left(1\right)}^{-} = 10 - \frac{5}{2} 0 = 7.5$$

  Primer elemento utilizando la teoría de flexión según Mindlin
• Segundo elemento:
  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 5 - 5 \mathrm{x}$$

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} + {\mathrm{m}}_{\mathrm{x},0} = 5 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(1.5) = 9.167 = {\int }_0^{1}5 \mathrm{x} - \frac{5}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{5}{2} {\mathrm{x}}^2 - \frac{5}{6} {\mathrm{x}}^3 + {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 1.667 + {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}(1{)}^{+} = 9.167 - 1.667 = 7.5$$

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{{\mathrm{m}}_{\mathrm{x}}(1{)}^{-} + {\mathrm{m}}_{\mathrm{x}}(1{)}^{+}}{2} = \frac{7.5 + 7.5}{2} = 7.5$$

  Segundo elemento utilizando la teoría de flexión según Mindlin

3 Distribución suavizada de los esfuerzos internos v-z y m-x según Mindlin

Se obtienen los siguientes valores para la teoría de flexión según Kirchhoff.

• Primer elemento:
  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 10.038 - 5.025 \mathrm{x}$$

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} {\mathrm{m}}_{\mathrm{x},0} = 10.038 \mathrm{x} - \frac{5.025}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(0.5) = 4.167 = {\int }_0^{1}10.038 \mathrm{x} - \frac{5.025}{2} {\mathrm{x}}^2 {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{10.038}{2} {\mathrm{x}}^2 - \frac{5.025}{6} {\mathrm{x}}^3 {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 4.182 {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}\left(0\right) = 4.167 - 4.182 = - 0.015 {\mathrm{m}}_{\mathrm{x}}{\left(1\right)}^{-} = 10.038 - \frac{5.025}{2} - 0.015 = 7.511$$

  Primer elemento utilizando la teoría de flexión según Kirchhoff
• Segundo elemento:
  $${\mathrm{v}}_{\mathrm{x}}\left(\mathrm{x}\right) = 5.013 - 5.013 \mathrm{x}$$

  $${\mathrm{m}}_{\mathrm{x}}\left(\mathrm{x}\right) = \int {\mathrm{v}}_{\mathrm{x}}\mathrm{dx} + {\mathrm{m}}_{\mathrm{x},0} = 5.013 \mathrm{x} - \frac{5.013}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},\mathrm{C}} = {\mathrm{m}}_{\mathrm{x}}(1.5) = 9.167 = {\int }_0^{1}5.013 \mathrm{x} - \frac{5.013}{2} {\mathrm{x}}^2 + {\mathrm{m}}_{\mathrm{x},0} = {\left[\frac{5.013}{2} {\mathrm{x}}^2 - \frac{5.013}{6} {\mathrm{x}}^3 + {\mathrm{m}}_{\mathrm{x},0}\mathrm{x}\right]}_0^1 = 1.671 + {\mathrm{m}}_{\mathrm{x},0}$$

  $${\mathrm{m}}_{\mathrm{x},0} = {\mathrm{m}}_{\mathrm{x}}(1{)}^{+} = 9.167 - 1.671 = 7.496$$

  $${\mathrm{m}}_{\mathrm{x}}\left(1\right) = \frac{{\mathrm{m}}_{\mathrm{x}}(1{)}^{-} + {\mathrm{m}}_{\mathrm{x}}(1{)}^{+}}{2} = \frac{7.511 + 7.496}{2} = 7.503$$

  Segundo elemento utilizando la teoría de flexión según Kirchhoff

4 Diagramas suavizados de los esfuerzos internos v-x y m-x según Kirchhoff

En el proceso de suavizado, los pasos descritos anteriormente se realizan gradualmente en el programa. Los valores suavizados se pueden mostrar gráficamente con la opción de visualización "Continuo dentro de las superficies".

Base de datos de conocimientos

KB 001606 | Suavizado de esfuerzos internos en superficies en RFEM