Metodos Energeticos Investigacion

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 22

Universidad Central de Venezuela

Facultad de Ingeniería
Escuela de Ingeniería Mecánica
Depto. De Diseño
Mecánica de Solidos II

Métodos Energéticos

Prof. José De Freitas Alumno: Leonardo Castillo


C.I.: 23.696.858

Caracas, agosto de 2020


Introducción
Como bien es sabido en el basto mundo del diseño e ingeniería, es imprescindible
artilugios que nos permitan de alguna u otra forma “predecir” el futuro, por medio de
cálculos, para obtener la mayor seguridad posible al momento de consumar un
proyecto, de lo contrario pudieran ocurrir accidentes que incluso, atentaran contra la
vida humana, dichos artilugios se presentan en el presente trabajo, pues por medio de
los Teoremas de Castigiliano y del trabajo mínimo, podemos dar un claro avance al
momento de diseñar, pues estos métodos son versátiles y muy útiles al momento de
estimar deformaciones en vigas, estructuras, armaduras, Etc. tanto isostáticas como
hiperestáticas lo cual es una poderosa arma al momento de diseñar, nosotros como
ingenieros.
Como se dijo anteriormente, estos teoremas surgen de una necesidad, necesidad de
métodos alternativos pero versátiles para el ingeniero estimar fácilmente deflexiones o
cargas, sin necesidad de llevar a cabo complejos análisis, dando como resultado una
optimización del tiempo y esfuerzo.
Índice
Marco teórico .............................................................................. 1
1-Energía de deformación elástica............................................... 1
2-Cargas de Impacto y Súbitas ..................................................... 2
3- Influencia de la forma de un sólido en su eficiencia para el
soporte de cargas de impacto...................................................... 2
4- Primer teorema de Castigliano ................................................ 2
5- Segundo teorema de Castigliano ............................................. 4
6- Teorema del trabajo mínimo ................................................... 4
7-Distribución de esfuerzos normales en discos giratorios y en
cilindros de pared gruesa y espesor constante. ........................... 4
Ejercicio 1- 2do Teorema de Castigliano ..................................... 7
Ejercicio 2- Deformación Elástica ............................................... 11
Ejercicio 3- Teorema del trabajo mínimo ................................... 13
Aplicación laboral. ..................................................................... 17
Conclusión ................................................................................. 18
Bibliografía................................................................................. 19
Marco teórico

1-Energía de deformación elástica


La energía de deformación es el aumento de energía interna acumulado en el interior de un
sólido deformable como resultado del trabajo realizado por las fuerzas que provocan la
deformación

La energía de deformación es igual al trabajo realizado por una carga, la cual se incrementa
lentamente aplicada al elemento.
𝑥1
𝑈 = ∫ 𝑃𝑑𝑥
0

[𝑈] = 𝐽𝑂𝑈𝐿𝐸𝑆 𝑜 𝐹𝑡 ∗ 𝑙𝑏 𝐸𝑛 𝑐𝑎𝑠𝑜 𝑑𝑒 𝑡𝑟𝑎𝑏𝑎𝑗𝑎𝑟 𝑒𝑛 𝑈𝐼

Ahora bien, el valor de la energía de deformación U de un cuerpo sometido a esfuerzos normales


uniaxiales seria:
𝑥1
∇𝑥 2 𝑃2 𝐿
𝑈=∫ 𝑑𝑣 → 𝑈 =
0 2𝐸 2𝐴𝐸

Donde:

P= Carga a la cual se esta sometiendo el elemento.

A= área transversal del elemento.

E= Coeficiente de elasticidad.

Esta expresión es válida solo para deformaciones elásticas y se conoce como “energía de
deformación elástica de un cuerpo”

También se tiene que si un cuerpo es sometido a un esfuerzo cortante plano, su energía de


deformación se estima por la siguiente expresión:
𝑥1
ţ𝑥𝑌 2 𝑇 2𝐿
𝑈=∫ 𝑑𝑣 → 𝑈 =
0 2𝐺 2𝐴𝐺

Donde:

T= torque o momento cortante

A= Área transversal del elemento

G= módulo de rigidez

Nota: El módulo de rigidez se estima dada la expresión:


𝐸
𝐺=
2(1 + 𝑣)
Donde v= es el módulo de poisson.

