Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 24 vistas

    Laboratorio N°6 Marcos - Levano - Hernandez - Carlos

    Cargado por

    samuel
    Este informe presenta los resultados de un análisis de elementos finitos de un marco bidimensional realizado en MATLAB y ANSYS. Se analizó el marco con diferentes números de elementos finitos, y se compararon las fuerzas de reacción, diagramas de esfuerzos y deformaciones. Se concluyó que a mayor número de elementos, mayor precisión en los resultados.

    Copyright:

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

    Laboratorio N°6 Marcos - Levano - Hernandez - Carlos

    Cargado por

    samuel
    24 vistas18 páginas
    Este informe presenta los resultados de un análisis de elementos finitos de un marco bidimensional realizado en MATLAB y ANSYS. Se analizó el marco con diferentes números de elementos finitos, y se compararon las fuerzas de reacción, diagramas de esfuerzos y deformaciones. Se concluyó que a mayor número de elementos, mayor precisión en los resultados.

    Descripción original:

    Título original

    LABORATORIO N°6 MARCOS_LEVANO_HERNANDEZ_CARLOS

    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 informe presenta los resultados de un análisis de elementos finitos de un marco bidimensional realizado en MATLAB y ANSYS. Se analizó el marco con diferentes números de elementos finitos, y se compararon las fuerzas de reacción, diagramas de esfuerzos y deformaciones. Se concluyó que a mayor número de elementos, mayor precisión en los resultados.

    Copyright:

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

    Laboratorio N°6 Marcos - Levano - Hernandez - Carlos

    Cargado por

    samuel
    Este informe presenta los resultados de un análisis de elementos finitos de un marco bidimensional realizado en MATLAB y ANSYS. Se analizó el marco con diferentes números de elementos finitos, y se compararon las fuerzas de reacción, diagramas de esfuerzos y deformaciones. Se concluyó que a mayor número de elementos, mayor precisión en los resultados.

    Copyright:

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

    UNIVERSIDAD NACIONAL DE INGENIERÍA

    FACULTAD DE INGENIERÍA MECÁNICA

    INFORME DEL LABORATORIO N°6


    “MARCOS”

    Curso : Cálculo Por Elementos Finitos – MC516


    Sección : “F “
    Profesor: Nain Máximo Ramos Álvarez
    Alumno : Lévano Hernández Carlos Samuel- 20135001j
    Fecha de entrega:
    07 de enero del 2023
    Periodo académico:

    2022-II

    LIMA - PERÚ
    DIGITACIÓN Y EJECUCIÓN EN MATLAB

    % MARCOS
    clearvars, clc;

    E_vec=200000*[1,1,1];
    A_vec=2343.75*[1,1,1];
    I_vec=7446289.0625*[1,1,1];

    % generation of coordinates and connectivities


    % x1=3000*(1+cos(pi/4))
    nodeCoordinates=[0 0; 1500 3000;4500 3000; 6000 0];
    elementNodes=[1 2; 2 3; 3 4];
    numberNodes=size(nodeCoordinates,1);

    % GDof: global number of degrees of freedom


    GDof=3*numberNodes;

    % Assembly stiffness matrix


    K_Assembly=formStiffness2Dframe(GDof,size(elementNodes,1),elementNodes,nodeCoor
    dinates,E_vec,I_vec,A_vec);

    %force vector
    F_col=nan(GDof,1);
    F_col(4)=0;
    F_col(5)=-52974;
    F_col(6)=0;
    F_col(7)=0;
    F_col(8)=-52974;
    F_col(9)=0;

    %displacement vector
    prescribedDof=[1 2 3 10 11 12];
    D_col=nan(GDof,1);
    D_col(prescribedDof)=0;

    % solution
    [D_col,F_col]=solution(prescribedDof,K_Assembly,D_col,F_col);

    %IMPRESION DE RESULTADOS
    i=1;
    for j=1:numberNodes
    fprintf('La deformacion en mm en el nodo %d es:\n',j);
    fprintf('\tEn el eje X: %f\n',D_col(i));
    fprintf('\tEn el eje Y: %f\n',D_col(i+1));
    fprintf('\tRotacion(°): %f\n',D_col(i+2));
    i=i+3;
    end
    fprintf('\n')
    i=1;
    for j=1:numberNodes
    fprintf('Las fuerzas en Newton en el nodo %d es:\n',j);
    fprintf('\tEn el eje X: %f\n',F_col(i));
    fprintf('\tEn el eje Y: %f\n',F_col(i+1));
    fprintf('\tMomento Flector(N.mm): %f\n',D_col(i+2));
    i=i+3;
    end

    function
    K_asembly=formStiffness2Dframe(GDof,numberElements,elementNodes,nodeCoordinates
    ,E_vec,I_vec,A_vec)
    %This corrects the strange node ordering of Ferreira's book

    N_node_DOF=3; %2D problem


    K_asembly=zeros(GDof);

    for iElement=1:numberElements
    iNodes=elementNodes(iElement,:);
    x_local=nodeCoordinates(iNodes(2),:)-nodeCoordinates(iNodes(1),:);
    L=norm(x_local);

    E=E_vec(iElement);

    if iElement==2
    A=2343.75/2;
    II=3051757.8125;
    else
    A=2343.75;
    II=7446289.0625;
    end

    k_local=E*[ A/L,0,0,-A/L,0,0
    0,12*II/L^3,6*II/L^2,0,-12*II/L^3,6*II/L^2
    0,6*II/L^2,4*II/L,0,-6*II/L^2,2*II/L
    -A/L,0,0,A/L,0,0
    0,-12*II/L^3,-6*II/L^2,0,12*II/L^3,-6*II/L^2
    0,6*II/L^2,2*II/L,0,-6*II/L^2,4*II/L];

    l=x_local(1)/L;
    m=x_local(2)/L;

    T_node=[l,m,0
    -m,l,0
    0,0,1];
    T_element=[T_node,zeros(N_node_DOF);zeros(N_node_DOF),T_node];

    k_global=T_element.'*k_local*T_element;

    elementDof=[N_node_DOF*iNodes(1)+(-N_node_DOF:-1)+1,N_node_DOF*iNodes(2)+(-
    N_node_DOF:-1)+1];
    K_asembly(elementDof,elementDof)=K_asembly(elementDof,elementDof)+k_global;
    end

    function [D_vec,F_vec]=solution(prescribedDof,K_assembly,D_vec,F_vec,F_eq_vec)
    if nargin<5
    F_eq_vec=zeros(size(F_vec));
    end

    GDof=length(D_vec);

    freeDof=setdiff(1:GDof,prescri

    bedDof);

    D_vec(freeDof)=K_assembly(freeDof,freeDof)\(F_vec(freeDof)+F_eq_vec(free
    Dof));

    %nonZeroDof=find(D_vec~=0);
    nonZeroDof=union(freeDof,prescribedDof(D_vec(prescribedDo
    f)~=0));

    F_vec(prescribedDof)=K_assembly(prescribedDof,nonZeroDof)*D_vec(nonZ
    eroDof)-F_eq_vec(prescribedDof);

    RESULTADOS DEL SOFTWARE MATLAB

    1
    DESARROLLO EN ANSYS

    GEOMETRY

    MODEL CON 10 ELEMENTOS FINITOS

    2
    MODEL CON 5 ELEMENTOS FINITOS

    MODEL CON 3 ELEMENTOS FINITOS

    MODEL CON 1 ELEMENTOS FINITOS

    3
    FUERZA DE REACCIÓN ENLOS APOYOS
    FUERZA DE REACCION CON 10 ELEMENTOS FINITOS

    FUERZA DE REACCION CON 5 ELEMENTOS FINITOS

    4
    FUERZA DE REACCION CON 3 ELEMENTOS FINITOS

    FUERZA DE REACCION CON 1 ELEMENTOS FINITOS

    5
    DIAGRAMA DE ESFUERZOS

    DIAGRAMA DE ESFUERZOS CON 10 ELEMENTOS FINITOS

    DIAGRAMA DE ESFUERZOS CON 5 ELEMENTOS FINITOS

    6
    DIAGRAMA DE ESFUERZOS CON 3 ELEMENTOS FINITOS

    DIAGRAMA DE ESFUERZOS CON 1 ELEMENTOS FINITOS

    7
    FORMA DEFORMADA DEL MARCO

    FORMA DEFORMADA CON 10 ELEMENTOS FINITOS

    FORMA DEFORMADA CON 5 ELEMENTOS FINITOS

    8
    FORMA DEFORMADA CON 3 ELEMENTOS FINITOS

    FORMA DEFORMADA CON 1 ELEMENTOS FINITOS

    9
    LA DEFORMADA DE LA VIGA

    DEFORMADA DE LA VIGA CON 10 ELEMENTOS FINITOS

    DEFORMADA DE LA VIGA CON 5 ELEMENTOS FINITOS

    DEFORMADA DE LA VIGA CON 3 ELEMENTOS FINITOS

    1
    0
    DEFORMADA DE LA VIGA CON 1 ELEMENTOS FINITOS

    1
    1
    CONCLUSIONES DEL
    TRABAJO
    1. Es mejor trabajar con un mayor número de elementos finitos para mejor precisión enlos resultados

    a obtener

    2. Los valores de las fuerzas resultantes difieren debido a que mientras se use un menornúmero de

    elementos finitos se tendrá menor exactitud

    3. Lo mismo ocurre con el grafico de diagrama de esfuerzos donde se puede observar que mientras

    se use menos elementos finitos, las gráficas que se generan tienden a serpoco exactas

    4. El mismo caso ocurre para el marco deformado donde se observa que al momento desometerlo a

    un solo elemento finito no ocurre un cambio significativo en la viga.

    1
    2
    1
    3

    14

    También podría gustarte