13.08.2013 Views

Rock Mechanics.pdf - Mining and Blasting

Rock Mechanics.pdf - Mining and Blasting

Rock Mechanics.pdf - Mining and Blasting

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

METHODS OF STRESS ANALYSIS<br />

The discretised form of equation 6.27 can be written<br />

n<br />

<br />

j=1<br />

n<br />

<br />

j=1<br />

se<br />

<br />

1<br />

qx<br />

−1<br />

se<br />

se<br />

qx(S)T x<br />

xi<br />

qx(S)T x<br />

yi<br />

y <br />

+ qy(S)T xi dS = txi<br />

y <br />

+ qy(S)T yi dS = tyi<br />

(6.29)<br />

(6.30)<br />

where n is the number of boundary elements, <strong>and</strong> each surface integral is evaluated<br />

over the range S j<br />

e of each boundary element j. Considering a particular element, one<br />

of the surface integrals can be expressed by<br />

<br />

<br />

qx(S)T x<br />

y <br />

xi + qy(S)T xi dS =<br />

N()T x<br />

xi ()dS<br />

d d + 1<br />

qy N()T<br />

−1<br />

y<br />

xi ()dS d (6.31)<br />

d<br />

The integrals of the interpolation function (N)–kernel (T ) products defined in equation<br />

6.31 can be evaluated readily by st<strong>and</strong>ard Gaussian quadrature methods. When all<br />

components of equations 6.29 <strong>and</strong> 6.30 have been calculated using the procedure<br />

defined in equation 6.31, it is found that for the m boundary nodes<br />

m<br />

j=1<br />

m<br />

j=1<br />

qxjT x∗<br />

x<br />

qxjT x∗<br />

x<br />

y∗<br />

+ qyjTx = txi<br />

y∗<br />

+ qyjTy = tyi<br />

(6.32)<br />

where T x∗<br />

x , etc., are the results of the various interpolation function–kernel integrations<br />

<strong>and</strong>, for the end nodes of each element, a summation with the appropriate integral for<br />

the adjacent element. When equations similar to 6.32 have been established for each<br />

of the m boundary nodes, they may be recast in the form<br />

[T ∗ ][q] = [t] (6.33)<br />

Equation 6.33 represents a set of 2m simultaneous equations in 2m unknowns, which<br />

are the nodal values of fictitious boundary load intensity.<br />

Once equation 6.33 has been solved for the vector [q] of nodal load intensities,<br />

all other problem unknowns can be calculated readily. For example, nodal displacements,<br />

or displacements at an internal point i in the medium, can be determined<br />

from<br />

182<br />

<br />

<br />

uxi = qx(S)U<br />

s<br />

x xi<br />

<br />

<br />

u yi = qx(S)U x yi<br />

s<br />

y <br />

+ qy(S)U xi dS<br />

y <br />

+ qy(S)U yi dS<br />

(6.34)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!