1
Si el elemento se encuentra bajo un momento flector, el valor de energía de deformación
almacenada viene dado por:
𝑥1
𝑀2
𝑈=∫ 𝑑𝑣
0 2𝐸𝐼

Donde:

M=momento flector

E= módulo de elasticidad

I= momento de inercia

2-Cargas de Impacto y Súbitas


Se le conoce como carga de impacto al efecto dinámico sobre una estructura, móvil o estática,
que tiene una carga aplicada de corta duración debido a su movimiento. También llamada carga
móvil. Un ejemplo muy común es el choque de dos objetos, estos sufrirán el efecto de la carga
por breve periodo de tiempo, muy distintas a las cargas “fluctuantes” cuyos efectos persisten a
lo largo del tiempo y su magnitud varia constantemente.

3- Influencia de la forma de un sólido en su eficiencia para el


soporte de cargas de impacto
Si se considera solamente la forma, un sólido puede soportar cargas de impacto con mayor
eficiencia mientras no tenga concentradores de esfuerzos, es decir curvas muy pronunciados,
por ejemplo las muescas son concentradores de esfuerzos debido a que poseen una menor área
los esfuerzos son mayores, en consecuencia, el material puede fracturarse con mayor facilidad
al ser sometido a una carga de impacto. Por el contrario, si la forma de la pieza o el sólido a
someter al ensayo de impacto tiene bordes mas suavizados, tiene mayor capacidad de absorber
energía en forma de deformación elástica sin llegar a fracturarse y/o deformarse plásticamente.

4- Primer teorema de Castigliano


Un método alternativo, no tan mencionado e incluso mas optimizado para analizar la deflexión
es posible mediante el llamado teorema de Castigliano. Esta herramienta constituye una forma
única de analizar deflexiones e incluso es útil para encontrar las reacciones de estructuras
indeterminadas sin necesidad de llevar a cabo extensos análisis.

El teorema de Castigliano estipula que cuando actúan fuerzas sobre sistemas elásticos sujetos
a desplazamientos pequeños, el desplazamiento correspondiente a cualquier fuerza, colineal
con la fuerza, es igual a la derivada parcial de la energía de deformación total con respecto a esa
fuerza.

2
Para determinar deformaciones (desplazamientos) en un cuerpo debido a la aplicación de
fuerzas externas. Considérese un cuerpo en equilibrio estático sometido a un conjunto de
fuerzas puntuales. Donde UT es la energía total acumulada por el cuerpo debido a todas las
cargas externas, podemos ver que:

UT = UT (F1, F2, ..., Fi, ..., Fn)

-es decir UT es función de todas estas fuerzas externas-

Ahora bien, asumamos ahora que la fuerza Fi sufre un incremento muy pequeño en su magnitud
dFi, luego usando una aproximación con series de Taylor tenemos
𝜕𝑈
UT (F1, F2, . . . , Fi + dFi, . . . , Fn) ≈ UT (F1, F2, . . . , Fi, . . . , Fn) + 𝜕𝐹𝑡 𝑑𝐹𝑖 (1)
𝑖

𝜕𝑈𝑡
Donde , está siendo evaluado con (F1, F2, ..., Fi, ..., Fn).
𝜕𝐹𝑖

En un cuerpo elástico la energía acumulada por el cuerpo es igual al trabajo externo hecho sobre
él, de modo que el trabajo inicial (antes de incrementar Fi) es igual a UT (F1, F2, ..., Fi, ..., Fn).

Por otra parte, si Fi se incrementa en dFi, podemos calcular el incremento en el trabajo de la


siguiente forma, si en un cuerpo en donde particularmente podemos apreciar la fuerza Fi, la
que en ese punto de aplicación genera un desplazamiento δi. Si la fuerza se incrementa en dFi
el incremento en trabajo en el punto de aplicación de Fi mecánico es aproximadamente igual a:

𝛿𝑖 𝑑𝐹𝑖
El trabajo final total al incrementar la fuerza es entonces igual a:

𝑈𝑇 (F1, F2, . . . , Fi, . . . , Fn) + 𝛿𝑖 𝑑𝐹𝑖 (2)

Como el trabajo final es igual a la energía UT (F1, F2, ..., Fi + dFi, ..., Fn) (debido a que el cuerpo
es elástico) de (1) y (2) tenemos:
𝜕𝑈𝑡
dFi = 𝛿𝑖 𝑑𝐹𝑖 .
𝜕𝐹𝑖

