| Line 3,082: | Line 3,082: | ||
<div id="cite-1"></div> | <div id="cite-1"></div> | ||
| − | '''[[#citeF-1|[1]]]''' Challis, Lawrie and Sheard, Fred. The green of Green functions, | + | '''[[#citeF-1|[1]]]''' Challis, Lawrie and Sheard, Fred. The green of Green functions, Physics Today, 56(12):41-46,2003. |
<div id="cite-2"></div> | <div id="cite-2"></div> | ||
| Line 3,091: | Line 3,091: | ||
<div id="cite-4"></div> | <div id="cite-4"></div> | ||
| − | '''[[#citeF-4|[4]]]''' Sánchez-Sesma, F. J. and Ramos-Martínez, J. and Campillo, M. An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves | + | '''[[#citeF-4|[4]]]''' Sánchez-Sesma, F. J. and Ramos-Martínez, J. and Campillo, M. An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves. Earthquake engineering & structural dynamics, 22(4):279-295, 1993. |
<div id="cite-5"></div> | <div id="cite-5"></div> | ||
| − | '''[[#citeF-5|[5]]]''' Fairweather, G. and Karageorghis, A. and Martin, P.A. The method of fundamental solutions for scattering and radiation problems, | + | '''[[#citeF-5|[5]]]''' Fairweather, G. and Karageorghis, A. and Martin, P.A. The method of fundamental solutions for scattering and radiation problems, Engineering Analysis with Boundary Elements, 27(7) :759-769, 2003. |
<div id="cite-6"></div> | <div id="cite-6"></div> | ||
| − | '''[[#citeF-6|[6]]]''' Thomson, William T. Transmission of elastic waves through a stratified solid medium | + | '''[[#citeF-6|[6]]]''' Thomson, William T. Transmission of elastic waves through a stratified solid medium. Journal of applied Physics, 21(2):89-93, 1950. |
<div id="cite-7"></div> | <div id="cite-7"></div> | ||
| − | '''[[#citeF-7|[7]]]''' Boussinesq, Joseph. Application des potentiels a l'étude de l'équilibre et du mouvement des solides élastiques. Gauthier-Villars | + | '''[[#citeF-7|[7]]]''' Boussinesq, Joseph. Application des potentiels a l'étude de l'équilibre et du mouvement des solides élastiques. Gauthier-Villars, 1885. |
<div id="cite-8"></div> | <div id="cite-8"></div> | ||
| Line 3,106: | Line 3,106: | ||
<div id="cite-9"></div> | <div id="cite-9"></div> | ||
| − | '''[[#citeF-9|[9]]]''' Mindlin, R.D. Force at a Point in the Interior of a Semi-Infinite Solid, | + | '''[[#citeF-9|[9]]]''' Mindlin, R.D. Force at a Point in the Interior of a Semi-Infinite Solid. Physics, 7(5):195-202, 1936. |
<div id="cite-10"></div> | <div id="cite-10"></div> | ||
| − | '''[[#citeF-10|[10]]]''' George Gebriel Stokes. On dynamical theory of diffraction | + | '''[[#citeF-10|[10]]]''' George Gebriel Stokes. On dynamical theory of diffraction. Transactions of the Cambridge Philosophical Society, 9:1-48, 1849. |
<div id="cite-11"></div> | <div id="cite-11"></div> | ||
| − | '''[[#citeF-11|[11]]]''' H. Lamb. On the propagation of tremors over the surface of an elastic solid | + | '''[[#citeF-11|[11]]]''' H. Lamb. On the propagation of tremors over the surface of an elastic solid. Philosophical Transactions of the Royal Society of London, 203:1-42, 1904. |
<div id="cite-12"></div> | <div id="cite-12"></div> | ||
| − | '''[[#citeF-12|[12]]]''' Chao, C.C. Dynamical response of an elastic half-space to tangential surface loadings | + | '''[[#citeF-12|[12]]]''' Chao, C.C. Dynamical response of an elastic half-space to tangential surface loadings. Journal of Applied Mechanics, 27:559, 1960. |
<div id="cite-13"></div> | <div id="cite-13"></div> | ||
| Line 3,142: | Line 3,142: | ||
<div id="cite-21"></div> | <div id="cite-21"></div> | ||
| − | '''[[#citeF-21|[21]]]''' M. Eisemberger and D.Z. Yankelevsky. Exact stiffness matrix for beams on elastic foundation | + | '''[[#citeF-21|[21]]]''' M. Eisemberger and D.Z. Yankelevsky. Exact stiffness matrix for beams on elastic foundation. Computer & Structures, 21(6):1355-1359, 1985. |
Juan Camilo Molina-Villegasa,b, Harold Nolberto Diaz Giraldob, Andrés Felipe Acosta Ochoab
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.
Las funciones de Green (F.G.) se definen como la respuesta de un medio ante una carga puntual unitaria y son usadas ampliamente para la solución de problemas de valores en la frontera. Desafortunadamente, en el análisis estructural su uso es limitado y solo se emplean de forma indirecta y con otro nombre en el cálculo de lineas de influencia y en la formulación del método del trabajo virtual. En este artículo se presenta una metodología para obtener la respuesta analítica o exacta de estructuras reticulares planas, la cual mezcla al método de rigidez y a las funciones de Green, estas últimas empleadas para el cálculo de los campos de desplazamiento. En particular se realizará la formulación para elementos tipo barra (sometidos a fuerza axial), viga (sometidos a fuerza cortante y momento flector), viga sobre fundación flexible (sometidos a fuerza cortante y momento flector) y pórtico plano (sometidos a fuerza axial, fuerza cortante y momento flector). Esta formulación tiene como propiedad principal que puede ser empleada para cualquier distribución de carga externa y minimiza el número de elementos a emplear en las discretizaciones. Además se presenta la equivalencia de esta formulación con aquella obtenida mediante una implementación “exacta” del método de elementos finitos.
Green functions (F.G.) are defined as the response of a medium to a unit point load and are widely used to solve boundary value problems. Unfortunately, in structural analysis its use is limited and they are only used indirectly and with another name in the calculation of influence lines and in the formulation of the virtual work method. This article presents a methodology to obtain the analytical or exact response of two dimensional frames, which mixes the stiffnes method and the Green functions, the latter used for the calculation of displacement fields. In particular, the formulation will be carried out for bar elements (subjected to axial force), beam elements (subjected to shear force and bending moment), beam over flexible foundation elements (subjected to shear force and bending moment) and two dimensional frames (subjected to axial force, cutting force and bending moment). This formulation has as its main property that it can be used for any external load distribution and minimizes the number of elements to be used in discretizations. In addition, the equivalence of this formulation with that obtained by an “exact” implementation of the finite element method is presented.
keywords Funciones de Green, Método de rigidez, Método de elementos finitos, Pórticos planos.
Las funciones de Green (F.G.) son la respuesta de un medio ante la acción de una fuente o fuerza puntual unitaria y son muy importantes en la solución de problemas de valor en la frontera (P.V.F.) de fenómenos físicos e ingenieriles, cuyas ecuaciones diferenciales (E.D.) gobernantes son lineales [1]. Esto es debido a que ellas por si solas son solución de problemas fundamentales, pueden ser empleadas para resolver problemas con fuentes distribuidas [2] o de forma indirecta se emplean en métodos numéricos de frontera o contorno, como el método directo de elementos de frontera [3], el método indirecto de elementos de frontera [4] o el método de las soluciones fundamentales [5].
En geotecnia las principales F.G. son: la respuesta de un espacio completo o infinito [6], la respuesta de un semi-espacio ante la acción de un fuerza en superficie [7], [8] y en el interior [9]. En sismología las principales F.G. son: la respuesta de un espacio completo [10],la respuesta de un semi-espacio debida a la acción de una fuerza normal [11] y tangencial [12], para el lector interesado en este tema se le recomiendan los libros [13] y [14]. Pese a lo anteriormente mencionado, es desafortunado el poco uso que se les da a las F.G. en el análisis estructural, donde su uso (con otro nombre) se limita al cálculo de lineas de influencia y a la formulación del método del trabajo virtual.
Para el caso particular de la elasticidad y análisis estructural las funciones de Green se expresan como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle G_{ij}(\mathbf{x},\boldsymbol \xi )} , la cual corresponde al desplazamiento en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i}
del punto Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{x}}
, debido a una fuerza puntual unitaria en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle j}
aplicada en el punto Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol \xi }
.
En la actualidad los métodos matriciales son los más empleados para la solución de problemas de análisis estructural, entre estos destacan el método de rigidez y el método de elementos finitos (M.E.F.). El primero de estos es mas complejo de implementar y busca la solución exacta del modelo estructural en estudio, mientras que el segundo es mas general y simple de implementar y busca una solución aproximada. Por esta razón es usual emplear el método de rigidez para solucionar problemas “sencillos” que involucran modelos lineales con elementos de sección transversal constante y cargas simples, como puntuales o distribuidas con variación lineal. Mientras que el M.E.F. suele emplearse para estructuras de sección transversal constante o variable pero sometidas a cargas cuya definición es compleja, o para problemas no lineales. Pese a lo anterior, ambos métodos tienen pequeñas falencias intrínsecas, en el método de rigidez hay gran dificultad para manejo de cargas genéricas y no es usual el cálculo del campo de desplazamientos ([15], [16] y [17]) mientras que en el M.E.F. la solución es aproximada, lo cual lleva contradicciones a la hora de calcular las fuerzas internas a partir de los campos de desplazamientos ([18] y [19]).
En este artículo se presenta una formulación particular del método de rigidez para obtener la solución analítica de problemas de análisis estructural formados por estructuras reticulares planas (las cuales pueden estar formadas por elementos tipo barra, viga, viga sobre fundación flexible o pórtico), la cual emplea a las funciones de Green para el cálculo de los campos de desplazamiento. Esta metodología busca combinar las principales fortalezas del método de rigidez (obtención de soluciones analíticas) y del M.E.F. (posibilidad de manejar cargas complejas), para la solución total de estas estructuras, es decir, para el calculo de sus reacciones, campos de desplazamiento y campos de fuerzas internas. Al final de este artículo se presentan tres apéndices donde se presenta la equivalencia de la actual metodología (la cual parte de las ecuaciones diferenciales gobernantes de cada tipo de elemento) y una “exacta” por el M.E.F. (la cual parte de la forma débil de las ecuaciones diferenciales en lugar de estas).
A continuación se comenzarán presentando las ecuaciones diferenciales gobernantes de cada uno de los cuatro tipos de elementos a analizar (barra, viga, pórtico y viga sobre fundación flexible), luego se describirá en detalle la metodología de análisis propuesta en este artículo, seguido a esto se realizará la formulación “exacta” del método de rigidez para cada uno de los tipos de elemento estudiados y se realizaran tres ejemplos ilustrativos.
Se define como elemento tipo barra, a aquel cuya única fuerza interna es axial y se encuentra sometido a una fuerza externa distribuida por unidad de longitud (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p(x)} ), la cual actúa en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
y su valor es positivo si se dirige en la dirección de dicho eje local (ver Figura 2). Para el caso que el elemento sea de sección transversal constante con área de la sección transversal Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
, material elástico lineal con módulo de elasticidad Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E} , su ecuación diferencial gobernante es:
|
(1) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u(x)}
es el desplazamiento en dirección axial del elemento (eje local Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
).
| Figura 2: Elemento tipo barra y el sistema coordenado local empleado para la definición de su E.D. gobernante. |
Mientras que a partir de la ley de Hooke unidimensional su fuerza interna axial se calcula como:
|
(2) |
El elemento tipo viga está sometido tanto a fuerza cortante en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
como a momento flector alrededor del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle z}
(perpendicular a los ejes Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
, ver Figura 3). Sus cargas externas se definen en términos de la carga por unidad de longitud en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
(Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(x)}
), la cual es positiva en la dirección de este eje (ver Figura 3). Si el elemento es de sección transversal constante, con momento de inercia Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I}
y material elástico lineal con módulo de elasticidad Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E}
, su ecuación diferencial gobernante es:
|
(3) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v(x)}
es el desplazamiento en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
, es decir, perpendicular al eje longitudinal del elemento.
| Figura 3: Elemento tipo viga y el sistema coordenado local empleado para la definición de su E.D. gobernante. |
Mientras que las fuerzas internas cortante y momento flector se calculan a partir del campo de desplazamiento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v(x)} , respectivamente como:
|
El elemento tipo pórtico plano se define como la superposición de un elemento tipo barra y uno tipo viga, por lo cual sus fuerzas internas son la fuerza axial, la fuerza cortante y el momento flector respecto al eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle z} . Sus ecuaciones diferenciales gobernantes son 1 para el desplazamiento axial y 3 para el desplazamiento perpendicular al eje longitudinal del elemento (eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y} ) y sus fuerzas internas se calculan a partir de los campos de desplazamiento empleando 2 y 2.2.
El elemento viga sobre fundación flexible se define como un elemento tipo viga que esta apoyado sobre un medio elástico, el cual tiene una rigidez por unidad de longitud Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
y genera una fuerza por unidad de longitud sobre la viga igual al Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -kv(x)}
(ver Figura 4). Sus fuerzas internas son iguales a las de la viga, es decir, fuerza cortante y momento flector, y su ecuación diferencial gobernante es ([20]):
|
(5) |
| Figura 4: Elemento tipo viga sobre fundación flexible y el sistema coordenado local empleado para la definición de su E.D. gobernante. |
Al igual que para los elementos tipo viga, las fuerzas internas se calculan a partir del campo de desplazamientos empleando 2.2.
En la formulación del método de rigidez para elementos prismáticos se expresan las fuerzas y momentos en los extremos de estos en función de los desplazamiento y rotaciones en esos mismos punto (lo cual comúnmente se conoce como grados de libertad). Esta formulación se expresa matricialmente en un sistema local de coordenadas como:
|
(6) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \{ {F'}_E \} }
es el vector de fuerzas en los extremos del elemento en coordenadas locales, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {[K'}_E]}
es la matriz de rigidez en coordenadas locales, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \{ \Delta' _E \} }
es el vector de desplazamientos en los extremos del elemento en coordenadas locales y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \{ \left.{F'}_E \right.^f \} }
es el vector de fuerzas de empotramiento en coordenadas locales.
De analizar los dos términos del lado derecho de 6 es evidente que el problema se puede descomponer como la superposición o suma de dos problemas fundamentales. El primero se llamará problema homogéneo (la justificación de esto se presentará mas adelante) y es el encargado de la aparición del término que incluye a la matriz de rigidez, este problema solo depende de los desplazamientos de los nodos (nodales) y no de las cargas externas. El segundo problema se llamará problema empotrado (la justificación para su nombre también se dará mas adelante) y es el encargado de la aparición del vector de fuerzas de empotramiento y no depende de los desplazamientos nodales, pero si de las fuerzas externas. A partir de esta idea, las componentes de los campos de desplazamientos se expresan como:
|
(7) |
Donde:
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d(x)}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d_h(x)}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d_f(x)}
A continuación, a partir de 7, se definirán y solucionarán los problemas de valor en la frontera (P.V.F.) empleados para la formulación del método de rigidez de cada uno de los cuatro tipos de elementos a analizar, así como los P.V.F. que dan lugar a los campos homogéneos y empotrados de estos mismos.
El P.V.F. que gobierna la formulación del método de rigidez para elementos tipo barra de sección transversal constante y longitud Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
es (ver la Figura 2):
|
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u(0)=u_i}
es el desplazamiento en dirección axial (eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
) en el extremo inicial del elemento y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u(L)=u_j}
es el desplazamiento en la misma dirección pero en el extremo final del elemento.
Para este caso la particularización de 7 es:
|
(9) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_h(x)}
es el campo de desplazamiento homogéneo (sección 4.1), mientras que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_f(x)} es el campo empotrado (sección 4.2).
El P.V.F. que gobierna al campo homogéneo se presenta en 4.1 y su nombre se debe a que 10.a es una ecuación diferencial homogénea.
|
Es fácil probar que la solución de 4.1 es:
|
(11) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
son conocidas como las funciones de forma “exactas” del problema y en este caso son iguales a aquellas que suelen emplearse en el M.E.F. para interpolar el campo de desplazamientos. En este punto es importante resaltar que en el M.E.F. es usual emplear 11 para aproximar el campo total de desplazamientos, lo cual, a partir de 9, se observa que es inexacto pues no incluye al campo empotrado.
Por su parte, empleando la ley de Hooke uniaxial, se tiene que la fuerza axial homogénea se puede calcular como:
|
(12) |
Mientras que las fuerzas en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
en los extremos de la barra son:
|
El problema de P.V.F. que gobierna al campo empotrado se presenta en 4.2. Su nombre se debe a que representa la respuesta de una barra doblemente empotrada sometida a una carga externa por unidad de longitud genérica definida por la función Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p(x)}
(ver Figura 5a).
|
Para la solución del P.V.F. 4.2 se empleará la función de Green de un elemento tipo barra de sección transversal constante y material elástico lineal (ver Figura 5b), la cual corresponde a la respuesta (campo de desplazamiento axial) debido a la aplicación de una carga puntual unitaria ubicada a una distancia Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi }
del origen de coordenadas y está gobernada por el siguiente P.V.F (notar que las dos últimas ecuaciones corresponden a la continuidad de desplazamiento y equilibrio en el punto de aplicación de la fuerza puntual unitaria Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (x=\xi )}
):
|
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta (x)}
es la función delta de Dirac.
La solución de 4.2 es:
|
(16) |
La cual por conveniencia se escribirá de forma compacta como:
|
(17) |
En este punto es importante resaltar la propiedad de simetría que tiene 16 y en general todas las funciones de Green, la cual se expresa como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle G_{xx}(\xi ,x)=G_{xx}(x,\xi )} .
Multiplicando a ambos lados de 15.a por Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p(\xi )}
e integrando respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi }
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
se obtiene:
|
(18) |
Donde se empleó la siguiente propiedad de la función delta Dirac: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p(x)=\int _0^L \delta (x-\xi ) p(\xi ) d\xi } .
De comparar 18 y 14.a se concluye que:
|
(19) |
En donde es importante resaltar que 19 permite calcular el campo de desplazamiento empotrado ante cualquier carga externa.
Empleando la ley de Hooke uniaxial, a partir de 19 el campo de fuerza axial empotrado se calcula como:
|
(20) |
Y las reacciones o fuerzas de empotramiento son:
|
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
son las funciones de forma para el elemento tipo barra que se definieron en la sección 4.1. Es de destacar que la forma que expresar las fuerzas de empotramiento en función de la funciones de forma es coherente con aquella presentada para la formulación “exacta” del M.E.F. presentada en el Apéndice 11.
Reemplazando 11 y 19 en 9 se tiene que el campo de desplazamiento en el elemento, expresado en función de los desplazamientos nodales y las cargas externas es:
|
(22) |
Mientras que a partir de sumar 4.1 y 4.2, se tiene que las fuerzas en los extremos del elemento son:
|
(23) |
La cual es la formulación del método de rigidez para un elemento tipo barra y es equivalente a la presentada en el Apéndice 11 para el caso de la formulación “exacta” por el método de elementos finitos.
Para finalizar es importante resaltar que las definiciones de los P.V.F. 4.1 y 4.2 no han sido arbitrarias y son totalmente coherentes con 4 y 9. Por ejemplo si se reemplaza 9 en 4 y se emplean 4.2, se obtiene directamente el P.V.F. 4.1, lo cual también ocurre para los demás tipos de elementos analizados en este artículo.
El P.V.F. que gobierna la formulación del método de rigidez de un elemento tipo viga de sección transversal constante, material elástico lineal y longitud Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
es (ver Figura 3):
|
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v(0)=v_i}
es el desplazamiento en dirección del eje local Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
en extremo inicial del elemento, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \dfrac{dv}{dx}(0)=\theta _i}
es la pendiente en ese mismo punto, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v(L)=v_j}
es el desplazamiento en dirección del eje local Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
en extremo final del elemento y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \dfrac{dv}{dx}(L)=\theta _j}
es la rotación en ese último punto.
Para este caso la particularización de 7 es:
|
(25) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_h(x)}
es el campo de desplazamiento homogéneo (sección 5.1), mientras que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x)} es el campo empotrado (sección 5.2).
Para este caso el P.V.F. gobernante es:
|
Cuya solución, es el siguiente polinomio de orden tres:
|
(27) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _6(x)}
se conocen coma las funciones de forma “exactas” de este problema y tienen el siguiente valor:
|
Al igual que para el caso de las barras, las funciones de forma “exactas” presentadas en 5.1 son iguales a las empleadas en el M.E.F. para aproximar o interpolar el campo total de desplazamiento, lo cual de nuevo es inexacto pues se omite el campo empotrado.
A partir de 27, las fuerzas internas homogéneas son:
|
Mientras que las reacciones homogéneas son:
|
El problema de P.V.F. que gobierna al campo empotrado se presenta en 5.2 y corresponde a la respuesta de una viga doblemente empotrada sometida a una carga externa genérica definida por la función Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(x)}
(ver Figura 6a).
|
Para solucionar el P.V.F. 5.2 se definirá la función que Green asociada con este problema, es decir, la respuesta de una viga doblemente empotrada de sección transversal constante y material elástico lineal, sometida a una fuerza puntual unitaria ubicada a una distancia Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi }
del extremo izquierdo del elemento (ver Figura 6b). Esta función de Green está gobernada por el siguiente P.V.F. que contiene el cumplimiento de la ecuación diferencial gobernante, condiciones de frontera, condiciones de continuidad de desplazamientos y rotaciones en (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=\xi }
) y equilibrio vertical y rotacional en este mismo punto:
|
La solución de 5.2, es decir, la función de Green para la viga doblemente empotrada presentada en la Figura 6b es:
|
(33) |
La cual, por comodidad se expresa de forma compacta como:
|
(34) |
Multiplicando a ambos lados de 32.a por Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(\xi )}
e integrando respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi }
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
se obtiene:
|
(35) |
Donde se empleó que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(x)=\int _0^L \delta (x-\xi ) q(\xi ) d\xi } .
De comparar 18 y 31.a se concluye que:
|
(36) |
Mientras que los campos de fuerzas internas se calculan como:
|
(37) |
|
(38) |
Y las reacciones o fuerzas de empotramiento son:
|
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _6(x)}
son las funciones de forma para un elemento tipo viga y se definieron en la ecuación 5.1.
Reemplazando 27 y 36 en 25 se obtiene que el campo de desplazamientos puede expresarse en función de los desplazamientos nodales y las fuerzas externas como:
|
(40) |
Mientras que las fuerzas en los extremos del elemento se calculan a partir de la suma de 5.1 y 5.2, dando como resultado:
|
(41) |
La cual es la formulación del método de rigidez para un elemento tipo viga y es equivalente a la formulación “exacta” del M.E.F. presentada en el Apéndice 12.
Como se mencionó en la sección 2.3 el elemento tipo pórtico plano es simplemente la superposición de un elemento tipo barra y uno tipo viga, con lo cual a partir de lo presentado en 23 y 41 se tiene que la formulación del método de rigidez para un elemento tipo pórtico plano en coordenadas locales es:
|
(42) |
Mientras que los campos de desplazamiento en dirección de los ejes locales Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y}
, se calculan a partir de las ecuaciones 22 y 40 respectivamente.
El P.V.F. que gobierna la formulación del método de rigidez de un elemento tipo viga sobre fundación flexible de sección transversal constante, material elástico lineal, rigidez del suelo por unidad de longitud contante e igual Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
y longitud Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
es (ver Figura 4):
|
Donde se ha empleado como ecuación diferencial una versión equivalente de 5, en la cual se define Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda =\sqrt=4}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_i}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \theta _i} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_j}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \theta _j}
tienen el mismo significado que para la viga (sección 5).
Para este caso la particularización de 7 es:
|
(44) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_h(x)}
es el campo de desplazamiento homogéneo (sección 7.1), mientras que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x)} es el campo empotrado (sección 7.2).
Para este caso el P.V.F. gobernante es:
|
Cuya solución se expresa como:
|
(46) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _6(x)}
son las funciones de forma “exactas” de este problema y tienen el siguiente valor:
|
Y para simplificar la escritura se ha definido: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hbox{s}=\sin (\lambda L)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hbox{c}=\cos (\lambda L)}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hbox{sh}=\sinh (\lambda L)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hbox{ch}=\cosh (\lambda L)}
.
A diferencia de lo que ocurre para el caso la barra y la viga, las funciones de forma “exactas” presentadas en 7.1 no son aquellas que suelen emplearse en la formulación del M.E.F. para vigas sobre fundación flexible, en el cual se emplean las mismas funciones de forma que en este artículo se definieron para la viga. Como se presentará en el apéndice 13, si en la formulación del M.E.F. se emplean 7.1, se obtiene una formulación matricial equivalente a la presentada en el presente método.
Para este caso el P.V.F. que define al problema empotrado es (ver Figura 7a):
|
Mientras que la función de Green (7b) se define por medio del siguiente P.V.F (ver Figura 7b):
|
Donde se ha usado como convención emplear una linea horizontal encima de la letra Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle G}
para diferenciar esta función de Green de aquella de las vigas.
La solución de 7.2 se expresa en forma compacta como:
|
(50) |
Donde:
|
Y
|
|
(53) |
|
(54) |
|
(55) |
|
(56) |
|
(57) |
|
(58) |
|
(59) |
|
(60) |
Multiplicando a ambos lados de 49.a por Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(\xi )}
e integrando respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi }
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
se obtiene:
|
(61) |
Donde se empleó que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q(x)=\int _0^L \delta (x-\xi ) q(\xi ) d\xi } .
De comparar 18 y 48.a se concluye que:
|
(62) |
Mientras que las fuerzas internas empotradas se calculan como:
|
Y las reacciones o fuerzas de empotramiento se calculan como:
|
Reemplazando 46 y 62 en 44, el campo de desplazamiento total se expresa en función de los desplazamientos de los extremos del elemento y de la carga externa como:
|
(65) |
|
(66) |
Donde los términos de la matriz de rigidez son ([21]):
|
Es importante resaltar que 66 corresponde a la formulación analítica del método de rigidez para una viga sobre fundación flexible y es equivalente a aquella formulación “exacta” del M.F.F. que se presenta en el Apéndice 13.
| (a) Problema a resolver. | (b) Discretización. |
| Figura 8: Viga con articulación interior y discretización empleada para su solución. | |
Discretización
La discretización a emplear en la solución de este problema se presenta en la Figura 8b y consta de solo dos elementos. El elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
tiene una carga externa definida en tres tramos, mientras que el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
tiene una carga externa lineal sobre toda su longitud.
Definición de las funciones de carga externa
Para el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
su carga externa se define como:
|
(68) |
O en forma compacta como:
|
(69) |
Por su parte la carga externa en el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
es:
|
(70) |
Cálculo de las fuerzas de empotramiento
Según lo presentado en 5.2, las fuerzas de empotramiento del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
se calculan como:
|
Mientras que las del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
como:
|
Sistema matricial de ecuaciones de cada elemento
A partir de lo presentado en 41, la formulación del método de rigidez para el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es:
|
(73) |
Mientras que para elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
este es:
|
(74) |
Cálculo de los desplazamientos nodales desconocidos
A partir de las fuerzas externas e internas conocidas en los nodos, y empleando continuidad de desplazamientos y equilibrio en estos, se obtiene el siguiente sistema lineal de tres ecuaciones con igual número de incógnitas:
|
(75) |
Cuya solución es:
|
(76) |
Cálculo del campo de desplazamiento del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): A
Reemplazando los desplazamientos nodales del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
en 27, se obtiene que el campo de desplazamiento homogéneo para este elemento es:
|
(77) |
Mientras que debido a que la carga externa en el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
tiene tres tramos, el campo de desplazamientos empotrado también los tendrá, es decir:
|
(78) |
Donde, a partir de 36, el campo de desplazamiento empotrado en cada uno de los tres tramos se calcula como:
|
(79) |
|
(80) |
|
(81) |
Ahora sumando 77 y 78 se obtiene que el campo de desplazamiento total en el elemento es:
|
(82) |
Donde:
|
(83) |
|
(84) |
|
(85) |
Cálculo del campo de desplazamiento del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B
Procediendo de igual manera a como se hizo para el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A} , se calculan los campos de desplazamiento homogéneo, empotrado y total del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
|
|
(87) |
| Figura 9: Diagrama del campo de desplazamientos. |
Cálculos de los campos de fuerzas internas
A partir de 82 y empleando 2.2, se obtiene que las fuerzas internas para el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
son:
|
(88) |
|
(89) |
De forma similar, a partir de 87, se tiene que los campos de fuerzas internas en el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
son:
|
| (a) Fuerza cortante. | (b) Momento flector. |
| Figura 10: Diagramas de fuerzas internas. | |
Cálculo de las reacciones
Aunque es posible calcular las reacciones mediante el cálculo de un sistema lineal de ecuaciones, donde las incógnitas sean las cuatro reacciones, a continuación se calcularán a partir de los campo de fuerzas internas calculados anteriormente:
|
Para revisar la validez de los resultados anteriores, a continuación se realizará el equilibrio vertical y rotacional respecto al nudo 1 de toda la viga:
|
(92) |
|
(93) |
En el siguiente link se encuentra el código de Python empleado para la solución de este ejemplo: https://drive.google.com/open?id=1-rqbQquE5YmvSKxRtQW2z7Lpm1O8C8-Y
| (a) Problema a resolver. | (b) Discretización. |
| Figura 11: Viga sobre fundación flexible y discretización empleada para su solución. | |
Discretización
Para este problema se empleará la discretización de un solo elemento presentada en la Figura 11b, la cual muestra una de las principales bondades de la presente metodología que consiste en minimizar el número de elementos a emplear, pues usualmente en la solución de este problema se emplearía una discretización con tres elementos para obtener la respuesta exacta.
Definición de la carga externa
La carga externa en este problema es:
|
(94) |
Formulación del sistema de ecuaciones de cada elemento
A partir de 66 se tiene que la formulación del método de rigidez para el único elemento de la discretización es:
|
(95) |
Obtención de los desplazamientos nodales
La solución de 95 es:
|
(96) |
Cálculo del campo de desplazamientos
Para el cálculo del campo de desplazamientos se realizará la siguiente superposición o descomposición:
|
(97) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_h(x)}
es el campo de desplazamientos homogéneo y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x)}
es el campo de desplazamientos empotrado. A partir de 46 y empleando los desplazamientos y rotaciones nodales, se obtiene que el campo de desplazamiento homogéneo es:
|
(98) |
Mientras que el campo de desplazamiento empotrado se calcula a partir de 62 y se expresa como:
|
(99) |
Donde:
|
(100) |
|
(101) |
|
(102) |
Con lo cual el campo de desplazamiento total es:
|
(103) |
|
(104) |
|
(105) |
|
(106) |
En la Figura 9 se presenta de forma gráfica el campo de desplazamiento vertical en toda la viga.
Cálculo de las fuerzas que el suelo le hace a la viga
La fuerza que el suelo ejerce sobre la viga se calcula a partir del campo de desplazamientos de la siguiente forma:
|
(107) |
Donde:
|
(108) |
|
(109) |
|
(110) |
| (a) Campo de desplazamientos. | (b) Campo de la fuerza que el suelo ejerce sobre la viga. |
| Figura 12: Campos de desplazamiento en la viga y fuerza que el suelo ejerce sobre la viga. | |
Como revisión las fuerzas que el suelo realiza sobre la viga y por ende del campo de desplazamientos, a continuación se revisará el equilibrio vertical y rotacional al rededor del punto 1 para toda la viga:
|
Cálculo de las fuerzas internas
A partir del campo de desplazamientos es posible obtener el campo de momento flector como:
|
(112) |
Donde:
|
(113) |
|
(114) |
|
(115) |
Mientras que el campo de fuerza cortante se calcula como:
|
(116) |
Donde:
|
(117) |
|
(118) |
|
(119) |
| (a) Fuerza cortante. | (b) Momento flector. |
| Figura 13: Campos de fuerzas internas. | |
En el siguiente link se encuentra el código de Python empleado para la solución de este ejemplo: https://drive.google.com/open?id=1yq_G2WYO3X6oUW3Nnlq4jb9ATUd_WvlP
Resolver el pórtico plano presentado en la Figura 14a cuyos elementos son rectangulares de base y altura iguales a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L/20}
y cuyo módulo de elasticidad es Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E}
.
| (a) Problema a resolver. | (b) Discretización. |
| Figura 14: Pórtico plano sometido a cargas externas oblicuas y discretización empleada para su solución. | |
Discretización
En la Figura 14b se presenta la discretización y ejes locales a emplear en la solución de este ejercicio.
Definición de la carga externa en coordenadas locales
A partir de la Figura 14a, es posible obtener el valor de la carga externa en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_A}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
|
(120.a) |
Mientras que la carga externa en dirección del eje Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {y'}_A}
en el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es:
|
(121.a) |
De forma similar, se tiene que las cargas externas en dirección de los ejes Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_B}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {y'}_B}
en el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
son respectivamente:
|
Cálculo del sistema de ecuaciones de cada elemento en coordenadas locales
A partir de lo presentado en 42 se tiene que las ecuaciones del método de rigidez en coordenadas locales para los elementos Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
son respectivamente:
|
(123) |
|
(124) |
Cálculo del sistema de ecuaciones de cada elemento en coordenadas globales
Las ecuaciones 123 y 123 se transforman a coordenadas globales como:
|
(125) |
|
(126) |
Cálculo de los desplazamientos nodales desconocidos
A partir del equilibrio de los nudos en las direcciones donde las fuerzas externas o internas son conocidas, se obtiene el siguiente sistema lineal de cindo ecuaciones con igual número de incógnitas:
|
(127) |
Cuya solución es:
|
(128) |
cálculo de las reacciones
Del equilibrio de los nodos en las direcciones donde se desconocen las reacciones se obtiene el siguiente sistema lineal de cuatro ecuaciones:
|
(129) |
Donde se emplearon los valores de los desplazamientos nodales presentados en 128.
Como revisión de las reacciones obtenidas anteriormente se realizará el equilibrio de toda la estructura:
|
Cálculo del campo de desplazamiento del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): A
Para poder calcular el campo de desplazamiento homogéneo del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es necesario haber transformado los desplazamientos en los extremos de este elementos de coordenadas globales a locales. Una vez realizado esto, se tiene que el campo de desplazamiento homogéneo en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_A}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es:
|
(131) |
Mientras que a partir de lo presentado en 19 y 8.3, se tiene que el campo de desplazamiento empotrado en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_A}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es:
|
(132) |
Con lo cual el campo de desplazamiento total en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_A}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
es:
|
(133) |
Procediendo de forma similar a como se realizó en el ejemplo 8.1, se tiene que los campos homogéneo, empotrado y total en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {y'}_A}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
son respectivamente:
|
(134) |
|
(135) |
|
(136) |
Cálculo del campo de desplazamiento del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B
Procediendo de forma similar a como se realizó con el elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A} , se tiene que los campos de desplazamiento homogéneo, empotrado y total en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {x'}_B}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
son respectivamente:
|
(137) |
|
(138) |
|
(139) |
Mientras que los campos de desplazamiento homogéneo, empotrado y total en dirección Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {y'}_B}
del elemento Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
son respectivamente:
|
(140) |
|
(141) |
|
(142) |
Como síntesis de los campos de desplazamiento calculados anteriormente en la Figura 15d se presenta la configuración deformada de la estructura.
Cálculo de las fuerzas internas
A partir de la derivación de los campos de desplazamientos en coordenadas locales del ambos elementos se obtiene las fuerzas internas en estos:
|
(143) |
|
(144) |
|
(145) |
|
En el siguiente link se encuentra el código de Python empleado para la solución de este ejemplo: https://drive.google.com/open?id=1WNnbqC-09Odkh3GajgONgi8Al1sKhRz_
Las ideas fundamentales de este artículo han nacido del los cursos de análisis estructural y mecánica estructural que el autor Juan Camilo Molina-Villegas a dictado en lo últimos años en la Universidad de Medellín y en la Universidad Nacional de Colombia - Sede Medellín, y de los cuales han sido estudiantes los otros dos coatores. Los comentarios de muchos de los estudiantes de dichos cursos han contribuido al nivel actual de maduración de las ideas presentadas en este artículo, por lo cual los autores les agradecen infinitamente.
A continuación se presentará una formulación “exacta” del M.E.F. la cual es equivalente a la presentada en la sección 4, la cual también se basa en la solución del P.V.F 4 pero ahora a partir de la llamada forma débil de 8.a en lugar de esta directamente. De forma similar, el campo de desplazamientos se expresara como:
|
(A.1) |
Donde Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
son las funciones de forma presentadas en 11 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_f(x)}
se calcula según lo presentado en 19.
Además, es fácil probar que las funciones de forma Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
cumplen la siguiente propiedad:
|
(A.2) |
Multiplicando a ambos lados de 8.a por una función de peso o ponderación Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)}
se obtiene:
|
(A.3) |
Ahora, integrando a ambos lados de A.3 con respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
se obtiene:
|
(A.4) |
Integrando por partes en el lado izquierdo de la anterior ecuación se obtiene:
|
(A.5) |
La cual luego de el primer término del lado izquierdo de la anterior ecuación y reordenándola, da lugar la siguiente ecuación:
|
(A.6) |
Ahora si se emplean como las funciones de ponderación o peso (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)} ) a las funciones de forma (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
) se tiene que:
|
(A.7) |
Lo anterior debido a que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u_f(0)=u_f(L)=0}
y a A.2.
Con lo cual, si como caso particular primero se emplea Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _1(x)}
en lugar de Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)}
y luego Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _4(x)}
, a partir de A-7 se obtiene respectivamente:
|
La cual es equivalente a la formulación del método de rigidez para elementos tipo barra de sección constante y material elástico lineal presentada en 23.
Ahora, de forma similar a como se realizó en el apéndice 11 para el elemento tipo barra, se procederá a realizar la formulación “exacta” del M.E.F. para un elemento tipo viga de sección transversal constante. Al igual que lo presentado en la sección 5, el objetivo es realizar esta formulación mediante la solución del P.V.F. 5 pero ahora empleando la forma de débil de la E.D. 24.a en lugar de esta directamente. A diferencial de la formulación tradicional del M.E.F., en esta alternativa el campo de desplazamientos no solo dependerá de los desplazamientos y rotaciones de los extremos del elemento (campo homogéneo) sino también de la carga externa (campo empotrado), es decir:
|
(B.1) |
Donde una propiedad importante de las funciones de forma Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _6(x)}
es que estas cumplen la siguiente ecuación diferencial:
|
(B.2) |
Multiplicando por una función de peso o ponderación Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)}
a ambos lados de 24.a se obtiene:
|
(B.3) |
E integrando con respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
a ambos lados de B.3 se obtiene:
|
(B.4) |
Resolviendo con integración por partes la integral del lado izquierdo de la anterior ecuación se obtiene:
|
(B.5) |
Ahora, teniendo en cuenta que la fuerza cortante se obtiene a partir del campo de desplazamiento vertical como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle V(x)=-EI \dfrac{d^3v}{dx^3}(x)} , se obtiene:
|
(B.6) |
Mientras que debido a que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle FY_i=-V(0)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle FY_j=V(L)}
, la anterior ecuación se reescribe como:
|
(B.7) |
Integrando por partes de nuevo, ahora en el tercer término del lado izquierdo de la ecuación anterior, se obtiene:
|
(B.8) |
Por último, teniendo en cuenta que el campo de momentos interno se calcula como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M(x)=EI\dfrac{d^2 v}{dx^2}(x)}
y que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_i=-M(0)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_j=M(L)}
, la ecuación anterior se reescribe como:
|
(B.9) |
Reemplazando B.1 en el lado derecho de B.9, se obtiene:
|
(B.10) |
Ahora, si como funciones de ponderación se emplean las funciones de forma, empleando integración por partes se prueba lo siguiente:
|
(B.11) |
Donde se ha empleado B.2 y los valores de Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \dfrac{dv_f}{dx}(x)}
son iguales a cero en Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=0}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
.
Con lo cual, si como caso particular primero se emplea Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi _6(x)}
, así como empleando B.11, a partir de B.10 se obtienen las siguientes cuatro ecuaciones:
|
Que son equivalentes a las presentadas en 41.
Para este caso el objetivo es resolver el P.V.F. 7 pero ahora empleando la forma débil de 43.a en lugar de esta directamente.
Como se presentó en 66, la solución de 43.a se expresa como la suma de un campo homogéneo Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_h(x)}
y un campo empotrado (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x}
), es decir:
|
(C.1) |
Donde el campo homogéneo se expresa en función de los desplazamientos y rotaciones en los extremos del elemento y de las funciones de forma Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _2(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _6(x)}
de la siguiente manera:
|
(C.2) |
Con lo cual, el campo de desplazamiento total, se expresa ahora como:
|
(C.3) |
Es importante resaltar que debido a que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_h(x)}
es una combinación lineal de las funciones de forma Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _2(x)}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _3(x)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _5(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Psi _6(x)}
, estas últimas cumplen la propiedad de ser solución de la ecuación diferencial homogénea asociada con 43.a, es decir:
|
(C.4) |
Multiplicando por una función de peso o ponderación Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)}
a ambos lados de 43.a se obtiene:
|
(C.5) |
E integrando con respecto a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x}
entre 0 y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
a ambos lados de C.5:
|
(C.6) |
Resolviendo con integración por partes la integral del lado izquierdo de la anterior ecuación, da como resultado:
|
(C.7) |
Ahora teniendo en cuenta que la fuerza cortante se obtiene a partir del campo de desplazamiento vertical como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle V(x)=-EI \dfrac{d^3v}{dx^3}(x)} , se obtiene:
|
(C.8) |
Teniendo en cuenta que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle FY_i=-V(0)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle FY_j=V(L)}
, la anterior ecuación se reescribe como:
|
(C.9) |
Integrando por partes de nuevo, ahora en el tercer término del lado izquierdo de la ecuación anterior, se obtiene:
|
(C.10) |
Por último, teniendo en cuenta que el campo de momentos interno se calcula como Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M(x)=EI\dfrac{d^2 v}{dx^2}(x)}
y que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_i=-M(0)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_j=M(L)}
, la ecuación anterior se reescribe como:
|
(C.11) |
Reemplazando 155 en el lado derecho de 163, se obtiene:
|
(C.12) |
Ahora, si como funciones de ponderación se emplean las funciones de forma, empleando integración por partes se prueba lo siguiente:
|
(C.13) |
Donde se ha empleado C.4 y que los valores de Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_f(x)}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \dfrac{dv_f}{dx}(x)}
son iguales a cero en Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=0}
y Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
.
Empleando C.13 en C.12 esta última toma la siguiente forma (donde para ser coherentes con C.13, se ha tomado que Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle w(x)=\Psi _l(x)} ):
|
(C.14) |
Reemplazando Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l}
por 2,3,5 y 6 en C.14 se obtienen cuatro ecuaciones, las cuales se expresan en forma matricial como:
|
(C.15) |
Donde:
|
(C.16) |
En este punto es importante resaltar que C.15 es equivalente a lo presentado en 66.
[1] Challis, Lawrie and Sheard, Fred. The green of Green functions, Physics Today, 56(12):41-46,2003.
[2] Duffy, Dean G. Green’s functions with applications. Chapman and Hall/CRC, 2015.
[3] Banerjee, PK and Butterfield, R. Boundary Element Methods in Engineering Science. McGraw-Hill, New York, 1981.
[4] Sánchez-Sesma, F. J. and Ramos-Martínez, J. and Campillo, M. An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves. Earthquake engineering & structural dynamics, 22(4):279-295, 1993.
[5] Fairweather, G. and Karageorghis, A. and Martin, P.A. The method of fundamental solutions for scattering and radiation problems, Engineering Analysis with Boundary Elements, 27(7) :759-769, 2003.
[6] Thomson, William T. Transmission of elastic waves through a stratified solid medium. Journal of applied Physics, 21(2):89-93, 1950.
[7] Boussinesq, Joseph. Application des potentiels a l'étude de l'équilibre et du mouvement des solides élastiques. Gauthier-Villars, 1885.
[8] Cerruti, Valentino. Ricerche intorno all'equilibrio de'corpi elastici isotropi: memoria. Coi tipi del Salviucci, 1882.
[9] Mindlin, R.D. Force at a Point in the Interior of a Semi-Infinite Solid. Physics, 7(5):195-202, 1936.
[10] George Gebriel Stokes. On dynamical theory of diffraction. Transactions of the Cambridge Philosophical Society, 9:1-48, 1849.
[11] H. Lamb. On the propagation of tremors over the surface of an elastic solid. Philosophical Transactions of the Royal Society of London, 203:1-42, 1904.
[12] Chao, C.C. Dynamical response of an elastic half-space to tangential surface loadings. Journal of Applied Mechanics, 27:559, 1960.
[13] Kausel, E. Fundamental solutions in elastodynamics: a compendium. Cambridge University Press, 2006.
[14] Aki, K. and Richards, P.G. Quantitative seismology. Univ Science Books, 2002.
[15] Colunga, Arturo Tena. Análisis de estructuras con métodos matriciales. Limusa, 2007.
[16] McCormac, Jack C. Structural Analysis: using classical and matrix methods. Wiley Hoboken, NJ, 2007.
[17] Kassimali, Aslam. Matrix analysis of structures SI version. Cengage Learning, 2012.
[18] J.N. Reddy. An introduction to the finite element method. McGrawHill, 2006.
[19] K.J. Bathe. Finite element procedures. McGrawHill, 2006.
[20] Hetényi, Miklós. Beams on elastic foundation: theory with applications in the fields of civil and mechanical engineering. University of Michigan, 1971.
[21] M. Eisemberger and D.Z. Yankelevsky. Exact stiffness matrix for beams on elastic foundation. Computer & Structures, 21(6):1355-1359, 1985.
Published on 23/09/20
Accepted on 14/09/20
Submitted on 01/01/20
Volume 36, Issue 3, 2020
DOI: 10.23967/j.rimni.2020.09.004
Licence: CC BY-NC-SA license
Are you one of the authors of this document?