UNI-VERSIDAD DE SANTANDE R

- ESCUELA TECNJCA SUPERIOR DE INGENIEROS -

- - - DE CAMINOS, CANALES Y PUERTOS - -

DEPARTAMUTO DE ANALISIS DE LAS ESTRUCTURAS
tATfDRA DE CALCULO DE ESTIUCTURAS

CALCUL O MATRICI AL DE

ESTRUC TURAS

problemas

AVELINO SAMARTIN QUIROGA

J.R. GONZALEZ DE CANGAS
LUIS MORENO GARCIA
PUBLICACION : AE - 80.1 JAVIER TORRES RUIZ

1980
 indice

1. Planteamiento general.
2. Matriz de rgidez de elementos.

3. Matriz de rgidez de estructuras.

4. Clculo de estructuras. Condiciones de borde.

5. Problemas especiales. Nudos de dimensi6n finita.

6. Clculo del grado de h i pe res ta t i s mo

7 . Sistemas bsicos.
8. Cambio de bases.
9 . Elementos finitos triangulares y cuadrangulares eo
1 o. Elementos finitos triangulares y cuadrqngulares e .
 presentacin

El Clculo matricial de estructuras corresponde a un

planteamiento moderno y original-en su acepci5n mis genuina
del tE!"rmino, de vuelta a los orgenes- del Analisis de Estruc
turas. Su desarrollo y por lo tanto aprendizaje ms eficaz, se
realiza con. la ayuda de un computador. Sin embargo, re.legar al
estudiante al simple papel de usuario de programas generales
de clculo de estructur~s por computador, que puede utilizar
1

como una misteriosa caja negra, sin comprensi5n de sus funda-

mentos, parece una penosa y poco formativa perspectiva educa-
cional. Por otra parte, el desarrollo por el alumno de peque-
os programas de ciculo, permite una adecuada y natural asi-
milaci6n de algunos conceptos del clculo.matrlcial de estruc
turas, si bien exige un previo conocimiento de procedimientos
nmericos .e informativos, y particularmente, un esfuerzo, a
veces, desproporcionado para un estudiante que s5lo desea com
prender pero no d~sarrollar nuevos programas de cl~ulo. En -
esta publicaci5n, se intenta un camino intermedio, en donde
se da especial ~nfasis a los conceptos fundamentales de matriz
de rigidez y cargas equivalentes de elementos y subestructu--
ras. Conceptos que son ll'tiles, y prc!cticos, tanto para el fu-
t u r o u s u a r i o ' e o mo a 1 r e a l. i z a d o r d e p r o g r a ma s Es t e h e e _h o s e
puede comprobar en los modernos programas generales donde exis-
ten facilidades de subestructuraci5n a diferentes niveles. Ade-
ms se incluyen aqu simples problemas de anlisis matricial de
estructuras que se plantean de un modo manual y que no precisan
para su resoluci6n necesariamente un computador. Por ltimo, se
muestran ejemplos de uso de un programa general de computador -
implmentado 'en.el Centro de Clculo de la U,niversidad de San--
tander.
Con esta publicaci5n c~eemos, que se intenta cubrir
un importante vaco en el campo del anlisis matricial de ~s

tructuras,: en donde existen excelentes 1 ibros de exposici5n

te5rica de losmtodos -flexibilidad, rigidez, mixtos, hbri-
dos, etc.- pero que carecen de esa intencin de acercar el
frecuentemente distante planteamiento general matricial al a-
lumno.
 Los problemas presentados son fruto, de la experie~

cia docent en la Catedra de Clculo d~ Estructuras, Y. corre~

ponden a prcticas, exmen~s o explicaciones en clases. No s&
han tratado, dada la 1 imitacin del curso otros mtodos dis--
tintos al mtodo del equilibrio o de los movimientos. Sin em-
bargo se esperan introdu~ir algunos ejemplos, de estos mto--
dos. Esta publicacin, corresponde a una continuacin del es-
fuerzo que la Ctedra esta realizando con objeto de poner a -
disposi~in de los alumnos de. la f.T.s. de lngenferos de Caml
nos, Canales y Puertos de Santander de aqu~11as ayudas que --
les permitan una ms feil y eficaz asimilacin de los canee~
tos mis importantes del Clculo matricial de.estructuras. Ese
ha sido nuestro propsito, y nuestra mejor recompensa serTa -
el haberlo lograda. En cualquier caso, se desea, y~ desde aho
ra, recoger en futuras edicio~es, aquellas posibles correccio
nes y mejora~ que a lo largo de los sucesivos cursos aparez--
can.

A. Samartfn.
 - 4 -

Ejercicio E 1.1.
Numerar las estructuras de la ffgura E l.la de mo-
que el ancho de banda sea mTnimo. Se supone estructuras de
barras.

(1 )

F ig. E 1. 1 a
 - 5 -

Si A-1 es la mayor. diferencia

. entrri dos nudos conti-
guos en cada una de las estructuras se tiene lo que se indi-
ca en 1 a f i g u r a E1 1 b

S 9 13 17 21 25 3
2
1 1

6 10 , 1 18 22 2
2 ---

1 AmS 13 7 11 15 19 23 27 4

4 8 12 16 20 24 28

(1 ) ( 2)

3
nL
6 9 10 11 14

18

19 22
. 1
12
11
10 . 13~~15

18
16

,/l;}j 9 13

' / _/,21

23. /
1 A: t. 1 =S

.
V/ (3

24

A . 3,::7
1

(:g E',,_
.. ' ~)

El ancho de banda par.a cada una de las estructuras

 - 6 -

anteriores es,

Estructura 1 : (2x5n-1= lOn-1)

. ll z: (2x3n-1= 6m-1)
11
3: (2x7n-1= 14n-1)
11
4: (2x5n-1= lOn-1)
siendo n el nmero de grados de 1 i bertad por nudo.

~esde el punto de vista de computacin, se aprove-

cha el hecho de que la matriz es simtrica, stendo preciso
en este caso cons)derar el semiancho de banda a efectos de
almacenamiento en la memoria del computador. Es decir el se-
miancho de banda es An.
 - 7 -

Ejercicio E 1.2.

-
Determinar el ancho de banda real (nOmero de diago
nales paralelas a la principal con trminos no nulos en la
matriz de rigidez) de la estructura de la figura E 1.2a.

1
\ (f)
1
1
1

Ftg. E i. 2a
 -8 -

Se considera en cada barra los sentidos indicados

en la figura E1.2b.

Fig. E 1. 2b
 - 9-

Con este supuesto se tiene la matriz siguiente:

1 2 3 4 5

(K22)a+ .
1
(K11)b+ o (K12)h (K12)b o
(K11)h
(K22)f+
o (K22)c+ o (K21)c (K12)g
( K11) g
..
(K22)h+
3
(K21)h o (Kll)i+ o (K 2) i
.
~

(K22) j
(K22)b+ .
(K22)d+

4 o (K22)e+
( K2 1 ) b (K12)c (K12)m
(Kll)c+

(K11 )m
(K22)g+
. (K22)m+
S o (K21)g ( K2 1 ) i ( K2 1) m (K22)i+
(K22)k+
.
( K2 2) 1

Suponiendo que hay 6 grados de 1 lbertad por nudo

(tres desplazamientos y tres giros). Se tiene la siguiente

matriz:
 - 10 -

23
...

ro-
. oooooo.
2 3 4 1 5
X lOOI.llll .XXX XXX X ltXXX X ~~o ooo
X o o -gooo
X o oo too
1 ll o ooo ~o
X o ooo o D
X o oooooo_.,
o XXX XXX
o
2 o
o
o
o
24 lllUllllll
Jt
l
X
lt
X
XXX XXX
1
llllllliXll
4 XXliXlllC
XXXXlllt
xxxxxx
"--- X 11 lt ltl.ll
!'.o o ooa
o' -
S
o~
o
' ' '".,,
>

La diferenci.a m~xima entre nudos es:.{ 5-2

_ ,.. 3} luego
4 1 3
A= 4, pero el ancho ser1: ~(~amero de gra~os de 1 ibertad por
nudo)= 4.6= 24. Ahora bien, ese sera el ancho reducido, es-
to es, si la matriz fuese simtrica. El ancho total es
2x24-1= 48-1= 47. (Se le resta 1 para no contar dos veces la
diagonal principal).
As pues la frmula es 2.A.n-1., donde A-1 es la
diferencia m(xima entre dos nudos y n el nQmero de grados
de libertad por nudo (gdl).
 - 11 -

Ejercicio E1.3.
Numerar los nudos de la estructura de la figura
E1.3a., de modo que el ancho de banda de la matriz de ri-

-
gidez total sea mnimo. Si N es el nGmero de ~rados de 1 i
bertad por nudo de la estructura, determinar el ancho de
banda anterior en los dos casos siguientes:
a) Estructura de barras.
b) Estructura de elementos finitos cuadrilteros.

Fig. E1. 3 a
 - 12 -

a) 6
- Ancho de banda 2.6.N-1
b) 7
- Ancho de banda 2.7.N-1
 - 13 -

Ejercicio E2.1.
Determinar la matriz de rigidez de la barra rectd
de la figura E2.1a. que tiene area variable linealmente en-
tre sus extremos y est sometida nicamente a esfuerzos
axiles (barra de estructura celosa).- Efectuar el clculo
mediante:
a) -Mtodo exacto.
b) -Ecuacin diferen-
cia 1.
e) -Elementos finitos
con funciones de

t X
L

Fg. E 2. 1a
.1 forma lneales
d) -Comparaci6n de
los mtodos para
A.= 1,10 y 100.

Caractersticas: E A2. =A.

'~
 - 14 -

a) Mtodo general~
Se considera la figura E2.1b.

-ti.-f
1~--~-------------------
J ~
. S

Fig; E2. lb

La energa_elstica de deformacin es:

U= .!_L N2 ds
2 o "EA
N= S (cons.tan.te) .,

s=x
2
Por lo tanto u~
s
----
L dx
A-l
ln A
2EA1 o 1+ L X --r:T

Segn el teorema de Castigl.iano, se obtiene:

au
a= ~S , es decir a= L
---EA ln A S
a 1 --r:T

La matriz de flexibilidad es:

ln A -1 A -1
F= --r:T ; F ... K... TiiT

La matriz de equilibrio es H=1, por lo que la

matriz de rigidez del elemento est dada en la frmula:
 - 15 -

b) Ecuacin diferencial.
La ecuacin gneral de la extensin de una barra
es:

d
dX {EA ~}
dx + pz O {2.1a.}

con las condiciones de contorno pertinentes en x=O yx=l

estticas (EA*= valor especificado) cinemticas (u= va
lor esficificado).
En este caso se tiene:
- En la viga, la carga actuante es p=O.
-Las ecuciones de borde son:
en x=O EA~=
lax
en x=l EA~=
2ax

La solucin general de la ecuacin {2.1a.} es:

- Adu
E Tx .. e,
l+lf:l X

es decir
A-1
u= EA (A _ 11
1
{e 1 1 n ( 1+ -r- x) + e2 }
1

Siendo e 1 y c 2 dos constantes arbitrarias.

La consideracin de las condiciones cinemticas con
ducen al siguiente sistema:

L''
= EA l (X -1 ) {~nA
En forma compacta se puede escribir
d= G e
-:d Las condiciones estgticas llevan a la ecuacin:

o bien matricialmente
2.= - e
Gp -
 - 16 -

Eliminando las constante~ c 1 y c (o el vector ~)

2
entre las ecuaciones anteriores s~ obti~ne la matriz de r-
9 i dez._
Es decir:
.e.= --
k d
G
con k= -p
-1
- i;
EA 1
con Gd =
-1
- L
A.-1
T'A{ lnk
-t 1
o}

Por lo tanto
!;Al A.-1 1
Ji= L TA {_,

e) M'todo de los elimentos finitos.

Se adoptan coordenadas adicionales ~= fy las fun
clones de forma lTneales. Es deci~
d
u(x)= (N 1 , N ) .(dl)
2
2

La relacin deformacin-desplazamientos es
d
e(x)= dx {u(x)}

es decir
e= (8p 8 )
2
(~1)
2
aN 1 aN
con B= ax- = Lar-
En este caso N1= --
1
L
N
2
= -.
L

La matriz de rgidez es

k . . = !L 8: D 8 . d x (i, j=l, 2)
1J o 1 J

con D= EA(x)= EA {1+(A-1)~}que es la matriz constitutiva del

1
materia 1
 - 17 -

Resulta, entonces

La matriz de rigidez es entonces

EA 1 1+A. { 1
k= -,---r -1

d) La comparacin de ambas matrices se resume en la siguien-

te tabla de valores de k, siendo
1
k= k - 1 } , que indica la influencia de la varia--
e i n de 1a se ce i n.

Tabla E2.1, Valor-es de k

Mtodo exacto elementos finitos

A.=1_,0
(seccin constan
te) - 1, o 1,0

A.-=10 3,91 5,5

A.=100 21 , 5o 50,5

Es decir la aproximacin se deteriora a medida que

le variacin de la seccin es ms fuerte, a no ser que se --
utilizen otras funciones de formas mas ricas.
 - ~8 -

Ejercicio E2,2.
1. Determinar .las funciones de forma o de interpolaci6n
del elemento finito monodfmensional (1-D) de la fig~
ra E2.2a. 11
barra recta de tres nudos 11 sometido a un
esfuerzo axil, Se supone un desarrollo pol in6mico p~

ra el desplazamiento longitudinal u(x).

2. Dibcjar dichas funciones de forma.
3. Indicar sus propiedades ms importantes,

t=
1

L/2 _ _
3

Ir--.- -
2

L 12--1

Fig. E2. 2a

4. Conocidas las funciones de forma N. del elemento fi-

nito 1-0 11
barra recta de tres nudos' 11 aterior, deter
minar la matriz de rgidez ~ Hallar la matriz de r
gidez para el caso en que se consideren solamente --
dos nudos y la matriz de rgidez exacta.
5. Si actOa una carga longitudinal uniforme en toda la
longitud de la barra, se desea conocer la expresi6n
de las cargas concentradas equivalentes en los nudos
(soluci6n inicial).
6. Mediante la tcnica de condensacin esttica, obte-
ner los resultados de los dos apartados anteriores,
para la barra recta con dos nudos 1-2, tras la el i-
minaci6n del grado de 1 ibertad del nudo intermedio 3.
j. Comparar estos resultados con los obtenidos en un
clculo exacto sin la aproximaci6n del mtodo.de los
element9s finitos,
M6dulo de elasticidad: E.
Area de la secci5n Ax variable 1 inealmente de~de un
valor A1 en el nudo 1, al valor A2 en el nudo 2.
 - 19 -

1) Considerando la figura E2.2b. y que:

~
l:"a TX
Ax= ( 1+~~ ) A1 con ~- 'A2-A1
A

F1g. E2 2 b

y sabiendo que las funciones de forma son nulas en todos los

nudos menos en uno de ellos se tiene que las funciones de in
terpolacin son:

N1= (2~-1) (~-1)

2) Los gr~ficos de las ~unciones de forma son:

'~-
, 1
....L.
1
.
.
.
L ,-- ' 1
t_j_
-~~
..............7. . . 2 ....-...-.-_-..-:;A:--:tl"r;___2
1-x::- 3 2
__ g

Fig E2. 2 e

3) Las funciones de forma cumplen:

N +N +N
1 2 3
= 1

N (x.)
1 J
= o..
1J
 - 20 -

4) La relacin deformacin-desplazamientos es:

aN 1 'J
E= B u
-- con B= -
ax = !:.. r e;

s,= -(E;-
L
4
t) B2= -(t;-
L
4
*) B3= -(1-2t;)
L
4

La matriz de rgidez es:

1
k .. = E[L B. A B dx= LE! B. A B. dE;
1J o 1 X J o 1 X J

k=
L+
3
1
2 ll
1 + 1
3 'b ll - !8 - 3
2
ll EA
1
-:-
l.+ b1 ll
3 L3 11
+ b""ll - 38 - 2ll

-- -r8 2
ll - 38 - 2ll 16 + 8
3 r 7.1 ll

11
(1) -Matriz de rgidez del elemento 1-D barra recta cvn dos
nudos de seccin variable.

2 = E;
N=1-E; N
1

k= EA f 1 -1 } ( 1+~)
L" -1 1 2

La matriz de rgidez exacta es:

k EA 1 ll { 1 - 1} -
-exacta=-:- tn(1+l.!J -1 1

L
5) ( p1 ' p 2 , p 3) = f o (N 1 ' N 2 ' N 3) q d x=
l
q Lf { ( 2 t; - 1) ( t; - 1 ) ~ ( 2 ~ - 1) 4~ ( 1 - ~ ) } d E;... q L ( 1/6., - 1 /6 , 2/3 )
0

p 1 ... f., p2= f., 1

p3 = LcL!:.
3

( 1) ( p ( 1)
1 '
PJ1) )= 1
f o (N 1 N2)q dx= ql[ o {(1-t;) t;'J d <;1:" =
l

1.
= ql <-} 'i")

p ( 1) = .9..!:. s.h.
1 2 2

p ( 1 1) L -1 1
1exacta=. tn (1+..t):- "C" q ' p 2 exacta = (rn{l+lll + 1+.~)
,.. q.l-

 - 21 -

6) Condensacin esttica
4+31J
u= 4 + 11 ul + 4(2+ ) u2
3 ij{2+1J) . lJ

EA 1 { -1} (p 2 +61J+6)
1
k*= T -1. 1 3(2+1J)

Condensacin esttica:

p = .;Lb.
-1
* p1-k13 k33
P1= 3 o
+ 4+1J ~
4 (2+1J) 3 -.9.1..
-3
3+p
2+1J

p
~~~

2
. -1
= p2-k23 k33 p =
3
tl + 2i(2+1J)
4+3p ~
3
=.9...!.. 3+2u
3. 2+u

k k k ( 1)
7) Valores lJ aerox. exacto aprox.
o 1 1
1 1,444 1,443 1 , 5 00
1o 4,611 4, 1 7 o 6.000
100 34,660 21,668 51.000

( 1)
Valores P. aprox. P cxacta P 1 aprox.
lJ 1 1
o 0.500 0.500 0.5
0,444 0,443 0.5
1o 0,361 0,317 0.5
100 0,337 0,207 0.5
 - 22 -

Ejercicio E2.3.
Determinar la matriz de rgidez de la viga de sec-
cin constante y longitud L mediante la resolucin de la ecua
cin diferencial de la flexin. (Ver figura E2.3a.).

Ct 1)
F1g. E 2 3a
 ,.. 23 -

Se procede a part1r de la soluci6n tomplementaria,

como es usual en otros problemas estructurales,
Partiendo de la ecuaci6n diferencial de flexi6n:

la deformada es del tipo:

Por definici6n se tienen los movimientos en los ex-

tremos de la viga:
w 1=A
1
L8 =A
1 2
w2 =A 1+A +A +A 4
2 3
L8 2 =A +2A +3A 4
2 3

es decir
d= w, = o o r = A
ol 1 Al ..d
Le o o

~J
1 A2
1
w2 1 1 A3
L'e
2 o 2 3 A4
L

Anlogamente los esfuerzos correspondientes son:

Q,= iL 6A
L3 4
El
LM 1 = 3 2A3
L
Q
2
=--
El
L3
6A 4

LM = l l (2A +6A 4 )
2 L3 3
o bien
El
.e.= Q1 = o o o 6 Al = G A
LM
~ o o -2 o
-p
1 A2
Q2 o o o -6 A3
LM
2 o o 2 6 A4
 - 24 -

y 1a matriz de r1gidez es dire.ctamente

-1
.e.= k d con k= G ..d
-p

En este caso

-1
..d = r, o e o y k=
El
73
12 6 -12 6
o o o 6 4 - 6 2
-3 -2 3 -1 -12 -6 12 -6
2 -2 6 2 - 6 4

o bien S i se cons1. d era:E. 1 = Q1 y dl= w1

t-11 e,
~ w2
M2 62
L
1 kl d
con a=
1 El
1a matriz de rigidez es :k = 6L
12 6 L -12
L3
6L 4L 2 -6L 2L
-12 -6L 12 -6L
pL 2L -6L 4L 2
 - 25 -

Ejercicio E2.4.
Determinar ia matriz dt~ rgdez del .~iem.:1nto en-
tramado plano de la figura E2.4a., mediante !0: p~ocedi~ie~

tos siguientes:
a) Mtodo general matricial.
b) Ecuacin diferencial.
e) Condensacin esttica.
d) Mtodo de los elementos finitos en desplazamientos.
(Funciones de forma hermticas).
e) Comparar los resultados de este timo mtodo para
1=1 ; A=lO ; A=lOD.

r- l
r L
---4
,
~
...,1
... i
~
:3 L
-"--
l,
:; ,,
.J
!.2
t ../
1, Modulo d. P.lastcidad E

Fig. E 2. 4 a
 - 26 -

a) Mtodo general matricial.

Se considera la figura E2.4b.

Fig. E2. 4b

1 1
f-2-l] = [ H F- HT - HF_- .]
. L.2.z - F -: 1 HT . F- 1 J
-1 T -1
.!?..2 = -F H i 1+ F . iz y si se hace 1
~ = 0:

Teniendo en cuenta la figur E2.4c. se tiene:

,I. ,.._ __ __r=

__ 1 . 3

L~.)M-(x)-2) G2
X t
Q(x) Pz

F(g. E 2. 4 e

2 2
M(x)= G2+P 2 x 2() l.L(G2.+P2.x~ dx+ 2
J. L.(G2.+P2.x) dx
Q( x ) = p
1
} E= 2 f ~ l ( ~ J .d X= 2 o E 1 2. 2 1 E1 1
2
2
aE 6 2= L Gz+P2.x dx + L Gz+P2X dx
aG2. =. . ~ 11
o ""X 1 E1
1
 - 27 -

G2L PzL 2 G2L P2 2 .J_ 1 .L.2 1

8 2=rrr + i1ET +Er + 'i'ET" 3 L =
. 1 1 1 1
-rr--<r
1
1 )G2+rrr(r +3)Pz
1

G2 2 p2 3 G2 2 P2 3 L 2 1 L3 1
d2 2f-EI1 L+ 3Xf11.L + '2E'T1',3L + m;7L = '2E'Tl'(t'"3)G2+'3'"E'i,(;::-+"7)P2

82] L [(a+l) L
[dz E'i7" f(a+3) F=-
E 11

L2 L
E1 3(a+7) .-~{a+ 3) 12 E 1 L2
-1 1 12 1 -(a+?) ""f(a+ 3
F = -e- X 2
L {a: 2+14a+1} -j{a+ 3) a+1
=
L3 {a 2+14a+l} 3
L
2(a.+3) a.+ l

L2
r<7>-+1)
L
-i<3>-+ nl
--(3>-+1) f-+ 1 1
2
J
Se considera la figura E2.4d.

La matriz H es:

"'-
s"= P~
[ G J-H [P~
G J = fL-.2Pz
G 2LP 2J [1
=-o 2L
1
J f_
1L
__ Gp~l +H=[o1 2!.... ] ,HT=[zll
1
ol
11
~
 - 28 -

-1 12A.EI
1
2L 2
-HF = -r<>-+1) --t<>-+3)
x
L3{ 2 + 1 4>-+,}
t<3>-+ 1) -(A.+1.)

2
-l T 12A.EI 1 ll..(>.+ 1) ~(3>-+1)
-F H =3 . 3
L ( A. 2 + 1 4+ 1 ) : L .
-2(>-+3) -(A.+ 1)

2
_ T 12A.EI f-<>-+7) f(A.+ 3)
1 1
HF . H = !..3 (A. 2+ 1 4A+ 1 )
t< )..+ 3) )..+ 1

2 2
j-<>..+7) t<~+3) -4-<>..+ 1) -t(>..+3)

te>..+ 3) )..+1 t<3>..+1) -(A.+1)

12).. E 1
1
K= 2 2
- L3 (A. 2 + 14>..+ 1) 2 L .
3()..+1) f(3A.+ 1). t-<7A.+1) -~(3>-+1)
-~(/.,+3) -.(>..+ 1) -}( 3>..+ 1) A.+ 1

4 (>..+7) 6(>..+3) 8 (>..+ 1) -6(>..+3)

e L2 L L2
6(>..+3) 12 (A.+ U 6(3)..+1) -12(>..+1}
L2 L3 L2 L3
).. E 1
1
= 8 (>..+1) 6(3>..+1) 4(7>..+1) -6(31-+1)
>.. i+ 14A.+ 1 t L2 l2 L2
-6(>..+3) -12()..+1)-6(3)..~1) 12(1-+1)
L2 L3 L2 L3

b) Ecuacin diferencial.
Se considera la figura E2.4e.
La ecuacin diferencial de una viga recta a flexin
es:
 - 29 -

r- L
.1
L

-l
1

