0% encontró este documento útil (0 votos)
4 vistas19 páginas

Simetrias

El documento aborda el análisis de estructuras simétricas y sus simplificaciones en el cálculo de esfuerzos y deformaciones. Se detalla cómo las cargas pueden ser simétricas, antisimétricas o arbitrarias, y cómo estas afectan los resultados en estructuras planas. Además, se presentan consideraciones sobre los efectos estáticos y cinemáticos en función de la simetría de la carga aplicada.

Cargado por

juantelo.09
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
4 vistas19 páginas

Simetrias

El documento aborda el análisis de estructuras simétricas y sus simplificaciones en el cálculo de esfuerzos y deformaciones. Se detalla cómo las cargas pueden ser simétricas, antisimétricas o arbitrarias, y cómo estas afectan los resultados en estructuras planas. Además, se presentan consideraciones sobre los efectos estáticos y cinemáticos en función de la simetría de la carga aplicada.

Cargado por

juantelo.09
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

4.

7 ESTRUCTURAS SIMÉTRICAS

Es relativamente frecuente que las estructuras presenten simetrías


respecto a uno o varios planos, lo que permitirá buscar, aprovechando esta
circunstancia, simplificaciones para su cálculo. Incluso resulta también habitual
que se presenten simetrías en la carga, como por ejemplo puede ocurrir con la
carga de peso propio en una estructura simétrica, pero aunque esto no resulte tan
habitual, una carga cualquiera se podrá tratar mediante la superposición de un
estado con la simétrica y otro con la carga antisimétrica, tal y como como se verá
más adelante. El hecho pues de que una estructura simétrica con carga simétrica
o antisimétrica, permita conocer a priori los valores de los esfuerzos y
deformaciones en algunos puntos de la misma, permitirá simplificar su cálculo y
de ahí el enorme interés que tiene el estudio de este tema. Debido al carácter
básico del curso, en adelante se centrará el análisis en el caso de estructuras
planas, en las que los planos de simetría serán ejes (intersección del plano de la
estructura con el plano de simetría).

4.7.1 ESTRUCTURAS PLANAS SIMÉTRICAS

Una estructura plana es simétrica geométricamente respecto a un eje, si


son simétricas respecto a dicho eje, la forma de la misma, las dimensiones de las
secciones, las propiedades elásticas del material y las condiciones de contorno
(cuestión que habrá de matizarse siguiendo el ejemplo 4.9). A su vez, la carga que
actúa sobre una estructura simétrica, podrá ser simétrica, es decir una imagen
especular respecto al eje de simetría (figura 4.7.1.a), antisimétrica, al invertir el
sentido de la parte simétrica (figura 4.7.1.b) o, como será más usual, arbitraria,
que siempre se podrá estudiar mediante la superposición de dos estados, uno
con carga simétrica y otro con la carga antisimétrica (figura 4.7.1.c).
Por otra parte, una carga P aplicada en la misma dirección del eje de simetría,
se considera equivalente a dos cargas de valor P/2 con el mismo sentido y
situadas a ambos lados del eje de simetría separadas una distancia infinitesimal,
lo que constituye un estado de carga simétrico. Si la dirección de una carga con
su punto de aplicación situado en el eje de simetría es perpendicular a dicho eje,
se considerará equivalente a dos cargas P/2 con el mismo sentido y situadas a
ambos lados del eje de simetría situadas a una distancia infinitesimal, lo que
constituye un estado de carga antisimétrico. Un momento aplicado sobre el eje
de simetría, constituye un estado de carga antisimétrico con dos momentos M/2
situados a una distancia infinitesimal a ambos lados del eje de simetría.
Por último, se podrá dar el caso de que una estructura plana presente una
doble simetría, tal y como se puede ver en la figura 4E7.1.

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

q q/2 q/2 q/2 q/2


M
M/2 M/2 M/2 M/2

∆T2 ∆T1 ∆T1 /2 ∆T2 /2 ∆T2 /2 ∆T1 /2 -∆T1 /2 -∆T2 /2 ∆T2 /2 ∆T1 /2

P P/2 P/2 P/2 P/2

= +

c.- Estructura simétrica con carga general Estado simétrico Estado antisimétrico

FIGURA 4.7.1

4.7.2 CONSIDERACIÓN DE LA SIMETRÍA

En una estructura geométricamente simétrica con carga simétrica, los efectos


