Flecha máxima viga carga

Created on: January 23, 2025

Answered using O1 by Chat01

Question

Calcula la flecha máxima de una viga biapoyada con forma de pentágono, de acero, con longitud L y bajo aplicacion de una carga que corresponda a la función F(x)=log(x^3) donde x es la posición en metros desde uno de los extremos de la viga

Thought

Calculating maximum deflection

Investigating the maximum deflection of a pentagon-shaped steel beam under a load function, leaning on the connections between cross-sectional geometry and deflection.

Assessing beam dynamics

OK, let me see. Integrating a beam's deflection equation four times under a load function. The undefined at x=0 complicates this.

Considering load implications

Reflecting on the non-symmetrical load and singularity at x=0, which likely complicates determining the maximum deflection. This makes the problem quite intricate.

Revisiting the coordinate system

I’m reconciling the problem, possibly adjusting coordinates to avoid singularities and shift the focus to a simpler approximation or integral approach for deflection analysis.

Integrating step-by-step

First, antiderivate 3lnx3 \ln x. Then, by parts, tackle xlnxx \ln x for the second derivative. This method carefully maps out the process.

Tackling integration

OK, let me see. I’m pulling together the second derivative d2ydx2\frac{d^2 y}{dx^2} by combining terms and integrating each term to simplify the function for further examination.

Breaking down the integral

