Teorema de Castigliano.

En este tema se limita las consideraciones a deformaciones pequeñas en materiales

linealmente elásticos, entonces el resultado de la relación carga – deformación para
cualquier miembro o estructura es lineal y en estas condiciones la energía de deformación
U es igual a la energía complementaria C y el principio de superposición tiene aplicabilidad.
Entonces:

“si se apoya un sistema elástico para que se eviten los desplazamientos de cuerpo rígido del
sistema, y si ciertas fuerzas concentradas de magnitudes F1, F2, …. P actúan sobre el
sistema, además de las cargas distribuidas y las tensiones térmicas, el componente de
desplazamiento Δi; del punto de aplicación de la fuerza Fi, viene determinado por la
ecuación.”

𝜕𝐶
∆𝑖 = , 𝑖 = 1,2,3 ….
𝜕𝐹𝑖

C= área sobre la curva = energía complementaria.

U= área bajo la curva = energía de deformación.

dC=e*dp

e= dC/dP

Para el caso de comportamiento lineal C=U

∂U
∆i = ∂Fi, i=1,2,,….n (5.4)

∂U
𝜃i = ∂Mi, i=1,2,,….n (5.5)

FIGURE 5.1 Nonlinear elastic load-elongation curve. (Boresi, 2003)

Donde U = U(F1, F2,…..M1, M2,…..)

La energía de deformación U es:

𝑈 = ∫ 𝑈0 𝑑𝑣 (5.6)

Donde U0 es la densidad de energía de deformación, la cual para el caso de materiales con

comportamiento linear elástico, isotrópicos y homogéneos se expresa como:
 1 2 2 2 𝜗 12 2
𝑈0 = ((𝜎𝑥𝑥 + 𝜎𝑦𝑦 + 𝜎𝑧𝑧 ) + (𝜎𝑥𝑥 𝜎𝑦𝑦 + 𝜎𝑦𝑦 𝜎𝑧𝑧 + 𝜎𝑧𝑧 𝜎𝑥𝑥 )) + (𝜎𝑥𝑦 + 𝜎𝑧𝑥 +
2𝐸 𝐸 2𝐺
2
𝜎𝑥𝑦 ) (5.7)
U0 puede ser expresado en función de las cargas que actúan sobre la estructura; en estas condiciones la
ecuación 5.6 expresa U como una función de las cargas. Las ecuaciones 5.4 y 5.5 pueden entonces
utilizarse para encontrar los desplazamientos y rotaciones en los puntos de aplicación de cargas
concentradas y momentos puntuales en la dirección de dichas acciones. En el presente artículo se
consideran tres tipos de esfuerzos: flexión, cortante y torsión.

5.3.1 Energía en elementos sometidos a carga axial. (UP)

En este caso se presenta esfuerzo únicamente en la dirección del eje del elemento (z) con lo que la
energía de deformación se obtiene de integrar sobre el volumen del elemento la densidad de energía.

𝑧𝑧𝜎2
𝑈𝑃 = ∫ 𝑈𝑜 𝑑𝑣 = ∫ ∫ ∫ 2𝐸 𝑑𝑥𝑑𝑦𝑑𝑧

Para la sección transversal (ejes x,y) el esfuerzo (P/A=σzz) se considera constante, con lo que el resultado
de las dos primeras integrales corresponde con el área de la sección transversal, reemplazando se llega
a la expresión que representa la energía de deformación por efecto de una carga axial.
𝑃2
𝑈𝑃 = ∫ 2𝐴2 𝐸 [∫ ∫ 𝑑𝑥𝑑𝑦]𝑑𝑧

𝑃2
𝑈𝑃 = ∫ 2𝐴𝐸 𝑑𝑧 (5.8)

Para un resorte la energía de deformación se calcula mediante la expresión

𝑘
𝑈 = 2 𝛿2 (5.12)

5.3.2 Energías de deformación por flexión y por cortante en vigas. (UV)

Considere una viga de sección uniforme o con una variación moderada de su sección
transversal, definiendo el sistema de ejes coordenados con Z en la dirección del eje principal
de la viga, y en sentido vertical positivo hacia abajo y x saliendo del plano del dibujo, de tal
forma que la sección transversal se ubica en el plano xy.

El esfuerzo generado por flexión es:

𝑀𝑥 𝑦
𝜎𝑧𝑧 = 𝐼𝑥

La densidad de energía de deformación es entonces

1 2
𝑈0 = 2𝐸 ∗ 𝜎𝑥𝑥
La energía de deformación se obtiene entonces de
2 𝑦2
𝑀𝑥𝑥
𝑈𝑀 = ∫ 𝑈𝑜 𝑑𝑣 = ∫𝑧 ∫𝑦 ∫𝑥 2𝐸𝐼𝑥2
𝑑𝑥𝑑𝑦𝑑𝑧
Teniendo en cuenta que dA = dx*dy, Mxx es constante en la sección transversal (x,y) y que
la integral de 𝑦 2 𝑑𝐴 𝑒𝑠 𝐼𝑥 se llega a la expresión que representa la energía de deformación
por flexión.

𝑀2
𝑥𝑥
𝑈𝑀 = ∫ 2𝐸𝐼 𝑑𝑧 (5.13)

En relación con el cortante, el esfuerzo de cortante se evalúa mediante la expresión

𝑉∗𝑄
𝜎𝑥𝑧 = 𝐼∗𝑏

Este esfuerzo se asume constante en el eje x, es variable en Y y en Z.

Realizando cálculos para una sección rectangular con las dimensiones indicadas en la figura:

𝑏ℎ2 𝑏𝑦2
Q =b*(h/2-y)*[y+(h/2-y)/2] = 8
− 2
h/2

𝑏ℎ3
I=
y

EN 12

Reemplazando se llega a:

𝑏ℎ2 𝑏𝑦 2
12𝑉( − )
b
𝜎𝑥𝑧 = 8 2
𝑏 2 ℎ3

2 1 2 1
Recordando que UV=∫ 𝑈0 𝑑𝑣 = ∫ 2𝐺 𝜎𝑧𝑦 𝑑𝑣 = ∫ ∫ ∫ 2𝐺 𝜎𝑧𝑦 𝑑𝑦𝑑𝑥𝑑𝑧

Integrando inicialmente con respecto a y (2 veces entre 0 y h/2 para abarcar toda la sección)

ℎ/2 1 3𝑉 2
Se obtiene 2 ∗ ∫0 𝜎 *dy= 5𝑏2 𝐺ℎ
2𝐺 𝑧𝑦

Integrando ahora con respecto a X, entre 0 b queda entonces la siguiente expresión:

𝑏 ℎ/2 1 3𝑉 2
∫0 [2 ∗ ∫0 (2𝐺 𝜎𝑥𝑧 )dy]dx= 5𝑏𝐺ℎ, afectando por 2 tanto el numerador como el denominador

se llega a:

6𝑉 2 6 𝑉2 𝑉2
UV=∫ 5∗(2∗𝐺)∗𝑏∗ℎ 𝑑𝑧 = 5 ∫ (2∗𝐺)∗𝐴 𝑑𝑧 = 𝑘 ∗ ∫ (2∗𝐺)∗𝐴 𝑑𝑧 con k=1.2

Determinar este valor de K resulta en un cálculo complejo para algunas secciones, algunos
autores recomiendan tomar a k como la relación entre 𝜎𝑥𝑧 máximo y 𝜎𝑥𝑧 promedio, esto
puede generar un error, pero en no es muy significativo específicamente si se tomara esta
decisión para el caso del rectángulo se tendría para k un valor de 1.5 en vez de 1.2.

Se han desarrollado valores para diferentes secciones los cuales se relacionan a

continuación.

Sección Valor de K
Rectángulo 1.2 (Valor exacto)
Círculo sólido 1.33 Calculado como la relación del
esfuerzo en el eje neutro con relación al
esfuerzo promedio.
Sección circular de pared delgada 2.0 Calculado como la relación del esfuerzo
en el eje neutro con relación al esfuerzo
promedio
Sección Canal, Cajón 1.0 El área para cualquiera de estas
secciones es el área del alma b*h donde h
es la altura de la sección y b el ancho del
alma.

5.3.3 Energía de deformación generada por Torsión.

La energía de deformación causada en un elemento circular por efectos de torsión UT puede

ser evaluada como se presenta a continuación.

Los esfuerzos generados en la sección transversal por efectos de torsión corresponden a

