A Meshless Method for Retrieving Nonlinear Large External Forces on Euler-Bernoulli Beams

Chih-Wen Chang

Department of Mechanical Engineering,National United University,Miaoli,360302,Taiwan

Abstract:We retrieve unknown nonlinear large space-time dependent forces burdened with the vibrating nonlinear Euler-Bernoulli beams under varied boundary data,comprising two-end fixed,cantilevered,clamped-hinged,and simply supported conditions in this study.Even though some researchers used several schemes to overcome these forward problems of Euler-Bernoulli beams;however,an effective numerical algorithm to solve these inverse problems is still not available.We cope with the homogeneous boundary conditions,initial data,and final time datum for each type of nonlinear beam by employing a variety of boundary shape functions.The unknown nonlinear large external force can be recuperated via back-substitution of the solution into the nonlinear Euler-Bernoulli beam equation when we acquire the solution by utilizing the boundary shape function scheme and deal with a smallscale linear system to gratify an additional right-side boundary data.For the robustness and accuracy,we reveal that the current schemes are substantiated by comparing the recuperated numerical results of four instances to the exact forces,even though a large level of noise up to 50% is burdened with the overspecified conditions.The current method can be employed in the online real-time computation of unknown force functions in space-time for varied boundary supports of the vibrating nonlinear beam.

Keywords:Inverse problems;nonlinear Euler-Bernoulli beams;ill-posed problems;nonlinear space-time dependent force;boundary shape functions

1 Introduction

As we all know,retrieving external forces of Euler-Bernoulli beams plays very important roles in many engineering and scientific areas.These equations occur in the vibration of a structure,the cutting process in engineering,the sandwiches beams,the cable-stayed beams,the rotating beams,the aircraft engineering,the design of mechanical cutting tools,the nondestructive testing and so forth.

For linear external forces of linear Euler-Bernoulli equations,Han et al.[1]tested four approximate models for a transversely vibrating beam:the Euler-Bernoulli,Rayleigh,shear,and Timoshenko models.For each model,the orthogonality data were identified,and the forced response was acquired employing the approach of eigenfunction expansion.Nevertheless,they only displayed one example to discuss the second frequency spectrum.Later,expressions for the maximum cutting force,variation of the cutting force in one complete rotation,surface quality and many other characteristics of the cutting force components can be derived from the proposed cutting force model[2].The results can be shown in charts to permit machining operators to easily choose the cutting data.Abu-Hilal[3]employed Green functions to decide the dynamic response of damped Euler-Bernoulli beams and demonstrate the dynamic behavior of single and multi-span beams,single and multi-loaded beams;however,he did not show the robustness of the proposed method.Then,Andrén et al.[4]demonstrated that the clamped boring bar has non-linear dynamic properties.They also showed that non-linearities can be reduced by modifying the clamping of the boring bar.Nevertheless,they did not display the comparisons with the other literature.After that,Yoon et al.[5]used a mechanistic cutting force model to prophesy the cutting process.The experimental coefficient modelling scheme was developed for the formulation of theoretical cutting force by pondering the specific cutting force coefficient.Later,Gradiˇsek et al.[6]presented expressions for semi-empirical mechanistic identification of specific cutting and edge force coefficients for a general helical end mill from milling tests at an arbitrary radial immersion.Apart from that,Nicaise et al.[7]pondered two inverse issues of deciding point sources in vibrating beams by boundary measurements.They displayed that the boundary observation at one extremity of the area decides uniquely the sources for an arbitrarily small time of observation.The output-feedback controller for an undamped shear beam was proposed by Krstic et al.[8,9].Then,Hasanov[10]formulated new classes of inverse source problems for vibrating cantilevered elastic beams,within the range of Euler-Bernoulli beam theory.This result showed the iteration parameter via the Lipschitz constant in gradient-type methods,which were mostly employed in the numerical implementation of inverse issues.Later,Huang et al.[11]dealt with an inverse forced vibration issue,based on the conjugate gradient method (CGM),which was tested in this study to evaluate the unknown spatial and temporal-dependent external forces for the cutting tools by using the simulated beam displacement measurements.The numerical experiments were performed to test the validity of the CGM by utilizing varied types of external forces and measurement errors.After that,Liu[12]recovered an unknown space-time-dependent force in an Euler-Bernoulli beam vibration equation by an effective combination of the Lie-group adaptive approach and the differential quadrature scheme through a few iterations.Then,Kawano[13]analyzed the issue of the identification of the distribution of asynchronous vibration sources and rigidity perturbations in Euler-Bernoulli beams.Hasanov et al.[14]presented that the collocation method combined with the truncated singular value decomposition was used to estimate the degree of ill-posedness of the pondered inverse source issue.The numerical results illustrated bounds of applicability of the proposed algorithm,also its efficiency and accuracy.Later,Hasanov et al.[15]solved these two inverse issues related to asynchronous load identification that have significant engineering applications.Note that when in some inverse issues the use of boundary data was desirable or a necessity,as in the chief example of electric impedance tomography.After that,an adjoint problem approach was used for a class of inverse issues related to the identification of temporal and spatial load distributions in the Euler-Bernoulli beam equation[16];furthermore,they claimed that the constructed iterative algorithm was robust,which allowed the use of random noisy data up to 10% noise level.However,the numerical results with noisy data are not good.Maciag et al.[17]addressed an approximate scheme of solving direct and inverse problems described by Bernoulli-Euler inhomogeneous equation of vibrations of a beam;nevertheless,they did not show the robustness of their algorithm.

For the free vibration of composite beams and non-uniform beams,Liu et al.[18]proposed a new upper bound theory to estimate the first few natural frequencies.They addressed the inverse problems of composite beam equations,where they used the orthogonal system of boundary functions as bases to expand the unknown functions and derived linear algebraic equations to decide the expansion coefficients.The robustness of the current inversion approaches was shown by numerical instances.Then,when adjoint eigenfunctions were adopted as the test functions in Green’s second identity for the Euler-Bernoulli beam equation,they can develop a simple noniterative numerical method to retrieve an unknown space-dependent external force exerted on the beam[19].Therefore,they had a noniterative algorithm to retrieve the unknown force supplemented by the noisy final time displacement data.Later,Liu et al.[20]resolved the inverse source issue of a nonlinear wave equation,developing a family ofm-order homogenization functions.This method did not require the iteration to resolve nonlinear equations,which was accurate for simultaneously solving the solution Then,Liu et al.[21]addressed the higher-dimensional inverse heat source issues of nonlinear convectiondiffusion reaction equations in 2-D rectangles and 3-D cuboids,of which the final time data and the Neumann boundary conditions on one-side were over-specified.Numerical exams found that the new approach was very accurate to reveal the solution After that,Bajkowski et al.[22]contemplated a theoretical analysis and experimental examination of a sandwich beam,with a core layer made of controllable material that can change its properties over time.Numerical simulations were performed to study the possibility of shifting beam vibration frequency towards ranges distant from resonance;however,they did not display the robustness of their scheme.Apart from that,Liu et al.[23]developed a simple and effective numerical skill,which aimed to accurately and quickly deal with thin plate bending issues.Note that the proposed algorithm is quite accurate for the thin plate,the clamped plate,and the simply supported plate problems.The nonlinear primary resonance in the vibration control of cable-stayed beam with time delay feedback was studied by[24].The numerical simulation was also given to determine the optimal value for control gain and time delays that can improve the vibrations suppression efficiency;nevertheless,this approach is complex and does not discuss the noisy effect.Later,Li et al.[25]mentioned that very few results were reported on the vibration issues of some novel multilayer sandwich structures with lattice truss cores in the literature.Contributions of this literature lie in the development of the new deformation relations of multilayer sandwich beams,the construction of the dynamic models,and the systematic analysis of vibration characteristics with both the numerical and experimental schemes.Apart from that,the nonlinear free vibration and principal parametric resonance of rotating beams were investigated taking into account the lagging-axial coupling motion because of Coriolis force[26].The nonlinear equations of motion were obtained through a direct Lagrangian formulation.Later,Liu et al.[27]recovered unknown space-time-dependent forces imposed on the vibrating Euler-Bernoulli beams under varied boundary conditions.The accuracy and robustness of the present schemes were confirmed by comparing the retrieved results of several instances to the exact forces,even though considerable noise was shown in the overspecified data.

