Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 620 vistas

    Tea 3 2019-02

    Cargado por

    mayli mayhua paucar
    Este documento presenta un proyecto para resolver una viga mediante diferentes métodos computacionales. Primero, se realiza un cálculo estático manual para hallar las reacciones, diagramas de corte y momento. Luego, se implementan los métodos de Bisección, Regula Falsi, Secante y Newton-Raphson en MATLAB para encontrar donde el momento es cero, y se grafican los resultados. Finalmente, se desarrolla una interfaz gráfica para mostrar los cálculos.

    Copyright:

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

    Tea 3 2019-02

    Cargado por

    mayli mayhua paucar
    620 vistas16 páginas
    Este documento presenta un proyecto para resolver una viga mediante diferentes métodos computacionales. Primero, se realiza un cálculo estático manual para hallar las reacciones, diagramas de corte y momento. Luego, se implementan los métodos de Bisección, Regula Falsi, Secante y Newton-Raphson en MATLAB para encontrar donde el momento es cero, y se grafican los resultados. Finalmente, se desarrolla una interfaz gráfica para mostrar los cálculos.

    Descripción original:

    trabajo IMC

    Título original

    TEA 3 2019-02

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta un proyecto para resolver una viga mediante diferentes métodos computacionales. Primero, se realiza un cálculo estático manual para hallar las reacciones, diagramas de corte y momento. Luego, se implementan los métodos de Bisección, Regula Falsi, Secante y Newton-Raphson en MATLAB para encontrar donde el momento es cero, y se grafican los resultados. Finalmente, se desarrolla una interfaz gráfica para mostrar los cálculos.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    620 vistas16 páginas

    Tea 3 2019-02

    Cargado por

    mayli mayhua paucar
    Este documento presenta un proyecto para resolver una viga mediante diferentes métodos computacionales. Primero, se realiza un cálculo estático manual para hallar las reacciones, diagramas de corte y momento. Luego, se implementan los métodos de Bisección, Regula Falsi, Secante y Newton-Raphson en MATLAB para encontrar donde el momento es cero, y se grafican los resultados. Finalmente, se desarrolla una interfaz gráfica para mostrar los cálculos.

    Copyright:

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

    CURSO:

    · INTRODUCCION A LOS METODOS COMPUTACIONALES

    TRABAJO:
    . “ELABORAR UNA PROGRAMACIÓN QUE NOS PERMITA RESOLVER
    UNA VIGA MEDIANTE LOS METODOS DE BISECCION, REGULA FALSI,
    SECANTE Y NEWTON RAPSHON”
    ALUMNO:

    · Mayhua Paucar, Boris Zósimo.

    PROFESOR:

    · Cuba Samaniego, Lisbeth Mónica.

    SECCIÓN:

    · CV64

    LIMA 2019-II
    INTRODUCCIÓN

    Este trabajo tiene como finalidad elaborar una programación que nos permita
    resolver una viga con la ayuda del programa MATLAB a través de diversos métodos
    vistos en el curso Introducción a los métodos computacionales. Desarrollaremos
    una programación para un caso en específico, hallando así los parámetros más
    importantes como son las reacciones, diagrama de fuerza cortante, momento
    flector, puntos donde el esfuerzo cortante y el momento sean cero. En primer lugar,
    elaboraremos una solución manual aplicando la estática para hallar las reacciones
    y los momentos en cada punto específico. Así mismo, graficaremos manualmente
    los diagramas de fuerza cortante y el diagrama de momento flector, hallando así
    todos los parámetros necesarios para resolver una viga. En segundo lugar,
    hallaremos los resultados obtenidos con el programa MATLAB, aplicaremos
    métodos de sistemas de ecuaciones no lineales para hallar el punto donde el
    momento flector sea cero, seguidamente utilizaremos el comando Guide para
    resolver los momentos y gráficas de esfuerzo cortante y momento flector.
    ANTECEDENTES

    Matlab (Matrix Laboratory, “laboratorios de matrices”). Es un Software


    matemático con entorno de desarrollo integrado (IDE) que tiene un lenguaje de
    programación propio (Lenguaje M) y es multiplataforma (Unix, Windows y Apple Mac
    Os X). Software de un gran uso en Centros de Investigación y Desarrollo, así como
    en universidades. En 1984 surge la primera versión, creado por Cleve Moler, con la
    idea de usar paquetes de subrutinas escritas en Fortran en los cursos de Análisis
    Numérico y Algebra Lineal, sin tener una necesidad de programas que usaran este
    lenguaje. El lenguaje de programación M se creó en 1970 proporcionando un
    acceso sencillo al Software de matrices LINPACK y EISPACK sin tener que hacer
    uso del lenguaje Fortran. Ya en 2004 se apreciaba que MATLAB era usado por
    aproximadamente más de un millón de personas, tanto académicos como
    empresarios.

    Funciones

    Dentro de sus principales funciones se encuentran:

    • Manipulación de Matrices.

    • La representación de datos y funciones.

    • Implementación de algoritmos.

    • Creación de interfaces de usuario (GUI).

    • Comunicación con programas en otros lenguajes y con otros dispositivos


    Hardware.

    Herramientas Adicionales

     Simulink (plataforma de simulación multidominio).


     GUIDE (editor de interfaces de usuario - GUI).
    DESARROLLO DEL TRABAJO

    1. PROBLEMA ASIGNADO:
    Se carga una viga de a manera que se aprecia en la figura, Emplee los
    métodos de Bisección, Regula falsi, Secante y Newton Rapshon para
    resolver la posición de la viga donde no hay momento.

    Unidades:
    Ib=libras
    ft=pies

    2. CALCULO ESTATICO DE LA VIGA:


    a) Calculo de reacciones:
     ∑ma=0
    −𝟏𝟓𝟎(𝟐) − 𝟑𝟎𝟎(𝟒. 𝟓) + 𝐁𝐯(𝟏𝟎) − 𝟏𝟎𝟎(𝟏𝟐) = 𝟎
    𝐁𝐯 = 𝟐𝟖𝟓 𝐈𝐛 (⤒)
     ∑fx=0
    𝐀𝐡 = 𝟎
     ∑fy=0
    𝐀𝐯 − 𝟏𝟓𝟎 − 𝟑𝟎𝟎 + 𝟐𝟖𝟓 − 𝟏𝟎𝟎 = 𝟎
    𝐀𝐯 = 𝟐𝟔𝟓 𝐈𝐛 (⤒)

    b) Análisis de ecuaciones de cortantes y momentos. (Método Cortes)


    TRAMO 1: (𝟎 < 𝐗 < 𝟑)

     ∑m(x)=0
    𝟓𝟎 𝐱 𝟑
    𝐌(𝐱) = − + 𝟐𝟔𝟓𝐱
    𝟗
    𝐌(𝟎) = 𝟎 𝐈𝐛. 𝐟𝐭
    𝐌(𝟑) = 𝟔𝟒𝟓 𝐈𝐛. 𝐟𝐭

     ∑fx=0
    𝐍=𝟎
     ∑fy=0
    𝟓𝟎 𝐱 𝟐
    𝐕(𝐱) = 𝟐𝟔𝟓 −
    𝟑
    𝐕(𝟎) = 𝟐𝟔𝟓 𝐈𝐛
    𝐕(𝟑) = 𝟏𝟏𝟓 𝐈𝐛

    TRAMO 2: (𝟑 < 𝐗 < 𝟔)


     ∑m(x)=0
    𝐌(𝐱) = −𝟓𝟎𝐱 𝟐 + 𝟒𝟏𝟓𝐱 − 𝟏𝟓𝟎
    𝐌(𝟑) = 𝟔𝟒𝟓 𝐈𝐛. 𝐟𝐭
    𝐌(𝟔) = 𝟓𝟒𝟎 𝐈𝐛. 𝐟𝐭

     ∑fx=0
    𝐍=𝟎
     ∑fy=0
    𝐕(𝐱) = 𝟐𝟔𝟓 − 𝟏𝟓𝟎 − 𝟏𝟎𝟎(𝐱 − 𝟑)
    𝐕(𝟑) = 𝟏𝟏𝟓 𝐈𝐛
    𝐕(𝟔) = −𝟏𝟖𝟓 𝐈𝐛

    TRAMO 3: (𝟔 < 𝐗 < 𝟏𝟎)

     ∑m(x)=0
    𝐌(𝐱) = −𝟏𝟖𝟓𝐱 + 𝟏𝟔𝟓𝟎
    𝐌(𝟔) = 𝟓𝟒𝟎 𝐈𝐛. 𝐟𝐭
    𝐌(𝟏𝟎) = −𝟐𝟎𝟎 𝐈𝐛. 𝐟𝐭

     ∑fx=0
    𝐍=𝟎
     ∑fy=0
    𝐕 = −𝟏𝟖𝟓 𝐈𝐛

    TRAMO 3: (𝟏𝟎 < 𝐗 < 𝟏𝟐)

     ∑m(x)=0
    𝐌(𝐱) = 𝟏𝟎𝟎𝐱 − 𝟏𝟐𝟎𝟎
    𝐌(𝟏𝟎) = −𝟐𝟎𝟎 𝐈𝐛. 𝐟𝐭
    𝐌(𝟏𝟐) = 𝟎 𝐈𝐛. 𝐟𝐭

     ∑fx=0
    𝐍=𝟎

     ∑fy=0
    𝐕 = 𝟏𝟎𝟎 𝐈𝐛
    c) Diagramas de fuerzas cortantes y momentos flector
    DIAGRAMA DE FUERZA CORTANTE

    DIAGRAMA DE MOMENTO FLECTOR


    d) Ecuaciones de momentos en cada tramo:
    50 x3
    f(x) = − + 265x, (0 ≤ X ≤ 3)
    9

    f(x) = −50x 2 + 415x − 150, (3 < X ≤ 6)


    f(x) = −185x + 1650, (6 < X ≤ 10)
    f(x) = 100x − 1200, (10 < X ≤ 12)

    e) Determinando la posición en el cual el momento es cero:


    Remplazar la función cuando f(x)=0 , para calcular la posición en la cual
    el momento es cero.
    Como evaluamos cada tramo y obtuvimos la gráfica del momento flector,
    podemos ubicar en que función evaluar para determinar la posición donde
    el momento es cero:

    𝐟(𝐱) = −𝟏𝟖𝟓𝐱 + 𝟏𝟔𝟓𝟎, (𝟔 < 𝐗 ≤ 𝟏𝟎)

    𝟎 = −𝟏𝟖𝟓𝐱 + 𝟏𝟔𝟓𝟎
    𝐱 = 𝟖. 𝟗𝟏𝟖𝟗 𝐟𝐭
    3. CALCULO DE MOMENTO UTILIZANDO PROGRAMA MATLAB:
    a) Graficando:
     Graficando método 1

    x=linspace(0,12,100);
    y1=((-50*x.^3)./9+265*x).*(0<=x & x<=3);
    y2=(-50*x.^2+415*x-150).*(3<x & x<=6);
    y3=(-185*x+1650).*(6<x & x<=10);
    y4=(100*x-1200).*(10<x & x<=12);
    f=y1+y2+y3+y4;
    plot(x,f)
    grid on
     Grafica método 2
    x1=0:0.1:3;
    x2=3:0.1:6;
    x3=6:0.1:6;
    x4=10:0.1:12;
    Mx1=(-50*x1.^3)./9+265*x1;
    Mx2=-50*x2.^2+415*x2-150;
    Mx3=-185*x3+1650;
    Mx4=100*x4-1200;
    x=[x1 x2 x3 x4];
    Mx=[Mx1 Mx2 Mx3 Mx4];
    plot(x,Mx)
    grid on
    b) Calculando la posición del momento cuando es cero por diferentes
    métodos:
     Método Bisección
    function [n,m] = bisec(a,b,er2)

    er1=100;
    m0=100;
    for n=0:100
    m=(a+b)/2;
    %Fa=-185*a+1650
    Fa=viga1(a);
    %Fm=-185*m+1650
    Fm=viga1(m);

    if Fa*Fm<0
    b=m;
    else
    a=m;
    end
    er1=(abs(m-m0)/m)*100;
    m0=m;
    if er1<=er2
    break;
    end
    end

     Método Regula Falsi


    function [n,m] = RegulaFalsi(a,b,er2)
    er1=100;
    m0=100;
    for n=0:100
    %Fa=-185*a+1650;
    Fa=viga1(a);
    %Fb=-185*b+1650;
    Fb=viga1(b);
    m=(a*Fb-b*Fa)/(Fb-Fa);
    %Fm=-185*m+1650;
    Fm=viga1(m);

    if Fa*Fm<0
    b=m;
    else
    a=m;
    end

    er1= (abs(m-m0)/m)*100;
    m0=m;
    if er1<=er2
    break;
    end
    end
     Método Secante
    function [n ,m ] = secante(x0,x1,er2)
    er1=100;
    m0=100;
    for n=0:100
    Fx0=viga1(x0);
    Fx1=viga1(x1);
    m=(x0*Fx1-x1*Fx0)/(Fx1-Fx0);
    x0=x1;
    x1=m;
    er1=abs((m-m0)/m)*100;
    m0=m;
    if er1<=er2
    break
    end
    end
    end

     Método Newton Rapshon


    function [n,m] = newton100(x0,er1)
    er2=0;
    m0=0;
    for n=0:100
    Fx0=-185*x0+1650;
    dFx0=-185;
    m=x0-Fx0/dFx0;
    x0=m;
    er2= abs((m-m0)/m)*100;
    m0=m;
    if er2<=er1
    break
    end
    end
    4. INTERFAZ GUIDE PROGRAMACION MATLAB
     IMPORTAR GRAFICA

    x=linspace(2,10,100)
    axes(handles.axes1);
    imshow(imread('viga1.png'));

     GRAFICA DE MOMENTO FLECTOR

    x1=0:0.1:3;
    x2=3:0.1:6;
    x3=6:0.1:6;
    x4=10:0.1:12;
    Mx1=(-50*x1.^3)./9+265*x1;
    Mx2=-50*x2.^2+415*x2-150;
    Mx3=-185*x3+1650;
    Mx4=100*x4-1200;
    x=[x1 x2 x3 x4];
    f=[Mx1 Mx2 Mx3 Mx4];
    axes(handles.axes2);
    plot(x,f)
    grid on
    title('DIAGRAMA DE MOMENTO FLECTOR')

     RESULTADOS DE LOS METODOS CODIFICACIÓN

    function popupmenu1_Callback(hObject, eventdata, handles)


    % hObject handle to popupmenu1 (see GCBO)
    % eventdata reserved - to be defined in a future version of MATLAB
    % handles structure with handles and user data (see GUIDATA)

    global a;
    global b;

    a=eval(get(handles.edit1,'string'));
    b=eval(get(handles.edit2,'string'));
    er2=eval(get(handles.edit3,'string'));

    contenido=get(hObject,'string');
    a1=get(hObject,'value');
    tipo=contenido(a1);

    er1=100;
    m0=100;
    switch cell2mat(tipo)

    case 'BISECCION'
    [n,m0] = bisec(a,b,er2);
    datos1=[n m0];
    set(handles.uitable1,'data',datos1);

    case 'REGULA FALSI'


    [n,m] = RegulaFalsi(a,b,er2);
    datos2=[0 0; n m0];
    set(handles.uitable1,'data',datos2);

    case 'SECANTE'
    x0=a;
    x1=b;
    [n ,m0 ] = secante(x0,x1,er1);
    datos3=[0 0; 0 0; n m0];
    set(handles.uitable1,'data',datos3);

    case 'NEWTON RAPSHON'


    x0=a;
    [n,m] = newton100(x0,er1)
    datos4=[0 0; 0 0; n m0];
    set(handles.uitable1,'data',datos4);
    end
    CONCLUSIONES
    • Se concluye que el uso de Matlab ayudo de forma óptima a desarrollar nuestro trabajo,
    ya que el proceso de programación coincidió a los resultados esperados que resolvimos
    manualmente. Resultó una buena herramienta de comprobación.
    • El programa nos ayudara a reducir el tiempo empleado en los cálculos de vigas
    hiperestáticas y al posible error que pueda surgir en las gráficas.
    • El programa servirá como una herramienta de corroboración para el estudio de los
    cursos de estructura.

    BIBLIOGRAFÍA
     Matlab para ingenieros. Manipulación de matrices Matlab. Disponible en:(
    http://dea.unsj.edu.ar/control2/matlab%20para%20ingenieros.pdf ). (Fecha de consulta: 6
    de mayo del 2019).
     Mathworks.Math.Graphics.Programming. Disponible en:(
    https://www.mathworks.com/products/matlab.html?s_tid=hp_products_matlab). (Fecha de
    consulta: 2 de mayo del 2019).
     Teoría y práctica del modelado de procesos mediante diagramas de flujos. Reglas
    de diseño. Disponible en (
    https://www.jramonet.com/sites/default/files/adjuntos/diagramas_flujo_jrf_v2013 .pdf ).
    (Fecha de consulta: 6 de mayo del 2019)

    También podría gustarte