FACULTAD DE INGENIERÍA

CARRERA DE INGENIERÍA CIVIL

CURSO:

METODOS NUMERICOS EN INGENIERIA

TEMA:
DESPLAZAMIENTO EN UNA BARRA CONICA
UTILIZANDO EL METODO DEMIEMBROS CARGADOS
AXIALMENTE Y ELEMENTO DE BARRA LINEAL

DOCENTE:

Ing. HUAMAN ROJAS, HEVER

INTEGRANTES:
 CAMPOS LEZAMA, DIEGO
DIAZ DIAZ, DENIN RIJKAARD
MALCA HERNANDEZ, JHAN DEYBIS
RAMOS BECERRA, KEVIN
VASQUEZ TONGO, DENYS IVAN
Cajamarca, 08 de Junio de 2017
 METODOS NUMERICOS EN INGENIERIA

RESUMEN

En el presente trabajo se hablara de como una fuerza puntual influye en el

cambio de la longitud de un cono, en el cual se pretende hallar dicho cambio de
longitud utilizando las fórmulas de resistencia de materiales concerniente al tema
(cambios en la longitud de barras no uniformes), este cono además será
analizado en el programa Matlab en donde se colocaran nueve nodos para
determinar su variación de longitud en cada uno de ellos , los resultados
obtenidos serán comparados por ambos métodos con la finalidad de poder
contrastar la variación de longitud que experimenta dicho elemento al ser
sometido a una fuerza.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

I. INTRODUCCIÓN

En ingeniería se necesita saber cómo responden los materiales sólidos a fuerzas

externas como la tensión, la compresión, la torsión, la flexión o la cizalladura. Los
materiales sólidos responden a dichas fuerzas con una deformación elástica (en
la que el material vuelve a su tamaño y forma originales cuando se elimina la
fuerza externa), una deformación permanente o una fractura.
es por ello que interesados por en conocer cómo es que varia la deformación o
el cambio de longitud en un cuerpo solido , el grupo de trabajo a tenido a bien
proponer un ejemplo en el cual el objeto a analizar es un tronco conico al cual se
le aplicara una determinada carga con el objetivo de observar el cambio de
longitud que experimenta este , dichos cálculos a realizar se harán a través del
método Miembros Cargados Axialmente , en el cual se realiza toda una
secuencia de análisis a fin de conocer cual es la deformación para un tronco
cónico , este tema es referido al curso de resistencia de materiales y finalmente
se hará uso del método elemento de barra lineal mediante el programa de
Matlab , lo que se pretende con el trabajo de investigación es calcular los valores
por ambos métodos , para luego compararlos y así de este modo observar si
existe similitud en el cambio de longitud entre cada método.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

II. KEYWORDS

 Deformación Axial. Es el cambio de longitud que sufre un cuerpo debido

a una carga axial, también se le conoce como deformación Total, y se
representa por δ (delta).

 Carga . Fuerza que actúa a lo largo del eje longitudinal de un miembro

estructural aplicada al centroide de la sección transversal del mismo
produciendo un esfuerzo uniforme. También llamada fuerza axial.

 Sección trasversal . Proyección de una sección de un objeto que se ha

realizado mediante un corte perpendicular al eje largo del mismo.

 Barra Ahusada. Barra que tiene un área transversal continuamente

variable.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

III. OBJETIVOS

1. OBJETIVO GENERAL

 Calcular el desplazamiento de la barra cónica al someterla a una fuerza

utilizando el método de Elemento de Barra Lineal y Miembros
Cargados Axialmente.

2. OBJETIVOS ESPECIFICOS

 Determinar los valores de los desplazamientos obtenidos usando el

método de elemento de barra lineal

 Determinar los valores de los desplazamientos obtenidos usando el

método de miembros cargados axialmente.

 Comparar los resultados obtenidos por ambos métodos para los

desplazamientos de la barra cónica sometidas a una misma fuerza.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

IV. MARCO TEORICO

3. ELEMENTO DE LA BARRA LINEAL

A. ECUACIONES BÁSICAS

El elemento de barra lineal es un elemento finito unidimensional en el que los

elementos locales y globales Coordenadas coinciden. Se caracteriza por
funciones de forma lineal y es idéntica Al elemento de resorte excepto que la
rigidez de la barra no se da directamente. los Lineal tiene módulo de elasticidad E,
área de sección A y longitud L. Cada elemento de barra lineal tiene nodos. En
este caso, el elemento Matriz de rigidez está dada por .