esfuerzos cortantes, tomando un sistema de coordenadas similar al utilizado en los casos
anteriores es decir con el eje Z en el sentido longitudinal del elemento, el eje Y vertical
positivo hacia abajo y el eje X perpendicular al plano de, se tiene que los esfuerzos por
torsión actuarían sobre la cara XY.

1
De acuerdo con la ecuación 5.7 Uo=2𝐺σxy2

se sabe que: σxy = T*r/J. con lo que el esfuerzo es proporcional a la distancia al radio.

𝜋∗𝑅 4
Utilizando coordenadas polares con dA=r*dθ*dr y J= donde R corresponde al radio
2
máximo y r a cualquier radio entre 0 y R.
Se integra entonces sobre el volumen del elemento.

UT=∫ 𝑈0 𝑑𝑣 = ∫𝑧 ∫𝐴 𝑈0 𝑑𝐴𝑑𝑧.
 𝑅 2𝜋
𝑟𝑇 2
𝑟( )
𝜋𝑅4
2 𝑇2
Realizando la integral sobre el área se tiene: ∫ 𝑑𝜃 𝑑𝑟=𝐺𝜋𝑅4
2𝐺

∫ 0 0
Con lo que la energía de deformación por torsión se puede expresar como:

𝑇2
UT= ∫ 2𝐺𝐽 𝑑𝑧 15.16

Resumiendo la energía de deformación se puede evaluar mediante la suma de los aportes

de los diferentes esfuerzos causados por las acciones sobre la estructura.

𝑈 = 𝑈𝑃 + 𝑈𝑀 + 𝑈𝑉 + 𝑈𝑇

A partir de estas expresiones se puede entonces utilizar el teorema de Castigliano para

determinar la deflexión o la rotación de la estructura en el punto de aplicación de una carga
o un momento puntual (Fi , Mi)

Reemplazando en las ecuaciones 5.4 y 5.5.

𝑑𝑈 𝑃𝐽 𝜕𝑃𝐽 𝐾 𝑉 𝜕𝑉 𝑀𝐽 𝜕𝑀𝐽 𝑇𝐽 𝜕𝑇𝐽

