Serie de Taylor: la aproximación de ángulos pequeños escondida dentro del análisis de segundo orden

La serie de Taylor que conociste en Cálculo II suele llegar como una manera de aproximar una función difícil con un polinomio: escribir sin, o cos, o 1/(1 menos x) como una suma infinita de derivadas, y detenerse cuando estás lo bastante cerca. Parece una herramienta para gráficas y cotas de error, lejos de una columna o de una placa base. No lo es. Cada vez que un ingeniero escribe el equilibrio sobre la forma no deformada de una estructura, está truncando en silencio una serie de Taylor después del primer término, y cada vez que corre un análisis de segundo orden, está conservando el siguiente.

Esa es toda la promesa de este artículo. Tomamos una columna de acero real, la inclinamos bajo carga, y mostramos que toda la escalera desde el primer orden hasta el segundo orden y hasta el análisis completo de grandes desplazamientos es un único desarrollo de Taylor, leído a tres profundidades de truncamiento distintas.

• El análisis de primer orden conserva solo el primer término: sin theta es theta, cos theta es 1, y la geometría nunca se mueve.
• El análisis de segundo orden (P-Delta) conserva el siguiente término: la carga axial actuando a través de la deflexión lateral, el momento P por Delta.
• El famoso amplificador 1/(1 menos P/Pcr) resulta ser una serie geométrica, un término de Taylor por cada pasada de la carga a través de la deflexión creciente.

Cada número aquí está ajustado contra el motor FEM de CalcSteel hasta el decimal. Ese es el sentido de tener un motor real detrás de las palabras: el cálculo a mano y la solución de la máquina tienen que encontrarse, y lo hacen. Todo corre gratis en el navegador, así que puedes reconstruir la columna tú mismo mientras lees.

Tres lectores sacarán algo distinto de esto. El estudiante de ingeniería obtiene una razón concreta por la que la serie de Taylor se gana su lugar en el temario de Cálculo I y II. El calculista en ejercicio obtiene un modelo mental limpio de dónde viene realmente el efecto P-Delta. El calculista independiente obtiene una forma gratuita y nativa del navegador de leer una carga crítica Pcr y formar el amplificador sin comprar una licencia de nada.

Una serie de Taylor reconstruye una función suave cerca de un punto a a partir de sus derivadas en ese único punto. Escríbela como f(x) = suma sobre n de f⁽ⁿ⁾(a)/n! por (x menos a)ⁿ: el valor, más la primera derivada por (x menos a), más la segunda derivada sobre 2! por (x menos a)², y así sucesivamente. Cuando el punto base es a = 0 la serie recibe su propio nombre, la serie de Maclaurin, y esa es la versión que usaremos para sin, cos y 1/(1 menos x).

La imagen que hay que conservar es esta: cada término añadido abraza la curva verdadera sobre un rango más amplio. Conserva un término y tienes una recta horizontal que toca la curva en a; añade el siguiente y tienes la tangente; añade el siguiente y el polinomio empieza a curvarse con la función. Detenerse pronto te da una aproximación deliberada, no un error.

Lo que te dice qué tan lejos estás es el resto, R_n. La forma de Lagrange acota el error de detenerse después de n términos por el tamaño del primer término que descartaste, evaluado en algún punto del intervalo. Ese resto es la contabilidad honesta de un truncamiento: dice exactamente cuándo conservar un término es seguro y cuándo no. Guárdalo, porque el mismo resto luego le dirá a un ingeniero estructural cuándo una aproximación de ángulos pequeños deja de ser inofensiva.

Esta es la secuela natural de Cálculo I. Allí conociste la derivada ordinaria, el único número f'(a). La serie de Taylor simplemente apila todas las derivadas en a en una sola expresión que reconstruye la función alrededor de ese punto. No se te pide nada nuevo, solo que conserves más de las derivadas que ya sabes calcular.

Apunta la máquina de Maclaurin a las funciones trigonométricas y las famosas reglas de ángulos pequeños salen de inmediato. Las series son:

• sin θ = θ − θ³/6 + θ&sup5;/120 − ..., así que para θ pequeño el primer término domina y sin θ ≈ θ.
• cos θ = 1 − θ²/2 + θ⁴/24 − ..., así que cos θ ≈ 1 (o, un término mejor, 1 − θ²/2).
• tan θ = θ + θ³/3 + 2θ⁵/15 + ..., así que tan θ ≈ θ.

