Microsoft Word - numero_58_art_19_3231.docx M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 254 Focussed on Steels and Composites for Engineering Structures Modelling of low-velocity impacts on composite beams in large displacement M. Utzeri*, M. Sasso Università Politecnica delle Marche - Department of Industrial Engineering and Mathematics, Italy m.utzeri@pm.univpm.it, m.sasso@univpm.it G. Chiappini Università degli Studi eCampus, Italy gianluca.chiappini@uniecampus.it S. Lenci Università Politecnica delle Marche - Department of Civil and Building Engineering, and Architecture, Italy s.lenci@univpm.it ABSTRACT. The paper provides an evaluation of the nonlinear dynamic response of a cantilever beam made of composite material subjected to low-velocity impacts. The structure is assumed to respond in a quasi-static manner and modelled by a continuous beam in large displacement with a lumped mass attached. First, an analytical model was developed to study the free vibrations of a beam, taking into account the nonlinearities due to large displacements and inertia. Then, the analytical findings were compared with experimental test data. The vibration of a real composite beam has been acquired through high-speed imaging technique. The displacements of the beam were extracted by digital image analysis; then, the nonlinear parameters of the analytical model were determined by the Fitting Time History technique. The results obtained by the analytical model and the experimental test are compared with numerical analysis. The validated analytical model was adapted to study a low-velocity impact; the lumped mass was associated with a rigid projectile, whose initial speed represents the impact velocity. KEYWORDS. Impact Dynamics; Composite Beam; Nonlinear Frequencies; Nonlinear Dynamics. Citation: Utzeri, M, Sasso, M., Chiappini, G., Lenci, S.., Modelling of low-velocity impacts on composite beams in large displacement, Frattura ed Integrità Strutturale, 58 (2021) 254-271. Received: 18.08.2021 Accepted: 01.09.2021 Published: 01.10.2021 Copyright: © 2021 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION n the last few years, composite materials (e.g. Carbon Fiber Reinforced Polymer, Glass Fiber Reinforced Polymer) are used extensively in widely impact engineering applications. The mains composite materials peculiarities are the high specific strength and flexibility. Therefore, a composite structure can undergo large deflection without breaks when a foreign object hits it. The theory behind the low-velocity impacts has been hugely discussed by Abrate with several works I https://youtu.be/nK2Pw5Hy7-4 M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 255 on this topic [1, 2]. The aim of his work is to extend the known model to predict the contact force history and the overall response of the structure in large displacements. The relevant models used to study the low-velocity impact are spring-mass and energy balance models [3]. In this study the spring-mass approach was used to describe the overall response of a slender cantilever beam which undergoes the impact of a projectile. Starting from the Abrate theory, the spring associated with the structure is analytically modelled as a continuous Euler-Bernoulli beam. The indentation phenomenon was neglected because the overall deflection of the beam is much larger than local. In this way, the mass of projectile can be added to the continuous cantilever beam. However, the slender beam reaches large deflection to absorb the projectile energy during the impact so only the nonlinear vibration theory allows establishing the correct dynamical behaviour. Many theoretical and experimental investigations of nonlinear vibrations of beams have appeared over the years [4, 5, 6, 7]. In general, the nonlinearity may be attributed to geometry (stiffness and inertia), and material (constitutive equation and damping) [8]. Firstly, geometric nonlinearities are caused by large displacements and slopes. Consequently, it is not possible to use the small-angle assumption, which would reduce the nonlinear curvature to the simple linear form [9, 10, 11]. For a deeper discussion on the nonlinear curvature the reader is referred to [12, 13]. Nonlinear inertia effects may be caused by the presence of concentrated or distributed masses. Actually, for transversal vibrations of an initially straight beam there are transversal and axial (linear and nonlinear) inertial effects; while the former can never be neglected, in certain cases (for example for axially immovable boundary conditions), the latter are negligible [5, 14]. According to Hamdan [15, 16, 17], the nonlinear stiffness induces nonlinear features of the hardening type, i.e., the vibration frequencies increase by increasing the amplitude of motion; on the contrary, the inertial nonlinearity, which is related to the inextensibility condition and to the lumped mass, determines a softening behaviour of the so-called “backbone curve”. Consequently, the nonlinear effects which modify the frequencies act also on the time in which the impact occurs. Therefore, the impact time is strictly connected with the amplitude reached during the impact. For instance, a hardening behaviour on the nonlinear frequency clearly reduces the impact time and the maximum amplitude reaches a lower amplitude than linear form. Accordingly, in this paper, the effects of nonlinearity in a cantilever beam were addressed further elaborating and extending the analytical approach of the works [16, 15, 18]. Analytical model takes into account the nonlinearity derived by large amplitude vibration and inertia. Starting from the energy definition, the equation of motion is determined including nonlinear terms associated to geometrical nonlinearity up to the fifth order. Hence, enforcing stationarity of the lagrangian of the system yields to the unimodal equation, which can be solved by Multiple Scale Method [19]. This analytical approach gives the approximate solution of the transversal displacement and the first nonlinear frequency where the nonlinear terms are up to fifth order. Then, the results of an experimental tests were compared with the the analytical findings to validate the nonlinear model. Free vibration experimental tests were done on a rectangular beam made of twill laminated carbon fiber reinforced polymer, whose vibration displacement amplitude was measured by high-speed imaging; indeed, the image analysis allows measuring the vibration displacement amplitude of any point of the beam. The obtained results were post-processed by the Fitting Time History (FTH) technique, which is based on the least square approximation of general damped sinusoidal function with the measured free damped vibrations of the specimen. This technique allows to determine the main natural frequencies, the related modal damping coefficients, and damping and frequency nonlinearities [20, 21]. Finally, the projectile impact against a cantilevered beam has been simulated through the commercial Finite Element software Abaqus ©. Exploiting the nonlinear Euler-Bernoulli beam coupled with implicit solver, the data about speed and displacement were exported and compared with the analytical findings. ANALYTICAL MODEL n according to Abrate [3], impacts can be modelled through spring-mass models. They provide accurate solution for low-velocity impacts. The most complete model consist in two-degree of freedom model which the projectile and the hit beam are associated with masses. The spring links each other represents the nonlinear contact stiffness. Instead, the second spring is associated with the linear and nonlinear stiffness of the structure. This kind of model can be exploited to describe the impact dynamics of projectile against a slender beam. The stiffness of the structure can be associated with the bending stiffness of Euler-Bernoulli beam. However, the hit beam undergoes a large deflection because of the slenderness. Therefore, the indentation is smaller than the overall deformation of the structure and can be neglected. The models assumes the single-degree of freedom form which the masses of the projectile and the beam are added together, as shown in Fig. (2). The low-velocity impacts models exist only when the waves phenomena inside the beam due to the hitting are neglected. In this case, flexural waves appear from the bending action of the impacted projectile. The Euler–Bernoulli theory of bending assumes that plane sections remain plane as the bending deformation. Shear deformation and rotary I M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 256 inertia effects are ignored. Considering a uniform beam of mass per unit length m, the elastic modulus E and the area moment of inertia I that undergoes a flexural wave motion with vertical displacement, the flexural waves speed is given by           1 4 F EI c m (1) where γ is the wavenumber and ω is the frequency [22, 23]. From the Eqn. (1) can be deduce the flexural waves speed is function of frequency thus they are also called dispersive waves. The phenomenon of flexural wave propagation can be neglected when the interested time is bigger than the time of flexural wave takes to travel inside the beam, FWt =  Fc L . Fig. (1) shows a simplified motion of slender cantilever beam during an impact of projectile. When the projectile hit the beam Fig. (1b), flexural waves start to travel inside the beam and have a velocity Fc . Fc is not constant as in axial waves but it is depended to the impact. More rapidly is the impulse and more is the waves speed. After the early stages, multiple reflections occur, and the bending motion is established. As the beam is bending, flexural waves keep traveling but their effects can be neglected. Consequently, the beam motion can be completely described through multi-modal superposition technique as in Fig. (1c). (a) (b) (c) Figure 1: Simplified motion of slender cantilever beam during an impact of projectile. (a) Before the impact. (b) After the impact and  FWt t (Flexural waves predominant). (c)After the impact and  FWt t (Multi-modal superposition) Multi-Modal Nonlinear Response The beam under consideration is shown in Fig. (2). It is homogeneous, i.e. the cross-section, the density, and all other properties are constant along the beam axis. The lumped mass is fixed symmetrically with respect to the beam center line to maintain the symmetry with respect to the 𝑋 reference line. Figure 2: Sketch of cantilever beam with intermediate mass M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 257 In order to describe the nonlinear response of the impacted beam in large deflection, the equations of motion are determined starting from the definition of the energies. The formulation takes into account large vibrations, the inextensibility of the beam, and the lumped mass inertial terms. As proposed by Hamdan [17], the kinetic (𝐸 ) and potential (𝐸 ) energy of the beam are defined as                    3 1 2 2 2 2 2 0 | | 2 k mL E x y d x y j (2)    1 2 0 ( , ) 2 p EI E k t d L (3) where dot means derivative with respect to time,   S L is the dimensionless distance from the left boundary,  X x L and  Y y L .   ,X t and   ,Y t are the deformed position of the generic point   , 0 in the reference straight configuration.   D L is the given dimensionless position of the lumped mass.   M mL and  3 J j mL are the inertia dimensionless parameters of the lumped mass. The term   ,k t is the curvature of the neutral axis of the beam, and is given by:  'k (4) Here prime means derivative with respect to  , and  is the rotation (equal to the beam slope, since we are considering an Euler-Bernoulli beam) [24],             2 2 2 ' '' ' ' ''1 tan     tan ' '     ' ' cos( ) ' ' y y x y x x x y (5) In this work, as in [6], the axis of the beam is considered inextensible. This implies that:  2 2' ' 1x y (6) namely,         2 4 6 2 ' ' '' 1 ' 1 ... 2 8 16 y y y x y (7) where in the last expression a Taylor approximation is considered. By means of Eqn.(7) we obtain, up to the due order,      3 5' 3 ' ' ... 6 40 y y y (8)       2 2 2 4' 1 ' ' ...y y y (9) (use is made of the boundary condition   0 0x )            2 4 2 '' ' 3 ' '' 1 ... 2 81 ' y y y k y y (10) M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 258 Note that the higher order terms are positive, meaning that the curvature increases (with respect to the linear case) for increasing amplitudes. This is the motivation for the hardening behaviour due to the geometric stiffness nonlinearity. The expression of the kinetic energy of the system then becomes:                                                            2 3 53 1 2 0 0 2 3 5 2 0 2 2 4 ' 3 ' ' ' 2 2 8 ' 3 ' ' ' 2 8 ' 1 ' ' k y ymL E y y y d d y y y y y d j y y y (11) Likewise, the potential energy becomes:           22 4 1 2 0 ' 3 ' ' ' 1 2 2 8 p y yEI E y d L (12) The vibration motion is assumed to have linear modes. The vibration of the beam with a concentrate intermediate mass is governed by the equation         ¨ '''' 1 0 m y y EI (13) Assuming a solution of Eqn.(13) in the form      , Φ i ty t e (14) Substituting the Eqn.(14) into Eqn.(13) gives               2 Φ '''' 1 Φ 0 m EI (15) and the boundary condition at  0x and  1x are           Φ 0 0    Φ 0 ' 0     Φ 1 '' 0     Φ 1 ''' 0EI EI (16) Taking the Laplace transformation of Eqn. (15), using the boundary conditions (16) and then taking the inverse Laplace transformation yields the eigenfunction of the beam with an intermediate concentrated mass                                                  2 3 2 Φ 0 '' Φ 0 ''' Φ cosh cos sinh sin 2 2 Φ sinh sin 2 k kL kL k kL kL kL H kL kL (17) where   2 4 mk EI . Hence the unknown constants 0Φ'' and 0Φ''' can be determined imposing the boundary condition in Eqn.(17). Substituting 0Φ'' and 0Φ''' into Eqn. (17), and then letting   to eliminate  Φ , gives the characteristic equation M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 259                                                1 cos cosh ( cos sinh cosh sin cos cosh 1 2 cos sinh cosh sin cos cosh 1 ) 0 kL kL kL kL kL kL kL bk bk bk bk bk bk kL kL (18) where   1b L [25]. Note that the  Φ guarantees that   Φ 1max . The Eqn.(18) can be simplified assuming   0 to obtain the characteristic equation of clamped-free beam case      cosh cosh 1 0kL kL (19) We consider that the arbitrary motion of a beam undergoing free vibration may be expressed as a superposition of its free vibration mode shapes (Φ 𝜁 ), each undergoing simple harmonic motion with frequency 𝜔 , namely we assume that: 𝑦 𝜁,𝑡 ∑ 𝑢 𝑡 Φ 𝜁 (20) where Φ 𝜁 is a given function and 𝑢 𝑡 is the unknown. The lagrangian of the system is defined as ℒ 𝐸 𝐸 . Inserting the Eqn.(20) in ℒ and enforcing stationariety, yields to the equation of motion of each mode defined as ℒ ℒ . After some computations we get [16, 26]: 𝑢 𝛼 𝛼 𝑢 𝛼 𝑢 𝑢 𝛼 𝑢 2𝛼 𝑢 𝛽 𝛼 𝑢 2𝛼 𝑢 3𝛼 𝑢 0 (21) where 𝛼 Φ 𝜁 𝑑𝜁 𝜇Φ 𝜂 𝑗Φ 𝜂 , 𝛼 Φ 𝜒 𝑑𝜒 𝑑𝜁 𝜇 Φ 𝜒 𝑑𝜒 𝑗Φ 𝜂 , 𝛼 Φ 𝜒 𝑑𝜒 Φ 𝜒 𝑑𝜒 𝑑𝜁, 𝜇 Φ 𝜒 𝑑𝜒 Φ 𝜒 𝑑𝜒 𝑗Φ 𝜂 , 𝛼 Φ 𝜁 𝑑𝜁, 𝛼 Φ 𝜁 Φ 𝜁 𝑑𝜁, 𝛼 Φ 𝜁 Φ 𝜁 𝑑𝜁 𝛽 (22) The nonlinear inertial terms are clearly visible. The accuracy of the proposed approximate solution is expected to diminish when the mass ratio  exceeds the value of 3 [16], since the spatial shapes  Φn is only considered in linear form. The Eqn.(21) is solved by the Multiple Scale Method [19], and the approximate solution is given by                           35 2 1 2 4 1 1 sin sin 3 ... 16 n n nl n n nl nu t A t A t (23) The nl is the nonlinear frequency, and can be written in the form       2 40 2 4 ...nl n n n n nA A (24) M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 260 where:                                          4 0 1 0 5 2 2 4 1 2 2 0 5 2 5 6 32 4 2 2 1 4 4 14 1 , 3 , 4   15 18 9 60 12 64   n n n n n (25) Eqn.(24) is the analytical expression of the backbone curve, approximated up to the fifth order. Considering only the first term in Eqn.(23), we remind that the physical tranvsersal position of the beam, which is needed for comparison with numerical and experimental results, is given by                       1 1 , Φ s in Φn n n nl n n n n Y t L u t L A t (26) In low-velocity impact dynamics the analytical model provides two boundary condition. The initial displacement of the system is zero. Secondly, the initial velocity of the system coincides with the impact velocity 0v . In a continuous Eulero- Bernoulli beams with a intermediate lumped mass the boundary conditions can be imposed through the Dirac function. Thereby the initial velocity along the beam is zero except in the point where the impact occurs.                 0, 0 0     , 0y y g v (27) If the impact occurs in a portion of the beam, the  g can be described through boxcar function or rectangular function as                                                    0 0Π 2 2 Imp Imp Imp g v v H H (28) where the  Imp is the adimensional portion of the beam impacted. Then, substituting the Eqn.(20) in Eqn.(27) yields             1 1 , 0 s in Φ 0,n n n n y A (29)                     0 1 1 , 0 c os Φn n n n n y v A (30) The initial displaced shape have to be equal to zero so the elementary solution is  1 0n . Thus, in this work the nA is only function of impact velocity. The magnitudes represents the strengths in which the various modes participate in the subsequent motion [27]. Multiplying both sides of Eqn.(30) by the m -th eigenfunction, integrating both sides over the length of the beam, and taking advantage of the well-known orthogonality property of the eigenfunctions yields to nA , the amplitude coefficients of the n -th mode                             1 1 00 0 1 12 2 0 0 Φ Φ Φ ( ) Φ ( ) n n n n n n n g d v d A d d (31) M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 261 EXPERIMENTAL METHOD he aim of the experimental test is to validate the analytical model proposed in the previous section. The model can be validated through a comparison with the vibration of cantilever beam with an intermediate fixed mass. In this way a more complex experimental test as gas gun test is avoided [28, 29]. The specimen is subjected to an impulse at its free end and the history of the oscillation amplitude is measured. During beam oscillations, the longitudinal section was recorded by a high speed camera (model FastCam Photron  SA4) at frame rate of 10000 Hz. By means of the Matlab Image Processing toolbox, each photo has binarized and processed to obtain the profile of beam longitudinal section, as shown in Fig.3. More detail of the experimental test can be found in Utzeri et al. work [15]. The specimen is a rectangular beam made of carbon fibre reinforced polymer pre-preg T700 TWILL. The linear density of the beam m is 0.0681 Kg m . The flexural Young’s modulus, i.e. E , are determined experimentally through the standard test method for flexural properties of polymer matrix composite materials (ASTM D7264/D7264M), and it is 44951 2 N mm . The beam cross-section is rectangular, where the base is 25.05 mm and the height is 1.85 mm. Figure 3: Image processing which shows the transverse vibration of the beam. The blue lines represents the upper and lower profiles of the beam. The circle line and the red line can be associated with the the neutral axis and the trajectory of one beam point, respectively. NONLINEARITIES IDENTIFICATION ollowing [20, 21], the results were post-processed by the Fitting Time History (FTH) technique, which is based on the least square approxi-mation of the measured free damped vibrations. FTH technique permits to reveal how the natural frequencies as well as the modal damping coefficients depend nonlinearly on the excitation amplitude in each i -th mode. Consequently, this technique is coherent with the analytical method involved in Analytical section, i.e. multi- modal approach. According to Eqn.(23), and retaining only the term proportional to A, during the free oscillations the amplitude iq of the i -th mode is given by              0, s in Φi il i n n i q t L A t (32) where     i i ti iA t B e represents the oscillation amplitude, which is assumed to vary in time because of the damping.     1il i i is the damped frequency, i is the natural frequency, i is the damping coefficient and iB and i are the starting amplitude and the phase delay, respectively, that depend on the initial conditions. The linear frequency is initially T F M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 262 computed by the Fast Fourier Transform (FFT). The nonlinear frequency is a function of the amplitude thus, according to Eqn.(24), we assume:           2 40 2 4il i i i i it A t A t (33) where  0i is the linear (natural) frequency and ij are the nonlinear correction terms of j -th order for the i -th mode. Consistently with the above explained approach, the nonlinear damping has been assumed as:            20 1 2il i i i i it A t A t (34) FTH optimization needs a cost or error function to be minimized. Therefore, the normal root means square deviation (NRSMD) between analytical function and experimental measures is used:           21 ( ) , , ,   ˆ N k kk ij ij i i q y NNRSMD B max y (35) where N is the number of acquired samples, ky is the k -th experimental observation, kq is the value of the analytical function at the same time. The NRSMD provides the optimal estimation of the vibrational properties of the system (linear,  0i and  0i , and nonlinear,  2i ,  4i ,  1i and  2i ). The minimization algorithm was created in-house with Matlab  . Finally, the nonlinear damping an frequency coefficients were progressively introduced into the Matlab optimization algorithm: first, only linear terms where considered, then a second term was added, and so on. The performance of the optimization algorithm was observed to improve increasingly. NUMERICAL MODEL he aim of the numerical modelling is to validate the analytical model proposed in the analytical section. The model can be validated through a comparison with the outcomes obtained by the impact of concentrate mass against the cantilever beam. The FEM analysis was carried out by means ABAQUS  software. The elements such as B23, B23H, B33, and B33H allow to model an Euler-Bernoulli beam in ABQUS. So, they do not allow for transverse shear deformation; plane sections initially normal to the beam’s axis remain plane (if there is no warping) and normal to the beam axis. In this study, 10000 B23 BEAM was used to model the length of the beam. The concentrate mass is considered rigid with specific initial condition. The mass has a initial speed 0v and is located transversely at the beam at D distance from the clamped end. The nonlinear transient analysis coupled with implicit solver has been used to compute the time history of the entire beam deflection. Two load step was used and a very small time step, i.e. 0.00001 s, was imposed in the FEM simulation. The first step defines the initial condition of the beam and the mass. In the second step the condition on the mass speed was deleted. So, the mass can impact against the beam and loose its velocity during the bending motion of the beam. To avoid the indentation effect and rebounds the mass was imposed to be bonded to the beam. After the simulation the data about speed and displacement were exported and compared with the analytical findings. RESULTS AND DISCUSSIONS he results obtained by the analytical model were compared with experimental data to validate the nonlinear model. After all, the boundary conditions associated with low-velocity impact were imposed to the model. The analytical findings were discussed and compared to the numerical solutions. Experimental Validation The FTH technique has been applied to the displacement signals measured in the experiment described in Experimental Method Section. A concentrate mass with a mass ratio   0.5 was fixed along to the beam at   0.95 . The vibration response is analyzed through FTH technique. T T M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 263 Unimodal linear Unimodal first order nonlinearities Unimodal second order nonlinearities Multimodal NRMSD 4.14% 2.64% 1.07% 0.76% 10 149.205 149.145 148.954 148.898 12 0 4.125 4.815 4.232 14 0 0 15.5231 12.6589 10 0.0091 0.0088 0.0063 0.0076 11 0 0.1910 0.3431 0.1793 12 0 0 -6.0519 -4.0536 1B 0.2426 0.2433 0.2451 0.2475 1 1.3219 1.3184 1.3170 1.2907 20 0 0 0 1296.5 20 0 0 0 0.0041 2B 0 0 0 -0.0012 2 0 0 0 3.0167 Table 1: FTH fitting coefficients for the signal at   0.95 of the beam length The coefficients of the FTH fittings, which synthetize the experimental results and thus are the most important results of this part, are reported in Tab. 1, together with the achieved NRSMD values. It can be noted that, as expected, the NRSMD values reduces as the modeling is enhanced from linear unimodal to multimodal nonlinear. The  2 2 1 1, , ,B B terms do not have a role in the system characterization. They compare only in the minimization algorithm, and are reported for completeness of information. Figure 4: Comparison of backbone curve of the first nonlinear frequency in analytical and experimental results. It can be noted that the backbone curve obtained by experimental method is very similar to that obtained with the analytical formulation. Moreover, the multi-modal identification increases the performance of the FTH. Indeed, inserting the second mode with its relative nonlinearities, the backbone curve of first frequency approaches that obtained analytically. The analytical model was able to describes correctly the first nonlinear frequency and therefore it was validated. Further information about the efficiency of FTH technique in experimental tests can be found in these works [15, 14]. M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 264 Analytical Findings First of all, the backbone curve of the first three modes are computed by means of Eqn. (24). The Fig. (5) shows the first nonlinear frequency and the influence of impacted distance  of a projectile with a mass ratio  of 50% and 200%. In both of cases, the nonlinear behaviour is hardening so the frequency increase with the vibration amplitude. The impacted distance  changes the nonlinear hardening trend. Especially, the hardening trend decrease up to the impact distance reaches the half of the beam. Instead, the increase of mass leads to an further hardening behaviour in all of case. Figure 5: Backbone curve of first nonlinear frequency. (a) Influence of impacted distance  of a projectile with a mass ratio  of 50%. (b)Influence of impacted distance  of a projectile with a mass ratio  of 200%. On the contrary the backbone curve of second and third modes are softening type, as shown in Fig. (6) and Fig. (7) respectively. Even in this cases were evaluated the influence of impacted distance  of a projectile with a mass ratio  of 50% and 200%. The impacted distance  changes the nonlinear trend. Figure 6: Backbone curve of second nonlinear frequency. (a) Influence of impacted distance  of a projectile with a mass ratio  of 50%. (b)Influence of impacted distance  of a projectile with a mass ratio  of 200% M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 265 Figure 7: Backbone curve of third nonlinear frequency. (a) Influence of impacted distance  of a projectile with a mass ratio  of 50%. (b)Influence of impacted distance  of a projectile with a mass ratio  of 200% Figure 8: Trend of adimensional impact duration in function of mass ratio  . The impact occurs at   0.75 The nonlinear behaviour of the frequency affects the response of the beam during the impacts. Following the spring-mass theory the projectile has the initial velocity and the beam is relaxed. Then, the projectile hits the beam and reduces its velocity up to the beam reaches the maximum deflection. The potential energy of the beam is equal to the initial kinetic energy of the projectile. After the maximum deflection point the projectile is pushed backward. The impact ends when the projectile overcomes the line of the clamped beam end. In large displacement the time of the impact can be reduced or augmented in function of the hardening or softening trend of nonlinear frequency. As the increase of the deflection of the beam this phenomenon becomes more and more relevant. The first nonlinear vibration is the most involved vibration when the impact occurs close to the free end of the beam. The amplitude due to the projectile impact with a velocity 0v can be found through the Eqn.(31). The amplitude associated with the first modes increase as the velocity impact increasing. Therefore, the impact time is affected by the nonlinear frequency which is function of the amplitude. Starting from the Eqn.(24) and the assumption of first vibration response, the impact time is defined as M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 266            2 41 10 12 1 14 1 Imp l T A A (36) Fig. (8) shows how the adimensional impact duration depends by the mass ratio  . Surely, the impact time is bigger as the increase of projectile mass. Moreover, the influence of nonlinearities is more relevant at high mass ratio with the same impact velocity. The explanation is related to the amplitude reached during the impact. The initial kinetic energy is strictly dependent to the mass ratio, so inevitably the beam has to deflect more at high mass ratio to absorb the impact at the same impact velocity. Likewise, can be evaluated the impact velocity in function of the maximum impact deflection. Starting from the Eqn.(31) and the assumption of first vibration response, the impact velocity can be computed as follows                 1 2 10 0 1 101 10 Φ ( )      , Φ d v A linear d (37)                       1 2 1 2 40 0 1 10 12 1 14 11 10 Φ ( )      Φ d v A A A nonlinear d (38) Figure 9: Impact velocity in function of the maximum deflection of the beam. (a) The projectile has a mass ratio  of 50%. (b) The projectile has a mass ratio  of 200%. Fig. (9) shows how the impact velocity depends by the maximum deflection of the beam with two different mass ratio  . In large displacement the impact velocity is not linear but follow the nonlinear frequency trend. Therefore, the hardening behaviour of the first frequency leads to an increase of impact velocity to obtain the predicted amplitude in the linear theory. As the mass ratio increase the hardening behaviour persists and the nonlinear terms have more relevance, as shown in Fig. (9)b). The changing of the impact distance  leads to an increase in slope of the impact velocity-amplitude relationship. Even in these case, the hardening behaviour persists but the nonlinear terms have less relevance. In impact dynamics the absorption energy on the impact is the key information as well as the displacements, and reactions forces. All information can be obtained by means of analytical model but the mono-modal approach is not able to describe the real time history of the system. For instance, the boundary conditions in Eqn.(30) impose the displacement is zero and the velocity is 0v only in the position where the impact occurs. Taking into account only the first vibration mode is not possible to have this trend M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 267 because the speed is function of only first mode shape. Following multi-modal superposition approach, a trend quite similar to the dirac function can be achieved. Therefore, a superposition of 10 modes is computed where only the first three modes takes into account nonlinear behaviour because of computational efforts. Fig. (10) shows the energy time history in a impact of projectile with   0.9 ,   0.5 , and 0 5 m v s . Figure 10: Energy time history. Figure 11: Deflection history of the free-end beam during the impact of projectile at   0.95 . (a) Influence of projectile mass with an initial speed of 5 m s . (b) Influence of projectile speed with a mass ratio  of 50%. In order to confirm the validation of the model, the kinetic energy of the system is split between the beam and the projectile. At  0t the kinetic energy of the projectile corresponds to the kinetic energy of free projectile with speed 0v . At the same time, the kinetic energy of the beam is very close to zero because it is relaxed. Increasing the numbers of modes into the model the beam reaches zero. In early stages, the beam undergoes the influence of the projectile velocity and starts to deflect. At the end of the impact, the beam continues to oscillate. The sum of the kinetic and potential energies represents the total energy of the system. The Fig. (10) shows the average of the total energy is constant. However, there are some fluctuations due to the numerical approximation in multi-modal approach and nonlinear equations. The effect of geometrical M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 268 nonlinearities are noticeable in the displacement history as shown in Fig. (11). The displacement in linear and nonlinear response are quite similar in the trend. However, the time of impact is clearly reduced due to the hardening behaviour in the first nonlinear frequency. The effect also involves the peak of deflection. In the nonlinear model, the maximum amplitude reaches a lower amplitude than linear form. It is reduced because of the increase of stiffness in large displacement or rather the hardening behaviour in the first nonlinear frequency. Numerical Results First of all, the natural frequency obtained by the analytical model was compared with the numerical model. Tab.(2) shows this comparison in which some combination of  and  were set. The frequencies in the tables confirm the accuracy of the analytical model. The first frequencies are the best approximated by the analytical model. Instead the second and third frequencies are slightly greater than numerical ones.   Mode n Analytical [Hz] FEM [Hz] 0 0 1° 39.333 39.331 - - 2° 246.39 246.39 - - 3° 690.20 690.20 0.75 0.5 1° 28.771 28.770 - - 2° 243.98 243.94 - - 3° 583.38 582.65 0.5 1 1° 32.328 32.326 - - 2° 159.04 158.99 - - 3° 690.03 689.30 0.75 2 1° 18.575 18.572 - - 2° 242.1731 242.01 - - 3° 531.9391 531.21 1 3 1° 11.697 11.697 - - 2° 176.49 176.43 - - 3° 563.31 562.80 Table 2: Comparison of natural frequency between analytical and numerical model The numerical analysis was carried out on the impact of a mass with   0.5 at the impact distance   0.9 with a initial velocity 0 5 m v s . Firstly, the early stages of the impact was evaluated. The flexural waves start to travel inside the beam after the impact of the projectile. The flexural waves is showed in Fig. (12) at two instant of time, 0.00006 second and 0.00009 s. In this time the flexural wave travels along to the beam and runs a distance of 40.58 mm, Δ WaveX in Fig. (12). The analytical flexural waves speed along the beam can be computed following the Eqn.(1). The impact occurs in a time about 0.00001 s. Hence, the flexural wave speed is 1362.362 m s and the theoretical Δ WaveX is 40.87 mm. Consequently, the flexural wave speed is verified. It is greater than the impact velocity and therefore the low-velocity assumption to model the impact is confirmed. After the early stages, the projectile and the beam start to motion following the multi-modal theory. In order to assess the compatibility of the analytical and numerical models, a comparison between the models results was evaluated. Fig. (13) shows the displacement and velocity histories of the impacted distance. The deflection computed with the analytical model overlaps the numerical ones, Fig. (13a). The impact time is perfectly described as well as the maximum deflection. Also, the analytical trend of velocity is consistent with the numerical one. However, the analytical signal appears smoother than the numerical solution, Fig. (13b). M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 269 Figure 12: Flexural waves propagation in early stages of impact ( Δ 0.00003 t s and scale factor 100). Figure 13: Impact of a mass with   0.5 at the impact distance   0.9 with a initial velocity 0 5 m v s ; Comparison of FEM and Analytical results.(a) Deflection history. (b) Velocity history CONCLUSION his paper investigated the nonlinear impact response in large displacements of a slender cantilever beam which undergoes the impact of a projectile. The spring-mass approach was used to describe the overall response and the indentation phenomenon was neglected because the overall deflection of the beam is much larger than local. The analysis compared the results coming from analytical approach with numerical model and experimental observations. As the nonlinear frequencies behaviour is function of the amplitude, the impact time is affected by the velocity impact. The hardening behaviour of the first frequency leads to a decrease of impact time as the velocity impact increase. Moreover, the hardening behaviour leads to an increase of impact velocity to obtain the same predicted amplitude in the linear theory. The experimental tests confirm the hardening behaviour of the first backbone curve. So, the analytical model was validated by the experimental findings. The numerical analysis also confirm the analytical outcomes. Consequently, the low-velocity assumption to model the impact is confirmed also in continuous beam hit by a projectile. REFERENCES [1] Beléndez T., Neipp C., Beléndez A., 2002. Large and small deflections of a cantilever beam. European Journal of Physics 23, 371. DOI: 10.1088/0143-0807/23/3/317. T M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 270 [2] Zheng Y., Shabana A. A., Zhang D., 2018. Curvature expressions for the large displacement analysis of planar beam motions. Journal of Computational and Nonlinear Dynamics 13. DOI: 10.1115/1.4037226. [3] Abrate, S., 1998. Low-velocity impact damage, in: Impact on Composite Structures. Cambridge University Press, pp. 135–160. DOI: 10.1017/cbo9780511574504.005. [4] Abrate, S., 2001. Modeling of impacts on composite structures. Composite Structures 51, 129–138. DOI: 10.1016/s0263-8223(00)00138-0. [5] Abrate, S., 2011. Impact Engineering of Composite Structures. Springer Vienna. DOI: 10.1007/978-3-7091-0523-8.  [6] Anderson T. J., Nayfeh A. H., Balachandran B., 1996. Experimental verification of the importance of the nonlinear curvature in the response of a cantilever beam. Journal of Vibration and Acoustics 118, 21–27. DOI: 10.1115/1.2889630.  [7] Awrejcewwicz J., Kryosko A. V., Soldatov V., 2012. Analysis of the nonlinear dynamics of the timoshenko flexible beams using wavelets. Journal of Computational and Nonlinear Dynamics 7. DOI: 10.1115/1.4004376.  [8] Babilio E., Lenci S., 2017. On the notion of curvature and its mechanical meaning in a geometrically exact plane beam theory. International Journal of Mechanical Science 128-129, 277–293. DOI: 10.1016/j.ijmecsci.2017.03.031.  [9] Crespo da Silva M. R. M., Glynn C. C., 1978. Nonlinear flexural-torsional dynamics of inextensional beams II. Forced motions. Journal of Structural Mechanics 6, 449–461. DOI: 10.1080/03601217808907349.  [10] Eisley J. G., 1964. Nonlinear deformation of elastic beams, rings and strings. Wear 7, 225–312. DOI: 10.1016/0043- 1648(64)90095-x. [11] Goel, R.P., 1973. Vibrations of a beam carrying a concentrated mass. Journal of Applied Mechanics 40, 821–822. DOI: 10.1115/1.3423102.  [12] Graff, K., 1991. Wave motion in elastic solids. Dover Publications, New York.   [13] Hamdan M. N., Dado M. H. F., 1997. Large amplitude free vibrations of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia. Journal of Sound and Vibration 206, 151–168. DOI: 10.1006/jsvi.1997.1081.  [14] Hinnant H. E., Hodges D. H., 1988. Nonlinear analysis of a cantilever beam. AIAA Journal 26, 1521–1537. DOI: 10.2514/3.10072.  [15] Leissa, A., Sonalla, M., 1991. Vibrations of cantilever beams with various initial conditions. Journal of Sound and Vibration 150, 83–99. DOI: https://doi.org/10.1016/0022-460X(91)90403-7.  [16] Lenci, S., Clementi, F., 2019. Flexural wave propagation in infinite beams on a unilateral elastic foundation. Nonlinear Dynamics 99, 721–735. DOI: 10.1007/s11071-019-04944-4.  [17] Lenci S., Clementi F., Rega G., 2016. A comprehensive analysis of hardening/softening behavior of shearable planar beams with whatever axial boundary constraint. Meccanica 51, 2589–2606. DOI: doi.org/10.1007/s11012-016-0614-9.  [18] Lenci S., Consolini L., Clementi F., 2017a. The use of the Fitting Time Histories method to detect the nonlinear behaviour of laminated glass. Journal of Vibration Testing and System Dynamics 1, 1–14. DOI: 10.5890/jvtsd.2017.03.001.  [19] Lenci S., Consolini L., Clementi F., Cocchi G., 2017b. Revealing nonlinear dynamical behaviour of laminated glass. Procedia Engineering 199, 1454–1459. DOI: 10.1016/j.proeng.2017.09.394.  [20] Navaee S., Elling R. E., 1991. Large deflections of cantilever beams. Transactions of the Canadian Society for Mechanical Engineering 15, 91–107. DOI: 10.1139/tcsme-1991-0005.  [21] Nayfeh A. H., 2000. Perturbation methods. John Wiley & Sons.  [22] Schwalbe D., Wagon S., 1997. The Duffing equation. John Wiley & Sons.  [23] Hamdan M. N.and Shabaneh N. H., 1997. On the large amplitude free vibrations of a restrained uniform beam carrying an intermediate lumped mass. Journal of Sound and Vibration 199, 711–736. DOI: 10.1006/jsvi.1996.0672.  [24] Utzeri, M., Bhagavatam, A., Mancini, E., Dinda, G., Sasso, M., Newaz, G., 2021a. Quasi- Static and Dynamic Behavior of Inconel 625 Obtained by Laser Metal Deposition: Experimental Characterization and Constitutive Modeling. Journal of Engineering Materials and Technology 143. DOI: 10.1115/1.4051087.  [25] Utzeri, M., Farotti, E., Coccia, M., Mancini, E., Sasso, M., 2021b. High strain rate compression behaviour of 3d printed carbon-PA. Journal of Materials Research DOI: 10.1557/s43578-021-00248-9.  [26] Utzeri, M., Sasso, M., Chiappini, G., Lenci, S., 2020. Nonlinear Vibrations of a Composite Beam in Large Displacements: Analytical, Numerical, and Experimental Approaches. Journal of Computational and Nonlinear Dynamics 16. DOI: 10.1115/1.4048913.  [27] Wagner H., 1964. Large-amplitude free vibrations of a beam. Journal of Applied Mechanics 32, 887–892. DOI: 10.1115/1.3627331.  M. Utzeri et alii, Frattura ed Integrità Strutturale, 58 (2021) 254-271; DOI: 10.3221/IGF-ESIS.58.19 271 [28] Warminski, J., Kloda, L., Lenci, S., 2020. Nonlinear vibrations of an extensional beam with tip mass in slewing motion. Meccanica 55, 2311–2335. DOI: 10.1007/s11012-020-01236-9.  [29] Zavodney L. D., Nayfeh A. H., 1989. The non-linear response of a slender beam carrying a lumped mass to a principal parametric excitation: theory and experiment. International Journal of Non-Linear Mechanics 24, 105–125. DOI: 10.1016/0020-7462(89)90003-6.