01.03.2013 Views

ACME 2011 Proceedings of the 19 UK National Conference of the ...

ACME 2011 Proceedings of the 19 UK National Conference of the ...

ACME 2011 Proceedings of the 19 UK National Conference of the ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

where<br />

and where a Kelvin unit is characterized by its shear modulus G1 and viscosity 1, and a Maxwell unit is<br />

characterized by its shear modulus G2 and viscosity 2. These creep model parameters can be obtained<br />

from laboratory creep experiments.<br />

The differential equations for <strong>the</strong> rock displacements can be expressed in <strong>the</strong> Laplace domain and <strong>the</strong><br />

solution can <strong>the</strong>n be back-transformed into <strong>the</strong> time domain.<br />

2 ANALYSIS<br />

We consider a tunnel <strong>of</strong> circular cross-section and radius r0 (Figure 1(a)). We focus here on <strong>the</strong> prediction<br />

<strong>of</strong> creep behaviour after <strong>the</strong> creation <strong>of</strong> <strong>the</strong> cylindrical cavity, thus <strong>the</strong> excavation <strong>of</strong> <strong>the</strong> tunnel is assumed<br />

to take place instantaneously. The rock surrounding <strong>the</strong> tunnel is divided into three zones (Figure 1(a)).<br />

The displacement field around <strong>the</strong> cavity is expressed as <strong>the</strong> product <strong>of</strong> two separate variables r and z<br />

(Table 1), where r u and u z are <strong>the</strong> radial and <strong>the</strong> longitudinal displacements, respectively. The terms<br />

are functions describing <strong>the</strong> variation <strong>of</strong> <strong>the</strong> radial and <strong>the</strong> longitudinal displacement in <strong>the</strong> z-direction.<br />

The terms are functions describing <strong>the</strong> attenuation <strong>of</strong> rock displacement away from <strong>the</strong> tunnel axis. As<br />

z→-∞ <strong>the</strong> plane-strain condition can be assumed, <strong>the</strong>refore d r 1 =0 and<br />

dz<br />

z 1(<br />

z)<br />

=0. At a far distance ahead<br />

<strong>of</strong> <strong>the</strong> tunnel face <strong>the</strong> displacement reduces to zero, <strong>the</strong>refore r 2 ( z)<br />

=0 and z 2 ( z)<br />

=0 as z→∞. It is<br />

assumed that 1 ( r)<br />

=1 and 2 ( r)<br />

=1 at r=r0 and 1 ( r)<br />

=0 and 2 ( r)<br />

=0 as r→∞ (to ensure that rock<br />

displacements decrease with radial distance away from <strong>the</strong> tunnel wall). Finally, it is assumed that<br />

3 ( r)<br />

=0 at r=0 and 3 ( r)<br />

=1 at r=r0 (to ensure compatibility ahead <strong>of</strong> <strong>the</strong> tunnel face).<br />

zone bounds ( r,<br />

z)<br />

I<br />

II<br />

III<br />

r 0<br />

r 0<br />

0<br />

0<br />

0<br />

z<br />

r<br />

r<br />

r<br />

z<br />

z<br />

0 r<br />

0<br />

fˆ ( s)<br />

(a) (b)<br />

Figure 1: (a) Model geometry and location <strong>of</strong> zones in which energy is dissipated; (b) Applied pressure.<br />

u r u z ( r,<br />

z)<br />

1 ( r)<br />

r1(<br />

z)<br />

2 ( r)<br />

z1(<br />

z)<br />

1 ( r)<br />

r2<br />

( z)<br />

2 ( r)<br />

z2<br />

( z)<br />

3 ( r)<br />

r2<br />

( z)<br />

z 2 ( z)<br />

Table 1: Displacements ur and uz for each <strong>of</strong> <strong>the</strong> rock zones.<br />

0<br />

f ( t)<br />

e<br />

st dt<br />

34<br />

(2)

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

Saved successfully!

Ooh no, something went wrong!