Estas no son reglas empíricas sacadas de una tabla. Son series de Taylor con todos los términos después del primero descartados, y el ángulo θ está en radianes, que es la única razón por la que el coeficiente principal es un limpio 1.

¿Qué tan bueno es el truncamiento? Aquí está el error de descartar todo lo que pasa del primer término, como porcentaje, en las rotaciones que un pórtico real realmente ve:

Lee la primera columna. En 5° la aproximación sin θ ≈ θ se desvía en cerca de una parte en mil, e incluso en 10° es medio por ciento. El pórtico en servicio rota unos pocos grados a lo sumo bajo carga de servicio, así que θ genuinamente es sin θ hasta una fracción de por ciento, y el coseno, conservado con un término más, es aún mejor. Ese error diminuto es exactamente el margen en el que vive el análisis estructural de primer orden.

Ahora el eje sobre el que gira todo el artículo. Cuando escribes las ecuaciones de equilibrio de una estructura, suma de momentos y suma de fuerzas, sobre su geometría original, no deformada, estás haciendo dos suposiciones silenciosas sin nombrarlas. Ambas son truncamientos de ángulos pequeños de una serie de Taylor.

Primero, tomas cos θ ≈ 1. Físicamente esto dice que la carga axial permanece vertical y que su brazo de palanca hasta el nudo no cambia mientras el miembro se inclina, porque la geometría sobre la que escribes el equilibrio no se ha movido. Segundo, tomas sin θ ≈ θ dondequiera que una rotación entre en un balance de fuerzas. Juntas conservan solo los términos de Taylor cero y primero de la trigonometría y descartan todo lo demás. La estructura en las ecuaciones está congelada en su pose descargada.

Ese es el significado preciso del análisis de primer orden: equilibrio sobre la forma no deformada, que es equilibrio con la serie de ángulos pequeños truncada después del primer término. La deflexión que calculas después es una consecuencia de las cargas, pero nunca se realimenta de vuelta a la geometría sobre la que se escribió el equilibrio. La carga y la forma quedan desacopladas por el truncamiento.

Lo importante es que el análisis de primer orden no está equivocado. Es un truncamiento deliberado cuyo error puedes acotar con el mismo término de resto que acota cualquier serie de Taylor, y para las pequeñas rotaciones de un pórtico rígido y bien arriostrado ese error es despreciable. El problema empieza solo cuando los términos descartados crecen lo suficiente como para importar, que es exactamente el caso para el que está construido el análisis de segundo orden. Así que la manera honesta de describir el análisis de primer orden no es "el aproximado" sino "el que conserva el primer término de Taylor".

El mismo truncamiento se esconde un nivel más profundo, dentro de la propia ecuación de la viga. La curvatura exacta de un miembro flexionado no es la segunda derivada y″. Es

κ = y″ / (1 + y′²)^3/2

donde y′ es la pendiente de la forma deflectada. Ese denominador es lo que una serie de Taylor puede aplanar. Desarrolla (1 + y′²)^−3/2 como una serie binomial, que es simplemente la serie de Taylor de (1 + u)^−3/2 en u = y′²:

(1 + y′²)^−3/2 = 1 − (3/2) y′² + (15/8) y′⁴ − ...

Para un miembro cuya pendiente y′ es pequeña, es decir, cuya rotación es un ángulo pequeño, cada término después del primero es diminuto, y el factor entero es ≈ 1. La curvatura exacta colapsa a

κ ≈ y″

que es la relación lineal de Euler-Bernoulli M = EI y″ que toda fórmula de deflexión, toda matriz de rigidez y toda verificación de norma asume en silencio. Es el mismo truncamiento de ángulos pequeños que el del equilibrio en la sección anterior, de la misma familia de series de Taylor, aplicado a una cantidad distinta. Una vez que ves que M = EI y″ es una serie truncada, el resto de la escalera es fácil de subir: conserva más términos y obtienes una estructura más honesta.

