CALCULO IV Examen Resuelto
CALCULO IV Examen Resuelto
CALCULO IV Examen Resuelto
FACULTAD DE INGENIERIA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERIA CIVIL
CALCULO 4
TITULO:
EXAMEN T2 (METODO DE ELEMENTOS FINITOS)
DOCENTE
INTEGRANTES
: CHACHA TANTA, Julio Milton ABANTO CABELLOS, Oswaldo David CAMPOS OCAS, Roco Maribel SANGAY CUSQUISIBAN, Milton Oswaldo AGUILAR OLIVARES, Karen Liseth
Datos: k1 = 200 kN/m k2 = 180 kN/m k3 = 210 kN/m k4 = 160 kN/m k5 = 220 kN/m k6 = 250 kN/m k7 = 150 kN/m P = 75 kN SOLUCIN
k2
Utilizamos los seis pasos para resolver este problema usando el elemento de resorte. 1.1. Discretizando el dominio El dominio se subdivide en siete elementos y seis nodos. La Tabla 1 muestra el elemento de conexin para este ejemplo. Tabla 1. Conectividad de elementos. N Elemento 1 2 3 4 5 6 7 Nodo i 1 3 3 4 4 6 5 Nodo j 3 5 4 5 6 5 2
1.2. Escribimos las matrices de rigidez del elemento Las siete matrices de rigidez del elemento k1, k2, k3, k4, k5, k6 y k7 se obtienen al cargar la funcin SpringElementStiffness en MATLAB. Cada matriz tiene un tamao de 2 2.
1.3. Ensamblaje de la matriz de rigidez global Dado que el sistema de resorte tiene seis nodos, el tamao de la matriz de rigidez global es 6 6. Por lo tanto, para obtener K se configura por primera vez una matriz cero de tamao 6 6 luego hacer siete llamadas a la funcin de MATLAB SpringAssemble ya que tenemos siete elementos elsticos en el sistema. Cada llamada a la funcin reunir a un elemento.
U1 U2 U3 U4 U5 U6
F1 F2 F3 F4 F5 F6
Las condiciones de contorno para este problema se dan como: U2 = 0, F1 = -75 kN, F3 = 0, F4 = 0, F5 = 0, F6 = 0 Insertando las condiciones anteriores en la matriz global se obtiene:
U1 0 U3 U4 U5 U6 -75 F2 0 0 0 0
1.5. Resolviendo las ecuaciones Resolviendo el sistema de ecuaciones se llevar a cabo mediante la particin (manualmente) y la eliminacin de Gauss (with MATLAB). En primer lugar generamos la submatriz de las filas 1, 3, 4, 5 y 6 y las columnas 1, 3, 4, 5 y 6. Por lo tanto se obtiene:
U1 U3 U4 U5 U6 -75 0 0 0 0
Se obtiene la solucin del sistema anterior usando MATLAB de la siguiente manera. Teniendo en cuenta que el operador barra invertida "\" se utiliza para la eliminacin de Gauss.
1.6. Post proceso En este paso, se obtienen las reacciones en los nodos y la fuerza en cada muelle, usando MATLAB de la siguiente manera. Primero creamos el vector de desplazamiento U, a continuacin, se calcula el vector de fuerza F.
Finalmente creamos el elemento vector de desplazamiento nodales U1, U2, U3, U4, U5, U6 y U7, a continuacin, se calcula el elemento vector de fuerzas f1, f2, f3, f4, f5, f6, f7; cargando la funcin en MATLAB SpringElementForces.
2.
E = 70 GPa A = 200 mm2 = 0.0002 m2 k = 2000 kN/m P = 8 kN SOLUCIN 2.1. Discretizando el dominio Tabla 2. Conectividad de elementos. N Elemento 1 2 3 Nodo i 1 2 3 Nodo j 2 3 4
2.2. Escribimos las matrices de rigidez del elemento Las siete matrices de rigidez del elemento k1, k2 y k3 se obtienen al cargar la funcin LinearBarElementStiffness y SpringElementStiffness en MATLAB. Cada matriz tiene un tamao de 2 2.
Las condiciones de contorno para este problema se dan como: U1= 0, U4= 0, F2 = 8 kN, F3 = 0 Insertando las condiciones anteriores en la matriz global se obtiene:
0 U2 U3 0
F1 8 0 F4
U2 U3
8 0
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: