(84 intermediate revisions by 4 users not shown) | |||
Line 1: | Line 1: | ||
− | |||
− | |||
− | |||
− | |||
==Resumen== | ==Resumen== | ||
− | 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 | + | 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 líneas de influencia y en algunos casos esporádicos en la formulación del método del trabajo virtual. En este artículo se presenta el método de rigidez con funciones de Green (Green functions stiffness method) el cual es una metodología novedosa para obtener la respuesta analítica o exacta de estructuras reticulares planas, y como su nombre lo indica mezcla al método de rigidez con las funciones de Green. 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 obtener la respuesta analítica ante 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. |
+ | '''Palabras clave''': Método de rigidez con funciones de Green, Funciones de Green, Método de rigidez, Método de elementos finitos, Vigas de Bernoulli, Pórticos planos, Vigas sobre fundación flexible, Funciones de forma | ||
+ | |||
+ | ==Abstract== | ||
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 the Green functions stiffness method, which is a novel 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 to compute the analytic reponse 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. | 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 the Green functions stiffness method, which is a novel 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 to compute the analytic reponse 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''': Green function stiffness method, Green functions, stiffness method, finite element method, Bernoulli beam, plane frames, beams on elastic foundations, shape functions |
− | ==1 Introducción== | + | ==1. Introducción== |
− | Las funciones de Green (F.G.) son la respuesta de un medio o elemento 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 sean lineales [ | + | Las funciones de Green (F.G.) son la respuesta de un medio o elemento 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 sean lineales <span id='citeF-1'></span>[[#cite-1|[1]]]. Esto es debido a que ellas por sí solas son solución de problemas fundamentales, pueden ser empleadas para resolver problemas con fuentes o fuerzas distribuidas <span id='citeF-2'></span>[[#cite-2|[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 <span id='citeF-3'></span>[[#cite-3|[3]]], el método indirecto de elementos de frontera <span id='citeF-4'></span>[[#cite-4|[4]]] o el método de las soluciones fundamentales <span id='citeF-5'></span>[[#cite-5|[5]]]. En geotecnia las principales F.G. son: la respuesta de un espacio completo o infinito <span id='citeF-6'></span>[[#cite-6|[6]]], la respuesta de un semi-espacio ante la acción de un fuerza en superficie <span id='citeF-7'></span><span id='citeF-8'></span>[[#cite-7|[7,8]]] y en el interior <span id='citeF-9'></span>[[#cite-9|[9]]]. En sismología las principales F.G. son: la respuesta de un espacio completo <span id='citeF-10'></span>[[#cite-10|[10]]], la respuesta de un semi-espacio debida a la acción de una fuerza normal <span id='citeF-11'></span>[[#cite-11|[11]]] y tangencial <span id='citeF-12'></span>[[#cite-12|[12]]], para el lector interesado en este tema se le recomiendan los libros <span id='citeF-13'></span><span id='citeF-14'></span>[[#cite-13|[13,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 líneas de influencia y en algunos casos esporádicos 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 <math display="inline">G_{ij}(\mathbf{x},\boldsymbol \xi )</math>, la cual corresponde al desplazamiento en dirección <math display="inline">i</math> del punto <math display="inline">\mathbf{x}</math>, debido a una fuerza puntual unitaria en dirección <math display="inline">j</math> aplicada en el punto <math display="inline">\boldsymbol \xi </math>. | Para el caso particular de la elasticidad y análisis estructural, las funciones de Green se expresan como <math display="inline">G_{ij}(\mathbf{x},\boldsymbol \xi )</math>, la cual corresponde al desplazamiento en dirección <math display="inline">i</math> del punto <math display="inline">\mathbf{x}</math>, debido a una fuerza puntual unitaria en dirección <math display="inline">j</math> aplicada en el punto <math display="inline">\boldsymbol \xi </math>. | ||
− | 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 más complejo de implementar y busca la solución exacta del modelo estructural en estudio, mientras que el segundo es más 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 de geometrías simples, como puntuales o distribuidas con variación lineal. Mientras que el M.E.F. suele emplearse para el análisis de 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 ([ | + | 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 más complejo de implementar y busca la solución exacta del modelo estructural en estudio, mientras que el segundo es más 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 de geometrías simples, como puntuales o distribuidas con variación lineal. Mientras que el M.E.F. suele emplearse para el análisis de 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 (<span id='citeF-15'></span><span id='citeF-16'></span><span id='citeF-17'></span>[[#cite-15|[15,16,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 <span id='citeF-18'></span><span id='citeF-19'></span>[[#cite-18|[18,19]]]. |
− | En este artículo se presenta una formulación particular del método de rigidez para obtener la respuesta analítica de problemas de análisis estructural formados por estructuras reticulares planas (pudiendo estar formadas por elementos tipo barra, viga, viga sobre fundación flexible o pórtico plano), la cual se llamará el método de rigidez con funciones de Green (Green functions stiffness method), y emplea a las funciones de Green para realizar el cálculo de los desplazamientos al interior de los elementos. 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 | + | En este artículo se presenta una formulación particular del método de rigidez para obtener la respuesta analítica de problemas de análisis estructural formados por estructuras reticulares planas (pudiendo estar formadas por elementos tipo barra, viga, viga sobre fundación flexible o pórtico plano), la cual se llamará el método de rigidez con funciones de Green (Green functions stiffness method), y emplea a las funciones de Green para realizar el cálculo de los desplazamientos al interior de los elementos. 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 cálculo 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 plano 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 concluirá con la realización tres ejemplos ilustrativos (para una viga, una viga sobre fundación flexible y un pórtico plano respectivamente). | 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 plano 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 concluirá con la realización tres ejemplos ilustrativos (para una viga, una viga sobre fundación flexible y un pórtico plano respectivamente). | ||
− | ==2 Ecuaciones diferenciales gobernantes== | + | ==2. Ecuaciones diferenciales gobernantes== |
− | A continuación se definirán las propiedades de los elementos tipo barra, viga, pórtico plano y viga sobre fundación flexible a estudiar en este documento, al igual que sus ecuaciones diferenciales gobernantes. Como generalidad, para todos estos tipos elementos se empleará la convención de fuerzas internas positiva presentada en la | + | A continuación se definirán las propiedades de los elementos tipo barra, viga, pórtico plano y viga sobre fundación flexible a estudiar en este documento, al igual que sus ecuaciones diferenciales gobernantes. Como generalidad, para todos estos tipos elementos se empleará la convención de fuerzas internas positiva presentada en la [[#img-1|Figura 1]], la cual se inspira en aquella usada para la definición de los esfuerzos positivos en la mecánica del medio continuo o de sólidos. |
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | |
+ | <div id='img-1'></div> | ||
+ | |||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;" | ||
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"| [[Image:Molina-Villegas_et_al_2020b-Figura01.png|420px|Convención positiva para las fuerzas internas. P(x) es la fuerza axial (dirección eje local x), V(x) es la fuerza cortante (dirección eje local y) y M(x) es el momento flector (dirección eje local z, perpendicular tanto a x como a y.)]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 1 | + | | colspan="1" style="padding:10px;"| '''Figura 1'''. Convención positiva para las fuerzas internas. <math>P(x)</math> es la fuerza axial (dirección eje local <math>x</math>), <math>V(x)</math> es la fuerza cortante (dirección eje local <math>y</math>) y <math>M(x)</math> es el momento flector (dirección eje local <math>z</math>, perpendicular tanto a <math>x</math> como a <math>y</math>) |
|} | |} | ||
+ | |||
===2.1 Elemento tipo barra=== | ===2.1 Elemento tipo barra=== | ||
− | 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 (<math display="inline">p(x)</math>), la cual actúa en dirección del eje <math display="inline">x</math>, y su valor es positivo si se dirige en la dirección de dicho eje local ( | + | 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 (<math display="inline">p(x)</math>), la cual actúa en dirección del eje <math display="inline">x</math>, y su valor es positivo si se dirige en la dirección de dicho eje local ([[#img-2|Figura 2]]). Para el caso que el elemento sea de sección transversal constante con área de la sección transversal <math display="inline">A</math>, material elástico lineal con módulo de elasticidad <math display="inline">E</math>, su ecuación diferencial gobernante es: |
<span id="eq-1"></span> | <span id="eq-1"></span> | ||
Line 47: | Line 50: | ||
|} | |} | ||
− | + | donde <math display="inline">u(x)</math> es el desplazamiento en dirección axial del elemento (eje local <math display="inline">x</math>). | |
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | |
+ | <div id='img-2'></div> | ||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;" | ||
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura02.png|420px|Elemento tipo barra y el sistema coordenado local empleado para la definición de su E.D. gobernante.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 2 | + | | colspan="1" style="padding:10px;"| '''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 <math display="inline">(P)</math> se calcula como: | Mientras que a partir de la ley de Hooke unidimensional su fuerza interna axial <math display="inline">(P)</math> se calcula como: | ||
Line 67: | Line 74: | ||
|} | |} | ||
− | + | donde <math display="inline">\sigma (x)</math> es el esfuerzo axial y <math display="inline">\epsilon (x)</math> es la deformación unitaria axial. | |
===2.2 Elemento tipo viga (viga de Euler)=== | ===2.2 Elemento tipo viga (viga de Euler)=== | ||
− | El elemento tipo viga está sometido tanto a fuerza cortante en dirección del eje local <math display="inline">y</math>, como a momento flector alrededor del eje local <math display="inline">z</math> (perpendicular a los ejes <math display="inline">x</math> y <math display="inline">y</math>, | + | El elemento tipo viga está sometido tanto a fuerza cortante en dirección del eje local <math display="inline">y</math>, como a momento flector alrededor del eje local <math display="inline">z</math> (perpendicular a los ejes <math display="inline">x</math> y <math display="inline">y</math>, [[#img-3|Figura 3]]). Sus cargas externas se definen en términos de la carga por unidad de longitud en dirección del eje local <math display="inline">y</math> (<math display="inline">q(x)</math>), la cual es positiva en la dirección de este eje ([[#img-3|Figura 3]]). Si el elemento es de sección transversal constante, con momento de inercia <math display="inline">I</math> y material elástico lineal con módulo de elasticidad <math display="inline">E</math>, su ecuación diferencial gobernante es: |
<span id="eq-3"></span> | <span id="eq-3"></span> | ||
Line 84: | Line 91: | ||
|} | |} | ||
− | + | donde <math display="inline">v(x)</math> es el desplazamiento en dirección del eje local <math display="inline">y</math>, es decir, perpendicular al eje longitudinal del elemento. | |
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | |
+ | |||
+ | <div id='img-3'></div> | ||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;" | ||
|- | |- | ||
− | |[[ | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura03.png|420px|Elemento tipo viga y el sistema coordenado local empleado para la definición de su E.D. gobernante.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 3 | + | | colspan="1" style="padding:10px;"| '''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 <math display="inline">v(x)</math>, respectivamente como: | + | |
+ | |||
+ | Mientras que las fuerzas internas cortante y momento flector se calculan a partir del campo de desplazamiento <math display="inline">v(x)</math>, respectivamente como: | ||
+ | |||
+ | <span id="eq-4"></span> | ||
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 99: | Line 113: | ||
|- | |- | ||
| style="text-align: center;" | <math>V(x)=-EI \dfrac{d^3 v}{dx^3}(x) </math> | | style="text-align: center;" | <math>V(x)=-EI \dfrac{d^3 v}{dx^3}(x) </math> | ||
+ | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (4.a) | | style="width: 5px;text-align: right;white-space: nowrap;" | (4.a) | ||
+ | |- | ||
+ | | | ||
+ | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
| style="text-align: center;" | <math> M(x)= EI \dfrac{d^2 v}{dx^2}(x) </math> | | style="text-align: center;" | <math> M(x)= EI \dfrac{d^2 v}{dx^2}(x) </math> | ||
+ | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (4.b) | | style="width: 5px;text-align: right;white-space: nowrap;" | (4.b) | ||
− | |||
|} | |} | ||
===2.3 Elemento tipo viga sobre fundación flexible=== | ===2.3 Elemento tipo viga sobre fundación flexible=== | ||
− | 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 <math display="inline">k</math> y genera una fuerza por unidad de longitud sobre la viga igual a <math display="inline">-kv(x)</math> ( | + | 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 <math display="inline">k</math> y genera una fuerza por unidad de longitud sobre la viga igual a <math display="inline">-kv(x)</math> ([[#img-4|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 <span id='citeF-20'></span>[[#cite-20|[20]]]: |
<span id="eq-5"></span> | <span id="eq-5"></span> | ||
Line 122: | Line 140: | ||
<div id='img-4'></div> | <div id='img-4'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura04.png|420px|Elemento tipo viga sobre fundación flexible y el sistema coordenado local empleado para la definición de su E.D. gobernante.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 4 | + | | colspan="1" style="padding:10px;"| '''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, para este elemento las fuerzas internas se calculan a partir del campo de desplazamientos empleando [[# | + | |
+ | |||
+ | Al igual que para los elementos tipo viga, para este elemento las fuerzas internas se calculan a partir del campo de desplazamientos empleando [[#eq-4|(4)]]. | ||
===2.4 Elemento tipo pórtico plano=== | ===2.4 Elemento tipo pórtico plano=== | ||
− | 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 local <math display="inline">z</math>. Sus ecuaciones diferenciales gobernantes son [[#eq-1|1]] para el desplazamiento axial y [[#eq-3|3]] para el desplazamiento perpendicular al eje longitudinal del elemento (eje <math display="inline">y</math>) y sus fuerzas internas se calculan a partir de los campos de desplazamiento empleando [[#eq-2|2]] y [[# | + | 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 local <math display="inline">z</math>. Sus ecuaciones diferenciales gobernantes son [[#eq-1|(1)]] para el desplazamiento axial y [[#eq-3|(3)]] para el desplazamiento perpendicular al eje longitudinal del elemento (eje <math display="inline">y</math>) y sus fuerzas internas se calculan a partir de los campos de desplazamiento empleando [[#eq-2|(2)]] y [[#eq-4|(4)]]. |
− | ==3 Metodología== | + | ==3. Metodología== |
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 puntos (lo cual comúnmente se conoce como grados de libertad). Esta formulación se expresa matricialmente en un sistema local de coordenadas como: | 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 puntos (lo cual comúnmente se conoce como grados de libertad). Esta formulación se expresa matricialmente en un sistema local de coordenadas como: | ||
Line 149: | Line 169: | ||
|} | |} | ||
− | + | donde <math display="inline">\{ {F'}_E \} </math> es el vector de fuerzas en los extremos del elemento en coordenadas locales, <math display="inline">{[K'}_E]</math> es la matriz de rigidez en coordenadas locales, <math display="inline">\{ \Delta' _E \} </math> es el vector de desplazamientos en los extremos del elemento en coordenadas locales y <math display="inline">\{ \left.{F'}_E \right.^f \} </math> es el vector de fuerzas de empotramiento en coordenadas locales. | |
− | De analizar los dos términos del lado derecho de [[#eq-6|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 para los elementos estudiados se expresan como: | + | De analizar los dos términos del lado derecho de [[#eq-6|(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 para los elementos estudiados se expresan como: |
<span id="eq-7"></span> | <span id="eq-7"></span> | ||
Line 164: | Line 184: | ||
|} | |} | ||
− | + | donde <math display="inline">d(x)</math> es el componente del campo de desplazamiento, <math display="inline">d_h(x)</math> el omponente del campo homogéneo de desplazamiento y | |
+ | <math display="inline">d_f(x)</math> el omponente del campo empotrado de desplazamiento. | ||
− | + | A continuación, a partir de [[#eq-7|(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 con funciones de Green 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. | |
− | + | ==4. Formulación analítica del método de rigidez con funciones de Green para elementos tipo barra== | |
− | + | 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 <math display="inline">L</math> es (Figura [[#img-2|(2)]]): | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | 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 <math display="inline">L</math> es ( | + | |
+ | <span id="eq-8"></span> | ||
<span id="eq-8.a"></span> | <span id="eq-8.a"></span> | ||
<span id="eq-8.b"></span> | <span id="eq-8.b"></span> | ||
Line 197: | Line 213: | ||
|} | |} | ||
− | + | donde <math display="inline">u(0)=u_i</math> es el desplazamiento en dirección axial (eje local <math display="inline">x</math>) en el extremo inicial del elemento y <math display="inline">u(L)=u_j</math> es el desplazamiento en la misma dirección pero en el extremo final del elemento. | |
− | Para este caso la particularización de [[#eq-7|7]] es: | + | Para este caso la particularización de [[#eq-7|(7)]] es: |
<span id="eq-9"></span> | <span id="eq-9"></span> | ||
Line 212: | Line 228: | ||
|} | |} | ||
− | + | donde <math display="inline">u_h(x)</math> es el campo de desplazamiento homogéneo (sección [[#4.1 Solución del P.V.F. homogéneo|4.1]]), mientras que <math display="inline">u_f(x)</math> es el campo empotrado (sección [[#4.2 Solución del P.V.F. empotrado|4.2]]). | |
===4.1 Solución del P.V.F. homogéneo=== | ===4.1 Solución del P.V.F. homogéneo=== | ||
− | El P.V.F. que gobierna al campo homogéneo se presenta en [[# | + | El P.V.F. que gobierna al campo homogéneo se presenta en [[#eq-10|(10)]] y su nombre se debe a que [[#eq-10.a|(10.a)]] es una ecuación diferencial homogénea. <span id="eq-10"></span> |
<span id="eq-10.a"></span> | <span id="eq-10.a"></span> | ||
Line 237: | Line 253: | ||
|} | |} | ||
− | Es fácil probar que la solución de [[# | + | Es fácil probar que la solución de [[#eq-10|(10)]] es el siguiente polinomio de primer grado: |
<span id="eq-11"></span> | <span id="eq-11"></span> | ||
Line 250: | Line 266: | ||
|} | |} | ||
− | + | donde: <span id="eq-12"></span> | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 265: | Line 281: | ||
|} | |} | ||
− | <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math> son conocidas como las funciones de forma “exactas” de este problema y en este caso son iguales a aquellas que suelen emplearse en el M.E.F. para interpolar el campo de desplazamientos en un elemento tipo barra. Un par de propiedades muy importantes de <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math>, las cuales nos permitirán expresar de forma compacta algunos cálculos posteriores, son: | + | <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math> son conocidas como las funciones de forma “exactas” de este problema y en este caso son iguales a aquellas que suelen emplearse en el M.E.F. para interpolar el campo de desplazamientos en un elemento tipo barra. Un par de propiedades muy importantes de <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math>, las cuales nos permitirán expresar de forma compacta algunos cálculos posteriores, son: <span id="eq-13"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 280: | Line 296: | ||
|} | |} | ||
− | En este punto es importante resaltar que en el M.E.F. es usual emplear [[#eq-11|11]] para aproximar el campo total de desplazamientos, lo cual, a partir de [[#eq-9|9]], se observa que es inexacto pues no incluye al campo empotrado. | + | En este punto es importante resaltar que en el M.E.F. es usual emplear [[#eq-11|(11)]] para aproximar el campo total de desplazamientos, lo cual, a partir de [[#eq-9|(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: | Por su parte, empleando la ley de Hooke uniaxial, se tiene que la fuerza axial homogénea se puede calcular como: | ||
Line 295: | Line 311: | ||
|} | |} | ||
− | Mientras que las fuerzas en dirección del eje local <math display="inline">x</math> en los extremos de la barra son: | + | Mientras que las fuerzas en dirección del eje local <math display="inline">x</math> en los extremos de la barra son: <span id="eq-15"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 312: | Line 328: | ||
===4.2 Solución del P.V.F. empotrado=== | ===4.2 Solución del P.V.F. empotrado=== | ||
− | El problema de P.V.F. que gobierna al campo empotrado se presenta en [[# | + | El problema de P.V.F. que gobierna al campo empotrado se presenta en [[#eq-16|(16)]]. Su nombre se debe a que representa la respuesta de una barra doblemente empotrada sometida a una carga externa por unidad de longitud igual a la del elemento en estudio ([[#img-5a|Figura 5a]]) <span id="eq-16"></span> |
<span id="eq-16.a"></span> | <span id="eq-16.a"></span> | ||
Line 336: | Line 352: | ||
<div id='img-5b'></div> | <div id='img-5b'></div> | ||
<div id='img-5'></div> | <div id='img-5'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura05a.png|360px|Carga distribuida arbitraria.]] |
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura05c.png|360px|Fuerza puntual unitaria.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
| (a) Carga distribuida arbitraria. | | (a) Carga distribuida arbitraria. | ||
| (b) Fuerza puntual unitaria. | | (b) Fuerza puntual unitaria. | ||
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 5 | + | | colspan="2" style="padding:10px;"| '''Figura 5'''. Barras de sección transversal constante doblemente empotradas, sometidas a una carga distribuida arbitraria y una carga puntual unitaria |
|} | |} | ||
− | Para la solución del P.V.F. [[# | + | |
+ | |||
+ | Para la solución del P.V.F. [[#eq-16|(16)]] se empleará la función de Green de un elemento tipo barra doblemente empotrada, que es de sección transversal constante y material elástico lineal ([[#img-5b|Figura 5b]]). Esta corresponde a la respuesta (campo de desplazamiento axial) debido a la aplicación de una carga puntual unitaria ubicada a una distancia <math display="inline">\xi </math> 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 <math display="inline">(x=\xi )</math>): <span id="eq-17"></span> | ||
<span id="eq-17.a"></span> | <span id="eq-17.a"></span> | ||
Line 371: | Line 389: | ||
|} | |} | ||
− | + | donde <math display="inline">\delta (x)</math> es la función delta de Dirac y lo superíndices - y + significan a la izquierda y a la derecha del punto en cuestión. | |
− | La solución de [[# | + | La solución de [[#eq-17|(17)]] expresada en términos de las funciones de forma [[#eq-12|(12)]], se puede escribir de forma compacta como: |
<span id="eq-18"></span> | <span id="eq-18"></span> | ||
Line 386: | Line 404: | ||
|} | |} | ||
− | + | donde, para su escritura se han explotado las propiedades de simetría de las funciones de forma presentadas en [[#eq-13|(13)]], así como las siguientes propiedades de simetría la función de Green: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 414: | Line 432: | ||
|} | |} | ||
− | Multiplicando a ambos lados de [[#eq-17.a|17.a]] por <math display="inline">p(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: | + | Multiplicando a ambos lados de [[#eq-17.a|(17.a)]] por <math display="inline">p(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: |
<span id="eq-21"></span> | <span id="eq-21"></span> | ||
Line 429: | Line 447: | ||
|} | |} | ||
− | + | donde se empleó la siguiente propiedad de la función delta de Dirac: <math display="inline">p(x)=\int _0^L \delta (x-\xi ) p(\xi ) d\xi </math> (<math display="inline">0<\xi{<}L</math>). | |
− | De comparar [[#eq-21|21]] y [[#eq-16.a|16.a]], así como empleando [[#eq-20|20]] se concluye que: | + | De comparar [[#eq-21|(21)]] y [[#eq-16.a|(16.a)]], así como empleando [[#eq-20|(20)]] se concluye que: |
<span id="eq-22"></span> | <span id="eq-22"></span> | ||
Line 444: | Line 462: | ||
|} | |} | ||
− | La cual también puede interpretarse mediante el principio de superposición pues el desplazamiento en cada punto <math display="inline">x</math> de la barra doblemente empotrada es la suma de los desplazamientos causados en ese mismo punto por infinitas cargas puntuales ubicadas en <math display="inline">\xi </math> desde 0 hasta <math display="inline">L</math> y cuyo valor es <math display="inline">q(\xi )d\xi </math>, donde los limites de ambas integrales se explican por [[#eq-20|20]]. | + | La cual también puede interpretarse mediante el principio de superposición pues el desplazamiento en cada punto <math display="inline">x</math> de la barra doblemente empotrada es la suma de los desplazamientos causados en ese mismo punto por infinitas cargas puntuales ubicadas en <math display="inline">\xi </math> desde 0 hasta <math display="inline">L</math> y cuyo valor es <math display="inline">q(\xi )d\xi </math>, donde los limites de ambas integrales se explican por [[#eq-20|(20)]]. |
− | Empleando la ley de Hooke uniaxial, a partir de [[#eq-22|22]] el campo de fuerza axial empotrado se calcula como: | + | Empleando la ley de Hooke uniaxial, a partir de [[#eq-22|(22)]] el campo de fuerza axial empotrado se calcula como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 464: | Line 482: | ||
|} | |} | ||
− | Y las reacciones o fuerzas de empotramiento en dirección del eje local <math display="inline">x</math> son: | + | Y las reacciones o fuerzas de empotramiento en dirección del eje local <math display="inline">x</math> son: <span id="eq-24"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 479: | Line 497: | ||
|} | |} | ||
− | Es de destacar que la anterior forma de 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., la cual se encuentra en el Apéndice [[# | + | Es de destacar que la anterior forma de 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., la cual se encuentra en el Apéndice [[#A Formulación exacta del M.E.F. para un elemento tipo barra|A]]. |
===4.3 Superposición del campo homogéneo y el campo empotrado=== | ===4.3 Superposición del campo homogéneo y el campo empotrado=== | ||
− | Reemplazando [[#eq-11|11]] y [[#eq-22|22]] en [[#eq-9|9]] se tiene que el campo de desplazamiento en el elemento, expresado en función de los desplazamientos nodales y las cargas externas es: | + | Reemplazando [[#eq-11|(11)]] y [[#eq-22|(22)]] en [[#eq-9|(9)]] se tiene que el campo de desplazamiento en el elemento, expresado en función de los desplazamientos nodales y las cargas externas es: |
<span id="eq-25"></span> | <span id="eq-25"></span> | ||
Line 496: | Line 514: | ||
|} | |} | ||
− | Mientras que a partir de sumar [[# | + | Mientras que a partir de sumar [[#eq-15|(15)]] y [[#eq-24|(24)]], se tiene que las fuerzas en dirección del eje local <math display="inline">x</math> en los extremos del elemento, se expresan en términos de los desplazamientos en esos mismos puntos y en la misma dirección axial: |
<span id="eq-26"></span> | <span id="eq-26"></span> | ||
Line 509: | Line 527: | ||
|} | |} | ||
− | La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo barra y es equivalente a la presentada en el Apéndice [[# | + | La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo barra y es equivalente a la presentada en el Apéndice [[#A Formulación exacta del M.E.F. para un elemento tipo barra|A]] 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. [[# | + | Para finalizar es importante resaltar que las definiciones de los P.V.F. [[#eq-10|(10)]] y [[#eq-16|(16)]] no han sido arbitrarias y son totalmente coherentes con [[#eq-8|(8)]] y [[#eq-9|(9)]]. Por ejemplo si se reemplaza [[#eq-9|(9)]] en [[#eq-8|(8)]] y se emplean [[#eq-16|(16)]], se obtiene directamente el P.V.F. [[#eq-10|(10)]], lo cual también ocurre para los demás tipos de elementos analizados en este artículo. |
− | ==5 Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga== | + | ==5. Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga== |
− | El P.V.F. que gobierna la formulación del método de rigidez de un elemento tipo viga de sección transversal constante con momento de inercia <math display="inline">I</math>, material elástico lineal homogéneo con módulo de elasticidad <math display="inline">E</math> y longitud <math display="inline">L</math> es ( | + | El P.V.F. que gobierna la formulación del método de rigidez de un elemento tipo viga de sección transversal constante con momento de inercia <math display="inline">I</math>, material elástico lineal homogéneo con módulo de elasticidad <math display="inline">E</math> y longitud <math display="inline">L</math> es (Figura [[#img-3|3]]): <span id="eq-27"></span> |
<span id="eq-27.a"></span> | <span id="eq-27.a"></span> | ||
Line 544: | Line 562: | ||
|} | |} | ||
− | + | donde <math display="inline">v(0)=v_i</math> es el desplazamiento en dirección del eje local <math display="inline">y</math> en extremo inicial del elemento, <math display="inline">\dfrac{dv}{dx}(0)=\theta _i</math> es la pendiente en ese mismo punto, <math display="inline">v(L)=v_j</math> es el desplazamiento en dirección del eje local <math display="inline">y</math> en extremo final del elemento y <math display="inline">\dfrac{dv}{dx}(L)=\theta _j</math> es la rotación en ese último punto. | |
− | Para este caso la particularización de [[#eq-7|7]] es: | + | Para este caso la particularización de [[#eq-7|(7)]] es: |
<span id="eq-28"></span> | <span id="eq-28"></span> | ||
Line 559: | Line 577: | ||
|} | |} | ||
− | + | donde <math display="inline">v_h(x)</math> es el campo de desplazamiento homogéneo (sección [[#5.1 Solución del P.V.F. homogéneo|5.1]]), mientras que <math display="inline">v_f(x)</math> es el campo empotrado (sección [[#5.2 Solución del P.V.F. empotrado|5.2]]). | |
===5.1 Solución del P.V.F. homogéneo=== | ===5.1 Solución del P.V.F. homogéneo=== | ||
− | Para este caso el P.V.F. gobernante es: | + | Para este caso el P.V.F. gobernante es: <span id="eq-29"></span> |
<span id="eq-29.a"></span> | <span id="eq-29.a"></span> | ||
Line 605: | Line 623: | ||
|} | |} | ||
− | + | donde <math display="inline">\psi _2(x)</math>, <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _6(x)</math> se conocen coma las funciones de forma “exactas” de este problema y tienen el siguiente valor: <span id="eq-31"></span> | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 626: | Line 644: | ||
|} | |} | ||
− | Fuera de lo anterior, las siguientes propiedades de simetría de las funciones de forma <math display="inline">\psi _2(x)</math>, <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _2(x)</math> son muy importantes de resaltar: | + | Fuera de lo anterior, las siguientes propiedades de simetría de las funciones de forma <math display="inline">\psi _2(x)</math>, <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _2(x)</math> son muy importantes de resaltar: <span id="eq-32"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 647: | Line 665: | ||
|} | |} | ||
− | Al igual que para el caso de las barras, las funciones de forma “exactas” presentadas en [[# | + | Al igual que para el caso de las barras, las funciones de forma “exactas” presentadas en [[#eq-31|(31)]] 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 [[#eq-28|(28)]]. |
− | A partir de [[#eq-30|30]] y [[# | + | A partir de [[#eq-30|(30)]] y [[#eq-4|(4)]], las fuerzas internas homogéneas son: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 664: | Line 682: | ||
|} | |} | ||
− | Mientras que las fuerzas homogéneas en los extremos del elemento en dirección de los ejes locales <math display="inline">y</math>, <math display="inline">z</math> son: | + | Mientras que las fuerzas homogéneas en los extremos del elemento en dirección de los ejes locales <math display="inline">y</math>, <math display="inline">z</math> son: <span id="eq-34"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 687: | Line 705: | ||
===5.2 Solución del P.V.F. empotrado=== | ===5.2 Solución del P.V.F. empotrado=== | ||
− | El problema de P.V.F. que gobierna al campo empotrado se presenta en [[# | + | El problema de P.V.F. que gobierna al campo empotrado se presenta en [[#eq-35|(35)]] y corresponde a la respuesta de una viga doblemente empotrada sometida a una carga externa genérica definida por la función <math display="inline">q(x)</math> ([[#img-6a|Figura 6a]]) |
+ | |||
+ | <span id="eq-35"></span> | ||
<span id="eq-35.a"></span> | <span id="eq-35.a"></span> | ||
Line 719: | Line 739: | ||
<div id='img-6b'></div> | <div id='img-6b'></div> | ||
<div id='img-6'></div> | <div id='img-6'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura06a.png|360px|Carga distribuida arbitraria.]] |
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura06c.png|360px|Fuerza puntual unitaria.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
| (a) Carga distribuida arbitraria. | | (a) Carga distribuida arbitraria. | ||
| (b) Fuerza puntual unitaria. | | (b) Fuerza puntual unitaria. | ||
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 6 | + | | colspan="2" style="padding:10px;"| '''Figura 6'''. Vigas de sección transversal constante doblemente empotradas, sometidas a una carga distribuida arbitraria y una carga puntual unitaria |
|} | |} | ||
− | Para solucionar el P.V.F. [[# | + | |
+ | |||
+ | Para solucionar el P.V.F. [[#eq-35|(35)]] se definirá la función de 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 <math display="inline">\xi </math> del extremo izquierdo del elemento ([[#img-6b|Figura 6b]]). Esta función de Green está gobernada por el P.V.F. [[#eq-36|(36)]], el cual comprende el cumplimiento de la ecuación diferencial gobernante, condiciones de frontera, condiciones de continuidad de desplazamientos y rotaciones en (<math display="inline">x=\xi </math>) y equilibrio vertical y rotacional en este mismo punto: <span id="eq-36"></span> | ||
<span id="eq-36.a"></span> | <span id="eq-36.a"></span> | ||
Line 766: | Line 788: | ||
|} | |} | ||
− | La solución de [[# | + | La solución de [[#eq-36|(36)]], es decir, la función de Green para la viga doblemente empotrada presentada en la [[#img-6b|Figura 6b]] es: |
<span id="eq-37"></span> | <span id="eq-37"></span> | ||
Line 774: | Line 796: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>G_{yy}(x,\xi )= \begin{cases}\dfrac{L^3}{6EI} \left[-\left(\dfrac{x}{L} \right)^3 \psi _2(\xi ) +3 \left(\dfrac{x}{L}\right)^2 \dfrac{\psi _3(\xi )}{L} \right] & 0 < x \leq \xi \\[0.3cm] \dfrac{L^3}{6EI} \left[-\left(1-\dfrac{x}{L} \right)^3 \psi _5(\xi ) -3 \left(1-\dfrac{x}{L}\right)^2 \dfrac{\psi _6(\xi )}{L} \right]\qquad & \xi \leq x | + | | style="text-align: center;" | <math>G_{yy}(x,\xi )= \begin{cases}\dfrac{L^3}{6EI} \left[-\left(\dfrac{x}{L} \right)^3 \psi _2(\xi ) +3 \left(\dfrac{x}{L}\right)^2 \dfrac{\psi _3(\xi )}{L} \right] & 0 < x \leq \xi \\[0.3cm] \dfrac{L^3}{6EI} \left[-\left(1-\dfrac{x}{L} \right)^3 \psi _5(\xi ) -3 \left(1-\dfrac{x}{L}\right)^2 \dfrac{\psi _6(\xi )}{L} \right]\qquad & \xi \leq x < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (37) | | style="width: 5px;text-align: right;white-space: nowrap;" | (37) | ||
|} | |} | ||
− | La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en [[# | + | La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en [[#eq-32|(32)]] y las siguientes propiedades de simetría de la función de Green: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 820: | Line 842: | ||
|} | |} | ||
− | Multiplicando a ambos lados de [[#eq-36.a|36.a]] por <math display="inline">q(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: | + | Multiplicando a ambos lados de [[#eq-36.a|(36.a)]] por <math display="inline">q(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: |
<span id="eq-41"></span> | <span id="eq-41"></span> | ||
Line 835: | Line 857: | ||
|} | |} | ||
− | + | donde se empleó que <math display="inline">q(x)=\int _0^L \delta (x-\xi ) q(\xi ) d\xi </math> <math display="inline">(0<\xi{<}L)</math>. | |
− | De comparar [[#eq-21|21]] y [[#eq-35.a|35.a]], así como empleando [[#eq-40|40]], se concluye que: | + | De comparar [[#eq-21|(21)]] y [[#eq-35.a|(35.a)]], así como empleando [[#eq-40|(40)]], se concluye que: |
<span id="eq-42"></span> | <span id="eq-42"></span> | ||
Line 850: | Line 872: | ||
|} | |} | ||
− | Al igual que para el caso de la barra, el campo empotrado presentado en [[#eq-42|42]] puede interpretarse a partir del principio de superposición como la respuesta en un punto <math display="inline">x</math> debido a la suma de los aportes de los diferenciales de fuerza <math display="inline">q(\xi )d\xi </math> ubicados a lo largo de toda la viga, para lo cual se emplea [[#eq-40|40]]. | + | Al igual que para el caso de la barra, el campo empotrado presentado en [[#eq-42|(42)]] puede interpretarse a partir del principio de superposición como la respuesta en un punto <math display="inline">x</math> debido a la suma de los aportes de los diferenciales de fuerza <math display="inline">q(\xi )d\xi </math> ubicados a lo largo de toda la viga, para lo cual se emplea [[#eq-40|(40)]]. |
Mientras que los campos de fuerzas internas se calculan como: | Mientras que los campos de fuerzas internas se calculan como: | ||
Line 886: | Line 908: | ||
|} | |} | ||
− | Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales <math display="inline">y</math> y <math display="inline">z</math> son: | + | Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales <math display="inline">y</math> y <math display="inline">z</math> son: <span id="eq-45"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 907: | Line 929: | ||
|} | |} | ||
− | Por su parte para el caso en que se presente un momento puntual, pese a que este se puede expresar como una función <math display="inline">q(x)</math> empleando el límite de dos funciones delta de Dirac, por comodidad a continuación se presenta la función de Green para este caso ( | + | Por su parte para el caso en que se presente un momento puntual, pese a que este se puede expresar como una función <math display="inline">q(x)</math> empleando el límite de dos funciones delta de Dirac, por comodidad a continuación se presenta la función de Green para este caso ([[#img-7|Figura 7]]): |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 919: | Line 941: | ||
|} | |} | ||
− | + | con lo cual: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 953: | Line 975: | ||
|} | |} | ||
− | + | con lo cual las fuerzas de empotramiento para este caso son: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 975: | Line 997: | ||
<div id='img-7'></div> | <div id='img-7'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura11.png|420px|Viga de sección transversal constante doblemente empotrada, sometidas a un momento puntual unitario.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 7 | + | | colspan="1" style="padding:10px;"| '''Figura 7'''. Viga de sección transversal constante doblemente empotrada, sometidas a un momento puntual unitario |
|} | |} | ||
+ | |||
+ | |||
Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sinó de valor <math display="inline">M</math>, los resultados anteriores se deberán simplemente multiplicar por <math display="inline">M</math>, y en caso de que sea un momento por unidad de longitud distribuido, este se deberá integrar de forma similar a como se realizó con la carga <math display="inline">q(\xi )</math> pero ahora empleando <math display="inline">G_{y\theta }(x,\xi )</math> en lugar de <math display="inline">G_{yy}(x,\xi )</math>. | Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sinó de valor <math display="inline">M</math>, los resultados anteriores se deberán simplemente multiplicar por <math display="inline">M</math>, y en caso de que sea un momento por unidad de longitud distribuido, este se deberá integrar de forma similar a como se realizó con la carga <math display="inline">q(\xi )</math> pero ahora empleando <math display="inline">G_{y\theta }(x,\xi )</math> en lugar de <math display="inline">G_{yy}(x,\xi )</math>. | ||
===5.3 Superposición del campo homogéneo y el campo empotrado=== | ===5.3 Superposición del campo homogéneo y el campo empotrado=== | ||
− | Reemplazando [[#eq-30|30]] y [[#eq-42|42]] en [[#eq-28|28]] se obtiene que el campo de desplazamientos puede expresarse en función de los desplazamientos nodales y las fuerzas externas como: | + | Reemplazando [[#eq-30|(30)]] y [[#eq-42|(42)]] en [[#eq-28|(28)]] se obtiene que el campo de desplazamientos puede expresarse en función de los desplazamientos nodales y las fuerzas externas como: |
<span id="eq-51"></span> | <span id="eq-51"></span> | ||
Line 998: | Line 1,022: | ||
|} | |} | ||
− | Mientras que las fuerzas en los extremos del elemento se expresan en términos de los desplazamiento y rotaciones en esos mismos puntos a partir de la suma de [[# | + | Mientras que las fuerzas en los extremos del elemento se expresan en términos de los desplazamiento y rotaciones en esos mismos puntos a partir de la suma de [[#eq-34|(34)]] y [[#eq-45|(45)]]: |
<span id="eq-52"></span> | <span id="eq-52"></span> | ||
Line 1,011: | Line 1,035: | ||
|} | |} | ||
− | La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo viga y es equivalente a la formulación “exacta” del M.E.F. presentada en el Apéndice [[# | + | La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo viga y es equivalente a la formulación “exacta” del M.E.F. presentada en el Apéndice [[#B Formulación exacta del M.E.F. para un elemento tipo viga|B]]. |
− | ==6 Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga sobre fundación flexible== | + | ==6. Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga sobre fundación flexible== |
− | El P.V.F. que gobierna la formulación del método de rigidez en coordenadas locales 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 | + | El P.V.F. que gobierna la formulación del método de rigidez en coordenadas locales 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 constante e igual <math display="inline">k</math> y longitud <math display="inline">L</math> es (Figura [[#img-4|4]]): <span id="eq-53"></span> |
<span id="eq-53.a"></span> | <span id="eq-53.a"></span> | ||
Line 1,044: | Line 1,068: | ||
|} | |} | ||
− | + | donde se ha empleado como ecuación diferencial una versión equivalente de [[#eq-5|(5)]], en la cual se define <math display="inline">\lambda =\sqrt[4]{\dfrac{k}{4EI}}</math> y <math display="inline">v_i</math>, <math display="inline">\theta _i</math>, <math display="inline">v_j</math> y <math display="inline">\theta _j</math> tienen el mismo significado que para la viga (sección [[#5 Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga|5]]) | |
− | Para este caso la particularización de [[#eq-7|7]] es: | + | Para este caso la particularización de [[#eq-7|(7)]] es: |
<span id="eq-54"></span> | <span id="eq-54"></span> | ||
Line 1,059: | Line 1,083: | ||
|} | |} | ||
− | + | donde <math display="inline">v_h(x)</math> es el campo de desplazamiento homogéneo (sección [[#6.1 Solución del P.V.F. homogéneo|6.1]]), mientras que <math display="inline">v_f(x)</math> es el campo empotrado (sección [[#6.2 Solución del P.V.F. empotrado|6.2]]). | |
===6.1 Solución del P.V.F. homogéneo=== | ===6.1 Solución del P.V.F. homogéneo=== | ||
− | Para este caso el P.V.F. gobernante es: | + | Para este caso el P.V.F. gobernante es: <span id="eq-55"></span> |
<span id="eq-55.a"></span> | <span id="eq-55.a"></span> | ||
Line 1,092: | Line 1,116: | ||
|} | |} | ||
− | + | cuya solución se expresa como: | |
<span id="eq-56"></span> | <span id="eq-56"></span> | ||
Line 1,105: | Line 1,129: | ||
|} | |} | ||
− | + | donde <math display="inline">\Psi _2(x)</math>, <math display="inline">\Psi _3(x)</math>, <math display="inline">\Psi _5(x)</math> y <math display="inline">\Psi _6(x)</math> son las funciones de forma “exactas” de este problema y tienen el siguiente valor: | |
<span id="eq-57"></span> | <span id="eq-57"></span> | ||
Line 1,153: | Line 1,177: | ||
|} | |} | ||
− | Y para simplificar la escritura se ha definido | + | Y para simplificar la escritura, se ha definido <math display="inline">\hbox{s}=\sin (\lambda L)</math>, <math display="inline">\hbox{c}=\cos (\lambda L)</math>, <math display="inline">\hbox{sh}=\sinh (\lambda L)</math> y <math display="inline">\hbox{ch}=\cosh (\lambda L)</math>. |
− | Al igual que para el caso de la viga, las funciones de forma para la viga sobre fundación flexible también poseen las siguientes propiedades de simetría: | + | Al igual que para el caso de la viga, las funciones de forma para la viga sobre fundación flexible también poseen las siguientes propiedades de simetría:<span id="eq-61"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,176: | Line 1,200: | ||
|} | |} | ||
− | A diferencia de lo que ocurre para el caso la barra y la viga, las funciones de forma “exactas” presentadas en [[#eq-57|57]] a [[#eq-60|60]] 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 (ecuaciones [[# | + | A diferencia de lo que ocurre para el caso la barra y la viga, las funciones de forma “exactas” presentadas en [[#eq-57|(57)]] a [[#eq-60|(60)]] 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 (ecuaciones [[#eq-31|(31)]]). Como se presentará en el apéndice [[#C Formulación exacta del M.E.F. para un elemento tipo viga sobre fundación flexible|C]], si en la formulación del M.E.F. para una viga sobre fundación flexible se emplean las funciones de forma [[#eq-57|(57)]] a [[#eq-60|(60)]], se obtiene una formulación matricial “exacta” y una matriz de rigidez y fuerzas de empotramientos iguales a los obtenidos en el método presentado en este artículo (método de rigidez con funciones de Green), siempre y cuando se sumen tanto el campo de desplazamiento homogéneo como el empotrado. |
===6.2 Solución del P.V.F. empotrado=== | ===6.2 Solución del P.V.F. empotrado=== | ||
− | Para este caso el P.V.F. que define al problema empotrado es ( | + | Para este caso el P.V.F. que define al problema empotrado es ([[#img-8a|Figura 8a]]): <span id="eq-62"></span> |
<span id="eq-62.a"></span> | <span id="eq-62.a"></span> | ||
Line 1,212: | Line 1,236: | ||
<div id='img-8b'></div> | <div id='img-8b'></div> | ||
<div id='img-8'></div> | <div id='img-8'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura07a.png|360px|Carga distribuida arbitraria.]] |
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura07c.png|360px|Fuerza puntual unitaria.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Carga distribuida arbitraria | + | | (a) Carga distribuida arbitraria |
− | | (b) Fuerza puntual unitaria | + | | (b) Fuerza puntual unitaria |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 8 | + | | colspan="2" style="padding:10px;"| '''Figura 8'''. Vigas sobre fundación flexible de sección transversal constante doblemente empotradas, sometidas a una carga distribuida arbitraria y una carga puntual unitaria |
|} | |} | ||
− | Mientras que la función de Green se define por medio del siguiente P.V.F ( | + | |
+ | |||
+ | Mientras que la función de Green se define por medio del siguiente P.V.F ([[#img-8b|Figura 8b]]): <span id="eq-63"></span> | ||
<span id="eq-63.a"></span> | <span id="eq-63.a"></span> | ||
Line 1,259: | Line 1,285: | ||
|} | |} | ||
− | + | donde se ha usado como convención emplear una linea horizontal encima de la letra <math display="inline">G</math> para diferenciar esta función de Green de aquella de las vigas. | |
− | La solución de [[# | + | La solución de [[#eq-63|(63)]] se expresa en forma compacta empleando las funciones de forma [[#eq-57|(57)]] a [[#eq-60|(60)]], como: |
<span id="eq-64"></span> | <span id="eq-64"></span> | ||
Line 1,269: | Line 1,295: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\bar{G}_{yy}(x,\xi )= \begin{cases}\bar{G}_{yy}^I(x,\xi ) & 0 | + | | style="text-align: center;" | <math>\bar{G}_{yy}(x,\xi )= \begin{cases}\bar{G}_{yy}^I(x,\xi ) & 0 < x \leq \xi \\ \bar{G}_{yy}^{II}(x,\xi ) \qquad & \xi \leq x < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (64) | | style="width: 5px;text-align: right;white-space: nowrap;" | (64) | ||
|} | |} | ||
− | + | donde: | |
<span id="eq-65"></span> | <span id="eq-65"></span> | ||
Line 1,300: | Line 1,326: | ||
|} | |} | ||
− | La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en [[# | + | La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en [[#eq-61|(61)]] y las siguientes propiedades de simetría de la misma función de Green: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,315: | Line 1,341: | ||
|} | |} | ||
− | Multiplicando a ambos lados de [[#eq-63.a|63.a]] por <math display="inline">q(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: | + | Multiplicando a ambos lados de [[#eq-63.a|(63.a)]] por <math display="inline">q(\xi )</math> e integrando respecto a <math display="inline">\xi </math> entre 0 y <math display="inline">L</math> se obtiene: |
<span id="eq-68"></span> | <span id="eq-68"></span> | ||
Line 1,330: | Line 1,356: | ||
|} | |} | ||
− | + | donde se empleó que <math display="inline">q(x)=\int _0^L \delta (x-\xi ) q(\xi ) d\xi </math>. | |
− | De comparar [[#eq-21|21]] y [[#eq-62.a|62.a]] se concluye que: | + | De comparar [[#eq-21|(21)]] y [[#eq-62.a|(62.a)]] se concluye que: |
<span id="eq-69"></span> | <span id="eq-69"></span> | ||
Line 1,345: | Line 1,371: | ||
|} | |} | ||
− | A partir de [[#eq-69|69]] y [[# | + | A partir de [[#eq-69|(69)]] y [[#eq-4|(4)]], se obtiene que las fuerzas internas se pueden calcular como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,371: | Line 1,397: | ||
|} | |} | ||
− | Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales se calculan como: | + | Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales se calculan como: <span id="eq-72"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,392: | Line 1,418: | ||
|} | |} | ||
− | De forma similar a como se realizó para el caso de la viga, la función de Green debida a un momento puntual unitario se puede obtener directamente de aquella debida a una fuerza puntual mediante la derivada de la misma respecto a <math display="inline">\xi </math> ( | + | De forma similar a como se realizó para el caso de la viga, la función de Green debida a un momento puntual unitario se puede obtener directamente de aquella debida a una fuerza puntual mediante la derivada de la misma respecto a <math display="inline">\xi </math> ([[#img-9|Figura 9]]), es decir: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,404: | Line 1,430: | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,441: | Line 1,467: | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,478: | Line 1,504: | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,502: | Line 1,528: | ||
|} | |} | ||
− | Con lo cual, a partir de [[#eq-76|76]] y [[#eq-79|79]] se tiene que las fuerzas de empotramiento son: | + | Con lo cual, a partir de [[#eq-76|(76)]] y [[#eq-79|(79)]] se tiene que las fuerzas de empotramiento son: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,523: | Line 1,549: | ||
|} | |} | ||
− | Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sino de valor <math display="inline">M</math>, los resultados anteriores se deberán simplemente multiplicar por <math display="inline">M</math>. Mientras que si se trata de un momento distruido, debe realizarse la integración del campo debido al momento puntual multiplicado por la función que define al momento distribuido. <div id='img-9'></div> | + | Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sino de valor <math display="inline">M</math>, los resultados anteriores se deberán simplemente multiplicar por <math display="inline">M</math>. Mientras que si se trata de un momento distruido, debe realizarse la integración del campo debido al momento puntual multiplicado por la función que define al momento distribuido. |
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | |
+ | <div id='img-9'></div> | ||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;" | ||
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;"|[[Image:Molina-Villegas_et_al_2020b-Figura12.png|420px|Viga sobre funcdación flexible doblemente empotrada y de de sección transversal constante, sometidas a un momento puntual unitario.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 9 | + | | colspan="1" style="padding:10px;"| '''Figura 9.''' Viga sobre funcdación flexible doblemente empotrada y de de sección transversal constante, sometidas a un momento puntual unitario |
|} | |} | ||
===6.3 Superposición del campo homogéneo y el campo empotrado=== | ===6.3 Superposición del campo homogéneo y el campo empotrado=== | ||
− | Reemplazando [[#eq-56|56]] y [[#eq-69|69]] en [[#eq-54|54]], el campo de desplazamiento total se expresa en función de los desplazamientos de los extremos del elemento y de la carga externa como: | + | Reemplazando [[#eq-56|(56)]] y [[#eq-69|(69)]] en [[#eq-54|(54)]], el campo de desplazamiento total se expresa en función de los desplazamientos de los extremos del elemento y de la carga externa como: |
<span id="eq-83"></span> | <span id="eq-83"></span> | ||
Line 1,554: | Line 1,582: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\begin{Bmatrix}FY_i \\ M_i \\ FY_j \\ M_j \end{Bmatrix}= \begin{Bmatrix}FY_i^h \\ M_i^h \\ FY_j^h \\ M_j^h \end{Bmatrix}+ \begin{Bmatrix}FY_i^f \\ M_i^f \\ FY_j^f \\ M_j^f \end{Bmatrix}= \begin{bmatrix}k_{22} & k_{23} & k_{25} & k_{26} \\ k_{32} & k_{33} & k_{35} & k_{36} \\ k_{52} & k_{53} & k_{55} & k_{56} \\ k_{62} & k_{63} & k_{65} & k_{66} \end{bmatrix} \begin{Bmatrix}v_i \\ \theta _i \\ v_j \\ \theta _j \end{Bmatrix}- \begin{Bmatrix}\int _0^L \Psi _2(x) q(x)dx \\ \int _0^L \Psi _3(x) q(x)dx \\ \int _0^L \Psi _5(x) q(x)dx \\ \int _0^L \Psi _6(x) q(x)dx \end{Bmatrix} </math> | + | | style="text-align: center;" | <math>\begin{Bmatrix}FY_i \\ M_i \\ FY_j \\ M_j \end{Bmatrix}= \begin{Bmatrix}FY_i^h \\ M_i^h \\ FY_j^h \\ M_j^h \end{Bmatrix}+ \begin{Bmatrix}FY_i^f \\ M_i^f \\ FY_j^f \\ M_j^f \end{Bmatrix}= \begin{bmatrix}k_{22} & k_{23} & k_{25} & k_{26} \\ k_{32} & k_{33} & k_{35} & k_{36} \\ k_{52} & k_{53} & k_{55} & k_{56} \\ k_{62} & k_{63} & k_{65} & k_{66} \end{bmatrix} \begin{Bmatrix}v_i \\ \theta _i \\ v_j \\ \theta _j \end{Bmatrix}- \begin{Bmatrix}\displaystyle\int _0^L \Psi _2(x) q(x)dx \\ \displaystyle\int _0^L \Psi _3(x) q(x)dx \\ \displaystyle\int _0^L \Psi _5(x) q(x)dx \\ \displaystyle\int _0^L \Psi _6(x) q(x)dx \end{Bmatrix} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (84) | | style="width: 5px;text-align: right;white-space: nowrap;" | (84) | ||
|} | |} | ||
− | + | donde los términos de la matriz de rigidez son (<span id='citeF-21'></span>[[#cite-21|[21]]]): <span id="eq-85"></span> | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,586: | Line 1,614: | ||
|} | |} | ||
− | Es importante resaltar que [[#eq-84|84]] corresponde a la formulación analítica del método de rigidez en coordenadas locales 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 [[# | + | Es importante resaltar que [[#eq-84|(84)]] corresponde a la formulación analítica del método de rigidez en coordenadas locales 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 [[#C Formulación exacta del M.E.F. para un elemento tipo viga sobre fundación flexible|C]]. |
− | ==7 Formulación analítica del método de rigidez con funciones de Green para elementos tipo pórtico plano== | + | ==7. Formulación analítica del método de rigidez con funciones de Green para elementos tipo pórtico plano== |
− | Como se mencionó en la sección [[#2.4 Elemento tipo pórtico plano|2.4]] 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 [[#eq-26|26]] y [[#eq-52|52]] se tiene que la formulación del método de rigidez con | + | Como se mencionó en la sección [[#2.4 Elemento tipo pórtico plano|2.4]] 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 [[#eq-26|(26)]] y [[#eq-52|(52)]] se tiene que la formulación del método de rigidez con funciones de Green en coordenadas locales para un elemento tipo pórtico plano es: |
<span id="eq-86"></span> | <span id="eq-86"></span> | ||
Line 1,598: | Line 1,626: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\begin{Bmatrix}FX_i \\ FY_i \\ M_i \\ FX_j \\ FY_j \\ M_j \end{Bmatrix}= \begin{bmatrix}\dfrac{AE}{L} & 0 & 0 & -\dfrac{AE}{L} & 0 & 0 \\ 0 & \dfrac{12EI}{L^3} & \dfrac{6EI}{L^2} & 0 &-\dfrac{12EI}{L^3} & \dfrac{6EI}{L^2} \\ 0 & \dfrac{6EI}{L^2} & \dfrac{4EI}{L} & 0 & -\dfrac{6EI}{L^2} & \dfrac{2EI}{L} \\ -\dfrac{AE}{L} & 0 & 0 & \dfrac{AE}{L} & 0 & 0 \\ 0 & -\dfrac{12EI}{L^3} & -\dfrac{6EI}{L^2} & 0 & \dfrac{12EI}{L^3} & -\dfrac{6EI}{L^2} \\ 0 & \dfrac{6EI}{L^2} & \dfrac{2EI}{L} & 0 & -\dfrac{6EI}{L^2} & \dfrac{4EI}{L} \end{bmatrix} \begin{Bmatrix}u_i \\ v_i \\ \theta _i \\ u_j \\ v_j \\ \theta _j \end{Bmatrix}- \begin{ | + | | style="text-align: center;" | <math>\begin{Bmatrix}FX_i \\ FY_i \\ M_i \\ FX_j \\ FY_j \\ M_j \end{Bmatrix}= \begin{bmatrix}\dfrac{AE}{L} & 0 & 0 & -\dfrac{AE}{L} & 0 & 0 \\ 0 & \dfrac{12EI}{L^3} & \dfrac{6EI}{L^2} & 0 &-\dfrac{12EI}{L^3} & \dfrac{6EI}{L^2} \\ 0 & \dfrac{6EI}{L^2} & \dfrac{4EI}{L} & 0 & -\dfrac{6EI}{L^2} & \dfrac{2EI}{L} \\ -\dfrac{AE}{L} & 0 & 0 & \dfrac{AE}{L} & 0 & 0 \\ 0 & -\dfrac{12EI}{L^3} & -\dfrac{6EI}{L^2} & 0 & \dfrac{12EI}{L^3} & -\dfrac{6EI}{L^2} \\ 0 & \dfrac{6EI}{L^2} & \dfrac{2EI}{L} & 0 & -\dfrac{6EI}{L^2} & \dfrac{4EI}{L} \end{bmatrix} \begin{Bmatrix}u_i \\ v_i \\ \theta _i \\ u_j \\ v_j \\ \theta _j \end{Bmatrix}- \left\{\begin{matrix}\displaystyle\int _0^L \psi _1(x) p(x)dx \\ \displaystyle\int _0^L \psi _2(x) q(x)dx \\ \displaystyle\int _0^L \psi _3(x) q(x)dx \\ \displaystyle\int _0^L \psi _4(x) p(x)dx \\ \displaystyle\int _0^L \psi _5(x) q(x)dx \\ \displaystyle\int _0^L \psi _6(x) q(x)dx \end{matrix}\right\} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (86) | | style="width: 5px;text-align: right;white-space: nowrap;" | (86) | ||
|} | |} | ||
− | Mientras que los campos de desplazamiento en dirección de los ejes locales <math display="inline">x</math> y <math display="inline">y</math>, se calculan a partir de las ecuaciones [[#eq-25|25]] y [[#eq-51|51]] respectivamente. | + | Mientras que los campos de desplazamiento en dirección de los ejes locales <math display="inline">x</math> y <math display="inline">y</math>, se calculan a partir de las ecuaciones [[#eq-25|(25)]] y [[#eq-51|(51)]] respectivamente. |
− | ==8 Ejemplos== | + | ==8. Ejemplos== |
===8.1 Viga=== | ===8.1 Viga=== | ||
− | Calcular la respuesta (reacciones, campos de desplazamientos y campos de fuerzas internas) de la viga presentada en la | + | Calcular la respuesta (reacciones, campos de desplazamientos y campos de fuerzas internas) de la viga presentada en la [[#img-10a|Figura 10a]]. |
+ | |||
+ | <div id='img-10a'></div> | ||
<div id='img-10b'></div> | <div id='img-10b'></div> | ||
<div id='img-10'></div> | <div id='img-10'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 85%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura08a.png|390px|Problema a resolver.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura08b.png|390px|Discretización.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Problema a resolver | + | | (a) Problema a resolver |
− | | (b) Discretización | + | | (b) Discretización |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 10 | + | | colspan="2" style="padding:10px;" | '''Figura 10'''. Viga con articulación interior y discretización empleada para su solución |
|} | |} | ||
+ | |||
+ | |||
<u>Discretización</u> | <u>Discretización</u> | ||
− | La discretización a emplear en la solución de este problema se presenta en la | + | La discretización a emplear en la solución de este problema se presenta en la [[#img-10b|Figura 10b]] y consta de solo dos elementos. El elemento <math display="inline">A</math> tiene una carga externa definida en tres tramos, mientras que el elemento <math display="inline">B</math> tiene una carga externa lineal sobre toda su longitud. |
<u>Definición de las funciones de carga externa</u> | <u>Definición de las funciones de carga externa</u> | ||
Line 1,647: | Line 1,679: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{q_A(x'}_A)= \begin{cases}q_A^{I}{(x'} A) \qquad & | + | | style="text-align: center;" | <math>{q_A(x'}_A)= \begin{cases}q_A^{I}{(x'} A) \qquad &{ 0 < x'}_A < \dfrac{L}{3} \\[0.3cm] q_A^{II}{(x'} A) & \dfrac{L}{3}{ < x'}_A < \dfrac{2L}{3} \\[0.3cm] q_A^{III}{(x'} A) & \dfrac{2L}{3}{ < x'}_A < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (88) | | style="width: 5px;text-align: right;white-space: nowrap;" | (88) | ||
Line 1,659: | Line 1,691: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{q_B(x'}_B)=Q\left(-2+2\frac{{x'}_B}{L} \right) </math> | + | | style="text-align: center;" | <math>{q_B(x'}_B)=Q\left(-2+2\frac{{x'}_B}{L} \right)\qquad {0<x'}_B<L </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (89) | | style="width: 5px;text-align: right;white-space: nowrap;" | (89) | ||
Line 1,666: | Line 1,698: | ||
<u>Cálculo de las fuerzas de empotramiento</u> | <u>Cálculo de las fuerzas de empotramiento</u> | ||
− | Según lo presentado en [[# | + | Según lo presentado en [[#eq-45|(45)]], las fuerzas de empotramiento del elemento <math display="inline">A</math> se calculan como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,673: | Line 1,705: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>FY_1^{Af}=-\int _0^L \psi | + | | style="text-align: center;" | <math>FY_1^{Af}=-\int _0^L \psi _2(x'_A)q_A(x'_A)dx'_A=\dfrac{40}{81}QL </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (90.a) | | style="width: 5px;text-align: right;white-space: nowrap;" | (90.a) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> M_1 ^{Af}=-\int _0^L \psi | + | | style="text-align: center;" | <math> M_1 ^{Af}=-\int _0^L \psi _3(x'_A)q_A(x'_A)dx'_A=\dfrac{71}{1215}QL^2 </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (90.b) | | style="width: 5px;text-align: right;white-space: nowrap;" | (90.b) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> FY_2^{Af}=-\int _0^L \psi | + | | style="text-align: center;" | <math> FY_2^{Af}=-\int _0^L \psi _5(x'_A)q_A(x'_A)dx'_A=\dfrac{40}{81}QL </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (90.c) | | style="width: 5px;text-align: right;white-space: nowrap;" | (90.c) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> M_2 ^{Af}=-\int _0^L \psi | + | | style="text-align: center;" | <math> M_2 ^{Af}=-\int _0^L \psi _6(x'_A)q_A(x'_A)dx'_A=-\dfrac{71}{1215}QL^2 </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (90.d) | | style="width: 5px;text-align: right;white-space: nowrap;" | (90.d) | ||
|} | |} | ||
Line 1,694: | Line 1,726: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>FY_2^{Bf}=-\int _0^L \psi | + | | style="text-align: center;" | <math>FY_2^{Bf}=-\int _0^L \psi _2(x'_B)q_B(x'_B)dx'_B=\dfrac{7}{10}QL </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (91.a) | | style="width: 5px;text-align: right;white-space: nowrap;" | (91.a) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> M_2 ^{Bf}=-\int _0^L \psi | + | | style="text-align: center;" | <math> M_2 ^{Bf}=-\int _0^L \psi _3(x'_B)q_B(x'_B)dx'_B=\dfrac{1}{10}QL^2 </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (91.b) | | style="width: 5px;text-align: right;white-space: nowrap;" | (91.b) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> FY_3^{Bf}=-\int _0^L \psi | + | | style="text-align: center;" | <math> FY_3^{Bf}=-\int _0^L \psi _5(x'_B)q_B(x'_B)dx'_B=\dfrac{3}{10}QL </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (91.c) | | style="width: 5px;text-align: right;white-space: nowrap;" | (91.c) | ||
|- | |- | ||
− | | style="text-align: center;" | <math> M_3 ^{Bf}=-\int _0^L \psi | + | | style="text-align: center;" | <math> M_3 ^{Bf}=-\int _0^L \psi _6(x'_B)q_B(x'_B)dx'_B=-\dfrac{1}{15}QL^2 </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (91.d) | | style="width: 5px;text-align: right;white-space: nowrap;" | (91.d) | ||
|} | |} | ||
Line 1,710: | Line 1,742: | ||
<u>Sistema matricial de ecuaciones de cada elemento</u> | <u>Sistema matricial de ecuaciones de cada elemento</u> | ||
− | A partir de lo presentado en [[#eq-52|52]], la formulación del método de rigidez para el elemento <math display="inline">A</math> es: | + | A partir de lo presentado en [[#eq-52|(52)]], la formulación del método de rigidez para el elemento <math display="inline">A</math> es: |
<span id="eq-92"></span> | <span id="eq-92"></span> | ||
Line 1,764: | Line 1,796: | ||
<u>Cálculo del campo de desplazamiento del elemento <math>A</math></u> | <u>Cálculo del campo de desplazamiento del elemento <math>A</math></u> | ||
− | Reemplazando los desplazamientos nodales del elemento <math display="inline">A</math> en [[#eq-30|30]], se obtiene que el campo de desplazamiento homogéneo para este elemento es: | + | Reemplazando los desplazamientos nodales del elemento <math display="inline">A</math> en [[#eq-30|(30)]], se obtiene que el campo de desplazamiento homogéneo para este elemento es: |
<span id="eq-96"></span> | <span id="eq-96"></span> | ||
Line 1,785: | Line 1,817: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{ v_A^f(x'}_A)= \begin{cases}v_A^{f-1}{(x'} A) \qquad & | + | | style="text-align: center;" | <math>{ v_A^f(x'}_A)= \begin{cases}v_A^{f-1}{(x'} A) \qquad &{ 0 < x'}_A \leq \dfrac{L}{3} \\[0.3cm] v_A^{f-2}{(x'} A) \qquad & \dfrac{L}{3} \leq {x'}_A \leq \dfrac{2L}{3} \\[0.3cm] v_A^{f-3}{(x'} A) & \dfrac{2L}{3} \leq {x'} < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (97) | | style="width: 5px;text-align: right;white-space: nowrap;" | (97) | ||
|} | |} | ||
− | + | donde, a partir de [[#eq-42|(42)]], el campo de desplazamiento empotrado en cada uno de los tres tramos se calcula como: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,828: | Line 1,860: | ||
|} | |} | ||
− | Ahora sumando [[#eq-96|96]] y [[#eq-97|97]] se obtiene que el campo de desplazamiento total en el elemento es: | + | Ahora sumando [[#eq-96|(96)]] y [[#eq-97|(97)]] se obtiene que el campo de desplazamiento total en el elemento es: |
<span id="eq-101"></span> | <span id="eq-101"></span> | ||
Line 1,836: | Line 1,868: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{ v_A(x'}_A)= \begin{cases}{v_A^I(x'}_A) & 0 | + | | style="text-align: center;" | <math>{ v_A(x'}_A)= \begin{cases}{v_A^I(x'}_A) &{ 0 < x'}_A \leq \dfrac{L}{3}\\[0.5cm] v_A^{II}{(x'} A) & \dfrac{L}{3} \leq {x'}_A \leq \dfrac{2}{3}L \\[0.5cm] v_A^{III}{(x'} A) \qquad & \dfrac{2}{3}L \leq {x'}_< L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (101) | | style="width: 5px;text-align: right;white-space: nowrap;" | (101) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,903: | Line 1,935: | ||
|} | |} | ||
− | Como resumen, en la | + | Como resumen, en la [[#img-11|Figura 11]] se presenta el campo de desplazamientos en toda la viga. |
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | |
+ | <div id='img-11'></div> | ||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;" | ||
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo1CampoDesplazamientoDefinitivo.png|450px|Diagrama del campo de desplazamientos.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" | '''Figura 11 | + | | colspan="1" style="padding:10px;" | '''Figura 11'''. Diagrama del campo de desplazamientos |
|} | |} | ||
+ | |||
+ | |||
<u>Cálculos de los campos de fuerzas internas</u> | <u>Cálculos de los campos de fuerzas internas</u> | ||
− | A partir de [[#eq-101|101]] y empleando [[# | + | A partir de [[#eq-101|(101)]] y empleando [[#eq-4|(4)]], se obtiene que las fuerzas internas para el elemento <math display="inline">A</math> son: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,919: | Line 1,955: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{V_A(x'}_A)= \begin{cases}QL \left[-\dfrac{3433}{3240} + 2\dfrac{{x'}_A}{L} - 2\left(\dfrac{{x'}_A}{L} \right)^2+\dfrac{4}{3}\left(\dfrac{{x'}_A}{L} \right)^3 \right]\qquad & 0 | + | | style="text-align: center;" | <math>{V_A(x'}_A)= \begin{cases}QL \left[-\dfrac{3433}{3240} + 2\dfrac{{x'}_A}{L} - 2\left(\dfrac{{x'}_A}{L} \right)^2+\dfrac{4}{3}\left(\dfrac{{x'}_A}{L} \right)^3 \right]\qquad &{ 0 < x'}_A \leq \dfrac{L}{3}\\[0.5cm] -\dfrac{611}{1080}QL & \dfrac{L}{3} \leq {x'}_A \leq \dfrac{2}{3}L \\[0.5cm] QL\left[-\dfrac{4553}{3240}+2\dfrac{{x'}_A}{L}-2\left(\dfrac{{x'}_A}{L} \right)^2+\dfrac{4}{3} \left(\dfrac{{x'}_A}{L} \right)^3 \right]& \dfrac{2}{3}L \leq {x'}_< L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (107) | | style="width: 5px;text-align: right;white-space: nowrap;" | (107) | ||
Line 1,934: | Line 1,970: | ||
|} | |} | ||
− | De forma similar, a partir de [[#eq-106|106]], se tiene que los campos de fuerzas internas en el elemento <math display="inline">B</math> son: | + | De forma similar, a partir de [[#eq-106|(106)]], se tiene que los campos de fuerzas internas en el elemento <math display="inline">B</math> son: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 1,941: | Line 1,977: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{V_B(x'}_B)=QL \left[-\dfrac{233}{3240} + 2\dfrac{{x'}_B}{L} - \left(\dfrac{{x'}_B}{L} \right)^2 \right] | + | | style="text-align: center;" | <math>{V_B(x'}_B)=QL \left[-\dfrac{233}{3240} + 2\dfrac{{x'}_B}{L} - \left(\dfrac{{x'}_B}{L} \right)^2 \right]{ 0 < x'}_B < L </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (109.a) | | style="width: 5px;text-align: right;white-space: nowrap;" | (109.a) | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{ M_B(x'}_B)=QL^2 \left[\frac{233}{3240}{x'}_B-\left(\frac{{x'}_B}{L} \right)^2 + \frac{1}{3}\left(\frac{{x'}_B}{L} \right)^3 \right] | + | | style="text-align: center;" | <math>{ M_B(x'}_B)=QL^2 \left[\frac{233}{3240}{x'}_B-\left(\frac{{x'}_B}{L} \right)^2 + \frac{1}{3}\left(\frac{{x'}_B}{L} \right)^3 \right]{ 0 < x'}_B < L </math> |
| style="width: 5px;text-align: right;white-space: nowrap;" | (109.b) | | style="width: 5px;text-align: right;white-space: nowrap;" | (109.b) | ||
|} | |} | ||
|} | |} | ||
− | Como resumen de los anteriores resultados, en la | + | Como resumen de los anteriores resultados, en la [[#img-12|Figura 12]] se presentan las fuerzas internas en toda la viga. |
+ | |||
+ | <div id='img-12a'></div> | ||
<div id='img-12b'></div> | <div id='img-12b'></div> | ||
<div id='img-12'></div> | <div id='img-12'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 90%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo1CortanteDefinitivo.png|480px|Fuerza cortante.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo1MomentoDefinitivo.png|480px|Momento flector.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Fuerza cortante | + | | (a) Fuerza cortante |
− | | (b) Momento flector | + | | (b) Momento flector |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 12 | + | | colspan="2" style="padding:10px;" | '''Figura 12'''. Diagramas de fuerzas internas |
|} | |} | ||
+ | |||
+ | |||
<u>Cálculo de las reacciones</u> | <u>Cálculo de las reacciones</u> | ||
− | Aunque es posible calcular las reacciones mediante | + | Aunque es posible calcular las reacciones mediante 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: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,004: | Line 2,044: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\sum M_1=M_1+M_3+2FY_3+\int _0^{L/3} Q\left[-2+4\frac{{x'}_A}{L}-4\left(\frac{{x'}_A}{L} \right)^2 \right]x'_A dx'_A+\int _{2L/3}^{L}Q\left[-2+4\frac{{x'}_A}{L}-4\left(\frac{{x'}_A}{L} \right)^2 \right]x'_A dx'_A </math> | + | | style="text-align: center;" | <math>\sum M_1=M_1+M_3+2FY_3+\int _0^{L/3} Q\left[-2+4\frac{{x'}_A}{L}-4\left(\frac{{x'}_A}{L} \right)^2 \right]{x'}_A {dx'}_A +\int _{2L/3}^{L}Q\left[-2+4\frac{{x'}_A}{L}-4\left(\frac{{x'}_A}{L} \right)^2 \right]{x'}_A {dx'}_A </math> |
|- | |- | ||
− | | style="text-align: center;" | <math> +\int _0^L Q \left(-2+2\frac{{x'}_B}{L} \right)(L+x'_B)dx'_B=0 </math> | + | | style="text-align: center;" | <math> +\int _0^L Q \left(-2+2\frac{{x'}_B}{L} \right){(L+x'}_B)dx'_B=0 </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (112) | | style="width: 5px;text-align: right;white-space: nowrap;" | (112) | ||
Line 2,017: | Line 2,057: | ||
===8.2 Viga sobre fundación flexible=== | ===8.2 Viga sobre fundación flexible=== | ||
− | Resolver la viga sobre fundación flexible presentada en la | + | Resolver la viga sobre fundación flexible presentada en la [[#img-13a|Figura 13a]]. |
+ | |||
+ | <div id='img-13a'></div> | ||
<div id='img-13b'></div> | <div id='img-13b'></div> | ||
<div id='img-13'></div> | <div id='img-13'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 85%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura09a.png|390px|Problema a resolver.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura09b.png|390px|Discretización.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Problema a resolver | + | | (a) Problema a resolver |
− | | (b) Discretización | + | | (b) Discretización |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 13 | + | | colspan="2" style="padding:10px;" | '''Figura 13'''. Viga sobre fundación flexible y discretización empleada para su solución |
|} | |} | ||
+ | |||
+ | |||
<u>Discretización</u> | <u>Discretización</u> | ||
− | Para este problema se empleará la discretización de un solo elemento presentada en la | + | Para este problema se empleará la discretización de un solo elemento presentada en la [[#img-13b|Figura 13b]], la cual muestra una de las principales bondades de la presente metodología. Consistente 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 analítica. |
<u>Preliminares</u> | <u>Preliminares</u> | ||
Line 2,043: | Line 2,087: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\lambda =\sqrt | + | | style="text-align: center;" | <math>\lambda =\sqrt[4]{\dfrac{k}{4EI}}=1.535260 </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (113) | | style="width: 5px;text-align: right;white-space: nowrap;" | (113) | ||
Line 2,100: | Line 2,144: | ||
|} | |} | ||
− | Por su parte la función de Green de la viga doblemente empotrada asociada con este problema se obtiene a partir de [[#eq-65|65]] y [[#eq-66|66]]: | + | Por su parte la función de Green de la viga doblemente empotrada asociada con este problema se obtiene a partir de [[#eq-65|(65)]] y [[#eq-66|(66)]]: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,107: | Line 2,151: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\bar{G}_{yy}(x,\xi )= \begin{cases}\bar{G}_{yy}^I (x,\xi ) & 0 | + | | style="text-align: center;" | <math>\bar{G}_{yy}(x,\xi )= \begin{cases}\bar{G}_{yy}^I (x,\xi ) & 0 < x \leq \xi \\ \bar{G}_{yy}^{II}(x,\xi ) & \xi \leq x < 10\hbox{ m} \end{cases}= \begin{cases}C(x)\Psi _2(\xi )+D(x)\Psi _3(\xi ) & 0 < x \leq \xi \\ A(x)\Psi _5(\xi )+B(x)\Psi _6(\xi ) & \xi \leq x < 10\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (118) | | style="width: 5px;text-align: right;white-space: nowrap;" | (118) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,147: | Line 2,191: | ||
|} | |} | ||
− | Con lo cual las fuerzas de empotramiento se calculan a partir de las funciones de forma presentadas en [[#eq-114|114]] a [[#eq-117|117]], como: | + | Con lo cual las fuerzas de empotramiento se calculan a partir de las funciones de forma presentadas en [[#eq-114|(114)]] a [[#eq-117|(117)]], como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,170: | Line 2,214: | ||
<u>Formulación del sistema de ecuaciones de cada elemento</u> | <u>Formulación del sistema de ecuaciones de cada elemento</u> | ||
− | A partir de [[#eq-84|84]] se tiene que la formulación del método de rigidez para el único elemento de la discretización es: | + | A partir de [[#eq-84|(84)]] se tiene que la formulación del método de rigidez para el único elemento de la discretización es: |
<span id="eq-122"></span> | <span id="eq-122"></span> | ||
Line 2,185: | Line 2,229: | ||
<u>Obtención de los desplazamientos nodales</u> | <u>Obtención de los desplazamientos nodales</u> | ||
− | La solución de [[#eq-122|122]] es: | + | La solución de [[#eq-122|(122)]] es: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,211: | Line 2,255: | ||
|} | |} | ||
− | + | donde <math display="inline">v_h(x)</math> es el campo de desplazamientos homogéneo y <math display="inline">v_f(x)</math> es el campo de desplazamientos empotrado. A partir de [[#eq-56|(56)]] y empleando los desplazamientos y rotaciones nodales, se obtiene que el campo de desplazamiento homogéneo es: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,225: | Line 2,269: | ||
|} | |} | ||
− | Mientras que el campo de desplazamiento empotrado se calcula a partir de [[#eq-69|69]] y se expresa como: | + | Mientras que el campo de desplazamiento empotrado se calcula a partir de [[#eq-69|(69)]] y se expresa como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,232: | Line 2,276: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>v_f(x)= \begin{cases}v_f^I(x) & 0 | + | | style="text-align: center;" | <math>v_f(x)= \begin{cases}v_f^I(x) & 0 < x \leq 3\hbox{ m} \\ v_f^{II}(x) & 3\hbox{ m} \leq x \leq 4\hbox{ m} \\ v_f^{III}(x) & 4\hbox{ m} \leq x < 5\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (126) | | style="width: 5px;text-align: right;white-space: nowrap;" | (126) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,246: | Line 2,290: | ||
| style="text-align: center;" | <math>v_f^I(x)=\int _3^4 \bar{G}_{yy}^I(x,\xi )q(\xi )d\xi </math> | | style="text-align: center;" | <math>v_f^I(x)=\int _3^4 \bar{G}_{yy}^I(x,\xi )q(\xi )d\xi </math> | ||
|- | |- | ||
− | | style="text-align: center;" | <math> =1.297541x10^{-8}\sin (\lambda x)\sinh (\lambda x) -3.344800x10^{-9}\sin (\lambda x)\cosh (\lambda x) | + | | style="text-align: center;" | <math> =1.297541x10^{-8}\sin (\lambda x)\sinh (\lambda x) -3.344800x10^{-9}\sin (\lambda x)\cosh (\lambda x)+3.344800x10^{-9}\cos (\lambda x)\sinh (\lambda x) </math> |
− | + | ||
− | + | ||
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (127) | | style="width: 5px;text-align: right;white-space: nowrap;" | (127) | ||
Line 2,288: | Line 2,330: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>v(x)= \begin{cases}v^I(x) & 0 | + | | style="text-align: center;" | <math>v(x)= \begin{cases}v^I(x) & 0 < x \leq 3\hbox{ m} \\ v^{II}(x) & 3\hbox{ m} \leq x \leq 4\hbox{ m} \\ v^{III}(x) & 4\hbox{ m} \leq x < 5\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (130) | | style="width: 5px;text-align: right;white-space: nowrap;" | (130) | ||
Line 2,333: | Line 2,375: | ||
|} | |} | ||
− | En la | + | En la [[#img-11|Figura 11]] se presenta de forma gráfica el campo de desplazamiento vertical en toda la viga. |
<u>Cálculo de las fuerzas que el suelo le hace a la viga</u> | <u>Cálculo de las fuerzas que el suelo le hace a la viga</u> | ||
Line 2,344: | Line 2,386: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>f_S(x)=-kv(x)= \begin{cases}f_S^I(x) & 0 | + | | style="text-align: center;" | <math>f_S(x)=-kv(x)= \begin{cases}f_S^I(x) & 0 < x \leq 3\hbox{ m} \\ f_S^{II}(x) & 3\hbox{ m} \leq x \leq 4\hbox{ m} \\ f_S^{III}(x) & 4\hbox{ m} \leq x < 5\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (134) | | style="width: 5px;text-align: right;white-space: nowrap;" | (134) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,393: | Line 2,435: | ||
|} | |} | ||
− | Como resumen de las anteriores tres ecuaciones, en la | + | Como resumen de las anteriores tres ecuaciones, en la [[#img-14b|Figura 14b]] se presenta la fuerza el suelo ejerce sobre la viga. |
+ | |||
+ | <div id='img-14a'></div> | ||
<div id='img-14b'></div> | <div id='img-14b'></div> | ||
<div id='img-14'></div> | <div id='img-14'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 90%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo2CampoDesplazamientoDefinitiva.png|450px|Campo de desplazamientos.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo2CampoFuerzaSueloVigaDefinitiva.png|450px|Campo de la fuerza que el suelo ejerce sobre la viga.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Campo de desplazamientos | + | | (a) Campo de desplazamientos |
− | | (b) Campo de la fuerza que el suelo ejerce sobre la viga | + | | (b) Campo de la fuerza que el suelo ejerce sobre la viga |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 14 | + | | colspan="2" style="padding:10px;" | '''Figura 14'''. Campos de desplazamiento en la viga y fuerza que el suelo ejerce sobre esta |
|} | |} | ||
+ | |||
+ | |||
Como revisión de 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: | Como revisión de 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: | ||
Line 2,430: | Line 2,476: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>M(x)=EI \frac{d^2 v}{dx^2}(x)= \begin{cases}M^I(x) & 0 | + | | style="text-align: center;" | <math>M(x)=EI \frac{d^2 v}{dx^2}(x)= \begin{cases}M^I(x) & 0 < x \leq 3\hbox{ m} \\ M^{II}(x) & 3\hbox{ m} \leq x \leq 4\hbox{ m} \\ M^{III}(x) & 4\hbox{ m} \leq x < 5\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (139) | | style="width: 5px;text-align: right;white-space: nowrap;" | (139) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,486: | Line 2,532: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>V(x)=-EI \frac{d^3 v}{dx^3}(x)= \begin{cases}V^I(x) & 0 | + | | style="text-align: center;" | <math>V(x)=-EI \frac{d^3 v}{dx^3}(x)= \begin{cases}V^I(x) & 0 < x \leq 3\hbox{ m} \\ V^{II}(x) & 3\hbox{ m} \leq x \leq 4\hbox{ m} \\ V^{III}(x) & 4\hbox{ m} \leq x < 5\hbox{ m} \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (143) | | style="width: 5px;text-align: right;white-space: nowrap;" | (143) | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,535: | Line 2,581: | ||
|} | |} | ||
− | En las | + | En las [[#img-12a|Figuras 12a]] y [[#img-12b|12b]] se presentan respectivamente los campos de fuerza cortante y momento flector para toda la viga. |
+ | |||
+ | <div id='img-15a'></div> | ||
<div id='img-15b'></div> | <div id='img-15b'></div> | ||
<div id='img-15'></div> | <div id='img-15'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 90%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo2CampoCortanteDefinitiva.png|450px|Fuerza cortante.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo2CampoMomentoFlectorDefinitiva.png|450px|Momento flector.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Fuerza cortante | + | | (a) Fuerza cortante |
− | | (b) Momento flector | + | | (b) Momento flector |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 15 | + | | colspan="2" style="padding:10px;" | '''Figura 15'''. Campos de fuerzas internas |
|} | |} | ||
En el siguiente enlace se encuentra el código de Python empleado para la solución de este ejemplo: | En el siguiente enlace se encuentra el código de Python empleado para la solución de este ejemplo: | ||
Line 2,554: | Line 2,602: | ||
===8.3 Pórtico plano=== | ===8.3 Pórtico plano=== | ||
− | Resolver el pórtico plano presentado en la | + | Resolver el pórtico plano presentado en la [[#img-16a|Figura 16a]] cuyos elementos son rectangulares de base y altura iguales a <math display="inline">L/20</math> y cuyo módulo de elasticidad es <math display="inline">E</math>. |
+ | |||
+ | <div id='img-16a'></div> | ||
<div id='img-16b'></div> | <div id='img-16b'></div> | ||
<div id='img-16'></div> | <div id='img-16'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 90%;" |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura10a.png|390px|Problema a resolver.]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-Figura10b.png|390px|Discretización.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | (a) Problema a resolver | + | | (a) Problema a resolver |
− | | (b) Discretización | + | | (b) Discretización |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura 16 | + | | colspan="2" style="padding:10px;" | '''Figura 16'''. Pórtico plano sometido a cargas externas oblicuas y discretización empleada para su solución |
|} | |} | ||
+ | |||
+ | |||
<u>Discretización</u> | <u>Discretización</u> | ||
− | En la | + | En la [[#img-16b|Figura 16b]] se presenta la discretización y ejes locales a emplear en la solución de este ejercicio. |
<u>Definición de la carga externa en coordenadas locales</u> | <u>Definición de la carga externa en coordenadas locales</u> | ||
− | A partir de la | + | A partir de la [[#img-16a|Figura 16a]], es posible obtener el valor de la carga externa en dirección <math display="inline">{x'}_A</math> del elemento <math display="inline">A</math>: |
<span id="eq-147"></span> | <span id="eq-147"></span> | ||
Line 2,581: | Line 2,633: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{ p_A(x'}_A)= \begin{cases}\dfrac{24}{25}Q\left(-\dfrac{1}{2}+\dfrac{{x'}_A}{L}\right) \qquad & | + | | style="text-align: center;" | <math>{ p_A(x'}_A)= \begin{cases}\dfrac{24}{25}Q\left(-\dfrac{1}{2}+\dfrac{{x'}_A}{L}\right) \qquad &{ 0 < x'}_A < \dfrac{L}{2} \\[0.3cm] \dfrac{24}{25}Q\left(-1+\dfrac{{x'}_A}{L}\right) & \dfrac{L}{2}{ < x'}_A < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (147) | | style="width: 5px;text-align: right;white-space: nowrap;" | (147) | ||
Line 2,594: | Line 2,646: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{ q_A(x'}_A)= \begin{cases}\dfrac{32}{25}Q\left(-\dfrac{1}{2}+\dfrac{{x'}_A}{L}\right)\qquad & | + | | style="text-align: center;" | <math>{ q_A(x'}_A)= \begin{cases}\dfrac{32}{25}Q\left(-\dfrac{1}{2}+\dfrac{{x'}_A}{L}\right)\qquad &{ 0 < x'}_A < \dfrac{L}{2} \\[0.3cm] \dfrac{32}{25}Q\left(-1+\dfrac{{x'}_A}{L}\right) & \dfrac{L}{2}{ < x'}_A < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (148) | | style="width: 5px;text-align: right;white-space: nowrap;" | (148) | ||
Line 2,616: | Line 2,668: | ||
<u>Cálculo del sistema de ecuaciones de cada elemento en coordenadas locales</u> | <u>Cálculo del sistema de ecuaciones de cada elemento en coordenadas locales</u> | ||
− | A partir de lo presentado en [[#eq-86|86]] se tiene que las ecuaciones del método de rigidez en coordenadas locales para los elementos <math display="inline">A</math> y <math display="inline">B</math> son respectivamente: | + | A partir de lo presentado en [[#eq-86|(86)]] se tiene que las ecuaciones del método de rigidez en coordenadas locales para los elementos <math display="inline">A</math> y <math display="inline">B</math> son respectivamente: |
<span id="eq-150"></span> | <span id="eq-150"></span> | ||
Line 2,642: | Line 2,694: | ||
<u>Cálculo del sistema de ecuaciones de cada elemento en coordenadas globales</u> | <u>Cálculo del sistema de ecuaciones de cada elemento en coordenadas globales</u> | ||
− | Las ecuaciones [[#eq-150|150]] y [[#eq-151|151]] se transforman a coordenadas globales como: | + | Las ecuaciones [[#eq-150|(150)]] y [[#eq-151|(151)]] se transforman a coordenadas globales como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,705: | Line 2,757: | ||
|} | |} | ||
− | + | donde se emplearon los valores de los desplazamientos nodales presentados en [[#eq-155|(155)]]. | |
Como revisión de las reacciones obtenidas anteriormente se realizará el equilibrio de toda la estructura: | Como revisión de las reacciones obtenidas anteriormente se realizará el equilibrio de toda la estructura: | ||
Line 2,739: | Line 2,791: | ||
|} | |} | ||
− | Mientras que a partir de lo presentado en [[#eq-22|22]] y [[#eq-147|147]], se tiene que el campo de desplazamiento empotrado en dirección <math display="inline">{x'}_A</math> del elemento <math display="inline">A</math> es: | + | Mientras que a partir de lo presentado en [[#eq-22|(22)]] y [[#eq-147|(147)]], se tiene que el campo de desplazamiento empotrado en dirección <math display="inline">{x'}_A</math> del elemento <math display="inline">A</math> es: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,758: | Line 2,810: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{u_A(x'}_A)= \begin{cases}\dfrac{Q}{E} \left[-\dfrac{1454}{9} \dfrac{{x'}_A}{L}+96 \left(\dfrac{{x'}_A}{L} \right)^2 -64\left(\dfrac{{x'}_A}{L} \right)^3 \right]& 0 | + | | style="text-align: center;" | <math>{u_A(x'}_A)= \begin{cases}\dfrac{Q}{E} \left[-\dfrac{1454}{9} \dfrac{{x'}_A}{L}+96 \left(\dfrac{{x'}_A}{L} \right)^2 -64\left(\dfrac{{x'}_A}{L} \right)^3 \right]&{ 0 < x'}_A \leq \dfrac{L}{2} \\[0.4cm] \dfrac{Q}{E} \left[24-\dfrac{2318}{9}\dfrac{{x'}_A}{L}+192\left(\dfrac{{x'}_A}{L} \right)^2 -64\left(\dfrac{{x'}_A}{L} \right)^3 \right]& \dfrac{L}{2} \leq {x'} < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (160) | | style="width: 5px;text-align: right;white-space: nowrap;" | (160) | ||
Line 2,865: | Line 2,917: | ||
<u>Cálculo de las fuerzas internas</u> | <u>Cálculo de las fuerzas internas</u> | ||
− | A partir de la derivación de los campos de desplazamientos en coordenadas locales de ambos elementos se obtiene las fuerzas internas en estos (ecuaciones [[#eq-2|2]] y [[# | + | A partir de la derivación de los campos de desplazamientos en coordenadas locales de ambos elementos se obtiene las fuerzas internas en estos (ecuaciones [[#eq-2|(2)]] y [[#eq-4|(4)]]): |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,872: | Line 2,924: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{P_A(x'}_A)=AE \dfrac{{du'}_A}{{dx'}_A}{(x'}_A)= \begin{cases}QL \left[-\dfrac{727}{1800} + \dfrac{12}{25}\dfrac{{x'}_A}{L}-\dfrac{12}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] \qquad & | + | | style="text-align: center;" | <math>{P_A(x'}_A)=AE \dfrac{{du'}_A}{{dx'}_A}{(x'}_A)= \begin{cases}QL \left[-\dfrac{727}{1800} + \dfrac{12}{25}\dfrac{{x'}_A}{L}-\dfrac{12}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] \qquad &{ 0 < x'}_A < \dfrac{L}{2} \\[0.4cm] QL \left[-\dfrac{1159}{1800} + \dfrac{24}{25}\dfrac{{x'}_A}{L}-\dfrac{12}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] & \dfrac{L}{2}{ < x'}L < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (170) | | style="width: 5px;text-align: right;white-space: nowrap;" | (170) | ||
Line 2,882: | Line 2,934: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{V_A(x'}I )=-EI \dfrac{{d^3 v'}_A}{\left.{dx'}_A \right.^3}{(x'}_A)= \begin{cases}QL \left[-\dfrac{14}{75} + \dfrac{16}{25}\dfrac{{x'}_A}{L}-\dfrac{16}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] \qquad & | + | | style="text-align: center;" | <math>{V_A(x'}I )=-EI \dfrac{{d^3 v'}_A}{\left.{dx'}_A \right.^3}{(x'}_A)= \begin{cases}QL \left[-\dfrac{14}{75} + \dfrac{16}{25}\dfrac{{x'}_A}{L}-\dfrac{16}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] \qquad &{ 0 < x'}_A \leq \dfrac{L}{2} \\[0.4cm] QL \left[-\dfrac{38}{75} + \dfrac{32}{25}\dfrac{{x'}_A}{L}-\dfrac{16}{25} \left(\dfrac{{x'}_A}{L}\right)^2 \right] & \dfrac{L}{2} \leq {x'} < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (171) | | style="width: 5px;text-align: right;white-space: nowrap;" | (171) | ||
Line 2,892: | Line 2,944: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{M_A(x'}_A)=EI \dfrac{{d^2 v'}_A}{\left.{dx'}_A \right.^2}{(x'}_A) \begin{cases}QL^2 \left[\dfrac{14}{75}\dfrac{{x'}_A}{L} - \dfrac{8}{25}\left(\dfrac{{x'}_A}{L}\right)^2 +\dfrac{16}{75} \left(\dfrac{{x'}_A}{L}\right)^3 \right]\qquad & | + | | style="text-align: center;" | <math>{M_A(x'}_A)=EI \dfrac{{d^2 v'}_A}{\left.{dx'}_A \right.^2}{(x'}_A) \begin{cases}QL^2 \left[\dfrac{14}{75}\dfrac{{x'}_A}{L} - \dfrac{8}{25}\left(\dfrac{{x'}_A}{L}\right)^2 +\dfrac{16}{75} \left(\dfrac{{x'}_A}{L}\right)^3 \right]\qquad &{ 0 < x'}_A \leq \dfrac{L}{2} \\[0.3cm] QL^2 \left[-\dfrac{2}{25}+\dfrac{38}{75}\dfrac{{x'}_A}{L} - \dfrac{16}{25}\left(\dfrac{{x'}_A}{L}\right)^2 +\dfrac{16}{75} \left(\dfrac{{x'}_A}{L}\right)^3 \right]& \dfrac{L}{2} \leq {x'} < L \end{cases} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (172) | | style="width: 5px;text-align: right;white-space: nowrap;" | (172) | ||
Line 2,913: | Line 2,965: | ||
|} | |} | ||
− | Como resumen en las | + | Como resumen en las [[#img-17a|Figuras 17a]] a [[#img-17c|17c]] se presentan los diagramas de fuerzas internas en los dos elementos que componen la estructura. |
− | + | ||
− | + | ||
<div id='img-17a'></div> | <div id='img-17a'></div> | ||
<div id='img-17b'></div> | <div id='img-17b'></div> | ||
− | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: | + | <div id='img-17c'></div> |
+ | <div id='img-17d'></div> | ||
+ | <div id='img-17'></div> | ||
+ | {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 85%;" | ||
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo3CampoFuerzaAxialDefinitivo.png|456px|Fuerza axial dividida QL, (\dfracPQL ).]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo3CampoFuerzaCortanteDefinitivo.png|456px|Fuerza cortante dividida QL, (\dfracVQL ).]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | ( | + | | (a) Fuerza axial dividida <math display="inline">QL</math>, <math display="inline">\left(\dfrac{P}{QL} \right)</math> |
− | | ( | + | | (b) Fuerza cortante dividida <math display="inline">QL</math>, <math display="inline">\left(\dfrac{V}{QL} \right)</math> |
|- | |- | ||
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo3CampoMomentoFlectorDefinitivo.png|456px|Momento flector dividido QL², (\dfracMQL² ).]] |
− | |[[Image: | + | |style="padding:10px;" |[[Image:Molina-Villegas_et_al_2020b-FiguraEjemplo3CampoDesplazamientoDefinitivo.png|456px|Configuración deformada de la estructura, los desplazamientos adimensionalizados han sido escalados 0.00002 veces.]] |
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | ( | + | | (c) Momento flector dividido <math display="inline">QL^2</math>, <math display="inline">\left(\dfrac{M}{QL^2} \right)</math> |
+ | | (d) Configuración deformada de la estructura, los desplazamientos adimensionalizados han sido escalados 0.00002 veces | ||
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="2" | '''Figura | + | | colspan="2" style="padding:10px;" | '''Figura 17'''. Diagramas de fuerza internas adimensionalizadas y configuración deformada de la estructura |
|} | |} | ||
+ | |||
+ | |||
En el siguiente enlace se encuentra el código de Python empleado para la solución de este ejemplo: | En el siguiente enlace se encuentra el código de Python empleado para la solución de este ejemplo: | ||
https://drive.google.com/open?id=1WNnbqC-09Odkh3GajgONgi8Al1sKhRz_ | https://drive.google.com/open?id=1WNnbqC-09Odkh3GajgONgi8Al1sKhRz_ | ||
− | ==9 Conclusiones== | + | ==9. Conclusiones== |
<ol> | <ol> | ||
<li>Se ha presentado una extensión del método de rigidez para estructuras reticulares planas, la cual puede ser aplicable para obtener la respuesta total de este tipo de estructural ante cualquier carga externa. Esta metodología ha sido nombrada el método de rigidez con funciones de Green (Green function stiffness method) y a partir de sus ideas básicas, esta metodología puede expandirse a otras configuraciones estructurales. </li> | <li>Se ha presentado una extensión del método de rigidez para estructuras reticulares planas, la cual puede ser aplicable para obtener la respuesta total de este tipo de estructural ante cualquier carga externa. Esta metodología ha sido nombrada el método de rigidez con funciones de Green (Green function stiffness method) y a partir de sus ideas básicas, esta metodología puede expandirse a otras configuraciones estructurales. </li> | ||
− | <li>La metodología presentada, al ser definida de forma | + | <li>La metodología presentada, al ser definida de forma explícita para cualquier carga externa, permite minimizar el número de elementos que se emplearían respecto a aquellos usando el método de rigidez o el M.E.F. Como caso particular se tiene que en los textos clásicos de análisis matricial de estructuras <span id='citeF-15'></span><span id='citeF-16'></span><span id='citeF-17'></span>[[#cite-15|[15,16,17]]] es usual subdividir los elementos si las cargas externas son complejas y estas no aparecen en sus tablas de fuerzas de empotramiento, lo cual no es necesario con la presente metodología. </li> |
<li>Se extiende el concepto de fuerzas de empotramiento típico del método de rigidez al concepto de campo de desplazamientos empotrado, el cual, al emplear las funciones de Green permite el cálculo de los campos de desplazamientos incluso cuando las fuerzas externas son complejas. </li> | <li>Se extiende el concepto de fuerzas de empotramiento típico del método de rigidez al concepto de campo de desplazamientos empotrado, el cual, al emplear las funciones de Green permite el cálculo de los campos de desplazamientos incluso cuando las fuerzas externas son complejas. </li> | ||
<li>En la metodología presentada se da una gran importancia al cálculo de los campos de desplazamientos. Esto se debe a que los desplazamientos son las variables dependientes principales de los problemas de análisis estructural e incluso a partir de estos es posible calcular los campos de fuerzas internas. </li> | <li>En la metodología presentada se da una gran importancia al cálculo de los campos de desplazamientos. Esto se debe a que los desplazamientos son las variables dependientes principales de los problemas de análisis estructural e incluso a partir de estos es posible calcular los campos de fuerzas internas. </li> | ||
<li>Se presenta una forma novedosa de expresar a las funciones de Green en términos de las funciones de forma “exactas” para los diferentes tipos de elementos estructurales estudiados. Esta se base propiedades de simetría observadas tanto en las funciones de forma como en las funciones de Green. </li> | <li>Se presenta una forma novedosa de expresar a las funciones de Green en términos de las funciones de forma “exactas” para los diferentes tipos de elementos estructurales estudiados. Esta se base propiedades de simetría observadas tanto en las funciones de forma como en las funciones de Green. </li> | ||
− | <li>Se hace | + | <li>Se hace explícita la gran importancia que tienen las funciones de forma “exactas” en la respuesta estructural, ya que estas no solo son la base para el cálculo del campo de desplazamiento homogéneo sinó también del empotrado pues aparecen de forma natural en las funciones de Green y en las fuerzas de empotramiento. </li> |
<li>En los apéndices se presenta la equivalencia entre la metodología presentada en este artículo, la cual parte de las ecuaciones diferenciales gobernantes para cuatro tipos de elementos diferentes y aquella obtenida por medio de una formulación “exacta” del M.E.F. (la cual parte de la forma débil de dichas ecuaciones diferenciales). </li> | <li>En los apéndices se presenta la equivalencia entre la metodología presentada en este artículo, la cual parte de las ecuaciones diferenciales gobernantes para cuatro tipos de elementos diferentes y aquella obtenida por medio de una formulación “exacta” del M.E.F. (la cual parte de la forma débil de dichas ecuaciones diferenciales). </li> | ||
</ol> | </ol> | ||
− | ==10 Agradecimientos== | + | ==10. Agradecimientos== |
− | 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 | + | 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 ha 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 autores. 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. |
==Apéndices== | ==Apéndices== | ||
− | + | ===A. Formulación exacta del M.E.F. para un elemento tipo barra=== | |
− | + | ||
− | == | + | |
− | 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 Formulación analítica del método de rigidez con funciones de Green para elementos tipo barra|4]], la cual también se basa en la solución del P.V.F [[# | + | 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 Formulación analítica del método de rigidez con funciones de Green para elementos tipo barra|4]], la cual también se basa en la solución del P.V.F [[#eq-8|(8)]] pero ahora a partir de la llamada forma débil de [[#eq-8.a|(8.a)]] en lugar de esta directamente. De forma similar a como se indicó para el método de rigidez con funciones de Green, el campo de desplazamientos se expresará como: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 2,973: | Line 3,028: | ||
|} | |} | ||
− | + | donde <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math> son las funciones de forma presentadas en [[#eq-12|(12)]] y <math display="inline">u_f(x)</math> se calcula según lo presentado en [[#eq-22|(22)]]. | |
− | Además, es fácil probar que las funciones de forma <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math> cumplen la siguiente propiedad asociada con [[#eq-10.a|10.a]] (ecuación diferencial del campo homogeneo): | + | Además, es fácil probar que las funciones de forma <math display="inline">\psi _1(x)</math> y <math display="inline">\psi _4(x)</math> cumplen la siguiente propiedad asociada con [[#eq-10.a|(10.a)]] (ecuación diferencial del campo homogeneo): |
<span id="eq-175"></span> | <span id="eq-175"></span> | ||
Line 2,983: | Line 3,038: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\dfrac{d^2 \psi _i}{dx^2}(x)=0 \ | + | | style="text-align: center;" | <math>\dfrac{d^2 \psi _i}{dx^2}(x)=0 \quad i=1,4 </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (175) | | style="width: 5px;text-align: right;white-space: nowrap;" | (175) | ||
|} | |} | ||
− | Multiplicando a ambos lados de [[#eq-8.a|8.a]] por una función de peso o ponderación <math display="inline">w(x)</math> se obtiene: | + | Multiplicando a ambos lados de [[#eq-8.a|(8.a)]] por una función de peso o ponderación <math display="inline">w(x)</math> se obtiene: |
<span id="eq-176"></span> | <span id="eq-176"></span> | ||
Line 3,001: | Line 3,056: | ||
|} | |} | ||
− | Ahora, integrando a ambos lados de [[#eq-176|176]] con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> se obtiene: | + | Ahora, integrando a ambos lados de [[#eq-176|(176)]] con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> se obtiene: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,050: | Line 3,105: | ||
|} | |} | ||
− | Lo anterior debido a que <math display="inline">u_f(0)=u_f(L)=0</math> y a [[#eq-175|175]]. | + | Lo anterior debido a que <math display="inline">u_f(0)=u_f(L)=0</math> y a [[#eq-175|(175)]]. |
− | Con lo cual [[#eq-179|179]] se puede reescribir como: | + | Con lo cual [[#eq-179|(179)]] se puede reescribir como: |
<span id="eq-181"></span> | <span id="eq-181"></span> | ||
Line 3,065: | Line 3,120: | ||
|} | |} | ||
− | Ahora, si como caso particular primero se emplea <math display="inline">\psi _1(x)</math> en lugar de <math display="inline">w(x)</math> y luego <math display="inline">\psi _4(x)</math>, a partir de [[#eq-181|181]] se obtiene respectivamente: | + | Ahora, si como caso particular primero se emplea <math display="inline">\psi _1(x)</math> en lugar de <math display="inline">w(x)</math> y luego <math display="inline">\psi _4(x)</math>, a partir de [[#eq-181|(181)]] se obtiene respectivamente: <span id="eq-182"></span> |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,080: | Line 3,135: | ||
|} | |} | ||
− | 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 [[#eq-26|26]]. | + | 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 [[#eq-26|(26)]]. |
− | == | + | ===B. Formulación exacta del M.E.F. para un elemento tipo viga=== |
− | Ahora, de forma similar a como se realizó en el apéndice [[# | + | Ahora, de forma similar a como se realizó en el apéndice [[#A Formulación exacta del M.E.F. para un elemento tipo barra|A]] 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 Formulación analítica del método de rigidez con funciones de Green para elementos tipo viga|5]], el objetivo es realizar esta formulación mediante la solución del P.V.F. [[#eq-27|(27)]] pero ahora empleando la forma de débil de la E.D. [[#eq-27.a|(27.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: |
<span id="eq-183"></span> | <span id="eq-183"></span> | ||
Line 3,097: | Line 3,152: | ||
|} | |} | ||
− | + | donde una propiedad importante de las funciones de forma <math display="inline">\psi _2(x)</math>, <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _6(x)</math> es que estas cumplen la siguiente ecuación diferencial: | |
<span id="eq-184"></span> | <span id="eq-184"></span> | ||
Line 3,105: | Line 3,160: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\dfrac{d^4 \psi _i}{dx^4}(x)=0 | + | | style="text-align: center;" | <math>\dfrac{d^4 \psi _i}{dx^4}(x)=0 i=2,3,5,6 </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (184) | | style="width: 5px;text-align: right;white-space: nowrap;" | (184) | ||
|} | |} | ||
− | Multiplicando por una función de peso o ponderación <math display="inline">w(x)</math> a ambos lados de [[#eq-27.a|27.a]] se obtiene: | + | Multiplicando por una función de peso o ponderación <math display="inline">w(x)</math> a ambos lados de [[#eq-27.a|(27.a)]] se obtiene: |
<span id="eq-185"></span> | <span id="eq-185"></span> | ||
Line 3,123: | Line 3,178: | ||
|} | |} | ||
− | E integrando con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> a ambos lados de [[#eq-185|185]] se obtiene: | + | E integrando con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> a ambos lados de [[#eq-185|(185)]] se obtiene: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,196: | Line 3,251: | ||
|} | |} | ||
− | Reemplazando [[#eq-183|183]] en el lado derecho de [[#eq-191|191]], se obtiene: | + | Reemplazando [[#eq-183|(183)]] en el lado derecho de [[#eq-191|(191)]], se obtiene: |
<span id="eq-192"></span> | <span id="eq-192"></span> | ||
Line 3,232: | Line 3,287: | ||
|} | |} | ||
− | + | donde se ha empleado [[#eq-184|(184)]] y que los valores de <math display="inline">v_f(x)</math> y <math display="inline">\dfrac{dv_f}{dx}(x)</math> son iguales a cero en <math display="inline">x=0</math> y <math display="inline">x=L</math>. | |
− | Con lo cual, a partir de [[#eq-193|193]], la ecuación [[#eq-192|192]] se puede reescribir como: | + | Con lo cual, a partir de [[#eq-193|(193)]], la ecuación [[#eq-192|(192)]] se puede reescribir como: |
<span id="eq-194"></span> | <span id="eq-194"></span> | ||
Line 3,249: | Line 3,304: | ||
|} | |} | ||
− | Si como caso particular primero se emplea <math display="inline">\psi _2(x)</math>, luego <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _6(x)</math>, a partir de [[#eq-194|194]] se obtienen las siguientes cuatro ecuaciones: | + | Si como caso particular primero se emplea <math display="inline">\psi _2(x)</math>, luego <math display="inline">\psi _3(x)</math>, <math display="inline">\psi _5(x)</math> y <math display="inline">\psi _6(x)</math>, a partir de [[#eq-194|(194)]] se obtienen las siguientes cuatro ecuaciones: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,270: | Line 3,325: | ||
|} | |} | ||
− | Que son equivalentes al sistema de ecuaciones presentado en [[#eq-52|52]]. | + | Que son equivalentes al sistema de ecuaciones presentado en [[#eq-52|(52)]]. |
− | == | + | ===C. Formulación exacta del M.E.F. para un elemento tipo viga sobre fundación flexible=== |
− | Para este caso el objetivo es resolver el P.V.F. [[# | + | Para este caso el objetivo es resolver el P.V.F. [[#eq-53|(53)]] pero ahora empleando la forma débil de [[#eq-53.a|(53.a)]] en lugar de esta directamente. |
− | Como se presentó en [[#eq-83|83]], la solución de [[#eq-53.a|53.a]] se expresa como la suma de un campo homogéneo (<math display="inline">v_h(x)</math>) y un campo empotrado (<math display="inline">v_f(x)</math>), es decir: | + | Como se presentó en [[#eq-83|(83)]], la solución de [[#eq-53.a|(53.a)]] se expresa como la suma de un campo homogéneo (<math display="inline">v_h(x)</math>) y un campo empotrado (<math display="inline">v_f(x)</math>), es decir: |
<span id="eq-196"></span> | <span id="eq-196"></span> | ||
Line 3,289: | Line 3,344: | ||
|} | |} | ||
− | + | 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 <math display="inline">\Psi _2(x)</math>, <math display="inline">\Psi _3(x)</math>, <math display="inline">\Psi _5(x)</math> y <math display="inline">\Psi _6(x)</math> de la siguiente manera: | |
<span id="eq-197"></span> | <span id="eq-197"></span> | ||
Line 3,317: | Line 3,372: | ||
|} | |} | ||
− | Es importante resaltar que debido a que <math display="inline">v_h(x)</math> es una combinación lineal de las funciones de forma <math display="inline">\Psi _2(x)</math>, <math display="inline">\Psi _3(x)</math>, <math display="inline">\Psi _5(x)</math> y <math display="inline">\Psi _6(x)</math>, estas últimas cumplen la propiedad de ser solución de la ecuación diferencial homogénea asociada con [[#eq-53.a|53.a]], es decir: | + | Es importante resaltar que debido a que <math display="inline">v_h(x)</math> es una combinación lineal de las funciones de forma <math display="inline">\Psi _2(x)</math>, <math display="inline">\Psi _3(x)</math>, <math display="inline">\Psi _5(x)</math> y <math display="inline">\Psi _6(x)</math>, estas últimas cumplen la propiedad de ser solución de la ecuación diferencial homogénea asociada con [[#eq-53.a|(53.a)]], es decir: |
<span id="eq-199"></span> | <span id="eq-199"></span> | ||
Line 3,325: | Line 3,380: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>EI\dfrac{d^4 \Psi _i}{d x^4}(x)+k\Psi _i(x)=0 \ | + | | style="text-align: center;" | <math>EI\dfrac{d^4 \Psi _i}{d x^4}(x)+k\Psi _i(x)=0 \quad i=2,3,5,6 </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (199) | | style="width: 5px;text-align: right;white-space: nowrap;" | (199) | ||
|} | |} | ||
− | Multiplicando por una función de peso o ponderación <math display="inline">w(x)</math> a ambos lados de [[#eq-5|5]] (o [[#eq-53.a|53.a]]) se obtiene: | + | Multiplicando por una función de peso o ponderación <math display="inline">w(x)</math> a ambos lados de [[#eq-5|(5)]] (o [[#eq-53.a|(53.a)]]) se obtiene: |
<span id="eq-200"></span> | <span id="eq-200"></span> | ||
Line 3,343: | Line 3,398: | ||
|} | |} | ||
− | E integrando con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> a ambos lados de [[#eq-200|200]] se obtiene: | + | E integrando con respecto a <math display="inline">x</math> entre 0 y <math display="inline">L</math> a ambos lados de [[#eq-200|(200)]] se obtiene: |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,420: | Line 3,475: | ||
|} | |} | ||
− | Reemplazando [[#eq-198|198]] en el lado derecho de [[#eq-206|206]], se obtiene: | + | Reemplazando [[#eq-198|(198)]] en el lado derecho de [[#eq-206|(206)]], se obtiene: |
<span id="eq-207"></span> | <span id="eq-207"></span> | ||
Line 3,460: | Line 3,515: | ||
|} | |} | ||
− | + | donde se ha empleado [[#eq-199|(199)]] y que los valores de <math display="inline">v_f(x)</math> y <math display="inline">\dfrac{dv_f}{dx}(x)</math> son iguales a cero en <math display="inline">x=0</math> y <math display="inline">x=L</math>. | |
− | Empleando [[#eq-208|208]] en [[#eq-207|207]] esta última toma la siguiente forma (donde para ser coherentes con [[#eq-208|208]], se ha tomado que <math display="inline">w(x)=\Psi _l(x)</math>): | + | Empleando [[#eq-208|(208)]] en [[#eq-207|(207)]] esta última toma la siguiente forma (donde para ser coherentes con [[#eq-208|(208)]], se ha tomado que <math display="inline">w(x)=\Psi _l(x)</math>): |
<span id="eq-209"></span> | <span id="eq-209"></span> | ||
Line 3,479: | Line 3,534: | ||
|} | |} | ||
− | Reemplazando <math display="inline">l</math> por 2,3,5 y 6 en [[#eq-209|209]] se obtienen cuatro ecuaciones, las cuales se expresan en forma matricial como: | + | Reemplazando <math display="inline">l</math> por 2,3,5 y 6 en [[#eq-209|(209)]] se obtienen cuatro ecuaciones, las cuales se expresan en forma matricial como: |
<span id="eq-210"></span> | <span id="eq-210"></span> | ||
Line 3,492: | Line 3,547: | ||
|} | |} | ||
− | + | donde: | |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 3,504: | Line 3,559: | ||
|} | |} | ||
− | Lo cual es equivalente a lo presentado en [[#eq-84|84]]. | + | Lo cual es equivalente a lo presentado en [[#eq-84|(84)]]. |
− | === | + | ==Referencias== |
+ | <div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;"> | ||
<div id="cite-1"></div> | <div id="cite-1"></div> | ||
− | + | [[#citeF-1|[1]]] Challis L., Sheard F. The green of Green functions. American Institute of Physics, Physics Today, 56(12):41-46, 2003. | |
<div id="cite-2"></div> | <div id="cite-2"></div> | ||
− | + | [[#citeF-2|[2]]] Duffy D.G. Green’s functions with applications. Chapman and Hall/CRC, 2015. | |
<div id="cite-3"></div> | <div id="cite-3"></div> | ||
− | + | [[#citeF-3|[3]]] Banerjee P.K., Butterfield R. Boundary element methods in engineering science. McGraw-Hill, New York, 1981. | |
<div id="cite-4"></div> | <div id="cite-4"></div> | ||
− | + | [[#citeF-4|[4]]] Sánchez-Sesma F.J., Ramos-Martínez J., 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., Karageorghis A., 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 W.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 J. Application des potentiels a l'étude de l'équilibre et du mouvement des solides élastiques. Gauthier-Villars, Imprimeur-Libraire, 1885. | |
<div id="cite-8"></div> | <div id="cite-8"></div> | ||
− | + | [[#citeF-8|[8]]] Cerruti V. Ricerche intorno all'equilibrio de'corpi elastici isotropi: memoria. Coi tipi del Salviucci, 1882. | |
<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. AIP Physics, 7(5):195-202, 19336. | |
<div id="cite-10"></div> | <div id="cite-10"></div> | ||
− | + | [[#citeF-10|[10]]] Stokes G.G. 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]]] Lamb H. 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. Journal of Applied Mechanics, 27:559, 1960. | |
<div id="cite-13"></div> | <div id="cite-13"></div> | ||
− | + | [[#citeF-13|[13]]] Kausel E. Fundamental solutions in elastodynamics: a compendium. Cambridge University Press, 2006. | |
<div id="cite-14"></div> | <div id="cite-14"></div> | ||
− | + | [[#citeF-14|[14]]] Aki K., Richards P.G. Quantitative seismology. Univ. Science Books, 700 pp., 2002. | |
<div id="cite-15"></div> | <div id="cite-15"></div> | ||
− | + | [[#citeF-15|[15]]] Colunga A.T. Análisis de estructuras con métodos matriciales. Limusa, 560 pp., 2007. | |
<div id="cite-16"></div> | <div id="cite-16"></div> | ||
− | + | [[#citeF-16|[16]]] McCormac J.C. Structural Analysis: using classical and matrix methods. Wiley Hoboken, NJ, 2007. | |
<div id="cite-17"></div> | <div id="cite-17"></div> | ||
− | + | [[#citeF-17|[17]]] Kassimali A. Matrix analysis of structures SI version. Cengage Learning, 640 pp., 2012. | |
<div id="cite-18"></div> | <div id="cite-18"></div> | ||
− | + | [[#citeF-18|[18]]] Reddy J.N. An introduction to the finite element method. McGraw-Hill, 2006. | |
<div id="cite-19"></div> | <div id="cite-19"></div> | ||
− | + | [[#citeF-19|[19]]] Bathe K.J. Finite element procedures. McGraw-Hill, 2006. | |
<div id="cite-20"></div> | <div id="cite-20"></div> | ||
− | + | [[#citeF-20|[20]]] Hetényi M. Beams on elastic foundation: theory with applications in the fields of civil and mechanical engineering. University of Michigan, 1971. | |
<div id="cite-21"></div> | <div id="cite-21"></div> | ||
− | + | [[#citeF-21|[21]]] Eisemberger M., Yankelevsky D.Z. Exact stiffness matrix for beams on elastic foundation. Computer & Structures, 21(6):1355-1359, 1985. | |
+ | </div> |
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 líneas de influencia y en algunos casos esporádicos en la formulación del método del trabajo virtual. En este artículo se presenta el método de rigidez con funciones de Green (Green functions stiffness method) el cual es una metodología novedosa para obtener la respuesta analítica o exacta de estructuras reticulares planas, y como su nombre lo indica mezcla al método de rigidez con las funciones de Green. 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 obtener la respuesta analítica ante 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.
Palabras clave: Método de rigidez con funciones de Green, Funciones de Green, Método de rigidez, Método de elementos finitos, Vigas de Bernoulli, Pórticos planos, Vigas sobre fundación flexible, Funciones de forma
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 the Green functions stiffness method, which is a novel 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 to compute the analytic reponse 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: Green function stiffness method, Green functions, stiffness method, finite element method, Bernoulli beam, plane frames, beams on elastic foundations, shape functions
Las funciones de Green (F.G.) son la respuesta de un medio o elemento 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 sean lineales [1]. Esto es debido a que ellas por sí solas son solución de problemas fundamentales, pueden ser empleadas para resolver problemas con fuentes o fuerzas 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,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 líneas de influencia y en algunos casos esporádicos 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 , la cual corresponde al desplazamiento en dirección del punto , debido a una fuerza puntual unitaria en dirección aplicada en el punto .
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 más complejo de implementar y busca la solución exacta del modelo estructural en estudio, mientras que el segundo es más 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 de geometrías simples, como puntuales o distribuidas con variación lineal. Mientras que el M.E.F. suele emplearse para el análisis de 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,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,19].
En este artículo se presenta una formulación particular del método de rigidez para obtener la respuesta analítica de problemas de análisis estructural formados por estructuras reticulares planas (pudiendo estar formadas por elementos tipo barra, viga, viga sobre fundación flexible o pórtico plano), la cual se llamará el método de rigidez con funciones de Green (Green functions stiffness method), y emplea a las funciones de Green para realizar el cálculo de los desplazamientos al interior de los elementos. 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 cálculo 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 plano 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 concluirá con la realización tres ejemplos ilustrativos (para una viga, una viga sobre fundación flexible y un pórtico plano respectivamente).
A continuación se definirán las propiedades de los elementos tipo barra, viga, pórtico plano y viga sobre fundación flexible a estudiar en este documento, al igual que sus ecuaciones diferenciales gobernantes. Como generalidad, para todos estos tipos elementos se empleará la convención de fuerzas internas positiva presentada en la Figura 1, la cual se inspira en aquella usada para la definición de los esfuerzos positivos en la mecánica del medio continuo o de sólidos.
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 (), la cual actúa en dirección del eje , y su valor es positivo si se dirige en la dirección de dicho eje local (Figura 2). Para el caso que el elemento sea de sección transversal constante con área de la sección transversal , material elástico lineal con módulo de elasticidad , su ecuación diferencial gobernante es:
|
(1) |
donde es el desplazamiento en dirección axial del elemento (eje local ).
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) |
donde es el esfuerzo axial y es la deformación unitaria axial.
El elemento tipo viga está sometido tanto a fuerza cortante en dirección del eje local , como a momento flector alrededor del eje local (perpendicular a los ejes y , Figura 3). Sus cargas externas se definen en términos de la carga por unidad de longitud en dirección del eje local (), la cual es positiva en la dirección de este eje (Figura 3). Si el elemento es de sección transversal constante, con momento de inercia y material elástico lineal con módulo de elasticidad , su ecuación diferencial gobernante es:
|
(3) |
donde es el desplazamiento en dirección del eje local , 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 , respectivamente como:
|
(4.a) | |
|
(4.b) |
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 y genera una fuerza por unidad de longitud sobre la viga igual a (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, para este elemento las fuerzas internas se calculan a partir del campo de desplazamientos empleando (4).
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 local . Sus ecuaciones diferenciales gobernantes son (1) para el desplazamiento axial y (3) para el desplazamiento perpendicular al eje longitudinal del elemento (eje ) y sus fuerzas internas se calculan a partir de los campos de desplazamiento empleando (2) y (4).
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 puntos (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 es el vector de fuerzas en los extremos del elemento en coordenadas locales, es la matriz de rigidez en coordenadas locales, es el vector de desplazamientos en los extremos del elemento en coordenadas locales y 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 para los elementos estudiados se expresan como:
|
(7) |
donde es el componente del campo de desplazamiento, el omponente del campo homogéneo de desplazamiento y el omponente del campo empotrado de desplazamiento.
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 con funciones de Green 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 es (Figura (2)):
|
donde es el desplazamiento en dirección axial (eje local ) en el extremo inicial del elemento y 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 es el campo de desplazamiento homogéneo (sección 4.1), mientras que es el campo empotrado (sección 4.2).
El P.V.F. que gobierna al campo homogéneo se presenta en (10) 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 (10) es el siguiente polinomio de primer grado:
|
(11) |
donde:
|
y son conocidas como las funciones de forma “exactas” de este problema y en este caso son iguales a aquellas que suelen emplearse en el M.E.F. para interpolar el campo de desplazamientos en un elemento tipo barra. Un par de propiedades muy importantes de y , las cuales nos permitirán expresar de forma compacta algunos cálculos posteriores, son:
|
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:
|
(14) |
Mientras que las fuerzas en dirección del eje local en los extremos de la barra son:
|
El problema de P.V.F. que gobierna al campo empotrado se presenta en (16). Su nombre se debe a que representa la respuesta de una barra doblemente empotrada sometida a una carga externa por unidad de longitud igual a la del elemento en estudio (Figura 5a)
|
Para la solución del P.V.F. (16) se empleará la función de Green de un elemento tipo barra doblemente empotrada, que es de sección transversal constante y material elástico lineal (Figura 5b). Esta corresponde a la respuesta (campo de desplazamiento axial) debido a la aplicación de una carga puntual unitaria ubicada a una distancia 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 ):
|
donde es la función delta de Dirac y lo superíndices - y + significan a la izquierda y a la derecha del punto en cuestión.
La solución de (17) expresada en términos de las funciones de forma (12), se puede escribir de forma compacta como:
|
(18) |
donde, para su escritura se han explotado las propiedades de simetría de las funciones de forma presentadas en (13), así como las siguientes propiedades de simetría la función de Green:
|
Como se verá más adelante, es importante reescribir a en los intervalos y , lo cual da como resultado:
|
(20) |
Multiplicando a ambos lados de (17.a) por e integrando respecto a entre 0 y se obtiene:
|
(21) |
donde se empleó la siguiente propiedad de la función delta de Dirac: ().
De comparar (21) y (16.a), así como empleando (20) se concluye que:
|
(22) |
La cual también puede interpretarse mediante el principio de superposición pues el desplazamiento en cada punto de la barra doblemente empotrada es la suma de los desplazamientos causados en ese mismo punto por infinitas cargas puntuales ubicadas en desde 0 hasta y cuyo valor es , donde los limites de ambas integrales se explican por (20).
Empleando la ley de Hooke uniaxial, a partir de (22) el campo de fuerza axial empotrado se calcula como:
|
(23) |
Y las reacciones o fuerzas de empotramiento en dirección del eje local son:
|
Es de destacar que la anterior forma de 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., la cual se encuentra en el Apéndice A.
Reemplazando (11) y (22) 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:
|
(25) |
Mientras que a partir de sumar (15) y (24), se tiene que las fuerzas en dirección del eje local en los extremos del elemento, se expresan en términos de los desplazamientos en esos mismos puntos y en la misma dirección axial:
|
(26) |
La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo barra y es equivalente a la presentada en el Apéndice A 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. (10) y (16) no han sido arbitrarias y son totalmente coherentes con (8) y (9). Por ejemplo si se reemplaza (9) en (8) y se emplean (16), se obtiene directamente el P.V.F. (10), 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 con momento de inercia , material elástico lineal homogéneo con módulo de elasticidad y longitud es (Figura 3):
|
donde es el desplazamiento en dirección del eje local en extremo inicial del elemento, es la pendiente en ese mismo punto, es el desplazamiento en dirección del eje local en extremo final del elemento y es la rotación en ese último punto.
Para este caso la particularización de (7) es:
|
(28) |
donde es el campo de desplazamiento homogéneo (sección 5.1), mientras que 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:
|
(30) |
donde , , y se conocen coma las funciones de forma “exactas” de este problema y tienen el siguiente valor:
|
Fuera de lo anterior, las siguientes propiedades de simetría de las funciones de forma , , y son muy importantes de resaltar:
|
Al igual que para el caso de las barras, las funciones de forma “exactas” presentadas en (31) 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 (28).
A partir de (30) y (4), las fuerzas internas homogéneas son:
|
Mientras que las fuerzas homogéneas en los extremos del elemento en dirección de los ejes locales , son:
|
El problema de P.V.F. que gobierna al campo empotrado se presenta en (35) y corresponde a la respuesta de una viga doblemente empotrada sometida a una carga externa genérica definida por la función (Figura 6a)
|
Para solucionar el P.V.F. (35) se definirá la función de 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 del extremo izquierdo del elemento (Figura 6b). Esta función de Green está gobernada por el P.V.F. (36), el cual comprende el cumplimiento de la ecuación diferencial gobernante, condiciones de frontera, condiciones de continuidad de desplazamientos y rotaciones en () y equilibrio vertical y rotacional en este mismo punto:
|
La solución de (36), es decir, la función de Green para la viga doblemente empotrada presentada en la Figura 6b es:
|
(37) |
La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en (32) y las siguientes propiedades de simetría de la función de Green:
|
La cual, por simplicidad en los próximos cálculos se reescribirá de forma compacta como:
|
(39) |
O de forma alternativa como:
|
(40) |
Multiplicando a ambos lados de (36.a) por e integrando respecto a entre 0 y se obtiene:
|
(41) |
donde se empleó que .
De comparar (21) y (35.a), así como empleando (40), se concluye que:
|
(42) |
Al igual que para el caso de la barra, el campo empotrado presentado en (42) puede interpretarse a partir del principio de superposición como la respuesta en un punto debido a la suma de los aportes de los diferenciales de fuerza ubicados a lo largo de toda la viga, para lo cual se emplea (40).
Mientras que los campos de fuerzas internas se calculan como:
|
(43) |
|
(44) |
Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales y son:
|
Por su parte para el caso en que se presente un momento puntual, pese a que este se puede expresar como una función empleando el límite de dos funciones delta de Dirac, por comodidad a continuación se presenta la función de Green para este caso (Figura 7):
|
(46) |
con lo cual:
|
(47) |
Mientras que los campos de fuerzas internas debidas al momento puntual son:
|
(48) |
|
(49) |
con lo cual las fuerzas de empotramiento para este caso son:
|
Figura 7. Viga de sección transversal constante doblemente empotrada, sometidas a un momento puntual unitario |
Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sinó de valor , los resultados anteriores se deberán simplemente multiplicar por , y en caso de que sea un momento por unidad de longitud distribuido, este se deberá integrar de forma similar a como se realizó con la carga pero ahora empleando en lugar de .
Reemplazando (30) y (42) en (28) se obtiene que el campo de desplazamientos puede expresarse en función de los desplazamientos nodales y las fuerzas externas como:
|
(51) |
Mientras que las fuerzas en los extremos del elemento se expresan en términos de los desplazamiento y rotaciones en esos mismos puntos a partir de la suma de (34) y (45):
|
(52) |
La cual es la formulación en coordenadas locales del método de rigidez con funciones de Green para un elemento tipo viga y es equivalente a la formulación “exacta” del M.E.F. presentada en el Apéndice B.
El P.V.F. que gobierna la formulación del método de rigidez en coordenadas locales 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 constante e igual y longitud es (Figura 4):
|
donde se ha empleado como ecuación diferencial una versión equivalente de (5), en la cual se define y , , y tienen el mismo significado que para la viga (sección 5)
Para este caso la particularización de (7) es:
|
(54) |
donde es el campo de desplazamiento homogéneo (sección 6.1), mientras que es el campo empotrado (sección 6.2).
Para este caso el P.V.F. gobernante es:
|
cuya solución se expresa como:
|
(56) |
donde , , y son las funciones de forma “exactas” de este problema y tienen el siguiente valor:
|
(57) |
|
(58) |
|
(59) |
|
(60) |
Y para simplificar la escritura, se ha definido , , y .
Al igual que para el caso de la viga, las funciones de forma para la viga sobre fundación flexible también poseen las siguientes propiedades de simetría:
|
A diferencia de lo que ocurre para el caso la barra y la viga, las funciones de forma “exactas” presentadas en (57) a (60) 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 (ecuaciones (31)). Como se presentará en el apéndice C, si en la formulación del M.E.F. para una viga sobre fundación flexible se emplean las funciones de forma (57) a (60), se obtiene una formulación matricial “exacta” y una matriz de rigidez y fuerzas de empotramientos iguales a los obtenidos en el método presentado en este artículo (método de rigidez con funciones de Green), siempre y cuando se sumen tanto el campo de desplazamiento homogéneo como el empotrado.
Para este caso el P.V.F. que define al problema empotrado es (Figura 8a):
|
Mientras que la función de Green se define por medio del siguiente P.V.F (Figura 8b):
|
donde se ha usado como convención emplear una linea horizontal encima de la letra para diferenciar esta función de Green de aquella de las vigas.
La solución de (63) se expresa en forma compacta empleando las funciones de forma (57) a (60), como:
|
(64) |
donde:
|
(65) |
|
(66) |
La cual se ha podido expresar de forma tan compacta empleando las propiedades se simetría de las funciones de forma presentadas en (61) y las siguientes propiedades de simetría de la misma función de Green:
|
Multiplicando a ambos lados de (63.a) por e integrando respecto a entre 0 y se obtiene:
|
(68) |
donde se empleó que .
De comparar (21) y (62.a) se concluye que:
|
(69) |
A partir de (69) y (4), se obtiene que las fuerzas internas se pueden calcular como:
|
(70) |
|
(71) |
Y las reacciones o fuerzas de empotramiento en dirección de los ejes locales se calculan como:
|
De forma similar a como se realizó para el caso de la viga, la función de Green debida a un momento puntual unitario se puede obtener directamente de aquella debida a una fuerza puntual mediante la derivada de la misma respecto a (Figura 9), es decir:
|
(73) |
donde:
|
(74) |
|
(75) |
Mientras que el campo de momento flector es:
|
(76) |
donde:
|
(77) |
|
(78) |
Y el campo de fuerza cortante es:
|
(79) |
donde:
|
(80) |
|
(81) |
Con lo cual, a partir de (76) y (79) se tiene que las fuerzas de empotramiento son:
|
Debido a que el problema estudiado es lineal, en caso de que el momento aplicado no sea unitario sino de valor , los resultados anteriores se deberán simplemente multiplicar por . Mientras que si se trata de un momento distruido, debe realizarse la integración del campo debido al momento puntual multiplicado por la función que define al momento distribuido.
Figura 9. Viga sobre funcdación flexible doblemente empotrada y de de sección transversal constante, sometidas a un momento puntual unitario |
Reemplazando (56) y (69) en (54), el campo de desplazamiento total se expresa en función de los desplazamientos de los extremos del elemento y de la carga externa como:
|
(83) |
Mientras que las fuerzas en los extremos del elemento se expresan en función de los desplazamientos en estos mismos puntos y de las cargas externas como:
|
(84) |
donde los términos de la matriz de rigidez son ([21]):
|
Es importante resaltar que (84) corresponde a la formulación analítica del método de rigidez en coordenadas locales 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 C.
Como se mencionó en la sección 2.4 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 (26) y (52) se tiene que la formulación del método de rigidez con funciones de Green en coordenadas locales para un elemento tipo pórtico plano es:
|
(86) |
Mientras que los campos de desplazamiento en dirección de los ejes locales y , se calculan a partir de las ecuaciones (25) y (51) respectivamente.
Calcular la respuesta (reacciones, campos de desplazamientos y campos de fuerzas internas) de la viga presentada en la Figura 10a.
(a) Problema a resolver | (b) Discretización |
Figura 10. 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 10b y consta de solo dos elementos. El elemento tiene una carga externa definida en tres tramos, mientras que el elemento tiene una carga externa lineal sobre toda su longitud.
Definición de las funciones de carga externa
Para el elemento su carga externa se define como:
|
(87) |
O en forma compacta como:
|
(88) |
Por su parte la carga externa en el elemento es:
|
(89) |
Cálculo de las fuerzas de empotramiento
Según lo presentado en (45), las fuerzas de empotramiento del elemento se calculan como:
|
Mientras que las del elemento como:
|
Sistema matricial de ecuaciones de cada elemento
A partir de lo presentado en (52), la formulación del método de rigidez para el elemento es:
|
(92) |
Mientras que para elemento este es:
|
(93) |
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:
|
(94) |
Cuya solución es:
|
(95) |
Cálculo del campo de desplazamiento del elemento
Reemplazando los desplazamientos nodales del elemento en (30), se obtiene que el campo de desplazamiento homogéneo para este elemento es:
|
(96) |
Mientras que debido a que la carga externa en el elemento tiene tres tramos, el campo de desplazamientos empotrado también los tendrá, es decir:
|
(97) |
donde, a partir de (42), el campo de desplazamiento empotrado en cada uno de los tres tramos se calcula como:
|
(98) |
|
(99) |
|
(100) |
Ahora sumando (96) y (97) se obtiene que el campo de desplazamiento total en el elemento es:
|
(101) |
donde:
|
(102) |
|
(103) |
|
(104) |
Cálculo del campo de desplazamiento del elemento
Procediendo de igual manera a como se hizo para el elemento , se calculan los campos de desplazamiento homogéneo, empotrado y total del elemento :
|
|
(106) |
Como resumen, en la Figura 11 se presenta el campo de desplazamientos en toda la viga.
Figura 11. Diagrama del campo de desplazamientos |
Cálculos de los campos de fuerzas internas
A partir de (101) y empleando (4), se obtiene que las fuerzas internas para el elemento son:
|
(107) |
|
(108) |
De forma similar, a partir de (106), se tiene que los campos de fuerzas internas en el elemento son:
|
Como resumen de los anteriores resultados, en la Figura 12 se presentan las fuerzas internas en toda la viga.
(a) Fuerza cortante | (b) Momento flector |
Figura 12. Diagramas de fuerzas internas |
Cálculo de las reacciones
Aunque es posible calcular las reacciones mediante 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:
|
(111) |
|
(112) |
En el siguiente enlace se encuentra el código de Python empleado para la solución de este ejemplo:
https://drive.google.com/open?id=1-rqbQquE5YmvSKxRtQW2z7Lpm1O8C8-Y
Resolver la viga sobre fundación flexible presentada en la Figura 13a.
(a) Problema a resolver | (b) Discretización |
Figura 13. 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 13b, la cual muestra una de las principales bondades de la presente metodología. Consistente 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 analítica.
Preliminares
El coeficiente para este problema es:
|
(113) |
Mientras que las funciones de forma son:
|
(114) |
|
(115) |
|
(116) |
|
(117) |
Por su parte la función de Green de la viga doblemente empotrada asociada con este problema se obtiene a partir de (65) y (66):
|
(118) |
donde:
|
Cálculo de las fuerzas de empotramiento
Para este problema la carga externa se expresa por la siguiente función que posee tres tramos:
|
(120) |
Con lo cual las fuerzas de empotramiento se calculan a partir de las funciones de forma presentadas en (114) a (117), como:
|
Formulación del sistema de ecuaciones de cada elemento
A partir de (84) se tiene que la formulación del método de rigidez para el único elemento de la discretización es:
|
(122) |
Obtención de los desplazamientos nodales
La solución de (122) es:
|
(123) |
Cálculo del campo de desplazamientos
Para el cálculo del campo de desplazamientos se realizará la siguiente superposición o descomposición:
|
(124) |
donde es el campo de desplazamientos homogéneo y es el campo de desplazamientos empotrado. A partir de (56) y empleando los desplazamientos y rotaciones nodales, se obtiene que el campo de desplazamiento homogéneo es:
|
(125) |
Mientras que el campo de desplazamiento empotrado se calcula a partir de (69) y se expresa como:
|
(126) |
donde:
|
(127) |
|
(128) |
|
(129) |
Con lo cual el campo de desplazamiento total es:
|
(130) |
|
(131) |
|
(132) |
|
(133) |
En la Figura 11 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:
|
(134) |
donde:
|
(135) |
|
(136) |
|
(137) |
Como resumen de las anteriores tres ecuaciones, en la Figura 14b se presenta la fuerza el suelo ejerce sobre la viga.
(a) Campo de desplazamientos | (b) Campo de la fuerza que el suelo ejerce sobre la viga |
Figura 14. Campos de desplazamiento en la viga y fuerza que el suelo ejerce sobre esta |
Como revisión de 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:
|
(139) |
donde:
|
(140) |
|
(141) |
|
(142) |
Mientras que el campo de fuerza cortante se calcula como:
|
(143) |
donde:
|
(144) |
|
(145) |
|
(146) |
En las Figuras 12a y 12b se presentan respectivamente los campos de fuerza cortante y momento flector para toda la viga.
(a) Fuerza cortante | (b) Momento flector |
Figura 15. Campos de fuerzas internas |
En el siguiente enlace 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 16a cuyos elementos son rectangulares de base y altura iguales a y cuyo módulo de elasticidad es .
(a) Problema a resolver | (b) Discretización |
Figura 16. Pórtico plano sometido a cargas externas oblicuas y discretización empleada para su solución |
Discretización
En la Figura 16b 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 16a, es posible obtener el valor de la carga externa en dirección del elemento :
|
(147) |
Mientras que la carga externa en dirección del eje en el elemento es:
|
(148) |
De forma similar, se tiene que las cargas externas en dirección de los ejes y en el elemento son respectivamente:
|
Cálculo del sistema de ecuaciones de cada elemento en coordenadas locales
A partir de lo presentado en (86) se tiene que las ecuaciones del método de rigidez en coordenadas locales para los elementos y son respectivamente:
|
(150) |
|
(151) |
Cálculo del sistema de ecuaciones de cada elemento en coordenadas globales
Las ecuaciones (150) y (151) se transforman a coordenadas globales como:
|
(152) |
|
(153) |
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 cinco ecuaciones con igual número de incógnitas:
|
(154) |
Cuya solución es:
|
(155) |
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:
|
(156) |
donde se emplearon los valores de los desplazamientos nodales presentados en (155).
Como revisión de las reacciones obtenidas anteriormente se realizará el equilibrio de toda la estructura:
|
Cálculo del campo de desplazamiento del elemento
Para poder calcular el campo de desplazamiento homogéneo del elemento es necesario haber transformado los desplazamientos en los extremos de este elemento de coordenadas globales a locales. Una vez realizado esto, se tiene que el campo de desplazamiento homogéneo en dirección del elemento es:
|
(158) |
Mientras que a partir de lo presentado en (22) y (147), se tiene que el campo de desplazamiento empotrado en dirección del elemento es:
|
(159) |
Con lo cual el campo total de desplazamiento en dirección del elemento es:
|
(160) |
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 del elemento son respectivamente:
|
(161) |
|
(162) |
|
(163) |
Cálculo del campo de desplazamiento del elemento
De igual forma a como se realizó con el elemento , se tiene que los campos de desplazamiento homogéneo, empotrado y total en dirección del elemento son respectivamente:
|
(164) |
|
(165) |
|
(166) |
Mientras que los campos de desplazamiento homogéneo, empotrado y total en dirección del elemento son respectivamente:
|
(167) |
|
(168) |
|
(169) |
Como síntesis de los campos de desplazamiento calculados anteriormente en la Figura 17d 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 de ambos elementos se obtiene las fuerzas internas en estos (ecuaciones (2) y (4)):
|
(170) |
|
(171) |
|
(172) |
|
Como resumen en las Figuras 17a a 17c se presentan los diagramas de fuerzas internas en los dos elementos que componen la estructura.
En el siguiente enlace 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 ha 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 autores. 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 (8) pero ahora a partir de la llamada forma débil de (8.a) en lugar de esta directamente. De forma similar a como se indicó para el método de rigidez con funciones de Green, el campo de desplazamientos se expresará como:
|
(174) |
donde y son las funciones de forma presentadas en (12) y se calcula según lo presentado en (22).
Además, es fácil probar que las funciones de forma y cumplen la siguiente propiedad asociada con (10.a) (ecuación diferencial del campo homogeneo):
|
(175) |
Multiplicando a ambos lados de (8.a) por una función de peso o ponderación se obtiene:
|
(176) |
Ahora, integrando a ambos lados de (176) con respecto a entre 0 y se obtiene:
|
(177) |
E integrando por partes en el lado izquierdo de la anterior ecuación el resultado es:
|
(178) |
La cual luego de evaluar el primer y reordenarse, da lugar a:
|
(179) |
Ahora si se emplean como las funciones de ponderación o peso () a las funciones de forma ( y ) se tiene que:
|
(180) |
Lo anterior debido a que y a (175).
Con lo cual (179) se puede reescribir como:
|
(181) |
Ahora, si como caso particular primero se emplea en lugar de y luego , a partir de (181) 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 (26).
Ahora, de forma similar a como se realizó en el apéndice A 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. (27) pero ahora empleando la forma de débil de la E.D. (27.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:
|
(183) |
donde una propiedad importante de las funciones de forma , , y es que estas cumplen la siguiente ecuación diferencial:
|
(184) |
Multiplicando por una función de peso o ponderación a ambos lados de (27.a) se obtiene:
|
(185) |
E integrando con respecto a entre 0 y a ambos lados de (185) se obtiene:
|
(186) |
Resolviendo con integración por partes la integral del lado izquierdo de la anterior ecuación se obtiene:
|
(187) |
Ahora, teniendo en cuenta que la fuerza cortante se obtiene a partir del campo de desplazamiento vertical como , se obtiene:
|
(188) |
Mientras que debido a que y , la anterior ecuación se reescribe como:
|
(189) |
Integrando por partes de nuevo, ahora en el tercer término del lado izquierdo de la ecuación anterior, se obtiene:
|
(190) |
Por último, teniendo en cuenta que el campo de momentos interno se calcula como y que y , la ecuación anterior se reescribe como:
|
(191) |
Reemplazando (183) en el lado derecho de (191), se obtiene:
|
(192) |
Ahora, si como funciones de ponderación se emplean las funciones de forma, empleando integración por partes se prueba lo siguiente:
|
(193) |
donde se ha empleado (184) y que los valores de y son iguales a cero en y .
Con lo cual, a partir de (193), la ecuación (192) se puede reescribir como:
|
(194) |
Si como caso particular primero se emplea , luego , y , a partir de (194) se obtienen las siguientes cuatro ecuaciones:
|
Que son equivalentes al sistema de ecuaciones presentado en (52).
Para este caso el objetivo es resolver el P.V.F. (53) pero ahora empleando la forma débil de (53.a) en lugar de esta directamente.
Como se presentó en (83), la solución de (53.a) se expresa como la suma de un campo homogéneo () y un campo empotrado (), es decir:
|
(196) |
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 , , y de la siguiente manera:
|
(197) |
Mientras que el campo empotrado se calcula solo empleando la carga externa y las funciones de Green.
Con lo cual, el campo de desplazamiento total, se expresa ahora como:
|
(198) |
Es importante resaltar que debido a que es una combinación lineal de las funciones de forma , , y , estas últimas cumplen la propiedad de ser solución de la ecuación diferencial homogénea asociada con (53.a), es decir:
|
(199) |
Multiplicando por una función de peso o ponderación a ambos lados de (5) (o (53.a)) se obtiene:
|
(200) |
E integrando con respecto a entre 0 y a ambos lados de (200) se obtiene:
|
(201) |
Resolviendo con integración por partes la integral del lado izquierdo de la anterior ecuación, da como resultado:
|
(202) |
Ahora teniendo en cuenta que la fuerza cortante se obtiene a partir del campo de desplazamiento vertical como , se obtiene:
|
(203) |
Teniendo en cuenta que y , la anterior ecuación toma la siguiente forma:
|
(204) |
Integrando por partes de nuevo, ahora en el tercer término del lado izquierdo de la ecuación anterior, se obtiene:
|
(205) |
Por último, teniendo en cuenta que el campo de momentos interno se calcula como y que y , la ecuación anterior se reescribe como:
|
(206) |
Reemplazando (198) en el lado derecho de (206), se obtiene:
|
(207) |
Ahora, si como funciones de ponderación se emplean las funciones de forma, empleando integración por partes se prueba lo siguiente:
|
(208) |
donde se ha empleado (199) y que los valores de y son iguales a cero en y .
Empleando (208) en (207) esta última toma la siguiente forma (donde para ser coherentes con (208), se ha tomado que ):
|
(209) |
Reemplazando por 2,3,5 y 6 en (209) se obtienen cuatro ecuaciones, las cuales se expresan en forma matricial como:
|
(210) |
donde:
|
(211) |
Lo cual es equivalente a lo presentado en (84).
[1] Challis L., Sheard F. The green of Green functions. American Institute of Physics, Physics Today, 56(12):41-46, 2003.
[2] Duffy D.G. Green’s functions with applications. Chapman and Hall/CRC, 2015.
[3] Banerjee P.K., Butterfield R. Boundary element methods in engineering science. McGraw-Hill, New York, 1981.
[4] Sánchez-Sesma F.J., Ramos-Martínez J., 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., Karageorghis A., 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 W.T. Transmission of elastic waves through a stratified solid medium. Journal of Applied Physics, 21(2):89-93, 1950.
[7] Boussinesq J. Application des potentiels a l'étude de l'équilibre et du mouvement des solides élastiques. Gauthier-Villars, Imprimeur-Libraire, 1885.
[8] Cerruti V. 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. AIP Physics, 7(5):195-202, 19336.
[10] Stokes G.G. On dynamical theory of diffraction. Transactions of the Cambridge Philosophical Society, 9:1-48, 1849.
[11] Lamb H. 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., Richards P.G. Quantitative seismology. Univ. Science Books, 700 pp., 2002.
[15] Colunga A.T. Análisis de estructuras con métodos matriciales. Limusa, 560 pp., 2007.
[16] McCormac J.C. Structural Analysis: using classical and matrix methods. Wiley Hoboken, NJ, 2007.
[17] Kassimali A. Matrix analysis of structures SI version. Cengage Learning, 640 pp., 2012.
[18] Reddy J.N. An introduction to the finite element method. McGraw-Hill, 2006.
[19] Bathe K.J. Finite element procedures. McGraw-Hill, 2006.
[20] Hetényi M. Beams on elastic foundation: theory with applications in the fields of civil and mechanical engineering. University of Michigan, 1971.
[21] Eisemberger M., Yankelevsky D.Z. 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?