estáticos, es decir las reacciones y esfuerzos, así como los cinemáticos (los
desplazamientos y deformaciones) deben ser simétricos, mientras que si la carga
es antisimétrica, deberán ser antisimétricos. Esto supone que con carga simétrica
los diagramas de esfuerzos axiles y momentos flectores, serán simétricos y los de
esfuerzos cortantes para ser simétricos deben tener signo contrario y, por tanto,
su diagrama será antisimétrico (figura 4.7.2.a ). En el caso de carga antisimétrica
ocurrirá lo contrario, tal y como se puede ver en la figura 4.7.2.b.

2
= = = =

a.- Diagramas de esfuerzos y deformada en b.- Diagramas de esfuerzos y deformada en


una viga simétrica con carga simétrica una viga simétrica con carga antisimétrica

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

Esfuerzo cortante (V) V=0 V≠0

Esfuerzo axial (N) N≠0 N=0

Desplazamiento horizontal (u) u=0 u≠0

Desplazamiento vertical (v) v≠0 v=0

Giro ϑ ϑ=0 ϑ≠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).

M=0 A/2 N=0 I/2


A
V=0
I

= =

a.-Estructura real b.-Estructura simétrica b.-Estructura simétrica


con carga simétrica con carga antisimétrica

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 + 𝛽𝑗𝑖 ) 𝐿

Por simetría de la viga:

k 𝑖𝑗 = k𝑗𝑖 = 𝑘; β𝑖𝑗 = β𝑗𝑖 = 𝛽 (4.7.2)


Para el caso de que sobre la viga actúe una carga simétrica, de acuerdo con las
consideraciones realizadas anteriormente, las condiciones que se deben cumplir son:
ϑ𝑖𝑗 = −ϑ𝑗𝑖 ; u𝑖 = −u𝑗 ; v𝑖 = v𝑗 (δ𝑖𝑗 = 0)
0 (6.2.3)
M𝑖𝑗 = −M𝑗𝑖0
Por tanto, al sustituir las condiciones expresadas en (4.7.2) y (4.7.3) en la primera
ecuación de (4.7.1), se tiene:
Mij = Mij0 + k(1 − β)ϑij (4.7.4)
Es decir que la rigidez al giro que se debe considerar al aplicar la simetría, es:
k′𝑖𝑗 = k(1 − β) (4.7.5)
que para el caso habitual de piezas de sección constante (𝛽 = 1⁄2)esta rigidez es
𝑘
k′𝑖𝑗 = 2
Para el caso de carga antisimétrica, las condiciones que se deben cumplir son:
ϑ𝑖𝑗 = ϑ𝑗𝑖 ; u𝑖 = u𝑗 ; v𝑖 = −v𝑗 (δ𝑖𝑗 = 2δ𝑖 ) (4.7.6)
Al sustituir las condiciones (4.7.2) y (4.7.4) en la primera ecuación de (4.7.1), se
tiene:
δ
Mij = Mij0 + + k(1 + β)ϑij − 2k(1 + β) Li (4.7.7)
que para el caso de piezas de sección constante (𝛽 = 1⁄2) supone una rigidez al giro
3𝑘 6𝐸𝐼
(segundo término del segundo miembro de (4.7.7)) 𝑘′𝑖𝑗 = = y una rigidez
2 𝐿
respecto al desplazamiento perpendicular a la directriz de la viga (tercer término del
12𝐸𝐼
segundo miembro de (4.7.7)) 𝑘′𝑖𝑗 = 2𝑘 = 𝐿3 .
vi vj
θi ui θj uj
i j

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

Se trata de una estructura de nudos articulados simétrica con carga simétrica


situada en el eje de simetría (eje vertical que pasa por el nudo B). De acuerdo con lo
indicado anteriormente, al aplicar las condiciones de simetría, la estructura simplificada
será la indicada en la figura 4E6.2, donde además se ha eliminado la barra AD ya que no
trabaja al estar entre dos nudos fijos, se han indicado con flechas las direcciones y sentidos
positivos de los ejes x locales de cada una de las barras y se han indicado los ejes globales.

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:

𝑲 = [𝒌𝑎𝟐𝟐 + 𝒌𝑏𝟏𝟏 ] (4E7.1)

Siendo las matrices elementales en coordenadas globales:


𝐸𝐴
0 0
𝐿 𝐸𝐴
3𝐸𝐼 3𝐸𝐼 0 0
𝒌𝑎22 = 0 − 2 = [ 𝐿 ]
𝐿3 𝐿 0 150 −150
3𝐸𝐼 3𝐸𝐼 0 −150 150
[ 0 − 𝐿2 𝐿 ]

