Article ID Journal Published Year Pages File Type
695044 Automatica 2016 12 Pages PDF
Abstract

Source identification problems in a system governed by Euler–Bernoulli beam equation ρ(x)utt+(r(x)uxx)xx=F(x)H(t)ρ(x)utt+(r(x)uxx)xx=F(x)H(t), (x,t)∈(0,l)×(0,T)(x,t)∈(0,l)×(0,T), from available boundary observation (measured data), namely, from measured slope θ(t):=ux(0,t)θ(t):=ux(0,t) at x=0x=0, are considered. We propose a new approach to identifying the unknown temporal (H(t)H(t)) and spatial (F(x)F(x)) loads. This novel approach is based on weak solution theory for PDEs and quasi-solution method for inverse problems combined with the adjoint method. It allows to construct not only a mathematical theory of inverse source problems for Euler–Bernoulli beam, but also an effective numerical algorithm for reconstruction of unknown loads. Introducing the input–output operators (ΦH)(t):=ux(0,t;H)(ΦH)(t):=ux(0,t;H) and (ΨF)(t):=ux(0,t;F)(ΨF)(t):=ux(0,t;F), t∈(0,T)t∈(0,T), we show that both operators are compact. Based on this result and general regularization theory, we prove an existence of unique solutions of the regularized normal equations (Φ∗Φ+αI)Hα=Φ∗θ(Φ∗Φ+αI)Hα=Φ∗θ and (Ψ∗Ψ+αI)Fα=Ψ∗θ(Ψ∗Ψ+αI)Fα=Ψ∗θ. Then we develop the adjoint problem approach to prove Fréchet differentiability of the corresponding cost functionals and Lipschitz continuity of the Fréchet gradients. Derived explicit gradient formulas via the adjoint problem solution and known load, allow use of gradient type convergent iterative algorithms. Results of numerical simulations for benchmark problems illustrate robustness and high accuracy of the algorithm based on the proposed approach.

Related Topics
Physical Sciences and Engineering Engineering Control and Systems Engineering
Authors
, ,