Hasta ahora cada truncamiento ha descartado los términos superiores de Taylor. El análisis de segundo orden hace lo contrario: conserva el siguiente. El truco es simplemente dejar que la geometría se mueva antes de escribir el equilibrio.

Imagina nuestra columna en ménsula. Una carga lateral empuja la punta libre hacia el lado una cantidad Δ. En el análisis de primer orden la carga axial P se desliza por un miembro que, hasta donde saben las ecuaciones, sigue perfectamente vertical, así que P no aporta ningún momento en la base. Pero la columna real se ha inclinado. La carga axial P ahora cuelga por fuera de la base una distancia igual a la deflexión Δ, y añade un momento

P × Δ

en la base empotrada, encima de lo que la carga lateral ya produjo. Ese momento extra es precisamente el término que el análisis de primer orden descartó cuando congeló la geometría. Es el siguiente término en el desarrollo del equilibrio alrededor de la forma deformada: en vez de detenerte en la pose no deformada, llevas el balance un paso de Taylor más allá y dejas que la carga actúe a través del desplazamiento que causó.

El nombre P-Delta viene directamente de este producto. La carga P por el desplazamiento lateral Δ es el momento de segundo orden, y cualquier análisis que lo incluya es un análisis de segundo orden. Intuitivamente, la ménsula inclinada es la imagen más limpia posible: cuanto más se inclina, más sobresale la carga axial de la base, más momento genera, lo que la inclina más. Ese lazo de realimentación es toda la historia, y la siguiente sección lo convierte en una serie que ya conoces.

Sigue el lazo de realimentación una pasada a la vez y el amplificador se escribe solo. Aplica la carga lateral H y la columna se deflecta a cierto valor de primer orden Δ₀. Ahora la carga axial P, colgando por fuera de la base una distancia Δ₀, añade un momento P por Δ₀ que flexiona la columna aún más, en un extra Δ₁ = Δ₀ × ρ, donde ρ = P/Pcr es la fracción del camino hasta la carga crítica. Pero esa nueva deflexión deja que P sobresalga aún más, añadiendo Δ₂ = Δ₀ × ρ², luego Δ₃ = Δ₀ × ρ³, y así sin fin.

La deflexión total es la suma de cada pasada:

Δ = Δ₀ (1 + ρ + ρ² + ρ³ + ...)

Ese paréntesis es una serie geométrica, la serie de Maclaurin de 1/(1 menos x) que conociste en Cálculo II, evaluada en x = ρ. Para ρ menor que 1 suma a una forma cerrada limpia, 1/(1 menos ρ), así que

Δ = Δ₀ / (1 − ρ)

y el momento se amplifica por el mismísimo factor. El amplificador es

B = 1 / (1 − P/Pcr)

y ahora puedes leerlo por lo que es: cada pasada de la carga axial a través de la deflexión creciente es un término más de la serie de Taylor. El análisis de primer orden conserva solo el 1 principal. El análisis de segundo orden, en su forma de primer orden amplificado, suma la serie entera. El Pcr en el denominador es la carga crítica de Euler de la columna: la carga en la que ρ alcanza 1, la serie deja de converger, y el amplificador estalla. Esa divergencia no es un error del álgebra, es el pandeo, contado en el lenguaje de una serie de Taylor que se ha quedado sin margen.

Aquí está la estructura que llevaremos por el resto del artículo, con números reales y un motor real detrás de cada uno de ellos. Es una columna en ménsula de acero, perfil HEB 200, altura L = 4.0 m, empotrada en la base y libre en la punta. La flexión y el pandeo gobiernan ambos alrededor del eje débil, cuyo segundo momento de área es I_y = 2001.54 cm⁴, y el acero tiene E = 20000 kN/cm².

Paso 1, la respuesta de primer orden. Empuja la punta libre de lado con una carga lateral H = 10 kN. El momento en la base es simplemente H por la altura:

M₀ = H · L = 10 × 4.0 = 40.0 kNm.

El motor FEM de CalcSteel, el mismo que corre dentro de la app del navegador, devuelve 40.0000 kNm, una diferencia del 0.000%. La deflexión de primer orden de la punta de una ménsula es la clásica cúbica:

Δ₁ = H · L³ / (3 E I) = 53.29 mm,

y el motor devuelve 53.29 mm hasta cada dígito, porque el elemento de viga cúbico es exacto para esta carga.