t T 1 .
2
=-=:f
)( 1 x
1

F1g. E2 <t e

con las condiciones de apoyo en xO , x=L:

- Cinemticas: w ~~ son valores especificados.
2
d w d 3w
- Estticas: -El --- El ---son valores especificados.
dx2 dx 3

En este ca~o, se tiene:

'E 1
1
= para el tramo 1-3

El
2
= par~ el tramo 3-2

La solucin general es:

Las condiciones de continuidad en 3 impl ica.n, para

esa seccin:
.d.WA dWB 2 : d 2W d 3W d 3 WB
-E1
d WA . - B - A
----r; ;;;r
WA =WB; 'd")( == dx 1 ' 1 :;7"'
= .El 2
dx 1
E1 1 =-E 1
2 crx;--
Es decir: A1=B 1 ; A~=-B
2; ~A
3 =B 3 ;aA 4 =-B 4
con a= r1 .
 - 30 -

Po r o t r a p a r te , 1 a s e o n d i e i o n e s e i n e m t i e a s de b o r -
de son:
JwA
x= L , WA= w,= d2= d yt -dx
-= e,= d1
'
dw_JJ
X ' =L;
1
we= w2 = d4= d y2 dx' = 82= d3

Resulta, entonces:
Ld 1=-A 2 -2A -3A4 o bi~n 1 en forma matricial:
3
d2= A1+A2+A3+A4
d= Gd.A
con-!~:{L; 1 , d , Ld , d }
2 3 4
Ld =-A +2aA -3aA 4 y
3 2 3 ~ ={ A1 , A2 , A , A4 }
3
d 4 = A' 1 -A 2+aA -aA 4
3

~
-1
1
-Gd= o -1
1 -1

An~logamente, las condiciones est&ticas de borde se

imponen como sigue:

Obteniendose: P
1
1'1
V o -2 -6
L P2 o o o -6
p3 o o 2 -6
L P4 o o o 6

3
Es.decir:
~,-[=
__i_p-.-~~ (E2.4b.) x....--i
De (E2.4a.) se deduce: 3 2
&a,.r= 1
~ l--y
y e n t r a n do en ( E2 4 b ) : E..= ~p Gd d
Fig. 2. 4 f
 - 31 -

Es decir: a= - -d,
K
1

lnversi6n de ~d: Gd= O -1 -2 -3

o - 1 2a. -3a.
1 -1 a. -a.

1
l
X X X X

X X X

1-5a. 3(1-a.) 5-cx 3(a.-1)

X
-(3a.+1) -2(a.+1) -(a.+ 3) 2(a.+1)

1
Como es a.=-
. A' resulta:

r X

X
X

X
X

X
X

X
Gd1 =--A._2_.... =
A5 3(A-1) 5A-1 3(1-A)
2
A +14A+1 ""'A X --x-: X
_(3+A) -2(1+/..) _ (1+3A) 2(1+A)
A A A X

X X X X

=---- A
X

A-5 3(A-1)
X

5A-1
X

3(1-A)
X

-(A+3) -~(A+l) - ( 3 A,- 1 ) 2(1+i\)

1
Por lo tanto, multiplicando Gp.Gd , se obtiene:
--
pl

LP2 +4(A+7) +6(A+3)

-
+8(/..+i) -6(/..+3)
l
/..E 1 . Ld
1
1 6(/..+3) 12(/..+1) . 6(3/..+1) -12(/..+1) . .d 2
p3 =
2
(/.. +14/..+l)L
2 +8 (/..+1) +6(3/..+1) +4(7.~+1) -6{3/..+1) Ld
LP.4 3
-6(/..+3) -12(/..+1) -6(3/..+1) +12 (/..+1) d4
 - 32 -

Queda, por tanto la matriz . de rigidez:

+ 4(A+7) + 6(A+3) +8(A+1) -6(A+3)

L L2 L L2
6(A+3) 12(A+1) 6 (3A+1) -12(A+1)
L2 :3 L2 L3
AE 1 1 +8(>.+1) +6 (3),+1) +4(7A+1) -6(3A+1)
-K= A2+14A+1 L L2 L L2
-6(A+3) -12(A+1) ... 6(3A+l) +12(A+1)
L2 L3 L2 L3
 - 33 -

C) Condensacin esttica.

11 3
12
2

(1 1 1
1
1)
-@ .,..---- @

Fig. E 2 4g

(1<12)a
(K22)a + (K11) b X

( K21 )b

s, 82

4E 1
-:- -
6Ei
L2 -
2E 1
L
-6EI
7 - 2 E1
L
6E 1
"2
L
6E 1 12EI K12= 6EI -12EI K21= -6E 1 -12E!
-r
L 7 7 ~
L
--"2
L L3

r 4E 1
L'
-6E.I
-----r
L
12 E 1
LlliL2 L3

4 E1
1
-6 E 1
1
l 4 E1
2
6EI
2 4L
2
-6L l
-r L.2 -r- -r E 11
A= -6EI
- L2
1
.12 E 1
t3
1
+
6E1
t

2
12EI
2
= L'3 -6L 12 J +

L2 L3
 - 34 -

6AL
1
E1 ~(1+A) L2 G(A.-1)LJ
12A =~ 6(A-1) L 12(:\+1)

2 2 1
1:!.1 2 L {A + 1 4 A+ 1 } ; A- = - - - L - : : - - - -
2 [
12(:\+1) -6(A.-.l)L J
12EI1 (A +14A.+1) -6(A.-1) L 4(/..+1) L2

L [12(:\+1) -6(t..-1)L ]-! [ 2L 2 6L

12(:\2+14:\+1) -6(A.-1)L 4 (A+ l ) L 2 L - 6L -1 2 J=

12L
2
(5~\-1) 1 44L A.
[
-12L 3 (3A+1) -12L 2 (7A.+1) J

L(SA.-1) 12:\ .,
B A -1 (K ) _: L - 2
-1 21 a L(A2+14A.+1) [-L (3/..+1) -L(7/..+1)J

2
-1 L[1 2 { A.+ 1 ) -6 (A. -1.) L ] [ 2L -6L]
J!.2=A (K12)b=
12(A. 2+14A.+1) - 6 (A.-1)L 4 (A.+ 1 ) L
2
~3 6L -12 =

2
-1 2L (t..- 5) -14 4L ]
[ 2
12L 3 (A.+3) -12L (A.+7)

-12 ]
-L (A+7)

-3LJ L fl(SA.-1)
2
1 2 A. J
-6 X . (A. 2+ 1 4A+ l) l-L (3:\+ 1) -L(7A.+1) =

2(/A+1)L3 2
3(11A+1)L ]
2 E 11
[3(11A+l)L 2 6(13:\+l)L
 - 35 -

-3 1] . A [L (5- A)
2
-12 1
-6 L(A2+14A+1) L (A+3) -L(A+7~=

-4(A+1)L 3 2
3.(A+3)L ]
2A E 1 1 [ 2
-3(3A+1)L 6(A+1)L

3LJ 1 [ L(SA-1) 12A 1

-6 -L2(A2+ 14 A+ 1 )
2
. -L {3A+1) -L(7~+1~=
-4(A+1)L 3 2
-3(3A+1)L ]
=: . 2AEI, [ 2
3(A+3)L 6(A+1)L
L (A +14~+1)
4 2

B =2AEIJ [L 2 3L J A [L (5-A) -12 J =

e 4 = <K21) b -2 3 -3L -6 -L (A+7)
L(A2+14A+1) L2(A+3)

3 2
[ 2(A+7)L -3(A+11)L ]
- 4 2
2
2A E 1 3
L (A +14A+1)
.
-3(A+11)L
2
6(A+13)L

2
2L
(K,, >r= (K11)a -e = ~ 3L]
6 -
2 E IJ

1 L3 [3L L3(A 2 14A+1)

2 4(A+7) 6 (A+3)
[2(7A+l)L 3(11A+1)L1 L
... AE 1 l L
3
' 3(11A+1)L 6(13A+1) 2
A +14A+1
6(A+3) 12 (A+l)
L2 L3
B(A+l) -6(A+3)
L . 2
(K ) -=-e= AEI, L
12 T 2 A2+l 4 A+ 1
6(3A+1)
L2
-12(A+1)
-
L3

8(A+l) 6(3A+1)
-L-
(K21 )r=-C3= AE 1l L2
2
_A +14A+1
-6(A+3) -12(A+l)
L"2 L3
 - 36 --

-3L]
6L

2
2 El 1 [2/.(t.+7)L -3/.(/.+11)LJ
L3(;>..2+l4/.+l) -3/.(/.+ll)L 6;>..(/.+13)

4 (71.+1) -6(3/.+1)
E1 L L2
1
-6(3/.+1) 12(/.+1)
L2 L3
Es decir:
4 (1.+7) 6(/.+3) 8().+1) -6(1.+3)
L L2 L L2
6().+3) 12(/.+1) 6 (3/.+1) -12(/.+1)
L2 L3 L2 L3

8 (/.+1) 4(71.+1)
e L

-12(/..+1) 6(3/.+1) 12(/.+1)

L3 L2 L3
 - 37 -

d) Mtodo de los elementos finitos.

t2 d,= e, -+ w1 (o)

~e
X
'1 )3 d 2 = d y1 -+ w (O)
d = 62 -+ W1 ( L)
3
fw1 x l d4= dy2-+ W(L)
Fig. E2 4h
Se toma L =2 L y se rea 1 iza el ca m
1
X
bio ;=r
1
d1
d2
W={cj>
1 4>2 cj> 3~ cj> ;} d3
d4

Las funciones de forma han de

cumplir:

' cP,. cfl2 ~3 cj>4

;=o o 1 o o
;:;:1 o o o 1
4>1 4>2 4>3 el> J.

;=o 1 o o o
;=1 o o 1 o

Y quedan como sigue:

cP
1
= ;(1-;) 2 .2L
2
cP = i. 1 - 3 ;" + 2; 3
2
cj>3=-; (1-;).2!..
$4= 3;2-2~3

!= ~ 2 { 2 L (- 4 ~ 6 ; ) , - 6+ 1 2 ; , 2 L ( - 2 + 6 ~ ) , 6 - 12 ; }
4
 - ~8 -

1
B="2{ (-4+6t,:) L, -3+6t,:, 1
(-2+6t,:) L, 3-6t,:}= .L 2 !*
- 2L 2

r:i;l (i;) dl;a1 1Jt l;di;+AI 1J}I;di;al 1 x'} +AI 1 x a i-(1+3A)

1 2 1/2 2 1 2 1 7 . 11 .
1 0 t,: 1 ( t,: ) d t,: =! 1 J 0. t,: d ~+A l 1 t,: d t,: = 1 x_ 'i1f +A 1 . IT = i1i" ( 1+ 7 A)
1 1 1 1
2
de donde se obtiene:

E 1 -* ...
K. , = - ! B 1 ( t,:) B: d t,:
IJ 2L 3 o 1 J

K13= ~~lf}(l+A)- 4(1+3A)+~(1+7A)}= ~~ 1 (A+1)

K ={....LJ.. -!2(1+A)+42(1+3A).:.36(1+7A)}= -3EIJ (A+3)
14 2L2 ... 6 22i' 8L2

K2 2 = . 9{
El 9 36 36
2 ( 1+A ) - ~ ( 1+ 3 A) + 21+ { 1+ 7 A) } = .ll.!.J.. (A+ 1 )
2L 4L 3

K23
1
E . 1 { .?_Z ( 1+A) -:-~0 ( 1+ 3 A.) +~
2 (1 + 7 A) } =
6 3 E 1 1 ( 3 A+ 1 )
= 2L 2 "'t BL 2

K24~ !1..J..{..:2.(1+A)+36(1+3A)-36(1+7A)}=

*(
-3EI, (A+1)
2L3 2 8 IT lLJ

K3 3 = ft-1 {~ ( 1+A ) - ~ ( 1+ 3 A) + 1 + 7 .\ ) } = ftl-{ 7 A+ 1 }

 - 39 -

K4 4 = --+{ ~ (
E1

2L ...
1+A ) - 3 ~ ( 1+ 3 A) +~ ( 1+ 7 A) } =---+ (
3E1
4L
;>. + 1 )

Es decir:

A+7 3(A+3) )..+1 -3()..+3)

4T 8L
2 2T 8L
2

3 (;>..+1) 3 (3)..+1) -3()..+1)

4L 3 8L 2 4L 3
K= El 1
7A+1 -3(3"-+1)
s,i-,., --r-
~tr
'~o
8 L2

3(A+1)
4L 3
J
 - 40 -

e) Comparacin de resultados.

Se va a prescindir del factor comn E1 y de los

1
2 3
valores L, L y L , comparando solamente los coeficientes
numricos.

1..=1

2 1, 5 -1 , 5 2 1, 5 -1 , 5
1, 5 1 '5 -1 , 5 1 '5 1 '5 -1 , 5

2 -1 '5 2 -1 '5
. 1' 5 1, 5

K exacta K elementos finitos

Exacta Elementos finitos

 - 41 -

Se observa que_ para A=l la soluci6n exacta y la

obtenida por el mtodo de elementos finitos coinciden. A
medida que A crece los resultados calculados por este se-
gundo mtodo son ms inexactos.
 - 42 -

Ejer-cicio E2.5.
Hallar la matriz de rfgidez de una barra recta de sec
ci6n circular variable, con un radio R en el extremo 1 y r=AR
e]extremo 2, considerada. como un elemento de un em~arrillado
plano. (No se consi.dera alabeo de la secci6n).

-r
R ~
_Ll --~-----
--rc-r
r = :\ R

__. , X

1

Fig. EZ. 5 a
 - 43 -

En una seccin generica x, se tiene:

R-r
RX = r + - x= r{1+
L
k.!.
L
x} con x= X
L

4
= ~~ {1+().-1)x} 4 1
X
= -2 J
X

La energTa de deformacin ~s:

2 ~
u=lLtl_ dx+!L_ dx
2 oE1 o G.J

/Bz,Mz

2J~ 6 ;Mx

dy, Qy .

Fig. E 2 5 b

Teniendo en cuenta la figura E2.5b. se tiene:

M= QY(1-x)+M 2 = (1-x,l)(~y)
z
T= MX = 1M
X

La aplicacin dei teorema de Castigliano, permite

escribir:
edey r =..{ 1 Lo eL- x >-l
El
<1 - x , 1 >d x 1 r~Y}
z
2 1

eX = {/Lo ~}M
GJ X

Las integrales a resolver son:

1
dx 2
. _ 1 { 1 1 } 1= 1 < 1 A +A+1
o{ 1+ (A -1) x} ij - A=T -3 ( 1+ (A -1 ) x) 3 O 3 (A- 1 ) 1 ->:3) = 3 A3
 - 44 -

1 A. (u-11du = 1 {-_!_ +_l_}A.=A.3-3A.+2 =~

2 2
{A.-1) 1
u (A.-1) 2u 2 3~ 3 1 6A. 3 (A.-1) 2 6A. 3
-2 - .
rA.(u-1) 2 du
1
x dx
4 . =
__1 ~{ -.!. + 1
o { 1+ CA. - 1 ) x} 4 = 3 -2
J
1 u (A. -1 ) u u

Por lo tanto se tiene:

2
l =f L ( L-x) 2 ( -) 2
dx =!L L 1-x dx = 4L 3 1 (1-x) - 2 dx
1 o E 1X -::-rf . 4
oE T{
1+ <A- 1 >X"!~-~--E-1T-~ --o { 1+ <A- 1 >x} .

1 4L 3
r 1 =3Tx 4
1T r E

l =fL ( L-x) dx=JL ~(1-x) dx

2 o El 0
X
Ei 4.{ 1+ ( A - 1 ) x} 4

1 = 3.l:!:.!. X 4L2
2 6A.2 ~

1
L d L dx 4L dx
'"1"3=!0 ~ =! 0 - - , . . - - - - - - = f =
El x 4 E'Ir r 4 0{ 1+ ( A- 1 ) x} 4
ET{ 1+ (A. -1 ) x}
.

A2+A.+ 1 4 L
13= X --
3A.3 E7T r
.......... _
L dx A.2+A.+1 4L
14=!0
rr-X = 6A.3
X
:---Ir
G7Tr

l
La. matriz de flexibilidad es:
eX A2+A.+1 1
x-
o ~1
G "'
~

dy
!J. 6A3 1 L3 o
!l!.lx!:.2
=; o n xr
7Tr GA2 E Qy
ez 2A+1 L 2
o --zx- A -FA+1 x- LJ M
z
6A E 3 A3 E

y la matriz de rfgidez es inversa de la anterlor, obtenindose:

 - 45 ;..

4
M 1Tr G 6A3 o ex
X =lf r 2
A +A+1

Qy o E
~3
4A(A 2 +A+1) -- E 2
2A (2A+1) d

j
L2

M
z o -E
~2
2
2A (2A+1)' I
L
4A3 6:
La matriz de equilibrio es:

[~ ~l
o
H= 1
L
~

.Por 1 o tanto 1 a matriz de rigidez de 1 a pieza es:

4 A3
~22=
1Tr
-- -6GL A2+A+1 o o

o ll_A (A 2 +A 1) -lfA 2 (2A+ 1) .

L3 L
2 (2A+1) !!_A3
o -ll.A
L
L2

l. T
1Tr 6G A3 o o =k21
~12=-!i 122= -,- ---c 2
A + A+ 1

o _!..A.( A2+ A+ 1 ) 2E 2
A (2A+1)
L3 L2

o -!f;.. (A+ 2)
.L
2EA2
L.
4
k =H k
-11 - -22
'IJ'r
HT= -,-
-
6G
L 2
A3
A +A+1
o o

2
o !fA ( A + A+ 1.) ll,A (A+2)
L L2 .

o !!.;.. (A+2) !!!.;..

L2 L
 6G A.3 6G A.3
o o o o
L" A.2+A.+l --e- A.2+A.+l f.

'
o 2
!!!_A. {A. +A.+ 1 ) 2EA. {A.+2) o -llA. (A. 2+A.+ 1) ~A. 2 (2A.+l)
l3 7 L3 L

o 2EA (A.+2)
?"
4-A.
L
o -!!A (A.+2)
L2
,..
2EA.2

K= [Kll Kl ~n
6G A3 \ 6G A3
--r- A.2+A.+1
o o :T A2+A.+l o o
K2t Kzz . l ~
~
a-.

o _4EA(A2+A+l) 2E { ) o 4 EA (A 2+A.+ 1 ) ~qA. 2 (2A.+ 1)

L3 . : -LiA. A.+2
i
i? . L

o f ~A 2 (2A+ 1) 2EA.2 o -!!_A 2 ( 2 A+ 1) 4EA3

L2 'L L2 r
 - 47 -

Ejercicio E2.6.
Deducir la matriz de rfgidez de la pieza recta de
seccin constante (E,I,A), plana, adoptando como elementos
isosttico~ las estructuras de la figura E2.6a.

2
1

(11 12)

Fig. E2. 6 a
 - 48 -

Considerando los grados de 1 i.bertad del estado (1)

de la figura E2.6a. las fuerza~ que actGan se representan en
1a f i g u r a E2 6 b

1 2 __

A--_--,,--o;A-, E. 1 J A.
'-./Gz 1 ~Gz 2

l L 1
Fig.E 2 6b

El equilibrio se puede plantear como sigue:

p -1 o o o o
x1 Px2l
P= p
y1 = o
1
rl r1
1

GzlJ
; H= o --1 --L
L
1

p
y2 o . --L --L Gz2 o r
1 1
'L

La energa u de deforl'(lacin (no se considera la de

formacin por cortante)es:

2
U=l.2 1M &_
E1
+ l. ! N2
2
2J:..
EA
G +G
l f(
___ z2 zlx-G )2 dx 1 2 dL
U -E 1 + - p x.2 -EA
2 . L z1 2

au .
aPx2 :=

G _L_ L G
zl 3EI 'bET z2
G +G x3 2
=-. -1 { ::2 Z 1
-G
X } L

El JL2 .zlit' 0

L G .
- Gz1 biT+ z2 mL
au
Es decir aP = F~, con F matriz de flexibilidad, cuya
 - 49 -

expresin es:

L
'EA o o

o -l
3EI
L
-...-
bE 1
.f=
o L L
-biT m
L3
-
siendoiFI=
12 EA (E 1)
2 1 a inversa es:

L2
o o o o
12(E1) 2 -
EA
L
L2 L2 4E 1 2E 1
o 3EAE 1 GEIEA
o -:- -:-
.s- 1=m
1 -1
;F
L2 L2 2 E1 4E 1
o GEIEA 3EAEJ.
o -r -r

--EA L
o o +g
"-
o o

o 6 E1 6 E1 o 12EI 12EI
+ --
-HF~l.., 7 7 -HF- 1 HT= L3 L3

12EI +12EI
6E 1 6EI o
o -7 -~
L
7 7J
Con esto, se obtiene la matriz ~' con el orden de
l~s grados de 1 ibertad (elementos de los vectores ~y !) que
se indican:
 - so -

p
xl
+EA
:- o o
--
EA
L
o o d
x1

+llll 12E.l 6E 1 6E 1
p o o d
yl L3 7 7 7 yl

llll +lill 6 E1 6 E1
py2 o o -7 -7 d
y2
L3 L3

px2
EA
-r o o -EA
L
o o d
x2

~ -- -4E
o 6EI 6E 1 o 1 2EI 6 z1
G
zl L2 L. -:-

zr -- -4EIL
G
z2 o 6E 1 6EI o 2E.I 6 z2
L2 -:-
 - 51 -

Apartado b)
Co n s i de r a n do e 1 e s t a do ( 2 ) d e ' a f i g u r a E2 6 a s e
tienen las fuerzas de la figura E2.6c:

S* =-HS

Fig. E 2. 6c

El equilibrio es:

o o o o
=- o o H= o o
o L o L

La energa de deformacin resulta:

lL. 2d lL 2 d
U= w-f
L o
(G '>+ P
z... y2
x) ...!.
E1
+
2f o
P
x2
x
EA

~ =!L =
dx L
)Px2 o Pxz'EA' EA' Px2

au ... L dx
- 1
... IT{ X
2 x3 L L3 L2
a'P'y2 o
(Gz 2 +P Y 2 x)x El Gz2 -2 +P y 2 J }o=m p y2+2ET Gz2

L
'EA o o
L3 L2
o m 211
Ls l
f= COI)Ifl= rr . 12(E1)
2

o TIT
L2
- El
L
 - 52 -

1:.4
2
o p
12(EI)

o
L2
EAE 1 -- L3
EA El

L4
o - L3
2EAE 1 EA3EI

es decir:
EA
r o o

!Szz=
o
-
.12 E 1
L3
--
6E 1
L2 - E.
-1

o --
6E 1
Lz
4E 1
-:-

HF- 1 = o
o
1
o
o
-1 .
F =
-
EA
L
12EI
o o
6E 1
o -7
o L ~
o 6 E1
-r
L
- ZE
L
1

EA
o o
:-
1 o o -
EA
L
o o
HF- 1 HT= o -
12 E 1
L
'3 -- 6 E1
Lz o 1 L. = o 1 2 E'l
--r
L
6E 1
7
o o 1
o 6EI
--z-
L
-- 2E 1
L o 6EI
~
4EI
-r
 - 53 -

La matriz de rigidez tiene la siguiente expresin:

EA
:- o o --
EA
L
o o

o 12 E 1 6 E1 o 12EI 6 E1
-~
~ 7 7
o 6E 1
7 -4EI
.L o -- 6 E1
L2
2 E1
-r-
K=
-
--EAL . o o EA
'L o o

o --
12 E
L3
1 6EI
-7 o 12 E 1
~
--
6E
L2
1

'

o 6EI
~ -2 E1
L o --
6E 1
L2
'4 E1
-r-
J
donde K es 1a relativa a 1 a
. ...
ecuac1on siguiente:

p
x1 ~ dxl
p
yl dyl
G
z1 =K ezl
p d x2
x2
py2 d y2
G
z2 9 z2
 - 54 -

Ejercicio E2.7.
Se conocen las caracteristicas siguientes de una
barra 1-2 de la figura E2.7a. de un entramado plano.

2 5

3)F..,..._1________Gi)_ 4
t 2

~
......
. - - - - l. -----.4
Fig. E2. 7a

Rigideces: R y R2
1
Area constante: A.
Coeficientes de transmisin:Y 12 , Y21
Se adopta como estructura b&sica la figura E2.7a~ corres
pon d i e nt e a 1o s g r a do s d e 1 i be r t a d ( 3 , 6 , 4 )
Se pide:
1.- Matriz de rigidez elemental ...... K
2.- Matrices de equi 1 ibrio y cinem&ticas !1 1 y !:f 2
3) Si se supone referida la barra a unos ejes generales
que forman un &ngulo a, se pide la matriz de rigidez
K1 (6 x 6). En estos nuevos ejes generales puede poner.
se en la forma: K'= ~T ~ ~. Se pide la expresin ex-
plcita de-
4) Si la pieza est& articulada en un extremo o en los dos
extremos, deducir las exprasiones de ! correspondien-
tes, adoptando siempre la expresin primitiva de ~
 - 55 -

1) Se considera la estructura bsica de la figura E2.7b.

Fig. E 2.7b . - Estructura bsica

La matriz de rigidez es la del sistema siguiente:

pl
p

p
p
xl

yl H K HT
--- -H K
--
d xl

d yl
d y2
d1
-
l
y2
= = = =K
Px2
-K HT K dx2
p2 Gzl d2
8 zi
Gz2
8 zz
... 1
-
donde K = F

Ha.c i en do. d 1=0: Pz=JS..dz es decir: Px2l a1 b1 e1 dx2

Gzlj = a2 bz- c2 ezl

Gz2 a3 b3 c3 8 z2

Para ha 11 a r los coeficientes se dan valores par-

ticulares a d 2 .
 1

!!,2=1 o1+
o
~
Fig. E2.7c
*~ ~
+~
1 = e 1= a 2
b1 a 3:--~

