土楔(土くさび)重量、重力式擁壁、段差有盛り土、CALCULATION OF WEIGHT OF SOIL WEDGE. ※RETAINING WALL,UNEVEN EMBANKMENT.

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

000269

Scilab
269sc01269sc02
269sc03269sc04

-->clear;
-->for www=1:4
-->    funcprot(0);
-->    equation='W=1/2*Gamma*((H+H1)^2*cosd(omega-alpha)/sind(omega)-H1^2*cosd(alpha-Beta)/sind(Beta))/cosd(alpha)';
-->    H1=1;//(m)
-->    H=10;//(m)
-->    Gamma=2;//(tf/m^3)
-->    xset("window",www);clf();xgrid();
-->    graphtitle='Weight of Soil wedge';
-->    if www==1 | www==2 then
-->        xtitle(graphtitle,'ω(°)','β(°)','W(tf/m)',boxed=1);
-->        alpha=20;//(°)
-->        fnt='[W]=func1(omega,Beta)';
-->        if www==1 then 
-->            omega=linspace(40,60);//(°)
-->            Beta=linspace(5,30);//(°)
-->            R=245;
-->            G=255;
-->            B=250;
-->            color(R,G,B);
-->            idcolor=color(R,G,B);
-->            deff(fnt,equation);
-->            fplot3d(omega,Beta,func1,flag=[idcolor , 2 , 4]);
-->        elseif www==2 then
-->            deff(fnt,equation);
-->            contour(omega,Beta,func1,30,flag=[2,2,4]); 
-->        end
-->    else
-->        xtitle(graphtitle,'ω(°)','α(°)','W(tf/m)',boxed=1);
-->        Beta=20;//°
-->        fnt='[W]=func2(omega,alpha)';
-->        if www==3 then 
-->            omega=linspace(40,60);//(°)
-->            alpha=linspace(5,30);//(°)
-->            R=216;
-->            G=191;
-->            B=216;
-->            color(R,G,B);
-->            idcolor=color(R,G,B);
-->            deff(fnt,equation);
-->            fplot3d(omega,alpha,func2,flag=[idcolor , 2 , 4]);
-->        elseif www==4 then
-->            deff(fnt,equation);
-->            contour(omega,alpha,func2,30,flag=[2,2,4]); 
-->        end
-->    end 
-->    currentaxes=get("current_axes");
-->    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  =
 
    161.92895  
-->func1(30,5)
 ans  =
 
    241.82452  
-->Beta=20;
-->func2(45,20)
 ans  =
 
    161.92895  
-->func2(30,5)
 ans  =
 
    217.32931  

SciNotes

clear;
for www=1:4
    funcprot(0);
    equation='W=1/2*Gamma*((H+H1)^2*cosd(omega-alpha)/sind(omega)-H1^2*cosd(alpha-Beta)/sind(Beta))/cosd(alpha)';
    H1=1;//(m)
    H=10;//(m)
    Gamma=2;//(tf/m^3)
    xset("window",www);clf();xgrid();
    graphtitle='Weight of Soil wedge';
    if www==1 | www==2 then
        xtitle(graphtitle,'ω(°)','β(°)','W(tf/m)',boxed=1);
        alpha=20;//(°)
        fnt='[W]=func1(omega,Beta)';
        if www==1 then 
            omega=linspace(40,60);//(°)
            Beta=linspace(5,30);//(°)
            R=245;
            G=255;
            B=250;
            color(R,G,B);
            idcolor=color(R,G,B);
            deff(fnt,equation);
            fplot3d(omega,Beta,func1,flag=[idcolor , 2 , 4]);
        elseif www==2 then
            deff(fnt,equation);
            contour(omega,Beta,func1,30,flag=[2,2,4]); 
        end
    else
        xtitle(graphtitle,'ω(°)','α(°)','W(tf/m)',boxed=1);
        Beta=20;//°
        fnt='[W]=func2(omega,alpha)';
        if www==3 then 
            omega=linspace(40,60);//(°)
            alpha=linspace(5,30);//(°)
            R=216;
            G=191;
            B=216;
            color(R,G,B);
            idcolor=color(R,G,B);
            deff(fnt,equation);
            fplot3d(omega,alpha,func2,flag=[idcolor , 2 , 4]);
        elseif www==4 then
            deff(fnt,equation);
            contour(omega,alpha,func2,30,flag=[2,2,4]); 
        end
    end 
    currentaxes=get("current_axes");
    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=20;
func2(45,20)
func2(30,5)

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

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