Métodos de integración

El método de integración predeterminado en RFEM es una cuadratura de Gauss-Lobatto con nueve puntos de integración. Esta configuración predeterminada es suficiente para la mayoría de los casos. Para representar diagramas tensión-deformación no lineales con una precisión suficiente (por ejemplo, en el caso de hormigón armado con fibras de acero), puede ser útil aumentar el número de puntos de integración. Por lo tanto, cuando se usa material no lineal , es posible personalizar el número de puntos de integración en cada capa entre tres y 99 puntos (consulte la imagen Ajustar los puntos de integración para las capas ). Tenga en cuenta, sin embargo, que un mayor número de puntos de integración da como resultado un tiempo de cálculo más largo.

Además, hay tres métodos de integración diferentes disponibles para su selección (consulte la imagen Especificación del método de integración ):

• Cuadratura de Gauss-Lobatto
• Regla de Simpson
• Regla trapezoidal

Estos métodos de integración, también conocidos como fórmulas de cuadratura, se explican a continuación. En general, se puede suponer que la cuadratura de Gauss-Lobatto con un número fijo de puntos de integración proporciona la mayor precisión sin aumentar significativamente el esfuerzo de cálculo. La regla trapezoidal o de Simpson solo conduce a un mejor resultado en casos especiales.

Fórmulas de cuadratura

El objetivo de una fórmula de cuadratura es calcular una aproximación numérica a una integral. Para esto, se seleccionan los puntos de apoyo en el dominio de integración y se agregan las evaluaciones de función ponderadas en estos puntos. El cálculo de un momento se usa como ejemplo para las siguientes explicaciones.

d es la altura de cualquier capa en una superficie multicapa. El programa intenta entonces aproximar la integral de la siguiente manera:

Para esto, n puntos de rejilla −d/2 ≤ z₁ < z₂ < . . . < z_n ≤ d/2 sobre la altura de la capa así como los pesos escalares ω₁,. . . , ω_n que son específicos de la fórmula de cuadratura respectiva. RFEM siempre usa un número impar n = 2k + 1 de puntos de rejilla y al menos la parte superior, media e inferior de la capa se selecciona como puntos de integración. Por lo tanto, da como resultado z₁ = -d/2, z_k = 0 y z_n = d/2. Esto es necesario para que los valores de tensión en estas posiciones se puedan calcular exactamente y luego mostrar en el programa. Si establece un número par de puntos de rejilla, se utiliza el siguiente número impar más alto para el cálculo. El grado de precisión de una fórmula de cuadratura se da en su orden ''p''. Este es el grado más alto de un polinomio, que está exactamente integrado. === Regla trapezoidal === Un ejemplo simple de una fórmula de cuadratura es la regla trapezoidal. En este caso, los dos bordes exteriores del dominio de integración se seleccionan como posiciones de apoyo; es decir, z1 = -d/₂ y z2 = d/₂. La regla trapezoidal es de primer orden, ya que las funciones lineales se integran exactamente, pero las funciones cuadráticas no.

En el caso de la variante de la regla trapezoidal acumulativa utilizada en RFEM, si selecciona ''n'' puntos de integración, el intervalo de integración se divide equidistantemente en n - 1 intervalos parciales y se aplica la regla trapezoidal a cada uno de ellos. La siguiente figura muestra un ejemplo de la regla trapezoidal resumida para cinco puntos de integración. === Regla de Simpson === La regla de Simpson' es una fórmula de cuadratura de tercer orden que utiliza tres puntos de integración, dos en los bordes y uno en el centro de la región de integración. Los pesos se seleccionan de tal manera que la solución aproximada corresponde a la integral sobre una parábola que pasa por estos tres puntos. Con z₁ = −d/2, z₂ = 0 y z₃ = d/2, obtenemos la forma: RFEM también usa aquí una cuadratura armada. Si hay un número impar ''n'' de puntos de apoyo, se obtienen (n − 1)/2 intervalos parciales igualmente largos a los que se puede aplicar la regla de Simpson'. La figura muestra esto de nuevo como un ejemplo para cinco puntos de integración. === Cuadratura de Gauss-Lobatto === En la cuadratura de Gauss-Lobatto, los dos puntos de contorno z₁ = −d/2 y z_n = d/2 siempre se especifican como posiciones de apoyo y todos los demás puntos de integración se seleccionan de tal manera que se logre el orden más alto posible. Esto es p = 2n − 3. Para n = 3, la cuadratura de Gauss-Lobatto es equivalente a la regla de Simpson'; para n = 5, por ejemplo: