Analysis of beams on elastic foundations using Green's functions

Análisis de vigas sobre fundación flexible empleando funciones de Green Analysis of beams on elastic foundations using Green's functions

Juan Camilo Molina-Villegas^a,b, Jorge Eliecer Ballesteros Ortega^b, Andrés Camilo Quintero Toro^a

a Universidad de Medellín - Facultad de ingenierías.

b Universidad Nacional de Colombia - Sede Medellín - Facultad de minas - Departamento de ingeniería civil.

Resumen

Las vigas sobre fundación flexible representan un modelo básico dentro del análisis estructural, las cueles suelen emplearse para modelar vigas de cimentación, pilas, muros de contención y estructuras más complejas que tengan algún tipo de estos elementos. Para el análisis de estas, usualmente se emplea el método de elementos finitos [1], el cual produce una solución aproximada del problema; y el método de rigidez con funciones de Green [2] y [3], el cual produce una solución exacta.

En este artículo se presenta una metodología 100 % basada en el empleo de las funciones de Green (respuesta ante una fuerza puntual unitaria), para obtener la respuesta exacta de las vigas sobre fundación flexible. La principal ventaja de esta formulación es que es mucho menos costosa computacionalmente que las citadas alternativas e incluso para gran cantidad de problemas, esta puede expresarse solo por medio de sumas e integrales, las cuales se pueden realizar fácilmente de forma numérica.

Por completez, también se presentan una gran variedad de funciones de Green para vigas sobre fundación flexible finitas e infinitas con diferentes condiciones de frontera, así como algunos ejemplos con la implementación de la metodología propuesta.

Beams on elastic foundation are basic elements within structural analysis, which are used to model foundation beams, foundation piles, retaining walls, and more complex structures that include some of these elements. For their analysis, the finite element method is usually used [1], which produces an approximate solution of the problem; and the Green's function stiffness method [2] and [3], which produces an exact solution. This article presents a methodology 100% based on the use of Green function's (response to a unit point force), to obtain the exact response of beams on flexible foundation. The main advantage of this formulation is its computational low cost compared to the aforementioned alternatives, and even for a large number of problems, it can be expressed only by means of sums and integrals, which can be easily performed numerically.

Also, a great variety of Green function's for finite and infinite beams on elastic foundations with different boundary conditions are also presented, as well as some examples with the implementation of the proposed methodology.

keywords Funciones de Green, vigas sobre fundación flexible, pilas, campos de desplazamiento.

1 Introducción

En mecánica, se define como función de Green, al campo de desplazamiento en un medio ante la acción de una fuerza puntual unitaria. En general, para casos tridimensionales, dicha fuerza unitaria puede estar ubicada en cualquier punto y dirección del medio, y el campo de desplazamiento será vectorial y tendrá tres componentes escalares. El concepto de función de Green puede fácilmente extenderse a cualquier problema físico en el cual la fuerza puntual será reemplazada por una fuente o acción puntual unitaria, mientras que el campo de desplazamiento será reemplazado por el campo de la variable dependiente principal para el problema. Por ejemplo, en transferencia de calor, la función de Green será el campo de temperatura debido a la acción de una fuente de calor puntual unitaria.

Por si solas las funciones de Green juegan un papel muy importante en la solución de problemas físicos puesto que presentan la solución a problemas fundamentales. Por ejemplo, en el famoso problema de Boussinesq, la función de Green es el campo de desplazamiento en un semiespacio (espacio semi-infinito) tridimensional ante la acción de una fuerza unitaria en su superficie. En textos como [4], [5] y [6] se presentan las principales funciones de Green empleadas en la física clásica, en [7] se presenta un compendio de las principales funciones de Green empleadas en la geotecnia, para sismología y elastodinámica sus principales funciones de Green son presentadas en [8], para transferencia de calor en [9], mientras que aquellas propias de los problemas de difusión se presentan en [10].

Además de lo anterior, las funciones de Green son la base de algunos métodos numéricos, los cuales permiten resolver problemas diferentes a aquellos para los cuales fueron definidas dichas funciones, es decir, problemas con diferentes condiciones de cargas o fuentes y de frontera. Entre estos métodos numéricos destacan el método directo de elementos de frontera [11], [12] y [13], el método indirecto de elementos de frontera [14] y [15], y el método de rigidez con funciones de Green [2] y [3].

Las vigas sobre fundación flexible representan un modelo básico dentro del análisis estructural, las cueles suelen emplearse para modelar vigas de cimentación, pilas, muros de contención y estructuras más complejas que tengan algún tipo de estos elementos. Pese a que en la actualidad existen diversas metodologías para el análisis de vigas sobre fundación flexible, las cuales van desde la solución directa de los problemas de valor en la frontera gobernantes empleando métodos clásicos de ecuaciones diferenciales [16] y [3], hasta el método de elementos finitos [1] o de rigidez con funciones de Green [2], es posible solo mediante el empleo de funciones de Green obtener su solución.

En este artículo se presenta una metodología alternativa basada 100 % en el empleo de las funciones de Green para el análisis de vigas sobre fundación flexible, la cual en muchos casos puede ser bastante menos costosa computacionalmente que las demás empleadas usualmente y que teóricamente obtiene la solución exacta del modelo, es decir, no se trata de una aproximación.

2 Ecuación diferencial gobernante y convenciones

La ecuación diferencial que gobierna el desplazamiento de una viga sobre fundación flexible (ver la figura 1), formada por una viga de material elástico lineal homogéneo y sección transversal constante, apoyada sobre un medio flexible con rigidez por unidad de longitud constante, es (ver [2] para una deducción paso a paso):

$EI{\frac {d^{4}v}{dx^{4}}}(x)+kv(x)=q(x)$

(1)

Donde:

${\textstyle E}$: Módulo de elasticidad del material de la viga.

${\textstyle I}$: Momento de inercia respecto al eje ${\textstyle z}$ de la sección transversal de la viga.

${\textstyle k}$: Constante de rigidez por unidad de longitud del suelo de soporte de la viga.

${\textstyle v(x)}$: Campo de desplazamiento en dirección del eje ${\textstyle y}$ de la viga.

${\textstyle q(x)}$: Fuerza externa por unidad de longitud sobre la viga.

Figura 1: Viga sobre fundación flexible sometida a carga externa distribuida.)

Por facilidad en su solución, es usual dividir (1) entre ${\textstyle EI}$, obteniéndose:

${\frac {d^{4}v}{dx^{4}}}(x)+4\lambda ^{4}v(x)={\frac {q(x)}{EI}}$

(2)

Donde: ${\textstyle \lambda ={\sqrt[{4}]{\dfrac {k}{4EI}}}}$.

Además de la fuerza externa por unidad de longitud ${\textstyle q(x)}$, la viga también se encuentra sometida a la fuerza distribuida que el suelo ejerce sobre esta, la cual es proporcional al desplazamiento en cada punto de la viga y se calcula como:

$f_{S}(x)=-kv(x)$

(3)

En este punto es importante resaltar que pese a que las ecuaciones (1) y (2) están escritas en términos de la fuerza externa por unidad de longitud ${\textstyle (q(x))}$ sobre la viga, es decir una carga distribuida, incluso fuerzas y momentos puntuales se pueden expresar de esta manera. Como ejemplo se tiene que la carga externa puntual aplicada a la viga sobre fundación flexible presentada en las figura 2a se definen por medio de la siguiente función:

$q(x,\xi )=F\cdot \delta (x-\xi )$

(4)

Donde ${\textstyle \delta ()}$ es la función delta de Dirac, cuyas principales propiedades se presentan en [5].

(a) Fuerza puntual.	(b) Momento puntual.
Figura 2: Vigas sobre fundación flexible sometidas a fuerzas y momentos puntuales.

Mientras que para el caso del momento externo puntual presentado en la figura (2b), la función ${\textstyle q(x)}$ se define a partir de (4), como (ver la figura (3) la una explicación):

$q(x,\xi )=\lim _{\Delta \xi \rightarrow 0}\left[{\dfrac {M}{\Delta \xi }}\delta (x-\xi {-\Delta }\xi )-{\dfrac {M}{\Delta \xi }}\delta (x-\xi )\right]=M\lim _{\Delta \xi \rightarrow 0}{\dfrac {\delta (x-\xi {-\Delta }\xi )-\delta (x-\xi )}{\Delta \xi }}$ $=M{\dfrac {\partial \delta (x-\xi )}{\partial \xi }}$

(5)

Figura 3: Definición de un momento puntual como el caso límite de dos fuerzas puntuales muy cercanas entre si.)

De otra parte, a partir del campo de desplazamiento, es posible obtener los campos de fuerzas internas en la viga como:

$M(x)=EI{\dfrac {d^{2}v}{dx^{2}}}(x)$ (6.a) $V(x)=-EI{\dfrac {d^{3}v}{dx^{3}}}(x)$ (6.b)

Donde:

${\textstyle M(x)}$: Campo de momento flector, cuya convención positiva presentada en la figura 4.

${\textstyle V(x)}$: Campo de fuerza cortante, cuya convención positiva presentada en la figura 4.

Figura 4: Convención positiva de las fuerzas internas.

3 Funciones de Green para vigas sobre fundación flexible

3.1 Fuerzas puntuales

En las figuras 5 y 6 se presentan los modelos de las vigas sobre fundación flexible finitas e infinitas para los cuales se definirán las funciones de Green a emplear en este artículo, todas las cuales son solución de la siguiente ecuación diferencial:

${\frac {\partial ^{4}G_{yy}}{\partial x^{4}}}(x,\xi )+4\lambda ^{4}G_{yy}(x,\xi )={\frac {\delta (x-\xi )}{EI}}$

(7)

Donde ${\textstyle G_{yy}(x,\xi )}$ es la función de Green o campo de desplazamiento del problema en particular.

La ecuación (7) es un caso particular de (2), en la cual en lugar de ${\textstyle v(x,\xi )}$ se ha empleado ${\textstyle G_{yy}(x,\xi )}$ y para definir la función de la carga externa se ha empleado (4) con ${\textstyle F=1}$.

(a) Libre - libre.	(b) Simplemente apoyada.
(c) Empotrada - libre.	(d) Empotrada - empotrada.
Figura 5: Vigas sobre fundación flexible finitas con diferentes condiciones de frontera y sometidas a fuerzas puntuales unitarias.

Figura 6: Viga sobre fundación flexible infinita sometida a una fuerza puntual unitaria.

3.1.1 Viga sobre fundación flexible finita libre - libre

La función de Green para la viga sobre fundación flexible finita, con ambos extremos libres y presentada en la figura 5a, es la solución del siguiente problema de valor en la frontera:

${\frac {\partial ^{4}G_{yy}^{LL}}{\partial x^{4}}}(x,\xi )+4\lambda ^{4}G_{yy}^{LL}(x,\xi )={\frac {\delta (x-\xi )}{EI}}$ (8.a) $EI{\dfrac {\partial ^{2}G_{yy}^{LL}}{\partial ^{2}x}}(0,\xi )=0$ (8.b) $-EI{\dfrac {\partial ^{3}G_{yy}^{LL}}{\partial ^{3}x}}(0,\xi )=0$ (8.c) $EI{\dfrac {\partial ^{2}G_{yy}^{LL}}{\partial ^{2}x}}(L,\xi )=0$ (8.d) $-EI{\dfrac {\partial ^{3}G_{yy}^{LL}}{\partial ^{3}x}}(L,\xi )=0$ (8.e)

Donde las dos ${\textstyle L}$ que acompañan al nombre de la función de Green hacen referencia a que cada uno de los extremos es ${\textstyle L}$ibre.

Cada uno de los tramos de esta función de Green se definen como:

$G_{yy}^{LL}(x,\xi )={\begin{cases}G_{yy}^{LL1}(x,\xi )&0<x\leq \xi \\G_{yy}^{LL2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(9)

Donde:

$G_{yy}^{LL1}(x,\xi )={\dfrac {A(x,\xi )}{4EI\lambda ^{3}\left[\sin ^{2}(\lambda L)-\sinh ^{2}(\lambda L)\right]}}$ (10.a) $G_{yy}^{LL2}(x,\xi )={\dfrac {B(x,\xi )}{4EI\lambda ^{3}\left[\sin ^{2}(\lambda L)-\sinh ^{2}(\lambda L)\right]}}$ (10.b)

$A(x,\xi )=2A_{1}(\xi )\cos(\lambda x)\cosh(\lambda x)+A_{2}(\xi )\left[\sin(\lambda x)\cosh(\lambda x)+\cos(\lambda x)\sinh(\lambda x)\right]$

(11)

$A_{1}(\xi )=\sin(\lambda L)\cosh(\lambda \xi )\cos[\lambda (L-\xi )]-\sinh(\lambda L)\cos(\lambda \xi )\cosh[\lambda (L-\xi )]$

(12)

$A_{2}(\xi )=\sin(\lambda L)\left\{\cosh(\lambda \xi )\sin[\lambda (L-\xi )]-\sinh(\lambda \xi )\cos[\lambda (L-\xi )]\right\}$ $+\sinh(\lambda L)\left\{\cos(\lambda \xi )\sinh[\lambda (L-\xi )](\lambda \xi )\cosh[\lambda (L-\xi )]\right\}$

$B(x,\xi )=A(\xi ,x)$

(13)

$G_{yy}^{LL2}(x,\xi )=G_{yy}^{LL1}(\xi ,x)$

(14)

3.1.2 Viga sobre fundación flexible finita simplemente apoyada

La función de Green para la viga sobre fundación flexible finita con ambos extremos libres, presentada en la figura 5b, es la solución del siguiente problema de valor en la frontera:

${\frac {\partial ^{4}G_{yy}^{SS}}{\partial x^{4}}}(x,\xi )+4\lambda ^{4}G_{yy}^{SS}(x,\xi )={\frac {\delta (x-\xi )}{EI}}$ (15.a) $G_{yy}^{SS}(0,\xi )=0$ (15.b) $EI{\dfrac {\partial ^{2}G_{yy}^{SS}}{\partial ^{2}x}}(0,\xi )=0$ (15.c) $G_{yy}^{SS}(L,\xi )=0$ (15.d) $EI{\dfrac {\partial ^{2}G_{yy}^{SS}}{\partial ^{2}x}}(L,\xi )=0$ (15.e)

Donde las dos ${\textstyle S}$ que acompañan al nombre de la función de Green hacen referencia a que cada uno de los extremos es un apoyo ${\textstyle S}$imple.

Cada uno de los tramos de esta función de Green se definen como:

$G_{yy}^{SS}(x,\xi )={\begin{cases}G_{yy}^{SS1}(x,\xi )&0<x\leq \xi \\G_{yy}^{SS2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(16)

Donde:

$G_{yy}^{SS1}(x,\xi )={\dfrac {C(x,\xi )}{4EI\lambda ^{3}\left[\cos(2\lambda L)-\cosh(2\lambda L)\right]}}$ (17.a) $G_{yy}^{SS2}(x,\xi )={\dfrac {D(x,\xi )}{4EI\lambda ^{3}\left[\cos(2\lambda L)-\cosh(2\lambda L)\right]}}$ (17.b)

$C(x,\xi )=\left\{\sin[\lambda (2L-\xi )]\sinh(\lambda \xi )-\sinh[\lambda (2L-\xi )]\sin(\lambda \xi )\right\}\left[\sin(\lambda x)\cosh(\lambda x)+\cos(\lambda x)\sinh(\lambda x)\right]$ $+\left\{\cosh[\lambda (2L-\xi )]\cos(\lambda \xi )-\cos[\lambda (2L-\xi )]\cosh(\lambda \xi )\right\}\left[\sin(\lambda x)\cosh(\lambda x)-\cos(\lambda x)\sinh(\lambda x)\right]$

(18)

$D(x,\xi )=C(\xi ,x)$

(19)

$G_{yy}^{SS2}(x,\xi )=G_{yy}^{SS1}(\xi ,x)$

(20)

3.1.3 Viga sobre fundación flexible finita empotrada - libre

La función de Green para la viga sobre fundación flexible finita con su extremo inicial ${\textstyle (x=0)}$ empotrado y su extremo final ${\textstyle (x=L)}$ libre, presentada en la figura 5c, es la solución del siguiente problema de valor en la frontera:

${\frac {\partial ^{4}G_{yy}^{EL}}{\partial x^{4}}}(x,\xi )+4\lambda ^{4}G_{yy}^{EL}(x,\xi )={\frac {\delta (x-\xi )}{EI}}$ (21.a) $G_{yy}^{EL}(0,\xi )=0$ (21.b) ${\dfrac {\partial G_{yy}^{EL}}{\partial x}}(0,\xi )=0$ (21.c) $EI{\dfrac {\partial ^{2}G_{yy}^{EL}}{\partial ^{2}x}}(L,\xi )=0$ (21.d) $-EI{\dfrac {\partial ^{3}G_{yy}^{EL}}{\partial ^{3}x}}(L,\xi )=0$ (21.e)

Donde la ${\textstyle E}$ del primer supeindice del nombre hace referencia a que el extremo inicial es ${\textstyle E}$mpotrado y la ${\textstyle L}$ a que el extremo final es ${\textstyle L}$ibre.

Cada uno de los tramos de esta función de Green se definen como:

$G_{yy}^{EL}(x,\xi )={\begin{cases}G_{yy}^{EL1}(x,\xi )&0<x\leq \xi \\G_{yy}^{EL2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(22)

Donde:

$G_{yy}^{EL1}(x,\xi )={\dfrac {E(x,\xi )}{4EI\lambda ^{3}\left[\cos ^{2}(\lambda L)+\cosh ^{2}(\lambda L)\right]}}$ (23.a) $G_{yy}^{EL2}(x,\xi )={\dfrac {F(x,\xi )}{4EI\lambda ^{3}\left[\cos ^{2}(\lambda L)+\cosh ^{2}(\lambda L)\right]}}$ (23.b)

$E(x,\xi )=2E_{1}(\xi )\sin(\lambda x)\sinh(\lambda x)+E_{2}(\xi )\left[\sin(\lambda x)\cosh(\lambda x)-\cos(\lambda x)\sinh(\lambda x)\right]$

(24)

$E_{1}(\xi )=\sinh(\lambda L)\sin(\lambda \xi )\sinh[\lambda (L-\xi )]-\sin(\lambda L)\sinh(\lambda \xi )\sin[\lambda (L-\xi )]$ $+\sin(\lambda \xi )\cosh(\lambda \xi )+\cos(\lambda \xi )\sinh(\lambda \xi )$

(25)

$E_{2}(\xi )=\sin(\lambda L)\left\{\cosh(\lambda \xi )\sin[\lambda (L-\xi )+\sinh(\lambda \xi )\cos[\lambda (L-\xi )]\right\}$ $-\sinh(\lambda L)\left\{\sin(\lambda \xi )\cosh[\lambda (L-\xi )]+\cos(\lambda \xi )\sinh[\lambda (L-\xi )]\right\}-2\cos(\lambda \xi )\cosh(\lambda \xi )$

(26)

$F(x,\xi )=E(\xi ,x)$

(27)

$G_{yy}^{EL2}(x,\xi )=G_{yy}^{EL1}(\xi ,x)$

(28)

3.1.4 Viga sobre fundación flexible finita doblemente empotrada

La función de Green para la viga sobre fundación flexible finita doblemente empotrada presentada en la figura 5d, es la solución del siguiente problema de valor en la frontera:

${\frac {\partial ^{4}G_{yy}^{EE}}{\partial x^{4}}}(x,\xi )+4\lambda ^{4}G_{yy}^{EE}(x,\xi )={\frac {\delta (x-\xi )}{EI}}$ (29.a) $G_{yy}^{EE}(0,\xi )=0$ (29.b) ${\dfrac {\partial G_{yy}^{EE}}{\partial x}}(0,\xi )=0$ (29.c) $G_{yy}^{EE}(L,\xi )=0$ (29.d) ${\dfrac {\partial G_{yy}^{EE}}{\partial x}}(L,\xi )=0$ (29.e)

Donde las dos ${\textstyle E}$ que acompañan al nombre de la función de Green hacen referencia a que cada uno de los extremos es un apoyo ${\textstyle E}$mpotrado.

Cada uno de los tramos de esta función de Green se definen como:

$G_{yy}^{EE}(x,\xi )={\begin{cases}G_{yy}^{EE1}(x,\xi )&0<x\leq \xi \\G_{yy}^{EE2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(30)

Donde:

$G_{yy}^{EE1}(x,\xi )={\dfrac {H(x,\xi )}{4EI\lambda ^{3}[\sin ^{2}(\lambda L)-\sinh ^{2}(\lambda L)]}}$ (31.a) $G_{yy}^{EE2}(x,\xi )={\dfrac {J(x,\xi )}{4EI\lambda ^{3}[\sin ^{2}(\lambda L)-\sinh ^{2}(\lambda L)]}}$ (31.b)

$H(\xi )=2\sin(\lambda x)\sinh(\lambda x)H_{1}(\xi )+[\sin(\lambda x)\cosh(\lambda x)-\cos(\lambda x)\sinh(\lambda x)]H_{2}(\xi )$

(32)

$H_{1}(\xi )=\sin(\lambda L)\sinh(\lambda \xi )\sin[\lambda (L-\xi )]-\sin(\lambda \xi )\sinh(\lambda L)\sinh[\lambda (L-\xi )]$

(33)

$H_{2}(\xi )=-\sin(\lambda L)\{\cosh(\lambda \xi )\sin[\lambda (L-\xi )]+\sinh(\lambda \xi )\cos[\lambda (L-\xi )]\}$ $+\sinh(\lambda L)\{\sin(\lambda \xi )\cosh[\lambda (L-\xi )]+\cos(\lambda \xi )\sinh[\lambda (L-\xi )]\}$

(34)

$J(x,\xi )=H(\xi ,x)$

(35)

$G_{yy}^{EE2}(x,\xi )=G_{yy}^{EE1}(\xi ,x)$

(36)

3.1.5 Viga sobre fundación flexible infinita

La función de Green para la viga sobre fundación flexible infinita presentada en la figura 6, se define como:

$G_{yy}(x,\xi )={\begin{cases}G_{yy}^{1}(x,\xi )&-\infty <x\leq \xi \\G_{yy}^{2}(x,\xi )\qquad &\xi \leq x<+\infty \end{cases}}$

(37)

$G_{yy}^{1}(x,\xi )={\dfrac {\exp =+\lambda (x-\xi )\{\cos =\lambda (x-\xi )}{8EI\lambda ^{3}}}$

(38)

$G_{yy}^{2}(x,\xi )={\dfrac {\exp =-\lambda (x-\xi )\{\cos =\lambda (x-\xi )\}}{8EI\lambda ^{3}}}$

(39)

La cual se puede escribir en un solo tramo desde ${\textstyle x=-\infty }$ a ${\textstyle x=+\infty }$ como:

$G_{yy}(x,\xi )={\dfrac {\exp(-\lambda |x-\xi |)[\cos(\lambda |x-\xi |)+\sin(\lambda |x-\xi |)]}{8EI\lambda ^{3}}}\qquad -\infty {<}x<+\infty$

(40)

3.2 Momentos puntuales

En las figuras 7 y 8 se presentan las vigas finitas e infinitas para las cuales se definirán las funciones de Green debidas a la aplicación de un momento puntual unitario, todas las cuales se pueden obtener directamente a partir de su contraparte sometida a fuerza puntual unitaria presentada en la sección 3.1.

(a) Libre - libre finita.	(b) Simplemente apoyada.
(c) Empotrada - libre.	(d) Empotrada - libre infinita.
Figura 7: Vigas sobre fundación flexible finitas con diferentes condiciones de frontera y sometidas a momentos puntuales unitarios.

Figura 8: Viga sobre fundación flexible infinita sometida a un momento puntual unitario.

3.2.1 Viga sobre fundación flexible finita libre - libre

A partir de (9), la función de Green para la viga sobre fundación flexible finita con ambos extremos libres presentada en la figura 7a, se define como:

$G_{y\theta }^{LL}(x,\xi )={\dfrac {G_{yy}^{LL}}{\partial \xi }}(x,\xi )$

(41)

Y sus tramos se definen como:

$G_{y\theta }^{LL}(x,\xi )={\begin{cases}G_{y\theta }^{LL1}(x,\xi )={\dfrac {G_{yy}^{LL1}}{\partial \xi }}(x,\xi )&0<x\leq \xi \\[0.3cm]G_{y\theta }^{LL2}(x,\xi )={\dfrac {G_{yy}^{LL2}}{\partial \xi }}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(42)

3.2.2 Viga sobre fundación flexible finita simplemente apoyada

A partir de (16), la función de Green para la viga sobre fundación flexible finita simplemente apoyada presentada en la figura 7b, se define como:

$G_{y\theta }^{SS}(x,\xi )={\dfrac {G_{yy}^{SS}}{\partial \xi }}(x,\xi )$

(43)

Y sus tramos se definen como:

$G_{y\theta }^{SS}(x,\xi )={\begin{cases}G_{y\theta }^{SS1}(x,\xi )={\dfrac {G_{yy}^{SS1}}{\partial \xi }}(x,\xi )&0<x\leq \xi \\[0.3cm]G_{y\theta }^{SS2}(x,\xi )={\dfrac {G_{yy}^{SS2}}{\partial \xi }}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(44)

3.2.3 Viga sobre fundación flexible finita empotrada - libre

A partir de (22), la función de Green para la viga sobre fundación flexible finita simplemente apoyada presentada en la figura 7c, se define como:

$G_{y\theta }^{EL}(x,\xi )={\dfrac {G_{yy}^{EL}}{\partial \xi }}(x,\xi )$

(45)

Y sus tramos se definen como:

$G_{y\theta }^{EL}(x,\xi )={\begin{cases}G_{y\theta }^{EL1}(x,\xi )={\dfrac {G_{yy}^{EL1}}{\partial \xi }}(x,\xi )&0<x\leq \xi \\[0.3cm]G_{y\theta }^{EL2}(x,\xi )={\dfrac {G_{yy}^{EL2}}{\partial \xi }}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(46)

3.2.4 Viga sobre fundación flexible finita doblemente empotrada

A partir de (30), la función de Green para la viga sobre fundación flexible finita simplemente apoyada presentada en la figura 7d, se define como:

$G_{y\theta }^{EE}(x,\xi )={\dfrac {G_{yy}^{EE}}{\partial \xi }}(x,\xi )$

(47)

Y sus tramos se definen como:

$G_{y\theta }^{EE}(x,\xi )={\begin{cases}G_{y\theta }^{EE1}(x,\xi )={\dfrac {G_{yy}^{EE1}}{\partial \xi }}(x,\xi )&0<x\leq \xi \\[0.3cm]G_{y\theta }^{EE2}(x,\xi )={\dfrac {G_{yy}^{EE2}}{\partial \xi }}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(48)

3.2.5 Viga sobre fundación flexible infinita

A partir de (37), la función de Green para la viga sobre fundación flexible finita simplemente apoyada presentada en la figura 8, se define como:

$G_{y\theta }(x,\xi )={\dfrac {G_{yy}}{\partial \xi }}(x,\xi )$

(49)

Y sus tramos se definen como:

$G_{y\theta }(x,\xi )={\begin{cases}G_{y\theta }^{1}(x,\xi )={\dfrac {G_{yy}^{1}}{\partial \xi }}(x,\xi )&0<x\leq \xi \\[0.3cm]G_{y\theta }^{2}(x,\xi )={\dfrac {G_{yy}^{2}}{\partial \xi }}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

(50)

4 Solución de vigas sobre fundación flexible empleando funciones de Green

Antes de comenzar con la explicación de esta metodología es necesario tener claro que todas las cargas externas aplicadas en una viga sobre fundación flexible pueden obtenerse como un operador lineal de cargas puntuales unitarias, en particular sumas, integrales o derivadas ponderadas por un número o función. Para facilitar la explicación de este concepto, se empleará como ejemplo la viga sobre fundación flexible presentada en la figura 9.

Figura 9: Viga sobre fundación flexible sometida a cargas externas distribuidas y puntuales.

La carga puntual ${\textstyle F}$ se expresa como el producto de una fuerza puntual unitaria por el factor ${\textstyle F}$, es decir, se encuentra definida por la siguiente función (ver la ecuación (4) para una explicación):

$q_{1}(\xi )=F\cdot \delta (\xi {-}x_{4})$

(51)

Algo similar ocurre con las fuerzas puntuales debidas a las reacciones verticales que ocurren en los puntos 1 y 2, las cuales se definen respectivamente como:

$q_{2}(\xi )=FY_{1}\cdot \delta (\xi )$ (52.a) $q_{3}(\xi )=FY_{2}\cdot \delta (\xi {-}x_{2})$ (52.b)

Donde ${\textstyle FY_{1}}$ y ${\textstyle FY_{2}}$ son respectivamente los valores de las reacciones en los puntos 1 y 2, las cuales son positivas en dirección del eje ${\textstyle y}$.

Por su parte, los momentos puntuales debidos al momento externo ${\textstyle M}$ y la reacción a momento en el punto 1, se expresan respectivamente como:

$q_{4}(\xi )=M\cdot {\dfrac {\partial \delta (\xi {-}r)}{\partial r}}(r=x_{5})$ (53.a) $q_{5}(\xi )=M_{1}\cdot {\dfrac {\partial \delta (\xi {-}r)}{\partial r}}(r=0)$ (53.b)

Donde ${\textstyle M_{1}}$ es la reacción a momento en el punto 1 y ${\textstyle r}$ es una variable auxiliar muda o temporal.

Mientras que la carga distribuida se expresa como:

$q(\xi )=\int _{x_{1}}^{x_{3}}q(x)\delta (x-\xi )dx$

(54)

Ahora que es claro que cualquier carga puntual o distribuida puede expresarse como un operador lineal de cargas puntuales, se debe seleccionar una función de Green cuya viga sobre fundación flexible cumpla o pueda cumplir, por medio de la aplicación de fuerzas o momentos puntuales, todas condiciones de frontera y compatibilidad de la viga sobre fundación flexible en estudio. Para este problema la función de Green que más se acomoda es aquella de la viga empotrada libre presentada en la figura 5c. Una vez seleccionada la función de Green, la viga que la define debe hacerse equivalente a aquella en estudio, lo cual se logra a partir de la aplicación de cargas distribuidas y/o puntuales (todas ellas expresables como un operador lineal de cargas puntuales) y del cumplimiento de condiciones de frontera y/o compatibilidad. En este caso es claro que para que la estructura presentada en la figura (5c) sea equivalente a aquella de la figura 9 debe cumplirse que el desplazamiento en el punto 2 sea cero, es decir, ${\textstyle v_{2}=v(x_{2})=0}$, condición de compatibilidad empleada para obtener el valor de la reacción ${\textstyle FY_{2}}$ que es desconocido, la cual se expresa en términos de la función de Green ${\textstyle G^{EL}(z,\xi )}$, y empleando el principio de superposición, como:

$v_{2}=v(x_{2})=0=FY_{2}\cdot G_{yy}^{EL1}(x_{2},x_{2})+F\cdot G_{yy}^{EL1}(x_{2},x_{4})+M\cdot G_{y\theta }^{EL1}(x_{2},x_{5})$ $+\int _{x_{1}}^{x_{2}}q(\xi )G_{yy}^{EL2}(x_{2},\xi )d\xi {+\int }_{x_{2}}^{x_{3}}q(\xi )G_{yy}^{EL1}(x_{2},\xi )d\xi$

(55)

Cuyos términos se explican a continuación:

${\textstyle FY_{2}\cdot G_{yy}^{EL1}(x_{2},x_{2})}$: Desplazamiento en ${\textstyle x=x_{2}}$ debido a la fuerza puntual ${\textstyle FY_{2}}$ ubicada en ${\textstyle x=x_{2}}$.

${\textstyle F\cdot G_{yy}^{EL1}(x_{2},x_{4})}$: Desplazamiento en ${\textstyle x=x_{2}}$ debido a la fuerza puntual ${\textstyle F}$ ubicada en ${\textstyle x=x_{4}}$.

${\textstyle M_{2}\cdot G_{y\theta }^{EL1}(x_{2},x_{5})}$: Desplazamiento en ${\textstyle x=x_{2}}$ debido al momento puntual ${\textstyle M}$ ubicado en ${\textstyle x=x_{5}}$.

${\textstyle \int _{x_{1}}^{x_{2}}q(\xi )G_{yy}^{EL2}(x_{2},\xi )d\xi }$: Desplazamiento en ${\textstyle x=x_{2}}$ debido a la carga distribuida ${\textstyle q(\xi )}$ ubicada entre ${\textstyle x=x_{1}}$ y ${\textstyle x=x_{2}}$.

${\textstyle \int _{x_{2}}^{x_{3}}q(\xi )G_{yy}^{EL1}(x_{2},\xi )d\xi }$: Desplazamiento en ${\textstyle x=x_{2}}$ debido a la carga distribuida ${\textstyle q(\xi )}$ ubicada entre ${\textstyle x=x_{2}}$ y ${\textstyle x=x_{3}}$.

En este punto es importante recordarle al lector que ${\textstyle G_{yy}^{EL1}(x,\xi )}$ corresponde al desplazamiento antes de la fuerza puntual unitaria, mientras que ${\textstyle G_{yy}^{EL2}(x,\xi )}$ corresponde a aquel después de esta.

Como siguiente paso, una vez obtenidas todas las condiciones de compatibilidad, se resuelve el sistema lineal de ecuaciones que ellas cumplen, para así obtener el valor de las reacciones puntuales desconocidas, diferentes a aquellas que también ocurrirían en la función de Green. En este ejemplo esto equivale a obtener el único valor de ${\textstyle FY_{2}}$ que cumple (55).

Figura 10: Aplicación de las cargas externas de la viga sobre fundación flexible a analizar sobre aquella empleada para la definición de la función de Green a analizar.

Una vez obtenidas las reacciones que hacen que se cumplan las condiciones de compatibilidad y/o frontera, es necesario calcular el campo de desplazamiento en cada uno de los tramos de la viga sobre fundación flexible, lo cual se realiza de nuevo empleando el principio de superposición. Para este ejemplo el campo de desplazamiento tiene 6 tramos y se define como:

$v(x)={\begin{cases}v_{1}(x)&0<x\leq x_{1}\\v_{2}(x)&x_{1}\leq x\leq x_{2}\\v_{3}(x)&x_{2}\leq x\leq x_{3}\\v_{4}(x)&x_{3}\leq x\leq x_{4}\\v_{5}(x)&x_{4}\leq x\leq x_{5}\\v_{6}(x)&x_{5}\leq x<L\end{cases}}$

(56)

Mientras que cada uno de sus valores se calcula como:

$v_{1}(x)=FY_{2}\cdot G_{yy}^{EL1}(x,x_{2})+F\cdot G_{yy}^{EL1}(x,x_{4})+M\cdot G_{y\theta }^{EL1}(x,x_{5})+\int _{x_{1}}^{x_{3}}q(\xi )G_{yy}^{EL1}(x,\xi )d\xi$

(57)

$v_{2}(x)=FY_{2}\cdot G_{yy}^{EL1}(x,x_{2})+F\cdot G_{yy}^{EL1}(x,x_{4})+M\cdot G_{y\theta }^{EL1}(x,x_{5})$ $+\int _{x_{1}}^{x}q(\xi )G_{yy}^{EL2}(x,\xi )d\xi {+\int }_{x}^{x_{3}}q(\xi )G_{yy}^{EL1}(x,\xi )d\xi$

(58)

$v_{3}(x)=FY_{2}\cdot G_{yy}^{EL2}(x,x_{2})+F\cdot G_{yy}^{EL1}(x,x_{4})+M\cdot G_{y\theta }^{EL1}(x,x_{5})$ $+\int _{x_{1}}^{x}q(\xi )G_{yy}^{EL2}(x,\xi )d\xi {+\int }_{x}^{x_{3}}q(\xi )G_{yy}^{EL1}(x,\xi )d\xi$

(59)

$v_{4}(x)=FY_{2}\cdot G_{yy}^{EL2}(x,x_{2})+F\cdot G_{yy}^{EL1}(x,x_{4})+M\cdot G_{y\theta }^{EL1}(x,x_{5})+\int _{x_{1}}^{x_{3}}q(\xi )G_{yy}^{EL2}(x,\xi )d\xi$

(60)

$v_{5}(x)=FY_{2}\cdot G_{yy}^{EL2}(x,x_{2})+F\cdot G_{yy}^{EL2}(x,x_{4})+M\cdot G_{y\theta }^{EL1}(x,x_{5})+\int _{x_{1}}^{x_{3}}q(\xi )G_{yy}^{EL2}(x,\xi )d\xi$

(61)

$v_{6}(x)=FY_{2}\cdot G_{yy}^{EL2}(x,x_{2})+F\cdot G_{yy}^{EL2}(x,x_{4})+M\cdot G_{y\theta }^{EL2}(x,x_{5})+\int _{x_{1}}^{x_{3}}q(\xi )G_{yy}^{EL2}(x,\xi )d\xi$

(62)

Por último, a partir del campo de desplazamiento calculado anteriormente y empleando (3) y (6) se calculan los campos de la fuerza que el suelo ejerce sobre la viga, de momento flector y fuerza cortante respectivamente como:

$f_{S}(x)={\begin{cases}f_{S}^{1}(x)=-kv_{1}(x)&0<x\leq x_{1}\\f_{S}^{2}(x)=-kv_{2}(x)&x_{1}\leq x\leq x_{2}\\f_{S}^{3}(x)=-kv_{3}(x)&x_{2}\leq x\leq x_{3}\\f_{S}^{4}(x)=-kv_{4}(x)&x_{3}\leq x\leq x_{4}\\f_{S}^{5}(x)=-kv_{5}(x)&x_{4}\leq x\leq x_{5}\\f_{S}^{6}(x)=-kv_{6}(x)\qquad &x_{5}\leq x<L\end{cases}}$

(63)

$M(x)=EI{\dfrac {d^{2}v}{dx^{2}}}(x)={\begin{cases}M_{1}(x)=EI{\dfrac {d^{2}v_{1}}{dx^{2}}}(x)&0<x\leq x_{1}\\[0.2cm]M_{2}(x)=EI{\dfrac {d^{2}v_{2}}{dx^{2}}}(x)&x_{1}\leq x\leq x_{2}\\[0.2cm]M_{3}(x)=EI{\dfrac {d^{2}v_{3}}{dx^{2}}}(x)&x_{2}\leq x\leq x_{3}\\[0.2cm]M_{4}(x)=EI{\dfrac {d^{2}v_{4}}{dx^{2}}}(x)&x_{3}\leq x\leq x_{4}\\[0.2cm]M_{5}(x)=EI{\dfrac {d^{2}v_{5}}{dx^{2}}}(x)&x_{4}\leq x\leq x_{5}\\[0.2cm]M_{6}(x)=EI{\dfrac {d^{2}v_{6}}{dx^{2}}}(x)\qquad &x_{5}\leq x<L\end{cases}}$

(64)

$V(x)=-EI{\dfrac {d^{3}v}{dx^{3}}}(x)={\begin{cases}V_{1}(x)=-EI{\dfrac {d^{3}v_{1}}{dx^{3}}}(x)&0<x\leq x_{1}\\[0.2cm]V_{2}(x)=-EI{\dfrac {d^{3}v_{2}}{dx^{3}}}(x)&x_{1}\leq x\leq x_{2}\\[0.2cm]V_{3}(x)=-EI{\dfrac {d^{3}v_{3}}{dx^{3}}}(x)&x_{2}\leq x\leq x_{3}\\[0.2cm]V_{4}(x)=-EI{\dfrac {d^{3}v_{4}}{dx^{3}}}(x)&x_{3}\leq x\leq x_{4}\\[0.2cm]V_{5}(x)=-EI{\dfrac {d^{3}v_{5}}{dx^{3}}}(x)&x_{4}\leq x\leq x_{5}\\[0.2cm]V_{6}(x)=-EI{\dfrac {d^{3}v_{6}}{dx^{3}}}(x)\qquad &x_{5}\leq x<L\end{cases}}$

(65)

5 Ejemplos

5.1 Viga sobre fundación flexible simplemente apoyada sometida a una carga externa sinusoidal

Resolver la viga sobre fundación flexible presentada en la figura 11 empleando la función de Green de la figura 5b.

Figura 11: Viga sobre fundación flexible simplemente apoyada sometida a una carga externa sinusoidal.)

Solución

Funciones de Green a emplear

Como se indicó en el enunciado, la función de Green a emplear en la solución de este problema es ${\textstyle G_{yy}^{SS}(x,\xi )}$, la cual se presenta en la ecuación (16).

Cálculo del campo de desplazamiento

La carga externa distribuida sobre está definida por medio de la siguiente función:

$q(x)=-Q\sin \left(\pi {\dfrac {x}{L}}\right)\qquad 0<x<L$

(66)

Ahora, debido a que la viga sobre fundación flexible empleada para la función de Green tiene los mismos tramos, condiciones de frontera y de continuidad que aquella a analizar, su campo de desplazamiento se obtiene de forma inmediata solo mediante la siguiente integral:

$v(x)=\int _{0}^{L}G_{yy}^{SS}(x,\xi )q(\xi )d\xi =\int _{0}^{x}G_{yy}^{SS2}(x,\xi )q(\xi )d\xi {+\int }_{x}^{L}G_{yy}^{SS1}(x,\xi )q(\xi )d\xi$ $=-{\dfrac {QL^{4}}{EI}}{\dfrac {1}{4(\lambda L)^{4}+\pi ^{4}}}\sin \left(\pi {\dfrac {x}{L}}\right)\qquad 0<x<L$

(67)

La cual se puede entender como el cálculo del desplazamiento en un punto ${\textstyle x}$ de la viga sobre fundación flexible como la suma del desplazamiento en dicho punto debido a la carga externa ${\textstyle q(\xi )}$ aplicada en cada uno de sus puntos ${\textstyle (0<\xi {<}L)}$.

Cálculo del campo de la fuerza que el suelo ejerce sobre la viga y de fuerzas internas

A partir de (3) y (6) la fuerza que el suelo ejerce sobre la viga y las fuerzas internas en esta se calculan respectivamente como:

$f_{S}(x)=-kv(x)={\dfrac {kQL^{4}}{EI}}{\dfrac {1}{4(\lambda L)^{4}+\pi ^{4}}}\sin \left(\pi {\dfrac {x}{L}}\right)0<x<L$ (68.a) $M(x)=EI{\dfrac {d^{2}v}{dx^{2}}}(x)=QL^{2}{\dfrac {\pi ^{2}}{4(\lambda L)^{4}+\pi ^{4}}}\sin \left(\pi {\dfrac {x}{L}}\right)0<x<L$ $V(x)=-EI{\dfrac {d^{3}v}{dx^{3}}}(x)=-QL{\dfrac {\pi ^{3}}{4(\lambda L)^{4}+\pi ^{4}}}\cos \left(\pi {\dfrac {x}{L}}\right)0<x<L$

Desde el siguiente enlace se puede descargar el código de Python con la solución de este ejemplo:

https://drive.google.com/file/d/18ee0fR2uzRQFk79DYxLziZwbhIHquAgS

5.2 Viga sobre fundación flexible empotrada - libre sometida a una carga distribuida constante

Resolver la viga sobre fundación flexible presentada en la figura 12 empleando: a) la función de Green empotrada libre presentada en la sección 3.1.3 y b) la función de Green de libre-libre presentada en la sección 3.1.1.

Figura 12: Viga sobre fundación flexible en voladizo sometida a una carga externa distribuida constante.

Solución

Definición de la carga externa

La carga externa que actúa sobre la viga está definida por la siguiente función:

$q(x)=-Q\qquad 0<x<L$

(69)

Funciones de Green a emplear

Para esta alternativa la función de Green a emplear es ${\textstyle G_{yy}^{EL}(x,\xi )}$, la cual se presenta en la ecuación (22).

Cálculo del campo de desplazamiento

Para este caso el campo de desplazamiento se calcula como:

$v(x)=\int _{0}^{L}G_{yy}^{EL}(x,\xi )q(\xi )d\xi =\int _{0}^{x}G_{yy}^{EL2}(x,\xi )q(\xi )d\xi {+\int }_{x}^{L}G_{yy}^{EL1}(x,\xi )q(\xi )d\xi$ $=-Q\left[\int _{0}^{x}G_{yy}^{II}(x,\xi )d\xi {+\int }_{x}^{L}G_{yy}^{I}(x,\xi )d\xi \right]$

(70)

El cual luego de realizar las integrales respectivas, da como resultado:

$v(x)=-{\dfrac {Q}{4EI\lambda ^{4}}}{\dfrac {f(x)}{\cos(2\lambda L)+\cosh(2\lambda L)+2}}\qquad 0<x<L$

(71)

Donde:

$f(x)=\sin(\lambda x)\sinh[\lambda (2L-x)]-\sin[\lambda (2L-x)]\sinh(\lambda x)+2\cos(\lambda x)\cosh(\lambda x)+\cos(\lambda x)\cosh[\lambda (2L-x)]$ $+\cos[\lambda (2L-x)]\cosh(\lambda x)-2-\cos(2\lambda L)-\cosh(2\lambda L)$

(72)

Cálculo del campo de la fuerza que el suelo ejerce sobre la viga y de fuerzas internas

A partir de (3) y empleando (71) se tiene que la fuerza que el suelo hace sobre la viga es:

$f_{S}(x)={\dfrac {Q}{4EI\lambda ^{4}}}{\dfrac {k\cdot f(x)}{\cos(2\lambda L)+\cosh(2\lambda L)+2}}\qquad 0<x<L$

(73)

Mientras que reemplazando (71) en (6) se obtiene que los campos de momento flector y fuerza cortante son respectivamente:

$M(x)={\dfrac {g(x)}{2\lambda ^{2}[\cos(2\lambda L)+\cosh(2\lambda L)+2]}}0<x<L$ (74.a) $V(x)={\dfrac {h(x)}{\lambda [\cos(2\lambda L)+\cosh(2\lambda L)+2]}}0<x<L$ (74.b) (74.c)

Donde:

$g(x)=Q\{2\sin(\lambda L)\sin[\lambda (L-x)]\cosh(\lambda x)-2\sin(\lambda L)\cos[\lambda (L-x)]\sinh(\lambda x)+3\sin(\lambda x)\sinh(\lambda x)$ $-\sin(\lambda x)\sinh[\lambda (2L-x)]-\cos(\lambda x)\cosh(\lambda x)+\cos(\lambda x)\cosh[\lambda (2L-x)]\}$

(75)

Y:

$h(x)=Q\{-\sin(\lambda x)\cosh(\lambda x)+\sin[\lambda (2L-x)]\cosh(\lambda x)-\cos(\lambda x)\sinh(\lambda x)+\cos(\lambda x)\sinh[\lambda (2L-x)]$

(76)

Funciones de Green a emplear

Para la solución de este problema se emplearán las siguientes funciones de Green:

${\textstyle G_{yy}^{LL}(x,\xi )}$: Campo de desplazamientos de la viga sobre fundación flexible libre-libre sometida a una carga puntual unitaria presentada en la figura 5a y cuyo valor se presenta en la ecuación (16).

${\textstyle G_{y\theta }^{LL}(x,\xi )}$: Campo de desplazamientos de la viga sobre fundación flexible libre-libre sometida a un momento puntual unitario presentada en la figura 7a y cuyo valor se presenta en la ecuación (43).

${\textstyle G_{\theta y}^{LL}(x,\xi )}$: Campo de rotaciones de la viga sobre fundación flexible libre-libre sometida a una carga puntual unitaria presentada en la figura 5a.

${\textstyle G_{\theta \theta }(x,\xi )}$: Campo de rotaciones de la viga sobre fundación flexible libre-libre sometida a un momento puntual unitario presentada en la figura 7a.

En este punto es importante resaltar que debido a que ${\textstyle G_{\theta y}^{LL}(x,\xi )}$, ${\textstyle G_{y\theta }^{LL}(x,\xi )}$ y ${\textstyle G_{\theta \theta }^{LL}(x,\xi )}$ son derivadas de ${\textstyle G_{yy}^{LL}(x,\xi )}$, en realidad es posible indicar que lo es necesario emplear está última, mientras que los tramos de las primeras se definen como:

$G_{\theta y}^{LL}(x,\xi )={\dfrac {\partial G_{yy}^{LL}}{\partial x}}(x,\xi )={\begin{cases}G_{\theta y}^{LL1}(x,\xi )(x,\xi )&0<x\leq \xi \\G_{\theta y}^{LL2}(x,\xi )(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

$G_{y\theta }^{LL}(x,\xi )={\dfrac {\partial G_{yy}^{LL}}{\partial \xi }}(x,\xi )={\begin{cases}G_{y\theta }^{LL1}(x,\xi )&0<x\leq \xi \\G_{y\theta }^{LL2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

$G_{\theta \theta }^{LL}(x,\xi )={\dfrac {\partial ^{2}G_{yy}^{LL}}{\partial x\partial \xi }}(x,\xi )={\begin{cases}G_{\theta \theta }^{LL1}(x,\xi )&0<x\leq \xi \\G_{\theta \theta }^{LL2}(x,\xi )\qquad &\xi \leq x<L\end{cases}}$

Para resolver este problema se debe tener en cuenta que solo integrando la función de Green libre-libre multiplicada por la carga externa que actúa sobre la viga, se obtiene una solución que cumplirá las mismas condiciones de frontera homogéneas que la función de Green posee, es decir:

1. Segunda derivada derivada (o momento flector) igual a cero en ${\textstyle x=0}$.
2. Tercera derivada (o fuerza cortante) igual a cero en ${\textstyle x=0}$.
3. Segunda derivada derivada (o momento flector) igual a cero en ${\textstyle x=L}$.
4. Tercera derivada (o fuerza cortante) igual a cero en ${\textstyle x=L}$.

Las cuales cuales no coinciden con aquellas que posee la viga sobre fundación flexible a estudiar (ver figura 12).

Cálculo del campo de desplazamiento

Para poder hacer cumplir las condiciones de frontera del problema real, a la viga sobre fundación flexible libre-libre no solo se le debe agregar el efecto de la carga distribuida sino también aquel del momento y fuerza puntuales debidos al empotramiento en el extremo izquierdo de esta, así como el cumplimiento de las condiciones de desplazamiento y rotaciones nulos en ${\textstyle x=0}$ (ver figura 13).

Figura 13: Viga sobre fundación flexible en voladizo sometida a una carga externa distribuida constante.

Con base en lo anterior, las dos ecuaciones a emplear para obtener el valor de las reacciones de la viga sobre fundación flexible son:

$v(0)=0=FY_{1}G_{yy}^{LL2}(0,0)+M_{1}G_{y\theta }^{LL2}(0,0)-Q\int _{0}^{L}G_{yy}^{LL1}(0,\xi )d\xi$ (77.a) ${\dfrac {dv}{dx}}(0)=0=FY_{1}G_{\theta y}^{LL2}(0,0)+M_{1}G_{\theta \theta }^{LL2}(0,0)-Q\int _{0}^{L}G_{\theta y}^{LL1}(0,\xi )d\xi$ (77.b)

De cuya solución se obtiene:

$FY_{1}=-{\dfrac {Q}{\lambda }}{\dfrac {\sin(2\lambda L)+\sinh(2\lambda L)}{\cos(2\lambda L)+\cosh(2\lambda L)+2}}$ (78.a) $M_{1}={\dfrac {Q}{2\lambda ^{2}}}{\dfrac {\cos(2\lambda L)-\cosh(2\lambda L)}{\cos(2\lambda L)+\cosh(2\lambda L)+2}}$ (78.b)

Con lo cual, se tiene que el campo de desplazamiento ahora se calcula como:

$v(x)=FY_{1}G_{yy}^{LL2}(x,0)+M_{1}G_{y\theta }^{LL2}(x,0)-Q\left[\int _{0}^{x}G_{yy}^{LL2}(0,\xi )d\xi {+\int }_{x}^{L}G_{yy}^{LL1}(0,\xi )d\xi \right]$

(79)

El cual luego de realizar las integrales y usando (78), da como resultado:

$v(x)=-{\dfrac {Q}{4EI\lambda ^{4}}}{\dfrac {f(x)}{\cos(2\lambda L)+\cosh(2\lambda L)+2}}\qquad 0<x<L$

(80)

Donde:

$f(x)=\sin(\lambda x)\sinh[\lambda (2L-x)]-\sin[\lambda (2L-x)]\sinh(\lambda x)+2\cos(\lambda x)\cosh(\lambda x)+\cos(\lambda x)\cosh[\lambda (2L-x)]$ $+\cos[\lambda (2L-x)]\cosh(\lambda x)-2-\cos(2\lambda L)-\cosh(2\lambda L)$

(81)

Cálculo del campo de la fuerza que el suelo ejerce sobre la viga y de fuerzas internas

Dado que el campo de desplazamiento presentado en (81) es igual al presentado en (71) la fuerza que el suelo ejerce sobre la viga y las fuerzas internas en esta serán también iguales con ambas alternativas.

Desde el siguiente enlace se puede descargar el código de Python con la solución de este ejemplo:

https://drive.google.com/file/d/1_i-gzT1M8_PfHonABU7wzhmksbqhUeIi

5.3 Viga sobre fundación flexible libre - libre sometida a cargas puntuales y distribuidas así como a momento puntual

Resolver la viga sobre fundación flexible presentada en la figura 11 empleando la función de Green libre-libre presentada en la sección 3.1.1.

Figura 14: Viga sobre fundación flexible simplemente libre libre sometida a cargas puntuales y distribuidas así como a momento puntual.)

Funciones de Green a emplear

Para la solución de este ejercicio se emplearán las funciones de Green ${\textstyle G_{yy}^{LL}(x,\xi )}$ y ${\textstyle G_{y\theta }^{LL}(x,\xi )}$, las cuales se presentan en las ecuaciones (9) y (41) respectivamente.

Cálculo del campo de desplazamiento

El campo de desplazamiento de la viga sobre fundación flexible tiene 4 tramos, los cuales se definen como:

$v(x)={\begin{cases}v_{1}(x)&0<x\leq 1{\hbox{ m}}\\v_{2}(x)&1{\hbox{ m}}\leq x\leq 4{\hbox{ m}}\\v_{3}(x)&4{\hbox{ m}}\leq x\leq 5{\hbox{ m}}\\v_{4}(x)&5{\hbox{ m}}\leq x<10{\hbox{ m}}\end{cases}}$

(82)

Donde cada uno de estos se calcula como:

$v_{1}(x)=-250G_{yy}^{LL1}(x,1)-100G_{y\theta }^{LL1}(x,4)-200\int _{5}^{10}G_{yy}^{LL1}(x,\xi )d\xi$

(83)

$v_{1}(x)=3.399293\times 10^{-4}\sin(\lambda x)\cosh(\lambda x)+3.399293\times 10^{-4}\cos(\lambda x)\sinh(\lambda x)$ $-1.716465\times 10^{-3}\cos(\lambda x)\cosh(\lambda x)$

(84)

$v_{2}(x)=-250G_{yy}^{LL2}(x,1)-100G_{y\theta }^{LL1}(x,4)-200\int _{5}^{10}G_{yy}^{LL1}(x,\xi )d\xi$

(85)

$v_{2}(x)=1.815758\times 10^{-3}\sin(\lambda x)\sinh(\lambda x)-2.085686\times 10^{-3}\sin(\lambda x)\cosh(\lambda x)$ $+1.952791\times 10^{-3}\cos(\lambda x)\sinh(\lambda x)-1.595276\times 10^{-3}\cos(\lambda x)\cosh(\lambda x)$

(86)

$v_{3}(x)=-250G_{yy}^{LL2}(x,1)-100G_{y\theta }^{LL2}(x,4)-200\int _{5}^{10}G_{yy}^{LL1}(x,\xi )d\xi$

(87)

$v_{3}(x)=1.331976\times 10^{-3}\sin(\lambda x)\sinh(\lambda x)-1.628203\times 10^{-3}\sin(\lambda x)\cosh(\lambda x)$ $-2.305245\times 10^{-4}\cos(\lambda x)\sinh(\lambda x)+4.693530\times 10^{-4}\cos(\lambda x)\cosh(\lambda x)$

(88)

$v_{4}(x)=-250G_{yy}^{LL2}(x,1)-100G_{y\theta }^{LL2}(x,4)-200\int _{5}^{x}G_{yy}^{LL2}(x,\xi )d\xi -200\int _{x}^{10}G_{yy}^{LL1}(x,\xi )d\xi$

(89)

$v_{4}(x)=-9.182279\times 10^{-5}\sin[\lambda (x-5)]\sinh[\lambda (x-5)]+4.461867\times 10^{-5}\sin[\lambda (x-5)]\cosh[\lambda (x-5)]$ $-1.632428\times 10^{-3}\cos[\lambda (x-5)]\sinh[\lambda (x-5)]+1.708853\times 10^{-3}\cos[\lambda (x-5)]\cosh[\lambda (x-5)]$ $-3.636363\times 10^{-3}$

(90)

Como resumen de lo anterior, en la figura 15 se presenta la gráfica del campo de desplazamiento en toda la viga.

(a) Campo de desplazamiento.	(b) Campo de la fuerza que el suelo le ejerce a la viga.
Figura 15: Campos de desplazamiento y fuerza que el suelo ejerce sobre la viga.

Cálculo de la fuerza distribuida que el suelo realiza sobre la viga

A partir del campo de desplazamiento es fácil calcular la fuerza que el suelo le hace a la viga como ${\textstyle -kv(x)}$, es decir:

$f_{S}(x)={\begin{cases}f_{S}^{1}(x)&0<x\leq 1{\hbox{ m}}\\f_{S}^{2}(x)&1{\hbox{ m}}\leq x\leq 4{\hbox{ m}}\\f_{S}^{3}(x)&4{\hbox{ m}}\leq x\leq 5{\hbox{ m}}\\f_{S}^{4}(x)&5{\hbox{ m}}\leq x<10{\hbox{ m}}\end{cases}}$

(91)

Donde:

$f_{S}^{1}(x)=-18.696111\sin(\lambda x)\cosh(\lambda x)-18.696111\cos(\lambda x)\sinh(\lambda x)+94.405583\cos(\lambda x)\cosh(\lambda x)$

(92)

$f_{S}^{2}(x)=-99.866674\sin(\lambda x)\sinh(\lambda x)+114.712745\sin(\lambda x)\cosh(\lambda x)$ $-107.403483\cos(\lambda x)\sinh(\lambda x)+87.740186\cos(\lambda x)\cosh(\lambda x)$

(93)

$f_{S}^{3}(x)=-73.258680\sin(\lambda x)\sinh(\lambda x)+89.551174\sin(\lambda x)\cosh(\lambda x)$ $+12.678850\cos(\lambda x)\sinh(\lambda x)-25.814414\cos(\lambda x)\cosh(\lambda x)$

(94)

$f_{S}^{4}(x)=5.050253\sin[\lambda (x-5)]\sinh[\lambda (x-5)]-2.454027\sin[\lambda (x-5)]\cosh[\lambda (x-5)]$ $+89.783516\cos[\lambda (x-5)]\sinh[\lambda (x-5)]-93.986929\cos[\lambda (x-5)]\cosh[\lambda (x-5)]+200$

(95)

Como revisión de los anteriores resultados a continuación se realizará la revisión del equilibrio vertical y rotacional respecto al nodo 1 de toda la viga:

$\sum FY=-250-200\cdot 5+\int _{0}^{10}f_{S}(x)dx\approx 0$ (96.a) $\sum M_{1}=-250\cdot 1-100-200\cdot 7.5+\int _{0}^{10}f_{S}(x)\cdot xdx\approx 0$ (96.b)

Cálculo de los campos de fuerzas internas

A continuación se calculan las fuerzas internas en la viga a partir de su campo de desplazamiento.

Campo de momento flector

A partir del campo de desplazamiento y empleando (6.a), el campo de momento flector se calcula como:

$M(x)=EI{\dfrac {d^{2}v}{dx^{2}}}(x)={\begin{cases}M_{1}(x)&0<x\leq 1{\hbox{ m}}\\M_{2}(x)&1{\hbox{ m}}\leq x<4{\hbox{ m}}\\M_{3}(x)&4{\hbox{ m}}<x\leq 5{\hbox{ m}}\\M_{4}(x)&5{\hbox{ m}}\leq x<10{\hbox{ m}}\end{cases}}$

(97)

Donde:

$M_{1}(x)=236.013957\sin(\lambda x)\sinh(\lambda x)-46.740277\sin(\lambda x)\cosh(\lambda x)+46.740277\cos(\lambda x)\sinh(\lambda x)$

(98)

$M_{2}(x)=219.350465\sin(\lambda x)\sinh(\lambda x)-268.508707\sin(\lambda x)\cosh(\lambda x)-286.781862\cos(\lambda x)\sinh(\lambda x)$ $+249.66668430318\cos(\lambda x)\cosh(\lambda x)$

(99)

$M_{3}(x)=-64.536036\sin(\lambda x)\sinh(\lambda x)+31.697125\sin(\lambda x)\cosh(\lambda x)-223.877936\cos(\lambda x)\sinh(\lambda x)$ $+183.146701\cos(\lambda x)\cosh(\lambda x)$

(100)

$M_{4}(x)=-234.967314\sin =[\lambda (x-5)]\sinh[\lambda (x-5)]+224.458790\sin[\lambda (x-5)]\cosh[\lambda (x-5)]$ $+6.135067\cos[\lambda (x-5)]\sinh[\lambda (x-5)]-12.625633\cos[\lambda (x-5)]\cosh[\lambda (x-5)]$

(101)

De forma gráfica, el campo de momento flector se presenta en la figura 16FigMomentoFlector.

Campo de fuerza cortante

Mientras que el campo de fuerza cortante se calcula a partir del campo de desplazamiento como (ver (6.b)):

$V(x)=-EI{\dfrac {d^{3}v}{dx^{3}}}(x)={\begin{cases}V_{1}(x)&0<x<1{\hbox{ m}}\\V_{2}(x)&1{\hbox{ m}}<x\leq 4{\hbox{ m}}\\V_{3}(x)&4{\hbox{ m}}\leq x\leq 5{\hbox{ m}}\\V_{4}(x)&5{\hbox{ m}}\leq x<10{\hbox{ m}}\end{cases}}$

(102)

Donde:

$V_{1}(x)=41.805775\sin(\lambda x)\sinh(\lambda x)-105.548650\sin(\lambda x)\cosh(\lambda x)-105.548650\cos(\lambda x)\sinh(\lambda x)$

(103)

$V_{2}(x)=-8.172003\sin(\lambda x)\sinh(\lambda x)+13.557825\sin(\lambda x)\cosh(\lambda x)-209.750846\cos(\lambda x)\sinh(\lambda x)$ $+248.333492\cos(\lambda x)\cosh(\lambda x)$

(104)

$V_{3}(x)=-114.296642\sin(\lambda x)\sinh(\lambda x)+110.767087\sin(\lambda x)\cosh(\lambda x)-53.044302\cos(\lambda x)\sinh(\lambda x)$ $+85.945872\cos(\lambda x)\cosh(\lambda x)$

(105)

$V_{4}(x)=-97.637337\sin[\lambda (x-5)]\sinh[\lambda (x-5)]+99.434222\sin[\lambda (x-5)]\cosh[\lambda (x-5)]$ $+110.72693z\cos[\lambda (x-5)]\sinh[\lambda (x-5)]-103.124708\cos[\lambda (x-5)]\cosh[\lambda (x-5)]$

(106)

Como resumen de los anteriores resultados, en la figura 16 se presenta de forma gráfica el campo de fuerza cortante en toda la viga.

(a) Momento flector.	(b) Fuerza cortante.
Figura 16: Campos de fuerzas internas.

Desde el siguiente enlace se puede descargar el código de Python con la solución de este ejemplo:

https://drive.google.com/file/d/1OVWDyO10_n7tHrZyyISo-OOYUmVrB8bn

6 Conclusiones

Las principales conclusiones de este artículo son:

1. Se presenta una metodología basada 100% en las funciones de Green para el análisis de vigas sobre fundación flexible, con la cual se pueden calcular de forma exacta los campos de desplazamiento, momento flector, fuerza cortante y fuerza que el suelo ejerce sobre la viga.
2. Una vez implementado de forma numérica, el procedimiento presentado permite analizar de forma muy eficiente vigas sobre fundación flexible.
3. Pese a que las funciones de Green más eficientes para resolver un determinado problema son aquellas que comparten las mismas condiciones de frontera (homogéneas) del problema en estudio, es posible su empleo para resolver problemas con diferentes condiciones de frontera. Esto se demuestra en el ejemplo [[#5.2.2 Solución empleando la función de Green de la viga sobre fundación flexible libre libre (figura 5a)|5.2.2]] y le da mucha generalidad al procedimiento.
4. Se presentan tres ejemplos con la solución de igual número de vigas sobre fundación flexible empleando la metodología presentada en este artículo, las cuales son la respuesta exacta de cada problema pues cumplen sus ecuaciones diferenciales gobernantes (equilibrio de cada punto), condiciones de frontera y equilibrio de toda la viga.

BIBLIOGRAFÍA

[1] Reddy, JN. (2004) "An introduction to the finite element method", Volume 1221. McGraw-Hill New York

[2] Molina-Villegas, Juan Camilo and Giraldo, Harold Nolberto Diaz and Ochoa, Andrés Felipe Acosta. (2020) "Analytical formulation of the stiffness method for 2D reticular structures using Green functions", Volume 36. Revista Internacional de Metodos Numericos para Calculo y Diseno en Ingenieria 3

[3] Molina-Villegas, Juan Camilo. (2021) "Análisis estructural: Métodos clásicos y matriciales.". Ecoe Ediciones

[4] Stakgold, Ivar and Holst, Michael J. (2011) "Green's functions and boundary value problems", Volume 99. John Wiley & Sons

[5] Duffy, Dean G. (2015) "Green's functions with applications". CRC Press

[6] Rother, Tom. (2017) "Green's Functions in Classical Physics", Volume 938. Springer

[7] Podio-Guidugli, Paolo and Favata, Antonino. (2014) "Elasticity for Geotechnicians: A Modern Exposition of Kelvin, Boussinesq, Flamant, Cerruti, Melan, and Mindlin Problems", Volume 204. Springer

[8] Kausel, Eduardo. (2006) "Fundamental solutions in elastodynamics: a compendium". Cambridge University Press

[9] Cole, Kevin and Beck, James and Haji-Sheikh, A and Litkouhi, Bahman. (2010) "Heat Conduction Using Greens Functions". CRC Press

[10] Mandelis, Andreas. (2013) "Diffusion-wave fields: mathematical methods and Green functions". Springer Science & Business Media

[11] Beer, Gernot. (2000) "Programming the boundary element method". John Wiley & Sons, Inc.

[12] Brebbia, Carlos Alberto and Dominguez, Jose. (1994) "Boundary elements: an introductory course". WIT press

[13] Katsikadelis, John T. (2002) "Boundary elements: theory and applications". Elsevier

[14] Rodriguez-Castellanos, Alejandro and Flores, E and Sánchez-Sesma, Francisco José and Ortiz-Aleman, Carlos and Nava-Flores, Mauricio and Martin, Roland. (2011) "Indirect Boundary Element Method applied to fluid–solid interfaces", Volume 31. Elsevier. Soil Dynamics and Earthquake Engineering 3 470–477

[15] Sánchez-Sesma, FJ and Ramos-Martinez, J and Campillo, M. (1993) "An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves", Volume 22. Wiley Online Library. Earthquake engineering & structural dynamics 4 279–295

[16] Hetényi, Miklós. (1971) "Beams on elastic foundation: theory with applications in the fields of civil and mechanical engineering". University of Michigan