K=
[−EA/
EA/ L
L
−EA/ L
EA/ L ]

FIGURA 3.1: Elemento de la barra lineal

Obviamente, la matriz de rigidez del elemento para el elemento de barra lineal es

similar a la del elemento de muelle con la rigidez reemplazada por EA / L. Es
evidente que la Bar tiene sólo dos grados de libertad - uno en cada nodo. Por
consiguiente Para una estructura con n nodos, la matriz de rigidez global K será
de tamaño n × n (Ya que tenemos un grado de libertad en cada nodo). La matriz
de rigidez global K es ensamblado haciendo llamadas a la función MATLAB
LinearBarAssemble que es Escrito específicamente para este propósito.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

B. FUNCIONES DE MATLAB UTILIZADAS

Las cuatro funciones de MATLAB utilizadas para el elemento de barra lineal son:

 LinearBarElementStiffness (E, A, L) . Esta función calcula la rigidez del

elemento Matriz para cada barra lineal con módulo de elasticidad E, área
de sección A y Longitud L. Devuelve la matriz de rigidez de 2 × 2
elementos k .:

 LinearBarAssemble(K, k, i, j) . Esta función ensambla la matriz de rigidez

de elementos k de la barra lineal que une los nodos i (en el extremo
izquierdo) y j (en el extremo derecho) a la matriz de rigidez global K.
Devuelve la matriz de rigidez global n × n K cada vez que un elemento Está
montado.

 LinearBarElementForces (k, u). Esta función calcula el vector de fuerza

del elemento Utilizando la matriz de rigidez del elemento k y el vector de
desplazamiento del elemento u. Se devuelve el vector de fuerza de 2 × 1
elementos f.

 LinearBarElementStresses (k, u, A) - Esta función calcula el elemento de

tensión Vector utilizando la matriz de rigidez de elementos k, el vector de
desplazamiento de elementos u y El área de la sección transversal A.
Devuelve el vector de esfuerzo del elemento 2 × 1 sigma o s.

La siguiente es una lista del código fuente de MATLAB para cada función:

Esta función devuelve el elemento Matriz de rigidez para una barra lineal con
Módulo de elasticidad E, sección transversal Área A y longitud L. El tamaño
de la matriz de rigidez de los elementos es de 2 x 2.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Esta función ensambla la rigidez del elemento Matriz k de la barra lineal con
nodos i y j En la matriz de rigidez global K. Esta función devuelve la rigidez
global % Matriz K después de la matriz de rigidez del elemento K se ensambla.

Esta función devuelve el elemento nodal del vector de fuerza dada la rigidez del
elemento Matriz k y el elemento desplazamiento nodal Vector u.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Esta función devuelve el elemento nodal Del vector de fuerza dada la rigidez del
elemento Matriz k y el elemento desplazamiento nodal Vector u.

Esta función devuelve el elemento nodal del vector de esfuerzo dada la rigidez
del elemento Matriz k, el elemento desplazamiento nodal de vectores u, y el área
de sección transversal A.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

4. MIEMBROS CARGADOS AXIALMENTE

A. BARRAS CON CARGA O DIMENCIONES EN VARIACION CONTINUA

Algunas veces, la fuerza axial N y el área de la sección transversal varían

continuamente a lo largo del eje de la barra, como se ilustra por la barra ahusada
en la figura 2-11 a. esta barra tiene no solo un área transversal continuamente
variable , sino también una fuerza axial continuamente variable. En esta
ilustración, la carga tiene dos partes: una fuerza sencilla Pb que actúa sobre el
extremo B de la barra y las fuerzas distribuidas p(x ) que actúan a lo largo del eje.
(La fuerza distribuida tiene unidades de fuerza por unida de distancia, como libras
por pulgada o Newtons por metro). Una carga axial distribuida puede deberse a la
a fuerzas centrifugas, fuerzas de fricción o del peso de una barra en posición
vertical.
En estas condiciones, no podemos usar ya la ecuación para obtener el cambio de
longitud. Para ello debemos determinar el cambio de longitud de un elemento
deferencial de la barra y luego integrarlo la longitud de esta.
Seleccionamos un elemento diferencial a la distancia X desde el extremo
izquierdo de la barra (fig.2-11a). La fuerza axial interna N ( x ) que actúa en esta
sección (fig.2-11b) que puede determinarse por equilibrio usando el segmento AC
o el segmento CB como cuerpo libre. En general, esta fuerza es función de X.
Además, si conocemos las dimensiones de la barra, podemos expresar el área
transversal A(x) como función de X.