Paso 2, la carga crítica de Euler. El análisis de pandeo lineal del motor da Pcr = 617.32 kN. La forma cerrada para una columna empotrada-libre es Pcr = π² E I / (2L)², con factor de longitud efectiva K = 2 porque la ménsula libre pandea en un cuarto de onda senoidal. Esa fórmula da 617.32 kN. Motor y mano concuerdan al 0.000%.

Paso 3, añade carga axial y amplifica. Ahora presiona una carga axial de servicio P = 247 kN hacia abajo por la misma columna. La razón de carga es

ρ = P / Pcr = 247 / 617.32 = 0.400,

así que el amplificador es B = 1 / (1 − 0.400) = 1.667. Aplicándolo a los resultados de primer orden:

M₂ = B · M₀ = 1.667 × 40.0 = 66.68 kNm, y Δ₂ = B · Δ₁ = 88.84 mm.

Léelo con claridad: apenas al 40% del camino hasta Pcr, los efectos de segundo orden han añadido un 67% al momento de primer orden. Los 26.68 kNm extra (es decir, 66.68 menos 40.0) no son otra cosa que la suma de cada término de Taylor después del primero. Puedes ver las sumas parciales geométricas trepar hacia el amplificador, un término por cada pasada de la carga:

El análisis de primer orden conserva solo el 1.0000 principal. Cada fila extra es un término más de la misma serie convergiendo a 1.667.

Una palabra honesta sobre lo que el motor hizo y no hizo. El motor FEM da la respuesta exacta de primer orden y la Pcr de Euler exacta, ajustadas al decimal arriba. El amplificador es la clásica estimación de segundo orden en forma cerrada, un análisis elástico de primer orden amplificado (el método B1/B2 de AISC), no una resolución aparte. Es la serie geométrica en acción, y es precisa precisamente mientras la serie converge, que es el tema de una sección posterior.

Acabas de ver una serie geométrica convertirse en un momento real sobre el papel. Ahora hazlo tú mismo, sin tomar un lápiz. La calculadora incrustada justo aquí corre un análisis de pandeo sobre una columna que tú defines, así que el abstracto Pcr se vuelve un número que puedes mover a voluntad.

Introduce tus propias altura, perfil y condiciones de apoyo, lee la carga crítica Pcr directamente, luego elige una carga axial de servicio P y forma la razón ρ = P/Pcr y el amplificador B = 1/(1 − ρ) tú mismo. Cambia la altura, mira caer Pcr y trepar el amplificador. Es el mismo lazo que la columna resuelta: la serie geométrica 1 + ρ + ρ² + ... hecha visible.

Es genuinamente gratis, no hay registro para las cuentas, y no hay límite de cuántas veces la corras. Explórala como explorarías una serie en Cálculo II: cambia una entrada, mantén el resto fijo, y lee qué tan rápido crece el amplificador a medida que te acercas a Pcr.

Todo truncamiento tiene una etiqueta de advertencia, y para la serie de Taylor es el término de resto. Mira atrás a la tabla de ángulos pequeños: el error de sin θ ≈ θ era 0.127% en 5° pero había crecido a 4.720% en 30°. El resto crece con el ángulo, y también la pregunta honesta de si el truncamiento aún se sostiene. Hay dos maneras distintas en que puede fallar.

Falla uno, las rotaciones mismas se vuelven grandes. Cables esbeltos, miembros muy flexibles y problemas de búsqueda de forma pueden rotar lo suficiente como para que sin θ ya no sea θ: el primer término de Taylor descartado es ahora una fracción real de la respuesta. Aquí ningún amplificador te salva, porque la geometría linealizada sobre la que el amplificador está construido ha caducado ella misma. Necesitas un verdadero análisis de grandes desplazamientos, geométricamente no lineal que conserve sin θ completo en vez de truncado, y actualice la geometría a medida que la estructura se deforma. Este es el régimen donde la imagen de Taylor deja de ser un atajo y el motor tiene que cargar la trigonometría completa.

