Ejemplo 1.5 (Pgina 138).

THE FINITE ELEMENT METHOD IN ENGINEERING FOURTH EDITION Elsevier Science & Technology Books December 2004.

La viga dispone de una seccin rectangular con 1cm de ancho y 2cm de altura y la longitud es 60cm, sujeto a una carga vertical concentrado de 1000N como se muestra en la figura 1.12(a) (actuando e 20cm del nudo inicial). La viga est empotrada en ambos extremos. Obtener la tensin en la viga modelando mediante dos elementos. El mdulo de elasticidad del material es E=107 N/cm2.

Figura 1.12(a)

Solucin.
El fin principal del ejemplo, es mostrar una comparacin de los resultados que se obtiene de analizar la viga con sap2000 y los resultados de un proceso detallado mediante el mtodo de los elementos finitos, que para este caso est programado en matCAD. No se consideran los efectos por cortante en los elementos.

1. Para modelar la viga mediante elementos finitos se toma un sistema de de orientacin global, la misma que viene por defecto en Sap2000, para este caso se analizar en el plano X-Z como muestra la figura 1.12.1, que muestra la numeracin de los nudos y de los elementos, la conectividad entre elementos o la orientacin local para cada elemento.

Figura: 1.12.1: conectividad

Seguidamente, se ordena los argumentos

2. ARGUMENTOS
Para realizar el anlisis por el Mtodo de Elementos Finitos es necesario ordenar los argumentos de la siguiente manera

2.1 Nudos
Las coordenadas de los nudos considerados en la estructura columna 1: coordenada "x" del nudo en el sistema global columna 2: coordenada "y" del nudo en el sistema global NODE := 1 2 3 1 0 20 60 2 0 0 0

2.2 Propiedades de los elementos

Las propiedades de la seccin trasversal de los elementos y del material que est compuesto, cada fila representa una propiedad distinta y las columnas son: columna 1: momento de inercia columna 2: mdulo de elasticidad del material

PROP := 1

1 0.667

2 1107

2.3 Elementos
Cada fila representa una barra, contiene la informacin de la conectividad del elemento en el sistema, cada columna representa: columna 1: nudo inicial del elemento columna 2: nudo final del elemento columna 3: nmero de propiedad del elemento MEMB :=
1 2

1 1 2

2 2 3

3 1 1

2.4 Restricciones - Apoyos

Cada fila representa un apoyo de la estructura, las columnas informan el comportamiento para cada grado de libertad, la convencin es: "1" para los grados de libertad de desplazamiento restringido "0" para los grados de libertad donde existe desplazamiento libre. Cada columna representa: columna 1: nmero del nudo donde existe el apoyo columna 2: "uy?" informacin del desplazamiento en la direccin "y" global columna 3: "Rz?" iinformacin de la rotacin en el eje "z" global

SUPP :=
1 2

1 1 3

2 1 1

3 1 1

2.5 Cargas 2.5.1 cargas en nudos(NLF)

Se requiere las cargas en los nudos, si los hubiera, cada columna representa. Columna 1: Nmero del nudo Columna 2: carga puntual en la direccin global Y. Columna 3: momento concentrado en la direccin global z.

NLF :=
1

1 2

2 -1103

3 0

3. MATRIZ DE RIGIDEZ RESPECTO AL SISTEMA GLOBAL DE LOS

ELEMENTOS
Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

Obteniendo la matriz de rigidez respecto al sistema global para el elemento:

m := 1

3.1 Matriz de rigidez respecto al sistema local

Momento de inercia respecto al eje z

I( m) = 0.667 Mdulo de elasticidad del material

7

E( m) = 1 10

Longitud del elemento

(xf (m) xi( m))2 + (yf (m) yi( m))2

l ( m) = 20

( 20 0 ) + ( 0 0 ) = 20

reemplazando estos valores, resuta:

gdl ( m) = ( 1 2 3 4 ) 10 1 10 1 10 133.333 10 66.667 4 k e( m) = 10 1 10 10 1 10 66.667 10 133.333 ..... se procede de igual manera para cada elemento 1 2 gdl ( m) = 3 4

3.2 Matriz de rigidez total ensamblado

Seguidamente se ensambla en una sola matriz segn los grados de libertad, lo que resulta.

10 1 10 1 10 133.333 10 66.667 1.125 7.5 1 10 K= 10 66.667 7.5 200 0 0.125 2.5 0 2.5 33.333 0 0