o o o
~
o1 r ~ 1 1 a1

~2= l 1 1 .,A '

"../
1

1
R1 = o b2 c2
11 1 1 b2=R1
-
o1 ... o b3 c 11 O 1 1 b 3=y 12 R1
1 Y12R1 3
Fig. E2.ld
1

U1
O'

'Yic=Y., R, o o
o1
!!.2 = 1 o 1
)
o
Y12R1
=.
al
o b2 cz l: 1
1 e 2 =y 12 Rf

1 ~2 1 c3=R2
1
R2 o b3 e
3
11 1 1
F1g E 2. 7e
 - 57 -

La Matriz elemental K es:

EA o o
L

K= F - 1=
o R1 Y12R1

o Y12R1 R2

2) Se puede expresar el equilibrio siguiente:

p
xl

pyl =
=- p
G
x2
~+
l
Gz2
-rt~
Px1

py1 o
o

-1/L -1/L
o
f:~zl
x2

l
L =

p y2
Gz1
=---
L G~2 J py2 o r
1
r
1
Gz2

Entonces las matrices de equilibrio y cinemtica H1 y

o o o o
H =H=
-1 -
o -r1 -- 1
L
H =HT=
-2 -
O -1/L -1/L
1 1 O -1/L 1/L
o r r
o EA o o
-1 o :-

.
-HF
-1
= o 1/L. 1/L. o =

o -1/L -1/L o Y12R1 R2

EA o o
-:-
R1 ( 1+y 1 2) R2 ( 1+y 21)
o L . L
-R 1 (1 +y 1 2) -R2(l+y21)
o L L
 -. 5.8
'"'
-EA
L
o o o o
-R1 (-t+y12) - R 2 ( 1+-y 2 1 )

-- -
HF- 1 HT= o L L
o -- 1
L
1
r =

R1 R2
o :-( l+y 12) r< 1+y21) o --Ll 1
r

-EA
L
o o
Rl (1+2y12)+R2 . _R1 (1+2y12)+R2
= o
L2 L2
R1 (1+2y12)+R2 R 1 ( 1+ 2y 1 2 ) + R2
o
L2 L2

-F-lHT=
-
-EA
L

o -R
o o

-y12R1
l o -1/L
o o

1/L =
1

o -y12R1 -R
2
o -1/L 1/L

-
-EA
L
o o
R 1 ( 1+y 1 2) - R 1 ( 1+y 1 2)
= o L L

R 2 ( 1+y 21) -R2( 1+y21)

o L L
 - 59 - .

y 1 a matriz de rgidez K es:

l
L
EA o o --
EA
L
o o

Rl (1+2y12)+Rz R (1+2y )+R 2 R 1 ( l +y 1 2) R 2 ( 1+y 21)

1 12
o o
L 2. L2 L L

Rl (1+2y12)+Rz R1 (1+2y12)+Rz -Rl (1+y12) -R2( 1 +y21)_

o o L L
L2. L2

K=
-EA o o EA o o
'L :-
R 1 ( 1+y 1 2) -R1 (1+y12)
e L L
o R1 Y12R1

R 2 ( l +y 21) - R 2 ( 1 +y 21)
o L L
o Y12R1 R2

que corr.esponde a 1 a ecuacin:

p d
xl xl
p d
yl yl

py2 d
y2

Px2 d
= K x2
G 8
zl z1

Gz2 8
z2
 - 60 -

3) Suponiend o el cambio de ejes definido en la figura

E2.7f. se tiene:
P'x2=P x2 cosa.-P y2 sena.

G1 = Gz2
z2

P'xl =P xl cosa.-P yl sena.

Fig. E 2. 7f -'Cambio de ejes

dyl =-d'xl sena.+d'yl cosa.

pJ
xl cosa. -sena. o o o o p
xl
p1 sena. cosa. o o o o p
yl
yl
P'y2 = o o cosa. sena. o o p
y2
P'
x2 o o -senel cosa. o o Px2
'
Gzl o o o o o G
zl
G'z2 o o o o o G
z2
 - 61 -

d xl
coscx sencx o o o o d'
d x1
yl -sencx coscx o o o o d1
d y2 yl
o o coscx -sencx o o d'
y2
d
x2 = o o sencx coscx o o d'
ezl x2
o o o o o ez1
ez2 o o o o o ez2
coscx -sena o o o o
sena cosa o o o o
o o cosa sena o o
BT==
o o -sena cosa o o
o o o o o
o o o o o

4) Las nuevas fuerzas son funcin de las de la matriz

fila ]t= (Px 1 Pyl Py 2 Px 2 Gzl Gz 2 ). Para el caso articulado
e n e 1 ex t re mo 1 ( Gz 1== O) y e o n s i de r a n do . e j e s d e b a r r a s e t i e
ne (ver figura E2.7g.):

_ G,z 1 Gz 2 -Gz 1 -Y 1 ~ 1 G'\

Px1 -E ~::t------;G;.~w------tU:~r-r------t~:/ Pxz

Py21 1 t 1 Py
-Gz 1(l+ Y12 l
L

Fig. E2. 7 g
 - 62 -

p = p pJ o o o o o
x1 x1 xl P)d
( l+y 12) l+y12
P'yl = p 1- L Gz1 p
y1 o o o o p
y ( 1+y, 2) 1+yL12 yl
p =
y2 py2+ L Gz1 pt
y2
o o o L
o p
y2
p =
x2 px2 p
x2
= o o o o O. Px2
1

G'z1 = o G'z1 o o o o o o G
zl
G~z=-y12Gz1+Gz2 G'z2 o o o o -y 12 Gz2

o o o o o
l+y12
o o o L
o
l+y 12
o o o L
o
BT=
o o o o o
o o o o o o
o o o o -y12 1
 - 63 -

Para el caso biart_iculado (Gzl== O y Gz 2 = O) se

tiene(ver figura E2.7h~):

Fig. E2. 7h

p = p p o o o o o p
xl
xl x1 xl
1 1
p = p y 1 -LG z 1-LG z2 p o o o -1/L -1/L p
yl
yl yl
1
p = Py2+LG 1 p o o 1 o 1/L 1/L p
y2 z 1+rG z2 y2 y2

p = Px2 p = o o o o o p
x2 x2 x2
G = o G
zl
o o o o o o G
z1
z1
G = o G
z2
o o o o o o G
z2
Z2

o o o o o
o o o -1/L -1/L

o o o 1/L 1/L
BT= o
o o o 1 . o
o o o o o o
o o o o o o
 - 64 -

Eje re i e i o E2. 8.
Determinar la matriz de rigidez de una pieza recta
de secci6n constante, considerando la deformaci6n por cortan
te, mediante:
a) Teorema de Castigl iano. Soluci6n exacta.
b) M~todo_de los elementos finitos.

- e1-'-------..;;:.~ 1) --
Fis. E2 .. 8 a

A, 1, A' son el rea, inercia y rea reducida de cor--

tante.
E y G son los mGdulos de elasticidad y cortante.
 - 65 -

a) Se adopta el sistema isosttico bsico que se indica

en la figura E2.8b.

Mz 2

..1-x_ _ _ _ _ _ _
ot- 2~~ - Px,
,
- L -----~1 Y2

Fig.E2.8 b

Sea S= Px matriz: en esfuerzos qu~. actan en el nudo 2.

2

'
Los correspondientes movimientos en ese nud.o (defor
maciones) son:
-
e=~
as
dx ds
en donde U= .!..L M2 rr+ .!..L N2 2.+ .!..rL Q2
2 o 2 o En 2 o bA

con M= Mz + p y2 . (L-x)
N= p
x2
T
Q= p y e= (ex2' ey2' ee2)
y2

Entonces, llevando a cabo las integrales anteriores,

se obtiene:

e
y2
=

2
e = aU
62 aH
L
~ {M +P ( L-x) }=M 2
L
IT +Py 2
L
'ITT
z 2 o z 2 y 2 z
 - 66 -

es decir
L
e= F. s con F= EA' o o

o L3 L L2
jEl + bA' 2 E 1.

o L2 L
TIT TI
-1
S= F e= K!_.
EA o
con
'L o
-1
- '1 o L L2 2 L 2
K= F = Gt\1"" + TI'ET GAT + biT

o -.1
L
L2 I!.+
2 L
w+ m 4
GA 1 + 3 E 1
L2
-
Para cal e u 1a r H se tiene en cuenta la figura E2.8b.

: o o
p p p
x1 x2 x2
*
S = p
y1
=-11 p
y2
=- o p
y2.

Mz 1 M L M
z2 Lo z2

o o o o
o 1 o HT= o
!:!=
o L
- o o
1 L
1

-- EA
L
o o

o -1 +1
L3 2
~L 2 L
-HF- ~
1
GA 1 + TiET GA 1 +bET

- 1 _!J..+ L
o L
2
GAT +
L2
biT
4
GA 1 -
L2
+ 3' E 1
 - 67 -

o o

o -1 -1
L3 L2
-F-1HT =-
..L.+
GA 1 -
12 E 1 w+ m
2

_!J.. + L
o 2 L2 L 4 L2
GA 1 + bTi w+ m
J
EA o o
:-

1
o 2
L L3 2 L
lLE- 1 J:!:r= G"'A'i + i2'IT GA 1 + biT

L.
o !..!..+
L
2 L2 4 L2
'G"AT + rrr - -+
GA' m
Luego, la matriz de rfgidez es:

EA o o EA o o
L -r-

o L L3
-z- ---~----. -Lt- o
L
-1
L3 2
1
L2.
GA 1 + T'iET GA 1 + bET GA 1 + T2"ET GA 1 + bET

o 1
El + L -1 E1 L
o -- + -
2 L2 'L 4 L ; 2 L2 L 4 L2
GA 1 + m GA 1 + 3IT GA 1 + biT GA 1 + )TI ~
():)
11
K= 11
EA EA
--
L
o o L o o

-1 -1 -1
o 2
o 2
L L:f 2 L L L) 2 L
GA 1 + 1TET GA 1 + bET GA 1 + 1TET GA 1 + bET

-1 El L
o 2
-E-
1 + L -- o - + --
2 L L 4 L2 2 L2: L . 4 L2
GA 1 + bET GAT +m 'G'A'I + biT 'ff'AT'+ 3EI.
 - 69 -

b) Mtodo de los .elementos finitos.

Conviene redefinir ei vector deformacim y las varia
bles de campo, ya que ahora no es Qnicamente la flecha w, si-
no el ngulo de giro de la seccin, e.
SegQn la figura
E2.8c. se tiene:
dw
dX +a=e.
con w flecha
o. giro debido al
cortante.
e ngulo de giro
total de la seccin.
b) Movimiento As, pues, en una sec
sin deformacin
cin genrica de una
viga no son iguales e
dw .d
y dx'. s 1 se .cons: era
la deformacin por --
e l Deformo.cton cortante. El vector -
~~~...::.:;;:.-.. de. fiex in
de campo es pues
, e) r.
u=\w,
Las ecuaciones -
constitutivas son
1 de
d l Deformacin M= IT dx
de co1tanw
e= GA'1 o..
La energa
de deformacin de una
barra es entonces:

el Deformacin de
une rebanada

'
\~w
-.ai:_
, Cl /
de viga recta
{E2.8a.}
..... \ . ,..,.1 ..... ..... - Se adoptan como
funciones de forma pa
ra un elemento las 1 i
-
nea 1es, ya que e 1
-
pro

Fig. E2. ae
 -70 -

Los grados de 1 ibertad del elemento de flexin se

ven en la figura E2.8e.

wl Wz

e,( t })e,
f-~~
Fig. E2.8 d.- Elemento de
flexin

Siendo
X
N1 = 1--
L
y

El vector deformacin es:

d8
e:= [ dx J
CL .

y el correspondiente vector tensin cr es:

M.,
g_= [ 8 J
La relacin deformacin-desplazamientos, se obtiene
a partir de la expresin:
d L
=
-
e: = 8
dx w
-u
d
-crx e
Por 1o tanto
B d
-e: = {!1
 - 71 -

d.={w.,
con -1 1
e.} vector de de.splazamientos (inc5gnitas b~sicas)
. 1
en el nudo l .
dN.

-
1
B.= o CiX
d.N.
1
-crx N.
1

La matriz el~stica del material ~' es

D= diag.{EI, GA'}

Por lo tanto la matriz de rigidez, es:

k= f L BT D B d x.
-- o - ........
No obstante, convien.~, por razones de estabi 1 idad
nGmerica, separar las conttibucienes de las deformaciones por
cortante y flexi5n en la matriz de rigidez. Por lo que si se
~ust[tuyen las expresiones {~2.8b.} en {E2.8a.}, se obtiene
(se supone secci5n constante):

U= fL ~T e +

GA'
-r{~'
T er}
- {/~ N2
1,X
N
11X N2 ,x -Nl ,xN1 -N 1 , X N2
dx} [:
N N2,x -N 2,x N2
2,x -N2,xN1
-N . N1 -N 2,x N1 N2 N 1 N2
1t X 1
-N l,x N2 -N 2,x N2 N2
N1N2 2
aN.1
con N.
1 1 X
= a"X""""t 9=

[::] -
.y w=

[::] .,
como L N . 2 1 L N.
o 1, X
dx= r o
1, X
N.
J 'X
dx=--
L
i# j

/L N 1 L . 1
N. dx= N N d x=--
o 1, X 11 X 2 o 2,x 1 2
L N~ dx= -1 L N N dx= .b.b
o 1 3 o 1 J
 - 72 -

se deduce

GA 1 {w T , aT} -1 -.L/2
-~] a
l2 lL-aT [ - Lf2
- [~]
U= l + n
-1
-1 L/2 L/2
-~ 2
-- -L2
'
L L
2 2 J b
L2 L2
--2L -L
2 b 3

La matriz de rgidez es; evidentemente:

k2 -k
2
-k
2
L
2
L
-k 2 2L
L
1 L -
con k = E 1 y GA 1
k2~ 'L
-k 2 -k 2 k2 2 k2 2
k= L L2 L2
-k 2 'i"'L k2 2 k 1+k2 -k1+k2
32 b2
.;.k L L L L
2 2 k2 2 -k,+k2 b kl +k2 r
El problema fundamental con la consideraci5n del
,cortante en la matriz de rigidez, es la posible inestabil i-
dad numrica que aparece. Si k es el canto de la viga, se
h'
denomina factor de aspecto al valor t Las caractersticas
~. - 1 3-
mecan1cas de la viga son 1= T2 bh y A1 = k b h 1 (con k de-
pendiendo del tipo de secci6n) siendo b el ancho de la sec
ci5n' (se supone, .la secci5n re~tangular. per()es_aplicable.
a cualquier secci6n general).
L2
Para el valor k 1+k 2 ! =
b ~: L
2
{ k + ~: rt }:,
h
Se comprende que si el factor de aspecto L' es -
pequeo, en el elemento (caso normal de vigas delgadas), la
contribuci6n de la flexi5n es despreciable, es decir ~<1
y con mayor motivo tr<.~) 2 <.<k, y no puede ser recogida en e 1
c~lculo con un computador que trabaja con un nOmero finito

de cifras significativas.
Existe la tcnica de la integraci6n reducida, con
objeto de obviar el problema anterior. Se utiliza, un solo
punto de Gauss, para evaluar las integrales de contribuci5n
 - 73 -

h
del cortante, en los casos de vigas delgadas (t<<1). En otras
situaciones, se calcula exactamente.
 - 74 -

Ejercicio E2.9.
Determinar la'~atriz de rfgidez de una viga de ton
gttud infinita (viga semi-infinita) y secci6n constante Figu
. -
ra E2.9., considerada como un elemento de emparrillado plano.
Se supone que esta viga se encuentra inmersa en un medio
el&stico, de m6dulo de balasto a desplazamientQs verticales
C.k w} y a g i ro s d e t o r s i 6 n (k 6 ) La s e a r a e t e r s t i e a s e 1 a s t o - -
mecnicas de la barra son: El, GJ.
No se considera la deformaci6n por cortante ni el
alabeo ae la secci6n.

Fig. E2.9a. - V1ga infinita. Grals de lbertaa

 - 75 -

...
La ecuac1on de 1a flexin es:
4'
d w kw
El 4 ' +k v1= o S i e4=
liET se tiene
dx w
4
w IV + 46 w= o para X (o' oo)

Las condiciones de borde son:

= Qo } para x= O

valores finitos de w para x= oo

La solucin del problema.anterior es:

w(x)=(A cos6x+B senSx)e-Sx

Por lo tanto
=
[:~
wo
J _; o
1
J [A]
B

[:a
w'' o
r [2
o .: l [A
B
J
1 1 'W ,,_ 1
'Wb= ' 'Wo '=
1 . 1 1 1
con . -62w 'o' y o -
~o wo
S!
Eliminando
I: J se tiene:

[=~
-1

[~~:o ']
= 2 2 = 4
L, ~ l JG~lo
1
[o . -2 J wo J [_2

resulta por lo tanto:

= 2EI

La ecuacin diferencial de la tersin es:

S j
k8
GJ-
 - 76 -

es decir

con las condiciones de bord~

y GJ 2,i
dx
= T
O para x=O

valores finitos de e para x=oo

La solucin es:
e= A COSAX.

Por lo tanto
A=e y GJ). A= T
o 0
.Eliminando A, se obtiene la rgidez

0
T ~ GJ). e0
La matriz de rgidez pedida, es por consiguiente~

2 2
ey = 4E113 2EII3 o d
y
2
Gz 2EI13 2 E 113 D e
z
G
X
o o GJ). ex

4 kw 2 ke
con ~ = '4Ti y A =GJ
-

NOTA.

La matriz propuesta puede servir para modelar losas

muy largas mediante el m~todo del emparrillado. En esa d~rec

cin se puede aproximar los valores de B y A como sigue, ob--

tenidos al determinar el coeficiente de balasto,o relacin
carga/flecha, que al mov.imiento de la viga opone las vigas
transversales del emparrillado:

con L luz transversal de la losa

b ancho de la viga infinita
 - 77 -

f factor de la flecha producida por una carga unidad puntual

aplicada en el arranque de la viga, es decir

( P= 1 ) L 3
f= f El
Por ejemplo f= ~simplemente apoyada con carga en
el centro de la luz.
eo mo
t a n t e o s e p u e d e s u p o n e r L a-1 <;'
-1 L
t< . .,
L( e n 1 o s a s 'b"1a -

poyadas o empotradas) y S ~!(en voladizos), para el valor

mximo de S. Ver figura E2.9b.

empotramiento
/ / / /,/,// / / / / /

apoyos o empotramientos
~'""
_ __.___a_ -
~"'max
''-........ libre

F ig. E 2. 9 b.- Jnterp~lacin dei mdulo de balasto

Figura E2.9b. Interpolacin del mdulo de balasto

2) ,..-~--

A= v'( 1+~)6bL
siendo. b y L cantidades anteriormente definidas y el coe- e
ficiente adimensibnal del giro a producido por un momento -
un1 d a d ap 1 1ca
d o a 1 a v1ga,
es d ec1r
a= -a (M=El1 ) L Por ejem-
- 1 '
p1o 8=-- en una viga simplemente apoyada con carga en el
4-2
centro d la luz.
 - 78 -

Ejercicio E2.10.
Resolver idntico problema al 2.9., con~iderando la
existencia de un m6dulo de balasto al giro ~ongitudinal k~.
1

Fig. 2. 10 a.- Viga infinita. Gradas de libertad

 - 79 -

Se supone por lo tanto la actuacin de las siguie~

tes cargas repartidas (se tiene en cuenta solo la flexin, -
ya que la torsin est desacoplada y se resuelve igual que en
E2 9. ) :
Fuerza vertical -k w
w
dw
Momento -k,,,y=- k,,,-
'~" dx
'~"

Equ i 1 i br i o de una rebanada (ver figura E2.10b.):

-jdxr-
Fig. E2. 10b

d Q.
dx + p= 0 es decir

+.2!i + g+Q= o es decir dm + Q k... dw 0

dx -dx -.~<J Tx =

Relacin constitutiva es una rebanada

2
El~=
2
M
dx
Por lo tanto, se obtiene la ecuacin diferencial
general:
2 2
d
(El ~) -~(k ~) +k w= O
dx2 dx2 dx ~ dx w

o bien en este caso

4 d2w
El ~-k --- +k w= O en x (0, 1).
dx'+ 1/J dx 2 w
 - 80 -

Con las condiciones cinem~icas (esenciales):

dw . ,
w=w o dx = 'Po

y las condiciones estticas (naturales):

para la seccin origen x=O. Los valores anteriores deben de

ser finitos para x==.
Para la resolucin de la ecuacin diferencial, con
viene distinguir tres casos:
1) Debil amortiguamiento del balasto k~, es decir:
2
k - 4EJ k <O.
VI w
2) Fuerte amortiguamiento del balasto kl/i, es decir:
k2 - 4E 1 k >O.
tP w .
3) Valor crtico del mdu 1o de balasto kt/J, es decir:

klP
2
k =O.
w - 4 Et
Se estudia a continuacin cada uno de los casos
anteriores:
1) Caso: k~"- 4EI k <O
"' w
k k
Se denomina s 4= w y cosa.= P._

La ecuacin diferencial resulta:

4
d w 2
2
d w 4
~ -28 cosa.--r +8 = O
dx dx
. La ecuacin caracterfstica tiene como rafees:

t= S(cos~ isen~)

De stas, por la condicin de solucin finita en x=O,

solo se considera la correspondiente a la parte real nega-
tiva, es decir, la solucin es:
-rx
w= e (c cos sx +
1
c 2 sen sx)
 - 81 -

a a
con r= Bcos'i' s= ....>asen-
2

Las derivadas de w son (ver nota al final del

ejercicio):

e< 1 >=-(e cos a - e2 sen!) B a

e (1)=- (e cos2 + e2 sen~) S
1 1 2 2 1

e (2)
1 = (e 1 cosa - e2 sena) S
2 e(2)= (e 1 cosa +
2 e2 sena)B
2

"a
c<1 3 >=-(e 1 cos-' 2 - c2 sen~)B 3 c<z),..=(e e os~+ c2 senza
3 1
3 ) (3 3

w(n)= e-rx{ein) cos sx + e ( n) sen sx} n= o, 1, 2., 3

2
Las condiciones en el borde x=O, son:

\'1 o= e 1
a
sen-
2

1
;z Mo El{-e 1 cosa + e2 sena}

1
;3 Q =-EI{ e, cosa- e2 senCX} + 2E 1 cosa{e, e os.-
0 2
a
- e2 se ni}

es decir, en f o r ma ma t r i e i a 1 :
d =

~:: J
=

[-c:s~ s:n~J
.
[e,2
J c2
= G e
-d

e. =
1
-3QO
B
El cos~- -sen~] e, = G
-p
e
1 -cosa sena e2
-32 Mo

Eliminando e se obtiene:

-1 -d
p= Gp Gd
 - 82 -

siendo
a. a.
[ cosr J[ sen2a. o
2cos~]
sen 2
-cosa. senot cos~
2

La matriz de rigidez es k:

..e.= = El [za 3 cos 2a. a2 = k d

J J
[ Q
Mo a2 2Scos~] [ wo
VJo
 - 83 -

2) Caso k:-4EI
'1'
k w>O

Si se denominaS=
4 kw
ET y eh~=
kljJ
, la ecua
2
2S El
ci6n "diferencial se convierte en:

Las races de la ecuaci6n caracterstica son:

Por la condici6n de valor finito en x=oo, solo se co~

sideran las dos races negativas. La soluci6n general es en--
tonces:

w= e 1 e t 1 BX .
+
C
2
t Bx
e 2

con t 1 =(Ch~+Sh~) 112 y t

2
= (Ch~-Sh~) 112
La derivada n-sima, es:

Las condiciones de borde son:

o en forma matricial

[l: 1
w

S o
=
[_ t
1

1
-t
1

2
]
e =
[e: 1
G

-o
e

2 2
.e.= .!.!lo =E 1 t1t2 t1t2 e~ = G e
s3 -p
1 -t2 -t2
-2M O 1 2 c2
S
 - 84 -

Eliminando C se deduce

-1 -d
.e.= G
-p ~

con

La matriz de rigidez es, considerando que t +t = 2Ch~:

1 2

.e.=
 - 85 -

2
3) Caso crtico k -4EI k =O
'lJ! \>J

La ecuacin diferencial es:

siendo S~ kw y
IT

Las r a e es (do b 1es ) de 1 a e e u a e i n e a r a'c ter s t i e a son :

t= tf3
Solo se considera t=-S, por la condicin de valores fi
nitos en x=~ la solucin general es:
w= ( c sx+C ) e -Sx
1 2
La derivada n-sima de esta ecuacin se obtiene, segn
se indica en la nota final de este ejercicio:

con
e< n ) = (- S ) n e
1 1

c(n)= n(-S)n- 1 s~ 1+(-S)nc 2 =(-S)nc 2 =(-S)n{-nc 1+C 2 }

2

Es decir

c<o)= e
1 1

C(1)=-SC
1 1 ,
c(2)= o2c
1 ...> 1 , .

c<1 3 )=-s
.
3c
1 ,

Las condicones de borde, producen los resultados:

 - 86 -

1
-M2 o= E1{2e 1 -e 2 }
6
1
- 3 Qo= El{3e 1 -e 2 }-2EI{e 1 -e 2 }: El{e 1 +e 2 }
6
En forma matricial

d=

!'a
1
wo =
[: _:] [::l =
~ e

.e.= -3 Qo =E 1 e1 = G
"""!p
e
S
1
-2M O 2 -1 e2
e
Por lo tanto
-1
G ~
a= -p d

con
-1
G = El [~
-p Gd.
[~ - ~ J [ ~ ~ J= El
:J
La matriz de rigidez es

.e. =
[Qo J = El [2:: 2S
ez
J
( wo 1= k d

Mo $0
 - 87 -

Notas al ejercicio:
a) Derivada n-sima de la funcin:

w= eax { y ces bx+ e sen b x}

2

luego la derivada n-sima es:

(a+bi)x
e = e p n (e os nw+ i sen n w) (e os bx+ i sen b x) e a?<
1

luego .
w(n= eax(e(n)cos bx+e(n)sen bx)
1 2
con '
( n) 1 n
e 1 =(e cos nw+e sen nw)p
1 2

. e~n)=(-c 1 sen nW+C cos nw)pn

2
1
p=(a 2 +b 2 )"2
b
W=arctg-
a

b) Derivada n-sima de la funcin

y=C xeax
1

Se considera y= a! {c 1eax}. La derivada n-sima es

n .
a {a~ (e e ax 1}= -a ( e a n e ax ) =e { (na n-1 +a n x ,L e ax}
y( n ) =........_
ax" aa 1 aa 1 1
 - 88 -

Ejercicio EZ.ll.
Resolver el ejercicio 2.9 mediante el m~todo de los
elementos finitos. Adoptar como funciones de forma funciones
del tipo (ax+b )e-A.x.

. 3

Fig. 1:2 11 a - Viga mfimta: Grados de liberrad

 - 89 -

Energa potencial total en flexn es:

2
d w 2
1
V(w)= 2 r El (--z) dx+-
o
00

dx 2 r o kww
1 00 2
dx+(w 0 ,w 0 ) (Qo)
M
o
-J..x
Por otra parte w=(a x+b)e , debe satisfacer
dw
x=O w=w 0 CiX =w0 , es decir
b=w 0 a -J..b= 1J1 0

Por lo tanto w={ (A x+ 1 ) e -Ax , x e -J..x} = N d

La relacin deformacin-desplazamiento es:

e:= B d
2
d w
con e:= 2
dx

La matriz de rTgidez es de dimensin (2 x 2) y su

ex p r e s i n e s :

Las integrales se ~alculan elemento a elemento, Pa-

ra ello, se realizan los clculos previos siguientes:

rooDe ':'2Axd X"'fi1 y por derivacin sucesiva, se deduce:

roooe -2J..xd x=-1

4;.. 2
rooox 2 e -2Axd x=-1
4;.. 3
Entonces:
 - 90 -

kw
Por lo tanto i se denomina 134=
S
TI se tiene:
4 A.2 a4
k= E 1 1 A3 + sa
1i" ""U r+
2A2

A2 134 5 A + 134
r+
-;z 4 4A.3

El valor de A puede ser arbitrario. Para cada uno

de ellos se obtiene una matriz de rg~dez m~s o menos aprox!
mada. Si se adopta A. de modo que k ten~a mnima la traza se
obtiene, A.=S.
La matriz de rigidez, m~s adecuada, es entonces:

k= E 1 3 (33
2

2
(3 313
2
Si se adopta la misma definici5n de que en el a
ejercicio 2. 9. , es decir S= t~l , 1a matriz de rigidez an-
terior se convierte en la siguiente:

~Yz
k= El 133 2 3 2
= El [4,2436
26 ] 26 ]
2S
2 3 t;"s 213 2 2,121(3

que se acerca a la soluci5n exacta k= El 2

[4 6 3 26 ]
213 2 213
 - 91 -

Ejercicio E2.12.
Determinar las 1eacciones de empotramiento rgido
para una carga horizontal uniformemente repartida de valor
p, del elemento representado en la figura E2.12a. suponien-
do forma parte de una celosa.
Se utilizar:
a) Mtodo general energtico. Teorema segundo de Cas--
tigliano.
b) Ecuacin diferencial.
e) Elementos finitos.

kA

~- --------._A-
-
X
L

Fig.E2. 12a
 - 92 -

a) Mtodo general,
La matriz de rgidez exacta, calculada en el ejer
cicio 2.1 es:
k= EA (k -1 )
- L 1n k

El desplazamiento en la estructura b&sica se deteL

mina segn el teorema de Castigliano. (Figura E2.12b.).
p
~----~
Fig. E 2. 1b

con N 0 = p(L-x)
N = 1
1
A = A (k+
1
f x) .Y S
X
t=~, se puede poner:

e = f..b L ( 1 - t) d n
o EA o k+(l-k)t
Las reacciones de empotramiento rTgido son:

[::] -
es decir:

"' l+ EA(k-1)
p,=-p L ln k

"
p
2
=- EA (k -1 )
L ln k.

Se comprueba que p1 + p2 = pl.

 - 93 -

b) La ecuacin diferencial es:

_2. ( EA ..!iLY) = en x {o, 1)

dx dx P

n=O para x=O y x=L.

La solucin de la ecuacin, es directa:

EA .Tx
du
px+ e 1 P L
{2.12a}
k
px+e]p P.J:.2
1 e1 - 1-k }
d L.F 1-k - dx= EA {1-k + k+(l-k)~ d~
EA (k+-Lx)

es decir

oL
2
~ k 1
n= !t l-k *(e 1 -r=k)l-k ln{k+(l-k)~}
1

Las condiciones de contorno conducen a:

O= (C
1
- k
1-k).!.!!l
1-k
+e
2

1
O=-+
k-1

... e 1 e 1 k
2 =~ y l =
lnk - - -
L a so 1 uc1on es k-1

La ecu.acin {2.12a.} permite calcular las reaccio

nes de empotramiento:

EA d !JI - { 1nk
1
dx x=O- PL
 - 94 -

e) Mtodo de los elementos finitos.

Segn P.l ejercicio 2.1., se tienen las funciones de
forma:

X
N "" 1--
1 L
y N
2
=-XL

y las reacciones de empotramtento son, simplemente:

1
L.T L
S. 1 =-!~
o 1 p dx=-pLI o (1-~)d~=-p
2

Que son independientes de la variaci5n de la secci5n.

 - 95 -

Ejercicio E2.13.

f- L
1 L

-lr2
\j
2
-='A
1
12

Fig.E 2. 73 a

Hallar las reacciones de empotramiento rgido, en

la viga de la figura E2.13a. para los casos de carga de la
figura E2.13b.

p=lmxt
p=lt/m. m

I 111 11111111 n
~----------------------~2 ------------------~2
a} b}
F1g. E2 13 b

Se utilizarn los procedimientos siguientes:

a) Condensacin esttica.
b) Mtodo de los elementos finitos.
 - 96 -

a) CON DEN SAC 1 ON ESTATICA,

f.., k11 k12 k13 .2., .e.,

"' e
.1 A
3

Fig. E 2.13 e
')
--- e
f..z - k21 k22 k23 d
-2- + "
E.z P=k d +
'
p "'e..3
-3 k31 k32 k33 .2.3

d =d...2. =O--
Pero se cumple que -1

Luego:
A .J.
1'"= "' -1 "'
1 .e..,-k13'k33 .e.3

l
p=1 t /m

L~ e~ ll 11 u 1111 1 ) g
"' . -1 "
~2= E.z- k23 k33 .e..3
--

1~ ~ tl ~ t~
Fig. E2.13.d

Teniendo en cuenta el apartado e) del ejercicio2.4.,

