TIME-DEPENDENT AND NON-LINEAR PROBLEMS 163
Using weighted residuals as previously shown we may obtain the
transformed dynamic equivalent of the boundary integrals of Chapter
5. Note that U and P
are
the transformed displacements and tractions.
For an internal point,
u
j+
U,P*
Jt
άΓ
= f P,
1/J,άΓ
+ f ß
f
U'
Jt
dß (6.53)
Using (6.42) a similar relationship could be obtained for the
transformed stress tensor if
(σ
0
·).
If we now take the point under consideration to the boundary Γ, we
obtain the starting integral expression for the boundary element
discretisation,
cUj+j I/, PJ,άΓ=ί P
t
I/J,άΓ + f QtU),άΩ (6.54)
For a smooth boundary the constant c = .
Once (6.54) has been solved for a sufficient number of values of our
transform parameter κ a numerical inversion process must be
performed on I//*, k) to obtain the time-dependent displacements
Uj(x,
t). The procedure used by Cruse
1
for this inversion is given by
Bellman.
14
6.6 STEADY STATE ELASTODYNAMICS
We shall now examine the problem considered above from the point
of view of finding the steady state solution to the elastodynamics
problem. We saw in Section 6.3 that the application of the Fourier
transform to a differential equation resulted in a solution which was
independent of the initial conditions of the problem, the particular
integral.
If
we
start with the equation of motion of an isotropic elastic body,
equation (6.35),
(c
2
i
~cl)u
u u
+ cu
h u
+ bj = d
2
uj/dt
2
(6.55)
with boundary conditions
u
t
=
M,(X,
t) on Tj (6.56)
Pi = °ijnj = Pi(x, t) on Γ
2
(6.57)