Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 83 vistas

    Informe Truss 3d

    Cargado por

    Karla Villegas
    Este documento presenta el código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

    Copyright:

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

    Informe Truss 3d

    Cargado por

    Karla Villegas
    83 vistas4 páginas
    Este documento presenta el código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

    Descripción original:

    Código y resultados de un elemento en 3 dimensiones por elementos finitos.

    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 código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

    Copyright:

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

    Informe Truss 3d

    Cargado por

    Karla Villegas
    Este documento presenta el código para resolver un problema de elementos finitos de una estructura tridimensional de vigas (truss) utilizando el método de los elementos finitos. Se define la geometría, condiciones de contorno y cargas de la estructura, y se ensambla la matriz de rigidez global considerando las matrices de rigidez locales de cada elemento viga. Luego, se resuelve el sistema para obtener los desplazamientos nodales, y con esto los esfuerzos y fuerzas internas en cada elemento. Finalmente, se muestran los resultados obten

    Copyright:

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

    ESCUELA POLITÉCNICA NACIONAL

    FACULTAD DE INGENIERÍA MECÁNICA


    ELEMENTOS FINITOS

    Villegas Ponce Karla Elizabeth


    Velasco Carrera Xavier Efren
    Gr 2
    10/01/2018

    DEBER

    TEMA: Programación por elementos finitos para el elemento barra 3D (Truss 3D).

    PROBLEMA:

    CODIGO:

    clear
    clc
    %% Datos de armadura 3D
    coord=[0 0 2000;4000 0 0;0 0 -2000;0 3000 0];
    conect=[1 2;2 3;2 4];
    f=[0 0 0 0 -10000 -12000 0 0 0 0 0 0]';
    nfijos=[1 2 3 7 8 9 10 11 12]';
    U=[0 0 0 0 0 0 0 0 0 0 0 0]';
    area=[100;100;100];
    mod_elas=[200e3;200e3;200e3];
    k1=[1 -1;-1 1]; % Matriz de rigidez elemental en ejes locales

    %% Rigidez axial
    r_axial=area.*mod_elas;

    %% Ensamble Matriz Global


    nnodos=size(coord,1);
    num_elem=size(conect,1);
    GDL=3*nnodos;
    K=zeros(GDL);
    for ii=1:num_elem
    % Nodos del elementos
    i=conect(ii,1);
    j=conect(ii,2);
    ESCUELA POLITÉCNICA NACIONAL
    FACULTAD DE INGENIERÍA MECÁNICA
    ELEMENTOS FINITOS

    % coordenadas nodales del elemento


    xn=[coord(i,1) coord(j,1)]';
    yn=[coord(i,2) coord(j,2)]';
    zn=[coord(i,3) coord(j,3)]';

    % Coordenadas nodales del elemento


    xi=xn(1); yi=yn(1); zi=zn(1);
    xj=xn(2); yj=yn(2); zj=zn(2);
    L=sqrt((xj-xi)^2+(yj-yi)^2+(zj-zi)^2);

    cx=(xj-xi)/L; %Coseno director X


    cy=(yj-yi)/L; %Coseno director Y
    cz=(zj-zi)/L; %Coseno director Z

    % matrices de transformacion
    R=[cx cx cz 0 0 0; 0 0 0 cx cy cz];
    X=mod_elas(1)*[-1/L(1 ) 1/L(1 )]*R;

    EA=r_axial(ii);
    k(ii)=(EA/L);
    %matriz de rigidez del elemento en coordenadas globales
    ke=(R')*((EA/L)*k1*R);

    in=[3*i-2 3*i-1 3*i 3*j-2 3*j-1 3*j];


    K(in,in)=K(in,in)+ke;
    end
    disp(K)

    %% SOLUCION DEL SISTEMA


    %Aplicacion de condiciones de contorno
    ic=nfijos;
    id=[1:GDL];%Define los indices de los GDL desconocidos
    id(ic)=[];
    Kmod=K([id, ic'],[id, ic']);
    dim=size(id, 2);
    K11=Kmod([1:dim],[1:dim]);
    UD=K11\f(id);

    U(4)=UD(1);
    U(5)=UD(2);
    U(6)=UD(3);
    U
    Fuerzas=Kmod*U
    L=[447.21 500 447.21];

    %% Esfuerzos
    X1=[-40.0000 -40.0000 20.0000 40.0000 0 -20.0000];
    X2=[40.0000 40.0000 20.0000 -40.0000 0 -20.0000];
    X3=[32 32 0 -32 24 0];
    Um=U(1:6);
    Esf1=X1*Um
    Esf2=X2*Um
    Esf3=X3*Um

    %% R obtenidos de cada elemento


    R1=[0.8944 0.8944 -0.4472 0 0 0; 0 0 0 0.8944 0 -0.4472];
    R2=[-0.8944 -0.8944 -0.4472 0 0 0; 0 0 0 -0.8944 0 -0.4472];
    ESCUELA POLITÉCNICA NACIONAL
    FACULTAD DE INGENIERÍA MECÁNICA
    ELEMENTOS FINITOS

    R3=[-0.8000 -0.8000 0 0 0 0 ; 0 0 0 -0.8000 0.6000 0];

    %% Fuerzas elementales
    % Para:
    disp('FUERZAS ELEMENTALES')
    Elemento1=k(1)*R1*Um
    Elemento2=k(2)*R2*Um
    Elemento3=k(3)*R3*Um

    RESULTADOS

    Matriz de rigidez

    Desplazamientos, fuerzas y esfuerzos: [mm] [N] [Pa] respectivamente


    ESCUELA POLITÉCNICA NACIONAL
    FACULTAD DE INGENIERÍA MECÁNICA
    ELEMENTOS FINITOS

    También podría gustarte