For the difficult nonlinear Euler-Bernoulli beams,Barari et al.[28]proposed the variational iteration method and parametrized perturbation method to study the non-linear vibration of beams and nonlinear responses of a clamped-clamped buckled beam.However,they did not show the robustness of their approaches and mechanical affairs.Then,Weeger et al.[29]analyzed the vibrations of nonlinear structures by means of the new method of isogeometric finite elements.Nevertheless,they did not demonstrate real applications and the noisy effect of their schemes.Kitarovic[30]resolved a nonlinear kinematics of the 2-D,non-shear-deformable and extensible Euler-Bernoulli beam imposed with the planar flexure and/or lengthening/shortening.He discussed implications of the derived formulations pertinent to the progressive collapse analysis methods based on the Smith’s approach;however,the 2-D issue was merely considered.After that,Bagheri et al.[31]coped with nonlinear responses of a clamped-clamped buckled beam to certify the vibrational behaviors of beam.They claimed that comparing with numerical results,note that the approximate solutions were in good agreement with the analytical solutions.On the other hand,about the recent nonlinear dynamical mathematical modeling and its applications,Raza and his coworkers have proposed many schemes to tackle those issues,such as the nonlinear stochastic leprosy epidemic model[32],the cervical cancer epidemic model[33],the dynamics of the pneumonia-like infections of epidemic models[34],the cancer virotherapy model[35],the pine wilt epidemic model[36],and dynamical analysis of coronavirus disease with crowding effect,and vaccination[37].

This article is arranged as follows.Section 2 displays the nonlinear problem statement and constructs the new boundary shape function,and homogeneous boundary data of the beam.In Section 3,we acquire shape functions and introduce a free parameter into the boundary shape function,which leads to a variety of boundary shape functions.Four numerical examples of the nonlinear large external forces on vibrating nonlinear Euler-Bernoulli beams are shown in Section 4.At last,we display the conclusions in Section 5.

2 Problem Statement and Establishment of Boundary Shape Function

We deliberate an inverse source issue to reveal an unknown nonlinear force functionH(x,t)being burdened with a vibrating nonlinear Euler-Bernoulli beam with varied boundary supports.It is a severely ill-posed problem in the engineering field.The nonlinear beam with simply supported boundary data is utilized as the first significant instance,for which the pair of unknown functions＜v(x,t),H(x,t)＞concurrently gratifies

whereΓ:= {(x,t)|0＜x ＜,0＜t≤tf},β:=EI/(κA)is a constant with Young’s modulusE,Ithe moment of inertia,κis the material density,andAis the cross-sectional area.

It cannot be resolved forthrightly to revealv(x,t)since the problem (1)-(3) has an unknown nonlinear force functionH(x,t).To retrieveH(x,t)and deal withv(x,t)in the entire area,we clarify the additional conditions of

First of all,we deliberate a partial boundary shape function in the time orientation:

which gratifies the conditions(2)and(4):

3 Combination of Boundary Shape Functions

4 Numerical Scheme and Experiments

4.1 Example 1 of A Simply Supported Beam

4.2 Example 2 of A Clamped-hinged Beam

Figure 2:Example 1 resolved utilizing the BSFM and illustrating the maximum errors of v and the retrieved of H

The noise isr= 0.2,and we takeβ= 1.5.We employ the BSFM withs= 3 andw= 200.Comparing the numerically retrieved solution ofH(x,t)with the exact forceH(x,t),as displayed in Figs.3a and 3c.We also illustrate the MEs ofv(x,t)and the recovery ofH(x,t)under the noisy effect in Fig.4.Excellent results are acquired with the MEs over the plane[0,1]×(0,1]being 4.95 forv(x,t)and 15039.55 forH(x,t).We also reveal the maximum absolute value ofv(x,t)andH(x,t)over the plane[0,1]×(0,1]is 1063.10 and 1.83×106,respectively,ande(H)=7.11×10-3is small.

Deliberating a large relative noise of ξ=0.2 in Eq.(40),we use the BSFM withs=3 andw=200 to retrieveH(x,t),as shown in Fig.3b.The accuracy compared with the aforementioned results is close,where the ME forv(x,t)is 3.35,the ME forH(x,t)is 15225.75,ande(H)= 8.51×10-3.For this instance,the CPU time is also less than 1.0 s.

4.3 Example 3 of A Cantilevered Beam

Figure 3:For example 2 of the unknown nonlinear force function,(a)displaying exact H,(b)numerical H with small noise r=0.2,and(c)numerical recovery of H with large noise ξ=0.2

