Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 26 vistas

    DEBER1 Elementodebarra-Alfonso Jimenez

    Cargado por

    Alfonso Jimenez
    Este documento presenta el análisis de un problema de elementos finitos de barras con 3 elementos. Se define la geometría, las propiedades de los materiales, las cargas y las restricciones. Luego, se ensambla la matriz de rigidez global y se resuelve el sistema de ecuaciones para obtener los desplazamientos nodales. Finalmente, se calculan los esfuerzos nodales, fuerzas nodales y se verifican las fuerzas internas.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    DEBER1 Elementodebarra-Alfonso Jimenez

    Cargado por

    Alfonso Jimenez
    26 vistas6 páginas
    Este documento presenta el análisis de un problema de elementos finitos de barras con 3 elementos. Se define la geometría, las propiedades de los materiales, las cargas y las restricciones. Luego, se ensambla la matriz de rigidez global y se resuelve el sistema de ecuaciones para obtener los desplazamientos nodales. Finalmente, se calculan los esfuerzos nodales, fuerzas nodales y se verifican las fuerzas internas.

    Descripción original:

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta el análisis de un problema de elementos finitos de barras con 3 elementos. Se define la geometría, las propiedades de los materiales, las cargas y las restricciones. Luego, se ensambla la matriz de rigidez global y se resuelve el sistema de ecuaciones para obtener los desplazamientos nodales. Finalmente, se calculan los esfuerzos nodales, fuerzas nodales y se verifican las fuerzas internas.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    26 vistas6 páginas

    DEBER1 Elementodebarra-Alfonso Jimenez

    Cargado por

    Alfonso Jimenez
    Este documento presenta el análisis de un problema de elementos finitos de barras con 3 elementos. Se define la geometría, las propiedades de los materiales, las cargas y las restricciones. Luego, se ensambla la matriz de rigidez global y se resuelve el sistema de ecuaciones para obtener los desplazamientos nodales. Finalmente, se calculan los esfuerzos nodales, fuerzas nodales y se verifican las fuerzas internas.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 6

    UNIVERSIDAD POLITECNICA SALESIANA

    _______________________________________________________________

    Fundamentos de Elementos Finitos

    Nombre: Alfonso Jimenez

    DEBER 1
    Calcular los esfuerzo nodales, deplazamientos y fuerzas del siguiente problema planteado.

    clear ;
    close all;

    Limpiar variable Workspace

    clc;

    PRE-PROCESAMIENTO

    Parametros de modelo

    Numero de elementos

    numeroelementos=3

    numeroelementos = 3

    Areas de elementos

    Area=[[5e-4];[3e-4];[1e-4]];
    Area_1=Area(1)

    Area_1 = 5.0000e-04

    Area_2=Area(2)

    Area_2 = 3.0000e-04

    1
    Area_3=Area(3)

    Area_3 = 1.0000e-04

    Longitud de elementos

    L=[[0.2];[0.3];[0.4]];
    L1=L(1)

    L1 = 0.2000

    L2=L(2)

    L2 = 0.3000

    L3=L(3)

    L3 = 0.4000

    Modulo de elasticidad elementos

    Modulo=[[200e9];[70e9];[200e9]]

    Modulo = 3×1
    1011 ×
    2.0000
    0.7000
    2.0000

    Modulo_1=Modulo(1)

    Modulo_1 = 2.0000e+11

    Modulo_2=Modulo(2)

    Modulo_2 = 7.0000e+10

    Modulo_3=Modulo(3)

    Modulo_3 = 2.0000e+11

    Fuerzas nodales

    F2=0 % Ya que no se aplica ninguna fuerza externe en este nodo

    F2 = 0

    F3=100e3 %Fuerza en N corresponde a 100kN

    F3 = 100000

    Restricciones

    U1=0

    U1 = 0

    2
    U4=0

    U4 = 0

    Conectividad

    Nodos1=[1 2];
    Nodos2=[2 3];
    Nodos3=[3 4];

    Coordenadas

    Coordenadas=[0;0.2;0.3;0.4]

    Coordenadas = 4×1
    0
    0.2000
    0.3000
    0.4000

    Grados de libertad

    GradoLibertad=4

    GradoLibertad = 4

    Matriz de rigidez de cada uno de los elementos

    k=zeros(2);
    k=cell(numeroelementos,1); %Sirve para crear un varible con nombre dinamico
    for Elementok=1:numeroelementos

    k{Elementok}=LinearBarElementStiffness(Modulo(Elementok),Area(Elementok),L(Elementok));
    end
    k_1=k{1} ;
    k_2=k{2} ;
    k_3=k{3} ;

    Inicializo matriz de rigidez global

    K=zeros(GradoLibertad);

    Ensamble de matriz de rigidez

    K=LinearBarAssemble(K,k_1,[1 2]);
    K=LinearBarAssemble(K,k_2,[2 3]);
    K=LinearBarAssemble(K,k_3,[3 4]);

    Reducir sistema de ecuaciones

    3
    Fa=[0;F3]

    Fa = 2×1
    0
    100000

    ka=K(2:3,2:3);

    Resolucion del sistema de ecuaciones reducidas

    Ua=ka\Fa;

    Vectores de desplazamientos globales

    U2=Ua(1);
    U3=Ua(2);
    U=[U1;U2;U3;U4]

    U = 4×1
    10-3 ×
    0
    0.1102
    0.8976
    0

    Vectores fuerza

    F=K*U

    F = 4×1
    105 ×
    -0.5512
    -0.0000
    1.0000
    -0.4488

    POST-PROCESAMIENTO

    Desplazamientos

    i=Nodos1(1,1); %asigna a i el valor i=1 representa lado izquierdo nodo


    j=Nodos1(1,2); %asigna a i el valor j=2 representa lado derechonodo

    4
    u_1=[U(i);U(j)]%Imprime los desplazamientos en los nodo #1

    u_1 = 2×1
    10-3 ×
    0
    0.1102

    i=Nodos2(1,1);
    j=Nodos2(1,2);

    u_2=[U(i);U(j)]%Imprime los desplazamientos en los nodo #2

    u_2 = 2×1
    10-3 ×
    0.1102
    0.8976

    i=Nodos3(1,1);
    j=Nodos3(1,2);
    u_3=[U(i);U(j)] %Imprime los desplazamientos en los nodo #3

    u_3 = 2×1
    10-3 ×
    0.8976
    0

    u=[u_1;u_2;u_3]

    u = 6×1
    10-3 ×
    0
    0.1102
    0.1102
    0.8976
    0.8976
    0

    Esfuerzos nodales

    s_1=LinearBarElementStresses(k_1,u_1,Area_1)

    s_1 = 2×1
    108 ×
    -1.1024
    1.1024

    s_2=LinearBarElementStresses(k_2,u_2,Area_2)

    s_2 = 2×1
    108 ×
    -1.8373
    1.8373

    s_3=LinearBarElementStresses(k_3,u_3,Area_3)

    s_3 = 2×1

    5
    108 ×
    4.4882
    -4.4882

    Calculo de fuerza nodales

    f_1=LinearBarElementForces(k_1,u_1)

    f_1 = 2×1
    104 ×
    -5.5118
    5.5118

    f_2=LinearBarElementForces(k_2,u_2)

    f_2 = 2×1
    104 ×
    -5.5118
    5.5118

    f_3=LinearBarElementForces(k_3,u_3)

    f_3 = 2×1
    104 ×
    4.4882
    -4.4882

    También podría gustarte