Lampiran A

  • Uploaded by: Raga Candradimuka
  • 0
  • 0
  • May 2020
  • PDF

This document was uploaded by user and they confirmed that they have the permission to share it. If you are author or own the copyright of this book, please report to us by using this DMCA report form. Report DMCA


Overview

Download & View Lampiran A as PDF for free.

More details

  • Words: 1,385
  • Pages: 12
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

Related Documents

Lampiran A
May 2020 14
Lampiran A
December 2019 23
1108 Lampiran A
June 2020 2

More Documents from ""