resu.lta .. (Figura E2,13d.).
k13=(kl2)a k23=(k21)b k33=(k22)a +

+~
2 2
2L -3LJ [2AJ., +3ALJ =
[ -3L 6 L 3AL 6A
2
2 E 11 2L (1+A) 3L(A-:-1)]
1
L3 [ 3L(A-l) 6(A+1)
 - 97 -

2
L -3 L -3L(-1)]
r(+l)

3L -6
X
6EI
1
(~ 2 +14+1)-3L(-1) 2L
2
(!c+1)

2 2
6L. (Sic-1) -3L (31c+1)

2
36~1c -3L (7!c+1)

2
6 (lc+1) -3L (A-l) ]
L 3L]
[ -3L -6 [ ~3L(Ic-1) 2L 2 (!c+1)

-36L
 - 98 -

1~ Caso: Carga unidad.

2
... 3 L3 (3+ 1 )

n
6L (5A-1)
1 =
2
26 LA -3L (7A+1)

4
-3 L (3A+1)

2 2 -JJ(7A+l)
3 L ( A + 1 4 A+ 1 )

=
-36 L

.L2 2
2 A +50A+13
-L2 -.L (3A+1) -X
12 12 A 2 +14A+l-
1 P*=
1 2
~r= 2 L
X
A +28A+3
L A +14+l) -L(7A+1) 2 2
...
";) A +1 4-A+ 1

:.L2 2 2
--12
[ t (A+3)] L2
-rr x
13 +50A+1
2
p*= A . P''(=
A +14A+1
2 2 - L(A+7) ' 2
A +14A+l
L 2
-r , L' 3 +28A+ 1
-X
2 2
A +t4+l
 - 99 -

Comprobacin:

Equilibrio fuerzas verticales:

L
-X
2

Equilibrio de momentos:

2
13A +50A+1 + 2 L2-L2 3A~+28).+1 =
A2 +14A+1 A2+14A+l
13A 2 +50A+1+24A 2 +336A+~4-36A 2 336A-12
2
A. +14A+1
 - 100 -

22 CASO: MOMENTO UNIDAD.

p =1m x t/m
1 3 2
1)))) t)) ))_f
Fig. E2. 13e

Para obtener las reacciones de empotramiento rgi-

do debidas a un momento unidad uniformemente repartido se de
rivan las de la carga vertical unidad situada en un punto g~
nrico y se integra"el resultado a lo largo de toda la viga.

M_,
e~
1
X

t ~ )~2 M,( 1
1
X
...
:j ~)~
Qlt l
12 ' of
1 t02

Fig.E2.73f Fig.E2.13g

2
M1=-x(1-T:) 2 4 + -
M =-{1-2. 3x }
1 L L2
2
Q =-{1-~~(3-2~)} ~( 1- ~)
1 L2 L . L2 L
Derivando cada
2 ~reaccin respecto {E2.13a.}
M=~ (1-!.)
2 L L
de x, resulta: M
2
= -x(
L
3x)
2-
L
-

-6x x
Q = -( 1--)
2 L2 L

Integrando los resultados a lo largo de toda

la luz L, se deduce (figura E2.13h.);
 - 1 o1 -

Fig. ~2.13h

~ L 4x 3x 2 2
t1 -=-! {1-- + - )dx=-{L--
1 o L L2 L

~ L 6x 6x 2 6 L2
Q1 = ! 0 (-z- -y )dx=-r{~
L L L
A L 6 6 2
Q =-f X X d 1
2 O. L2 - ~ x=-
2 2
M2= rL(~- 3x )dx=!h
, o L L2 2L

Queda por lo tanto lo indicado

en la figura.E2.13i.
' ~

Fig. E2 13

Se observa que
p =O
3 A..J.-A
L uego P1=P 1
 - 102 -

b) METODO DE LOS ELEMENTOS FINITO S.

Tenien do en cuenta las funcion es de forma definid as

en el ejercic io 2.4. (aparta do d), resulta n las siguie ntes
fuerzas consis tentes en los nudos:
A

- t1
,..l
-Q,..1
-M,..2 .=! L cj>.b
T
dx
o - -
-Q 2

Caso de p=l ! o ~(1-~) 2 .2L d~

1

1
1
b --3L2
(1-3~ +2~ )d~
2 3 1 -L
! 2
p=2L!L !E_T (-1)d~=-2L o
1 =-2L L2
o -! ~ 2 (1-~).2L.: d~ L
o -b 3
/1 (3~2-2~3)d~ L
-b -L
o

!
1
(1-4~+3~ 2 )2L d~
_ro-- o
o
(-6~+6~ 2 )d~
Caso de m=l 1
/
o -1 -2L

o (2~-3~ ) .2L-
1
2
p=2L/
1
cj>'T.1.d~=2L
- o-
/ d~
.. 2L o =
o
1
! (6~-6~ 2 )d~ 2L
o
 - 103 -

Ejercicio E2.14.
Oe t e r 01 i n a r 1 a s ea r g a s e q u i va 1 e n t e s ( s o 1u e i n i n i e i a 1)
de la barra cnica de la figura E2.14a. sometida la torsin
pura.

Fig. 2 14 o
 - 104 -

La ecuacin diferencial es:

con K-=GJ
1
mdulo de torsin

y las condiciones de borde 8=0 para x=O y x= 1

El problema constituye un ejemplo de determinacin

de la funcin de Creen (o lnea de influencia). Se utiliza
un procedimiento tpico de clculo.

G1T{ - x }4 =
+R-r X ' R
K= -r con ~= y 1\=-
T 2 L L r

Integrando la ecuacin diferencial resulta:

K ~ =C
T dx 1

Por lo tanto
2C 1L 1 1
8=- - 4 (X-1) + c2
3G1Tr {1+().-1)~} 3
1
3.E.J.1.4 r:r
Para ~=O 8=- +C =O
6G1rr 2

luego
ZC 1.L
8=
3G1Tr 4 {).-1)
{ 1-
{1+(~-1)~}3} para

Anlogamente adoptando los ejes de la figura E2.14b.,

se tiene:

8= 41 {1- ,, ,-3 }. para ~'= rx'

3G'lTR (--1) {1+<r-,>~}
A
 -.105-

-
- Gx

~ ------------~-
1'-
)(
X~ 2
a --...f---- b

F1g. E2.14 b

En el punto x=a, se cumple el .equi 1 ibrio:

KT ~ .2,! 1 1 b=- Gx
dx 1 a + KTdx

es decir e1 + e1=-G X {E12.14a.}

Adems se debe satisfacer la compatibilidad o bien e= e

a b

. e1 r
1- . 1 . "1 1
4 { 1- 1 -X - } 3} , { E2. 1 4 b. }
r (A.-1) { 1+ ( ~.. - 1 ) a} 3 {1+ -x- b

La solucin de las ecuaciones {E2.14a.} y {E2.14b.}

es:

e=-{1- 1 } { 1 }- 1 G
1 (at..+ 'b) 3 1 --;:_- 3 . x

1..3
e 1= {
1-
(-
3
-)3}{1.. -1}
-1
Gx
aA.+b
 - 106 -

Los momentos de empotramiento rTgido (soluci6n

inicial) son:

3
0.3-1}-1
t1T1= {1- A } G
(al.+ 'b) 3 X

1
MT2= { 1 - 1 } {1--} G
CaJ.+'b) 3
X
1.3

Si J.-+-1, resultan expresiones indeterminadas. Sin embargo,

desarrollando en serie d J.-1 o bien por L 1 H8pita1 se obtie
ne:

y
 - 107 -

Ejercicio E2.15.
Determinar la matriz de rgidez de la barra curv~
de la figura E2.15a. cuya directriz es parablica de segundo
grado, suponiendo que forma parte de una estructura celosa
plana.
No se considera deformacin por cortante. A es el
~ ~. e 1 su inercia a flexin. Se supone
are a total de 1 a secc1on
que se cumple en cualquier seccin de 1a barra: A
- = Ao,
cosa
1 cosa= lo con Ao e constantes.
1
'o

f
h

Fig E 2 75 a
 - 108 -

La ecuaci6n de la directriz es

Se utilizar la tcnica de la matriz de flexib.ilidad para

obtener la rigidez de la barra. La fuerza ~ctuante en la
barra isosttica es S y su deformaci6n eficaz e.

Fig. E2.15 b

La energTa de deformacin de la barra sustentada

tomo se indica en la figura E2.15b., es:

2
U= lL ~ + l L N2 ~ =lL {S y) i_E~I + l L (S cosa) 2
2 o El 2 o EA 2 o 2 o EA d~

Con. L la longitud del arco.

Por el segundo teorema de Castigl iano, se puede

escribir
2
au .2,l + L cos a
e=-=
aH EA
E1 o

con dx=d~ cosa, la ecuacin anterior se transforma en la si-

guiente:

e =s{....!. '- 2 1 L dx}

E 1 o y dx+E'A o
0 0

es decir:
..,
. . 8h... 1
e=L{,_--SEI +EA }S= F.S
o o
 - 109 -

La matriz de equi 1 ibrio es H=(1), con lo que la ma

triz de rgidez resulta:

~
-1
K= [ HkH -Hk J con k=F
-kH k

Fe r 1o t a n t o :

E ( 1
K= -1
 - 11 o -

Ejercicio E2.16.
Determinar la matriz de rigidez del elemento barra
(Figura E2.16a.), sometida a torsin no uniforme, consideran
do la coaccin de alabeo.

GJ mdulo de torsin de St. Venant.

El mdulo de torsin por alabeo.

w

Fig. E 2. J6 a - ~rados de libertad de 1a barra a tors1n

no uniforme

Grados de libertad de la barra a torsin no unifor

me.
 - 1, 1 -

Se sabe que M = GJ ~ 1 -EI 9' 1 '= MST + Mw

X
X W X

con M momento torsional total

X
ST
MX momento torsional de St, Venant = GJ~'

t1w momento torsional de alabeo =-El w~

X .

La ecuacin diferencial general (seccin constan-

te) es:

IV
GJ9' '-El W<P =m
X

La matriz de rigidez se obtiene, de la ecuacin ho

mognea. La solucin general (complementa ria) de la ecuacin:

GJ~''-EI ~IV=O
w
es:
con p=~
Las ecuaciones de borde son:

cinemticas ~, 8= ~
estticas +GJ~'=MST -El = Mw
X' W X

En los extremos x=O y x=L.

d= ~, = o o Al = ~d A
e, p o p o Az
~2 S e pl A3
ez pe ps p o A4

-
E.= MX 1 = -GJ 9 o o = G A
2
All -p
M
w1 p E lw o o o A2

M
x2
cGJ sGJ o o A3
2 2
M
w2 -p e 1 lw -p sE lw o o A4
 - 112 -

siendo s=senpl y c=cospl.

Por 1 o t a n t o .e.= ~ ~ con k= G G-1
-p -d

E 1 resultado final es

k=
[~11
1.21
~12]
k
-22
siendo

~1 1
1
= "F.
r-ps -(c-1) J
l-(c-1) l(s-pLc)
p

k12=~~ 1 =
1 [ ps ~(c-1)]
A c-1 -( p- S)
p

1.22= A
1 [-ps . c-1
1 J
c-1 -(S-pLc)
p

A GJ
-2(c-i)-pLs

-ps -(c-1) ps -(c-1)

1 1
-(p-s)
-(c-1) -(s-pLc) e -1
1 p p
K=
A e -1 -ps c-1
ps
1 1
-(s-pLc)
-(c-1) -(p- S) c-1 p
p
 ':' 11 3 -

Ejercicio E2.17.
Determinar las funciones de forma del elemento mo-
nodimensional c 1 (viga recta a flexin) considerando:
a) dos nudos (figura E2.17a.).
b) tres nudos (figura E2.17b.).

Fig. E2. 17 a

o 2

1 ~-----

~,:+L/2-i
Fig. E2. 17 b
 - 11 4 -

a) Los grados de 1 ibertad en cada nudo son en trminos de des

plazamientos : la flecha o funcin (w) y la derivada de dicha
fun-cin o giro (ddw) particulariza das en el nudo.
. X

Se tiene:

4)
w(x)=(N , N , N , N w
1 2 3 0

eo

e,
S i'endo
w.= w (x) 1 x=x. y e.= ~(x) 1
(i=O,l)
1 1 dx x=x
1

con X =O y
o
Se comprende que N =N (x) tiene que cumplir las
1 1
condiciones:

N (x )=1 N (x )=0 N (x )=0 N (x )=0

1 1 2 1 3 1 4 1
dN dN dN dN
1 2
-1
dx x=x =o ~~x=x =l -1
dx
3
x=x =O -1
dx
4
x=x =O
o o o 0

N (~ )=0
1 2 N (x )=0
2 2
N (x )=1
3 2
N (x )=0
4 2

dN dN dN . dN
1 2
-1
dx x=x =O -1
dx x=x
=O -1
dx x=x
3
=O -1
dx
4
x=x
=1
1 1 1 1

que permiten determinar las funciones de forma N. (x) supo-

1
niendo que son funcionespo linmicas cbicas, es decir

(i=1,2,3,4)

Resulta:
N =
1
1-3~ 2 +2~ 3 N2 = L(1-2~+~ )~
2

X
N = 3t;2-2t;3 Nl= L(t;-1)~ con~='[
3

Estas funciones de forma constituyen las funciones

viga o polinomios hermticos.E stos polinomios se _definen del
 - 11 5 -

siguiente modo: (para m=1,2,3, ... n):

es un polinomio de orden 2n-1, que toma el valor

dkH
---( =1 para k=m y x=x.
1
dx

d KH
-dx K
=O para k1 rn y cualquier x=x.
J
(j~i)

o para k=n y x=x.

1

Se puede comprobar que

N1 = H(l)(x) N=H( 1 )(x)

00 2 10

N = H(o11. )(x) 1
Nl.= H(1 1 )(x)
3

b) De un modo anlogo se define:

w(x)=(N , N , N , N , N , N6 )
1 2 3 4 5

L
X =- y
1 2

Se puede demostrar que las funciones de forma sa

t(sfacen a las condiciones siguientes:

dN2i+1 o
N2 i + 1 (x i } ~1 dx 1 x=x.
_jK
o x.
1 J

N (x.} =O (con x=x. x.)

21 1 1 J
 - 11 6

dN dN.,.
2 -o
-dx1 x=x .=O
L
~ x=x~1
1 1 (i,j=0,1,2 ; jli)
J

Si se consideran polinomios qunticos, se obtiene:

N 1 =(1+6~) (4~ 4 -12~ 3 +13~ 2 -6~+1)=-24~ 5 -60~ 4 +66~ 3 -23~ 2 +1

N =
3
16~ 4 -32~ 3 +16~ 2

N 5 =-24~ 5 +52~ 4 -34~ 3 +7~ 2

N2 = L~(4~ -12~ 3 +13~ 2 -6~+1)

4

N4= L~(16~ -4o~3+32~ -a~)

4 2

.N
6 = L~ ( 4~ - 8~ +5~ - ~ )
4 3 2
siendo

se"comprueba que con la definicin de polinomios

hermticos anteriores se puede escribir:

N =H ( 2 ) (x)
2 10

N 3 -Ho 1
- (2)
(x) , N=H( 2 )(x)
4 11

N =H ( 2 ) (x) N =H ( 2 ) (x)
4 02 5 12

Se puede demostrar la siguien~e relacin existen-

te entre Jos pol inomicos de interpolaci6n de Lagrange y. los
h t: r mft i C O S :

1 -2 -{
d p n. (x) } (x- x . ) }
{ dx 1 x=x. 1
1

. (n)
H1 i = (x -X ) { p ? (X ) } Z
1 1

Observacin:
Es ms conveniente utilizar en el anlisis mediante
el m~todo de los elementos finitos sistemas de coordenadas con
 - 1 17 -

origen en el punto medio del eLemento y normalizado a longitud

2, es decir, en el caso de dos nudos t,; =-1 y t; =1. En el caso
1 2
de tres nudos t,; =-1, t,; =o y t,; =1. De esta forma se puede pro-
1 2 3
ceder a la integracin numrica de la matriz de rigidez) las
funciones de forma son ms simples. Se obtiene entonces:
Caso a):
( ~.+ 1 )
2
N
1
= (2 F,; - 1)

N = (~-1} (F,;+1) L
2
N = (2~+1)(1-F,;)
3
N
4
= (~+1)(F,;-1)L

Ca.so b):

2
N
4 =(1-~ )F,;

Ns= t<1+t,;)2t,;2(4-3t,;)

N6=-t(1+~)F,;2(1-~)
 - 11 8 -

Ejercicio E3.1.
Formar la matriz de rgidez de la estructura indi-
cada en la figura expresando por k~. un termino correspondien
lj -
te al elemento a. Se supone un s61o grado de 1 ibertad por nu-
do.
 - 1 19 -

<D
1
<Do

<Do

Fig. 3.1 b

Teniendo en cuenta la figura E3.1b. se tiene

1 2 3 4 5 6

1
e
K11+K22
e e
K1 a - K12
e
- -
2 e
Ka1
a e
K+Kaa Ka
14
a
K12+Ka2
e a
K a -
3 - K41
a Va
'44+ Kd2.2
a
K42
a
K4 a -
4 e a e a a b a b ,b
K2. 1 K2.1+K2a . K2 4 K2.2+K11 K2a+K1a 1'\.12
b
K22
5 - K:n
a Ka3 4 a
Ka2+K31
b a b
Kaa+Kaa
b
Ka2
'
6
- - - K2 1
b b
K2 ;
b f
K2 2+ K2 2
 - 120 -

Ejercicio E3.2.
Todas las barras del entramado plano de la figura
E3.2a. son inextensibles.
SE PIDE:
Obtener la matriz de rigidez total de la estructu-
ra. A efectos de este ejercicio se supone, por simplicidad,
que E1=1 y L=l para todas las barras.

TL

+l

ll
Fig. EJ. 2-a
 - 1 21 -

Los grados de 1 ibertad de la estructura se ven en

