Ejemplo 1.8 (Pgina 248).

INTRODUCCIN A ESTUDIO DEL ELEMENTO FINITO EN INGENIERA SEGUNDA EDICIN Tirupathi R. Chandrupatla Ashok D. Belegundu Prentice Hall, Mxico, 1999 Para la viga 1 y carga mostrada en la figura E8.1, determine los pendientes en 2 y 3 y la deflexin vertical en le punto medio de la carga distribuida. E=200GPa, I=4*106 mm4

Figura E-8.1

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.

Seccin trasversal del elemento es: 48mm*100mm

1 Convenciones 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.

Figura: E-8.1.1 numeracin de nudos y barras

Seguidamente se establece el nudo inicial y final para cada elemento con la cual queda implcitamente establecida la orientacin local.

Figura: E-8.1.2 conectividad de cada elemento

Luego, 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 1 2 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 410-6

2 21011

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 libre "Valor del desplazamiento" para los grados de libertad restringidos. 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 3

1 1 2 3

2 0 0 0

3 0 1 1

2.5 Cargas
Cargas\Cargas en nudos Cargas\Carga Puntual en Elemento Cargas\Carga Repartida Trapezoidal

2.5.3 Carga Repartida Trapezoidal

MLY := 0

Cada columna representa: Col_1: nmero de elemento Col_2: carga Inicial. Col_3: carga Final Col_4: distancia inicial Col_5: distancia final

MLY :=
1

1 2

2 -1.2104

3 -1.2104

4 0

5 1

Cargas\Carga Repartida Trapezoidal

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

6

I( m) = 4 10

Mdulo de elasticidad del material

11

E( m) = 2 10

Longitud del elemento

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

l( m) = 1

( 1 0) + ( 0 0) = 1

reemplazando estos valores, resuta:

gdl ( m) = ( 1 2 3 4 )

12 6 k e( m) = 12 6

4 6 2 5 10 8 6 12 6 2 6 4
6

12 6

1 2 gdl ( m) = 3 4

..... se procede de igual manera para cada elemento

3.4 Matriz de rigidez total ensamblado

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

12 6 12 K= 6 0 0

6 4 2 0 0

12 6 6 0 6 2 0 8 2

0 0 12 6 6

6 24

12 6 12

6 5 10 8 2 6 4
0 0

4. VECTOR DE FUERZAS NODALES EQUIVALENTES

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

4.1 Vector de carga nodal equivalente(de las cargas puntuales en nudos)

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 0 NLFnudos = 0 0 0

No existe cargas directamente en los nudos para la viga analizada.

4.2 Vector de carga nodal equivalente (de cargas que actuan en el elemento(puntual,trapezoidal))
En este caso se obtiene el momento de empotramiento perfecto para cada carga en todos los elementos y se ensanbla en un vector con signo cambiado. Para la viga existe una carga distribuida constante en el elemento 2. Las funciones de interpolacin se utiliza tambien para reducir las cargas a los grados de libertad de los nudos. Nmero del elemento: Longitud del elemento: m := 2 l( m) = 1 d1 := 0

Distancia a la carga desde el nudo inicial: Carga a la distancia d1: q1 := 12000 N m

Distancia a la carga desde el nudo final: N m

d2 := 1

Carga a la distancia d2:

q2 := 12000

con los arqumentos, la ecuacin de la carga es: q ( x , q1 , q2 , d1 , d2) = 12000 constante a lo largo del elemento

y los momentops de empotramiento perfecto resulta

6 103 3 1000 = 1 10 MEP_MLY( m , x , q1 , q2 , d1 , d2) 6000 3 1000 6 10 3 1 10

6000

3 4 gdl ( m) = 5 6

Estos resultados se ensambla en el vector de momentos de empotramiento perfecto de todas las cargas que actan en los elementos.

 0 0 6 103 MEPbarras = 1 103 3 6 10 3 1 10

4.3 Finalmente, el vector de carga nodal equivalente resulta