0 0 0.125 2.5 0.125 2.5

0 2.5 4 10 33.333 2.5 66.667

0

4. VECTOR DE FUERZAS NODALES EQUIVALENTES

Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

4.2 Vector de cargas puntuales en nudos

Existe carga puntual en el nudo 2, en el grado de libertad 3, las cargas en los nudos se ensamblan directamente en el vector de fuerzas expresados en coordenadas globales(para vigas el sistema de corrdenada global y el sistema de coordenadas locales coinciden)

0 0 1 103 FNLF = 0 0 0
4.3 Finalmente, el vector de carga nodal equivalente resulta
F=F.nld + CARGAS NODALES EQUIVALENTES

0 0 1 103 FNLF + 0 = 0 0 0 0 0 1 103 F= 0 0 0

5 DESPLAZAMIENTOS EN NUDOS
Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

5.1 Imponiendo las condiciones de contorno

La matriz de rigidez total(K) de la estructura obtenido anteriormente est libre en todos sus grados de libertad(desplazamiento libre), pero como la estructura cuenta con apoyos en las cuales los desplazameintos son nulos, se modifica "K" para tal efecto.

4 1 10307 1 105 1 10 1 105 1 10307 1 105 1 104 1 105 1.125 104 Km = 1 105 6.667 105 7.5 104 3 0 0 1.25 10 4 2.5 10 0 0

5 0 0 6.667 10 4 3 4 2.5 10 7.5 10 1.25 10 6 4 5 2 10 2.5 10 3.333 10 4 307 4 2.5 10 1 10 2.5 10 5 4 307 3.333 10 2.5 10 1 10
1 10
5

5.2 Matriz aumentada

Teniendo los coeficientes de las incgnitas(Km) y los trminos independientes(F) de las ecuaciones simultneas que se forman para cada grado de libertad en los nudos, existe muchas maneras para resolver la ecuacin. Formando la matriz aumentada

4 1 10307 1 105 1 10 1 105 1 10307 1 105 1 104 1 105 1.125 104 augment( Km , F) = 1 105 6.667 105 7.5 104 3 0 0 1.25 10 4 2.5 10 0 0

5 0 0 0 6.667 10 4 3 4 3 2.5 10 1 10 7.5 10 1.25 10 6 4 5 2 10 2.5 10 3.333 10 0 4 307 4 2.5 10 1 10 2.5 10 0 5 4 307 3.333 10 2.5 10 1 10 0
1 10
5

5.3 Resolviendo la ecuacicn

la ecuacin anterior en su forma escalonada reducida

1 0 0 rref ( augment( Km , F) ) = 0 0 0
5.4 Los desplazamientos
D = 0 0 0.119 4.444 10
T

0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1

0 0.119 3 4.444 10 0 0
0

0 0

Ordenando

0 0 1 NDIS = 2 0.118519 0.004444 0 0 3

Donde: Columna 1: Nmero del nudo Columna 2: Reaccin e Y(Sistema Global) Columna 3: Giro en

Para comparar, se muestra en el siguiente cuadro los desplazamientos obtenidos con un anlisis realizado en sap2000 12.0.0 Educacional.

Joint Text 1 2 3

OutputCase Text DEAD DEAD DEAD

CaseType Text LinStatic LinStatic LinStatic

TABLE: Joint Displacements U1 U2 U3 R1 R2 R3 cm cm cm Radians Radians Radians 0 0 0 0 0 0 0 0 -0.118519 0 0.004444 0 0 0 0 0 0 0

se obserava que coinciden ambos resultados, entonces puede ser que sap2000 12.0.0 utilice estas mismas formas de analizar la estructura.

6 REACCIONES EN LOS APOYOS

La matriz de rigidez total de la estructura es.
4 1 104 1 105 1 10 1 105 1.333 106 1 105 1 104 1 105 1.125 104 K= 1 105 6.667 105 7.5 104 3 0 1.25 10 0 4 2.5 10 0 0

5 0 0 6.667 10 4 3 4 2.5 10 7.5 10 1.25 10 6 4 5 2 10 2.5 10 3.333 10 4 3 4 2.5 10 1.25 10 2.5 10 5 4 5 3.333 10 2.5 10 6.667 10
1 10
5

el vector de desplazamientos:
T