Figure 4:Example 2 resolved utilizing the BSFM and illustrating the maximum errors of v and the retrieved of H

The noise isr=0.9,and we takeβ=5.We utilize the BSFM withs=2 andw=500.Comparing the numerically retrieved solution ofH(x,t)with the exact forceH(x,t),as exhibited in Figs.5a and 5c.We also display the MEs ofv(x,t)and the recovery ofH(x,t)under the noisy effect in Fig.6.Good results are acquired with the MEs over the plane[0,1]×(0,1]being 0.18 forv(x,t)and 103.27 forH(x,t).We also indicate the maximum absolute value ofv(x,t)andH(x,t)over the plane[0,1]×(0,1]is 25.94 and 15477.89,respectively,ande(H)=8.08×10-3is small.

Pondering a large relative noise of ξ=0.5 in Eq.(40),we use the BSFM withs=3 andw=200 to retrieveH(x,t),as shown in Fig.5b.The accuracy compared with the aforementioned results is close,where the ME forv(x,t)is 0.32,the ME forH(x,t)is 125.60,ande(H)=4.68×10-3.For this example,the CPU time is less than 1.0 s.

Figure 5:For example 3 of the unknown nonlinear force function,(a)displaying exact H,(b)numerical H with small noise r=0.9,and(c)numerical recovery of H with large noise ξ=0.5

Figure 6:Example 3 resolved utilizing the BSFM and illustrating the maximum errors of v and the retrieved of H

4.4 Example 4 of A T wo-end Fixed Beam

be an exact solution of the beam Eq.(1)under the boundary conditions in Eq.(50).The true functionH(x,t)can be obtained by introducing the above-mentionedv(x,t)into Eq.(1).

The noise isr=0.9,and we takeβ=0.2.We use the BSFM withs=4 andw=100.Comparing the numerically retrieved solution ofH(x,t)with the exact forceH(x,t),as shown in Figs.7a and 7c.We also demonstrate the MEs ofv(x,t)and the recovery ofH(x,t)under the noisy effect in Fig.8.Good results are acquired with the MEs over the plane[0,1]×(0,1]being 8.18×10-2forv(x,t)and 60.68 forH(x,t).We also exhibit the maximum absolute value ofv(x,t)andH(x,t)over the plane[0,1]×(0,1]is 7.78 and 23989.37,respectively,ande(H)=2.47×10-3is small.

Considering a large relative noise of ξ=0.01 in Eq.(40),we use the BSFM withs=8 andw=100 to retrieveH(x,t),as presented in Fig.7b.The accuracy compared with the above-mentioned results is close,in which the ME forv(x,t)is 8.53×10-2,the ME forH(x,t)is 79.82,ande(H)=2.72×10-3.For this instance,the CPU time is also less than 1.0 s.

Figure 7:(Continued)

Figure 7:For example 4 of the unknown nonlinear force function,(a)displaying exact H,(b)numerical H with small noise r=0.9,and(c)numerical recovery of H with large noise ξ=0.01

Figure 8:Example 4 resolved utilizing the BSFM and illustrating the maximum errors of v and the retrieved of H

5 Conclusions

We addressed the recovery issues of revealing unknown forces burdened with the nonlinear Euler-Bernoulli beams with four boundary-supported data,e.g.,cantilevered,simply supported,two-end fixed,and clamped-hinged beams,by utilizing the boundary shape functions method.The innovation of this scheme is the establishment of a variety of boundary shape functions with varied orders for each beam sort,such that we have influential foundations from which to extend the solution in the whole space-time realm.The proposed scheme can be utilized in the online real-time estimation of unknown force functions in space-time for varied boundary supports of the vibrating beam.On the basis of those numerical examples,we display that the proposed algorithm is applicable to the nonlinear external forces of nonlinear Euler-Bernoulli equations and pretty excellent computational efficiency,and even for adding the large random noise up to 50%.Furthermore,to the author’s best knowledge,there has no report in the literature that the numerical schemes for those four issues can offer more accurate results than the present results.The present approach can be extended to cope with the multidimensional inverse nonlinear transient PDEs and will be worked out in the future.

Funding Statement:This work was financially supported by the National United University[grant numbers 111-NUUPRJ-04].

Conflicts of Interest:The author declares that he has no conflicts of interest to report regarding the present study.