1 a figura E3.2b.
k .. , es la fuerza en el grado de 1 ibertad i cuan-
'J
do en el j se aplica desplazamiento unidad, siendo cero Jos
dems dep'Jazamientos.

__..,.
. 8
6El

- 7
4EI
L
cffi
l
Fig. E3. 2.b
Fig. E2.2c

Si en una barra se da una deformacin unidad se

tiene el resultado de la figura E3.2c.
E1 =l para todas las barras.
'L

Segn esto:

k11=12; k 21 = 2., k31=0 k41= 2.' k51 = O' k61=0; k 71 = O; ka,= 6

k22=16; k32=2 k42= O, ksz= 2.' k62=0; k72= o ,. kou2 = 6
k33=12; k43= O' k 53= o.' k63=2; k73= o., k83= 6

k44= 8 .' k 54= 2.' k64=0; k74=-6; ks4= 6

k55=12; k65=2; k 75---6, kas= 6
k66=8; k76=-6; ka6= 6

k77=72; k 8 7=-36.
k 88 =36
 - 122 -

La matriz de rgidez es:

12
2 16 S
11>-) El
t,..
o 2 12 le-
o
2 o o 8
o 2 o 2 12
o o 2 o 2 8
o o o -6 -6 -6 72
6 6 6 6 6 6 -36 36
 - 123 -

Ejercici-o E4.1.

Fig. Eit. 7 a

En la celosa plana de la figura E4.1a., todas las

2
barras son de seccin constante A=20 cm y mdulo de elasti-
6 -2
cidad E= 2,1x10 .kg,cm
La barra 2-3 sufre un incremento de longitud de 3
cm producido por la temperatura.
Se pide:
1).- Matriz de rgidez de la estructura.
2).- Movimientos de todos los nudos.
3).- Esfuerzos en todas las barras.
4).- Reacciones en los apoyos.
 - 124 -

El apoyo del nudo 4 no es corcondante con los ejes

generales (x 1 ,y 1 ) de la figura E4.1b. Se podra utilizar la
t~cnica, de considerar ejes de nudo. No obstante aqu, se ha

procedido a introducir una barra ficticia de rigidez a esfuer

zo axil K muy grande y nul a flexin.

Y'

Las barras b, e y .f no trabajan

Fig. E 4 1 b

Teniendo en cuenta la figura E4.1b. se tiene:

t<=

El cambio de ejes generales a particulares es

(figura E4.1c.):

Px= Px cosa
Py 1 = Px sena.
J= cosaj
JPx
ry ~ena.
. Px+ T= lcosa]
~ena

Y'li x
Fig. E4. "1 e
 - 125 -

Las matrices de rigidez de cada elemento son:

2
(K'?' 2) -[cosa]EA[
- s e na --L cosa sena]_EA
---L [cosa cosasena]
- a
senacosa sen 2 a

(K'
22 a
) =~
L [{ 1]
[~ ~]
1 -'{3/4
a= 1 2 o o (K '
EA
) =-
4
, "11 eL
-vr4 3/4

r1
O ~]

t: :J
[: :J
[ 1/4 -fi/4 ]
a= 1 2 Oo , K 11 1 ) g =
EA
(
L- -lfi/4 3/4

K = K ; a=60; K' = ~A [
1<.1/4 K./3/4]
K.l3/4 K. 3/4

La matriz de rigidez total es:

3/2 o -l
o 3/2

K'= EA
L -1 o
o o
 126 -

Teniendo en cuenta que el apoyo 2 no p~ede tener

1
desplazamiento vertical (d y =O) se tiene:
2 .

r-
" 3/2 - r-
d1
-
p1 o -1 o x2
x2

L p1
l--y2_ ~= EA
L
o 4K
2
o
K+ S
o d~2 1
p1
x4 -1 o -4- Y.3(K-1) di
x4
p1 o o YfcK-1) 1(K+1) di
y4 y4
~
- '- - '- -
La fila y columna recuadradas se pueden no tener
en cuenta en principio.
Siendo ~
2 _ 3 = 3 cm. de alargamiento
EA EA
F= r- ~2-3= r- .3

/
1 -

Estado final

Estado micial
V
Estado . modal
 - 127 -

Para la resolucin del problema el estado real o

final se descompone en los estados .indicados en la figura
E4.1d.
En la solucin modal se tiene:

p 1 =-1 EA
x2 2 L

de donde se deduce:

-1 o di
x2

EA EA K+5 di
r- o = L
-1 -,- Yf<K-1) x4

o o 1( K-1) ~
t< K+ 1) di
y4

invirtiendo queda:

~~~=~(4K+1)

di t(2K+1)
3
'lf(K-1)
-v'3:
lf:(K-1) -3/2
x2
di 8
- 3 (4K+l) 3 (K+ 1 ) i( K+ 1 ) .:1f( K-1) o
~K-1)
x4
di
y4
- 3 f{K-1) -,-
3 K+7 o

di
x2
-i-(2K+1) d 1 =- 3 2K+1
x2 '4'K'+T
1
0
x2
=-12
di 8 d1 K+ 1 e u a n do K-+co d1 __ 3
= -i-(K+1) --
x4- 34K+1 x4- 'lf
x4 3(4K+1)
di 3fi(K-1) d~4=
fi ffiT
K-1 diy4 =Yf
y4
8

2).- MOVIMIENTOS DE LOS NUDOS:

d~
d~
2 =-1.5 cm. 4 =-0.75 cm.
d 1
y2
=o d~ =
4 0.433 cm.
 - 128 -

3).- ESFUERZOS EN LAS BARRAS: Soluci6n final= Soluci6n inicial

+Solucin modal.
Hay que recordar que la solucin inicial solo va
a afectar a la barr~ 2-3.

dx

=78,75 ton.
~ =- 7 8 , 75 ton.
-fiJ [-3/2J=-3EA
2 0

..
P or 1 a so 1 uc1on 1n1.c1a 3 EA = 3 1 ~e .ton.
. 1: L

=-78,75ton .
EA[l 'V3J[-3/4]_3EA_
Barra g: P 19 = ( K 11 ) 9 d 19 =r- 2- 2 V3/4 --~ -- 7875
' ton.

Fijndose en los ejes locales de barra y llamando:

Traccin (+) } ,a la traccin positiva y a compresin negati-
Compresin (-) va se puede hace~ la tabla E4.1a.
BARRA SOLUCION SOLUCION
INICIAL MODAL ESFUERZO
1-2 o -7'0,75 -7'0,75
1 -3 o o o
2-3 -315 +78,75 ...;.236.25
Tabla E4.1a.
2-4 o +78,75 +78,75
3-4 o o o
3-~ o o o
4-5 o +78,75 +78,75
 - 129 -

A) Clculo de las reacciones (en toneladas).

En los ejes generales (x ,y): R.= reaccin en nudo

-1
i en ejes gener~

1 es.

~1= E.1a= Ta.(~12)a.!_! ~2=-~ v3t 1"] 3

[-1,5.10-2
o
=
4 4

1 -2
-4. 1 , 5. 1 o 3~,375
=1,05.10
fi -2
= 68,2
.- 4' 1 "'5. 1 o

EA
L
[1'I' + 41 +1
o ~l
2 -2 -
d' gL [: :] ~4=1,05.10
2 [1 o'5 OJ[1,5]
1, 5 o

2 [~ ~ J [-0,75 J
-1,05.10 0,433

Teniendo en cuenta, adems, la solucin inicial:

R = 315 r1~2] + 105 [-2,25+0,75] = [ o .]

[Vf
l .
[ l-
-:-2
o - 2 7 2 , es

S. in[ic~all
-2 S.modal
s.modal
[ -1/2 -1/2
3c 315
,!!_ Vf +E:Icc315 Vf +K2 1lc!!.2=315
1
-1 95[~13
. 1i
-v'3/4l
3/4 J
-1,5]
[ O
. . -1/2] -105 -0,375] = [118,125]
= 315
v';l/2 +o,6495 204,6
 - 130 -

r- 1 o J[- 1 5] +105 1 , 25
- V3~ 4 ] fo, J 0,37!:)]
~
75 =1 05 =
+1 05 L
0 o [
-V3 1+ l0,433 [ 0,6495
1+

[
3 9, 3 7 5 J
68,1975
La ventaja de operar con ejes de nudo estriba en
que, en el apoyo 4 (no conordante), la reaccin ya se ob-
tiene en dichos ejes y no en los generales tal como sucede
aqu.

.[1 14 . - fi 4] l r-0,75 J -0,375 J=

R.= p' =(K' ) ~d'=-105
-5 -2g 21 g -4 -4 3/4
0,433 =-lo!) [ 0,6495

[
39.375 J
-68.1975
Como comprobacin se puede real izar el clculo de
la figura E4.1e.
 - 131 -

272.8 ton. 236.2 S

/78.75
l \ o

?8 75 ton.
~-78.75
0*-0 ~236.25 ton.
k 1
78.?5
\
236.2'5
APOYO APOYO 2 APOYO 3

?8.?5 ton
78.75
. b_oo \
78.75 .,.._._
--lf
o
o-
\
78 75
~S ton.
6
APOYO 4 APOYO 5

Ftg E4 1 e

El resultado final se representa en la figura E4.1f.

+ 78.?5

- - - COMPRESION
TRACCION

Fig. E 4 1f
 - 132 -

Ejercicio E4.2.
Determinar las reacciones de apoyo del emparrilla-
do de la figura E4.2a. suponiendo que las barras son prism-
ticas y tales que EI/GJ= 10. en los siguientes casos:
1) Nudos y 3 coaccionados todos los. movimientos.
2) El nudo 3 se 1 iberal iza permitiendo la rotacin a la
torsin de la barra 2-3.

DATOS:
p=100 t L= 1O m

Q) 2
/
L

F19. E4.2a.
 - 133 -

a)-Matrices de rigidez en ejes locales.

La numeracin de los grados de l1bertad de las ba-
rras 1 y 2 se representan en la figura E4.2.b. (a) y (b).

~----------------~2
(j)

2
al. EJes locales b) Ejes tocales

e ) Ejes globales

Fig.E4.2b

Las matrices de r .g id ez son entonces:

GJ o o o o
:- o' 1
o -r6
(k22) =
~
12EI 6EI
-7 = - El
L
o 12
~
o _6EI
L2
4EI
--r o -r6 4

con i = 1 '2.
Se adoptan como ejes globales 1 1
1 2' y 31 1 los re-
presentados en 1 a figura. E4.2b. (e)
 - 134 -

Las matrices de transformaci n son:

[~ ~]
T
-1
= o -2
T = o o
o o
o -1 o o
Las matrices de rgidez en ejes globales son:

El
( ,, 1 )
:''22 1- (k22) =
-
L o, 1 o o
o 12
~
-r6
6
o L
4

(k 1
22 2
) = TT=
l(~22)2 -2
El
L
4 --6L o
6 12
o
L
~
o o o. 1

La ecuacin de equilibrio del nudo 2 es:

con p=
.:....2.
o y u=
-z eX
p d
y
o ez

es decir
El
L 4, 1 --6L o eX = o
-r6 24
~
-r6 d
y
p

o 4,1 o
--6L e2

WL 2
,
eX = 6 0,227 PL ...'>
26,4EI
d
y
4, 1 L 0,1S5L
ez 6 0,227

y las reacciones en el extremo 1 son:

 - 135 -

o bien
El p L2
M
X
= L -o' 1 o o 2.6,4[1
6 = -O, 023.PL
p 12 6
-o,:;ooP
o 4, 1 L
y - L~ r
Mz o --bL 2 6 -O, 4 77 PL

y las reacciones en el otro extremo se deduce mediante sime-

tra. En el caso numrico P=lOO t y L=10 m se obtiene:

1
eX = IT 2273 M
X
= -23
d
y
1553 p
y
-so
ez 2273 M2 -477
2 1

b) La matriz de rgidez de 1a barrq 2 (gdl activos)es:

El
L o' 1 -o' 1 . o o ex3 = o

-o' 1 o' 1 o o e x2 M
x2
.

o o 12
2L --6 L
d
y2 p y;

o o --6L 4 e z2 M
z2

y eliminando el gdl e por condicin esttica se deduce:

x3

El o o o o o o
'L e
X
= M
X (~22)2= El
L
12 12
o
L2 -r6 d
y
p
y
o
~
--bL
o --6L 4 ez
2
M
z
o --6L 4

en ejes generales

TT= El 6
(~22)2= !.2 (~22) i -2 L
4 --
L
o

--6
L
12
L2
o

o o o
 - 136 -

La ecuacin de equilibrio en el nudo 2 es:

WL 2
El
L 4,1 --6L o eX = o eX = IT 0,235

6 24
L ~
-r6 dy p ; d
y
0,161L

o --6L 4 ez 2
o e.,...
2
0,241

y 1 as reacciones son:
El
-o, 1 o PL 2
M
X
=
L
o IT 0,235 = -O, O2 3 PL
12 6
p
y o - L2 '[ 0,161L -0,48bP

M
z o --6L 2 o, 241 -0,484 Pl

.,
\n ..
p
y
=
El
'L o --12
L2
6
'[ o o -1 IT 0,235 =

M
z 2
o -r6 2 o o 0,161L
o o 0,241

PL o 12
- L2
6
L
-0,241
0,161L
= [0,522P
-O,
] 4 9 PL
o --6L 2 0,235

Los resultados numricos son:

:J
1
IT 235 M
X
= -23 y M
X
= o
1 61 o p
y
-48,6 p
y
-52,2

ezj 2 4.1 M
z
-484 M
z
-496
2
 - 137 -

Ejercicio E4.3.
Determinar los esfuerzos en todas las barras de la
estructura representada en la figura E4.3a.
E=3x1~ 6 t/m 2
Barras de hormig6n de 0,50 m. de canto y 0,20 m.
de ancho.

r2 :o S cm~
E0 : 2 x 10 7 t m~ 2

10 t/m

.r---s:oo
. 1P= 100 t.
ALZAD~
--- 5.00
.1
1

\
T o
o
tJi

PLANTA
l
Fig. E4 3a
 - 138 -

Considerando la figura E4.3b. se tiene:

"-..X..

Fig. E~.t 3 b

BARRA @ ~;r.1-----'-___;;;0---'----2 );,

z =z
Fig. E4.3c

r GJ
M
x2 'L o o X
' o o X

p
y2 = o 12EI
7 --
6EI
L2
y . = o o y

M
z2 o - 76EI 4EI
--:- z' o o z

(K22)a Ta

(K22)a
 - 139 -

Barra b
x o o X

y = o o y

BARRA @ ~:r- ___;;;.__:@~--2

3 _

Z E X'
;;.z z o o
----r--o
z

Fig. E 4.3 d

(KZ2)b= (K22)a

(K2_2 )b = T b' (K22)b .T~= --r

4EI
--
6EI
L2
o

6EI 12EI
o
-7 7
GJ
o o L

BaTra e ~tirante.

x o o X

=-1 o o y

BARRA
tirante
~4 -
2
z
};-;-y
=z'
z' o o
~

T
z

e
Fig. E 4.3e

Ean o o o o o
--r Ean
o o o (K2.2)c = Te' (K22)c. T~= O -:- o
o o o o o o
 - 140 -

La matriz de rigidez total es:

GJ + 4EI 6EI o
L L" -7
6EI E8 n 6EI
~+
K= -7 L3 L" -7
6EI g + 4EI
o -7 L L"

SUSTITUCION DE LOS DATOS

6
E
G--2(1+v) {E=3x10 G=l ,25x106 t/m2.
' v=0,2 +

?.
J= a~b{0,33-0,22~} con a<b a=0,20
{ b=0,50

Luego:
.
-2
t.m
L=5 m
1=20,38x10-l. m4 5,24 -1 '5 o
n=Sxlo- 4 m2 + K =1 o~ -1 '5 3,2 -1 '5
G=1,25xlo 6 tm- 2
o -1 '5
J=9,68x10-l. m4
E =2x107 tm- 2
a
 - 141 -

Las acciones externas se consideran como se indica

en la figura E4.3f.

10 T/m

--

.
~

Solucin inicial

Solucin modal

Fig. 4.3 f
 - 142 -

.
Por tanto, el vector de cargas , en la soluci n mo-
dal, es:

t=
.
[
M x2] = [ -125O
py2
Mzz 20,8
l P 1 = K 1 .d 1

. 0=5,24x -1 ,5y

-o, 125=10 3 {-1 ,5 e~ 2 +3,2d~ 2 -1 ,5 e~ 2 }-o, 125=1 ,Sx +3,2y-1 ,5z

20,8 =10 {
3 -1 ,5d~ 2 +5,24ez 2 } 0,0208= -1 ,5y+5,2 4z

x= e =-14 56x10-3 . r d.
x2 '
y= d y2 =-50 , 8Sx10- 3 m.
1

. -3 rd.
z= ez2=-1 0,59x10

LOS ESFUERZOS EN LAS BARRAS son:

GJ
-r- o o --
GJ
L
o o

!1z= o --
12El 6EI
L3 7 .!S.z1= o --
12EI 6EI
L3 -7

o - 76El 2EI
"T"
o
7
6EI
-2EI
L

p =
Soluci n modal -1 !12"~2
.E..2= K22 .!!2

p =10 3
-1
-~' 242
o
-El,6
-1 '5
o
2,5
1
1 , 5 x1 O-. 3
-14,56
] 3' 52]
-5O , 8 5 = 14,63
-10,59 49,80

0,242 o o -14,56 -3,52

p
-2.
= o 0,6 -1 '5 -50,85 = -14,63
o -1 '5 5 -10,59 23,33
 - 143 -

El resultado final se halla sumando las soluciones

inicial y modal (figura E4.3g.).

10 t/m 2
49.8 23.33

152-E-
1
o- ---------~ 2
...:. 3.52

t 14.63 @ 11..63 l
Soluctn tnictal Solucin modal

10 t/m

706 mt
C 52 ryH l J J JJJJJ l J J~t\2.53mxt
- t J
139.63 t 10.37 t

Fig E4. 3g

Barra b. Slo hay esfuerzos debido a la solucin modal

l -1]o [-14' 56]

0,242
. t o
f.1.= f1 2 T b ~l = O -0,6 1 -50,85 =
. o -1 '5 o o -10,59

-0,6
o -2,5b]
8,67
-1 '5
[39,88
o 1o,sgJ L,::;bl
0,6 -50,85 = -8,67
-1 '5 [ -14,56 [ 3,48
El resultado se expresa en la figura E4.3h.
 - 144 -

3.48 m e 2.S6mt 3o------------o2 2.56 mxt) 39.8 m xt

t8.67 t - @ 8.67 t l
F1g E4 3 h

Barra e e
.E.2 = ( K2 2 ) e ' T ~ '~2

[~ ~] [~ o] [5o, a:] = o o 7]
o o o -1 -14,56 o
1,
p =
-2
o o o -50,85 = o o -14,56
o o o o -10,59 o o -10,59 o

E 1 resultado se representa en la figura E4.2i.

fm1.7t
'

Fig. E4 .3 i

La comprobacin equilibrio fuerzas nudo 2 se hace en la figu-

ra E4.3j.

103719
/.
8.671/ ,'
Fig. E4.3j
 - 145 -

Ejercicio E4.4.
La estructura emparrillado plano de la figura E4.4a.
-1
estsometida a la carga vertical uniforme de 2 t.m so-
bre la barra 3-5. Se pide:
Calcular el movimiento del nudo 3.

FJg.E4. 4 a

Todas.las barras tienen igual luz= 6 m. y caracte

ristic~s elsticas:
E=3x10 6 t.m- 2

4
1=0,0050 m (i"nercia a flexin)

J=0,0030 m4 (inercia a torsin)

v=0,20
No se considera la deformacin de cortante.
Los nudos 2 y S estn totalmente empotrados y el
nudo 1 articulado a flexin (el giro a torsin y movimiento
vertical estn impedidos), y en el nudo 4 est impedido el
giro a flexin y el movimiento vertical.
Se utilizar en la resolucin de esta problema los
siguientes procedimientos:
a) Mtodo general matricial, con introduccin a posteri~
ri de las condiciones de apoyo.
b) Condensacin esttica de los g~ados de 1 ibertad no
activ6s a nivel de elemento.
 - 146 -

a) Mtodo general.
Se considera la nomenclatura indicada en la figura
E4.4b.

y___x ~1--~--~--~3----~~--~4
Z' @

Fig. E4. 4 b

El problema se resolver descomponindolo en dos

estados inicial y modal (figura E4.4c,).

Final In retal Modal

Fig.E4.4c

SOLUCION FINAL = SOLUCION INICIAL + SOLUCION MODAL

f.., (,!S.ll) a (JS.1z) a -o ~1
':_3 = (K' )
21 a C.!S.2z)a + (JS.ll ) e + (JS.1z) e ~3
(.!S.2 2) b + (.!5.1, )-d
4 b (JS.2 1 ) e (K' )
-22 e ~
 - 147 -

con d.=
-r
eX i y

d2i
e
yiJ

_1
tey, 2
t8y2
P.= G
-1 X i ex, dz 1 dz 2 ex2
p
zi Fig. Et. t.d
G
yi

-GJL o -GJ
'L o o

K11=. o 12EI
7
6:11 K12= o 12EI
--y 6EI
2
L 7
o
7
6EI .;1 J o --
6EI
L2
2EI
--:-

l
-GJ GJ
r o o
'L o o

6EI
K21= o -12EI K22= o
12EI -6EI
7 7 7 7
O
-6EI 2EI o -6EI 4EI
--:- --:-
7 7
E 3x10 6 3 6
G-2(1+v) = 2x 1 , 2 = 2,"1+ X 10

g= 3xto 6 3xto- 3 3 xto 3

L ~X 6 = P'
 - 1 48 .-

12El
--=
12x3x10 6 x5x10 -3 = ix 1 o3
L3
6x6xb 6

6EI
-=
6x3x10 6x5x10 -3 = .2.x10 3
L
2 6x6 2

2El
-=
2x3x10 6x5x10 -3 = 5x1o 3
L 6

4EI
-=
4x3x10 6x5x10 -3 = 10x10 3
L

Barras b y d:

[~ l~ l~ ~] [~:]
x=y' o
z=z' -+

y=-x' o

~~ ~] [ -~]
o o o
Tt = T = o
o 1 o

t
(JS.2z) b= .!. (!22) b Tt (JS.1, ) d= .!. (JS.,,) d".!.

-~] ~]
o o 4
no o o o o
(JS.22)b= o 10 3 5
--
5 o =
6 2
o o 5 1o -1 o
-2

o o -1 o o 3 1o 5 o
'4,1r 2
10 3 o o 5 5 o =10 3 5 5
2 b 2 b o
o o 5
-10 -2 o o o 3
'lf;lr

o o -1 3 o o o o -1
(JS.11) d= o o 10 3 no 5 5 o o =
o o b 2
o 5 1o -1 o o
2

o o -1 o o 3 1o 5 o
-2
10 3 o o 5 5 ~ 5
-2 b o =10 3 -2
5 o
o o 5
b
-10 o o o 3
2 '1f,lr
 - 149 -

La matriz de rgidez de 1 a estructura es:

G
x1 n
3 o o -
' ' .
o o o o
ol
DI
!..z 1 o 5/6 5/2 o -5/6 5/2 o o o ~1

G'y1 o 5/2 1o o -5/2 5 o o o

-3 o o 51 o o -3 o o
n n n
20
f.3
=103 o -5/6 -512 o 6 o o -5/6 5/2 ~3

o +5/2 o o 51 o -5/2
5
n 5

GK l. o o o -3 o o 3 o o
n n
P'
-z4
o o o o -51~ -5/2 o 5/6 -5/2 ~4

G'
y4
o o o o 5/2 5 o -5/2 1o

lntroducin de 1as condiciones de contorno

e x1 =o ex4f 0 (M x4=O)
NUDO 1 d
z1
=O NUDO l. d z4=0

e y 1f o (M y l =0.) ey 4 =o.

Se pueden uti 1 izar e 1 mtodo general pero en este

caso se consideran muelles ficticios para simular estas con
diciones de borde cinemticas. Sea K un nmero grande.
La matriz de rgidez queda:
 3 o o 3 o o o o o
o+K -4.lf
'
5
o 1 1 o b5 +K 2 o -5/6' 5/2 o o o 119..1

o 5 o o o o
z 10 -5/2 5

-3 o o 51 o o -3 o o
4,1f ~ 1i";1r

P 1= 1 o3
-3
o -5/6 -5/2 o 20
6 o o -!-j/6 5/2 u ~3 ..

o 5/2 o o 51 o -5/2 5 .....

5

o o o -3 o
2,4

o 3 o o
ll 1
V1
o

-q,lf n
o 1 1 o o o o -5/6 -5/2 . o t +K -S/ 2 11 ~4

o o o o 5/2 5 o -5/2 lO+K__, ,_ _,

los grados de 1 ibertad de los nudos 1 y 4 se pueden eliminar de forma particionada.

Se obtiene as:
-1
O= K11 ~1 + !12 ~3 ~1= -.!5_11".!5.12"~3

O= .!5_32"~3 + !33~4
} -1
~4= -.!5.33 .!5.32 ~3
} t3= ~21"~1 + .!5.22"~3 + ~23"~4
 - 1 51 -

3 +K o o
"li78
!,,- 10 3 o 5
b +K 5/2 .
'
t.,,= 10 3 (lOK 2 +, ... )

