04.07.2014 Views

download completo - SET - USP

download completo - SET - USP

download completo - SET - USP

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

A mixed BEM-FEM formulation for layered soil-superstructure interaction<br />

71<br />

Similarly, the surface forces can be expressed as:<br />

g = g ξ +<br />

(23)<br />

i 1 + g jξ2<br />

gkξ3<br />

After transforming the Cartesian coordinates into homogeneous coordinates -<br />

REDDY (1984) - and substituting equations (22) and (23) for Eq. (21), one has:<br />

T<br />

( g ξ + g ξ + g ξ )( ⋅ wξ<br />

+ w ξ + w )dA<br />

= ∫<br />

(24)<br />

e i 1 j 2 k 3 i 1 j 2 kξ3<br />

A<br />

Moreover, minimising the portion of potential energy due to external loading<br />

∫ 1 2,<br />

3<br />

A<br />

and keeping in mind that the integral f ( , ξ ξ ) dA<br />

ξ can be determined by:<br />

∫<br />

A<br />

ξ<br />

η1<br />

1<br />

ξ<br />

η2<br />

2<br />

ξ<br />

η3<br />

3<br />

dA = 2A<br />

η ! η ! η !<br />

( η + η + η + 2) !<br />

1<br />

1<br />

2<br />

2<br />

3<br />

3<br />

(25)<br />

the transverse nodal force vector is then given by:<br />

⎧ Fi<br />

⎫ ⎡2<br />

⎪ ⎪ A<br />

= ⋅<br />

⎢<br />

⎨F<br />

j ⎬ ⎢<br />

1<br />

⎪ ⎪ 12<br />

⎩ Fk<br />

⎭ ⎢⎣<br />

1<br />

1<br />

2<br />

1<br />

1⎤<br />

⎧ g<br />

⎪<br />

1<br />

⎥<br />

⎥<br />

⋅ ⎨g<br />

2⎥<br />

⎪<br />

⎦ ⎩g<br />

i<br />

k<br />

⎫<br />

⎪<br />

⎬<br />

⎪<br />

⎭<br />

j<br />

(26)<br />

By following the same procedure in the other directions, the relation between<br />

the nodal and surface forces for the case of the Free Formulation and DKT element<br />

can be written as:<br />

⎧Fi<br />

⎪<br />

⎨F<br />

j<br />

⎪<br />

⎩<br />

Fk<br />

l<br />

l<br />

l<br />

⎫ ⎧g<br />

⎪ ⎪<br />

⎬ = [ Q ] ⋅ ⎨g<br />

⎪ ⎪<br />

⎭ ⎩<br />

g<br />

l<br />

i<br />

l<br />

j<br />

l<br />

k<br />

⎫<br />

⎪<br />

⎬<br />

⎪<br />

⎭<br />

(27)<br />

where l = 1,2, 3 represent the three orthogonal directions of the element and the matrix<br />

Q given by:<br />

Cadernos de Engenharia de Estruturas, São Carlos, v.9, n. 38, p. 63-82, 2007

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

Saved successfully!

Ooh no, something went wrong!