Falla dos, la propia serie del amplificador deja de converger. A medida que la carga axial P trepa hacia Pcr, la razón ρ = P/Pcr se acerca a 1, y la serie geométrica 1 + ρ + ρ² + ... diverge: B = 1/(1 − ρ) se dispara al infinito. Mucho antes de que literalmente diverja, el atajo de primer orden amplificado pierde precisión, y por eso las normas lo acotan. AISC limita el uso del amplificador B1/B2, y el Eurocode EN 1993-1-1 pide un análisis completo de segundo orden siempre que la razón crítica elástica α_cr = Pcr/P caiga por debajo de 10. Esa α_cr es simplemente 1/ρ: α_cr ≥ 10 significa ρ ≤ 0.1, la carga está a lo sumo a una décima parte del camino hasta el pandeo, la serie converge rápido, y la imagen de ángulos pequeños está cómoda. Cuando α_cr < 10, la carga está lo bastante cerca de Pcr como para que los términos superiores importen y se requiera una resolución rigurosa de segundo orden. Nuestra columna resuelta se situaba en ρ = 0.400, así que α_cr = 2.5, bien dentro de la región donde un análisis riguroso de segundo orden, no un simple amplificador, es la elección honesta.

El término de resto es lo que hace ambos límites precisos en vez de una cuestión de gusto. Es la etiqueta de advertencia honesta impresa sobre la aproximación de ángulos pequeños: segura mientras los términos descartados sean diminutos, peligrosa en el momento en que no lo sean.

Da un paso atrás y todo el arco desde Cálculo I y II aterriza sobre acero real, una herramienta a la vez. La derivada ordinaria da la tasa de cambio a lo largo de una viga, dM/dx es igual a la fuerza cortante V. La segunda derivada da la curvatura de la línea elástica, M = EI y″. La integral definida recupera áreas y propiedades de sección. La derivada parcial da una deflexión directamente a partir de la energía almacenada, como muestra el teorema de Castigliano. Y ahora la serie de Taylor da el amplificador de segundo orden, y revela que el análisis de primer orden que subyace a todos los demás era un truncamiento de esa misma serie.

Es un único cálculo coherente, no un estante de trucos inconexos. La serie de Taylor es la herramienta que finalmente nombra la aproximación que el resto de ellas hacían en silencio, y te da el siguiente término cuando lo necesitas. Si quieres ver de dónde viene la carga crítica Pcr del denominador, nuestro artículo de pandeo de Euler la deriva de principio a fin. El mismo acero, un miembro más de la familia del cálculo.

El amplificador no es una curiosidad de aula. Está escrito, casi letra por letra, en las normas de diseño que un ingeniero firma. Cada una de ellas toma la misma serie geométrica y la convierte en una verificación obligatoria.

• AISC 360 permite un análisis elástico de primer orden amplificado con dos factores: B1 para el efecto P-δ del miembro (una carga axial actuando a través de la deflexión dentro de un miembro) y B2 para el efecto P-Δ del piso (la carga actuando a través del desplazamiento lateral del piso entero). Ambos son de la forma 1/(1 − α P_r/P_e), la misma serie 1/(1 − ρ) con P_e la carga de pandeo elástico. El Capítulo C establece el método de análisis directo y el Apéndice 8 la ruta B1/B2.
• Eurocode EN 1993-1-1 ancla toda la decisión a α_cr: si α_cr ≥ 10 (elástico) los efectos de segundo orden son lo bastante pequeños para ignorarse o manejarse por amplificación, y por debajo de eso se requiere un análisis de segundo orden más completo. Ese es el test de convergencia de la serie, ascendido a cláusula.
• ABNT NBR 8800 lleva sus propios coeficientes B1 y B2 con el mismo espíritu de primer orden amplificado, así que un calculista brasileño está sumando la idéntica serie de Taylor bajo notación brasileña.

Dos vecinos de la misma familia de normas valen un enlace. El desplazamiento lateral que B2 amplifica tiene que ser disparado en primer lugar, que es el trabajo de las cargas nocionales, pequeñas fuerzas laterales que exponen la inestabilidad de un pórtico imperfecto sin una corrida completa de P-Delta. Y el momento amplificado no actúa solo: se combina con la carga axial en la verificación de interacción de la flexión y axial combinadas. En cada caso el punto es el mismo: cuántos términos de Taylor conservas es una decisión de norma, no una corazonada. La norma te dice cuándo el primer término basta y cuándo debes sumar el resto.