Donde finalmente tenemos lo que se conoce como Teorema de Castigliano

𝜕𝑈𝑡
𝛿𝑖 =
𝜕𝐹𝑖
En el caso de desplazamiento rotacional, la ecuación puede escribirse como

∅ 𝜕𝑈
𝑖=
𝜕𝑀𝑖

3
5- Segundo teorema de Castigliano
Está dirigido para resolver estructuras hiperestáticas e isostáticas, donde se conoce
previamente, que su desplazamiento o su rotación es igual a 0.
𝑥1
𝜕𝑀𝑧 𝑑𝑠
∇𝑦 = ∫ 𝑀𝑧
0 𝜕𝑃 𝐸𝐼

Para un problema estáticamente indeterminado, el equilibrio aun constituye la base del


método usado, pero el número de incógnitas estáticas excedería el número de
disponibles de ecuaciones de equilibrio, como se vio en el primer teorema, estas
incógnitas adicionales, se denominan Redundantes. Una vez que se selecciona la o las
incógnitas redundantes, las ecuaciones de equilibrio se resuelven para las incógnitas
restantes en función de las redundantes y/o ficticias en cuyo caso que se quiera saber
el valor de la deflexión en un punto en específico.

6- Teorema del trabajo mínimo


Este teorema es consecuencia del teorema de Castigliano, en los apoyos donde se
cumpla que di=0,es decir, que los desplazamientos lineales sean nulos o los giros, ya sea
por empotramiento, apoyos fijos, etc.
𝑥1
𝜕𝑢 𝜕𝑀𝑧 𝑑𝑠
=∫ 𝑀𝑧
𝜕𝑅 0 𝜕𝑃 𝐸𝐼
Donde la derivada parcial del desplazamiento respecto a las reacciones es igual a la
derivada parcial del momento en z respecto a la carga, multiplicada a su vez por el
mismo momento en z por un diferencial de sección del elemento dividido entre el
modelo de elasticidad y el momento de inercia, como se ilustra en la formula.
Este teorema es beneficioso cuando se desea conocer las reacciones en una viga
estáticamente indeterminada por ejemplo, pues se conocen varios puntos donde el
desplazamiento vertical es 0 e incluso el angular en caso haber empotramiento,
facilitando asi la expresión obteniendo:
𝑥1
𝜕𝑀𝑧 𝑑𝑠
∫ 𝑀𝑧 =0
0 𝜕𝑃 𝐸𝐼

7-Distribución deesfuerzos normales en discos giratorios y en


cilindros de pared gruesa y espesor constante.
La naturaleza de las fuerzas que surgen en el tubo de paredes gruesas bajo la acción de
una presión y las que se originan en un disco que gira a gran velocidad es distinta. Sin
embargo ek cakcuki de estas dos piezas conduce a un esquema de calculo común el

4
esquema de “cuerpo de revolución. El análisis posterior demuestra también que son
idénticas las ecuaciones diferenciales que determinan los desplazamientos y las
tensiones (esfuerzos) en los dos casos. Es conveniente, por tanto, análisis ambos casos
simultáneamente.
Cada punto del cilindro al deformarse este, se desplaza, de las condiciones de simetría
se deduce que estos desplazamientos ocurren en los planos radiales. El punto puede
desplazarse en la dirección del radio y a lo largo de la generatriz correspondiente
Ahora bien, en las secciones axiales del cilindro (plano ABCD), debido a las condiciones
de simetría axial, son nulas. Las tensiones son inversamente proporcionales al cuadrado
del radio r, entonces en los puntos situados a esta distancia del eje, las tensiones
constituirán solamente 1/16 parte de las tensiones máximas. Es por ello que es muy
importante guir la variación de las tensiones y a medida que disminuye el espesor del
cilindro, admitamos que b=a+δ, siendo δ el espesor del cilindro. Entonces:

(𝑎 + 𝛿)2 + 𝑎2
𝜎𝑡 = 𝑝
𝛿(2𝑎 + 𝛿)
Cuando δ es pequeño,
𝑎
𝜎𝑡 ≈ 𝑝
𝛿
Las tensiones normales en las secciones transversales de un cilindro de paredes gruesas
pueden aparecer como consecuencia de la solicitación del cilindro por cargas orientadas
según el eje. Se supone que estas tensiones no varian tanto en dirección axial, como
radial del cilindro, la distribución seria:

5
Sección transversal

Incidencia de las cargas en un Diferencial infinitesimal

6
Ejercicio 1- 2do Teorema de Castigliano
Se tiene una viga empotrada en voladizo con una carga distribuida uniforme parcial en
el lado empotrado de 4kN/m. Se pide la deflexión en el punto “B” por el 2do teorema
de Castigliano y además su respectiva comprobación por otro método alternativo ya
estudiado al inicio del curso.

Sabemos que para conocer la deflexión en un punto cualquiera de la viga es necesario


ubicar sobre este una carga “ficticia” en el punto de interés que queramos conocer, con
el fin de derivar parcialmente en función de esa carga ficticia más adelante en el
problema como consecuencia del Teorema.
Como queremos conocer la deflexión en el punto B, procedemos a ubicar la carga ficticia
allí, como se muestra a continuación:

Procedemos a hacer el respectivo DCL de la viga para conocer sus reacciones y


expresarlas en función de nuestra carga fictcia:

7
5
∑ 𝑀𝑎 = 0 → − ∗ 20 − 5𝑃 + 𝑀𝑎 = 0
2
𝑀𝑎 = 5𝑃 + 50 (1)

∑ 𝐹𝑦 = 0 → 𝑅𝑎 − 𝑃 − 20 = 0

𝑅𝑎 = 𝑃 + 20 (2)

Ahora necesitamos saber los momentos internos para ello haremos 2 cortes en la viga
para determinar el valor de esos momentos:

8
Tramo 1-1 (0≤x≤5)

𝑥
∑ 𝑀1 = 0 → 𝑀1 + (4𝑥) + 𝑀𝑎 − 𝑥𝑅𝑎 = 0
2
𝑀1 = −2𝑥 2 − 𝑀𝑎 − 𝑥𝑅𝑎 (3)
Sustituyendo (1) y (2) en (3) tenemos
𝑀1 = −2𝑥 2 − (5𝑃 + 50) + 𝑥(20 + 𝑃) (4)
Derivando parcialmente con respecto a la carga ficticia se tiene:
𝜕𝑀1
= −5 + 𝑥 (5)
𝜕𝑃

Análogamente para el tramo 2-2 (5≤x≤10) se tiene:

∑ 𝑀2 = 0 → 𝑀2 + (𝑥 − 5)𝑃 + (𝑥 − 2,5)20 + 𝑀𝑎 − 𝑥𝑅𝑎 = 0

9
𝑀2 = −(𝑥 − 5)𝑃 − (𝑥 − 2,5)20 − (5𝑃 + 50) + 𝑥(20 + 𝑃)
𝜕𝑀2
= −5 + 𝑥 − (𝑥 − 5)
𝜕𝑃
Ahora bien por 2do teorema de castigliano tenemos:
𝑥1
𝜕𝑀𝑧 𝑑𝑠
∇𝑦 = ∫ 𝑀𝑧
0 𝜕𝑃 𝐸𝐼
Por consiguiente para saber la deflexión en el punto B, basta con sustituir las ecuaciones
(4) y (5) en el teorema teniendo en cuenta que la carga P es “Ficticia” por lo cual no se
toma en cuenta al momento de sustituir:

5
1
∇𝑦 𝐵 = [∫ (−2𝑥 2 − (50) + 𝑥(20))(−5 + 𝑥) dx ]
𝐸𝐼 0
5
1
= [∫ (−2𝑥 3 − 50𝑥 + 20𝑥 2 + 10𝑥 2 + 250 − 100𝑥) dx ]
𝐸𝐼 0
5
1
= [∫ (−2𝑥 3 + 30𝑥 2 − 150𝑥 + 250) dx ]
𝐸𝐼 0

1 −2𝑥 4 30𝑥 3 150𝑥 2


= [ + − + 250𝑥]
𝐸𝐼 4 3 2

Evaluado entre 0 y 5 se tiene que la deflexión en B es:


312,5
∇𝑦 𝐵 =
𝐸𝐼
Nota: como el valor es positivo quiere decir que el valor que escogimos de la deflexión,
hacia abajo, es correcto.
Ahora bien, para comprobar que efectivamente esa es la expresión correcta de la
deflexión en B, procedemos a analizar la misma viga por el teorema de superposición,
gracias a las tablas podemos observar que la expresión para la deflexión de una viga en
voladizo con carga uniforme parcial en el lado del empotrado se rige por:

