Taller2 Ecdi Grupo4

TALLER – ECUACIONES DIFERENCIALES-4
PRESENTADO POR:
BRAYAN FELIPE SALAZAR QUINTANA
MARIA CAMILA HERNANDEZ MORENO
JUAN SEBASTIAN SANCHEZ PATARROYO
ARLEX SANTIAGO FARFAN CARDENAS
PRESENTADO A: LIDA BUITRAGO GARCIA
UNIVERSIDAD: ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO

BOGOTÁ. D, C
2023
Para realizar el análisis de flexión de una placa circular, se considera a la placa en un sistema
coordenado polar, de tal manera que la función de flexión w(r, φ) se puede obtener al establecer el
equilibrio sobre un elemento diferencial en el sistema cilíndrico (r, φ, z) , considerando los esfuerzos,
momentos y cortantes referidos al sistema coordenado y las relaciones de equilibrio adecuadas. En
el análisis de la flexión de una placa circular uniformemente cargada, w(r) depende solo de la
distancia radial r desde cualquier punto de la placa a su centro, dando como resultado que la curva
de flexión de la placa está determinada por la ecuación.
(1)
Donde q y D son constantes
1. Encuentre la solución general de la ecuación, de manera analítica. Verifique su

resultado usando un asistente computacional.
Si se multiplica la ecuación por 𝑟 3 queda:
3
𝑑3 𝑤 2
𝑑2𝑤 𝑑𝑤 𝑞 4
𝑟 + 𝑟 − 𝑟 = 𝑟
𝑑𝑟 3 𝑑𝑟 2 𝑑𝑟 2𝐷
Ahora la ecuación tiene forma a una de Cauchy-Euler, así que resolviendo la parte
homogénea:
𝑑3𝑤 𝑑2𝑤 𝑑𝑤
𝑟3 3 + 𝑟2 2 − 𝑟 =0
𝑑𝑟 𝑑𝑟 𝑑𝑟
Suponiendo que una solución es 𝑤 = 𝑟 𝑚 , se reemplaza en la ecuación
𝑟 3 ∗ (𝑚(𝑚 − 1)(𝑚 − 2)𝑟 𝑚−3 ) + 𝑟 2 ∗ (𝑚(𝑚 − 1)𝑟 𝑚−2 ) − 𝑟 ∗ (𝑚𝑟 𝑚−1 ) = 0
𝑟 𝑚 (𝑚3 − 3𝑚2 + 2𝑚 + 𝑚2 − 𝑚 − 𝑚) = 0
𝑟 𝑚 (𝑚2 (𝑚 − 2)) = 0
Por lo tanto, las raíces son: m=0 (multiplicidad 2), m=2. Así que la solución de la parte
homogénea es:
𝑤ℎ = 𝑐1 + 𝑐2 ln 𝑟 + 𝑐3 𝑟 2
Para la parte no homogénea se encuentra una solución particular por el método de variación
de parámetros.
Suponiendo que la solución particular tiene forma: 𝑤𝑝 = 𝑤1 𝑢1 + 𝑤2 𝑢2 + 𝑤3 𝑢3
Donde
𝑊1 𝑊2 𝑊3
𝑢′1 = , 𝑢′2 = , 𝑢′3 =
𝑊 𝑊 𝑊
𝑞
Y 𝑓(𝑟) = 2𝐷 𝑟, ya que se usa la formula original para este método
1 ln 𝑟 𝑟2
1
| 0 2𝑟| 4
𝑊=| 𝑟 |=𝑟
1
0 2
𝑟2
0 ln 𝑟 𝑟2
1 2 2
0 2𝑟| 2𝑞𝑟 ∗ ln 𝑟 − 𝑞𝑟
𝑊1 = || 𝑟 | =
𝑞 1 2𝐷
𝑟 2
2𝐷 𝑟2
1 0 𝑟2 2
𝑊2 = |0 0 2𝑟| = − 𝑟 𝑞
𝑞 𝐷
0 𝑟 2
2𝐷
1 ln 𝑟 0
1
0 0 | 𝑞
𝑊3 = || 𝑟 =
| 2𝐷
1 𝑞
0 𝑟
𝑟2 2𝐷
2𝑞𝑟 2 ∗ ln 𝑟 − 𝑞𝑟 2
( ) 2𝑞𝑟 3 ∗ ln 𝑟 − 𝑞𝑟 3
2𝐷
𝑢 ′1 = =
4 8𝐷
(𝑟 )
𝑟 2𝑞
(− 𝐷 ) 𝑟 3𝑞
𝑢′ 2 = =−
4 4𝐷
(𝑟 )
𝑞
(2𝐷 ) 𝑞𝑟
𝑢′3 = =
4 8𝐷
(𝑟 )
Integrando cada expresión:

2𝑞 3
𝑞 3
𝑞𝑟 4 ∗ ln 𝑟 3𝑞𝑟 4
𝑢′1 = ∫ 𝑢′1 𝑑𝑟 = ∫ 𝑟 ∗ ln 𝑟 𝑑𝑟 − ∫ 𝑟 𝑑𝑟 = −
8𝐷 8𝐷 16𝐷 64𝐷
′ ′
𝑞 3
𝑞𝑟 4
𝑢 2 = ∫ 𝑢 2 𝑑𝑟 = − ∫ 𝑟 𝑑𝑟 = −
4𝐷 16𝐷
𝑞 𝑞𝑟 2
𝑢′3 = ∫ 𝑢′3 𝑑𝑟 = ∫ 𝑟 𝑑𝑟 =
8𝐷 16𝐷
La solución particular es:
𝑞𝑟 4 ∗ ln 𝑟 3𝑞𝑟 4 𝑞𝑟 4 2
𝑞𝑟 2
𝑤𝑝 = ( − ) + (ln 𝑟 ∗ − ) + (𝑟 ∗ )
16𝐷 64𝐷 16𝐷 16𝐷
𝑞𝑟 4
𝑤𝑝 =
64𝐷
Finalmente, la solución general de la ecuación es:
𝑤(𝑟) = 𝑤ℎ + 𝑤𝑝
𝑞𝑟 4 2
𝑤(𝑟) = 𝑐1 + 𝑐2 ln 𝑟 + 𝑐3 𝑟 + (3)
64𝐷
Verificando con Wolfram Mathematica:
Por lo que, la solución encontrada es correcta

2. Determine una solución que satisfaga las condiciones de frontera:
w′(0) = 0, w(a) = 0 y w′(a) = 0
donde a > 0 es el radio de la placa. [Sugerencia: La condición w′(0) = 0 es correcta]
1 𝑞𝑟 3
𝑤 ′ (𝑟) = 𝑐2 ∗ + 2𝑐3 𝑟 +
𝑟 16𝐷
A simple vista podemos ver que tanto para w(0) y w’(0) la solución se indetermina, pero
esto no puede ser así, debido a que, la placa está siendo uniformemente cargada, así que,
hasta en el radio 0 (en el origen) aunque no se flexione, no se puede descartar ya que, ahí
hay carga. En conclusión, 𝑐2 = 0.
𝑞𝑎4
𝑤(𝑎) = 𝑐1 + 𝑐3 𝑎2 + =0
64𝐷
′ (𝑎)
𝑞𝑎3 𝑞𝑎2
𝑤 = 2𝑐3 𝑎 + =0 → 𝑐3 = −
16𝐷 32𝐷
Reemplazando
𝑞𝑎2 2 𝑞𝑎4
𝑤(𝑎) = 𝑐1 + (− )𝑎 + =0
32𝐷 64𝐷
𝑎4 𝑞
𝑐1 =
64𝐷
La solución queda:
𝑎4 𝑞 𝑞𝑎2 2 𝑞𝑟 4
𝑤(𝑟) = + (− )𝑟 +
64𝐷 32𝐷 64𝐷
𝑎4 𝑞 + 𝑞𝑟 4 − 2𝑞𝑎2 𝑟 2 𝑞(𝑎2 − 𝑟 2 )2
𝑤(𝑟) = =
64𝐷 64𝐷
3. Use un asistente computacional para resolver el problema de valores en la

frontera para diferentes valores de a.
Tomando los valores a= 10,20,30,40,50

a=10
a=20
a=30
a=40
a=50
4. Grafique las diferentes curvas solución para los valores en la frontera y escriba
una breve conclusión.
En esta gráfica se puede observar como a medida que aumentamos el valor de “a”
radio de la placa) la flexión va a ser mayor en el centro de la placa (radio cero) y va
a ser menor en el borde de la placa (radio igual al de la placa). Los valores de q y D
se tomaron como 1 ya que son constantes proporcionales y la idea es ver el
comportamiento entre las gráficas modificando a “a”.
5. Verifique que la ecuación (1) puede escribirse en la forma alternativa