La idea es corta, que es exactamente por lo que los errores son fáciles de cometer. Aquí está la lista de verificación que mantiene honesta una estimación de segundo orden.

• "Ángulo pequeño" no es "deflexión pequeña". El truncamiento que hace el análisis de primer orden es sobre la rotación, no sobre el desplazamiento. Una ménsula esbelta puede mover su punta un buen trecho (nuestra columna resuelta alcanzó 88.84 mm) mientras sus miembros aún rotan solo unos pocos grados, así que sin θ ≈ θ sigue siendo válido aunque la deflexión sea visiblemente grande. Confundir las dos cosas lleva a la gente a desconfiar de un análisis perfectamente bueno, o a confiar en uno malo.
• Usar ρ = P/Pcr con la Pcr equivocada. Este es el error real habitual. La carga crítica depende del factor de longitud efectiva K y de cuál eje pandea; mete la Pcr del eje fuerte, o K = 1 para una columna que es en realidad empotrada-libre (K = 2), y tu amplificador no significa nada. Nota que nuestra columna resuelta pandeó alrededor del eje débil, con I_y = 2001.54 cm⁴, que es el eje que gobierna. Acierta con Pcr antes de formar ρ.
• El segundo orden no es lo mismo que la no linealidad del material. Todo el análisis en este artículo es elástico lineal: el acero nunca fluye, la ley de Hooke se cumple en todo momento. La no linealidad de segundo orden (geométrica) trata sobre el equilibrio en la forma deformada, algo distinto de una curva tensión-deformación no lineal. Un problema puede ser una cosa, la otra, ambas, o ninguna.
• B estalla cerca de Pcr por una razón real. Cuando alguien ve 1/(1 − ρ) dispararse al infinito a veces supone un fallo numérico. Es física: en ρ = 1 la carga ha alcanzado el pandeo de Euler, y la amplificación infinita es la cara matemática honesta de una columna inestable. Respeta la divergencia, no la recortes.
• El primer orden es un truncamiento, no un error. Mientras α_cr = Pcr/P sea cómoda (10 o más en términos del Eurocode), conservar solo el primer término de Taylor es una elección legítima y avalada por la norma. El análisis de primer orden se gana su lugar precisamente porque los términos descartados son demostrablemente diminutos allí.

Una columna inclinada, una serie de Taylor, leída a tres profundidades de truncamiento, y cada número de primer orden y de Euler ajustado a un motor FEM real hasta el decimal.

• El análisis de primer orden es un truncamiento de ángulos pequeños. El equilibrio sobre la forma no deformada mantiene sin θ como θ y cos θ como 1, solo el primer término de Taylor, y la geometría nunca se mueve.
• El amplificador es una serie geométrica. B = 1/(1 − P/Pcr) es exactamente 1 + ρ + ρ² + ρ³ + ..., un término por cada pasada de la carga axial a través de la deflexión creciente.
• La columna HEB 200 resuelta. El primer orden dio M₀ = 40.0 kNm y Δ₁ = 53.29 mm; en ρ = 0.400 el amplificador B = 1.667 los elevó a M₂ = 66.68 kNm y Δ₂ = 88.84 mm, 67% más momento apenas al 40% del camino hasta Pcr.
• El motor coincidió con el cálculo a mano hasta el decimal. Su análisis de pandeo devolvió Pcr = 617.32 kN, igual a la forma cerrada de Euler π²EI/(2L)² con K = 2 al 0.000%.
• El resto es la etiqueta de advertencia. Cerca de Pcr la serie diverge, y en rotaciones grandes sin θ deja de ser igual a θ. Allí el truncamiento falla y necesitas un análisis completo de segundo orden o de grandes desplazamientos.

Ahora ponle las manos encima. Construye la columna, lee Pcr, y forma el amplificador tú mismo en CalcSteel, una herramienta de diseño y análisis de acero estructural nativa del navegador con un motor FEM genuino. Las cuentas y las calculadoras son gratis, sin registro, ilimitadas, y si eres estudiante la oferta /education es gratis para ti. La prueba no es un ensayo con cuenta regresiva; es una herramienta que puedes abrir ahora mismo y verificar contra cada número de esta página.