o 5/2 1o

Dado que l(+oo slo interesan 1os menores que contengan trmi-
nos en K2
o o o
1 o- 3
.!S.,,- ;o- o o o ..

o o

3 o o
p
- 3 o 65 +K -5/2 ., = 1 o3. ( K2 + )
K
33 -1 O ,
!:J. .
'
33
o -5/2 1 O+ K

[~ ~]
o
-1 =10-3.438 o
. .!S.3 3
o

3
-n o o o o o -, o o

o -5/6 -5/2 o o o o -5/6 5/2 =

o 5/2 5 o o o -5/2 5

o o o -3 o o o o o
D"
= o o -512 o -5/6 5/2 = o 25 ~2
4 2

o o 5 o -5/2 5 o _25 25
2
 - 152 -

3 3
-n o o o o -4,lr o o

o -5/6 5/2 o o o o -5/6 -5/2 =

o -5/2 5 o o o o 5/2 5

3 o 3
o .-4."8 o o -L2 o o
-4."8
' ' 4,8
= o o o o -5/6 -5/2 = o o o

o o o o 5/2 5 o o o

o o o 51 o o
n
p
-3
= 10 3 -o' 1 o 25
4 --
25
2
o b
20 o

o 25 25 o o 51
-2 n
9 o o
~2
' t:S

--
4,8
3
o o o ~3
o O o

99 o o ex 1 3
Mx 1 3 4,lr
32,5
p
z'3 = 10 3 o 12
2' 5
2' d
z'3
{ 4. 4a. }

M o 2,5 45 ey , 3
y'3 2 2,5

Invir-ti-endo se obtiene:
eX 1 3 0,04849 o o M
z'3

d
z'3 = 10- 3 o 0,38095 -0,02540 p
z 1,3
e-y. 3 o -0,02540 0,05503 M
y'3
 - 153 -

Solucin final = Solucin inicial + Solucin modal

P=2x6_ 6 t M -6
x'3
2
2 T /m p
P= = 6
z'3
r:lllll lllllll l -2
M=2xo_ 6 mxt.
12 My'3 o
G \..._31 . ~S
~ 6.00------l

tp

e. 0,04849 o o -6 ex, =-0,29D94 rad.

x'3 3
.10- 3
10- 3 o m.lo- 3
dz'3 = 0,38095 -0,02540 6 dz'3 2,2857

6y'3 o -0,02540 0,05503 o ley'3=-0,1524 rad.

1 o- 3

(b) Condensacin esttica a nivel de elemento.

En general si en la matriz de una barra (elemento)
se p a r t i e i o na d e a e u ~ r do e o n 1o s g r a do s d e 1 i be r ta d a e t i v os
. (a) y b no activos (b), se puede escribir:
{4.4b.}
. [ f.a
p
J = [!.a
k
a ~abJ
-b- -ba ~bb
-1
Si !.b =O, se tiene ib=-kbb d
kba -a y 1 a matriz de r
gidez condensada es:
k*a a = k a a - k ab - 1
kbb kba {4.4c.}

En este caso la matriz de rigidez de las barras del

emp a r r i 1 1a do es :
Barra b y d:
5 3 -2.2.
(.!S,2 2) b = 103. 1o 2 o (!S_ll)d= 1 o 1o o

5 S o 5 5
...... o
2 6 -2 b
3 3
o o
O"
o o
n
 - 154 -

Barra a
r 3
( K')a= 10 3 p o o o o

o 5
~

b
5
2
o -i 5
2

o i
2
1o o --25 5

o o o o

o 5 -i2 o 5 _,2.
-b 6 2

5 o 5 1o
o 2 5 -2

En el nudo
de e~ta barra: (nudo 1 general):
9x 1 =0 y dz 1 =o. En cambio e giro de flexin no est
coaccion~do y My 1 =o. Se tiene pues en este caso, segn
la nomenglat ura de {4.4b}:

= 3 o o
k
-a a n ~bb= 10

o 5 5
b -2

o --5
2
1o

K
-ba= ( o ' -25 ' 5)= kT
-ab

Por lo tanto, aplicando la frmula de condensac in

esttica{4 .4c} se tiene:

k-;'t
-a a
=
n
3 o o o 1
lo
(o, --25 ,5)=
n
3 o o

o 5
6 --5
2
5
-2 o 5
24 -4
5

-4 .!.i
_,2. 1o o 5
o 2
5 ...
'
 - 155 -

Barra e
(K') =(Ki) Es una barra idntica a la a
- e a

En este e as o en e 1 nudo 2 (nudo 4 gen era l.) de 1a

barra, existen las siguientes coacciones al movimiento:
dz 4=o y ey 4=o. En cambio el giro de to~si6n no est~

coaccionado y Mz
4 =o. Por lo tanto, se tiene, en este caso:
3 o o
k =
-a a n
o 5 5
b 2

o 5 1o
2

= {- 3 o, O}= kT
1
~ba n ab
Po 1 o tanto:
3 3
k*a a =
n o o -n 32{
--r
3
-n o, O}= o o u

o b
5
--5 2
o o 5
b
5
-2

o 5
-2 1o o o --5
2
1o

La matriz de rigidez de la estructura, se obtiene.

mediante adici6n de las matrices condensadas anteriores_y re
su 1 ta:
K= 10+10+.
'
--
5 5
2 2
o =

---
5 5
2 2 i+t+.J+1
b 2 b

o 5 5 ~~.ll.r1o
2-14 !f,O !f,O 2

99 o o
q
= o 32,5 2;5
12 2'
o 2,5 45
2' n
 - 156 -

Esta matriz coincide con la obtenida en el mtodo

general. Frmula{4,4a.}
 - 157 -

Ejercicio E4.5.
En la estructura articulada de la figura ~4~5a.
todas las barras son de seccin constantes y tales que
En/L=103 tm-l. La barra 4-5 sufre un incremento de la lon
gitud de 5 cm. se pide:
1) Desplazamientos en todos los nudos referidos a los
ejes (x 1 , y 1 ) .
2) Esfuerzos en todas las barras.
3)Reacciones en los apoyos.

Y'

T,..,

1 X'

~-- 4.00
.. ---+-1---- 4. oo ---1
Fg.E4.5a
 - 158 -

El apoyo del nudo 5 no es concordante con los ej~s

generales(x',y') de la figura E4.5a . Se podra utilizar la

tcnica de considerar ejes de nudo. No obstante aqu se ha
procedido a introducir una barra ficticia de rigidez a es--
fuerzo axil K muy grande y nula a flexin.
Se conside~an los sentidos de la figura E4.5b.

Y'

~.
x

Fig.E45b

A efectos de trabajo quedan las barra~ de la

figura E4.5c.

(~l2)f
!~ + (~22) f+
(~22) b+
(~ l1 ) d

F1g. E4. 5 e

.Las matri~es de rigidez en ejes locales son:

. . EA
(f11) e= (~11) f= (~22) f= (~22) b= (~11) d=L

( .!S.1 2 ) f= <.!S.z 1 ) i== - ~

 - 159 -

x=xcosa
X'= X cosa y 1 =ysena

v~ Y'= Y cosa
x]= [cosa] x
[ y sena

X'
Fig. E4 5 d

Considerando el cambio de base de la figura E4.5d.

se tiene:

cosa] eos 2 a cosase~a J

-
T= [ sena ; .I.T=[c?sa sena]; K =K
_senacosa
(por ser K un nmero).
1
r sen a

Barra e:
-------
{ cosa=4/5
sena=- 315
(K )
-1 1 e
= 1 o3 ~
-E
12
12]
-E
2; .
Ba r r a f : { sena""' O( ~2 2 ) f= ( ~ l1 )
co~a=1
f= 1 0
3[O1
(~1 ~ ) f= ( ~ 1 ) f= 1 o3 [
-1
o

{coscx=O
Barra b :
sena=l <~22) b= 1 o3
[~ ~ J

Barra d:
4
cosa=-
(K 1
-11 d
) =10 3
_g]
~~ rr25
12
25
9
{ sena==3/5 -E

Barra m

{ cosa=4/5
sena=3/5
(K 1 ) =10 3
- m
[~K
25
J2K
25
;-;-K
12
9K
25
l
 - 160 -

41
E -12/25 -1 o

-12/25 9/25 o o
K'-103 41+16K
-1 o 25
.!l.(K-1)
25 ..

o o 12 34+9K
'25 (K-1) 25

r- - r-
12
- - -
P'
-x4 41/2, -E -1 o d'
x4
P' _.!2 ~K o o d~4
1 -y4 25 25 1
P' - 10, -1 o 41+1bK .!l.(K-1)
-x5 25 25 d~5
P' o o .!l.(K-1) 34+9K
-y5 25 25 d~5
.... - .....
PI
- .... -
-y 4 { -~ 1 +d 1 }
K 25 x4 y4

Cu a n d o K-+oo ..
d 1
y
4 =O 1
.

La fila y columna recuadradas se pueden no tener en

cuenta.
-Para la resolucin del problema se tienen en cuenta
los estados inicial y modal de la figura E4.5e.
 - 161 -

A= S cm.

Final

= S cm. S

Inicial. Modal

F19 E4. 5 e

Si la barra 4-5 aumenta 5 cm: EAi:,P

L -'

P= ~A .5= 10 3 ~ .5.10- 2 m~ 50 ton.

Luego, al aumentar la barra en 5 cm. de longitud,

aparecen unas fuerzas de empotramiento que se indican en la
figura E4.5e.
E1 es t a do moda 1 de 1 a f i g u r a E4 5 e p ropo re i o na.
'
el siguiente resultado:
 - 162 -

-.so 41 o
-1 d'x4
E
= 1 o3 -1 41+16K .ll.(K-1)
50 25 d'
25 x.S

o o ~; (K-1) 34+9K d'y5

25

2._ = i!._ ( 4 1+ 1 6 K) ( 3 4+ 9 K) - i!._ 12 2(K-l)2_3~;9K =

10 3 25 3 25 3

= ~{1394+144K 2 +931K-144K 2 +288K-144}--1- (21250+5b25K)=

25' 25 3
= ~~ {49241K+51250-5625K-21250}+
3 ~= 10 ~{43616K+30000}
L;;l 25
-
34+91<
:;_ ~ = ~{ 1 2 O1 K+ 1 2 5 O} 25 . 25
-ll_(K-1)
25
34+9K 492
2!) - 2(K-1)
25
:12(K-1)
25 ~(656K+1 025)
25
-
d'
x4 -.so
3
d' = 25 50
x5 3
10 {43616K+30000} ~A
d'
y5 o

Teniendo en cuenta que K+oo, se tiene:

25 3 1201K+1250 34+9K 1 3
d~4= 4361GK+30000{(-.SO)
. 25
2 + 5 2 5
25
2
1o
=- 2 8.10- m. }--r
25 3 34+9K . 41 (34+9) -3 -3
d~.s= 43616K+30000{(-50).25. 252 +50. 252 }10 =4,13.10 m.

3
d~5= 43616~~30000{(-50) (-12):~~ .25-50.49~~~-1)}10-3=-5,5.10-3m.

DESPLAZAMIENTOS EN LOS NUDOS

d' =-2 8 cm d' =O d~ =0,41. cm
x4 ' ' y4 5 d's=:..o,s.s
y
cm.
 - 1 6 3 -

Clculo de los esfuerzos en las barras:

( p 1 ) f = ( k 1 1 ) f d 1+ ( K2 1 ) f d 2 = \ l( 1 1 ) f T T d l + ( K2 1 ) f T T d 2_ =

o] 14 1 3 1 ~ ; 3] =-28-4,13=-32,13 t
l:-s,s.1o

-280.10-3]
( p 1 ) e= (K 11 ) e TT d =
1
1
10 3 [ 4
5
-_]
5 [ =-22,4 t

(P2)b= (K22)b T
T
d'= 10 3 [o
2
~ 4,13.10
1] 5,5.10-3-3] =-5,50 t

3
[ 4,13.10- ]
(P1)d= (K 11 )dT
T
d'= 10
1
3
[~ -3
5
J -5,5.10 -3 = 6,6 t

Se tienen as los esfuerzos en les barras que tra-

bajan, pero en la barra f hay que tener en cuenta el est~do
inicial que es una compresin de 50 t.
Luego:
50-32,13 = 17,87 t. de compresin.
Llamando positiva a la traccin y negativa a la com
presin se tiene el resultado de la figura E4.5f.

Fig. E4.5f.- Esfuerzos en las barras

 - 164 -

La obtencin de ls reacciones por tcnicas matriciales

(en ejes generales x' ,y'., y con resultados en toneladas)
es:
16 -12 -2,8.10 -2.
R =
-2 .E.ze + .E., b--(1/
1 1 1 ) d 1 +tK 1 )
_;21 e .::.L \-1.2 b -5-
d'--10 3 25 rr
-12 9 o
25 E

-1o 3 o~ 1. 1o: :] =- 1o [- 1' 7 92] 10 . o ] [1 2]7' 9

[ 00 01] [-0,55.10 1,344 - [ -0,55 = -7,94

~3= .e.~d= <!i2,)d.~s=-lo

-~~
1256 - 25
12]
2~
0,41
] =-1o ] [- l
0,5264
=
-5,264

[*
-0,55 -0,3948 - 3,94ti

5
R " modal=P'
!!.4 -1 e +P'
-1 f ={(K'
-1 1 ) e +(K'
- 1 1 ) f }d'+(K'
-4 -1 2 ). f .d'=10
-5 _ 12
-~~] - ]+
9
2 6
'

rr 25 ~

+10 [-: :] ~::::] =10 t~::::] +1{:l10 [1.::.]= ~:~:.] .

Teniendo en cuenta la' soluci.n inicial:

-50 J
R = 50
-4 [01J + [ 13,44 = L;,44 J

41
n o.1T= 10 T3.736l.
+1 o [- !l
25
[-o,ssj ~0,9448J

.Teniendo en cuenta la solucin inicial:

- 1 J [ 3 7 ' 3 6' J =
-12,64]
( -9,45 ..
~s= 5 0 [ o + -9,448
 - 165 -

Al operar con ejes de nudo esta ~S ya sale en dichos

ejes, lo cual es mucho ms conveniente.
En la figura E4.5g. se representan las reacciones y
los esfuerzos en las barras:

A"-..
Apoyo 1 ; REACCJON :: O
1 22.4 t.
5.5 t.

Apoyo 2 : REACCJON Apoyo 3: REACCJON

S.St

"l~m?i
6.6 t.
22~v
-..... t7.87t
Apoyo 4; R EACCJON Apoyo S : REACCION

Fig. E4.5g
 - 166 -

Ejercicio 4.6.
En la estructura entramado plano de la figura E4.6a.,
las barras lA y A2 son de hormig6n armado de m6dulo de elsti
cidad de 200 tcm- 2 , rea 1 m2 e inercia 1 m La barr~ 2-3 es
4
2
un tirante sin rigidez a flexi6n; de rea 10 cm y m6dulo de
2
elsticidad de 2.000 t.cm-
Se pide:
1) Matriz de rigidez de la pieza 1A2.
2) Desplazamiento del nudo 2.
3) Desplazamiento del nudo A
4) Esfuerzos del tirante.

1
o
-
o

+-1
A 2

!
10 t.
o
o
-'

_L
1 1.00 1.00 ---t
FJg E4 6a
 - 167 -

a) La matriz de la barra 1A2 se_ puede determinar siguiendo

varios pr~cedimientos (condensaci6n est~tica, teorema de
Castigliano, etc.). Aqu se indica el mtodo del centro
elstico.

H/2

-+
X

. L
Y=-X+217'2"

Fig. E4. 6b

Teniendo en cuenta 1a figura E4.6b., 1 as coo.rd en a-

das de 1 centro elstico 0 son:
1
,. L
J~Ly 1o =J-vri[y
dx dx
2
x=O .
J
y= V'!
L L
=
L /2
~
= L
m

J* f~
dx dx
r;;-

El centro elstico ~st situadu en el punto medio

de la altura del tringulo.
S~ adopta como origen de coordenadas el centro
elstico.
 - 168 -

borras
infinitamente
r igidas

Fig. Et,. 6 e

Teniendo en cuenta la figura E4.6c., se aplica una

fuerza p en el centro el~stico, el movimiento que se produce
-.x
en l es:

1o= 1 cose= lV2

-
2

Las integrales se obtienen como sigue:

Por lo tanto

dx= k dx
11
 - 169 -

Anlogamente con una carga p.y' 5e produce un movi-

miento vertical d , de valor:
y

En este caso se tiene:

L L 2
~ +w 1!L dx 1 = n{Wo
1 "l
1 L L 2.}
x~dx 1 +
~ o no
V! Vi
Por lo tanto p = E't/2
y , L2 2
L\31 +r-)
o o
Para un momento G en el centro elstico, e 1 giro
z
que se produce en 1 , es:

~
dx dx 2.L
e z = Gz ! = G = G
~ z ~ z VzE 1o
vz
El 0
es decir Gz = V2 2 --:-- e z = k33 ez

La matriz de rgidez, del centro elstico es:

S = .--
K e

con s= ( px, py' Gz )T , e= (d x' ' d y, e ) T

z

K= diag.(f 11 , f22' f 3 3)

La matriz de rgid~z general de la estructura refe

rida al punto o1 ,se calcula considerando la matriz de equili-
brio ~=l ,
3 (matriz unidad de dimensin.3):

k kT =
= -21 -H K = -K, K
-11 = -HT K H= K
~
12 - - -

Estas matrices, se modifican al cambiar los ejes del

nudo 0 1 a los nudos 1 y 2, siendo las matrices de transforma-
cin estticas:
 - 170 -

T
-1
= o o
-T,=
... o o
1 o o o
v1 ~
2 2

k~.= T. k
Las matrices nuevas son -1 T 1 (i,j=1,2).
J -1 _ j _j
Resultando, referido a los ejes (x,y)

.Yff
4 11

kl =kl T =- o
-1 2 -21

f-1{ o

112
-4 f-11

24 6 6 6
En este caso, se tiene f 11 = -y .10 , f 22 =7.10 y
6
f3 3 =1.10 y la matriz de rigidez queda:
 - 1 71 -

24 o ---=-
6Vz
--=
24 o 6Vz
7 7

o 6
7 7
3l2 o --76 +3v'2
7

6Vz
---:- 31Z 11. 6v'2 3V
---:- -1
7 7 7

k=
24
---: o 6V 24 o 6V2
---=-
7 7
,......
o 6
-y -- 31.12
7
o 6
7 --
3V2
7

6Vz 61/2
7
3Vz
7
-1 --
7 --
3v'2
7 7
13

2) Movimientos del nudo 2.

Fig. E4. 6 d

La matriz de rgidez de la barra 2-3, corresponde

a una barra de entramado con ausencia de rgidez de flexin,
.es decir 1=0

~]
o

n
Solo se necesita conocer =
(kll)b= o
o

~]o
o

H
6
1o o
o
 - 172 -

La ecuacin resultante en desplazamientos es:

ll+vz o --y-
6112 d = -5 V2
7 x2

o 6 _3112 d -5v2
7 7 y2

6v'2
-- --
3112 13
ez2 o
7 7 7
Resolviendo el sistema, resulta

d' 2,965.10- 6 metros

=-
x2
d' =-12 50 .1o- 6 metros
y2 '
ez2 =- 6 ' 02 .10
-6 radianes

3) Movimientos del nudo A.

Los esfuerzos en el nudo 2 son:

( Pza=
') (k'22a
) d '= 1o 6 .1o- 6
2
24
7 o --
61/2
7
-2,965 =

o 312
-- -12,50
7
13 -6,02
7

-2' 875]
-7,06
[o
Estos esfuerzos trasladados al punto A resultan:

~~.19]
~ 2' 14
Valores que se representan en la figura E4.6e.
 - 173 -

12 .1i.
(\
---2.19
A

Fig.E4.6e

1
Q=f
0
(
12,1 4 -2,19x )dx ___
IT =El1 ( 12 , 14 _~)=
2
11,045
El
-
8=5,52.10 -6 rd.

df= Mxdx
El f= f: (1 2 , 1 4 X- 2 , 1 9 X ) ~
2
7 =rt (1Z 21 4 ~)=5,34
3 El
f:;:2,67.10 -6 m.

-6
dxA=1,88.10 m;

dyA=-1,88,10
-6
m;

-
ezA~-5,52,10
-6 rd.

4) Esfuerzo en el tirante;
(P2) b= (k22) b d2 +

p~2b f2 o o -2,965 -4,19

= 10 6 ,1o- 6 o o = o
p~2b o

G~2b o o o o

Este resultado est en ejes generales, En ejes lo-

cales es Pxzb= 4,19 t.+4,19.
Esfuerzo en. e 1 t i r a n te : t rae e i n + 4 , 1 9 ton ,
 - 174 -

Ejercicio E4.7.
Obtener mediante la utilizacin de un programa de
c~lculo matricial los esfuerzos en todas las barras de la -
figura E4.7a., siendo EITL =Cte.

~3 T
~toc+t.ooi

r-3.00 1200-+-200---+-- 3.00 ---1

Fig. E 4. la

NOTA.- Las especificaciones de acceso al ordenador y caract~

rsticas de entrada de datos para los ejercicios E4.7,

E4.8. y E4.9. se indican en el apendice a estos ejer-
cicios.
 - 175 -

DATOS
NUMERO DE NUDOS 7
NUHERO DE BARRAS 11
NUMERO DE NUDOS CARGADOS: 3
NUMERO DE APOYOS 2
MODULO DE ELASTICIDAD 20000000.

COORDENADAS NODALES
NUDO X y

.oo .oo
2 3.00 3.00
3 4.00 1 50
4 . 5. 00 3.50
5 6.00 1 5o
6 7.00 3.00
7 10.00 .oo

CONEXION y PROPIEDADES DE LAS BARRAS

BARRA PRIMER NUDO SEGUNDO NUDO AREA
1 2 .00212
2 3 .00214
3 2 3 .00090
4 2 4 .00103
5 3 4 .00112
6 3 5 .00100
7 4 5 . 00112
8 4 6 .00103
9 5 6 .00090
. 10 5 7 .00214
11 6 7 .00212

CARGAS NODALES
NUDO PX PY
3 .00 -3.00
4 . 00 -4.00
6 2.00 .oo
 - 176 -

CONDICIONES DE CO~TORNO
CODIGO VALORES PREESTABLECIDOS
(O:PREESTAbLECID0,1:LIBRE)
NUDO u V u V
o o .0000 .0000
7 o o .oooo .0000

RESULTADOS
DESPLAZAMIENTOS NODALES
NUDO u V
1 .oooo .oooo
2 .0005 -.0010
3 .0003 -.o. o1 3
4 .0001 -.0011
5 .0002 -.0007
6 .0001 -.0003
7 .oooo .0000

REACCIONES NODALES
NUDO PX py

4.2588 3.2000
2 -.0000 -.0000
3 ':",0000 -3. oooo
4 .0000 -4.0000
5 -.0000 .oooo
6 2.0000 .0000
7 -6.2588 3.8000

ESFUERZOS EN BARRAS
BARRA ESFUERZO AXIL
-3.6270
2 -1.8094
3 1.9815
4 -3.7766
5 .8005
6 -.9530
7 -3.7994
8 -.1.6562
9 2.3112
1o -4.2017
1l "'3.2876
 - 177 -

Ejercicio E4.8.
En la estructura de la figura E4.aa. el nudo 4 su-
-2
fre un descenso de 1 cm. y un giro de 0,5.10 radianes en
sentido positivo, y la barra 2-3 sufre un incremento de tem-
peratura de 40C, Obtener mediante la u~il izacin de un pro-
grama de clculo matricial los siguientes esfuerzos y reaccio
nes.
a) Reacciones en 1.
b) Esfuerzos en A (punto medio del dintel 3-5).
6 2
E= 2.10 ton/m

ex= 10-5 oc-1.

2
Seccin de todas las barras: 0,30x0,40 m

o
o
O')

l
_j_
.~

i-- 6.00 --+-~- 6.00

1
!

----.,..--- 6.0C -- -j
Ftg. E 4. 8 a
 - 178 -

Se considerar el punto A como un nudQ ms, el 5 de

la figura E4.8b.

iJ e q = 2 4 1 o- 4 ~n =
6
2 1 o- 4 2 1 o . ! 2 1o
-2
= 8 t

El estado a resolver es el de la figura E4.8b.

//
)/
A 1

1
3
r - - - - - - - r- 81

Flg.E4.8 b

DATOS
NUMERO DE NUDOS T
NUMERO DE BARRAS 6
NUMERO DE NUDOS CARGADOS: 2
NUMERO DE APOYOS 4
MODULO DE ELASTICIDAD 2000000.

COORDENADAS NODALES
NUDO X V
.oo .oo
2 6.00 6.00
3 12.00 6.00
 - 1 7 9 -

4 12.00 -3.00
5 15.00 8.25
6 18.00 10.50
7 18.00 -3.00

CONEXION y PROPIEDADES DE LAS BARRAS

BARRA PRIMER NUDO SEGUNDO NUDO AREA .. ~1. DE INERCIA
1 1 2 .12000 .00160
2 2 3 .12000 .00160
3 3 4 .12000 .00160
4 3 5 .12000 00.160
5 5 6 .12000 .00160
6 6 7 .12000 .00160

CARGAS NODALES
NUDO PX PV MZ
2 -8.00 .oo .oo
3 8.00 .oo .00