𝑑 1𝑑 𝑑𝑤 𝑞
[ (𝑟 )] = 𝑟 (2)
𝑑𝑟 𝑟 𝑑𝑟 𝑑𝑟 2𝐷
Empezamos desde la ecuación (2) para ver si es una forma alternativa de (1)
[ (𝑟 )]= 𝑟
Realizamos la derivada que está entre paréntesis
𝑑 1 𝑑𝑤 𝑑2𝑤 𝑞
[ ( + 𝑟 2) ] = 𝑟
Multiplicamos ese resultado por 1/r
𝑑 1 𝑑𝑤 𝑑 2 𝑤 𝑞
[ + 2 ]= 𝑟
Realizamos la última derivada de esa suma.
1 𝑑𝑤 1 𝑑 2 𝑤 𝑑 3 𝑤 𝑞
− 2
+ + 2 = 𝑟 (1)
𝑟 𝑑𝑟 𝑟 𝑑𝑟 𝑑𝑟 2𝐷
Concluimos que (2) si es una forma de escribir la ecuación (1)
6. Resuelva la ecuación (2) usando sólo integración con respecto a r y compruebe

que el resultado obtenido es equivalente a la solución general obtenida
inicialmente.
[ (𝑟 )]= 𝑟
Integramos respecto a r en ambos lados
∫ [ (𝑟 ) ] 𝑑𝑟 = ∫ 𝑟 𝑑𝑟
1𝑑 𝑑𝑤 𝑞 2
[ (𝑟 )] = 𝑟 + 𝑐1
𝑟 𝑑𝑟 𝑑𝑟 4𝐷
Multiplicamos r en ambos lados de la igualdad
𝑑 𝑑𝑤 𝑞 3
[ (𝑟 )] = 𝑟 + 𝑟𝑐1
𝑑𝑟 𝑑𝑟 4𝐷
𝑑 𝑑𝑤 𝑞 3
∫[ (𝑟 ) ] 𝑑𝑟 = ∫ 𝑟 𝑑𝑟 + ∫ 𝑟𝑐1 𝑑𝑟
𝑑𝑟 𝑑𝑟 4𝐷
𝑑𝑤 𝑞 4 𝑟2
(𝑟 )= 𝑟 + 𝑐1 + 𝑐2
𝑑𝑟 16𝐷 2
𝑑𝑤 𝑟2 𝑞 4
(𝑟 ) − 𝑐1 + 𝑐2 = 𝑟
𝑑𝑟 2 16𝐷
Multiplicamos 1/r en ambos lados de la igualdad
𝑑𝑤 𝑟 𝑐2 𝑞 3
( ) − 𝑐1 + = 𝑟
𝑑𝑟 2 𝑟 16𝐷
𝑑𝑤 𝑟 𝑐2 𝑞 3
∫( ) 𝑑𝑟 − ∫ 𝑐1 𝑑𝑟 + ∫ 𝑑𝑟 = ∫ 𝑟 𝑑𝑟
𝑑𝑟 2 𝑟 16𝐷
𝑟2 𝑞 4
𝑤 − 𝑐1 + c2ln[𝑟] + 𝑐3 = 𝑟
4 64𝐷
𝑟2 𝑞 4
𝑤 = 𝑐3 + 𝑐2 ln[𝑟] − 𝑐1 + 𝑟
4 64𝐷
Como c1 es constante podemos reducir el coeficiente de -1/4 con c1
𝑞 4
𝑤 = 𝑐3 + 𝑐2 ln[𝑟] + 𝑐1𝑟 2 + 𝑟
64𝐷
Podemos observar que esta ecuación es equivalente a la obtenida en (3) dando a entender
que otra forma de resolver la ED (1) es reescribirla como (2) para luego integrar y obtener
el mismo resultado que si lo hiciéramos por otro método (en este caso Cauchy-Euler).

Taller2 Ecdi Grupo4

Cargado por

Copyright:

Formatos disponibles

Taller2 Ecdi Grupo4

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Taller2 Ecdi Grupo4

Cargado por

Copyright:

Formatos disponibles

TALLER – ECUACIONES DIFERENCIALES-4

PRESENTADO A: LIDA BUITRAGO GARCIA

UNIVERSIDAD: ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO

1. Encuentre la solución general de la ecuación, de manera analítica. Verifique su

Integrando cada expresión:

Por lo que, la solución encontrada es correcta

3. Use un asistente computacional para resolver el problema de valores en la

Tomando los valores a= 10,20,30,40,50

5. Verifique que la ecuación (1) puede escribirse en la forma alternativa

6. Resuelva la ecuación (2) usando sólo integración con respecto a r y compruebe

También podría gustarte