1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
//example 9.5
clc; funcprot(0);
Mmax=43.72;
sp=["PSA-31","PSA-23"];//sheet pile
H=[7.9,7.9]
I=[4.41e-6,5.63e-6];
p=[0.00466,0.00365];
S=[10.8e-5,12.8e-5];
Md=[18.63,22.08];
printf("SheetPile I(m^4/m) H(m) p\t Logp S(m^3/m) Md(kN-m/m) Md/Mmax \n")
for i=1:2
Logp(i)=log10(p(i));
k(i)=Md(i)/Mmax;
printf("%s\t %.2e %.2f %f %.2f %.2e %.2f\t %.3f\n",sp(i),I(i),H(i),p(i),Logp(i),S(i),Md(i),k(i));
end
plot(Logp,k);
xtitle("","LogP","Md/Mmax");
|