F= -NLFnudos - MEPbarras)

0 0 0 NLFnudos = 0 0 0

0 0 6 103 MEPbarras = 1 103 3 6 10 3 1 10

0 0 6 103 F= 1 103 3 6 10 3 1 10

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.

1 10307 4.8 106 9.6 106 Km = 4.8 106 0 0

307 6 6 0 0 1 10 4.8 10 1.6 10 6 307 6 6 9.6 10 4.8 10 4.8 10 1 10 0 6 6 6 6 1.6 10 0 6.4 10 4.8 10 1.6 10 6 6 307 6 0 4.8 10 9.6 10 4.8 10 1 10 6 6 6 6 0 4.8 10 1.6 10 4.8 10 3.2 10
4.8 10
6

9.6 10

4.8 10

Modificando el vector de fuerzas, segn a los desplazamientos en los apoyos.

0 0 0 Fm = 3 1 10 0 3 1 10

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

1 10307 4.8 106 9.6 106 augment( Km , Fm) = 4.8 106 0 0

307 6 6 0 0 0 1 10 4.8 10 1.6 10 6 307 6 6 9.6 10 4.8 10 0 4.8 10 1 10 0 6 6 6 6 3 6.4 10 4.8 10 1.6 10 1 10 1.6 10 0 6 6 307 6 0 4.8 10 0 9.6 10 4.8 10 1 10 6 6 6 6 3 4.8 10 1.6 10 4.8 10 3.2 10 0 1 10
4.8 10
6

9.6 10

4.8 10

5.3 Resolviendo la ecuacicn

la ecuacin anterior en su forma escalonada reducida

1 0 0 rref ( augment( Km , Fm) ) = 0 0 0

5.4 Los desplazamientos
D = 0 0 0 2.679 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 4 2.679 10 0 4 4.464 10
0 0

0 4.464 10

Ordenando

0 1 0 4 NDIS = 2 0 2.679 10 4 3 0 4.464 10

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.

en forma de tabla

0 1 0 2 0 0.000268 NDIS = 3 0 0.000446

Joint Text 1 2 3

U1 m

TABLE: Joint Displacements U2 U3 R1 R2 R3 m m Radians Radians Radians 0 0 0 0 0 0 0 0 0 0 0.000268 0 0 0 0 0 -0.000446 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.

9.6 106 4.8 106 9.6 106 K= 4.8 106 0 0

6 6 6 0 0 3.2 10 4.8 10 1.6 10 6 7 6 6 9.6 10 4.8 10 4.8 10 1.92 10 0 6 6 6 6 1.6 10 0 6.4 10 4.8 10 1.6 10 6 6 6 6 0 9.6 10 4.8 10 9.6 10 4.8 10 6 6 6 6 4.8 10 1.6 10 4.8 10 3.2 10 0
4.8 10
6

9.6 10

4.8 10

el vector de desplazamientos:
T

D = 0 0 0 2.679 10

0 4.464 10

1285.714 428.571 8142.857 K D + MEPbarras + NLFnudos = 0 5142.857 0

Ordenado :

nudo

Ry

Mz

1 1285.714286 428.571429 REAC = 2 8142.857143 0 0 3 5142.857143

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

Joint Text 1 2 3

F1 N 0 0 0

TABLE: Joint Reactions F2 F3 M1 M2 N N N-m N-m 0 -1285.71429 0 428.571429 0 8142.857143 0 0 0 5142.857143 0 0

M3 N-m 0 0 0

Se observa que los resultados son los mismos.

6 ESFUERZOS EN LOS EXTREMOS DE LOS ELEMENTOS

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

Obteniendo los resultados para el elemento

m := 2

6.1 Matriz de rigidez respecto al sistema local

12 6 k e( m) = 12 6

4 6 2 5 8 10 6 12 6 2 6 4
6

12 6

6.2 Desplazamientos en los nudos del elemento en el sistema local

0 2.679 10 4 De( m) = 0 4.464 10 4 Grados de libertad de los nudos