∆𝑖 = 𝑑𝐹 = ∑𝑚
𝑗=1(∫ 𝐴 𝑑𝑧 + ∫ 𝐺𝐽𝐴𝐽 𝜕𝐹𝐽 𝑑𝑧 + ∫ 𝐸 𝑑𝑧 + ∫ 𝐺 𝑑𝑧 (5.17)
𝑖 𝐽 𝐸𝐽 𝜕𝐹𝑖 𝐽 𝐽 𝑖 𝐽 𝐼𝐽 𝜕𝐹𝑖 𝐽 𝐽𝐽 𝜕𝐹𝑖

𝑑𝑈 𝑃𝐽 𝜕𝑃𝐽 𝐾 𝑉 𝜕𝑉 𝑀𝐽 𝜕𝑀𝐽 𝑇𝐽 𝜕𝑇𝐽

𝜃𝑖 = 𝑑𝑀 = ∑𝑚
𝑗=1(∫ 𝐴 𝑑𝑧 + ∫ 𝐺𝐽𝐴𝐽 𝜕𝑀𝐽 𝑑𝑧 + ∫ 𝐸 𝑑𝑧 + ∫ 𝐺 𝑑𝑧 (5.18)
𝑖 𝐽 𝐸𝐽 𝜕𝑀𝑖 𝐽 𝐽 𝑖 𝐽 𝐼𝐽 𝜕𝑀𝑖 𝐽 𝐽𝐽 𝜕𝑀𝑖

Donde m corresponde al numero de elementos que componen la estructura

El procedimiento a seguir para evaluar las deflexiones utilizando el teorema de Castigliano

se presenta a continuación:

1. Escriba las expresiones que representen de las fuerzas internas (Axial, momento,
cortante, torsión) en cada miembro de la estructura en función de las fuerzas
externas aplicadas.
2. Derive cada una de las expresiones anteriores con respecto a una acción puntual (Fi
o Mi ) aplicada en el punto y en la dirección en la cual se desea conocer el
desplazamiento.
3. Substituya cada una de las expresiones obtenidas en los pasos anteriores en la
ecuación 5.17 o 5.18 según sea el caso y realice las integrales correspondientes.
4. Obtenga la magnitud buscada mediante la substitución de los valores numéricos
teniendo cuidado de emplear unidades consistentes (las rotaciones se obtienen en
radianes, los desplazamientos en las unidades utilizadas).
 Ejemplos.

E5.1 determine para la viga mostrada la deflexión del extremo libre generada por la
aplicación de la carga P

P
P

M
Z
Z
L V

Tomando el origen de coordenadas en el extremo libre y haciendo un corte a una distancia Z, se

plantean entonces las ecuaciones de equilibrio y se encuentran las expresiones que representan las
fuerzas internas.

V=P
M=P*z

Se derivan estas expresiones con respecto a la carga puntual aplicada en el sitio donde se desea
conocer la deflexión.

𝜕𝑉
=1
𝜕𝑃

𝜕𝑀
=𝑧
𝜕𝑃

Nótese que las expresiones de las derivadas corresponden con las fuerzas internas para cuando la
carga aplicada (p) es una carga unitaria.

Substituyendo ahora estas expresiones en la Ecuación 5.17

𝑑𝑈 𝑃 𝑃∗𝑧
∆𝑖 = 𝑑𝐹 = ∫(𝐸𝐴 ∗ 1)𝑑𝑧 + ∫( 𝐸𝐼 ∗ 𝑧)𝑑𝑧
𝑖

𝑃𝐿 𝑃𝐿3
∆= 𝐴𝐸 + 3𝐸𝐼

El primer termino corresponde al aporte del cortante y el segundo al aporte de la flexión. En general
el aporte por cortante es poco representativo y en la mayoría de los casos puede no tenerse en
cuenta sin que esto genere un error representativo en la respuesta, específicamente para vigas con
relación L/h>5 el efecto del cortante puede despreciarse.

Para fines prácticos se considerará una viga de 25*35 (b*h) con una longitud de 2 m, para este caso
200 0.229𝑃
la relación L/h=200/35=5.71, ∆𝑉 = 𝑃 ∗ = , por efectos de flexión se tiene ∆𝑀 =
25∗35𝐸 𝐸
2003 89.6𝑃
𝑃 1 = , con lo que se puede evidenciar que el aporte del cortante es muy pequeño
𝐸∗25∗353 E
12
comparado con lo que genera la flexión.

E5.2 determine para la viga mostrada la deflexión del extremo libre generada por la
aplicación de la carga uniformemente distribuida W a lo largo de toda la viga.

W
Se hará uso de lo encontrado en le ejercicio
anterior, es decir a cambio de utilizar la derivada, se
Z utilizara el termino correspondiente a la fuerza
L interna producida por una carga unitaria aplicada en
el punto donde se quiere determinar el
desplazamiento.

W
A partir del corte se encontrarán las expresiones para
las fuerzas internas generadas por el sistema de
M cargas aplicado (M) y por una unitaria aplicada en el
Z
V extremo libre.

𝑉 = 𝑊∗𝑧 v=1

𝑧2
𝑀= 𝑊∗ 𝑚 = 1∗𝑧 =𝑧
2

Las expresiones a integrar son entonces:

𝑉∗𝑣 𝑀∗𝑚
∫ 𝑑𝑧 𝑦 ∫ 𝑑𝑧
𝐴𝐺 𝐸𝐼

𝐿
1
𝐿 ∫ 𝑊𝑧 2 𝑧 ⅆ𝑧
∫0 𝑊𝑧 ⅆ𝑧 0
2 𝑊𝐿2 𝑊𝐿4
+ = +
𝐴𝐺 EI 2𝐴𝐺 8EI

Igual que en el ejercicio anterior el primer termino corresponde al aporte del cortante y el segundo
al aporte de la flexión.
En caso que el objetivo sea conocer la rotación del extremo libre los valores correspondientes a la
carga unitaria se reemplazarán por los generados por un momento unitario aplicado igualmente en
el punto donde se desea conocer la rotación.

m=mi m mi=1
v=0. Z
V

Con lo cual se evidencia que el cortante no genera rotación en el extremo libre, este se produce
únicamente como resultado de la flexión.
𝐿 𝑀∗𝑚 𝐿 𝑊𝑧 2 𝑊∗𝐿3
𝜃 = ∫0 𝑑𝑧 = ∫0 ( ∗ 1) 𝑑𝑧 = en el mismo sentido de rotación del momento aplicado.
𝐸𝐼 2∗𝐸𝐼 6𝐸𝐼

E5.3 determine la deflexión que presenta el extremo libre de la viga curva presentada
en la figura por efecto de la carga aplicada.

considere: la sección transversal es cuadrada de 30 mm de lado,

radio de curvatura R=65 mm, la viga está construida de acero E=200
GPa, ν=0.29.

haciendo un corte sobre el eje girado un ángulo θ respecto de la

línea horizontal, se encuentran las expresiones para N, V y M que
se presentan a continuación.

N=P*cos(θ)
V=P*seno (θ)
M=P*R*(1-cos(θ))
Ahora, aplicando una carga unitaria en el extremo y que corresponde con el punto de aplicación de
la carga P, los valores de las fuerzas internas generadas por esta carga unitaria serán.

N=cos(θ)

V=Seno(θ)

M=R*(1-cos(θ))

Para efectos de hacer la integral a lo largo del elemento se utilizaran coordenadas polares (dz=rd θ).

En estas condiciones.

1 3𝜋/2
∆= ∫ (𝑃 ∗ cos(θ) ∗ cos(θ) ∗ R) 𝑑θ
𝐴𝐸 0
1 3𝜋/2 3𝜋/2
+ ∫ 𝑘 ∗ 𝑃 ∗ 𝑠𝑒𝑛(𝜃) ∗ sin(𝜃) ∗ 𝑅𝑑𝜃 + 1/(EI) ∫ 𝑃𝑅(1 − cos(𝜃)) ∗ 𝑅(1 − cos(𝜃) ∗ 𝑅𝑑𝜃
𝐴𝐺 0 0

Axial cortante Flexion

9𝜋 3
3∗𝑃∗𝜋∗𝑅 3𝑘𝑃𝜋𝑅 𝑃(2+ )𝑅
4
= + +
4AE 4AG EI

Reemplazando numéricamente se llega a:

.005105+.01580+1.107=1.127 mm

Se evidencia nuevamente que en este caso el mayor aporte es generado por la flexión sobre el
elemento como se muestra en la siguiente tabla.

Axial Cortante Flexión

0.005105 0.015675 1.107

0.453% 1.390% 98.157%

Con lo que puede concluirse que en casos como este es justificable utilizar únicamente el aporte
generado por la flexión para el calculo de la deflexión, condición que agiliza y simplifica los cálculos.
 E5.4 la viga semicircular en cantiléver mostrada en la figura soporta 2 cargas de
magnitud P aplicadas en los puntos B y C, la viga tiene una sección transversal circular de
diámetro d. determine la deflexión vertical del punto C.

En este caso deben considerarse 2 intervalos, uno

dependiente de Φ (0-π/2) y el otro dependiente de θ
(0-π/2). Igualmente se utilizará el sistema de
coordenadas polares para facilitar la solución.

Para el primer intervalo: CB

Por Cargas aplicadas Por carga Unitaria en C

N=P*cos(Φ) n=cos(Φ)

V=P*seno(Φ) v=Seno(Φ)

M=P*R(1-cos(Φ)) p=R(1-cos(Φ))

Para el segundo intervalo.

Por Cargas aplicadas Por carga Unitaria en C

N=(P+Q)*Seno(θ) n=Seno(θ)

V=(P+Q)*Cos(θ) v=Cos(θ)

M= P*R(1+seno (θ)+Q*R*Seno(θ) p=R(1+seno(θ))

Integrando cada uno de los términos en los intervalos correspondientes (teniendo en cuenta que
dz=RdΦ o Rdθ según el intervalo) se llega a:

Para el intervalo BC
3𝜋 3
𝑝𝜋𝑅 𝑘𝑝𝜋𝑅 𝑝(−2 + 4 )𝑅
∆𝐵𝐶 = + +
4𝐴𝐸 4𝐴𝐺 EI
Para el intervalo AB.

𝑝𝜋𝑅 𝑘𝑝𝜋𝑅 𝑝(3 + 𝜋)𝑅 3

∆𝐴𝐵 = + +
2𝐴𝐸 2𝐴𝐺 EI
Con lo que el desplazamiento total es:
 3𝑝𝜋𝑅 3𝑘𝑝𝜋𝑅 𝑝(4+7𝜋)𝑅3
∆= 4𝐴E
+ 4𝐴𝐺
+ 4EI

E5.5 encuentre la deflexión en el extremo de libre de la viga mostrada.

se trabaja igualmente en dos intervalos. Se utilizará

únicamente el aporte de flexión.

Para determinar la deflexión del punto A se tendrá en

cuenta la carga unitaria FA

MAB= 0 mAB = z.

MBC= P(z-L/2) mBC = z.

En estas condiciones se plantea la ecuación para deflexión del punto A

1 𝐿 5𝐿3 𝑃
Δ=𝐸𝐼 ∗ ∫𝐿⁄2 𝑃 ∗ (𝑧 − 𝐿⁄2) ∗ 𝑧 𝑑𝑧 = 48𝐸𝐼

Para determinar la rotación del punto A se tiene en cuenta un momento unitario MA.

MAB= 0 mAB = 1

MBC= P(z-L/2) mBC = MA =1

𝐿
1 𝐿2 𝑃
𝜃= ∗ ∫ 𝑃 ∗ (𝑧 − 𝐿⁄2) 𝑑𝑧 =
𝐸𝐼 𝐿⁄2 8𝐸𝐼
 E5.6 La viga semicircular de radio R mostrada en la figura tiene una sección
transversal circular de radio r, se encuentra fija en O y se ubica en la plano xy con el centro
de curvatura en C , una carga P paralela al eje z actúa en el punto A, determine la componente
z de la deflexión del punto B.

Se supone entonces una carga unitaria en dirección z aplicada en B

Considerando la sección D en la viga curva que se ubica a
un Angulo θ, medido desde el eje AC, se plantean
entonces las expresiones para determinar el momento
flexionante y el torque resultantes de las fuerza aplicadas
en A y B.

M=P*R*seno(θ) m=R*cos(θ)

T=P*R(1-cos(θ)) t=R*(1+seno (θ))

Como quiera que en el tramo AB no se presenta ni M ni

T no se tendrá en cuenta para la evaluación.

𝜋⁄2
Δ=1⁄(EI) ∫0 𝑃 ∗ 𝑅 ∗ Sin[𝜃] ∗ 𝑅 ∗ Cos[𝜃] ∗ 𝑅 𝑑𝜃 +
𝜋⁄2
1/(𝐽𝐺) ∫ 𝑃 ∗ 𝑅(1 − Cos[𝜃]) ∗ 𝑅 ∗ (1 + Sin[𝜃]) ∗ 𝑅 𝑑𝜃
0

𝑃𝑅3 𝑃𝑅3 (𝜋−1)

De donde se obtiene: Δ= +
2EI 2𝐺𝐽

𝐸 𝜋∗𝑟 4 𝜋∗𝑟 4
Si se reemplaza G por 2(1+𝜗) , I= 4
, 𝐽= 4
se obtiene:

2𝑃𝑅3
Δ= E𝜋𝑟4 (1 + (1 + 𝜗)(𝜋 − 1))

5.3.4 Deflexión en estructuras estáticamente indeterminados.

Para el caso de estructuras estáticamente indeterminadas se hace uso del principio de
superposición; una estructura no se altera si se retiran elementos redundantes o soportes y se
reemplazan estos por fuerzas y momentos externos que sean idénticos a los aportados por los éstos.
Estas fuerzas y momentos son inicialmente desconocidos, pero es posible plantear ecuaciones de
condición que garanticen que la estructura simplificada (estáticamente determinada) tenga el
comportamiento de la estructura original, y a partir de éstas determinar la magnitud de las fuerzas
y momentos aplicados.
 E5.7 Solucione la viga mostrada.
a) Determinar la magnitud de la reacción en B

A B A B B
W W W
M V
MA
RA RB RB
L L X

Ecuación de condición ΔB = 0

Atendiendo los lineamientos presentados, se elimina el apoyo en B, se reemplaza por una fuerza
puntual, y se genera una ecuación de condición que garantice el comportamiento de la estructura
de acuerdo con su condición inicial es decir que el desplazamiento del punto B sea nulo.

Se plantean las ecuaciones necesarias para determinar la deflexión en B, con lo cual se supondrá la
carga unitaria actuando en B.

Considerando únicamente el aporte de flexión se tiene:

M= Wx2/2-Rx m=x

Con lo que se obtiene la expresión que determina la deflexión.

1 𝐿 1 𝐿 𝑤𝑥 2 𝐿3 𝑅 𝐿4 𝑊
∆= 𝐸𝐼 ∫0 𝑀 ∗ 𝑚 𝑑𝑥 = 𝐸𝐼 ∫0 ( 2
− 𝑅𝑥) ∗ 𝑥 𝑑𝑥 = − 3
+ 8

Aplicando la ecuación de condición se obtiene la magnitud de R

𝐿3 𝑅 𝐿4 𝑊 3𝑊𝐿
− + =0 ⇒𝑅=
3 8 8
b) Determinar la deflexión en L/2.

Si se desea conocer la deflexión en un punto intermedio es posible aplicar nuevamente el

procedimiento aplicando la carga unitaria en el punto que se desea determinar.
1.0
A B A B B
W W W
M V
MA
RA
C 3/8*WL 3/8*WL
L L X

Se requieren ahora dos intervalos AC y CB.

Intervalo AC Intervalo CB
M m M m
W*x2/2-3/8*WL 1.0*(X-L/2)=X-L/2 W*x2/2-3/8*WL 0

Con lo que el tramo BC se anula y la integral se debe realizar únicamente en el tramo AC.
 𝐿
1 𝐿 𝑊𝑋 2 3𝐿𝑊𝑋 𝐿4 𝑊
∆𝐿/2 = ∫ (𝑋 − ) ( − ) 𝑑𝑋 =
𝐸𝐼 𝐿 2 2 8 192EI
2

E5.8 determine las reacciones en C para el elemento mostrado ABC y la deflexión del punto B en la
dirección de P. no considere los efectos de cortante y axial.

El apoyo C permite rotación pero restringe los

desplazamientos vertical y horizontal. Se procede
entonces a considerar como redundantes las
componentes de reacción en C. (Cz, Cy) y sen tendrán dos
ecuaciones de condición relacionadas con los
desplazamientos horizontal y vertical del punto C.

∆𝑐𝑧 = 0, ∆𝑐𝑦 = 0

Se considera entonces dos condiciones primero con una carga unitaria actuando verticalmente en
C, y a continuación con la carga unitaria actuando en C en dirección horizontal.

Carga unitaria actuando en dirección vertical aplicada en C

Tramo CB TRAMO BA
M m M m
-Cy*R*seno(θ)+Cz*R*(1-cos(θ)) -R* seno(θ) -(P-Cy)*s + Cz*2R s

La deflexión vertical del punto se obtiene mediante la siguiente ecuación.

1 𝜋 1 2𝑅
∫ (−cy ∗ 𝑅 ∗ Sin[θ] + cz ∗ 𝑅 ∗ (1 − Cos[θ])) ∗ (−𝑅 ∗ Sin[θ]) ∗ 𝑅 𝑑θ + ∫ (−(𝑃 − cy) ∗ 𝑠 + cz ∗ 2 ∗ 𝑅) ∗ 𝑠 𝑑𝑠
𝐸𝐼 0 𝐸𝐼 0

8cy𝑅3 8𝑝𝑅3 1
= + 4cz𝑅3 − + (−4cz + cy𝜋)𝑅3 = 0 (a)
3 3 2

Aplicando ahora la carga unitaria en C en dirección horizontal.

Tramo CB TRAMO BA
M m M m
-Cy*R*seno(θ)+Cz*R*(1-cos(θ)) R* (1-cos(θ)) -(P-Cy)*s+Cz*2R 2R

Se realiza la integral presentada a continuación.

1 𝜋
∫ (−cy ∗ 𝑅 ∗ Sin[θ] + cz ∗ 𝑅 ∗ (1 − Cos[θ])) ∗ 𝑅 ∗ (1 − Cos[θ]) ∗ 𝑅 𝑑θ
𝐸𝐼 0
1 2𝑅
+ ∫ (−(𝑃 − cy) ∗ 𝑠 + cz ∗ 2 ∗ 𝑅) ∗ 2𝑅 𝑑𝑠
𝐸𝐼 0

3cz𝜋
= 4cy𝑅3 + 8cz𝑅3 − 4𝑝𝑅3 + (−2cy + ) 𝑅3 = 0 (b)
2
Con lo que resulta un sistema de dos ecuaciones con dos incógnitas. (a) y (b).

Solucionando se llega a:

cz = 0.2329𝑝, cy = 0.5194𝑝
Para encontrar la deflexión en el punto B, el momento “m” se aplica la carga unitaria en B,
con lo que para el tramo BC esta magnitud es nula, entonces la deflexión se calcula únicamente
con el aporte del tramo BC y considerando la carga unitaria aplicada en B en dirección vertical
(hacia abajo).
1 2𝑅 0.35𝑃𝑅3
Deflexión en B=𝐸𝐼 ∫0 (−(𝑝 − 𝐶𝑦) ∗ 𝑠 + 𝐶𝑧 ∗ 2𝑅 ) ∗ −𝑠𝑑𝑠 = 𝐸𝐼
.

E5.8 solución armadura mixta.

El Rey invertido esta construido con una viga rectangular de acero ABC de 60*160 mm (b*h) Es=200
GPa, σy=240 MPa, una varilla de 15 mm de diámetro ADC, E =200 GPa, σy=500 MPa, y un puntal de
roble blanco BD de 40*40 mm EBD=12.4 GPa, YBD=29.6 MPa. Determine la máxima magnitud de la
carga P que puede ser aplicada a la estructura si todos sus elementos se diseñan usando un factor
de seguridad de 2.0 respecto de su esfuerzo de fluencia. No considere las concentraciones de
esfuerzos.

Se determinan inicialmente los esfuerzos admisibles y las propiedades geométricas de los diferentes
elementos.

Elemento Esfuerzo admisible (MPa)

VIGA ABC 120
VARILLA ADC 250
PUNTAL BD 14.8

En esta estructura, la viga trabaja a compresión, cortante y a flexión, el puntal trabaja en compresión
y las varillas a tracción. Para efectos de evaluar la energía de deformación se tendrán en cuenta los
aportes de: flexión en la viga y carga axial en el puntal y las varillas, los demás esfuerzos no generan
aportes significativos y pueden ser no considerados sin que esto genere un error apreciable en los
resultados.
Aabc = 9600 mm2, IABC=20.48*106 mm4

ABD= 1600 mm2 ADC= AAD 176.7 mm2

En este caso se aplica como redundante el puntal y de acuerdo con el teorema de Castigliano la
derivada de la energía interna con respecto a esta fuerza interna es igual a 0.

Analizando el nudo D:

Σ𝐹𝑦 = 0 − 𝑁𝐵𝐷 + 2 ∗ 𝑁𝐷𝐶 ∗ 𝑠𝑒𝑛𝑜(0.245) = 0

NDC = 2.061NDB (a)

Considerando como carga unitaria los efectos de NBD, en los Nudos B y D (Llegando al nudo en cada
caso). El elemento ABC está a flexión, el elemento BD a compresión y el elemento DC a tracción, los
termino que representa estas condiciones para cada elemento son:

Elemento Flexión Axial

M m N n
AB, CB (P/2-NBD/2)*s -s/2
DB NBD 1
DC DA 2.061 NBD 2.061

𝑑𝑈
Teniendo como ecuación de condición que 𝑑𝑁 =0
𝐵𝐷

1 2000 P NBD 1 500

Se llega a: 2 ∗ ∫ ( − )∗ s ∗ ( − s/2)ds + ∫ 𝑁𝐵𝐷 𝑑𝑦 +
𝐸𝐵𝐶 𝐼𝐵𝐶 0 2 2 𝐴𝐵𝐷 𝐸𝐵𝐷 0

1 2016
2∗𝐴 ∗ ∫0 2.061𝑁𝐵𝐷 ∗ 2.061 𝑑𝑥 =0 (Se trabajó en N y mm)
𝐷𝐶 𝐸𝐷𝐶

Solucionando esta ecuación se llega a: NBD=0.390P.

Y reemplazando en las otras ecuaciones se tiene:

NDC = 2.061NDB = 0.804P

Máximo momento en ABC (para s=2000) = 610P.

Verificando para cada elemento los esfuerzos admisibles:

BD 0.390P/1600=14.8 P=60700 N.

BC 0.804P/176.7=250 P= 55000 N.

ABC 610P*80/2048*106=120 P= 50400 N.

Con lo que la carga de diseño es 50.400 N.