%CÁLCULOMECÁNICO %TABLADETENDIDO %ALTURADELAINSTALACIÓN321-334(m) %ZONAA %LA_56 %VANO[637060]: clear alfa=1.91*1e-5;%Coeficientededilatación(ºC-1) D=9.45;%Diámetroaparenteen(mm) S=54.6;%Secciónen(mm²) P=0. P=0.18 186; 6; %P %Pes eso opr prop opio iou uni nita tari rio o(d (daN aN/m /m) ) Q=1629;%Cargaderotura(daN) E=7900; %M %Módulodeelasticidaden ódulodeelasticidaden(daN/mm²) (daN/mm²) t1=-5;%Temperaturainicial P1=P;%Enelvanoderegulación n=5;%C.S.minimo T1=Q/n;%Tensiónmáxima(daN) %SOBRECARGAPORACCIÓNDELVIENTO: ifD>16 K_viento=0.05; else K_viento=0.06; end Pv=K_viento*D;%Accióndelvientoen(daN/m) Pv_p=sqrt(Pv^2+P^2);%Accióndelvien Pv_p=sqrt(Pv^2+P^2); %Accióndelvientoypesopropio(daN/m) toypesopropio(daN/m) fprintf('\n') fprintf('Laresultantedelvientoyelpesopr fprintf('Laresultante delvientoyelpesopropioes%0.3fdaN/m\n',Pv opioes%0.3fdaN/m\n',Pv_p) _p) fprintf('\n') %VANODEREGULACIÓN: v=[637060];%Longituddelosvanosen( v=[637060];%L ongituddelosvanosen(m) m) v_r=sqrt(sum(v.^3)/sum(v));%Vanoid v_r=sqrt(sum(v.^3)/sum (v));%Vanoidealderegulación ealderegulación %v_r=round(v_r); fprintf('Lalongituddelvanoderegulaciónes fprintf('Lalongitudd elvanoderegulaciónesde%0.0fm\n',v_r) de%0.0fm\n',v_r) fprintf('\n') f=v_r^2*Pv_p/8/T1;%FlechadelvanoderegulacionconPv_p+t=-5ºC f=v_r^2*Pv_p/8/T1;%Flechadelvanoderegul acionconPv_p+t=-5ºC fprintf('Laflechaesiguala%0.2fm\n',f) fprintf('\n') %TRACCIÓNMÁXIMA(P+V;t=-5ºC): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=-5; B=(A+alfa*t)*E*S; C=v_r^2*Pv_p^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end
T2=max(resul);%Soluciónreal f=v_r^2*Pv_p/8/T2;%Flechacont=-5ºC fcat=(1+D^2/2/v_r^2)*(v_r^2*Pv_p/8/T2); H=T2/P; fprintf('TracciónMáxima(P+V;t=-5)\n') fprintf('Vanoderegulación(m)TracciónMáxima(daN)Flecha(m)FlechaCat(m)\n') fprintf('%9.0f%25.1f%15.2f%11.2f\n',v_r,T2,f,fcat) %CHS(P;t=-5ºC): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=-5; B=(A+alfa*t)*E*S; C=v_r^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal chs=T2*100/Q; f=v_r^2*P/8/T2;%Flechacont=15ºC f2=T2/P*(((exp(v_r*P/2/T2)+exp(-v_r*P/2/T2))/2)-1); fprintf('CHS(P;t=-5)\n') fprintf('Vanoderegulación(m)Tracción(daN)CHSF(m)Fcat(m)\n') fprintf('%9.0f%20.1f%12.2f%8.2f%7.2f\n',v_r,T2,chs,f,f2) %EDS(P;t=15): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=15; B=(A+alfa*t)*E*S; C=v_r^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal eds=T2*100/Q; f=v_r^2*P/8/T2;%Flechacont=15ºC f2=T2/P*(((exp(v_r*P/2/T2)+exp(-v_r*P/2/T2))/2)-1); fprintf('EDS(P;t=15)\n') fprintf('Vanoderegulación(m)Tracción(daN)EDSF(m)Fcat(m)\n') fprintf('%9.0f%20.1f%12.2f%8.2f%7.2f\n',v_r,T2,eds,f,f2) %FLECHAMÁXIMA(P+V;t=15): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=15; B=(A+alfa*t)*E*S; C=v_r^2*Pv_p^2*E*S/24;
bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal f=v_r^2*Pv_p/8/T2;%Flechacont=15ºC f2=T2/Pv_p*(((exp(v_r*Pv_p/2/T2)+exp(-v_r*Pv_p/2/T2))/2)-1); fprintf('Flechamáxima(P+V;t=15)\n') fprintf('Vanoderegulación(m)Flechamáxima(m)Tracción(daN)Fcat(m)\n') fprintf('%9.0f%25.2f%17.1f%14.2f\n',v_r,f,T2,f2) %FLECHAMÁXIMA(P;t=0): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=0; B=(A+alfa*t)*E*S; C=v_r^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal f=v_r^2*P/8/T2;%Flechacont=0ºC f2=T2/P*(((exp(v_r*P/2/T2)+exp(-v_r*P/2/T2))/2)-1); fprintf('Flechamáxima(P;t=0)\n') fprintf('Vanoderegulación(m)Flechamáxima(m)Tracción(daN)Fcat(m)\n') fprintf('%9.0f%25.2f%17.1f%14.2f\n',v_r,f,T2,f2) %FLECHAMÁXIMA(P;t=50): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S; t=50; B=(A+alfa*t)*E*S; C=v_r^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal f=v_r^2*P/8/T2;%Flechacont=50ºC f2=T2/P*(((exp(v_r*P/2/T2)+exp(-v_r*P/2/T2))/2)-1); fprintf('FlechaMáxima(P;t=50)\n') fprintf('Vanoderegulación(m)Flechamáxima(m)Tracción(daN)Fcat(m)\n') fprintf('%9.0f%25.2f%17.1f%14.2f\n',v_r,f,T2,f2) %FLECHAMÍNIMA(P;t=-5): A=v_r^2*P1^2/24/T1^2-alfa*t1-T1/E/S;
t=-5; B=(A+alfa*t)*E*S; C=v_r^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul);%Soluciónreal f=v_r^2*P/8/T2;%Flechacont=-5ºC f2=T2/P*(((exp(v_r*P/2/T2)+exp(-v_r*P/2/T2))/2)-1); fprintf('TracciónMáxima(P;t=-5)\n') fprintf('Vanoderegulación(m)Flechamáxima(m)Tracción(daN)Fcat(m)\n') fprintf('%9.0f%25.2f%17.1f%14.2f\n',v_r,f,T2,f2) %Tabladetendido: %Temperaturas->[-5:5:50] a=v_r; d=1;%Desnivelentreapoyosen(m) tr=-5:5:50; freg=[]; freg1=[]; freg2=[]; tens=[]; fork=1:length(tr) tf=tr(k); B=(A+alfa*tf)*E*S; C=a^2*P^2*E*S/24; bb=[1,B,0,-C]; solucion=roots(bb); resul=[]; fork=1:3 ifabs(solucion(k))==real(solucion(k)) resul=[resul,solucion(k)]; end end T2=max(resul); tens=[tens,T2]; fr=a^2*P/8/T2; f2=(1+d^2/2/a^2)*(a^2*P/8/T2); fr1=T2/P*(((exp(a*P/2/T2)+exp(-a*P/2/T2))/2)-1); freg=[freg,fr]; freg1=[freg1,f2]; freg2=[freg2,fr1]; end fprintf('Tabladetendido87\n') fprintf('Temperatura(ºC)Tension(kg)Flecha(m)Flecha2(m)Flechac at(m)\n') fork=1:length(tr) fprintf('%8.f%20.1f%11.2f%15.2f%15.2f\n',tr(k),tens(k),freg(k),freg1(k),fre g2(k)) end