D = 0 0 0.119 4.444 10

0 0

 740.741 8.889 103 13 1.137 10 K D F = 0 259.259 3 4.444 10

740.741 8.889 103 13 1.137 10 R= 0 259.259 3 4.444 10

Ordenado :

nudo REAC =

Ry

Mz

1 740.740741 8888.888889 3 259.259259 4444.444444

Con la finalidad de comparar, el siguiente cuadro mustra los resultados de analizar con sap2000 12.0.0.

Joint Text 1 3

OutputCase Text DEAD DEAD

F1 N

TABLE: Joint Reactions F2 F3 M1 M2 N N N-cm N-cm 0 0 740.740741 0 -8888.8889 0 0 259.259259 0 4444.444444

M3 N-cm 0 0

Se observa que los resultados son los mismos.

7 ESFUERZOS EN LOS EXTREMOS DE LOS ELEMENTOS

Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

Obteniendo los resultados para el elemento

m := 2

7.1 Matriz de rigidez respecto al sistema local

3 4 2.5 10 1.25 10 4 5 2.5 10 6.667 10 k e( m) = 3 4 1.25 10 2.5 10 2.5 104 3.333 105

1.25 10 2.5 10 1.25 10 2.5 10

4 3 4

5 3.333 10 4 2.5 10 5 6.667 10

2.5 10
4

7.2 Desplazamientos en los nudos del elemento en el sistema local

Grados de libertad de los nudos

0.119 3 4.444 10 De( m) = 0 0

3 4 gdl ( m) = 5 6

7.3 Esfuerzos en los extremos del elemento

259.259 5.926 103 k e( m) De( m) = 259.259 4.444 103 3 4 gdl ( m) = 5 6

............... de igual manera para cada elemento Resumiendo para todos los elementos, se tiene:

1 1 MFOR = 2 3 4 1 1 2 2

2 1 2 2 3

3 740.740741 -740.740741 -259.259259 259.259259

4 8888.888889 5925.925926 -5925.925926 -4444.444444

Donde: columna 1: nmero de elemento columna 2: nmero de nudo columna 3: exfeurzo axial columna 4: esfuerzo cortante columna 5: momento flector Los resultados Obtenidos con sap20000 12.0.0 se muestra en el cuadro.

Frame Text 1 1 2 2 1 2 2 3

Joint Text

TABLE: Element Joint Forces - Frames F1 F2 F3 M1 M2 N N N N-cm N-cm 0 0 740.740741 0 -8888.8889 0 0 -740.740741 0 -5925.9259 0 0 -259.259259 0 5925.925926 0 0 259.259259 0 4444.444444

M3 N-cm 0 0 0 0

los resultados son los mismos

8. LEYES DE DESPLAZAMIENTOS
Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

Se obtendrn las ecuaciones que gobiernan los desplazamientos en todo el elemento, por esta razn todo est referido de acuerdo su sistema de orientacin local. Con la finalidad de obtener el grfico desde pares de puntos, se evaluar las ecuaciones para 11 puntos sobre el elemento. Analizando para el elemento: m := 2

8.1 Grfico del desplazamiento lineal - en el eje "y" local

Desplazamientos en los extremos de los elementos los desplazamientos en los extremos del elemento, para los grados de libertad asociados a los nudos del elemento son Desplazamiento del nudo inicial u 1 := De( m)
1

1 := De( m) 2

Desplazamiento del nudo final u 2 := De( m)

3

0.119 3 4.444 10 De( m) = 0 0

3 4 gdl ( m) = 5 6

2 := De( m) 4

Longitud del elemento L := l( m) *

Ecuacin del desplazamiento lineal. En Sap2000 12.0.0, para una distancia de 20cm desde el nudo inicial, el desplazamiento, resulta.

( x ) :=

2 x3 3 x2 + 1 u 1 + 3 2 L L

2 3 x 2 x + x + L 2 1 L

3 x2 2 x3 u + 2 3 2 L L

x3 x2 L2 L 2

Evaluando en 11 puntos. x := 0 , l( m) 10 .. l ( m)

x =
0 4 8 12 16 20 24 28 32 36 40

( x ) =
-0.118519 -0.1296 -0.128948 -0.119052 -0.1024 -0.081481 -0.058785 -0.0368 -0.018015 -0.004919 0
0.15 0.1 ( x) 0.05 0

Desplazamiento lineal
10 20 30 40