10
𝑦𝑎𝑏=−𝑊𝑜(6𝑎2 𝑥 2 −4𝑎𝑥 3 +𝑥 4 )
24𝐸𝐼

Para nuestro problema los valores son los siguientes:


Wo=4; a=5; x=5

𝑦𝑎𝑏= −4 2 2 3 4
24𝐸𝐼 (6∗5 ∗5 −4∗5∗5 +5 )

𝑦 −312,5
𝑏= 𝐸𝐼
Nota: el resultado es negativo pues, la deflexión va hacia abajo.
Al sustituir esos valores en la formula obtenemos la misma expresión conseguida por el
2do teorema de castigliano, confirmando así la veracidad del procedimiento.

Ejercicio 2- Deformación Elástica


A continuación se presenta un ejercicio simple, de una viga apoyada a ambos extremos
con apoyos simples y una carga puntual de 10Kn en el medio de la viga, con el fin de
calcular su deformación elástica, aplicaremos el teorema de doble integración.

Haciendo cortes antes y después de la aplicación de la fuerza puntual, tenemos los


siguientes diagramas:

11
Analizando los diagramas sacamos la siguiente ecuación de momento:
𝑀 = 5𝑥 − (𝑥 − 4)10
Ahora bien por teorema de doble integración se tiene que:

𝑑2𝑦
𝐸𝐼 =𝑀
𝑑𝑥 2
𝑑𝑦
𝐸𝐼 = ∫ 5𝑥 − (𝑥 − 4)10
𝑑𝑥
𝑑𝑦 5𝑥 2 (𝑥 − 4)2
𝐸𝐼 = − ∗ 10 + 𝑐1
𝑑𝑥 2 2
5𝑥 3 (𝑥 − 4)3
𝐸𝐼𝑦 = − ∗ 10 + 𝑐1𝑥 + 𝑐2
6 6

12
Para saber cuales son las constantes de integración, apelamos por las condiciones de
contornos.
Para X=0; Y=0 C2=-106.66
Para x=8; y=0 C1=-26.66
5 ∗ 83 (8 − 4)3
0= − ∗ 10 + 𝑐1 ∗ 8 − 106.66
6 6
Por consiguiente tenemos pues la ecuación completa para la deflexión:

5𝑥 3 (𝑥 − 4)3
𝐸𝐼𝑦 = − ∗ 10 − 26.66𝑥 − 106.66
6 6
Para saber la deformación máxima que sufrirá la viga, basta con sustituir x=4, pues es
allí donde se producirá la max deformación, o curvatura de la viga.
1 5𝑥 3 (𝑥 − 4)3
𝑦= [ − ∗ 10 − 26.66𝑥 − 106.66]
𝐸𝐼 6 6

1 5(4)3 (4 − 4)3
𝑦= [ − ∗ 10 − 26.66(4) − 106.66]
𝐸𝐼 6 6

−𝟏𝟔𝟎
𝒚=
𝑬𝑰

Ejercicio 3- Teorema del trabajo mínimo


A continuación se presenta una viga estáticamente indeterminada, empotrada en el
extremo izquierdo y con un apoyo simple en el extremo derecha, sobre ella se aplica una
carga distribuida triangular de 4 kN/m, la longitud total de la viga es de 6m. Se pide Se
requiere saber el valor de las reacciones utilizando el Teorema del Trabajo mínimo.

13
Como se puede ver, estamos en presencia de un caso de una viga hiperestático o como
también se le conoce “estáticamente indeterminado”
Se determina le grado de indeterminación restando los grados de libertad con el numero
de incógnitas, tenemos en total 3 incógnitas, dos por empotramiento y una por apoyo
simple, sabiendo que tenemos 4GDL se tiene:
H=4-3=1
Esto quiere decir que necesitamos una ecuación más para poder determinar todo el
sistema.
Como dicta el Teorema del trabajo mínimo debemos escoger una incógnita
“Redundante” en este caso, por practicidad, será la reacción en B

Al realizar el corte tomando el lado de la viga mas conveniente, en este caso el lado
derecho, obtenemos:

