% Descripcion del programa: %Este programa determina las fuerzas que se generan en dos nodos. %Autor: Edwin Molina %Fecha: 20.06.2018 E=input('Ingrese el m�dulo de elasticidad del material [N/m^2]: '); disp('La seccion es rectangular') b=input('Ingrese la medida de la base (b) [mm]: '); h=input('Ingrese la medida de la altura (h) [mm]; '); I=(1/12)*(b/1000)*(h/1000)^3; L=input('Ingrese el valor de la longitud [mm]: '); L1=L/1000; keq=(3*E*I)/(L1)^3; u1=input('Ingrese el maximo desplazamiento del nodo 1: '); u2=input('Ingrese el maximo desplazamiento del nodo 2: '); %Definir la matriz de rigidez K K=[keq -keq;-keq keq]; disp('La matriz de rigidez es: ') disp(K) %Definir la matriz de desplazamiento U=[u1; u2]; disp('La matriz de desplazamiento es:') disp(U) %Calcular las fuerzas F1 y F2 F=K*U; disp('La matriz de fuerza es:') disp(F)