8.2 Desplazamiento angular. Conocido la ecuacin del desplazamiento lineal, el desplazamiento angular est dado por ( x ) := d ( x ) dx

Evaluando para 11 puntos.

x =
0 4 8 12 16 20 24 28 32 36 40

( x ) =
-0.004444 -0.0012 0.001422 0.003422 0.0048 0.005556 0.005689 0.0052 0.004089 0.002356 0
5 10
3

Desplazamiento angular
0.01

5 10 ( x)

10

20

30

40

9. LEYES DE TENSIONES Y ESFUERZOS

Reference:D:\FEM\Vigas2D\0 VIGA2D_FUNCIONES.xmcd

Se obtendrn las ecuaciones que gobiernan los tensiones y esfuerzos en todo el elemento, por esta razn todo est referido de acuerdo su sistema de orientacin local. Con la finalidad de obtener el grfico desde pares de puntos, se evaluar las ecuaciones para 11 puntos sobre el elemento. Analizando para el elemento: Resumiendo, el desplazamiento est dado por: 0.118519 0.004444 De( m) = 0 0 u 1 := De( m) u 2 := De( m) L := l( m) *
2 3 x 2 x + x + L 2 1 L 1 3

m := 2

1 := De( m) 2 2 := De( m) 4

( x ) :=

2 x3 3 x2 + 1 u 1 + 3 2 L L

3 x2 2 x3 u + 2 3 2 L L

x3 x2 L2 L 2

9.1 Tensin en la direccin local x' (Direccin local "1" en Sap2000 12.0.0)
La tensin para elementos viga(frame en Sap2000 12.0.0) est dado por. d
2 2

xx( x , y ) := y E( m)

( x )

dx

la tensin ser evauada en la fibra superior del elemento m = 2 , para lo cual y=1. Evaluando la funcin para obtener las tensiones mximas, que en este caso suceden cuando x=0(Nudo unicial) y cuando x=40(nudo final). . xx( 0 , 1 ) = 8888.888894

N/cm2

xx( 40 , 1 ) = 6666.666667

N/cm2

9.2 El momento flector(M33 en Sap2000 12.0.0)

El momento flector est establecido mediante. d
2 2

M ( x ) := E( m) I( m)

( x )

dx

Evaluando en 11 puntos. x := 0 , l( m) 10 .. l ( m)

x =
0 4 8 12 16 20 24 28 32 36 40

M( x) =
-5925.925929 -4888.888889 -3851.851852 -2814.814815 -1777.777778 -740.740741 296.296296 1333.333333 2370.37037 3407.407407 4444.444444
1 10
4

Momento Flector
5 10
3

0 M ( x) 5 10
3

10

20

30

40

Los resultados del Sap2000 12.0.0 se muestra en la siguiente cuadro..

Frame Text 2 2 2 2 2 2 2 2 2 2 2

Station cm

TABLE: Element Forces - Frames V2 M3 FrameElem ElemStation N N-cm Text cm 0 259.259259 5925.925926 2-1 0 4 259.259259 4888.888889 2-1 4 8 259.259259 3851.851852 2-1 8 12 259.259259 2814.814815 2-1 12 16 259.259259 1777.777778 2-1 16 20 259.259259 740.740741 2-1 20 24 259.259259 -296.296296 2-1 24 28 259.259259 -1333.33333 2-1 28 32 259.259259 -2370.37037 2-1 32 36 259.259259 -3407.4074 2-1 36 40 259.259259 -4444.4444 2-1 40

9.3 Fuerza cortante(V22 en Sap2000 12.0.0)

V( x ) := d M( x) dx

Evaluando para 11 puntos.

x =
0 4 8 12 16 20 24 28 32 36 40

V( x ) =
259.259259 259.259259 259.259259 259.259259 259.259259 259.259259 259.259259 259.259259 259.259259 259.259259 259.259259
250 0 10 260 V ( x) 270 280

Fuerza Cortante

20 x

30

40

Finalmente, el grfico generado por sap2000 de los esfuerzos para la viga modelada.

Conclusin.
Sap2000 v12.0.0 Educacional utiliza para el anlisis de objetos "Frame"(vigas - columnas) elementos con nudos en sus extremos, en este caso utiliza las mismas teoras que se encuentra en muchos libros que trata sobre FEM para elementos viga.

Comentarios Sugerencias