Civil Engineering Reference
In-Depth Information
% Programme File #3: linstoch.m
% This Matlab programme performs wavelet analysis of a
soil-tank-footing interaction problem
% where the system is subjected to strong ground shaking
(vide Section 4.7.1). Values of main
% parameters are:
% Shear wave velocity: vs = 200 m/s
% Tank heights ('ht') considered: ht = 3, 6 and 10 m
% Height-radius ratio of the tank assumed: hr = 3.0
(slender)
%
clear all;
close all;
%
zei=0.05; % damping ratio of impulsive
liquid mass
gr=9.81; % acceleration due to gravity
(SI unit)
est=2.1*(10.0^11.0); % Elastic modulus of steel tank
(N/m^2)
pst=7850.0; % mass density of steel tank
material
tsm=0.02;
bb=0.2;
g=9.81;
sig=2^0.25;
dw=0.0001;
%
ep=sqrt(est/pst);
cpsi=2.0*log(sig)/(2.0*(sig-1.0)*pi);
sk=2.0/(4.0*pi*cpsi)*(sig-1.0/sig);
skk=sk/((sig-1.0)*pi);
load(‚wvcfaccsq.mat'); % Squared wavelet coefficients
are read
wv=wvcfaccsq; % Squared wavelet coefficients
are stored in a variable
clear wvcfaccsq;
vs=200.0;
ht=[3 6 10];
for ikl=1:3
for ihr=1:3
if(ihr==1)
hr=0.5;
ci=0.0719;
elseif(ihr==2)
hr=1.0;
ci=0.0875;