3 4 gdl ( m) = 5 6

6.4 Esfuerzos en los extremos del elemento

857.143 142.857 k e( m) De( m) = 857.143 3 1 10

3 4 gdl ( m) = 5 6

6.5 Los momentos de empotramiento perfecto para el elemeto "m"

Los esfuerzos obtenidos en 6.4 debe ser corregido por lo momentos de empotramiento perfecto. sumando al anterior, resulta.

6 103 3 1 10 MEP( 2 ) = 3 6 10 1 103

6.857 103 857.143 k e( m) De( m) + MEP( m) = 5.143 103 0

............... 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 -1285.714286 1285.714286 6857.142857 5142.857143

4 -428.571429 -857.142857 857.142857 0

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-m N-m 0 0 -1285.71429 0 428.571429 0 0 1285.714286 0 857.142857 0 0 6857.142857 0 -857.142857 0 0 5142.857143 0 -2.274E-13

M3 N-m 0 0 0 0

los resultados son los mismos

7. LEYES DE ESFUERZOS Y DESPLAZAMIENTOS

Se obtendrn las ecuaciones que gobiernan los desplazamientos en todo el elemento, por esta razn todo est referido de acuerdo a su sistema de orientacin local. Con la finalidad de obtener el grfico desde pares de puntos, se evaluar las ecuaciones para "n" puntos sobre el elemento. Cada elemento se trata como a una estructura completa: con nudos en sus extremos que estn sometidos a desplazamientos iniciales, que son los desplazamientos que se obtuvieron en el proceso anterior; Con cargas puntuales actuando en sus apoyos, que son los esfuerzos en los extremos del elemento que se obtuvieron en el proceso anterior, y adems, todas las cargas que acta sobre el elemento. Analizando para el elemento: m=2 Se dividir el elemento 10 segmentos que hacen 11 nudos. 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 4 5 6 7 8 9 10 11

1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

2 0 0 0 0 0 0 0 0 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 410-6

2 21011

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 3 4 5 6 7 8 9 10

1 1 2 3 4 5 6 7 8 9 10

2 2 3 4 5 6 7 8 9 10 11

3 1 1 1 1 1 1 1 1 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 libre "Valor del desplazamiento" para los grados de libertad restringidos. 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 11

2 0 0

3 -2.67910-4 4.46410-4

2.5 Cargas
Cargas\Cargas en nudos

2.5.1 cargas en nudos(NLF)

NLF := 0

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 2
Cargas\Cargas en nudos Cargas\Carga Puntual en Elemento Cargas\Carga Repartida Trapezoidal

1 11

6.857103 5.143103

857.143 0

2.5.3 Carga Repartida Trapezoidal

MLY := 0

Cada columna representa: Col_1: nmero de elemento Col_2: carga Inicial. Col_3: carga Final Col_4: distancia inicial Col_5: distancia final

MLY :=
1 2 3 4 5 6 7 8 9 10

1 1 2 3 4 5 6 7 8 9 10

2 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104

3 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104 -1.2104

4 0 0 0 0 0 0 0 0 0 0

5 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1

Cargas\Carga Repartida Trapezoidal

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

7.1 leyes de desplazamientos

nudo

desplazamiento 2 1 2 3 4 5 6 7 8 9 10 11 0 -0.000031 -0.000065 -0.000095 -0.000117 -0.000128 -0.000126 -0.00011 -0.000082 -0.000044 0

Rotacin 3 -0.000268 -0.000335 -0.000331 -0.000271 -0.000171 -0.000045 0.000092 0.000225 0.000338 0.000417 0.000446

Puntos en que se evaluan los resultados a lo largo del elemento x := li 0 for i 1 .. 11 x

i, 1

1 1 2 3 4 NDIS = 5 6 7 8 9 10 11

li
1 11 1

li x 1 1 2 3 4 x= 5 6 7 8 9 10 11 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

+ li

