function [Wtps,Wstr,Wfuel,Wb,Ww,Wfinh,Wfinv,Wgear,Wthrua,Swfh,Swfv] = pesi_strutturali(Lb,ULF,Dbe,Qmax,Sbtot,Wgtot,WH2,mf,Sref,AR,...
    lambda,tsuc,lambda_mezzi,Ahfp,Avfp,Ttott,Wins,rho_f,rho_th,delta)
    
%calcolo pesi strutturali
    o=abs((Lb*ULF/Dbe)^0.15*Qmax^0.16*Sbtot^1.05);
    Wb=0.341*mf*o;
    Wfuel=Wgtot*WH2;
    Wemp=Wgtot-Wfuel;
    Wtnk=Wfuel/rho_f*rho_th;
    delta=1-delta;
    Ww=0.2958*mf*(((Wemp-Wtnk*delta)*ULF/1000)^0.52*Sref^0.7*AR^0.47...
       *((1+lambda)/tsuc)^0.4*(0.3+0.7/cos(lambda_mezzi)))^1.017;
    Swfh=Ahfp*Sref; %superficie tail orizzontale
    Wfinh=0.0035*(Wgtot/Sref)^0.6*(Swfh)^1.2*Qmax^0.8;
    Swfv=Avfp*Sref; %superficie tail verticale
    Wfinv=5*(Swfv)^1.09;
    Wgear=0.00916*Wgtot^1.124;
    Wthrua=0.00625*Ttott+69;
    Stb=Sbtot/2;

    Wtps=Wins*(Stb+Sref+Swfh);
    Wstr=Wb+Ww+Wfinh+Wfinv+Wtps+Wgear+Wthrua;
end