Abstract
An inverse problem of identifying an unknown shear force
{g(t)}
on the inaccessible boundary
{x=l}
in a system governed by the general form Euler–Bernoulli beam equation
\rho(x)u_{tt}+\mu(x)u_{t}+(r(x)u_{xx})_{xx}=0,\quad(x,t)\in(0,l)\times(0,T)
subject to the boundary conditions
u(0,t)=u_{x}(0,t)=0,\quad u_{xx}(l,t)|_{x=l}=0,\quad-(r(x)u_{xx}(x,t))_{x}|_{x%
=l}=g(t),
is studied.
The bending moment
{\mathtt{M}(t)\coloneq-r(0)u_{xx}(0,t)}
given at the accessible boundary
{x=0}
is assumed to be a measured output.
The Neumann-to-Neumann operator
\Phi[\,{\cdot}\,]\colon\mathcal{G}\subset H^{p}(0,l)\mapsto L^{2}(0,T),\quad(%
\Phi g)(t)\coloneq-r(0)u_{xx}(0,t_{g})
corresponding to this inverse problem is shown to be compact (
{p=3}
) and Lipschitz continuous (
{p=2}
).
These properties allow us to prove the existence of a solution of the minimization problem for the Tikhonov functional
{J(g)\coloneq\lVert\Phi g-\mathtt{M}\rVert^{2}_{L^{2}(0,T)}}
.
It is proved that this functional is Fréchet differentiable.
Furthermore, an explicit formula for the Fréchet gradient of this functional is derived by making use of the unique solution to corresponding adjoint problem.
A numerical method based on Hermitian finite elements and conjugate gradient algorithm is developed for the solution of the inverse boundary value problem.
Numerical examples with random noisy measured outputs are presented to illustrate the validity and effectiveness of the proposed approach.
Identification of the Transverse Distributed Load in the Euler-Bernoulli Beam Equation from Boundary Measurement