CONDICIONES DE CONTORNO
CODIGO VALORES PREESTABLECIDOS
(o: PREESTABLECIDO, 1:LIBRE)
NUDO u V RZ u V RZ
1 o o o .oooo .oooo .oooo
4 o o o .0000 -.0100 .0050
6 o 1 1 .0000 .oooo .oooo
7 o o .0000 0000. .oooo

RESULTADOS
DESPLAZAMIENTOS NODALES
NUDO u V RZ
.0000 .0000 .oooo
2 .0070 -.0069 -.0010
3 .0072 -.0099 -.0006
4 .oooo -.0100 .ooso
5 .0050 -.0070 .0021
6 .oooo -.0001 . oo1 9
7 .oooo .oooo -.0010
 - 180 -

REACCIONES NODALES
NUDO PX PY .MZ
-.8345 -.3259 1.8925
2 -8.0000 -.0000 00.00
3 ~.0000 .0000 .oooo
4 -1.4150 -2.1857 8.3697
5 -.0000 -.oooo -.0000
6 2.3516 .oooo -.0000
7 -.1020 2.5116 .oooo

ESFUERZOS EN BARRAS
BARRA . NUDO FX FY MZ
1 1 -.8205 .3597 1.8925
2 .8205 -.3597 1.1594
2 2 -8.8345 -.3259 -1.1594
3 8.8345 .3259 -.7960
3 3 -2.1857 L4150 4.3652
4 2.1857 -1.4150 8.3697'
4 3 -3.3066 -.6596 -3.5693
5 3.3066 .6596 1.0959
5 5' -3.3066 -.6596 -1.0959
6 3.3066 .6596 -1.3775
6 6 2.5116 . 1 02 o 1.3775
7 -2.5116 -.1020 .0000
 - 1 81 -

Ejercicio E4.9.
La estructura emparrillado plano de la figura E4.9a.
est sometida a una carga vertical uniforme de 2 t/m. sobre
la barra 3-5. Obtener los esfuerzos en todas las barras, median
te la utilizacin de un programa de clculo matricial.

Fig. E4.9 a

Todas las barras tienen igual luz = 6 m. y caracte-

rsticas elsticas:
6
Et= 3.10 t.m
-2

1= 0,0050 m4 (inercia a flexin)

4
Jc 0,0030 m (inercia a torsin)
v= o, 2 o
Los nudos 2 y 5 estn totaJmente empotrados y el nu-
do 1 articulado a flexin (el giro a torsin y movimiento ver-
tical estn impedidos), y el nudv 4 tiene impedidos el giro a
flexin y el movimiento vertical.
 - 182 -

Dado que el programa utilizado solo admite la posl

bi 1 idad de <..argas en los nudos, es preciso real izar la descom
posicin en estados de la figura E4.9b.

7:
-
1 A ~

l--x
Y' ...
;3
+ 1~---M~x~~3~~~~----0~y~3~----~
3

i
~
2 2
Estado i me 10 l E.stado modal

Fig E 4 9 b

Resultando:
Solucin final =solucin inicial +solucin modal.
La solucin inicial es inmediata y slo afecta .a la barra 3-S.
La solucin modal se obtiene mediante el programa de clculo
matricial.
DATOS
NUMERO DE NUDOS 5
NUMERO DE BARRAS 4
NUMERO DE NUDOS CARGADOS: 1
.NUMERO DE APOYOS 4
MODULO DE ELASTICIDAD 3000000.'
COEFICIENTE DE POISSON .200
 - 183 -

COORDENADAS NODALES
NUDO X z
1 .oo .oo
2 6.00 -6.00
3 6.00 .oo
l 12.00 .oo
5 6.00 6.00

CONEXION y PROPIEDADES DE LAS BARRAS

BARRA PRIMER NUDO SEGUNDO NUDO 1N. A FLEXION 1N. A TORSION
1 3 .00500 .00300
2 2 3 - .00500 .00300
3 3 l :oosoo .00300
l 3 5 .00500 .00300

CARGAS NODALES
NUDO MX PY MZ
3 -6.00 6.00 .oo

CONDICIONES DE CONTORNO
CODIGO VALORES PREESTABLECIDOS
(O:PREESTABLECIDO, 1 :LIBRE)
NUDO RX V RZ RX V RZ
.1 o o o .ooo .oooo .0000
2 o o o .oooo ;oooo .0000
l o o .oooo .0000 .0000
5 o o o .0000 .0000 .oooo

RESULTADOS
DESPLAZAMIENTOS NODALES
NUDO RX V RZ
1 .0000 .oooo .0006
2 .oooo .oooo .oooo
3 -.0003 .0023 -.ooo2
l -:.0003 .0000 .0000
5 .oooo .0000 .0000
 - 184 -

REACCIONES NODALES
NUDO HX py MZ
181 8 -.6667 -.0000
2 4.2597 -1.1775 .0952
3 -6.0000 6.0000 .0000
4 .0000 -1.5238 4.9524
5 -7.1688. -2.6320 .0952

ESFUERZOS EN BARRAS
BARRA NUDO MX FY MZ ..
1 81 8 -.6667 -.0000
3 -.1818 .6667 -4.0000
2 2 .0952 -1.1775 -4.2597
3 -.0952 1.1775 -2.8052
3 3 ,0000 1. 5238 4.1905
4 .oooo -1.5238 4.9524
4 3 -.0952 2. 6320 8.6234
5 .0952 -2.6320 7.1688

Los esfuerzos en 1a barra 3-5 sern, teniendo en

cuenta 1a figura E4.9c.:

Fig. E 4. 9 e

M =-6 + 8.6234 = 2.6234 m. t.

z35
M = 6 +
z53 7.168? =13.1688
p =-6 + 2.632 = 3.368 t.
y35
p
y 53
=-6 - 2.638 =-8.638
Todos los dems esfuerzos no varan.
 - 185 -

APENDICE A LOS EJERCICIOS 4.6, 4.7, 4.8 SOBRE. RESOLUCION DE

ESTRUCTURAS MEDIANTE PROGRAMAS GE CALCULO MATRICIAL.
1. GENERALIDADES
Los programas disponibles en el Centro de Clculo
de la Universidad de Santander son:
CELOSIA (Versin 1) para la resolucin de estructuras arti-
culadas planas, con las siguientes 1 imi
taciones:
Mximo n5mero de nudos: lOO
Mximo nmero de barras:lOO
Semiancho de banda mximo:20
Mximo nmero de apoyos:20
ENTRAMADO (Versin 1) para la resoluci~ de estructuras re-
ticuladas planas, siendo sus 1 imita--
cienes:
Mximo nmero de nudos: 50
Mximo nmero de barras:100
Mximo nmero de.apoyos:20
de banda mximo:21
-Semiancho
- - -- -
.. ..
EMPARRILLADO (Versin 1) para la resolucin de emparrilla--
dos, con las mismas limitactones-
que el programa ENTRAMA~O.
2. ACCESO AL ORDENADOR.
Para 1~ ejecucin de cualquiera de ~stos tres pr~
gramas deber disponerse la siguiente secuencia de tarjetas
( b: espacio en blanco; ~=cero).
TARJETA 1
JOBbCAMI .EST.~S

TARJETA 2
PASSWORDbESTRUCTURAS
TARJETA 3
DELETEbCEL.DAT si se trata de ejecutar el programa CELOSIA,
DELETEbENT.DAT si es en l caso del ENTRAMADO, o bien
DELETEbEMP.DAT si es en el EMPARRILLADO.
TARJETA 4
CREATE/IbCEL.DAT o bien
CREA TE 1 1b EIJT . OAT o
CREATE/IbEMP.DAT segn sea el caso.
 - 186 -

GRUPO DE TARJETAS 5
En este lugar ir una. serie de tarjetas con los da
tos del problema, segn estructura del apartado 3 del prese.!l
te apndice.
TARJETA 6
)
Esta tarjeta slo 1 levar un parntesis derecho en su prl
mera posi.cin.
TARJETA 7
XbC ELOS lA
XbENTRAMADO
XbEMPARRILLADO en cada caso.
TARJETA 8
RENAMEb CEL.OUTb (nombre del problema) o bien
RENAMEb 'ENT.OUTb (nombre del problema) o
RENAMEb EMP.OUTb (nombre del problema) segn los casos
TARJETA 9
QPRINTb (nombre.del problema)
TARJETA 10
DELETEb (nombre del problema)
TARJETA 11
BYE
TARJETA 12
END
3. ENTRADA DE DATOS.
La entrada de datos se real iza con formato 1 ibre por
lo que se puede emplear el nmero de decimales que se quiera,
debindose separar simplemente los nGmeros entre sf por una co
ma.
3.1. DATOS DEL PROGRAMA CELOSIA. (Versin 1)
Deber disponerse la siguiente secuencia de tarjetas:
1) TARJETA DE PARAMETROS BASICOS. Una tarjeta que contenga
el nmero de nudos, nmero de barras, nmero de nudos car
. -
gados, nmero de apoyos, y el mdulo de elsticidad del
material.
2) TARJETAS DE COOR~ENADAS NODALES. Habr tantas tarjetas. co -
mo nudos y cada una de ellas contendr el namero del nudo,
 - 187 -

su abscisa X y su ordenada Y.
3) TARJETAS DE CONEXION Y PROPIEDADES DE LAS BARRAS. Habr
tantas tarjetas como barras tenga la estructura, conte-
niendo cada una de el las el nmero de la barra, el nm~

ro de su primer nudo, el nmero de su segundo nudo, y

el rea de su seccin recta.
4) TARJETAS DE CARGAS NODALES. Habr que disponer tantas
tarjetas como nudos cargados haya. Cada tarjeta conten-
dr el nmero del nudo, la componente x de la carga y
la componente y, en ejes generales.
5) TARJETAS DE DATOS DE LOS APOYOS~ Tantas tarjetas como
apoyos existan. Cada tarjeta contendr el nmero del
nudo, el cdigo de la condicin de apoyo para la com-
ponente del desplazamiento en la direccin x, idem en
la direccin y, y los valores preestablecidos de las
componentes x e ydel desplazamiento.
El cdigo de condiciones de apoyos es el siguiente:
Cdigo indica que el desplazamiento correspondiente no
esta coaccionado, es decir, no es conocido.
Cdigo O ndica que el desplazamiento correspondiente esta
coaccionado y es conocido (Desplazamiento impues-
to o apoyo).
En el caso del cdigo O, se dar el valor del movi-
miento (O en el caso de apoyo fijo y.el valor especifi
o -

cado si existe un asiento).

En el caso del cdigo 1, cualquier valor del movi-
miento no se interpreta, ya que se supone desconocido.
Se aconseja suponer . el valor O al movimiento.
 - 188 -

3.1 .1.EJEMPLO
La secuencia de tarjetas de datos correspondientes
a la estructura E4.7. serfa (cada lTnea corresponde a una tar
jeta):
7 , 1 1 ' 3 ' 2 , 20000000. (1)

1, o. , o.
2, 3 , 3.
3 , 4. , 1.5
(2)
4, 5 . 3.5
t

5~ 6. , 1.5
6, 7 t 3.
7, 1 o. , o.

1, 1, 2t .00212
2, 1 , 3, .00214

3 ' 2, 3, .0009
4' 2, 4, .00103
5 ' 3, 4' .00112 ( 3)
6, 3 , 5' .001
7 , 4, 5 ' .00112
8 ' 4, 6, .00103
9, 5, 6, .0009
1 o. 5. 7. .00214
1 1 ' 6, 7 t .00212

3t o. , -3.
(4)
4' o ' .- 4.
6, 2. t o.

1' o' o, o . ' o. (5)

7 ' o' o' o. ' o.
 - 189 -

3.2. DATOS DEL PROGR-AMA ENTRAMADO (Versin 1)

La sec~encia de tarjetas ha d~ ser la siguiente:
1) TARJETA DE PARAMETROS BASICOS. No vara.

2) TARJETAS DE COORDENADAS NODALES. No vara.

3) TARJETAS DE CONEXION Y PROPIEDADES DE LAS BARRAS. Es ne

cesario dar el momento de inercia de cada barra. Habr
una tarjeta por cada barra, conteniendo el nmero de
la barra, el nmero de su primer nudo, el del segundo
nudo, el rea de su seccin recta, y el momento de
inercia.

4) TARJETAS DE C~RGAS NODALES. Para cada nudo cargado ha-

br una tarjeta que contenga el nmero del nudo, y las
cargas P , P- y G c~rrespondientes a ese nudo.
X y

5) TARJETAS DE DATOS DE LOS APOYOS. Es igual que para el

programa CELO S 1A pero e o n s i de r ando t res g r a dos de 1 i be!.
tad en lug~r de dos. Por cada apoyo habr una tarjeta
conteniendo el nmero del nudo, los cdigos de las in--
cgnitas u, v y e, y Jos valores preestablecidos, si
procede, de las mismas. El sistema de cdigos es simi-
lar al del programa CELOSIA.

3.2.1. EJEMPLO
Para la estructura del ejercicio E4.8. la entrada de
datos sera la siguiente:
7, 6, 2, 4, 2000000. ( 1)

1, o., o.
2, 6., 6.
3 ' 12 ' 6
4,12.,-3. ( 2)
5,15., 8.25
6,18.,10.5
7,18.,-;a.

1, 1, 2, . 12, .0016

2, 2~ 3, .12, .0016 (3)

 - 190 -

3' 3 ' 4 , 1 2 ' .0016

4 , 3, 5 ' 1 2 t 001 6
( 3)
5 ' 5 ' 6 ' 1 2 ' .0016
6' 6, 7, 1 2 ' oo1 6

2,-8., o., o.
3, 8., o., o. (4)

1, o' o' o, o . ' o . ' o .

4, o, o, o ' o.,-.01, .005
6, o, 1 , 1 , o . ' o. ' o (5)
7 ' o, o ' 1 ' o . ' o . ' o .

3.3. DATOS DEL PROGRAMA EMPARRILLADO. (Versin 1)

Las vartaciones con respecto al ~rograma ENTRAMADO
son mnimas.
1) TARJEtA DE PARAMETROS BASICOS. Es necesario dar el va
lor del coeficiente de Poisson. Esta tarjeta contendr el n-
mero de nudos, nmero de barras, nmero de nudos cargados, --
nmero de apoyos, mdulo de elasticidad, y coeficiente de Poi
sson.
2) TARJETAS DE COORDENADAS NODALES. No vara.
3} TARJETAS DE CONEXION Y PROPIEDADES DE LAS BARRAS. Habri
una tarjeta por cada barra conteniendo el nmero de la barra,
el nmero de un primer nudo, el del segundo nudo, su momento
de inercia 1, y su inercia a la torsin J.
4) TARJETAS DE CARGAS NODALE~. Una tarjeta por cada nudo
cargado, que contenga el nmero del nudo, el valor de G , el
X
de P y , y el de M
z
5) TARJETAS DE DATOS DE LOS APOYOS. Por cada apoyo habr
una tarjeta conteniendo el nmero del nudo, los cdigos de las
incgnitas 8 , V y 8 , y los Valores preestablecidos, S pro-
X y Z
cede, de las mismas. El sistema de cdigos es similar al del
programa CELOSIA.
 - 191 -

3.3.1. EJEMPLO
~ara
la estructura del ejercicio E4.9. la entrada
#
de datos ser1a la siguiente:

5, 4, 1, 4, 3000000., .2 (1)

1' o. ' o.
2' 6 ' -6.
3' 6. ' o.
4' 12 , o.
5' 6. ' 6. (2)

1' 1' 3 , .005, .003

2, 3 ' 4, .005, .003
3' 2 ' 3' .005, .003
4' 3' 5 ' oo5' .003 ( 3)

3, -6., 6., o. (4)

1' o, o' 1 ' o. ' o. ' o.

2' o, o ' o, o. ' o. ' o.
4, 1 , o , o, o. ' o. ' o.
5' o' o' o' o . ' o. ' o . (5)
 - l92 -

Ejercicio E5.1.
Determinar los movimientos en todos los nudos de la
estructura entramado plano de la figura E5.1a.
Todas las barras son de hormign con E=3.10 6 t.m- 2
y seccin 0,50.0,20 m.

Pieza rgida

o
o
~

~
.....~ - - 5 . 0 0 --+1.00+-- 5.00 ~

F1g. E 5 1 a
 - 193 -

Se tomarn como nudos los expresados en la figura

ES. lb.

Fg. ES./ b

Teniendo en cuenta la figura E5.1c. y planteando

el equilibrio.se tiene:

Gz 1c P'x2 = p xlc p~2

~
P'2 = p ylc P'y2 =

tPylc G'z2 = Gz1c

+P -P
ylc xlc G'z2
Gz 2
~ ~ ~~~~~~~ o o P- x 1 e

t Py2 o o p
y1c
-1 G
Fig. ES. 1 e
zlc

Teniendo en cuenta la figura E5.1d. la compatibill

dad se puede expresar como sigue:

1""'

xlc = d'x2 -e'z2 .1

d el
x le
d ylc = d~2+6~2.1 d ylc =

6 zlc = 6~2 6 z1c

 - 194 -

~ o -1 d1
x2

t = o d'y2

o o 1 e~2
*t
TN

F1g. E5.1d

Con estas matrices de transformacin se puede pa-

sar de ejes locales a generales. Para el ensamblaje d~ben de
observarse las relaciones.

p x 1 e= { K 1 1 } e d x 1 e+ { K 1 2 } e d x 2 e -+ [ ~
p 1 e] U OO2 ~ ~ { K 1 1 } e T ~ t, ~2 +
{K''" }
11 e

T~{K12}c.T~ ~J

{ Kl2} b o
{ Kll} b

{Kl2}b+
K { Kl2} b
{Kll}c

o
 - 195 -

3.10 6 .0,1 .0,2= b00.10

. 2
--
EA
L
6
3.10 .0,5.0,2
L
L5
{ L=6 10 6 .0,5.0,1 500.10 2

6 1 3
12EI 12:3.10 "f20,2.0,5 rL=5 3.0,2.10 -3 .10 6 =6.10 2
-=
0 3
L3 L3 L6 3g .1o 5 -3,47.1o 2

6 1
6.3.10 120,2.0,5 3 L=5 0,1.3.0,5.10 4 =15.10 2
6EI
-=
L2 L2 { L6 0,53 6 2
0,1 .3.-rb .10 =10,4.10

6 1
4. 3 1 o 12" o, 2. o, 5 3 L5 0,2.0,25.10 5=50.10 2
--
4EI
L L
{ L=6 0 53
3. .10 5=41,7.10 2

6 1 .. 25.10 2
2.3.10 "12"0,2!0,5 3 L5
--
~

2EI { L6
L L ~ - 20,8.10 2
600 o o o -1 o
{K22}a=102. o 6 -15 {K } .. 1 o o
22 a
o -15 50 o o 1

o 1 o 6 o 15
-1 o o 2
{K22}a .. 10 ' o 600 o
o o 15 o 50

{Kl1 }b{K11 }b= 1 02 [6~0 o

6 1:]
15 50

-6oo o o

~
{K~2}b{K12}b102! -6 15
l -15 25

t6~o
o
2
{K21}b{K21} b= 1 0 -6
15
-1:}
25
 - 196 -

-1~]
o.
{K22}b{K22}b= 102 [600~ 6
-15 50

o
{K } =10
[600
2 O 6 o]
15
11 e
o 15 50

o o o -1 600 o -600
2
{K 1
11
}
e
= o o ~ K11 }e o 1 = 1o . o 6 21
-1 1 o o -600 21 686

{K12}c=102. l6oo~ -6
-15
o
1
25
~ l
{Kl2}e=
-1
o
D

1 . 1
o
o { Kl 2} e o
o
o o

o
o =1 o. 2
1'
-600
o
600
o
-6 .
-21
~
140 l
r-600O
-1~]
o
{K21}e=10 2 -6
L O 15 25

~] {K21} e
o o
{K 1 } = o o -:]
1 1
.
=1 o .
2-
-600
o -6
600]
-21
21 e
o o o 1 o 15 40

[600~
-1~]
. o
{Kh}c={K22}c=102. 6
-15 50
 - 197 -

{K22}d=1Dz. [500~ -1~,4]

o
3,47
-10,4 41 , 7

[ ~. 47 1~ 4]
{K2zldc[~ ~]
-1 o
2
o
o
JK22}dt o
o
=10 .
1 o' 4
.500
o 41 , 7

606 o 15 -600 o o o o o
p =O o 606 15 o -6 15 o o o ~1
-1
15 15 100 o -15 25 o o o
2 -600 .:.o o 1200+K o -600 -600 o o
o 1 o . dx2
-Joo o -6 -15 o 12 6 o -6 15 d
y2
o o 15 25 -600 6 736 600 -21 40 ez2
o o o -600 o 600 603,5 o 1 o, 4
p =O
-3
o o o o -6 -21 o 506 -15 .2.3
o o o o 15 40 1 o, 4 -15 91 , 7

Aplicand o condensa cin esttica :

. O= Kll E., + K12 ;!.2 -1

d =-K 11
-+
-1 -
K12 d?...
-1
-1
O= K32 !2 + K33 -1 !.2={- K21 Kl 1 K12+K22-K23K33K32}d2
~ -+ ~3=-K33 K32 !z

En lugar de utilizar la constan te K para introduc ir

1as condicio nes de contorno , es mejor modific ar 1a matriz de
rgidez (dx2=0):
606 o 15 o o o o o o
o 606 1 5 o -6 15 o o .O
15 15 1oo o -1 5 25 o o o
o o o 1 o o o o o
2
1o o -6 -15 o 1 2 6 o -6 15
o 15 25 o 6 736 600 -21 40
o o o o o 6oo 603,5 o 1 o, 4 .
o o o o -6 -21 o 506 -15
o o o o 15 40 1 o' 4 -15 91 , 7
 - 198 -

.;..1 1
[ 99.63
0,37
0,37
15] --1
o o
[o
Kll = 60150 99,63 -1 5 ; -K21'Ktl'K12= ~ -2,2816 3,833
-15 -15 606 3,833 -6,48

-1 1
76,51 -0,258 -8,72
K33= 46o84,s. -0,258 91 1 52 1 5 -K 23 .K -1 .K 32
[o o
= O -2,483
o
-4,975
33
-8,72 1 5 506 o -4,975 -606,65

4,86 1d:2
o o d
o d 2 x2 =O
o
~
-too 7,235
o o 4,86 122,87 . ez2
J .
r-1]=[7,235
lo 4,86 . [:::]
7,235d y 2 + 4,86e z 2 =-1. d x2 = o
~

4,86d 2 + 1z2.87e = o
Y z2 d
y2 = -14,2 cm.
e z2 = 0,00562 rd

o -0,0037 0,00614 o 0,006

d 1=-K -1
11 K12 d 2 = o 0,0062 .-0,0186 -0,142 = -0,001
o 0,1496 -0,2481 0,00562 -0,0227

d x1 = 6,0 mm.
d y1 = -1 , o mm.
e z1 = -0,0227 rd.

[o o,oo28 -o, 9887 lo,

-1
d 3 =-K 33 K32 d 2 = O 0,0~7
o -0,1627
0,032
-0,3188
0.~14+
0,00562
oo6] 0,0008
0,0213.

d x3 = -6,0 mm.

dy3= -o,a mm.

e z3 = 0,0213rd.
 - 199 -

d
1e =
o
o
o
1
o
-:] d
-2
= o
o
o1
o
-1H
1
1
o
o.o1~2

o.oo562
l -0,0056l
= -0.1364
0,0056J

d
x1e
= -5 ' 6 mm.

d
Y1 e =-13.64
. cm.

e 1 e = o' oos 6 rd o
 - 200 -

Ejercicio E9.1.
Determinar las funciones de forma de los elementos
con continuidad C0 , de tipo rectangular de la figura E9.1a.

1) 2} 3)

~V +v V
1
!
2 1 2 1 2
!
u u

3 4 3 4 3 4 -
S S

4) 6}
~ V 4V
7 7
2 1 2 1 1

u 6 u 8 6

3 4 3 4 3 4
S S 5

Fig. E 9. 1 a
 - 201 -:-

1)
V
1
N1 = 'Zi(l+u) (l+v)
T2 1
1 1
N2 = 'Zi(l-u) (l+v)
+
1
u
N
3
1
= (1-u) (1-v)
-L3 4
J--1 +
Fig. E!J.lb
~
N4= j(1+u) (1-v)

t
1
u*

N2 =T1 (1 - u )( 1 V )

Fig. E!J. 7 e

1
N1 = 'Il(l+u) (l+v)

2)
Para el nudo 5(0,-1), los valores de
V
las funciones de forma anteriores son:
2 1 s_ 1.
NS-o
,- ,
N5-o
2- t
N5_1
3-2 y N4- 2
u
2
La funcin de forma N es: N = 1-u ,
5 5
3 l.
5 por lo tanto las nuevas funciones de
Fig. E9. ld forma son:
 - 202 -

1
N = 1"(1+u) (1+v) - N5 N = 1 (1+u) (l+v)
1 1 5 4

1 1
N = 1(1-u) (1+v)
2
- N5 N = 1"'(1-u)
2 5
(1+v)

3
1
N = 1"'(1-u) (1-v) - N53 1 1 2 1
N = -(1-u) (1-v)--(1-u )=-(1-u) (1+2u+v)
5 4 2

