Simetrias
Simetrias
7 ESTRUCTURAS SIMÉTRICAS
1
q q
∆T ∆T ∆T - ∆T
P P P P
a.-Estructura simétrica con carga simétrica b.-Estructura simétrica con carga antisimétrica
∆T2 ∆T1 ∆T1 /2 ∆T2 /2 ∆T2 /2 ∆T1 /2 -∆T1 /2 -∆T2 /2 ∆T2 /2 ∆T1 /2
= +
c.- Estructura simétrica con carga general Estado simétrico Estado antisimétrico
FIGURA 4.7.1
2
= = = =
FIGURA 4.7.2
Los esfuerzos y desplazamientos en un punto de una estructura simétrica
con carga simétrica situado en el eje de simetría, podrán existir si son simétricos
y serán nulos si son antisimétricos. En el caso de que la carga sobre una estructura
simétrica sea antisimétrica, en un punto del eje de simetría podrán existir los
movimientos y esfuerzos que sean antisimétricos, mientras que los simétricos
serán nulos. En la figura 4.7.3 se hace un resumen de los de los esfuerzos y
movimientos en una rebanada diferencial perteneciente a un eje de simetría
indicándose cual debe ser su tratamiento al hacer la simplificación por simetría.
En la figura 4.7.4 se indica la simplificación para los ejemplos mostrados en la
figura 4.7.1.
CARGA SIMÉTRICA CARGA
ANTISIMÉTRICA
Momento flector (M) M≠0 M=0
x, u
θ SIMPLIFICACIÓN
y, v
FIGURA 4.7.3
En el caso de que una barra de una estructura simétrica coincida con el eje de
simetría, las anteriores consideraciones son válidas para todos los puntos del eje de
simetría. Por tanto, si la carga es simétrica, el desplazamiento perpendicular al eje de
simetría y el giro son nulos, lo que supone que el momento flector y el esfuerzo cortante
3
que ellos originarían, serán nulos (figura 4.7.5.a). Al realizar la simplificación, el esfuerzo
axil deberá ser la mitad del que tendría la estructura completa, por lo que para que la
elástica sea real, se debe asignar a dicha barra en la simplificación un valor mitad del área
de la sección real (con la mitad de la fuerza axial se obtendrá la deformación real) (figura
4.7.5.b).
q q
∆T ∆T ∆T - ∆T
P P P P
q q
∆T ∆T
P P
a.-Estructura simétrica con carga simétrica b.-Estructura simétrica con carga antisimétrica
FIGURA 4.7.4
Al realizar la simplificación en el caso de que la carga sobre una estructura
simétrica sea antisimétrica, sobre una barra coincidente con el eje de simetría el esfuerzo
axil será nulo y el valor del momento flector el esfuerzo cortante será la mitad del de la
estructura completa. Por tanto, para que la elástica de la estructura simplificada sea igual
a la real, se debe asignar a dicha barra un momento de inercia cuyo valor sea la mitad del
de la sección real (ver figura 4.7.5.c).
= =
FIGURA 4.7.5
4
Las ecuaciones generales de comportamiento de la viga de la figura 4.7.6 son:
0 𝛿𝑖𝑗
𝑀𝑖𝑗 = 𝑀𝑖𝑗 + 𝑘𝑖𝑗 (𝜗𝑖𝑗 + 𝛽𝑖𝑗 𝜗𝑗𝑖 ) − 𝑘𝑖𝑗 (1 + 𝛽𝑖𝑗 ) 𝐿
𝛿𝑖𝑗
(4.7.1)
𝑀𝑗𝑖 = 𝑀𝑗𝑖0 + 𝑘𝑗𝑖 (𝜗𝑗𝑖 + 𝛽𝑗𝑖 𝜗𝑖𝑗 ) − 𝑘𝑖𝑗 (1 + 𝛽𝑗𝑖 ) 𝐿
FIGURA 4.7.6
5
EJEMPLO 4.6
Calcular los desplazamientos en todos los nudos de la estructura representada en
la figura 4E6.1, así como los esfuerzos en todas sus barras.
Datos: EA = 3x104 kN
D E
4m
16 T
A C
F B G 3m
3m 3m 3m
FIGURA 4E6.1
c d
8T
y
a b
x A
F B
FIGURA 4E6.2
La matriz de rigidez de la estructura, tras aplicar las condiciones de contorno en
los nudos A y D, es:
𝒌𝒂 + 𝒌𝒃𝟏𝟏 + 𝒌𝒄𝟏𝟏 𝒌𝒃𝟏𝟐
𝑲 = [ 𝟐𝟐 ] (4E6.1)
𝒌𝒃𝟐𝟏 𝒌𝒃𝟐𝟐 + 𝒌𝒅𝟏𝟏
Siendo las matrices elementales en coordenadas globales:
𝐸𝐴
0 10000 0
𝒌𝒂22 = [ 𝐿𝑎 ]=[ ] = 𝒌𝒃𝟏𝟏 = 𝒌𝒃𝟐𝟐 = −𝒌𝒃𝟏𝟐 = 𝒌𝒃𝟐𝟏
0 0 0 0
0 −1 7500 0 0 1 0 0 (4E6.2)
𝒌𝒄𝟏𝟏 = [ ][ ][ ]=[ ]
1 0 0 0 −1 0 0 7500
𝒄 −0.6 −0.8 7500 0 −0.6 0.8 2160 −2880
𝒌11 =[ ][ ][ ]=[ ]
0.8 −0.6 0 0 −0.8 −0.6 −2880 3840
6
Sustituyendo las matrices elementales en la matriz de la estructura (4E6.1):
0 20000 0 −10000 0 𝑢𝐹
{ 0 }=[ 0 7500 0 0 ] {𝑣𝐹 }
𝑋𝐵 −10000 0 12160 −2880 0
−8 0 0 −2880 3840 𝑣𝐵
0 20000 0 0 𝑢 𝐹
{ 0 }=[ 0 7500 0 ] {𝑣𝐹 }
−8 0 0 3840 𝑣𝐵
Que resuelto :
VB = -2.0833x10-3 m; uF = vF = uB = 0
Los esfuerzos son:
Na = N b = N c = 0
6000 0 −0.6 −0.8 0 10
𝑁𝑑 = [ ][ ]{ −3 } = { }kN
0 0 0.8 −0.6 −2.0833𝑥10 0
EJEMPLO 4.7
Calcular el movimiento del nudo A de la estructura de la figura 4E7.1, así como
las reacciones en C y E.
Datos: Para todas las barras: EI = 100 kNm2, EA = 105 kN, q = 10 kN/m
1m q
q A
B
D
q
1m q
1m 1m
FIGURA 4E7.1
Se trata de una estructura de nudos con dos ejes de simetría tal y como se indica
en la figura 4E7.1, por lo que en el nudo A los desplazamientos son nulos y solo quedará
el giro, resultando la estructura simplificada de la figura 4E7.2, en la que se debe
7
considerar la mitad de la rigidez a flexión (valor I/2) en las barras AB y AC, y el término
de la flexión axial no intervendrá.
C
b
q/2
q/2
B
a A
FIGURA 4E7.2
La matriz de rigidez de la estructura, tras aplicar las condiciones de contorno en
los nudos A y D, es:
𝐸𝐴
0 0
𝐿 600 0 −300
𝑐𝑜𝑠𝛼 −𝑠𝑒𝑛𝛼 0 12𝐸𝐼 6𝐸𝐼 𝑐𝑜𝑠𝛼 𝑠𝑒𝑛𝛼 0 𝐸𝐴
𝑏
𝒌11 = [𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 0] 0 [−𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 0] = [ 0 0 ]
0 0 1 𝐿3 𝐿2 0 0 1 𝐿
6𝐸𝐼 4𝐸𝐼 −300 0 200
[0 𝐿2 𝐿 ]
El vector de cargas:
𝑞𝐿2 5
𝐵𝑎𝑟𝑟𝑎 𝑎: 𝑀𝐴 = − =−
8 8 𝑀 = − 5 − 5 = − 25
𝐴
𝑞𝐿2 5 8 12 24
𝐵𝑎𝑟𝑟𝑎 𝑏: 𝑀𝐴 = − =− }
12 12
8
Sustituyendo las matrices elementales en la matriz de la estructura (4E7.1):
𝐸𝐴
𝑋𝐴 + 600 0 −300
𝐿 0
𝑌𝐴
{ }= 𝐸𝐴 {0}
25 0 + 150 −150 𝜗
− 𝐿 𝐴
24 [ −300 −150 350 ]
Que resuelto :
θA = -0.002976 rad
Las reacciones en C (Se recuerda que, como se verá a continuación, los valores
obtenidos de la estructura simplificada deben multiplicarse por dos):
Debidos al giro:
0
𝑹 𝑘𝑏 𝑏
𝑘12 𝜹𝑨
{ 𝐴 } = [ 11 𝑏 𝑏 0
𝑏 ] 𝑳 { 𝟎 } ⇒ 𝑹𝐶 = 𝑘21 𝐿𝜹𝑨 = 𝑘21 𝐿 { }
𝑹𝐶 𝑏
𝑘21 𝑘22 𝜗 𝐴
EJEMPLO 4.8
Calcular los movimientos de los nudos A y B de la estructura de la figura 4E8.1,
considerando para su cálculo la simetría de la misma.
Datos: Para todas las barras: EI = 4.5x104 kNm2, EA = 1.6x106 kN, q = 10kN/m
9
10kN/m
B C
2m
D 100 m kN E
A
FIGURA 4E8.1
y, v
= +
θ x, u 50 50 50
D E D A D A
A
Simétrica Antisimétrica
FIGURA 4E8.2
LBD = LBA = 2.5 m
Las ecuaciones matriciales, tras aplicar las condiciones de contorno del nudo E,
correspondientes a cada uno de estos estados dos estados son:
- Caso simétrico:
𝑝 𝐿2 10 𝑥 32
𝑀𝐵𝐶 = − =− = −7.5 𝑚𝑘𝑁
12 12
𝑝𝐿 10 𝑥 3
𝑅𝐵𝑦 =− =− = −15 𝑘𝑁
2 2
10
MBC = 50 mkN MB = MA/2 =25 mkN; RB x 2.5=75⇒ RB=30 kN, cuyas
componentes son: RBx= -30 x sen α= -30 x 0.8 = -24 kN; RBy= -30 x 0.6 = -18 kN
𝑢𝐵 0
𝑣𝐵 −15
𝒌𝑩𝑫 + 𝒌𝑨𝑩 𝑩𝑪′
𝟐𝟐 + 𝒌𝟏𝟏 𝒌𝑩𝑪′ 𝜗𝐵 −7.5
𝑲 = [ 𝟐𝟐 𝟏𝟐
𝑩𝑪′ ] = 𝑋𝐶′ (4E8.1)
𝑩𝑪′
𝒌𝟐𝟏 𝒌𝟐𝟐 0
𝑣𝐶′ 0
{0} { 𝑀𝐶′ }
Donde para la barra DB, se usará la matriz de rigidez elemental correspondientes
la barra articulada-empotrada y para las barras AB y BC’ las de la barra empotrada-
empotrada.
- Caso antisimétrico:
Las cargas debido al momento sobre la barra BA, son:
MAB = 50 mkN MBA = MA/2 =25 mkN; RB x 2.5=75⇒ RB=30 kN, cuyas
componentes son: RBx= -30 x sen α= -30 x 0.8 = -24 kN; RBy= -30 x 0.6 = -18 kN
𝑢𝐵 −24
𝑣𝐵 −18
𝒌𝑩𝑫 𝑨𝑩
𝟐𝟐 + 𝒌𝟐𝟐 + 𝒌𝟏𝟏
𝑩𝑪′
𝒌𝑩𝑪′ 𝜗𝐵
𝑲=[ 𝟏𝟐
𝑩𝑪′ ] 𝑢 = −25 (4E8.2)
𝑩𝑪′
𝒌𝟐𝟏 𝒌𝟐𝟐 𝐶′ 0
0 𝑌𝐶′
{𝜗𝐶′ } { 0 }
Donde para las barras DB y AB, se usarán las matrices de rigidez elementales
correspondientes la barra articulada-empotrada y para la barra BC’ la de la barra
empotrada-empotrada.
También se puede plantear:
0 𝑋𝐴
0 𝑌𝐴
𝜗𝐴 50
𝒌𝑨𝑩
𝟏𝟏 𝒌𝑨𝑩
𝟏𝟐 𝟎 𝑢𝐵 0
𝑲 = [𝒌𝑨𝑩
𝟐𝟏 𝒌 𝑩𝑫
𝟐𝟐 + 𝒌𝑨𝑩
𝟐𝟐 + 𝒌𝑩𝑪′
𝟏𝟏 𝒌𝑩𝑪′
𝟏𝟐 ] 𝑣𝐵 = 0 (4E8.3)
𝟎 𝒌𝑩𝑪′ 𝒌𝑩𝑪′ 𝜗𝐵 0
𝟐𝟏 𝟐𝟐
𝑢𝐶′ 0
0 𝑌𝐶′
{𝜗𝐶′ } {0}
Donde para la barra DB, se usará la matriz de rigidez elemental correspondientes
la barra articulada-empotrada y para las barras AB y BC’ las de la barra empotrada-
empotrada.
11
𝐸𝐴
0 0
𝐿
𝑐𝑜𝑠𝛼 −𝑠𝑒𝑛𝛼 0 3𝐸𝐼 3𝐸𝐼 𝑐𝑜𝑠𝛼 𝑠𝑒𝑛𝛼 0
𝒌𝐵𝐷
22
𝑇 𝐵𝐷
= 𝑳 𝒌′22 𝑳 = [𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 0] 0 − 2 [−𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 0]
0 0 1 𝐿3 𝐿 0 0 1
3𝐸𝐼 3𝐸𝐼
[0 − 2
𝐿 𝐿 ]
23.593 30.306 1.728
= 104 [30.306 41.271 −1.296]
1.728 −1.296 5.4
∗ 𝐶𝑎𝑠𝑜 𝑎𝑛𝑡𝑖𝑠𝑖𝑚é𝑡𝑟𝑖𝑐𝑜:
𝐸𝐴
0 0
𝐿
3𝐸𝐼 3𝐸𝐼 23.593 −30.306 1.728
𝒌𝐵𝐴
22 =𝑳 𝑇
0 4
− 2 𝑳 = 10 [−30.306 41.271 1.296]
𝐿3 𝐿 1.728 1.296 5.4
3𝐸𝐼 3𝐸𝐼
[0 − 2
𝐿 𝐿 ]
12
𝐸𝐴 𝐸𝐴
0 0 − 0 0
𝐿 𝐿
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
0 0 −
𝐿3 𝐿2 𝐿3 𝐿2
6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼
0 0 − 2
𝒌𝐵𝐶′
= 𝐿2 𝐿 𝐿 𝐿
𝐸𝐴 𝐸𝐴
− 0 0 0 0
𝐿 𝐿
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
0 3
−
− 2 0 3
− 2
𝐿 𝐿 𝐿 𝐿
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼
[ 0 𝐿2 𝐿
0 − 2
𝐿 𝐿 ]
106.667 0 0 −106.667 0 0
0 16 12 0 −16 12
= 0 12 12 0 −12 6
−106.667 0 0 106.667 0 0
0 −16 −12 0 16 −12
[ 0 12 6 0 −12 12 ]
- Caso antisimétrico:
Con el planteamieto dado en 4E8.2, se obtiene:
153.853 0 3.456 −106.667 0 0 𝑢𝐵 −24
0 98.542 12 0 −16 12 𝑣𝐵 −18
3.456 12 22.8 0 −12 6 𝜗 𝐵 −25
104 𝑢𝐶′ =
−106.667 0 0 106.667 0 0 0
0 −16 −12 0 16 −12 0 𝑌𝐶′
[ 0 12 6 0 −12 12 ] {𝜗𝐶′ } { 0 }
Que al ser resuelto da los siguientes resultados:
13
uB= -0.424x10-4m; vB=-0.129x10-4m; ϑB=-1.15 x10- 4rad; uC’=-0.424 x10-4m;
ϑC’=0.703 x10-4rad
El movimiento de B, será la suma del obtenido en los casos simétrico y
antisimétrico, lo que da los siguientes resultados para sus componentes:
uB= -0.407x10-4m; vB= -0.302x10-4m; ϑB= -1.622x10-4 rad
Para obtener el giro en A, solo se precisa el caso antisimétrico ya que este será nulo
en el simétrico y lo primero será obtener los movimientos de B en coordenadas locales.
𝑢′𝐵 𝑐𝑜𝑠𝛼 𝑠𝑒𝑛𝛼 0 𝑢𝐵
{𝑣′𝐵 } = [−𝑠𝑒𝑛𝛼 𝑐𝑜𝑠𝛼 0] {𝑣𝐵 } 𝑐𝑜𝑛 𝑠𝑒𝑛𝛼 = 0.8 𝑦 𝑐𝑜𝑠𝛼 = −0.6
𝜗′𝐵 0 0 1 𝜗𝐵
Obteniédose al sustituir, los siguientes movimientos en coordenadas locales:
u’B= 0.151x10-4m; v’B= 0.4166x10-4m; ϑ’B= ϑB = 0.151x10-4rad
El momento en la articulación A es nulo, por tanto:
6𝐸𝐼 4𝐸𝐼 ∗ 6𝐸𝐼 ′ 2𝐸𝐼 ′
𝑣𝐴 + 𝜗𝐴 − 𝑣 𝐵 + 𝜗 𝐵 = 0 ⇒ 𝜗𝐴∗ = 0.825𝑥10−4 𝑟𝑎𝑑
𝐿2𝐴𝐵 𝐿𝐴𝐵 𝐿2𝐴𝐵 𝐿𝐴𝐵
Por tanto el giro en A, es:
50 𝐿𝐴𝐵 50𝑥2.5
𝜗𝐴 = + 𝜗𝐴∗ = ( + 0.825) 10−4 = 7.769𝑥10−4 𝑟𝑎𝑑
4𝑥𝐸𝐼 4𝑥4.5
EJEMPLO 4.9
Calcular el desplazamiento en el nudo E de la estructura de la figura 4E9.1,
considerando la simetría de la misma.
Datos: P= 40 KNw
Para todas las barras el módulo de elasticidad es E=2 x 104MPa.
Para las barras AB, CB y EF, el área de la sección es A=0.15 m2 y el momento de
inercia I=30 x 10-4m4
Para las barras AC, CE, BD y DF, el área de la sección es A=0.30 m2 y el momento
de inercia I=60 x 10-4m4
14
P
A C
E
4m
B D F
4m 4m
FIGURA 4E9.1
= +
D
B F
(1) (2)
FIGURA 4E9.2
Estado (1):
La barra AB tiene un esfuerzo de tracción de valor P/2.
La barra CD tiene un esfuerzo de compresión de valor P.
La barra EF tiene un esfuerzo de compresión de valor P/2.
Estado (2):
Se trata de un estado simétrico respecto al eje vertical que pasa por C y D, que se
puede simplificar tal y como se puede ver en el primer paso presentado en la figura 4E9.3.
En este caso el esfuerzo axil en la barra CD es nulo, lo que se cumple al considerar los
empotramientos en C y D, siendo además correctos los desplazamientos relativos entre
los nudos de la estructura. Por otra parte, esto permite evitar la singularidad de la matriz
de coeficientes de rigidez de la ecuación matricial de la estructura que resultaría al
considerar deslizaderas verticales en C y D (por supuesto con cargas de valor P/2 hacia
arriba aplicadas sobre ambos nudos). La estructura obtenida es simétrica con carga
antisimétrica respecto al eje horizontal que pasa por su centro y se puede simplificar tal y
como se muestra en el segundo paso presentado en la figura 4E9.3.
B’
y a
P/2 x
P/2 P P/2
15
FIGURA 4E9.3
Si se consideran como sentidos positivos de los ejes locales los B’A (α=π/2) y AC
(α=0), la ecuación matricial para resolver el problema es:
𝑢𝐵′ = 0 𝑋𝐵′
𝑣𝐵′ 0
𝐵′𝐴 𝐵′𝐴
𝐾11 𝐾12 𝜗𝐵′ 0
[ 𝐵′𝐴 𝐵′𝐴 𝐴𝐶 ] = = (4E9.1)
𝐾21 𝐾11 + 𝐾11 𝑢𝐴 0
𝑣𝐴 −2 10−2
{ 𝜗𝐴 } { 0 }
Barra B’A:
𝑐𝑜𝑠 𝛼 −𝑠𝑒𝑛 𝛼 0
′ 𝑻
𝑲=𝑳𝑲 𝑳 ; 𝑳 = [𝑠𝑒𝑛 𝛼 𝑐𝑜𝑠 𝛼 0] ; 𝛼 = 𝜋/2
0 0 1
𝐸𝐴 𝐸𝐴
0 0 − 0 0
𝐿 𝐿
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
0 0 −
𝐿3 𝐿2 𝐿3 𝐿2
6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼
0 0 − 2
′
𝑲 = 𝐿2 𝐿 𝐿 𝐿
𝐸𝐴
0 0
𝐿
12𝐸𝐼 6𝐸𝐼
𝑆𝐼𝑀 0 −
𝐿3 𝐿2
6𝐸𝐼 4𝐸𝐼
[ 0 − 2
𝐿 𝐿 ]
90 0 −90
𝐵′𝐴
𝑲11 =[ 0 1500 0 ]
−90 0 120
−90 0 −90
𝐵′𝐴
𝐾12 =[ 0 −1500 0 ]
90 0 60
90 0 90
𝑲𝑩′𝑨
𝟐𝟐 =[0 150000 0 ]
90 0 120
Barra AC:
1500 0 0
𝑲𝑩′𝑨
𝟐𝟐 = [ 0 22.5 45 ]
0 45 120
16
90 0 −90 −90 0 −90 𝑢𝐵′ = 0 𝑋𝐵′
0 1500 0 0 −1500 0 𝑣 𝐵′ 0
−90 0 120 90 0 60 = 𝜗𝐵′ 0
= (4E9.2)
−90 0 90 1590 0 90 𝑢𝐴 0
0 −1500 0 0 1522.5 45 𝑣𝐴 −2 10−2
[−90 0 60 90 45 240] { 𝜗𝐴 } { 0 }
FIGURA 4E9.4
17
Las componentes radial, circunferencial y axial de una carga son:
∞ ∞
z ϑ r x
FIGURA 4.7.7
EJEMPLO 4.10
Dibujar las componentes radial axial y circunferencial del desarrollo en serie de
Fourier, para n=1, indicando si si son simétricas o antisimétricas respecto del plano ϑ=0
ya que el problema axilsimétrico se plantea sobre el plano xy (ver figura 4.7.7. Las
componentes de la fuerza se considerarán constantes e igual a uno, es decir:
Fr0 = Fy0 = Fϑ0 = Fr1 = F’r1 = Fy1 = F’y1 = Fϑ1 = F’ϑ1 = 1
Para n = 0 :
y y y
z x z x z x
FIGURA 4E10.1
Para n = 1 :
18
y y y
x x x
y y y
x x x
FIGURA 4E10.2
19