𝐸𝐴
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 𝜗 𝐴

Al haber solo un giro no nulo el vector de desplazamientos en locales y globales es


el mismo, por tanto:
𝐸𝐴
𝑋𝐶 0 0 0 0
𝐿
{ 𝑌𝐶 } = [ 0 −600 −300 ] { 0 } = { 0.8928𝑘𝑁 }
𝑀𝐶 0 300 100 𝜗𝐴 −0.2976 𝑘𝑁𝑚
El valor de la reacción es: (en ejes locales)
𝑅𝑥𝐶 = 0
𝑅𝑦𝐶 = −0.8928 𝑘𝑁
𝑀𝐶 = −0.2976 𝑘𝑁𝑚
Debidos a la carga: (en ejes globales)
5 𝐿2 5
𝑀𝐶0 = − =− = −0.417 𝑘𝑁𝑚
12 12
0
5𝐿2
𝑅𝑥𝐶 = = −2.5 𝑘𝑁
2
La reacción en C, es:
RxC = 2 (-0.8928-2.5) = -6.7856 kN
MC = 2 (-0.2976-0.417) = -1.4252 kNm

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

1.5m 1.5m 1.5m 1.5m

FIGURA 4E8.1

El cálculo se plantea mediante la superposición de dos estados uno simétrico y otro


antisimétrico, tal y como se puede ver en la figura 4E8.2. En dicha figura ya se hay
incluido las simplificaciones por simetría y antisimetría de la carga, así como los ejes x
locales de cada una de las barras y los ejes globales.
10 10 C’
B C B B C’

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:

Las cargas debido al momento sobre la barra BC’, son:

𝑝 𝐿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.

Siendo las matrices elementales en coordenadas globales:

- Barra BD (casos simétrico y antisimétrico): senα = 0.8; cosα = 0.6

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

- Barra AB: senα = 0.8; cosα = -0.6


∗ 𝐶𝑎𝑠𝑜 𝑠𝑖𝑚é𝑡𝑟𝑖𝑐𝑜 𝑦 𝑠𝑒𝑔𝑢𝑛𝑑𝑜 𝑝𝑙𝑎𝑛𝑡𝑎𝑚𝑖𝑒𝑛𝑡𝑜 𝑑𝑒𝑙 𝑎𝑛𝑡𝑖𝑠𝑖𝑚é𝑡𝑟𝑖𝑐𝑜:
𝐸𝐴
0 0
𝐿
12𝐸𝐼 6𝐸𝐼 25.252 −29.061 3.456
𝒌𝐵𝐴
22 = 𝑳
𝑇
0 − 2 𝑳 = 104
[ −29.061 42.204 2.592]
𝐿3 𝐿 3.456 2.592 7.2
6𝐸𝐼 4𝐸𝐼
[0 − 2
𝐿 𝐿 ]
𝐸𝐴
0 0
𝐿
12𝐸𝐼 6𝐸𝐼 25.252 −29.061 −3.456
𝐵𝐴
𝒌11 = 𝑳𝑇 0 𝑳 = 104
[ −29.061 42.204 −2.592]
𝐿3 𝐿2 −3.456 −2.592 7.2
6𝐸𝐼 4𝐸𝐼
[0 𝐿2 𝐿 ]
𝐸𝐴
− 0 0
𝐿
12𝐸𝐼 6𝐸𝐼 −25.252 29.061 3.456
𝐵𝐴
𝒌12 = 𝑳𝑇 0 − 𝑳 = 104
[ 29.061 −42.204 2.592]
𝐿3 𝐿2 −3.456 −2.592 3.6
6𝐸𝐼 4𝐸𝐼
[ 0 − 2
𝐿 𝐿 ]

∗ 𝐶𝑎𝑠𝑜 𝑎𝑛𝑡𝑖𝑠𝑖𝑚é𝑡𝑟𝑖𝑐𝑜:
𝐸𝐴
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
𝐿 𝐿 ]

- Barra BC’ (casos simétrico y antisimétrico):

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 ]

Sustituyendo valores, se obtienenlas ecuaciones matriciales que se incluyen a


continuación.
- Caso simétrico:
Al sustituir en 4E8.1, se obtiene:

155.512 1.244 5.184 −106.667 0 0 𝑢𝐵 0


1.244 99.475 13.296 0 −16 12 𝑣𝐵 −15
0 −12 6 𝜗𝐵 −7.5
104 5.184 13.296 24.6 = 𝑋𝐶′
−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:


uB= 0.0171x10-4m; vB=-0.173x10-4m; ϑB=-0.472 x10- 4rad; vC’=-0.5277 x10-4m

- 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

Con el planteamiento dado en 4E8.3, se obtiene:


7.2 3.456 2.592 3.6 0 0 𝜗𝐴 50
3.456 155.512 1.244 5.184 −106.667 0 𝑢𝐵 0
99.475 13.296 0 12 𝑣
104 2.592 1.244 𝐵
= 0
3.6 5.184 13.296 24.6 0 6 𝜗𝐵 0
0 −106.667 0 0 106.667 0 𝑢𝐶′ 0
[ 0 0 12 6 0 12 ] { 𝜗𝐶′ } { 0}
Obteniéndose los mismos resultados que con el anterior planteamiento.

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

Se calculan las reacciones y se plantea la estructura mediante la superposición de


dos estados tal y como se puede ver en la figura 4E9.2.
P P/2 P P/2 P/2 P P/2
A C
E

= +
D
B F

P 2P P/2 P P/2 P/2 P P/2

(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.

P/2 P P/2 P/2 A P/2


C
SIMETRÍA ANTISIMETRÍA
EJE VERTICAL EJE HORIZONTAL a/2
H

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

La ecuación matricial (4E9.1), resulta:

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 }

Que da como resultados:

𝑢𝐴 = −9.96 10−6 𝑚; 𝑣𝐴 = 𝑣𝐵′ = 0.00156 𝑚; 𝜗𝐴 = 0.000337 𝑟𝑎𝑑;


𝜗𝐵′ = −0.00015 𝑟𝑎𝑑

FIGURA 4E9.4

Como en la estructura real (figura 4E9.1) el desplazamiento vertical de B y D es


cero, para obtener el movimiento real de la estructura en E, habrá que girar la estructura
sobre D, un ángulo:
𝛿
𝑡𝑎𝑔 𝛼 ≈ 𝛼 = 𝑎 (4E9.3)
Esto supone que el punto E desciende:
𝛿𝐸 = 𝛿 + 𝛼 𝑎 = 2𝛿 = 0.0031 𝑚 (4E9.4)

4.7.3 SIMETRÍA AXIAL

Existen muchas estructuras, como torres de refrigeración, tanques o depósitos, en


los que la geometría y las propiedades mecánicas son independientes de la coordenada
circunferencial (figura 4.7.7). Si las cargas exteriores son también de revolución, su
cálculo se puede simplificar ya que se puede abordar de forma bidimensional en las
direcciones radial y axial, resultando su formulación fuera del alcance de este texto. No
obstante puede resultar de interés en este momento lo que se indica a continuación. Si las
cargas no son simétricas, para evitar un análisis tridimensional se puede hacer un
desarrollo en serie de Fourier en dirección circunferencial de la misma, hacer un análisis
bidimensional (en el plano x,y) para cada uno de los términos del desarrollo
separadamente, y a continuación una superposición de los resultados obtenidos.
Normalmente no se precisa más que hacer el cálculo con unos pocos harmónicos (diez a
lo sumo) para alcanzar una buena aproximación.

17
Las componentes radial, circunferencial y axial de una carga son:

∞ ∞

𝐹𝑟 = ∑ 𝐹𝑟𝑛 (cos 𝑛𝜗) + ∑ 𝐹′𝑟𝑛 (sen 𝑛𝜗)


𝑛=0 𝑛=0
∞ ∞

𝐹𝑦 = ∑ 𝐹𝑦𝑛 (cos 𝑛𝜗) + ∑ 𝐹′𝑦𝑛 (sen 𝑛𝜗)


𝑛=0 𝑛=0
∞ ∞

𝐹𝜗 = ∑ 𝐹𝜗𝑛 (cos 𝑛𝜗) + ∑ 𝐹′𝜗𝑛 (sen 𝑛𝜗)


𝑛=0 𝑛=0

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

Radial Fr simétrica Axial Fy simétrica Circunferencial Fϑ antisimétrica

FIGURA 4E10.1
Para n = 1 :

18
y y y

x x x

Radial Fr simétrica Axial Fy simétrica Circunferencial Fϑ antisimétrica

y y y

x x x

Radial F’r simétrica Axial F’y simétrica Circunferencial F’ϑ antisimétrica

FIGURA 4E10.2

19

También podría gustarte