OK, let me see. Substituting lnx\ln x and integrating x2x^2 gives uvvu2\frac{u v - v u'}{2}. This approach streamlines the integration process.

Piecing together derivatives

Blending terms and integrating constant multiples, I combine like terms of x3x^3 to simplify the first derivative dydx\frac{dy}{dx} step-by-step.

Pulling together the elements

OK, let me see. Integrating the remaining terms to find y(x)y(x), considering integer exponents and logarithms. This leads to a comprehensive expression for yy.

Integrating each term

I'm gathering all the ingredients for y(x)y(x), meticulously piecing together terms like x48lnx\frac{x^4}{8} \ln x and C1x36\frac{C_1 x^3}{6}, ensuring no detail is overlooked.

Layers of boundary conditions

I'm curious about setting y(0)=0y(0) = 0 and y(L)=0y(L) = 0 for deflection and y(0)=0y''(0) = 0 and y(L)=0y''(L) = 0 for moments, simplifying the analysis with these constraints.

Taking a closer look

Analyzing as xx nears 0, I notice x48lnx\frac{x^4}{8} \ln x trends to 0, leading to C4=0C_4 = 0. This simplifies the system significantly.

Breaking down the boundary

OK, let’s see. Evaluating y(x)y''(x) at x=0x=0 simplifies to C2C_2. So, C2C_2 must be zero to satisfy the boundary condition y(0)=0y''(0) = 0.

Pinpointing coefficients

I’m working through boundary conditions to pinpoint coefficients C1C_1 and C3C_3. This involves setting up equations for y(L)y(L) and y(L)y''(L) and solving for these unknowns.

Balancing terms

I'm aligning terms to find C1C_1 and C3C_3. By comparing expressions on both sides, I ensure the equation is balanced and consistent with boundary conditions.

Rearranging terms

I’m organizing terms to simplify the left side and adjust for coefficients. This creates a more balanced equation, hinting at solving for C1C_1 and C3C_3 by combining and isolating terms.

Rearranging the equation

I’m working through isolating C3C_3 by tweaking the equation to group terms with L3lnLL^3 \ln L. This approach clarifies the solution and simplifies the process.

Finalizing coefficients

I've pinpointed C1C_1, C2C_2, and C3C_3, concluding with the coefficients necessary for a precise deflection function for the simply supported beam.

Calculating the peak

Now I evaluate the deflection function y(x)y(x) by deriving it with respect to xx and identifying the xx value where y(x)y(x) likely reaches its maximum by setting the derivative equal to zero.

Mapping out the peak

I’m setting y(x)y'(x) to zero and factoring out x2x^2. This helps identify the critical points where the deflection function reaches its maximum, using the derived coefficients.

Taking a closer look

I’m considering a dimensionless approach by introducing α=x/L\alpha = x/L, rewriting the equation in terms of α\alpha. This might help in identifying the dimensionless quantity for maximum deflection.

Transforming the equation

I’m translating y(x)y'(x) into dimensionless terms using α=x/L\alpha = x/L. This approach simplifies the coefficients C1C_1 and C3C_3 for a clearer depiction of maximum deflection.

Adjusting my approach

I’m honing in on the terms with lnL\ln L, factoring them out methodically. This involves carefully balancing each term to ensure a precise dimensionless form for the deflection function.

Breaking down the equation

I’m fine-tuning the deflection function by isolating terms with lnL\ln L and those without. This approach links dimensionless α\alpha to the specified LL, aiming for a comprehensive solution.

Weighing deflection dynamics

I'm mapping out how load placement impacts deflection. Given the indefinite load near x=0x = 0, peak deflection may be near x=Lx = L.

Wrestling with boundaries

Figuring out peak deflection involves solving a complex boundary value problem, often leading to intricate solutions and typically lacking a straightforward close-form expression.

Thought about flecha máxima en viga pentagonal durante 2m 58s

Answer

Análisis General

Se desea calcular la deflexión (o “flecha”) de una viga simplemente apoyada en sus dos extremos x=0x=0 y x=Lx=L, de sección pentagonal de acero (lo cual afecta al momento de inercia II, pero no altera la forma básica de las ecuaciones de la flexión), sometida a una carga distribuida que depende de la posición como

w(x)  =  F(x)  =  log(x3)  =  3ln(x),w(x) \;=\; F(x)\;=\;\log\bigl(x^{3}\bigr)\;=\;3\,\ln(x),

donde (suponiendo) ln()\ln(\cdot) denota logaritmo natural y xx está en metros.

Nota sobre la singularidad en x=0x=0
Físicamente, ln(x)\ln(x) tiende a -\infty a medida que x0+x \to 0^{+}. Esto implica que, estrictamente, la función de carga w(x)=3ln(x)w(x) = 3\,\ln(x) no está bien definida justo en el apoyo izquierdo (x=0x=0). En la práctica, podría interpretarse que la carga real se aplica solo a partir de un cierto tramo (x>0x>0) o que se trata de un ejercicio meramente académico para ilustrar la integración con un tipo de carga “no usual”.
En lo que sigue, se procede a la integración formal y a la imposición de condiciones de apoyo en x=0x=0 y x=Lx=L. Sin embargo, téngase presente que, físicamente, habría que revisar cómo manejar la vecindad de x=0x=0.


1. Ecuación diferencial de la viga

Para una viga de longitud LL, apoyada en x=0x=0 y x=Lx=L, la ecuación de flexión (en hipótesis de Bernoulli–Euler) es:

EId4ydx4  =  w(x),E\,I\,\frac{d^{4}y}{dx^{4}} \;=\; w(x),

donde:

  • y(x)y(x) es la flecha (deflexión) vertical de la viga en la posición xx.
  • EE es el módulo de elasticidad (para acero).
  • II es el momento de inercia de la sección transversal respecto al eje neutro.
  • w(x)=3ln(x)w(x)=3\,\ln(x) en este problema (unidades de fuerza por unidad de longitud, p.e. N/m).

Por comodidad de notación, en muchos desarrollos uno escribe directamente

d4ydx4  =  w(x)EI  =  3ln(x)EI.\frac{d^{4}y}{dx^{4}} \;=\; \frac{w(x)}{E\,I} \;=\; \frac{3\,\ln(x)}{E\,I}.

Aquí mostraremos primero la forma general (omitiendo, de momento, la constante 1/(EI)1/(E\,I)), y al final se puede factorizar si se desea.


2. Integración paso a paso

Sea

y(x)  =  3ln(x).y''''(x) \;=\; 3\,\ln(x).

(Ignoramos la constante 1EI\frac{1}{E\,I} para mantener las integraciones más limpias; luego se reintroduce multiplicando donde corresponda, o bien se sobreentiende que cualquier “3” es en realidad 3/(EI)3/(E\,I).)

2.1. Primera integración

y(x)  =  3ln(x)  dx  =  3ln(x)  dx.y'''(x) \;=\; \int 3\,\ln(x)\;dx \;=\; 3\int \ln(x)\;dx.

Sabemos que

ln(x)dx  =  xln(x)    x  +  C,\int \ln(x)\,dx \;=\; x\,\ln(x) \;-\; x \;+\; C,

así que

y(x)  =  3[xln(x)    x]  +  C1  =  3xln(x)    3x  +  C1.y'''(x) \;=\; 3 \bigl[x\,\ln(x) \;-\; x\bigr] \;+\; C_{1} \;=\; 3\,x\,\ln(x) \;-\; 3\,x \;+\; C_{1}.

2.2. Segunda integración

y(x)  =  (3xln(x)3x+C1)dx.y''(x) \;=\; \int \bigl(3\,x\,\ln(x) - 3\,x + C_{1}\bigr)\,dx.

Separamos términos:

  1. 3xln(x)dx=3xln(x)dx.\displaystyle \int 3\,x\,\ln(x)\,dx = 3\int x\,\ln(x)\,dx.

    • Mediante partes,
      (\int x,\ln(x),dx = \tfrac{x^2}{2},\ln(x) - \tfrac{x^2}{4}.)

    • Por tanto,
      (\displaystyle 3 \int x,\ln(x),dx = 3\Bigl(\tfrac{x^2}{2},\ln(x) - \tfrac{x^2}{4}\Bigr) = \tfrac{3,x^2}{2},\ln(x) - \tfrac{3,x^2}{4}.)

  2. (\displaystyle \int \bigl(-3,x\bigr),dx = -3;\tfrac{x^2}{2} = -\tfrac{3,x^2}{2}.)

  3. C1  dx=C1x.\displaystyle \int C_{1}\;dx = C_{1}\,x.

Agrupando,

y(x)=3x22ln(x)    3x24    3x22  +  C1x  +  C2.y''(x) = \tfrac{3\,x^2}{2}\,\ln(x) \;-\; \tfrac{3\,x^2}{4} \;-\; \tfrac{3\,x^2}{2} \;+\; C_{1}\,x \;+\; C_{2}.

La parte puramente polinómica:

  3x24    3x22  =    3x24    6x24  =    9x24.-\;\tfrac{3\,x^2}{4} \;-\; \tfrac{3\,x^2}{2} \;=\; -\;\tfrac{3\,x^2}{4} \;-\; \tfrac{6\,x^2}{4} \;=\; -\;\tfrac{9\,x^2}{4}.

Así,

y(x)  =  3x22ln(x)    9x24  +  C1x  +  C2.y''(x) \;=\; \tfrac{3\,x^2}{2}\,\ln(x) \;-\; \tfrac{9\,x^2}{4} \;+\; C_{1}\,x \;+\; C_{2}.

2.3. Tercera integración

y(x)  =  y(x)dx  =  [3x22ln(x)9x24+C1x+C2]dx.y'(x) \;=\; \int y''(x)\,dx \;=\; \int \Bigl[\tfrac{3\,x^2}{2}\,\ln(x) - \tfrac{9\,x^2}{4} + C_{1}\,x + C_{2}\Bigr]\,dx.

Por partes:

  • x2ln(x)dx:\displaystyle \int x^2 \ln(x)\,dx:
    x2ln(x)dx=x33ln(x)x39.\int x^2 \ln(x)\,dx = \frac{x^3}{3}\,\ln(x) - \frac{x^3}{9}.

    Por tanto,
    (\tfrac{3}{2}\int x^2 \ln(x),dx = \tfrac{3}{2}\Bigl[\frac{x^3}{3},\ln(x) - \frac{x^3}{9}\Bigr] = \frac{x^3}{2},\ln(x) - \frac{x^3}{6}.)

  • (\displaystyle \int \bigl(- \tfrac{9,x^2}{4}\bigr),dx = -;\tfrac{9}{4};\frac{x^3}{3} = -;\tfrac{3,x^3}{4}.)

  • (\displaystyle \int C_{1},x,dx = C_{1},\tfrac{x^2}{2}.)

  • (\displaystyle \int C_{2},dx = C_{2},x.)

Sumando todo:

y(x)=(x32ln(x)x36)    3x34  +  C1x22  +  C2x  +  C3.y'(x) = \Bigl(\tfrac{x^3}{2}\,\ln(x) - \tfrac{x^3}{6}\Bigr) \;-\; \tfrac{3\,x^3}{4} \;+\; C_{1}\,\tfrac{x^2}{2} \;+\; C_{2}\,x \;+\; C_{3}.

Se combinan los términos en x3x^3:

  x36    3x34  =  x36    9x312  =  2x312    9x312  =  11x312.-\;\tfrac{x^3}{6} \;-\; \tfrac{3\,x^3}{4} \;=\; -\,\tfrac{x^3}{6} \;-\; \tfrac{9\,x^3}{12} \;=\; -\,\tfrac{2\,x^3}{12} \;-\; \tfrac{9\,x^3}{12} \;=\; -\,\tfrac{11\,x^3}{12}.

Por tanto,

y(x)=x32ln(x)    11x312  +  C1x22  +  C2x  +  C3.y'(x) = \tfrac{x^3}{2}\,\ln(x) \;-\; \tfrac{11\,x^3}{12} \;+\; \tfrac{C_{1}\,x^2}{2} \;+\; C_{2}\,x \;+\; C_{3}.

2.4. Cuarta (y última) integración

y(x)  =  y(x)dx  =  [x32ln(x)11x312+C1x22+C2x+C3]dx.y(x) \;=\; \int y'(x)\,dx \;=\; \int \Bigl[\tfrac{x^3}{2}\,\ln(x) - \tfrac{11\,x^3}{12} + \tfrac{C_{1}\,x^2}{2} + C_{2}\,x + C_{3}\Bigr]\,dx.
  • x3ln(x)dx:\displaystyle \int x^3\ln(x)\,dx:
    Por partes,
    (\int x^3\ln(x),dx = \frac{x^4}{4},\ln(x) - \frac{x^4}{16}.)

    Por tanto,
    (\tfrac{1}{2}\int x^3\ln(x),dx = \frac{1}{2}\Bigl[\frac{x^4}{4},\ln(x) - \frac{x^4}{16}\Bigr] = \frac{x^4}{8},\ln(x) - \frac{x^4}{32}.)

  • (\displaystyle \int \bigl(-\tfrac{11,x^3}{12}\bigr),dx = -;\tfrac{11}{12};\frac{x^4}{4} = -;\tfrac{11,x^4}{48}.)

  • (\displaystyle \int \tfrac{C_{1},x^2}{2},dx = \tfrac{C_{1}}{2};\frac{x^3}{3} = \tfrac{C_{1},x^3}{6}.)

  • (\displaystyle \int C_{2},x,dx = C_{2},\tfrac{x^2}{2}.)

  • (\displaystyle \int C_{3},dx = C_{3},x.)

Así,

y(x)=(x48ln(x)x432)    11x448  +  C1x36  +  C2x22  +  C3x  +  C4.y(x) = \Bigl(\frac{x^4}{8}\,\ln(x) - \frac{x^4}{32}\Bigr) \;-\; \frac{11\,x^4}{48} \;+\; \frac{C_{1}\,x^3}{6} \;+\; C_{2}\,\frac{x^2}{2} \;+\; C_{3}\,x \;+\; C_{4}.

Agrupando los términos en x4x^4:

  x432    11x448=  3x496    22x496=25x496.-\;\frac{x^4}{32} \;-\; \frac{11\,x^4}{48} = -\;\frac{3\,x^4}{96} \;-\; \frac{22\,x^4}{96} = -\,\frac{25\,x^4}{96}.

Por tanto,

y(x)=x48ln(x)    25x496  +  C1x36  +  C2x22  +  C3x  +  C4.\boxed{ y(x) = \frac{x^4}{8}\,\ln(x) \;-\; \frac{25\,x^4}{96} \;+\; \frac{C_{1}\,x^3}{6} \;+\; \frac{C_{2}\,x^2}{2} \;+\; C_{3}\,x \;+\; C_{4}. }

3. Condiciones de contorno (viga simplemente apoyada)

Para una viga biapoyada en x=0x=0 y x=Lx=L, las condiciones usuales son:

  1. Deflexión nula en los apoyos: y(0)  =  0yy(L)  =  0.y(0) \;=\; 0 \quad\text{y}\quad y(L) \;=\; 0.
  2. Momento flector nulo en los apoyos: M(x)=EIy(x)y(0)=0,  y(L)=0.M(x) = E\,I\,y''(x)\quad \Longrightarrow \quad y''(0)=0,\; y''(L)=0.

Estas condiciones son las apropiadas para una viga articulada (o simplemente apoyada) en los extremos.
Nótese la dificultad en x=0x=0 por la presencia de ln(x)\ln(x). Aun así, tomaremos el límite conforme x0+x\to 0^{+} y veremos que, gracias al factor xnx^n, varios términos pueden tender a 0.

3.1. Aplicando y(0)=0y(0)=0

Miremos la expresión de y(x)y(x) cuando x0x\to 0:

x48ln(x)    0(puesto que x40 maˊs raˊpido que ln(x)).\frac{x^4}{8}\ln(x) \;\to\; 0 \quad (\text{puesto que }x^4\to 0\text{ más rápido que }\ln(x)\to -\infty).

Los demás términos polinómicos en xx también van a 0 (pues cada uno tiene factor xnx^n, n1n\ge1), excepto C4C_{4}. Por tanto,

0=y(0)=limx0y(x)=C4        C4=0.0 = y(0) = \lim_{x\to 0}y(x) = C_{4} \;\;\Longrightarrow\;\; C_{4} = 0.

3.2. Aplicando y(0)=0y''(0)=0

Recordemos que

y(x)=3x22ln(x)    9x24  +  C1x  +  C2.y''(x) = \frac{3\,x^2}{2}\,\ln(x) \;-\; \frac{9\,x^2}{4} \;+\; C_{1}\,x \;+\; C_{2}.

Cuando x0x\to 0, el término x2ln(x)0x^2\ln(x)\to 0, el término 9x240-\tfrac{9\,x^2}{4}\to 0, y C1x0C_{1}\,x\to 0. El único que podría quedar es C2C_{2}. Por tanto,

0=y(0)=limx0y(x)=C2        C2=0.0 = y''(0) = \lim_{x\to 0} y''(x) = C_{2} \;\;\Longrightarrow\;\; C_{2} = 0.

3.3. Aplicando y(L)=0y(L)=0

Conocemos ya C4=0C_{4}=0, C2=0C_{2}=0. Entonces:

y(L)=L48ln(L)    25L496  +  C1L36  +  C3L  =  0.y(L) = \frac{L^4}{8}\,\ln(L) \;-\; \frac{25\,L^4}{96} \;+\; \frac{C_{1}\,L^3}{6} \;+\; C_{3}\,L \;=\; 0.

Reordenamos:

C1L36  +  C3L  =  25L496    L48ln(L).\frac{C_{1}\,L^3}{6} \;+\; C_{3}\,L \;=\; \frac{25\,L^4}{96} \;-\; \frac{L^4}{8}\,\ln(L).

Podemos factorizar LL a la izquierda:

L  (C1L26+C3)  =  L4  (259618ln(L)).L\;\biggl(\,\tfrac{C_{1}\,L^2}{6} + C_{3}\biggr) \;=\; L^4\;\Bigl(\tfrac{25}{96} - \tfrac{1}{8}\ln(L)\Bigr).

Dividiendo entre L0L\neq 0,

C1L26+C3  =  L3  (259618ln(L)).(1)\tfrac{C_{1}\,L^2}{6} + C_{3} \;=\; L^3\;\Bigl(\tfrac{25}{96} - \tfrac{1}{8}\ln(L)\Bigr). \tag{1}

3.4. Aplicando y(L)=0y''(L)=0

y(L)=3L22ln(L)    9L24  +  C1L  =  0.y''(L) = \frac{3\,L^2}{2}\,\ln(L) \;-\; \frac{9\,L^2}{4} \;+\; C_{1}\,L \;=\; 0.

Despejamos C1C_{1}:

C1L=9L24    3L22ln(L)        C1=9L4    3L2ln(L).(2)C_{1}\,L = \frac{9\,L^2}{4} \;-\; \frac{3\,L^2}{2}\,\ln(L) \;\;\Longrightarrow\;\; C_{1} = \frac{9\,L}{4} \;-\; \frac{3\,L}{2}\,\ln(L). \tag{2}

3.5. Hallando C3C_{3}

Sustituyendo C1\displaystyle C_{1} de (2)(2) en la ecuación (1)(1):

C1L26+C3  =  L3  (259618ln(L)).\frac{C_{1}\,L^2}{6} + C_{3} \;=\; L^3\;\Bigl(\tfrac{25}{96} - \tfrac{1}{8}\ln(L)\Bigr).

Con

C1=L(9432ln(L)),C_{1} = L\Bigl(\tfrac{9}{4} - \tfrac{3}{2}\ln(L)\Bigr),

tenemos

C1L26=L36  (9432ln(L))=L36  94    L36  32ln(L)=3L38    3L312ln(L).\frac{C_{1}\,L^2}{6} = \frac{L^3}{6}\;\Bigl(\tfrac{9}{4} - \tfrac{3}{2}\ln(L)\Bigr) = \frac{L^3}{6}\;\frac{9}{4} \;-\; \frac{L^3}{6}\;\frac{3}{2}\ln(L) = \frac{3\,L^3}{8} \;-\; \frac{3\,L^3}{12}\,\ln(L).

(Si se hace con más cuidado usando fracciones:
94÷6=94×16=924=38,\tfrac{9}{4}\div 6 = \tfrac{9}{4}\times \tfrac{1}{6} = \tfrac{9}{24} = \tfrac{3}{8},
32÷6=32×16=312.\tfrac{3}{2}\div 6 = \tfrac{3}{2}\times \tfrac{1}{6} = \tfrac{3}{12}.)

Por tanto,

C1L26  =  3L38    3L312ln(L).\frac{C_{1}\,L^2}{6} \;=\; \frac{3\,L^3}{8} \;-\; \frac{3\,L^3}{12}\,\ln(L).

Reemplazamos:

(3L383L312ln(L))  +  C3  =  L3(259618ln(L)).\Bigl(\frac{3\,L^3}{8} - \frac{3\,L^3}{12}\,\ln(L)\Bigr) \;+\; C_{3} \;=\; L^3\Bigl(\tfrac{25}{96} - \tfrac{1}{8}\ln(L)\Bigr).

Aislamos C3C_{3}:

C3=L3(259618ln(L))    3L38  +  3L312ln(L).C_{3} = L^3\Bigl(\tfrac{25}{96} - \tfrac{1}{8}\ln(L)\Bigr) \;-\; \frac{3\,L^3}{8} \;+\; \frac{3\,L^3}{12}\,\ln(L).

Agrupamos por un lado lo que depende de ln(L)\ln(L) y por otro lado los términos numéricos:

  • Términos con ln(L)\ln(L):

    (\displaystyle -,L^3;\tfrac{1}{8},\ln(L) ;+; \frac{3,L^3}{12},\ln(L)

    L^3,\ln(L)\Bigl(-\tfrac{1}{8} + \tfrac{1}{4}\Bigr)

    L^3,\ln(L),\tfrac{1}{8}. )

  • Términos sin ln(L)\ln(L):

    (\displaystyle L^3;\tfrac{25}{96} ;-; \frac{3,L^3}{8}

    L^3\Bigl(\tfrac{25}{96} - \tfrac{36}{96}\Bigr)

    L^3\Bigl(-\tfrac{11}{96}\Bigr)

    -,\frac{11,L^3}{96}. )

Por tanto,

C3=11L396  +  L38ln(L).\boxed{ C_{3} = -\,\frac{11\,L^3}{96} \;+\; \frac{L^3}{8}\,\ln(L). }

y

C1=9L4    frac3L2ln(L),C2=0,C4=0.\boxed{ C_{1} = \frac{9\,L}{4} \;-\; frac{3\,L}{2}\,\ln(L), \quad C_{2}=0, \quad C_{4}=0. }

4. Expresión final de y(x)y(x)

Sustituyendo las constantes halladas:

y(x)=x48ln(x)    25x496  +  (9L43L2ln(L))=C1  x36  +  (11L396+L38ln(L))=C3  x.y(x) = \frac{x^4}{8}\,\ln(x) \;-\; \frac{25\,x^4}{96} \;+\; \underbrace{\bigl(\tfrac{9\,L}{4} - \tfrac{3\,L}{2}\ln(L)\bigr)}_{=\,C_{1}} \;\frac{x^3}{6} \;+\; \underbrace{\bigl(-\tfrac{11\,L^3}{96} + \tfrac{L^3}{8}\ln(L)\bigr)}_{=\,C_{3}} \;x.

Para mayor claridad, puede factorizarse un LL, etc. Pero ya tenemos la función de deflexión que cumple las cuatro condiciones de apoyo. En notación más compacta:

y(x)=x48ln(x)    25x496  +  (9L43L2ln(L))x36  +  (11L396+L38ln(L))x.\boxed{ y(x) = \frac{x^4}{8}\,\ln(x) \;-\; \frac{25\,x^4}{96} \;+\; \Bigl(\frac{9\,L}{4} - \frac{3\,L}{2}\,\ln(L)\Bigr)\,\frac{x^3}{6} \;+\; \Bigl(- \frac{11\,L^3}{96} + \frac{L^3}{8}\,\ln(L)\Bigr)\,x. }

Si quisiéramos incluir el factor 1EI\tfrac{1}{E\,I} (que aparece en la ecuación diferencial real EIy(x)=w(x)\,E\,I\,y''''(x)=w(x)), habría que multiplicar todos los términos que vengan de integrar w(x)w(x) por 1EI\tfrac{1}{E\,I}. Las constantes de integración se ajustarían de igual manera para satisfacer las condiciones de contorno. En otras palabras, el resultado final sería

y(x)=1EI  [  (la misma expresioˊn obtenida arriba)].y(x) = \frac{1}{E\,I} \;\Bigl[\;\text{(la misma expresión obtenida arriba)}\Bigr].

5. Cálculo de la flecha máxima

Para hallar el máximo de la flecha en el intervalo (0,L)(0,L), se necesita:

  1. Calcular la derivada y(x)y'(x) (la pendiente de la elástica).
  2. Resolver la ecuación y(x)=0y'(x)=0 dentro de (0,L)(0,L).
  3. Verificar cuál raíz da el valor máximo (en lugar de un mínimo local, etc.), y comparar también con los extremos si fuera necesario (aunque en los extremos la flecha es cero).

La derivada que ya obtuvimos (con las constantes incluidas) es:

y(x)  =  x32ln(x)    11x312  +  C1x22  +  C3,y'(x) \;=\; \frac{x^3}{2}\,\ln(x) \;-\; \frac{11\,x^3}{12} \;+\; \frac{C_{1}\,x^2}{2} \;+\; C_{3},

con

C1=9L43L2ln(L),C3=11L396  +  L38ln(L).C_{1} = \frac{9\,L}{4} - \frac{3\,L}{2}\,\ln(L), \qquad C_{3} = -\,\frac{11\,L^3}{96} \;+\; \frac{L^3}{8}\,\ln(L).

Por tanto, la condición de extremal es

y(x)  =  0x32ln(x)    11x312  +  C1x22  +  C3  =  0.y'(x) \;=\; 0 \quad\Longrightarrow\quad \frac{x^3}{2}\,\ln(x) \;-\; \frac{11\,x^3}{12} \;+\; \frac{C_{1}\,x^2}{2} \;+\; C_{3} \;=\; 0.

Es una ecuación transcendental (contiene x3ln(x)x^3\ln(x), etc.). En la práctica, no existe una solución cerrada elemental. Lo usual es:

  1. No hay simetría (la carga w(x)=3ln(x)w(x)=3\ln(x) no es simétrica), así que no basta con mirar x=L/2x=L/2.
  2. Se introduce una variable adimensional α=x/L\alpha = x/L.
  3. Se reescribe la ecuación anterior en términos de α\alpha y ln(α)\ln(\alpha) (más ln(L)\ln(L)), y se resuelve numéricamente (p.ej. con Newton-Raphson) para α(0,1)\alpha\in(0,1).

Hecha esa búsqueda numérica, se hallará una αmax\alpha_{\mathrm{max}} que da el punto de máxima flecha. Finalmente, se evalúa y(αmaxL)y(\alpha_{\mathrm{max}}\,L) en la fórmula de y(x)y(x).

En muchos problemas de cargas “clásicas” (uniforme, triangular, etc.), se tienen fórmulas cerradas para la flecha máxima. En un caso como éste log(x3)\log(x^3), la forma de la carga obliga a resolver una ecuación trascendental.


6. Comentario sobre la sección pentagonal

El que la sección transversal sea un pentágono de acero afecta principalmente a:

  1. El cálculo del momento de inercia II, que se usaría en la fórmula.
  2. El módulo de Young EE.

La forma general de la solución de la viga no cambia. En otras palabras, toda la integración anterior es exactamente la misma:

  • Simplemente, en la expresión final de y(x)y(x) aparecerá un factor global 1EI\tfrac{1}{E\,I} (con II calculado específicamente para el pentágono).
  • Si uno quiere la flecha en unidades físicas (metros), habrá que incorporar 1EI\frac{1}{E\,I} y la magnitud real de la carga (en N/m, etc.).

7. Resumen de la Solución

  1. Ecuación de flexión:

    EIy(x)=3ln(x).E\,I\,y''''(x) = 3\,\ln(x).
  2. Solución general tras integrar 4 veces (omitiendo la constante 1/(EI)1/(E\,I) para simplificar en cada paso) y aplicar las condiciones de apoyo y(0)=y(L)=0,  y(0)=y(L)=0y(0)=y(L)=0,\;y''(0)=y''(L)=0, da:

    y(x)  =  x48ln(x)    25x496parte por la carga  +  (9L43L2ln(L))x36  +  (11L396+L38ln(L))xajuste para satisfacer BC.y(x) \;=\; \underbrace{\frac{x^4}{8}\,\ln(x) \;-\; \frac{25\,x^4}{96}}_{\text{parte por la carga}} \;+\; \underbrace{\Bigl(\frac{9\,L}{4} - \frac{3\,L}{2}\ln(L)\Bigr)\,\frac{x^3}{6} \;+\; \Bigl(- \frac{11\,L^3}{96} + \frac{L^3}{8}\,\ln(L)\Bigr)\,x}_{\text{ajuste para satisfacer BC}}.

    En rigor físico, la fórmula final completa sería

    y(x)=1EI[x48ln(x)25x496+(9L43L2ln(L))x36+(11L396+L38ln(L))x].\boxed{ y(x) = \frac{1}{E\,I} \,\Biggl[ \frac{x^4}{8}\,\ln(x) - \frac{25\,x^4}{96} + \Bigl(\frac{9\,L}{4} - \frac{3\,L}{2}\ln(L)\Bigr)\,\frac{x^3}{6} + \Bigl(- \frac{11\,L^3}{96} + \frac{L^3}{8}\,\ln(L)\Bigr)\,x \Biggr]. }
  3. Flecha máxima:

    • Se determina numéricamente, resolviendo y(x)=0y'(x)=0 en el intervalo (0,L)(0,L).
    • Se comprueba cuál raíz produce el máximo (y no un mínimo local).
    • Se evalúa y(x)y(x) en dicho valor.

No existe una expresión analítica sencilla (tipo “máximo en x=L/2x=L/2” o parecida) porque la carga w(x)=3ln(x)w(x)=3\,\ln(x) no es ni simétrica ni polinómica simple. Por tanto, el resultado final para la flecha máxima se obtiene de:

xmax  tal que  y(xmax)=0,0<xmax<L,(resolviendo numeˊricamente),x_{\max} \;\text{tal que}\; y'(x_{\max})=0, \quad 0 < x_{\max}< L, \quad \text{(resolviendo numéricamente)},

y luego

δmax  =  y(xmax).\delta_{\max} \;=\; y\bigl(x_{\max}\bigr).

Conclusión

  1. Se ha derivado la ecuación de la flecha elástica y(x)y(x) para la viga biapoyada con carga w(x)=3ln(x)w(x)=3\ln(x).
  2. Se han hallado las constantes de integración imponiendo las condiciones de contorno (flecha nula y momento nulo en x=0x=0 y x=Lx=L).
  3. La flecha máxima se localiza resolviendo y(x)=0y'(x)=0 numéricamente, pues la ecuación que surge es trascendental y no tiene solución cerrada simple.

El hecho de que la sección sea pentagonal de acero afecta a EE e II, pero no modifica la estructura de la solución; simplemente, en la respuesta final, se multiplica por el factor 1EI\tfrac{1}{E\,I} (donde II es el momento de inercia de la sección pentagonal) para obtener valores de flecha en unidades de longitud.

Share this Q&A