El alargamiento dδ del elemento diferencial (Fig. 2-11c) puede obtenerse con la

ecuación δ=PL/ EA sustituyendo N ( x ) por P, dx por L y A ( x ) por A, como sigue:

N ( x ) dx
dδ=
EA(x )

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

El alargamiento de la barra entera se obtiene integrando sobre la longitud:

2 2
N ( x ) dx
δ=∫ dδ=∫
0 0 EA(x )

Si las expresiones para N ( x ) y A ( x ) no son muy complicadas, la integral se

puede avaluar en forma analítica y obtener una fórmula para δ, como se ilustra
después del ejemplo 2-4. Sin embargo, si la integración formal es difícil o
imposible, deberá usarse un método numérico para evaluar la integral.
Una barra AB ligeramente ahusada de la sección circular sólida y longitud L
(Fig.2-13ª) esta fija en su extremo B y sometida a una carga P de tracción en su
extremo libre A. Los diámetros de la barra en los extremos A y B son de d A y d B ,
respectivamente.
Determinar el alargamiento de la barra debido a la carga P.

SOLUCION
La barra analizada en este ejemplo tiene una fuerza axial constante (igual a la
carga P) en toda su longitud, pero su área transversal varia continuamente de un
extremo a otro; por tanto debemos usar la integración (Ec.2-6) para determinar el
cambio de longitud.
Área transversal. El primer paso en la solución es obtener una expresión para el
área A(x ) transversal en cualquier sección transversal de la barra .Para esto,
debemos fijar un origen para la coordenada x. Una posibilidad es situar el origen
de coordenadas en el extremo libre A de la barra. Sin embargo, la integración por
efectuar se simplifica un tanto si localizamos el origen de coordenadas
extendiendo los lados de la barra ahusada hasta que se encuentren en el punto 0
como se muestra en la figura 2-13 b

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Las respectivas distancias L A y LBdel origen O a los extremos A y B están en la

proporción.
LA d A
= (a)
LB d B

Relación que se obtiene por triángulos semejantes en la figura 2-13b.A partir de

triángulos semejantes también podemos obtener la razón del diámetro d (x ) a la
distancia x desde el origen, el diámetro d A en el extremo pequeño de la barra:

d (x) x dA x
= O d ( x )= (b)
dA LA LA

Por tanto, el área transversal a la distancia x desde el origen es:

2
π [ d ( x )] π d A X
2 2

= 2
(c)
4 4 LA

Cambio de longitud sustituimos la expresión para A(x) en la ecuación (2-/)

Y obtenemos el alargamiento δ.
LB LB
N ( x ) dx Pdx(4 L2A ) 4 PL L2A dx
δ=∫ =∫ = ∫ (d)¿
EA(x ) L E (πd ¿ ¿ A 2 X 2) πEd 2A L x 2
A A

Al integrar (consulte las formular de integración en el apéndice C) y sustituir los

límites, obtenemos:

[ ]
2 2

( )
4 PL A −1 L B 4 PLA 1 1
δ=∫ 2
= 2
− (e)
πEd A x L A πEd A L A LB

Esta expresión para δ es entonces:

δ= ( )
4 PL L A
πEd 2A L B
(f )

Por último, sustituimos L A /L B=d A /d B (vea la Ec.a) y obtenemos:

4 PL
δ= (2.8)
πE d A d E

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

EJEMPLO DE APLICACION

Consideremos la barra cónica mostrada en la Fig. 3.3 con E = 210GPa y P =

18kN.Las áreas de sección transversal de la barra en los extremos izquierdo y
derecho son 0.002m2 y 0.012m2, respectivamente. Utilice elementos de barra
lineal para determinar el desplazamiento en el extremo libre de la barra cónica.

Barra Cónica

Una vez planteado el problema se procederá a realizar la solución del mismo por
dos métodos totalmente diferentes pero que en si guardan una relación muy
aparecida, ya que ambos nos proporcionaran una respuesta precisa, al cambio de
longitud que va a experimentar la barra cónica sometida a una fuerza P, como ya
se dijera en anteriores líneas el primer método consiste en desarrollar el problema
a través de fórmulas que tienen que ver con el curso de resistencia de materiales
mediante Miembros Cargados Axialmente, en el cual se determinara el cambio de
longitud total que experimenta dicho elemento y por último está el desarrollo del
problema usando para ello el programa Matlab mediante el tema Elemento de
Barra Lineal en el cual se determinara el cambio de longitud en nodos y puntos
diferentes . Los resultados obtenidos serán analizados para poder ver la similitud
que podemos encontrar utilizando estos dos métodos.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

