土楔(土くさび)重量、重力式擁壁、段差無盛り土、Calculation of Weight of Soil wedge. ※Retaining wall,Smooth Embankment.

Equation
W=\frac{1}{2}\gamma\cdot H^{2}\frac{\cos(\omega-\alpha)\cos(\alpha-\beta)}{\cos^{2}\alpha\sin(\omega-\beta)}

000268

Scilab
268_01268_02
268_03268_04

 -->clear;
-->for www=1:4;
-->    funcprot(0);
-->    equation='W=1/2*Gamma*H^2*cosd(Omega-Alpha)*cosd(Alpha-Beta)/(cosd(Alpha))^2/sind(Omega-Beta)';
-->    H=10;//(m)
-->    Gamma=2;//(tf/m^3)
-->    xset("window",www);clf();xgrid();
-->    graphtitle='Weight of Soil wedge';
-->    n=100;
-->    Omega=linspace(40,60,n);//(°)
-->    Beta=linspace(5,30,n);//(°)
-->    Alpha=linspace(5,30,n);//(°)
-->    if www==1 | www==2 then
-->        xtitle(graphtitle,'ω(°)','β(°)','W(tf/m)',boxed=1);
-->        Alpha=20;//(°)
-->        fnt='[W]=func1(Omega,Beta)';
-->        deff(fnt,equation);
-->        if www==1 then 
-->            R=245;
-->            G=255;
-->            B=250;
-->            color(R,G,B);
-->            idcolor=color(R,G,B);
-->            fplot3d(Omega,Beta,func1,flag=[idcolor , 2 , 4]);
-->        elseif www==2 then
-->            nz=30;
-->            contour(Omega,Beta,func1,nz,flag=[2,2,4]);
-->            ax=gca();
-->            for iii=1:nz;
-->                cnt=ax.children(iii).children;
-->                eee=length(cnt);
-->                if eee==2 then
-->                    cnt(2).font_size=4;
-->                end
-->            end
-->        end
-->    else
-->        xtitle(graphtitle,'ω(°)','α(°)','W(tf/m)',boxed=1);
-->        Beta=20;//°
-->        fnt='[W]=func2(Omega,Alpha)';
-->        deff(fnt,equation);
-->        if www==3 then 
-->            R=216;
-->            G=191;
-->            B=216;
-->            color(R,G,B);
-->            idcolor=color(R,G,B);
-->            fplot3d(Omega,Alpha,func2,flag=[idcolor , 2 , 4]);
-->        elseif www==4 then
-->            contour(Omega,Alpha,func2,30,flag=[2,2,4]); 
-->            ax=gca();
-->            for iii=1:nz;
-->                cnt=ax.children(iii).children;
-->                eee=length(cnt);
-->                if eee==2 then
-->                    cnt(2).font_size=4;
-->                end
-->            end
-->        end
-->    end 
-->    currentaxes=gca();
-->    currentaxes.font_size=6;
-->    currentaxes.title.font_size=6;
-->    currentaxes.x_label.font_size=6;
-->    currentaxes.y_label.font_size=6;
-->    currentaxes.z_label.font_size=6;
-->    currentaxes.tight_limits="on";
-->    currentaxes.cube_scaling="on";
-->    currentaxes.auto_clear="on";
-->    currentaxes.auto_scale="on";
-->end
-->Alpha=20;
-->func1(45,20)
 ans  =
 
    242.8599  
-->func1(30,5)
 ans  =
 
    254.90323  
-->Beta=15;
-->func2(45,20)
 ans  =
 
    204.49293  
-->func2(30,5)
 ans  =
 
    347.49014  

SciNotes

clear;
for www=1:4;
    funcprot(0);
    equation='W=1/2*Gamma*H^2*cosd(Omega-Alpha)*cosd(Alpha-Beta)/(cosd(Alpha))^2/sind(Omega-Beta)';
    H=10;//(m)
    Gamma=2;//(tf/m^3)
    xset("window",www);clf();xgrid();
    graphtitle='Weight of Soil wedge';
    n=100;
    Omega=linspace(40,60,n);//(°)
    Beta=linspace(5,30,n);//(°)
    Alpha=linspace(5,30,n);//(°)
    if www==1 | www==2 then
        xtitle(graphtitle,'ω(°)','β(°)','W(tf/m)',boxed=1);
        Alpha=20;//(°)
        fnt='[W]=func1(Omega,Beta)';
        deff(fnt,equation);
        if www==1 then 
            R=245;
            G=255;
            B=250;
            color(R,G,B);
            idcolor=color(R,G,B);
            fplot3d(Omega,Beta,func1,flag=[idcolor , 2 , 4]);
        elseif www==2 then
            nz=30;
            contour(Omega,Beta,func1,nz,flag=[2,2,4]);
            ax=gca();
            for iii=1:nz;
                cnt=ax.children(iii).children;
                eee=length(cnt);
                if eee==2 then
                    cnt(2).font_size=4;
                end
            end
        end
    else
        xtitle(graphtitle,'ω(°)','α(°)','W(tf/m)',boxed=1);
        Beta=20;//°
        fnt='[W]=func2(Omega,Alpha)';
        deff(fnt,equation);
        if www==3 then 
            R=216;
            G=191;
            B=216;
            color(R,G,B);
            idcolor=color(R,G,B);
            fplot3d(Omega,Alpha,func2,flag=[idcolor , 2 , 4]);
        elseif www==4 then
            contour(Omega,Alpha,func2,30,flag=[2,2,4]); 
            ax=gca();
            for iii=1:nz;
                cnt=ax.children(iii).children;
                eee=length(cnt);
                if eee==2 then
                    cnt(2).font_size=4;
                end
            end
        end
    end 
    currentaxes=gca();
    currentaxes.font_size=6;
    currentaxes.title.font_size=6;
    currentaxes.x_label.font_size=6;
    currentaxes.y_label.font_size=6;
    currentaxes.z_label.font_size=6;
    currentaxes.tight_limits="on";
    currentaxes.cube_scaling="on";
    currentaxes.auto_clear="on";
    currentaxes.auto_scale="on";
end
Alpha=20;
func1(45,20)
func1(30,5)
Beta=15;
func2(45,20)
func2(30,5)

参考文献
右城猛(1989).『新道路土工指針による擁壁の設計法と計算例』.理工図書.226pp.

アプリケーション
URL Scilab http://www.scilab.org/