1 1 1 2 1
N4= 'Ii(1+u) (1-v) - N5 N = (1+u)
4 5
(1-v)-2(1-u )=-(1+u) (1-2u+v)

2
N = 1-u
5

Fig. E9.1 ~

. 3)

V Para el nudo 6(1,0), los valores de la

2 1 funcin de forma anteriores son:
6 1 6 6 6 1 6
N1=z , N2 =0 , N =0 , N4 =2 , NS=O
6 3
u
La funci6n de forma nueva N6 , es: N6 =1-v 2 ,
3 l.
S por lo que se obtienen las nuevas funciones
Fig.E9. 1 f
de forma:
 - 203 -

1 6 1 1 2 1
N = (1+u) (1+v) - N N =(1+u) (l+v)-2(1-v )=-(1+v) (1+u-2v)
1 1 6

1 6 1
N = "4(1-u) (1+v) - N N =4(1-u) (1+v)
2 2 6

N 3 ~-.)
1 (1+u) (1+2u[l v)

N =-4(1+u) (1-2u++-v)
4
-

- N
N~ N 6 =-.)(1~u)
6 1 2
(1+2u+v)

2
4 N 6 =- 4 (3-u+v-2u +uv-v )

1
N3 --l1- uH1 2 u vl N = 1- v2
4 .
Fig. E9.1g

Para el nudo 7(0,1), los valores de las

.. 4)
funciones de forma anteriores son:
V

7
2 1

u La nueva funci6n de forma N es:

- 7
3 4 1 2
N =2 ( 1 - u ) ( 1+ v ) , p o r 1o q u e s e o b t i e n e n
5 7

Fig- E9.1 h
las nu~vas funciones de forma:
 - 204 -

1 7 1 1 2 1 .
N = 1f(1+u) (1+v) - N N =(1+u) (l+v)- (1-u ) (1+v)='Iiu(1+u) (1+v)
1 1 7 4
1 7 1 1 2 1 .
N = 4(1-u) (1+v) - N N =4(1-u) (1+v)- (1-u ) (1+v)- u(1-u) (1+v)=
2 2 7 4 4

- 41 u ( 1.- u) ( 1 - v)

1
N =-4(1-u) (1+2u+v) -
7 N =-4(1-u)
N
1 1 2 1
(1+2u+v)+4(1-u ) (1+v)=- u (1-u) (1-v)
3 3 7 4

N S~ 1 - u ~ ( 1- u2 )
2 2 2
- N; N = 1- u - ( 1 + v) =t ( 1 - u ) ( 1 - v )
7

1 2
N
7
= 2 (1-u ) (1+v)

-r 3
1
..L

Para el nudo 6(1,0), los valores de las

5) funciones de forma anteriores son:
V
? 6 1
2 1
N4=2
6 u
La funcin de forma nueva N
6 es:
31 L.
S
por lo tanto. las nuevas funciones de for-
F1g. E9. 1
ma son:
 - 205 -

1 6 1 1 21 2
N = 'Iiu ( 1 +U ) ( 1 +V ) - N N ='Iiu ( 1 + U) ( 1 +V ) - l ( 1 - V ) ='Ii ( 1 +V ) ( U +U - 2 + 2 V )
1 1 6

N =-i-u (1-u) (1+v)

2
- N~ N =-i-u (1-u) (1+v)
6

N3 =-tu(1-u) (1-v) -N; N6 =-i-u(1-u) (1-v)

N = tu ( 1 +U) ( 1 -V)
4 - N~ N b = tu ( 1 +U) ( 1 -V) - ~ ( 1 -V 2 ) =t ( 1 -V) 2
(u"!'" u - 2-2 V)

N =-11-ul(l-vl
S 2
Fig.ESJ.1 k

Para el nudo 8(-1,0), los valores de

6) las anteriores funciones de forma son:
V
8 8 1 N85 =O
2
7
1 N1 =--21 , =-
N2 2 .,
8 6 u

3 4
S
La funcin de forma nueva N es:
8
Fig. E9.7t 1 2
N =2 ( 1 - u ) ( 1 - v )
8
 - 206 -

1 -2 8 1 2 1 2
N = 'Z"(l+v) (u+u -2+2v)-N N ='i(l+v) (u+u -2+2.v)+ (1-u) (1-v )=
1 1 8 4

1 2 1 .
q ( 1 +V ) (U - 1 +V+ UV ) =1" ( 1 +V ) ( 1 +U ) (U+ V- 1 )

1 8 1 1 2 1
N =-'Z"u{1-u) (1+v)-N N =-'liu{1-u) {l+v)-'1"(1-u) (1-v )=-'1"(1-u) (1+v) (u-v+l)
2 2 8

. 1 8 1 1 2 1
N =-'lfu(1-u) (1-v)-N N :=-'Z"u(1-u) (l~v)-4(1-u) (1-v )=- {1-u) (l+v) (u+v+l)
3 3 8 4

1 2 8 1 2 1 . 2
N4= 'lf ( 1 -V) (U+ U - 2-2 V) -N 4N =4 ( 1 -V) {U+ U -2- 2 V) +Lf { 1 -U ) ( 1 -V ) =
8

q1 ( 1 -V) (U 2 - 1 -V- UV) =4

1
( 1 -V) (U+ 1 ) (U- V -1 )

N
5
= 21 (1-u 2 ) (1-v)

2
N = 1-v -r(1-u)
6
(1-v 2 )=~(1+u) 2
(1-v )

N}= ~(1-u 2 )_(1+v)

1 2
Ns=.'2(1-u)(1-v ).

::=
Fig. E9_ 1m
 - 207 -

Ejercicio E9.2.
Determinar las funciones de forma de los elementos
con continuidad C0 , de tipo triangular de la figura E9.2a.
Utilizar coordenad~s naturales o triangulares.

Fig. ES. 2a
 - 208 -

1) Las funciones de forma son:

N1=L 1
1)
N2=L2

3~1
N3=L3

Fig. E9. 2 b

Fig. E9. 2 e

2) Para el nudo
2) 4(1/2,0,1/2), los
2
valores de las

3Dl 4
Fig. E9.2 d
funciones de for
ma anteriores son:
4 1
N=- 4
1 2 'N 2
4 1
=O ,N=-
32

La nueva funcin
de forma N4 , es:
N = 4 L L
4 1 3
 - 209 -

Por lo tanto
4
N1= L1 -N 1 N4= L1. - 2 L1 L = L1 ( 1 - 2 L )
3 3
4
N2= L2 -N2 N4= L2

4
N = L3 -N 3 Nl= L -2L 1 L = L (1-2L 1)
3 3 3 3

N4= 4L 1 L3

2 2

Fig. E!J. 2 t

Para el nudo 5( 1 , 1 , O), los valores de las

2 2
anteriores funciones de forma, son:
N5=l N5=l N5=0 Nl5 =o
12' 22' 3 '

La nueva funcin de forma N , es:

5
N5= lL1L2

Por lo tanto:

N1= L1 (1 -2 L ),:-N 51 NS= Ll ( 1 -2 L ) 2 L1 L2 = Ll ( 1 -2 L2 2 L )

3 3 3

N = L2 -N 52 N =L 2 -2L 1 L2 =L 2 (1-2L 1 )
2 5
 - 210 -

3)
2
N4= 4L 1 L3

N
5
= 4L L
1 2
3D, 4
Fig.E9.2f

N 1 :: L11 1- 2 L 2 -2 L 3 )

Fig. ~!). 2 g

4)

4) Para el nudo 6(0, 1/2, 1/2), los va-

2 lores de las anteriores funciones de

forma son:
6 61 61 6 6
N1=0 ; N~= : N3 =y ; N4=o ; N5=o la
3 . 1.. l nueva funcin de forma N6 , es: N6 =4L L
. 2 3
F1g. E!) 2 h Por lo tanto:

N1= L1 (1-2L 2 -2L )

3

N2 = L2 (1-2L 1 )-N 62 N6 =L 2 (1-2L 1 )-2L 2 L =L 2 (1-2L 1 -2L )

3 3

N = L (1-2L 1 )-N 6 N6 =L (1-2L 1 )-2L 2 L =L (1-2L 1 -2L 2 ).

3 3 3 3 3 3

N
4= 4L L
1 3

NS= 4L L 2
1

N
6= 4L L
2 3
 - 211 -

Ejercicio E9.3.
A) Deducir 1 os po 1 i nom i os de Lag range de orden 1,
2 y 3.
B) Aplicacin a la obtencin de las funciones de
forma de los elementos finitos tringulares e~

C) Comprobar en el elemento triangular de elastici

dad p 1 a na 1 os re que r i m i n tos de e o n ve r g_e n e i a

Fig.E9. 3a
 - 212 -

A) Una funcin polnomica de grado n, f(x) que toma

los valores f. para x=x., con i=1,2 n, es:
1 1

n
f(x)= .L P~(x) f(x)
1=o 1

siendo: donde:
p~(x)
n 1
P. (x)= ~n~-
1 P. (x.)
1 1

. (x-x )(x-x.)..~.{x-x.)
0
p ~ (x) = . ---- ___ L_ . _______ .L. --
1 x-x.
1
Si 1 os puntos x , x , x x n es t n a i n ter va 1 os
0 1 2
iguales, en un intervalo L se obtiene con x.=h.
1 n
1)1er d 1 x 1
- o r e n {lineal) f.o(x)=1-r =po

L--1
____
.......,.. ~
i
p 1 (x) =~ =p 1
-1 L 1
o
Fig.E9.3b

2) 2~ orden
(cuadrtico)

p
2
O
(X ) = {2 ~-
L
1 ) ( ~- 1 )=( 1 - 2 p l ) p 1
L 1 O
1--
1
l. --...;-- .h. ___j
2 2 -----

o .i 2
Fig. E 9. 3 e
p (X
2
2
) =~
L
( 2 ~- 1 ) = p ( 2 p -
L
1
. 1
1
1
1)

3) 3=..!:.. 'orden (cbico)

p 3 (X)
1
=--
9
2L
X ( X
3--2
L
) (X
L
9 1 ( 2-3p 1 ) p 1
--1 ) =-p
21 1 O

3( ) 9 X ( X
P 2 X 2L 3--1L
=-- ) ( X) 9 1 (
1--
L
=-p
21
3p 1 -1 ) .p 1
1 O
 - 213 -

B). Se deduce la aplicacin pedida de la tabla si-

guiente:
Polinomio de Lagrange. Funcin de forma: 3

1) Lineal
N =L
1 1

N2=L2

N3=L3

2) Cuadrtico Vertices Fig. ES. 3d

2 1 1
N ~(2L
p 0 =(1-2p 1 )p 0
1 1 -1)L 1
p2=4plp1 N =(2L -1)L 2
1 o 12 2
2 1 1
p2=(2p1-1)p1 N =(2L -1)L
3 3 3
Puntos medios de lados
N =4L L
4 1 2
N =4L 2 L
5 3
N6 =4L L
3 1 Fig. ES. 3e

3) Cbico
3 1 1 1 1
p o=2 ( 3 p 1 - 1 ) ( 3 p 1 - 2 ) p o
3 9 1 1 1
p 1 =2p 1 ( 3 p o- 1 ) p o

p ~ =~ p ~ (. 3 p ~ - 1 ) p ~
3 1 1( 1 2
p 3=2p 1 3 p 1 - 1 ) ( 3 p 1 - 2 )
Puntos en lados
9 L (3L -1)
N4 =rL 1 2 2
1
9
Ns=!LiL 2 (3L 2 -1)
N6 =fL 2 L (3L 2 -1)
3
N =~L 2 L 3 (3L 3 -1)
7 Fg. E9.3 f
Na=~L 3 L 1 (3L 3 -1)
. 9
N =rL L1 ( 3 L1 - 1 )
9 3
Nudo interior

Para este nudo interior no existe comparacin

directa con .Lagrange.

 - 214 -

C). u=(N 1 ,N 2 ,N ,----N 6 ) u

3 1
u2
u3
u4
u5
u6
v=(N ,N ,N ----N 6 ) v
1 2 3 1
v2
v3
v4
vs
v6
Los requereim ientos en convergen cia son:
1) Movimient os de slido rtgido a) u(x)=cons tante=u 0
b) u(x)=-e 0 y
e) v(x)=v
0
d) v(x)=e 0 x.
Po r 1o t a n t o a ) u 1 =u 2 =u ==u 4= u =u 6 =u 6
3 5 0
lo que i mp 1 i e a u ( x ) = (N l , N , N - - - -N ) . u O
2 3 6

6
Es decir l:N.=l
1. 1

b) v 1 =v 2 --------vi =v 0 Anlogame nte se deduce

6
l:N.=l
1 1

C) Ul = - e 0 y 1 ' U 2 = - e 0 y 2 - - - - - - - - . U6l:= ~ ~ 0 y 6
Por 1o tanto
y= N.y.
. -1
I
1 1
d) v,=eoxl-6 v2=eox2- -------- vo=eox6
es decir x= l: N x
i=1
1 1

A continuac in se com~rueban las igualdades an-

teriores (L 1+L +L =1):
2 3
 - 215 -

6 3 3 2 2 2
1=.L: N.=.L: (2L.-1)L.+4.r L.L.=2L +2L +2L -(L +L +L )+
1=1 1 1=1 1 1 1 =1 1 J 1 2 3 1 2 3
+4(L L +L L +L L1 )=2(L 1_L +L ) 2 1=1
1 2 2 3 3 2 3
6
y=i~1NY con

es decir:
3
yi~l (2L 1-1)Liyi+2L 1 L2 (y 1+y 2 )+2L 2 L (y 2 +y )+2L L1 (y 1+y ).
3 3 3 3

y= ( 2 L~- L 1j2 L1 L + 2 L1 L ) y + ( 2 L~ - L + 2 L L + 2 L L ) y +
1 3 3 1 2 2 1 3 2 2
2
+(2L -L +2L L +2L 1 L )y = L y +L y +L y =y
3 3 3 2 3 3 1 1 2 2 3 3
An~logamente se cumple:
x=!:N 1. X 1
2) Movimientos de deformaci6n constante
dU
e; = - =a
av au -av
e; = - =b y =-+-=e (a,b,c, constantes)
X dX Y ay iy ay ax

u=ax+w (y) v=by+lP, 2 (x)

1
lJi 1(y ) + lJi 2(x) =e , 1 u ego lJi 1(y ) =e o n s tan te ljJ 2(x) = -'Con s tan te
La expresi6n de los desplazamientos es:

u=a,x+f31y+y1
v=a.zx+Szy+y2
Los valores nodales son u
1 =a. 1 x_,+S 1 y 1+y 1 ,u 2 =a. 1 x 2 +B 1 y 2 +y 1 ,---
uc)=a.1x6+131y6+y1
u=~N
Sustituido en la expresi6n de 1 ,N
2 ,N 3 ,N 4 ,N 5 ,N 6 ) u
1
se obtiene
uz
u3
u4
u5
u6
G ~ 6 6 6
u=.~
1- 1 1 1 1- 1
N.u.=.~ N. (a x +s
1 1 2 1 2
y +y
1 )=a
1 .L:
1= 1 N.x.+s
1 1 1 .~ N.y.+y
1- 1 1 1 1 .I; N.=
1- 1 1

=a x+s y+y . como se quera demostrar

1 1 l
6
An~logamente se cumple
v=.r
1= 1
>N.v
1 1
.
 - 216 -

2). La continuidad C0 se satisface obviamente en

los nudos. A lo largo de un lado, las funciones de forma
son cuadrticas y estn definidos por tres valores. Luego
el elemento es conforme.
 - 217 -

Ejercicio E10.1.
El elemento finito no conforme de flexin (hipoele-
mento) de la figura E10.1a. tiene un campo de desplazamientos
que puede suponerse aproximado por una ley parablica de se--
gundo grado.

Fig. E 10 la

Se pide:
1) Determinar la matriz de rigidez~
2) Para una carga vertical uniforme en todo el ele
mento, hallar las cargas equlvalentes en los nu
dos.
 - 218 -

1) Se supone -la expresin de la flecha w parablica, es de-

cir, en coordenadas triangulares se ~iene:
2 2 2
w=B 1 L1 +a 2 L2 +a L +2S 1 L L +2S 2 L L +2S L L {ElO;la.}
3 3 2 3 1 3 3 1 2

Se formular~n las ecuaciones que siguen en un nudo.

En los restantes se obtendrn mediante permutacin circular
1-+2-+3-+1
Los coeficientes a 1 se determinan al imponer la CO,!l
dicin de flecha para el nudo i. Resulta para el nudo 1(1,0,0),
a 1 =w 1 Por lo tanto:

Para calcular los coeficientes B, es preciso conocer

antes la expresin de la derivada segn la normal a un lado.

b3 -1
;,g E 10 1b

Las frmulas de transformacin de coordenadas son:

(ver figura E1.0.lb.)

y anlogas
con a1=y2-y3
b 1 =x -x
3 2
c,=(x2y3-x3y2)
 - 219 -

x, yl
1
A=- x2 y2 = rea del tringulo 123
2
x3 y3

Entonces:
aw =(a aw + a2 aw + a3 aw )....!.
ax 1 "1:"1 '"'T., L 2A
3
'"

aw =(b aw aw aw ) 1
ay 1 '"'T, + b2 L2 + b3 L"'iA
3

. Se a .a. e 1 n g u 1 o q u e f o r ma e 1 v e e t o r 1 a do
3
IT en e 1
eje positivo de abscisas x. Su valor es:

b
cosa. =...1. con
3 .m3

los cosenos directores de la normal (dirigida hacia el interior

del tringulo) son:

, ~
m
3
Por lo tanto la derivada segn dicha normal es:
aw
-=- a3
~+ ~ ay
aw 1 a 1 a3 + bl b3
= 2A( Clw
an m
3 . 3
dX m3 m.3 ar,
2
a3. + b23
+ ~)
m3 ClL
3
es decir:
m
aw - ' 3 aw
an3- n(d3 rr, + e3

con
a3 al + b3 bl a3 a2 + b3 b2
d = e = y h3 L3 = 2A.
3 m2 3 m2
3 . 3
h3 es 1a altura de 1 tringulo que parte del nudo 3.
 - 220 -

Las derivadas de la flecha w respecto a las coorde

nadas triangulares son:
aw
aLl = 2 (al Ll + B2L3 + B3L2)

Al particularizar para el punto 6(1/2, 1/2, O) resulta:

-Por lo tanto, al considerar {E10.1b.} se obtiene:

Como d + e + 1 =,0 se puede_escribir

3 3

Resumlendo las ecuaciones an~logas en los nudos 4 y 5 se tiene:

-1 -1 s, -h
1 e nl o dl el w1
-1 -1 Bz = -h 2 e n2 + e2 o d2 w2
-1 -1 B3 -h e n3 d3 e o w3
3 3

Sistema cuya solucin es:

13 1 o 1 h 1 e nl o o d1 el w,

Bz
1
=-
2
o ti 2 e n2 --
1
.2
o e-2 o d2 w2

B3 o h3 e n3 o d3 e3 o w3
 - 221 -

o bien

s, 'o h 1 en 1 e2+d3 e3 d2 w,

82
1
=-
2
o h2 6 n2 --21 d3 dl+e3 e, w2

83 o h3 8 n3 e2. dl e 1+d 2 w3

Mul t ip'l icando por el vector f i 1a

(2L L , 2L L1 , 2L 1 L2 )
2 3 3

se obtienen las funciones de forma, al sumar~e al resultado

~ w1 L21 + w2 L22 + w L2 Es dec'tr
1 a expreston
3 3

Reduciendo se obtiene:

N1= L~ - e 2 L 2 (1-L~) - d 3 L3 (1-L 3 )= Nwl

N 4= L (1-e 1 )h 1 = N8 1
1 n1 _

Las matrices que relacionan deform~ciones

~=(hxx' hyy' hxy)T y los movimientos son:

T
B =(N N
- 1 1 ,xx 1 , yy

2A
2
Nl,yy=(b~ + e2 b~ + d
3
b~) = ~l

2A2 N1 ,xy={albl + e2 a2 b2 + d3 a3 ~3) = ""Z",

La matriz de rigidez se calcula a partir de la frmula:

T
k..=//A B. D B. dxdy
- J .-1 - -J
 - 222 -

DX D
1 o
con D= D1 Dy o matriz de elasticidad del material
o o Dxy

orttropo de 1 a losa.
Se obtiene por tanto:

k
11
= {-2
a D +
1 X za-1 b1 o, + -2 -2 1
bl Dy + el Dxy}3"
4A

k 12 =ra-, a2 DX + <a, b.,... + a2 r, >o, + 'b, ;-2 Dy + el c2

D }..J._
xy 4A3

T
i4=(N4 ,xx ,. N
4,yy N4,xy)

1 2
con N = - - hl al
4,xx 2 A2

1
N = - - hl b21
4,.yy 2A2

1
N =--
4,xy 2 A-2 hl al bl

Obtenindose:

hl - 2 b2 + a,2
K 4 =--{a a,
1 A3 1
oX + Ca1 1
r, >o 1 + b1 b2 D + el a 1 bl Dxy}
1 y
4

h2 . - 2 D + b2 + a22 'b,) o1 +
k
15
=---y{a 1 a2 X
Ca 1 2 r, b2 D +
2 y e 1 a2 b2 oxy }
4A
2
h1 4 2 2 2 4 2 b2
K.44= 3(a1 DX +(b21 al + al b 1 ) o1 + b 1' Dy +al 1 Dxy)
4A
hlh2 2 2 2 2 2
(al a2 Ox + (b~ a2 + al b 2 )o 1 + 1 2 oY+a1 bl a2 b2 Dxy)
K4' 5= 3 b2 b2
4A

y los demis elementos se hallan por permutacin circular que-

dando la matriz de rTgidez! de la forma siguiente:
1 k ww k wa]
K=- donde:
- 4A 3 [ kew kee
 a~Dx+2a 1 b 1 D 1 +b~Dy a~a~Dx+(a 1 b2 +i 2 'b 1 )o 1 a 1a Dx+(a1'b +a 'b 1 )o 1
3 3 3
+e2o
l_xy +'b1'b2oy+ctc2oxy O +~~-O
+b- 1b ~--y . 1 ~~y ________ _

a22 ox+2a- 2 b 2 D1+ a2a3Dx+(a2b3+a3b2)01

k
ww -2 -2
+b D +e o
2 Y 2 xy +'b2 'b3 oY+c~'C3 oxy
.sfJ....,
"''etr. -2 -
'co a 0x+2a b D +
3 3 3 1
CV'\

"'
N +b-2 o +e
-2
o
3 y 3 xy

-h (- (~ 2 2- ) - h 3 Ca 1 a~ D (a 1p ~+a ~b 1)'o 1
2 (_,; 2 2- } 0 +
a a D~+ a b +a b h (- 2
1 1 1 1 1 1 1 1 - 2 a 1a 2 o.x + al b 2+a 1 b 1 O1
-. 2 -
~blblDy+~lalblDxy) +b- 1 b 22 Dy+c- 1 a 2 b 2 Dxy ) +b 1 b;DY+c 1 a 3 b~Oxy)
(- 2 - 2 2- . :... ~ . - . 2 2....4 ' (- 2 - 2 2- )D
-hl ala1Dx+(a2bl+alb2)01 -h2(a2a2Dx+(a2b2+a2b2).01 -h3 ~2a30x+a2b3+a3b2 1
k
w8
= 2
- 2
+b 2 b 2 Dy+c- 2 a 2 b 2 Dxy
1
- 2 - +b b -
+b 2 b 1 Dy+c 2 a 1 b 1 Dxy) 2 3Dy+c 2 a 3 b 3 Dxy)

<-
- h 1 a a 2O
3 1
.xt r-a 3 b z1+a 2- )
1 b 3 O1
(- 2 - 2 2- )
-h2 a3a2Dx+a3b2+a2b3 01
- 2 - 2 2- )
-h3(a3a3Dx+a3~3+a3b3 1

- 2 -
- b 2 DY+c- a b Dxy )
+b 3 1 3 1 1 +b- b 22 Dy+c- 3 a 2 b 2 Dxy) +b b 3 Dy+c
3 a b Dxy
3 3 3
3
 k = kT
ew. we

2 4 2 2 2 2 2 2 2 2 2 2
h 1 (a 1 Dx+2a 1b 1 o.1 + hlh2 ( ala2Dx+bla201+ hlh3 ( ala3Dx+bla301+

2 2 2 2 2 2 2 2
4 2 2 +a 1 b 2 D1 +b 1 b 2 DY+ +a b o +b 1 b Dy+
1 3 1 3
+b 1 Dy+a 1 b 1 oxy)

+al'bla2b2Dxy) +alb1a3b3Dxy)
~e e=
2 4 2 2 2 2 2
h 2 a 2 Dx+2h 2 a 2 b 2 o1 h2h3(a2a3Dx+ N
N
.,:-

. 2 4 2 2 2 2 2 2 2
+h 2 b 2 DY+h 2 a 2 b 2 Dxy +b2a3Dl+a2b301+
-t;.,
?ti'
........
"e-() +b 2 b 3 Dy+a 2 b 2 a 3 b 3 Dxy )
2 2

2 4 2 2 2
h a Dx+2h a b o 1+
3 3 3 3 3
2 4 2 2 2
+h b Dy+h 3 a 3 b 3 Pxy
3 3
 - 225 -

2) Las cargas equivalentes a una carga vertical uniforme en

todo el elemento es:

Se comprueba que