LIGO-T080003-x0
16 Jul 2008, 00:00 |
2) one can integrate the strain energy density (Ws_smaxi) over the subdomain to get the total strain energy due to the Levin applied Gaussian pressure on the HR face [Postprocessing, Subdomain integration ...]
3) However, one cannot use the [Postprocessing, Boundary integration ...] in order to integrate the strain energy density at the boundary.
4) Instead export the strain energy along the barrel to be integrated in Matlab, as follows:
a) Plot the strain energy density in the subdomain
b) Draw a line for a cross-section plot which is on the barrel edge (using the auto end-point selection feature). A figure with a line plot of strain energy density versus depth (z) along the barrel will be displayed.
c) Export the plot in an ASCII file. Select the "Write script that creates postprocessing data". The Matlab script file (*.m) created has the date in the following format:
clear pd
pd.p = [...
Z1 SE1 ;...
Z2 SE2 ;...
...
Zn SEn ;...
]';
pd = rmfield(pd,'p');
d) Edit the end of the file as follows:
]';
pd.x = pd.p(1,:);
pd.y = pd.p(2,:);
pd = rmfield(pd,'p');
plot(pd.x,pd.y);
z=trapz(fliplr(pd.x),fliplr(pd.y))
diameter=0.17;
SEbarrel=z*pi*2*diameter