Lampiran A. KODE PROGRAM JOZLOPE % membuat window win1=figure('units','points','position',[100 100 600 400],'color',[1 1 1],'menubar','none','resize','off',... 'numbertitle','on','name','JOZLOPE : Grid Search'); menu1=uimenu('parent', win1,'label','File'); menu1_1=uimenu('parent', menu1,'label','Save Figure as JPEG','callback','savemenu'); zoomm=uimenu('parent',menu1,'label','Zoom','callback','zoommenu'); panm=uimenu('parent',menu1,'label','Drag','callback','panmenu'); frame1=uicontrol('parent',win1,'units','points','position',[0 360 116 100],'backgroundcolor',[0.1 0 0.4],... 'style', 'Frame'); frame2=uicontrol('parent',win1,'units','points','position',[0 0 600 60],'backgroundcolor',[0.1 0 0.4],... 'style','Frame'); frame3=uicontrol('parent',win1,'units','points','position',[110 60 2 300],'backgroundcolor',[0.1 0 0.4],... 'style', 'Frame'); frame4=uicontrol('parent',win1,'units','points','position',[114 60 2 300],'backgroundcolor',[0.1 0 0.4],... 'style', 'Frame'); frame5=uicontrol('parent',win1,'units','points','position',[580 60 2 400],'backgroundcolor',[0.1 0 0.4],... 'style','Frame'); frame6=uicontrol('parent',win1,'units','points','position',[585 60 2 400],'backgroundcolor',[0.1 0 0.4],... 'style','Frame'); frame7=uicontrol('parent',win1,'units','points','position',[0 124 110 2],'backgroundcolor',[0.1 0 0.4],... 'style', 'Frame'); frame8=uicontrol('parent',win1,'units','points','position',[0 120 110 2],'backgroundcolor',[0.1 0 0.4],... 'style','Frame'); frame9=uicontrol('parent',win1,'units','points','position',[116 390 490 2],'backgroundcolor',[0.1 0 0.4],... 'style', 'Frame'); frame10=uicontrol('parent',win1,'units','points','position',[116 386 490 2],'backgroundcolor',[0.1 0 0.4],... 'style','Frame'); %Judul-judul INL=uicontrol('parent',win1,'units','points','position',[30 360 65 20],'backgroundcolor',[0.1 0 0.4],... 'style','Text','horizontalalignment','center','string','INPUTS','fontname','arial',... 'fontsize',14,'fontweight','bold','foregroundcolor',[1 1 1]); CL2=uicontrol('parent',win1,'units','points','position',[10 320 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','c (KN/m^2)','fontname','arial',... 'fontsize',10); HL4=uicontrol('parent',win1,'units','points','position',[10 280 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','H (m) ','fontname','arial',... 'fontsize',10); LebarL5=uicontrol('parent',win1,'units','points','position',[10 260 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','L (m) ','fontname','arial',...
73
'fontsize',10); LebarSAL6=uicontrol('parent',win1,'units','points','position',[10 240 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','b (m)','fontname','arial',... 'fontsize',10); TinggiSAL7=uicontrol('parent',win1,'units','points','position',[10 220 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','h (m)','fontname','arial',... 'fontsize',10); xcL8=uicontrol('parent',win1,'units','points','position',[10 200 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','x-start','fontname','arial',... 'fontsize',10); xpL9=uicontrol('parent',win1,'units','points','position',[10 180 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','xp','fontname','arial',... 'fontsize',10); D10=uicontrol('parent',win1,'units','points','position',[10 160 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','D (m)','fontname','arial',... 'fontsize',10); eo11=uicontrol('parent',win1,'units','points','position',[10 140 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','right','string','Incr.(m)','fontname','arial',... 'fontsize',10); %label output TimeL=uicontrol('parent',win1,'units','points','position',[0 40 60 10],'backgroundcolor',[0 0 0],... 'style','Text','string','Time Process :','fontname','arial',... 'fontsize',8,'foregroundcolor',[1 1 1]); Time2L=uicontrol('parent',win1,'units','points','position',[85 40 45 10],'backgroundcolor',[0 0 0],... 'style','Text','string','Minutes','fontname','arial',... 'fontsize',8,'foregroundcolor',[1 1 1]); TanggalL=uicontrol('parent',win1,'units','points','position',[0 5 30 14],'backgroundcolor',[0 0 0],... 'style','Text','string','Made:','fontname','arial',... 'fontsize',12,'foregroundcolor',[1 1 1]); trialL=uicontrol('parent',win1,'units','points','position',[150 40 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','N','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); xlowL=uicontrol('parent',win1,'units','points','position',[205 40 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','XL','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); xtopL=uicontrol('parent',win1,'units','points','position',[260 40 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','XH','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); MFSTL=uicontrol('parent',win1,'units','points','position',[320 40 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','FS MIN','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]);
74
xcenL=uicontrol('parent',win1,'units','points','position',[380 40 30 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','Xc','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); ycenL=uicontrol('parent',win1,'units','points','position',[415 40 30 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','Yc','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); RadL=uicontrol('parent',win1,'units','points','position',[450 40 50 15],'backgroundcolor',[0.9 0.9 0.9],... 'style','Text','horizontalalignment','center','string','R','fontname','arial',... 'fontsize',12,'fontweight','bold','foregroundcolor',[0 0 0]); %input edit1=uicontrol('parent',win1,'units','points','position',[65 340 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit2=uicontrol('parent',win1,'units','points','position',[65 320 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit3=uicontrol('parent',win1,'units','points','position',[65 300 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit4=uicontrol('parent',win1,'units','points','position',[65 280 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit5=uicontrol('parent',win1,'units','points','position',[65 260 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit6=uicontrol('parent',win1,'units','points','position',[65 240 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit7=uicontrol('parent',win1,'units','points','position',[65 220 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit8=uicontrol('parent',win1,'units','points','position',[65 200 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit9=uicontrol('parent',win1,'units','points','position',[65 180 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit10=uicontrol('parent',win1,'units','points','position',[65 160 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); edit11=uicontrol('parent',win1,'units','points','position',[65 140 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10); %edit output grafik1=axes('parent',win1,'units','points','position',[155 95 420 280],'xgrid','on',... 'ygrid','on','xcolor',[0 0 1],'ycolor',[0 0 1],'fontsize',8,'color',[1 1 1]); Time=uicontrol('parent',win1,'units','points','position',[60 40 30 10],'backgroundcolor',[0 0 0],... 'style','Edit','string','0','fontname','arial','fontsize',8,'foregroundcolor',[1 1 1]); Tanggal=uicontrol('parent',win1,'units','points','position',[35 5 100 14],'backgroundcolor',[0 0 0],... 'style','Edit','string','0','fontname','arial','fontsize',12,'foregroundcolor',[1 1 1]); trial=uicontrol('parent',win1,'units','points','position',[150 20 50 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]);
75
xlow=uicontrol('parent',win1,'units','points','position',[205 20 50 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); xtop=uicontrol('parent',win1,'units','points','position',[260 20 50 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); MFST=uicontrol('parent',win1,'units','points','position',[320 20 50 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); xcen=uicontrol('parent',win1,'units','points','position',[380 20 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); ycen=uicontrol('parent',win1,'units','points','position',[415 20 30 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); Rad=uicontrol('parent',win1,'units','points','position',[450 20 50 15],'backgroundcolor',[1 1 1],... 'style','Edit','string','0','fontname','arial','fontsize',10,'foregroundcolor',[1 0 0]); %tombol-tombol tomproses2=uicontrol('parent',win1,'units','points','position',[30 100 65 15],'backgroundcolor',[1 1 1],... 'style','pushbutton','callback','gs3','string','Preview','fontname','arial','fontsize',10); tomproses1=uicontrol('parent',win1,'units','points','position',[30 70 65 15],'backgroundcolor',[1 1 1],... 'style','pushbutton','callback','gs2','string','Compute','fontname','arial','fontsize',10); %logo img2=importdata('boy.jpg'); pbh2 = uicontrol('parent',win1,'Style','pushbutton',... 'Position',[700 0 82 77],... 'CData',img2); %label bj img3=importdata('bj2.JPG'); pbh3 = uicontrol('parent',win1,'Style','pushbutton',... 'Position',[14 450 67 26],... 'CData',img3); %label sudut geser img4=importdata('FA.JPG'); pbh4 = uicontrol('parent',win1,'Style','pushbutton',... 'Position',[14 400 66 24],... 'CData',img4); %PROGRAM MULAI tic %A-INPUT bJ=str2num(get(edit1,'string')); %mengambil data input berat jenis ch=str2num(get(edit2,'string'));% mengambil data input kohesi SDG=str2num(get(edit3,'string'));%mengambil data input sudut geser tanah H=str2num(get(edit4,'string'));%mengambil data input tinggi lereng L=str2num(get(edit5,'string'));% mengambil data input bentang lereng b=str2num(get(edit6,'string'));%mengambil data lebar area grid
76
h=str2num(get(edit7,'string'));%mengambil data tinggi area grid xc1=str2num(get(edit8,'string'));%mengambil data input absis awal area grid xp=str2num(get(edit9,'string'));%mengambil data input ab. awal area masuk D=str2num(get(edit10,'string'));%mengambil data input kedalaman tanah keras eo=str2num(get(edit11,'string'));%mengambil data input spasi antar grid %setelah input, dilakukan pengumpulan nilai FS, A, B, dan a untuk kemudian%dioptimasi %B-inisiasi yc1=H; xc2=xc1+b;%menent. ab. akhir area grid yc2=yc1+h;%menent. ord. akhir area grid xF=0;%inisiasi indeks baris matriks A, B, a, dan FS mo=atan(H/L)*180/pi;%menent. sudut kemiringan lereng xr=L;%menent. ab. akhir garis kemiringan lereng xe=0;%menent. ab. akhir area masuk ling. longsor potensial dx=1;%menent. spasi atau jarak antar pias nF=(fix(b/eo)+1)*(fix(h/eo)+1)*(-xp+1);%menentukan jumlah indeks kolom%matriks A,B,F %C-OPTIMASI %C1--pesan tempat buat FS dan koordinat FSI=zeros(nF,1);%membuat matriks nF x 1 untuk FS ASI=zeros(nF,1);%membuat matriks nF x 1 untuk A BSI=zeros(nF,1);%membuat matriks nF x 1 untuk B aSI=zeros(nF,1);%membuat matriks nF x 1 untuk a %C1a --- ITERASI wait=waitbar(0,'Please Wait...');%menampilkan waitbar for A=xc1:eo:xc2 %iterasi A for B=yc1:eo:yc2 %iterasi B for a=xp:dx:xe%iterasi a xF=xF+1;%menent. indeks matriks FS, A, B, a waitbar(xF/nF); %menent. lama proses iterasi pd waitbar %C1b --- SYARAT TIPE KERUNTUHAN R=sqrt(B^2+(a-A)^2);% menent. jari-jari ling. Y4O=B-(sqrt(R^2-(A)^2));%menent. ord. lingkaran di absis 0 Y41=B-(sqrt(R^2-(0.00001-A)^2));%menent. ord. ling. di ab. 0,00001 Y2O=(0.00001-xr)*tan(mo*pi/180)+H;%menent. ord. garis miring di ab. 0,00001 Y4D=B-R;%menent. ord. ling. di ab. A xH=A+(sqrt(R^2-(H-B)^2));%menent. ab. akhir pias if Y4O<=0 & Y41<=Y2O & Y4D>=-D & xH>=xr %seleksi tipe ling. %C1c ---REGENERASI FAKTOR AMAN m=fix(xH-1);%menent. m
77
n=m-a+1;%menent. n AI=zeros(1,n);% menent. matriks luasan 1 x n XS=zeros(1,n);% menent. matriks titik berat 1 x n bi=zeros(1,n);%menent. matriks busur lingk. 1 x n sdti=zeros(1,n);%menent. matriks sudut pias 1 x n Wi=zeros(1,n);%menent. matriks berat pias 1 x n CH=zeros(1,n);% menent. matriks gaya kohesi 1 x n mcos=zeros(1,n);%menent. bagian cos matriks persamaan sudut bishop 1 x n mtan=zeros(1,n);%menent. bagian tan matriks persamaan sudut bishop 1 x n TR=zeros(1,n);% menent. matriks momen penahan 1 x n TD=zeros(1,n);% menent. matriks momen pengguling 1 x n zi=0;%inisasi nomor indeks kolom matriks AI dkk for xi1=a:dx:m %iterasi titik masuk ling. longsor potensial xia=xi1;%menent. ab. awal pias xib=xia+dx;%menent. ab. akhir pias zi=zi+1;%menentukan nomor indeks kolom matriks AI dkk if xia<0 %seleksi untuk pias Tipe I Ai=-xib*B-1/4*(-2*xib+2*A)*(R^2-xib^2+2*xib*A-A^2)^(1/2)... +1/2*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^(1/2))... +xia*B+1/4*(-2*xia+2*A)*(R^2-xia^2+2*xia*A-A^2)^(1/2)... -1/2*R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^(1/2)); xsdti=(-1/3*(R^2-xib^2+2*xib*A-A^2)^(3/2)... +1/2*A*xib*(R^2-xib^2+2*xib*A-A^2)^(1/2)-1/2*... (R^2-xib^2+2*xib*A-A^2)^(1/2)*A^2+1/2*A*R^2*atan((xib-A)/... (R^2-xib^2+2*xib*A-A^2)^(1/2))-1/2*xib^2*B+1/3*... (R^2-xia^2+2*xia*A-A^2)^(3/2)-1/2*A*xia*... (R^2-xia^2+2*xia*A-A^2)^(1/2)+1/2*(R^2-xia^2+2*xia*A-A^2)^... (1/2)*A^2-1/2*A*R^2*atan((xia-A)/... (R^2-xia^2+2*xia*A-A^2)^(1/2))+1/2*xia^2*B)/Ai; elseif xia>=0 if xia<xr%seleksi untuk pias Tipe IIa if xia<m if xib<=xr Ai=1/2*tan(1/180*mo*pi)*xib^2-... tan(1/180*mo*pi)*xr*xib+H*xib-xib*B-1/4*(-2*xib+2*A)*(R^2-xib^2+2*xib*A-A^2)^... (1/2)+1/2*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^... (1/2))-1/2*tan(1/180*mo*pi)*xia^2+tan(1/180*mo*pi)*... xr*xia-H*xia+xia*B+1/4*(-2*xia+2*A)*R^2-xia^2+2*xia*A-A^2)^(1/2)-1/2*... R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^(1/2)); xsdti=(-1/3*(R^2-xib^2+2*xib*A-A^2)^... (3/2)+1/2*A*xib*(R^2-xib^2+2*xib*A-A^2)^...
78
(1/2)-1/2*(R^2-xib^2+2*xib*A-A^2)^(1/2)*... A^2+1/2*A*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^... (1/2))+1/3*tan(1/180*mo*pi)*xib^3-1/2*xib^2*... tan(1/180*mo*pi)*xr+1/2*xib^2*H-1/2*xib^2*B+1/3*... (R^2-xia^2+2*xia*A-A^2)^(3/2)-1/2*A*xia*... (R^2-xia^2+2*xia*A-A^2)^(1/2)+1/2*(R^2-xia^2+2*xia*A-A^2)^... (1/2)*A^2-1/2*A*R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^... (1/2))-1/3*tan(1/180*mo*pi)*xia^3+1/2*xia^2*tan(1/180*mo*pi)*... xr-1/2*xia^2*H+1/2*xia^2*B)/Ai; elseif xib>xr%seleksi untuk pias Tipe IIb xib=xr; Ai=1/2*tan(1/180*mo*pi)*xib^2-tan(1/180*mo*pi)*... xr*xib+H*xib-xib*B-1/4*(-2*xib+2*A)*... (R^2-xib^2+2*xib*A-A^2)^(1/2)+1/2*R^2*atan((xib-A)/... (R^2-xib^2+2*xib*A-A^2)^(1/2))-1/2*tan(1/180*mo*pi)*xia^2+... tan(1/180*mo*pi)*xr*xia-H*xia+xia*B+1/4*(-2*xia+2*A)*... (R^2-xia^2+2*xia*A-A^2)^(1/2)-1/2*R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^(1/2)); xsdti=(-1/3*(R^2-xib^2+2*xib*A-A^2)^... (3/2)+1/2*A*xib*(R^2-xib^2+2*xib*A-A^2)^... (1/2)-1/2*(R^2-xib^2+2*xib*A-A^2)^(1/2)*A^2+1/2*A*R^2*... atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^(1/2))+1/3*... tan(1/180*mo*pi)*xib^3-1/2*xib^2*tan(1/180*mo*pi)*... xr+1/2*xib^2*H-1/2*xib^2*B+1/3*(R^2-xia^2+2*xia*A-A^2)^... (3/2)-1/2*A*xia*(R^2-xia^2+2*xia*A-A^2)^(1/2)+1/2*... (R^2-xia^2+2*xia*A-A^2)^(1/2)*A^2-1/2*A*R^2*atan((xia-A)/... (R^2-xia^2+2*xia*A-A^2)^(1/2))-1/3*tan(1/180*mo*pi)*... xia^3+1/2*xia^2*tan(1/180*mo*pi)*xr-1/2*xia^2*H+1/2*xia^2*B)/Ai; end elseif xia==m%seleksi untuk pias Tipe IIc xia1=xia; xib1=xr; xia2=xr; xib2=xH; Ai11=1/2*tan(1/180*mo*pi)*xib1^2-tan(1/180*mo*pi)*... xr*xib1+H*xib1-xib1*B-1/4*(-2*xib1+2*A)*... (R^2-xib1^2+2*xib1*A-A^2)^(1/2)+1/2*R^2*atan((xib1-A)/... (R^2-xib1^2+2*xib1*A-A^2)^(1/2))-1/2*tan(1/180*mo*pi)*... xia1^2+tan(1/180*mo*pi)*xr*xia1-H*xia1+xia1*B+1/4*... (-2*xia1+2*A)*(R^2-xia1^2+2*xia1*A-A^2)^... (1/2)-1/2*R^2*atan((xia1-A)/(R^2-xia1^2+2*xia1*A-A^2)^(1/2)); Ai12=H*xib2-xib2*B-1/4*(-2*xib2+2*A)*...
79
(R^2-xib2^2+2*xib2*A-A^2)^(1/2)+1/2*R^2*atan((xib2-A)/(R^2-xib2^2+2*xib2*A-A^2)^... (1/2))-H*xia2+xia2*B+1/4*(-2*xia2+2*A)*(R^2-xia2^2+2*xia2*A-A^2)^(1/2)-1/2*R^2*... atan((xia2-A)/(R^2-xia2^2+2*xia2*A-A^2)^(1/2)); xsdti11=(-1/3*(R^2-xib1^2+2*xib1*A-A^2)^(3/2)+1/2*... A*xib1*(R^2-xib1^2+2*xib1*A-A^2)^(1/2)-1/2*... (R^2-xib1^2+2*xib1*A-A^2)^(1/2)*A^2+1/2*A*... R^2*atan((xib1-A)/(R^2-xib1^2+2*xib1*A-A^2)^... (1/2))+1/3*tan(1/180*mo*pi)*xib1^3-1/2*xib1^2*... tan(1/180*mo*pi)*xr+1/2*xib1^2*H-1/2*xib1^2*B+1/3*... (R^2-xia1^2+2*xia1*A-A^2)^(3/2)-1/2*A*xia1*(R^2-xia1^2+2*xia1*A-A^2)^(1/2)+1/2*... (R^2-xia1^2+2*xia1*A-A^2)^(1/2)*A^2-1/2*A*R^2*atan((xia1-A)/(R^2-xia1^2+2*xia1*A-A^2)^... (1/2))-1/3*tan(1/180*mo*pi)*xia1^3+1/2*xia1^2*... tan(1/180*mo*pi)*xr-1/2*xia1^2*H+1/2*xia1^2*B)/Ai11; xsdti12=(-1/3*(R^2-xib2^2+2*xib2*A-A^2)^... (3/2)+1/2*A*xib2*(R^2-xib2^2+2*xib2*A-A^2)^... (1/2)-1/2*(R^2-xib2^2+2*xib2*A-A^2)^(1/2)*... A^2+1/2*A*R^2*atan((xib2-A)/(R^2-xib2^2+2*xib2*A-A^2)^... (1/2))+1/2*xib2^2*H-1/2*xib2^2*B+1/3*... (R^2-xia2^2+2*xia2*A-A^2)^(3/2)-1/2*A*xia2*... (R^2-xia2^2+2*xia2*A-A^2)^(1/2)+1/2*(R^2-xia2^2+2*xia2*A-A^2)^... (1/2)*A^2-1/2*A*R^2*atan((xia2-A)/(R^2-xia2^2+2*xia2*A-A^2)^... (1/2))-1/2*xia2^2*H+1/2*xia2^2*B)/Ai12; Ai=Ai11+Ai12; xsdti=(xsdti11*Ai11+xsdti12*Ai12)/(Ai); end elseif xia>=xr%seleksi untuk pias Tipe IIIa if xia<m if (xia-xr)
80
elseif (xia-xr)>=dx%seleksi untuk pias Tipe IIIb Ai=H*xib-xib*B-1/4*(-2*xib+2*A)*(R^2-xib^2+2*xib*A-A^2)^... (1/2)+1/2*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^... (1/2))-H*xia+xia*B+1/4*(-2*xia+2*A)*(R^2-xia^2+2*xia*A-A^2)^... (1/2)-1/2*R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^(1/2)); xsdti=(-1/3*(R^2-xib^2+2*xib*A-A^2)^(3/2)+1/2*A*xib*... (R^2-xib^2+2*xib*A-A^2)^(1/2)-1/2*(R^2-xib^2+2*xib*A-A^2)^... (1/2)*A^2+1/2*A*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^... (1/2))+1/2*xib^2*H-1/2*xib^2*B+1/3*(R^2-xia^2+2*xia*A-A^2)^... (3/2)-1/2*A*xia*(R^2-xia^2+2*xia*A-A^2)^(1/2)+1/2*... (R^2-xia^2+2*xia*A-A^2)^(1/2)*A^2-1/2*A*R^2*atan((xia-A)/... (R^2-xia^2+2*xia*A-A^2)^(1/2))-1/2*xia^2*H+1/2*xia^2*B)/Ai; end elseif xia==m if (xia-xr)=dx%seleksi untuk pias Tipe IIId xib=xH; Ai=H*xib-xib*B-1/4*(-2*xib+2*A)*(R^2-xib^2+2*xib*A-A^2)^... (1/2)+1/2*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^(1/2))... -H*xia+xia*B+1/4*(-2*xia+2*A)*(R^2-xia^2+2*xia*A-A^2)^(1/2)... -1/2*R^2*atan((xia-A)/(R^2-xia^2+2*xia*A-A^2)^(1/2)); xsdti=(-1/3*(R^2-xib^2+2*xib*A-A^2)^(3/2)+1/2*A*xib*... (R^2-xib^2+2*xib*A-A^2)^(1/2)-1/2*(R^2-xib^2+2*xib*A-A^2)^... (1/2)*A^2+1/2*A*R^2*atan((xib-A)/(R^2-xib^2+2*xib*A-A^2)^... (1/2))+1/2*xib^2*H-1/2*xib^2*B+1/3*(R^2-xia^2+2*xia*A-A^2)^... (3/2)-1/2*A*xia*(R^2-xia^2+2*xia*A-A^2)^(1/2)+1/2*... (R^2-xia^2+2*xia*A-A^2)^(1/2)*A^2-1/2*A*R^2*atan((xia-A)/... (R^2-xia^2+2*xia*A-A^2)^(1/2))-1/2*xia^2*H+1/2*xia^2*B)/Ai;
81
end
end
end end ri=sqrt((xib-xia)^2+(-(R^2-xib^2+2*xib*A-A^2)^(1/2)+(R^2-xia^2+2*xia*A-A^2)^(1/2))^2);%menent. jari-jari tembereng sdARC=asin(ri/(2*R))*180/pi;%menent. sudut juring AI(1,zi)=Ai;%menent. luas pias ke-zi XS(1,zi)=(xsdti);%menent.titik berat pias ke-zi bi(1,zi)=2*sdARC*(2*pi*R)/360;%menent.busur pias ke-zi sdti(1,zi)=atan(-1/2/(R^2-XS(1,zi)^2+2*XS(1,zi)*A-A^2)^(1/2)*(-2*XS(1,zi)+2*A));%menent. sudut pias ke-zi Wi(1,zi) =AI(1,zi)*bJ;%menent. berat pias ke-zi CH(1,zi)=ch*bi(1,zi);%menent. gaya kohesi pias ke-zi mcos(1,zi)=cos(sdti(1,zi));%menent. bagian cos persamaan koefisien sudut bishop pias ke-zi mtan(1,zi)=tan(SDG*pi/180)*sin(sdti(1,zi));%menent. bagian tan persamaan koefisien sudut bishop pias zi TR(1,zi)= CH(1,zi)*cos(sdti(1,zi))+Wi(1,zi)*tan(SDG*pi/180); %menent. momen penahan pias ke-zi TD(1,zi)= Wi(1,zi)*sin(sdti(1,zi)); %menent. momen pengguling pias ke-zi end FS=1;%menentukan FS sumTD=sum(TD);%menjumlahkan seluruh momen pengguling for num=1:10 %iterasi FS sepuluh kali FS1=FS; mi=1./(mcos+mtan/FS1);%substitusi nilai FS ke mi TR2=TR.*mi;%menent. momen penahan sejati pias FS=sum(TR2)/sumTD;%menent. FS end %C1d --- PENGELOMPOKAN FAKTOR AMAN FSI(xF,1)=real(FS);%mengumpulkan ASI(xF,1)=A;%mengumpulkan A tiap BSI(xF,1)=B;%mengumpulkan B tiap aSI(xF,1)=a;%mengumpulkan a tiap end end
FS real tiap ling. ling. potensial ling. potensial ling. potensial
end end close(wait);%menutup waitbar %C2-- SELEKSI POPULASI UNTUK FAKTOR AMAN MINIMUM minFSI=min(nonzeros(FSI));%optimasi FS minaSI=max(((aSI-FSI)==(aSI-minFSI))*1000+aSI)-1000;%menent. a min
82
minASI=max(((ASI-FSI)==(ASI-minFSI))*1000+ASI)-1000;%menent. A min minBSI=max(((BSI-FSI)==(BSI-minFSI))*1000+BSI)-1000;%menent. B min %D - OUTPUT %D1 -- OUTPUT GAMBAR LINGKARAN KELONGSORAN KRITIS a=minaSI; A=minASI; B=minBSI; R=sqrt(B^2+(a-A)^2); xH=A+(sqrt(R^2-(H-B)^2)); m=fix(xH-1); xL=a:0.001:xH; YL=B-sqrt(R^2-(xL-A).^2); x5=a:dx:m; Y5=B-sqrt(R^2-(x5-A).^2); t4=[a 0;A B]; t5=[A B;xH H]; %D2 -- OUTPUT GRIDSEARCH t9=[ASI BSI]; hold on axis equal grid on set(win1,'CurrentAxes',grafik1); plot(xL,YL,'-r',A,B,'or',x5,Y5,'or',t4(:,1),t4(:,2),'-ok',t5(:,1),t5(:,2),'-ok',t9(:,1),t9(:,2),'xr'); %D3 -- OUTPUT PARAMETER FAKTOR AMAN set(trial,'String',num2str(nF)); set(xlow,'String',num2str(a)); set(xtop,'String',num2str(xH)); set(MFST,'String',num2str(minFSI)); set(xcen,'String',num2str(A)); set(ycen,'String',num2str(B)); set(Rad,'String',num2str(R)); %D4 -- OUTPUT WAKTU EKSESKUSI PROGRAM t=toc/60; set(Time,'string',num2str(t)); set(Tanggal,'string',datestr(now)); %PROGRAM SELESAI %Membuat "preview" H=str2num(get(edit4,'string')); L=str2num(get(edit5,'string')); b=str2num(get(edit6,'string')); h=str2num(get(edit7,'string'));
83
xc1=str2num(get(edit8,'string')); xp=str2num(get(edit9,'string')); D=str2num(get(edit10,'string')); eo=str2num(get(edit11,'string')); yc1=H; xc2=xc1+b; yc2=yc1+h; R=sqrt(xc2^2+yc2^2); xr=L; a=xp; xt=xc2+R+10; if xr<=50 xt2=xt; else xt2=xt+yc2; end %gambar space area t1=[xc1 yc1;xc2 yc1;xc1 yc1;xc1 yc2]; t2=[xc2 yc1;xc2 yc2;xc1 yc2;xc2 yc2]; t11a=a:1:0; t11b=zeros(1,-a+1); %gambar geometri lereng %gambar garis miring m=atan(H/L)*180/pi; x2=0:0.1:xr; Y2 =(x2-xr).*tan(m*pi/180)+H; %gambar permukaan tanah atas dan bawah t3=[-xt2 0;0 0;xr H;xt2 H]; t8=[-xt2 -D;xt2 -D]; t9=[-xt2 0;0 0;xr H;xt2 H;xt2 -D;-xt2 -D]; %Box Properties hold on axis equal grid on xlabel('Distance (m)'); ylabel('Elevation (m)'); axis ([-200 200 -70 200]); fill(t9(:,1),t9(:,2),'y'); plot(x2,Y2,'-k',t1(:,1),t1(:,2),'-r',t2(:,1),t2(:,2),'-r',t3(:,1),t3(:,2),... '-k',t8(:,1),t8(:,2),'-k',t11a,t11b,'xr'); %selesai
84