DIARIO DEL EXAMEN RESUELTO EN MATLAB

% Primero vamos a encontrar las matrices de rigidez en coordenadas locales% Para elmentos 1234 es la misma rigidez localhelp rlocald [k]=rlocal(g1,g2,g3,l,f,d,e,miu) la rigidez local se obtiene a partir de los siguientes datos g1,g2,g3 coeficientes de guldan 4,4,2 l= longitud de elemento f= fator de forma d= diámetro interno e=espesor sección [k1234]=rlocald(4,4,2,6,2,0.250,0.009,0.3);[k5]=rlocald(4,4,2,5,2,0.250,0.009,0.3);[k6]=rlocald(4,4,2,sqrt(73),2,0.250,0.009,0.3);help rlocalt [k]=rlocal(d,e,I) la rigidez local se obtiene a partir de los siguientes datos l= longitud de elemento d=diametro exterior e=espesor de seccion [k7]=rlocalt(0.08,0.006,sqrt(73); [k7]=rlocalt(0.08,0.006,sqrt(73); |{_Error: Unbalanced or unexpected parenthesis or bracket.}_ [k7]=rlocalt(0.08,0.006,sqrt(73));[k8]=rlocalt(0.08,0.006,sqrt(145));[k9]=rlocalt(0.08,0.006,sqrt(289));% Ahora encontramos las matricez de transformacion de coordenadas % para elementos 1234 es la misma matrizhelp transf [T]=transf(alfa,beta1,beta2,l,dx1,dy1,dx2,dy2) la matriz de transformacion cambia las acciones de ejes locales a coordenadas de nudo alfa=angulo de eje local beta1=angulo nudo inicial beta2=angulo nudo final l=longitud dx1,dy1,dx2,dy2, longitudes de elementos rigidos en ejes locales [T1234]=transf(0,0,0,6,0,0,0,0);[T5]=transf(90,0,0,5,0,0,0,0);[T6]=transf(110.55,0,0,sqrt(73),0,0,0,0);help transfarm [T]=transfarm(alfa,beta1,beta2,l,dx1,dy1,dx2,dy2) la matriz de transformacion cambia las acciones de ejes locales a coordenadas de nudo alfa=angulo de eje local beta1=angulo nudo inicial beta2=angulo nudo final

l=longitud dx1,dy1,dx2,dy2, longitudes de elementos rigidos en ejes locales [T7]=transfarm(69.44,0,0,sqrt(73),0,0,0,0);[T8]=transfarm(138.36,0,0,sqrt(145),0,0,0,0);[T9]=transfarm(151.92,0,0,sqrt(289),0,0,0,0);% Encontrar ahora las matrices de rigidez en coordenadas de NudoK1234=T1234*k1234*T1234';K5=T5*k5*T5';K6=T6*k6*T6';K7=T7*k7*T7';K8=T8*k8*T8';K9=T9*k9*T9';% Ensamble de la matriz de rigidez DE LA ESTRUCTURA% primero definimos los vectores Leelee1=[15 16 1 2 3 4]; lee2=[2 3 4 5 6 7]; lee3=[5 6 7 8 9 10]; lee4=[8 9 10 17 18 11];lee5=[19 20 21 2 3 4]; lee6=[2 3 4 12 13 14]; lee7=[15 16 1 12 13 14]; lee8=[5 6 12 13]; lee9=[8 9 12 13]; lee5=[19 20 21 2 3 4]; lee6=[2 3 4 12 13 14]; lee7=[15 16 12 13]; lee8=[5 6 12 13]; lee9=[8 9 12 13];% Incialisamos la matriz de rigidez de la estructura y el vector de cargasS=zeros(14); B=zeros(21,1);[S,B]=ensamble(lee1,K1234,14,S,B);[S,B]=ensamble(lee2,K1234,14,S,B);[S,B]=ensamble(lee3,K1234,14,S,B);[S,B]=ensamble(lee4,K1234,14,S,B);[S,B]=ensamble(lee5,K5,14,S,B);[S,B]=ensamble(lee6,K6,14,S,B);[S,B]=ensamble(lee7,K7,14,S,B);[S,B]=ensamble(lee8,K8,14,S,B);[S,B]=ensamble(lee9,K9,14,S,B);% Ahora definimos vector de acciones de empotramiento perfectoaep1234=[0 4.5 4.5 0 4.5 -4.5]';[Ta1234]=transfEMP(0,0,0,6,0,0,0,0); %Matriz de transformacion de coordenadas de acciones de empotramiento%Crear el vector AEP en coordenadas de nudoAEP1234=Ta1234*aep1234;% Encontramos el vector AD-ADLhelp acum [PR]=acum[lee,Ae,PR] [B]=acum[lee1,-AEP1,B]; [B]=acum[lee1,-AEP1,B]; |{_Error: Unbalanced or unexpected parenthesis or bracket.}_ [B]=acum(lee1,-AEP1,B);{_Undefined function or variable 'AEP1'.}_ [B]=acum(lee1,-AEP1234,B);[B]=acum(lee2,-AEP1234,B);[B]=acum(lee3,-AEP1234,B);[B]=acum(lee4,-AEP1234,B);% Desplazamientos en coordenadas de nudODelta=S\B(1:14);

% Ahora encontramos los desplazamientos de cada elemento en coordenadas de nudo para esto primero llenar el vector Delta(21)=0D1=Delta(lee1); D2=Delta(lee2); D3=Delta(lee3); D4=Delta(lee4); D5=Delta(lee5); D6=Delta(lee6); D7=Delta(lee7); D8=Delta(lee8);D= [D1 D2 D3 D4 D5 D6 D7 D8 D9];{_Undefined function or variable 'D9'.}_ D9=Delta(lee9);D= [D1 D2 D3 D4 D5 D6 D7 D8 D9];{_Error using <a href="matlab:helpUtils.errorDocCallback('horzcat')" style="font-weight:bold">horzcat</a>CAT arguments dimensions are not consistent.}_ % Encontrar ahora los desplazamientos en coordendas reducidasd1=T1234'*D1;d2=T1234'*D2;d3=T1234'*D3;d4=T1234'*D4;d5=T5'*D5;d6=T6'*D6;d7=T7'*D7;d8=T8'*D8;D=[D1 D2 D3 D4 D5 D6];D2=[D7 D8 D9];D20=[D7 D8 D9];D2=Delta(lee2);d2=T1234'*D2;% Acciones de elemento en coordenadas locales reducidasa1r=k1234*d1;a2r=k1234*d2;a3r=k1234*d3;a4r=k1234*d4;a5r=k5*d5;a6r=k6*d6;a7r=k7*d7;a8r=k8*d8;a9r=k9*d9;{_Undefined function or variable 'd9'.}_ d9=T9'*D9;a9r=k9*d9;% Ahora encontrar las acciones en coordenadas reducidad completas para esto utilizo una matriz de equilibrio[h]=mequilibriomp(6);a1c=h*a1r+aep1234;a2c=h*a2r+aep1234;a3c=h*a3r+aep1234;a4c=h*a4r+aep1234;[h]=mequilibriomp(5);a5c=h*a5r; %No hay acciones de empotramientos para elementos 56789[h]=mequilibriomp(sqrt(73); [h]=mequilibriomp(sqrt(73); |{_Error: Unbalanced or unexpected parenthesis or bracket.

}_ [h]=mequilibriomp(sqrt(73));a6c=h*a6r;[h]=mequilibriomp(sqrt(73));a7c=h*a7r;[h]=mequilibriomp(sqrt(145));a8c=h*a8r;[h]=mequilibriomp(sqrt(289));a9c=h*a9r;%Acciones finales en coordenadas de nudoA1=Ta1234*a1c;A2=Ta1234*a2c;A3=Ta1234*a3c;A4=Ta1234*a4c;[Ta5]=transfEMP(90,0,0,5,0,0,0,0);[Ta6]=transfEMP(110.55,0,0,sqrt(73),0,0,0,0);A5=Ta5*a5c;A6=Ta6*a6c;[Ta7]=[-cosd(69.44) -sind(69.44) cosd(69.44) sind(69.44)]';[Ta8]=[-cosd(138.36) -sind(138.36) cosd(138.36) sind(138.36)]';[Ta9]=[-cosd(151.92) -sind(151.92) cosd(151.92) sind(151.92)]';A7=Ta7*a7c;{_Error using <a href="matlab:helpUtils.errorDocCallback('mtimes')" style="font-weight:bold"> * </a>Inner matrix dimensions must agree.}_ [h]=mequilibrioar(sqrt(73));a7c=h*a7r;[h]=mequilibrioar(sqrt(145));a8c=h*a8r;[h]=mequilibriomp(sqrt(289));a9c=h*a9r;A7=Ta7*a7c;{_Error using <a href="matlab:helpUtils.errorDocCallback('mtimes')" style="font-weight:bold"> * </a>Inner matrix dimensions must agree.}_ A7=Ta7*a7c';A7=Ta7'*a7c;A8=Ta8'*a8c;A9=Ta9'*a9c;{_Error using <a href="matlab:helpUtils.errorDocCallback('mtimes')" style="font-weight:bold"> * </a>Inner matrix dimensions must agree.}_ [h]=mequilibrioar(sqrt(289));a9c=h*a9r;A9=Ta9'*a9c;a123456f=[a1c a2c a3c a4c a5c a6c]; a789f=[a7c a8c a9c];PR=zeros(21,1);% Ahora vamos a encontrar el vector de equilibro help acum [PR]=acum[lee,Ae,PR] [PR]=acum[lee1,A1,PR];

[PR]=acum[lee1,A1,PR]; |{_Error: Unbalanced or unexpected parenthesis or bracket.}_ [PR]=acum(lee1,A1,PR);[PR]=acum(lee2,A2,PR);[PR]=acum(lee3,A3,PR);[PR]=acum(lee4,A4,PR);[PR]=acum(lee5,A5,PR);[PR]=acum(lee6,A6,PR);[PR]=acum(lee7,A7,PR);{_Attempted to access Ae(2); index out of bounds because numel(Ae)=1. Error in <a href="matlab:helpUtils.errorDocCallback('acum', 'C:\Users\Usuario\Desktop\ANALISIS 2013 MATLAB\Examenes\acum.m', 6)" style="font-weight:bold">acum</a> (<a href="matlab: opentoline('C:\Users\Usuario\Desktop\ANALISIS 2013 MATLAB\Examenes\acum.m',6,0)">line 6</a>) PR(ii)=PR(ii)+Ae(i);}_ A7=T7*a7c;{_Error using <a href="matlab:helpUtils.errorDocCallback('mtimes')" style="font-weight:bold"> * </a>Inner matrix dimensions must agree.}_ diary off