Program Untuk Menghitung Waktu Layan Beton Bertulang (Matlab)

Wd(1) = 0;                                                                        % definisi vaiabel Wd sebelum terjadi korosi

k = 1.6;                                                                             % konstanta k

c = input('c = ');                                                                % input nilai variabel

Dc = input('Dc = ');                                                           % input nilai variabel

Do2 = input('Do2 = ');                                                       % input nilai variabel

Co2 = input('Co2 = ');                                                       % input nilai variabel

rho = input('rho = ');                                                         % input nilai variabel

phi = input('phi = ');                                                          % input nilai variabel

shi = input('shi = ');                                                           % input nilai variabel

ttz = input('ttz = ');                                                           % input nilai variabel

v = input('v = ');                                                                % input nilai variabel

Db = input('Db = ');                                                            % input nilai variabel

fc = input('fc (28 fc 100) = ');                                            % input nilai variabel

TI = (1/Dc)*(c/k)^2;                                                          % definisi rumus TI

r = (Do2/c)*Co2;                                                                % definisi rumus r

TII = (phi*ttz*v)/(1.5*r);                                                   % definisi rumus TII

alph = ((0.0007*Db - 0.04)*c) + (-0.066*Db + 5.92);      % definisi rumus alpha

beta = (0.0016*fc + 1.0037);                                            % definisi rumus beta

alph0 = (0.003*(c^1.5558));                                             % definisi rumus alpha0

beta0 = (3.8865*(Db^(-0.537)));                                     % definisi rumus beta0

gamma0 = (0.0049*fc) + 0.8539;                                     % definisi rumus gamma0

i=2;                                                                                     % variabel bantu i

Wd(i)=Wd(1)+0.1;                                                             % definisi Wd sesaat sebelum korosi

while (Wd(i)

       Wcorr(i)=(rho/(phi*(shi-1)))*((c+Db)/((3*Db)+(5*c)))*((alph0*beta0*gamma0)+(alph*beta*Wd(i)));

       M(i)=Wcorr(i);

       TIII(i)=(((M(i))*c)/(Do2*Co2))-TII;

       i = i + 1;

       Wd(i)=Wd(i-1) + 0.1;

end

Mplot = M([2:(i-1)]);

Wdplot = Wd([2:(i-1)]);

TIIIplot = TIII([2:(i-1)]);

figure, plot(Mplot, Wdplot);

figure, plot(TIIIplot, Mplot);

referensi : Sudjono, A.S., Prediksi Waktu Layan Bangunan Beton, Civil Engineering Dimension, Vol. 7, No. 1, 6–15, March 2005

---