Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 36 vistas

    Fexibilidad

    Cargado por

    Lyco Ceykhann Ramos Cavalcanti
    Este documento presenta el análisis estructural de una viga de 3 nudos y 3 barras mediante el método de los elementos finitos. Se define la matriz de incidencia, la matriz de equilibrio y la matriz de flexibilidad. Luego, se aplican cargas puntuales y distribuidas, y se resuelve el sistema mediante la matriz general para obtener los vectores de fuerzas y desplazamientos.

    Copyright:

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

    Fexibilidad

    Cargado por

    Lyco Ceykhann Ramos Cavalcanti
    36 vistas3 páginas
    Este documento presenta el análisis estructural de una viga de 3 nudos y 3 barras mediante el método de los elementos finitos. Se define la matriz de incidencia, la matriz de equilibrio y la matriz de flexibilidad. Luego, se aplican cargas puntuales y distribuidas, y se resuelve el sistema mediante la matriz general para obtener los vectores de fuerzas y desplazamientos.

    Descripción original:

    fexibilidad

    Título original

    fexibilidad

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

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

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta el análisis estructural de una viga de 3 nudos y 3 barras mediante el método de los elementos finitos. Se define la matriz de incidencia, la matriz de equilibrio y la matriz de flexibilidad. Luego, se aplican cargas puntuales y distribuidas, y se resuelve el sistema mediante la matriz general para obtener los vectores de fuerzas y desplazamientos.

    Copyright:

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

    Fexibilidad

    Cargado por

    Lyco Ceykhann Ramos Cavalcanti
    Este documento presenta el análisis estructural de una viga de 3 nudos y 3 barras mediante el método de los elementos finitos. Se define la matriz de incidencia, la matriz de equilibrio y la matriz de flexibilidad. Luego, se aplican cargas puntuales y distribuidas, y se resuelve el sistema mediante la matriz general para obtener los vectores de fuerzas y desplazamientos.

    Copyright:

    © All Rights Reserved
    Descargue como TXT, PDF, TXT o lea en línea desde Scribd
    Descargar como txt, pdf o txt
    Está en la página 1de 3

    %% 1.

    Datos
    E=1 % M�dulo de elasticidad
    I=1 % Momento de Inercia
    A=1 % Area de la secci�n

    N=4 % N�mero de nudos


    M=3 % N�mero de barras

    %% matriz de incidencia

    a=[ 1 2 3 4 5 6
    4 5 6 7 8 9
    7 8 9 10 11 12]

    %% vector de ordenamiento
    aa=zeros(M,3);

    aa=[ 1 2 3
    4 5 6
    7 8 9]

    %% calculamos la matriz de equilibrio

    C=zeros(3*N,3*M);

    [C1,L1,cs1,sn1]=axibend_C(0,0,2.5,5);
    [C2,L2,cs2,sn2]=axibend_C(0,5,5,5);
    [C3,L3,cs3,sn3]=axibend_C(5,5,5,0);

    C1;
    C2;
    C3;

    C(a(1,:)',aa(1,:))=C1;
    C(a(2,:)',aa(2,:))=C2 ;
    C(a(3,:)',aa(3,:))=C3;
    disp( 'matriz de equilibrio ensamblada')
    C % Matriz de equilibrio completo.
    % suprimimos las restricciones
    C([1,2,3,10,11,12],:)=[];

    disp( 'matriz de equilibrio reducida')


    C %Matriz de equilibrio con restricciones.

    %% calculamos la matriz de Flexibilidad


    %Longitud de la barra
    L1 %metros
    L2 %metros
    L3 %metros

    % Matriz de Flexibilidad
    J=zeros(3*M,3*M);

    J1 = A_B_J(L1,E,A,I)
    J2 = A_B_J(L2,E,A,I)
    J3= A_B_J(L3,E,A,I)
    J(aa(1,:)',aa(1,:))=J1;
    J(aa(2,:)',aa(2,:))=J2 ;
    J(aa(3,:)',aa(3,:))=J3;

    disp( 'matriz de Flexibilidad')


    J

    %% PRIMER CASO
    %definici�n del vector de cargas
    %
    Q=zeros(N*3,1);
    Ql=zeros(N*3,1);
    % % % CaRGA Local
    % pd=1 % Carga Distribuida ( kgf/m)
    % L=L3 % Longitud de la barra
    %
    % P=pd*L %P resultante de la carga normal
    % T=0; %T resultante de la carga tangencial
    % rp=0.5; %rP distancia de P al primer nudo, fracci�n de L
    %
    %
    % [QL] = axibend_Ql(L,cs3,sn3,P,T,rp);
    % Q(a(3,:),1)= Q(a(3,:),1)+(QL);

    % Carga en el nodo

    Q(4,1)=5;

    % Carga distribuidas en cada barra

    Ql(5,1)=2.5;
    Ql(8,1)=2.5;

    Q([1,2,3,10,11,12],:)=[]
    Ql([1,2,3,10,11,12],:)=[]

    Qr=Q-Ql

    %% Resolvemos por el m�todo general


    disp( 'Matriz General')
    G=[ J, -C'
    C, zeros(6,6)] % matriz General

    fl=zeros(M*3,1);

    fl(4,1)=5.208;
    fl(5,1)=-5.208;

    b1= [fl;Qr] % Matriz b

    x1= G\b1;

    % vector de fuerzas

    disp('=================================================')
    disp('Resultados')
    F1=x1(1:9)

    % vector de desplazamientos
    qr1=x1(10:15)

    %bPrint2('F1','qr1');

    También podría gustarte