14
Como podemos observar al tratar de puntualizar la carga triangular para poder efectuar
la suma de momentos, nos encontramos con que no sabemos a que altura en función
de x se encuentra, para ello utilizamos Semejanza de triángulos:

4 𝑦 2
= ↔𝑦= 𝑥
6 𝑥 3
𝑥2 𝑥
∑ 𝑀𝑧 = 0 → ( ) (− ) − 𝑀𝑧 + 𝑥𝑅𝑏 = 0
3 3

𝑥3
𝑀𝑧 = (− ) + 𝑥𝑅𝑏 (1)
9

Ahora derivamos parcialmente con respecto a la incógnita redundante que escogimos


en primera instancia (Rb)
𝜕𝑀𝑧
= 𝑥 (2)
𝜕𝑅𝑏

Ahora por teorema de trabajo mínimo tenemos:


𝑥1
𝜕𝑀𝑧 𝑑𝑠
∫ 𝑀𝑧 =0
0 𝜕𝑃 𝐸𝐼

Pues sabemos que la deflexión en B es igual a 0. Además sustituyendo (1) y (2) en el


teorema se tiene:

6
𝑥3 𝑑𝑥
∫ ((− ) + 𝑥𝑅𝑏)(𝑥) =0
0 9 𝐸𝐼

1 6 𝑥4
∫ ((− ) + 𝑅𝑏 ∗ 𝑥 2 )dx = 0
𝐸𝐼 0 9

15
1 𝑥5 𝑥 3 𝑅𝑏
[(− ) + ]=0
𝐸𝐼 45 3

Evaluando de 0 a 6 tenemos:

63 𝑅𝑏 65
− = 0 → 𝑅𝑏 = 2,4𝑘𝑁
3 45

Para obtener las otras incógnitas se recurre a las ecuaciones estáticas de equilibrio

DCL

∑ 𝑀𝑎 = 0 → −2(24) + 6(2,4) − 𝑀𝑎 = 0

𝑀𝑎 = −33,6𝑁𝑚
Da negativo, lo que quiere decir que la dirección es contraria a la escogida en primera
instancia es decir, la dirección correcta es en sentido anti horario.

∑ 𝐹𝑦 = 0 → 𝑅𝑎 + 𝑅𝑏 − 24 = 0

16
𝑅𝑎 = 21,6𝑘𝑁

Aplicación laboral.
Estos métodos energéticos alternativos son herramienta fundamental de los ingenieros
civiles al verse envueltos en armaduras y estructuras todo el tiempo este método puede
darle una precisa estimación en cuanto a cálculos de deflexiones y cargas se refiere al
momento de analizar vigas, estructuras y edificios. Otra aplicación que puede tener
además es en la construcción de “andamios” como bien sabemos, es la armadura más
común que se puede observar en el día a día, detrás de esta armadura hay cálculos que
determinan pues que cargas son permisibles de modo de garantizar la seguridad de cada
pieza para que esta no se deforme o falle y pueda poner en peligro la vida de su
operador, estos métodos energéticos pueden ayudar con la determinación de esas
cargas.

17
Conclusión
Como se pudo observar a lo largo del proyecto, los teoremas de Castigiliano y el
teorema del trabajo mínimo, demostraron ser muy versátiles al momento de diseñar y
estimar cargas y deformaciones desconocidas, pero también es conocido que estos
teoremas tienen un punto limitante y es que si no se conoce previamente, por ejemplo
en el caso del teorema del trabajo mínimo, la deformación vertical en algún punto, se
complicarían los cálculos al punto de seguir teniendo más incógnitas que ecuaciones lo
cual sería contraproducente, además en caso de tener múltiples cargas solicitadas en la
viga y múltiples reacciones, se tendría que trabajar con varios sistemas “Redundantes”
al mismo tiempo lo cual se convertiría en un artilugio poco práctico.

18
Bibliografía
Mott,R. (Ed).(2009). Resistencia de Materiales. México: Editorial: PEARSON EDUCATION.

Bickford, W.(Ed).(). Mecánica de sólidos. España: Editorial: McGraw-Hill/Irwin

Shigley, J. & C. Mischke. 2002. Diseño en Ingeniería Mecánica. McGraw-Hill/Irwin

Bustamante Plaza, R. Mecánica de Solidos.

19

También podría gustarte