METODO 1: MIEMBROS CARGADOS AXIALMENTE

Para realizar el cálculo de la deformación , se trabajara con diferentes distancias

al igual que el método de elemento de barra lineal , y para el cálculo del área en
diferentes puntos se usara la formula que se indica a continuación

1 DATOS

DESCRIPCION VALOR UNIDAD

Modulo de Elasticidad(E) 210 GPa
Fuerza (P) 18 KN
Área Mayor 0.012 m2
Área Menor 0.002 m2
Diámetro Mayor (db) 0.124 m
Diámetro Menor (da) 0.050 m
Longitud Barra 3 m

2. SOLUCION

formula para encontrar el área a una distancia X

 Distancias tomadas

DESCRIPCION VALOR UNIDAD

X1: 0.3 m
X2: 0.9 m
X3: 1.5 m
X4: 2.1 m
X5: 2.7 m

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

 Área a la distancia X tomada

VALOR UNIDAD
DIAMETRO
A1 0.003 m2
A2 0.005 m2
A3 0.007 m2
A4 0.009 m2
A5 0.011 m2

 Diámetro de cada área calculada

DIAMETRO VALOR UNIDAD

d1 0.062 m
d2 0.080 m
d3 0.094 m
d4 0.107 m
d5 0.118 m

 Calculo de la deformación según longitud tomada

Formula deducida en el marco teórico para cambiar la variación de la longitud

según las distancias tomadas

DEFORMACION VALOR UNIDAD

δ1 0.0043 m
δ2 0.0244 m
δ3 0.0344 m
δ4 0.0424 m
δ5 0.0493 m
δLT= 0.0525 m

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

METODO 2: ELEMENTO DE BARRA LINEAL

Paso 1 - Discretizar el Dominio:

Utilizaremos cinco elementos de barra lineales en este problema. Por lo tanto, la
barra cónica discretizado en cinco elementos y seis nodos como se muestra en la
figura . Para mayor precisión, los Resultados pueden utilizarse más elementos
(véase el problema). La Tabla muestra el elemento Conectividad para este
ejemplo. En este ejemplo, la longitud de cada elemento será 0,6 m. El área de la
sección transversal para cada elemento se tomará como la sección transversal En
el centro del elemento. Esto se puede calcular usando la siguiente fórmula Que da
el área de la sección transversal a cualquier distancia x desde el extremo
izquierdo.

Discretizacion del Ejemplo de la Barra cónica

Elemento de conectividad para el ejemplo

NUMERO DE NUD NUD

ELEMENTO Oi OJ
1 1 2
2 2 3
3 3 4
4 4 5
5 5 6

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Paso 2 - Escribir las Matrices de Rigidez del Elemento:

Las matrices de rigidez de cinco elementos k1, k2, k3, k4 y k5 se obtienen

haciendo llamadas a la función MATLAB LinearBarElementStiffness. Cada
matriz tiene tamaño 2 × 2.

E=210*10^6

E =

210000000

L=3/5

L =

0.6000

A1=0.002+(0.01*0.3/3)

A1 =

0.0030

A2=0.002+(0.01*0.9)/3

A2 =

0.0050

A3=0.002+(0.01*1.5/3)

A3 =

0.0070

A4=0.002+(0.01*2.1/3)

A4 =

0.0090

A5=0.002+(0.01*2.7/3)

A5 =

0.0110

k1=barras(E,A1,L)

k1 =

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

1050000 -1050000
-1050000 1050000

k2=barras(E,A2,L)

k2 =

1.0e+06 *

1.7500 -1.7500
-1.7500 1.7500

k3=barras(E,A3,L)

k3 =

2450000 -2450000
-2450000 2450000

k4=barras(E,A4,L)

k4 =

1.0e+06 *

3.1500 -3.1500
-3.1500 3.1500

k5=barras(E,A5,L)

k5 =

3850000 -3850000
-3850000 3850000

Paso 3 - Montaje de la Matriz Global de Rigidez:

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Dado que la estructura tiene seis nodos, el tamaño de la matriz de rigidez global
es 6 × 6. Por lo tanto, para obtener K primero establecemos una matriz cero de
tamaño 6 × 6 y luego hacer cinco llamadas A la función MATLAB LinearBarAssemble
ya que tenemos cinco elementos de barra lineal En la estructura. Cada llamada a
la función reunirá un elemento. El seguimiento Son los comandos MATLAB:

K=zeros(6,6)

K =

0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0

K=ensamble(K,k1,1,2)

K =

1050000 -1050000 0 0 0
0
-1050000 1050000 0 0 0
0
0 0 0 0 0
0
0 0 0 0 0
0
0 0 0 0 0
0
0 0 0 0 0
0

K=ensamble(K,k2,2,3)

K =

1.0e+06 *

1.0500 -1.0500 0 0 0 0
-1.0500 2.8000 -1.7500 0 0 0
0 -1.7500 1.7500 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0

K=ensamble(K,k3,3,4)

K =

1.0e+06 *

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

1.0500 -1.0500 0 0 0 0
-1.0500 2.8000 -1.7500 0 0 0
0 -1.7500 4.2000 -2.4500 0 0
0 0 -2.4500 2.4500 0 0
0 0 0 0 0 0
0 0 0 0 0 0

K=ensamble(K,k4,4,5)

K =

1.0e+06 *

1.0500 -1.0500 0 0 0 0
-1.0500 2.8000 -1.7500 0 0 0
0 -1.7500 4.2000 -2.4500 0 0
0 0 -2.4500 5.6000 -3.1500 0
0 0 0 -3.1500 3.1500 0
0 0 0 0 0 0

K=ensamble(K,k5,5,6)

K =

1.0e+06 *

1.0500 -1.0500 0 0 0 0
-1.0500 2.8000 -1.7500 0 0 0
0 -1.7500 4.2000 -2.4500 0 0
0 0 -2.4500 5.6000 -3.1500 0
0 0 0 -3.1500 7.0000 -3.8500
0 0 0 0 -3.8500 3.8500

Paso 4 - Aplicación de las Condiciones de Límite:

La matriz (3.2) para esta estructura se obtiene de la siguiente manera usando la

rigidez global Matriz obtenida en la etapa anterior:

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

Las condiciones de contorno para este problema se dan como:

F1 = −18, F2 = F3 = F4 = F5 = 0, U6 = 0

Paso 5 - Solución de las ecuaciones:

Resolviendo el sistema de ecuaciones en (3.11) se realizará particionando

(manualmente) Y la eliminación Gaussiana (con MATLAB). Primero particionamos

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

extrayendo la Submatrix en las filas 1 a 5 y columnas 1 a 5. Por lo tanto

obtenemos:

k=K(1:5,1:5)

k =

1.0e+06 *

1.0500 -1.0500 0 0 0
-1.0500 2.8000 -1.7500 0 0
0 -1.7500 4.2000 -2.4500 0
0 0 -2.4500 5.6000 -3.1500
0 0 0 -3.1500 7.0000

f=[-18 ; 0 ; 0 ; 0 ; 0 ]

f =

-18
0
0
0
0

u= k\f

u =

1.0e-04 *

-0.4517
-0.2802
-0.1774
-0.1039
-0.0468

V. ANALISIS DE RESULTADOS
Los resultados obtenidos muestran una similitud casi similar los cuales se
muestran a continuación en un cuadro

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

CUADRO COMPARATIVO
METODO
NODO
MIEMBROS CARCADOS AXIALMENTE ELEMENTO DE BARRA LINEAL

1 0.0043 0.0452
2 0.0100 0.0104
3 0.0140 0.0468
4 0.0173 0.0177
5 0.0201 0.0280

Cabe mencionar que los valores obtenidos en el matlab salen negativos porque
indica que la barra se desplaza hacia el lado izquierdo , asi como sus unidades
de medida al valor obtenido.

VI. CONCLUSIONES

 Se logro Calcular el desplazamiento de la barra cónica al someterla a

una fuerza utilizando el método de Elemento de Barra Lineal y
Miembros Cargados Axialmente.

 Se Logro determinar los valores de los desplazamientos obtenidos

usando el método de elemento de barra lineal

 Se logro determinar los valores de los desplazamientos obtenidos

usando el método de miembros cargados axialmente.

 Se compararon los resultados obtenidos por ambos métodos para los

desplazamientos de la barra cónica sometidas a una misma fuerza.

VII. BIBLIOGRAFIA

 Timoshenko.Resistencia de materiales 5ta° edición.

INGENIERIA CIVIL
 METODOS NUMERICOS EN INGENIERIA

 Petter Kattan.Guia de matlab.

 https://www.youtube.com/watch?v=gFf_YJxdM-c
 http://www.cartagena99.com/recursos/otros/apuntes/
Resistencia_de_Materiales_Aplicada.pdf

INGENIERIA CIVIL