DESPLAZAMIENTO LINEAL
0 0.2 0.4 0.6 0.8 1

Desplazamiento(m)

5 10 5 2 NDIS 1 10 4

1.5 10 4 x

Distancia(m)

6 10

DESPLAZAMIENTO ANGULAR

4 10

Rotacin(rad)

3 NDIS

2 10

0.2

0.4

0.6

0.8

2 10 4 4 10 4 x

distancia(m)

Comparando los resultados con los de SAP2000 12.0.0, se obtiene los mismos resultados.

7.2 Leyes de esfuerzos

Elemento
1 5 6 7 8 9 10 11 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10

Nudo
2 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11

Cortante
3 4456.8 -3256.8 3256.8 -2056.8 2056.8 -856.8 856.8 343.2 -343.2 1543.2 -1543.2 2743.2 -2743.2 3943.2 -3943.2 5143.2

Momento
4 -274.4 660.08 -660.08 925.76 -925.76 1071.44 -1071.44 1097.12 -1097.12 1002.8 -1002.8 788.48 -788.48 454.16 -454.16 ...

MFOR = 12
13 14 15 16 17 18 19 20

Esfuerzos internos ordenados

MFORo :=

a1 for i 1 , 3 .. rows( MFOR) MFORo MFORo MFORo

a, 1 a, 2 a, 3

a i, 3 i, 4

MFOR MFOR

aa+1 MFORo MFORo MFORo MFORo

rows( x) , 1 rows( x) , 2 rows( x) , 3

rows( x) i+ 1 , 3 i+ 1 , 4

MFOR MFOR

Comparando resultados V2 M3
3 -856.96 -231.28 274.4 660.08 925.76 1071.44 1097.12 1002.8 788.48 454.16 -0.16

2 2 2 2 2 2 2 2 2 2 2

TABLE: Element Forces - Frames Frame Station V2 M3 Text m N N-m 0 -6857.14 -857.14 0.1 -5657.14 -231.43 0.2 -4457.14 274.29 0.3 -3257.14 660 0.4 -2057.14 925.71 0.5 -857.14 1071.43 0.6 342.86 1097.14 0.7 1542.86 1002.86 0.8 2742.86 788.57 0.9 3942.86 454.29 1 5142.86 -7.745E-13

1 1 2 3 4 5 6 7 8 9 10 11 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

2 -6856.8 -5656.8 -4456.8 -3256.8 -2056.8 -856.8 343.2 1543.2 2743.2 3943.2 5143.2

MFORo =

ESFUERZO CORTANTE
1 10 4

5 10

Esfuerzo cortante(N)

2 MFORo

0.2

0.4

0.6

0.8

5 103

1 104 x

distancia(m)

Momento flector
1 10 3

Momento flector(N-m)

500

3 ( MFORo) 500

0.2

0.4

0.6

0.8

1 103

1.5 103 x

Distancia(m)

Conclusiones:
El mtodo de los elementos finitos, es ideal para analizar estructuras mediante programacin, esta vez se hizo detalladamente, con la finalidad de ver el proceso "paso a paso". Inicialmente se calcula la estructura tomando a todos los elementos con nudos en sus extremos, luego para obtener los grficos de desplazamientos y esfuerzos(las ecuaciones) se vuelve a subdividir a cada elemento, tomando como argumentos los resultados de la primera etapa y se procede de la misma manera. Los grficos no son mas que unin de pares de puntos mediante lneas de los resultados evaluados sobre el elemento, para el ejemplo se subdividi en 10 elemento(11 nudos). De esta opcin dispone sap2000 12.0.0, el Etabs 9.5, FEM4849 v5.3 sr6 con su mdulo "QUERY". Teniendo los valores para una serie de puntos sobre el elemento, una forma de obtener la ecuacin, es buscar su lnea de tendencia mediante regresin, el mejor ajuste!. Autor: Edmundo Canchari Gutirrez Comentarios y sugerencias: cgedmundo@gmail.com Visite: http://cgedmundo.wordpress,com