Mathematical modeling and experimental investigation of an embedded vibro-impact system

The subject of this work is the experimental investigation and the mathematical modeling of the impact force behavior in a vibro-impact system, where a hammer is mounted on a cart that imposes a prescribed displacement. By changing the hammer stiffness and the impact gap it is possible to investigate the impact force behavior under different excitation frequencies. The experimental data will be used to validate the mathematical model. The hammer behavior is studied in more detail using a nonlinear analysis, which shows the various responses of the hammer, such as dynamical jumps, bifurcations and chaos.     

Assessment of structural integrity and durability—a task of experimental mechanics

The necessity of health-monitoring and supervising structures will be justified under aspects of reliability and safety as well as with regard to economical reasons. Considering the achievements in measurement techniques combined with computer techniques, the requirements on the evolution of efficient monitoring systems will be indicated. The prerequisite will be pointed out, to conceive such systems in close co-ordination with the mathematical modelling of the structure. This is inalienable with concern to system-identification, as generally the control-parameters cannot be measured directly; they are to calculate on the basis of the mathematical model and the measurable structural response symptoms. This requires mathematical complicate solution of inverse problems. During service/operation many effects give rise for degradation of the structural resistance, reducing the safety and the life-time as well. The results of system identification enable the determination of damage indicators, which provide information on the scale of degradation in the course of time to estimate the limit of service-life and the residual life-time.     

An Exact Derivation from Three-Dimensions of the Linear Equations for the Bending of an Elastically Monoclinic Plate and Associated Three-Dimensional Solutions

The exact linear three-dimensional equations for a elastically monoclinic (13 constant) plate of constant thickness are reduced without approximation to a single 4th order differential equation for a thickness-weighted normal displacement plus two auxiliary equations for weighted thickness integrals of a stress function and the normal strain. The 4th order equation is of the same form as in classical (Kirchhoff) theory except the unknown is not the midsurface normal displacement. Assuming a solution of these plate equations, we construct so-called modified Saint-Venant solutions—“modified” because they involve non-zero body and surface loads. That is, solutions of the exact three-dimensional elasticity equations that exhibit no boundary layers and that are subject to a special set of body and surface loads that leave the analogous plate loads arbitrary.     

BEM–FEM coupling for the analysis of flexible propellers in non-uniform flows and validation with full-scale measurements

The first part of the paper presents a partitioned fluid–structure interaction (FSI) coupling for the non-uniform flow hydro-elastic analysis of highly flexible propellers in cavitating and non-cavitating conditions. The chosen fluid model is a potential flow solved with a boundary element method (BEM). The structural sub-problem has been modelled with a finite element method (FEM). In the present method, the fully partitioned framework allows one to use another flow or structural solver. An important feature of the present method is the time periodic way of solving the FSI problem. In a time periodic coupling, the coupling iterations are not performed per time step but on a periodic level, which is necessary for the present BEM–FEM coupling, but can also offer an improved convergence rate compared to a time step coupled method. Thus, it allows to solve the structural problem in the frequency domain, meaning that any transients, which slow down the convergence process, are not computed. As proposed in the method, the structural equations of motion can be solved in modal space, which allows for a model reduction by involving only a limited number of mode shapes.The second part of the paper includes a validation study on full-scale. For the full-scale validation study a purposely designed composite propeller with a diameter of 1 m has been manufactured. Also an underwater measurement set-up including a stereo camera system, remote control of the optics and illumination system has been developed. The propeller design and the underwater measurement set-up are described in the paper. During sea trials blade deflections have been measured in three different positions. A comparison between measured and calculated torque shows that the measured torque is much larger than computed. This is attributed to the differences between effective and nominal wakefields, where the latter one has been used for the calculations. To correct for the differences between measured and computed torque the calculated pressures have been amplified accordingly. In that way the deformations which have been computed with the BEM–FEM coupling for non-uniform flows became very similar to the measured results.     

Application of the method of direct separation of motions to the parametric stabilization of an elastic wire

The paper considers the application of the method of direct separation of motions to the investigation of distributed systems. An approach is proposed which allows one to apply the method directly to the initial equation of motion and to satisfy all boundary conditions, arising for both slow and fast components of motion. The methodology is demonstrated by means of a classical problem concerning the so-called Indian magic rope trick (Blekhman et al. in Selected topics in vibrational mechanics, vol. 11, pp. 139–149, [2004]; Champneys and Fraser in Proc. R. Soc. Lond. A 456:553–570, [2000]; in SIAM J. Appl. Math. 65(1):267–298, [2004]; Fraser and Champneys in Proc. R. Soc. Lond. A 458:1353–1373, [2002]; Galan et al. in J. Sound Vib. 280:359–377, [2005]), in which a wire with an unstable upper vertical position is stabilized due to vertical vibration of its bottom support point. The wire is modeled as a heavy Bernoulli–Euler beam with a vertically vibrating lower end. As a result of the treatment, an explicit formula is obtained for the vibrational correction to the critical flexural stiffness of the nonexcited system.     

Evidence for changes in the dynamics of Earth crust tilts caused by the large dam construction and reservoir filling at the Enguri dam international test area (Georgia)

This study is focused on the analysis of the dynamics of Earth crust tilting at the Enguri Dam International Test Area (EDITA). We aimed to perform a qualitative and quantitative evaluation of the dynamical properties of Earth tilt time series in the context of the development of a real-time telemetric monitoring system at the Enguri test area. The set of Earth tilt hourly time-series data, recorded during the construction of the Enguri power station and subsequent fill of reservoir (1971–1983), has been analyzed. Modern methods of linear and nonlinear time-series analysis have been used. It has been observed that Earth tilt sequences clearly reveal properties of low-dimensional dynamics. Quantitative changes in the dynamics of tilt generation connected with different stages of anthropogenic influence, such as the construction of the Enguri high dam and reservoir fill, have also been detected. It was shown that these changes are transient, since the character of Earth tilt dynamics has been restored at the end of observation.     

Soliton dynamics and interaction in the Bose–Einstein condensates with harmonic trapping potential and time-varying interatomic interaction

Investigated in this paper is the quasi-one-dimensional Gross–Pitaevskii equation, which describes the dynamics of the Bose–Einstein condensates with the harmonic trapping potential and time-varying interatomic interaction. Via the Horita method and symbolic computation, analytic bright N-soliton solution is obtained. One-, two- and three-soliton solutions are analyzed graphically. Based on the limit analysis on the one- and two-soliton solutions, the modulation on the speed of the matter-wave bright solitons is realized. Via the parameters, the interaction between the matter-wave solitons are adjustable. Furthermore, an approach to construct the interference between the matter-wave solitons has been proposed. Finally, investigation on the three-soliton solution verifies our conclusions drawn from the one and two solitons. Our conclusions might be useful in the fields of the control on the matter-wave solitons, atom lasers, and atomic accelerators.     

Nonlinear response and nonsmooth bifurcations of an unbalanced machine-tool spindle-bearing system

In this effort, a six-degree-of-freedom (DOF) model is presented for the study of a machine-tool spindle-bearing system. The dynamics of machine-tool spindle system supported by ball bearings can be described by a set of second order nonlinear differential equations with piecewise stiffness and damping due to the bearing clearance. To investigate the effect of bearing clearance, bifurcations and routes to chaos of this nonsmooth system, numerical simulation is carried out. Numerical results show when the inner race touches the bearing ball with a low speed, grazing bifurcation occurs. The solutions of this system evolve from quasi-periodic to chaotic orbit, from period doubled orbit to periodic orbit, and from periodic orbit to quasi-periodic orbit through grazing bifurcations. In addition, the tori doubling process to chaos which usually occurs in the impact system is also observed in this spindle-bearing system.     

Complete Characterization and Synthesis of the Response Function of Elastodynamic Networks

The response function of a network of springs and masses, an elastodynamic network, is the matrix valued function W(ω), depending on the frequency ω, mapping the displacements of some accessible or terminal nodes to the net forces at the terminals. We give necessary and sufficient conditions for a given function W(ω) to be the response function of an elastodynamic network, assuming there is no damping. In particular we construct an elastodynamic network that can mimic a suitable response in the frequency or time domain. Our characterization is valid for networks in three dimensions and also for planar networks, which are networks where all the elements, displacements and forces are in a plane. The network we design can fit within an arbitrarily small neighborhood of the convex hull of the terminal nodes, provided the springs and masses occupy an arbitrarily small volume. Additionally, we prove stability of the network response to small changes in the spring constants and/or addition of springs with small spring constants.     

Nonlinear modeling and dynamic analysis of the rotor-bearing system

To overcome the shortcomings of extreme time-consuming in solving the Reynolds equation, two efficient calculation methods, based on the free boundary theory and variational principles for the unsteady nonlinear Reynolds equation in the condition of Reynolds boundary, are presented in the paper. By employing the two mentioned methods, the nonlinear dynamic forces as well as their Jacobians of the journal bearing can be calculated saving time but with the same accuracy. Of these two methods, the one is called a Ritz model which manipulates the cavitation region by simply introducing a parameter to match the free boundary condition and, as a result, a very simple approximate formulae of oil-film pressure is being obtained. The other one is a one-dimensional FEM method which reduces the two-dimensional variational inequality to the one-dimensional algebraic complementary equations, and then a direct method is being used to solve these complementary equations, without the need of iterations, and the free boundary condition can be automatically satisfied. Meanwhile, a new order reduction method is contributed to reduce the degrees of freedom of a complex rotor-bearing system. Thus the nonlinear behavior analysis of the rotor-bearing system can be studied time-sparingly. The results in the paper show the high efficiency of the two methods as well as the abundant nonlinear phenomenon of the system, compared with the results obtained by the usual numerical solution of the Reynolds equation.     

A comparative experimental study on the use of three denoising methods for bearing defect detection

In the conditional maintenance of mechanical components by vibratory analysis, one distinguishes two types of analyses which are necessary for obtaining a reliable diagnosis. The first analysis lies in the detection of potential defects; there are currently various succeeded methods based on the treatment of the vibratory signals allowing the localization of a defect. One can quote among these methods the analysis of spectrum (with constant resolution (RC) or percentage of bands constant (PCB)), the analysis of envelope, the analysis cepstrale, the analysis time-frequencies or the analysis time-scales (wavelet). The second analysis is interested in the determination and the evaluation of severity of a defect detected to estimate the influence of this defect on the operation of a mechanism. The scalar indicators which make it possible to estimate the gravity of a defect are indicators known as total which are based on the statistical analysis of a temporal signal. However, the signals resulting from accelerometer sensors are the results of a mixture of sources of vibrations, sources being able to be allotted to one or more defects and are generally polluted by noise. This work presents the three principal methods of denoising and the study of their influence on the scalar parameters (kurtosis, factor peak, value rms) and this within the framework of the detection of defects of the chippings types of bearings.     

Results of an experimental and numerical study of the aerodynamic heating of the undersurface of delta wings with sharp leading edges at mach numbers M∞=6.1 and 8

The heat transfer between a supersonic flow and the undersurface of delta wings with leading-edge sweep angles x=65 and 70° is investigated in a shock tunnel at angles of attack 15°. The supersonic inviscid flow over these wings in regimes in which the bow shock is attached to the sharp leading edges is calculated numerically. The compressible boundary layer problem is solved for the calculated inviscid flow fields in the laminar, transition and turbulent flow zones. The calculations and experimental values of the heat flux on the surface of the wings are compared. The calculations are in satisfactory agreement with the experimental data in the laminar and transition zones, but diverge significantly (by up to 20%) in the turbulent zone.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 183–188, July–August, 1991.The authors wish to thank A. A. Golubinskii for assisting with the solution of the problem of supersonic inviscid gas flow over a wing.     

Bifurcation and route-to-chaos analyses for Mathieu–Duffing oscillator by the incremental harmonic balance method

Bifurcations and route to chaos of the Mathieu–Duffing oscillator are investigated by the incremental harmonic balance (IHB) procedure. A new scheme for selecting the initial value conditions is presented for predicting the higher order periodic solutions. A series of period-doubling bifurcation points and the threshold value of the control parameter at the onset of chaos can be calculated by the present procedure. A sequence of period-doubling bifurcation points of the oscillator are identified and found to obey the universal scale law approximately. The bifurcation diagram and phase portraits obtained by the IHB method are presented to confirm the period-doubling route-to-chaos qualitatively. It can also be noted that the phase portraits and bifurcation points agree well with those obtained by numerical time-integration.     

Autoparametric interaction of a liquid surface in a rectangular tank with an elastic support structure under 1:1 internal resonance

Autoparametric interaction of a liquid free surface in a rectangular tank with an elastic support structure, which is subjected to vertical excitation, is investigated. When the natural frequency of the structure is equal to the lowest natural frequency of liquid sloshing, this system is categorized as an autoparametric system with an internal resonance ratio 1:1. The structure is elastically supported so there is a higher possibility that the 1:1 internal resonance can be observed. The nonlinear theoretical analysis is conducted for a fluid assumed to be perfect in a tank with a finite liquid depth. The equations of motion for the first three sloshing modes are derived employing Galerkin’s technique and considering both the nonlinearity of the fluid motion, and the viscous damping effect. Then the theoretical frequency response curves for the harmonic oscillations of the structure and sloshing are determined using van der Pol’s method. The frequency response curves show that high amplitudes of the structure’s vibrations facilitate the liquid sloshing. Furthermore, the influence of the internal detuning parameter is investigated by showing the frequency response curves and bifurcation sets. Hopf bifurcations may occur followed by amplitude-modulated motions. The theoretical results are in quantitative agreement with the experimental data.     

A nonlinear visco-elastoplastic impact model and the coefficient of restitution

A visco-elastoplastic model for the impact between a compact body and a composite target is presented. The model is a combination of a nonlinear contact law that includes energy loss due to plastic deformation and a viscous element that accounts for energy losses due to wave propagation and/or damping. The governing nonlinear equations are solved numerically to obtain the response. A piecewise linear version of the model is also presented, which facilitates analytical solution. The model predictions are compared to those of the well-known and commonly used Hunt–Crossley model. The effects of the various impact parameters, such as impactor mass, velocity, plasticity, and damping, on the impact response and coefficient of restitution are investigated. The model appears to be suitable for a wide range of impact situations, with parameters that are well defined and easily calculated or measured. Furthermore, the resulting coefficient of restitution is shown to be a function of impact velocity and damping, as confirmed by published experimental data.     

Analytical and experimental investigations of the pulsed air–water jet

This paper describes a new way of generating pulsed air–water jet by entraining and mixing air into the cavity of a pulsed water jet nozzle. Based on the theory of hydro-acoustics and fluid dynamics, a theoretical model which describes the frequency characteristic of the pulsed air–water jet is outlined aimed at gaining a better understanding of this nozzle for generating pulses. The calculated result indicates that as the air hold-up increases, the jet oscillation frequency has an abrupt decrease firstly, and then reaches a minimum gradually at α (air hold-up)=0.5, finally it gets increased slightly. Furthermore, a vibration test was conducted to validate the present theoretical result. By this way, the jet oscillation frequency can be obtained by analyzing the vibration acceleration of the equal strength beam affected by the jet impinging. Thereby, it is found that the experimental result shows similar trend with the prediction of the present model. Also, the relationship between vibration acceleration and cavity length for the pulsed water jet follows a similar tendency in accord with the pulsed air–water jet, i.e. there exists a maximum for each curve and the maximum occurs at the ratio of L/d₁ (the ratio of cavity length and upstream nozzle diameter) =2.5 and 2.2, respectively. In addition, experimental results on specimens impinged by the pulsed water jet and pulsed air–water jet show that the erosion depth increases slightly with air addition within a certain range of cavity length. Further, this behavior is very close to the vibration test results. As for erosion volume, the air entrained into the cavity significantly affects the material removal rate.     

Statistical modeling and comparison with experimental data of tire–soil interaction for combined longitudinal and lateral slip

The interaction of a tire with a soft terrain has multiple sources of uncertainties such as the mechanical properties of the terrain, and the interfacial properties between the tire and the terrain. These uncertainties are best characterized using statistical methods such as the development of stochastic models of tire–soil interaction. The quality of the models can be assessed via statistical validation measures or metrics. Although validation of stochastic tire–soil interaction models has recently been reported with good results, it involves longitudinal slip only without considering lateral slip which can occur simultaneously with longitudinal motion. This paper presents results of the validation of a simple stochastic tire–soil interaction model for the more complicated case of combined slip. The statistical methods used for validation include the development of a Gaussian process metamodel, the calibration of model parameters using the approach of the maximum likelihood estimate in conjunction with new test data. The validation of the calibrated model, when compared with test data, is obtained using four validation metrics with good results.     

Gust response: Simulation of an aeroelastic experiment by a fluid–structure interaction method

In this paper fluid–structure interaction simulations regarding a gust generator experiment are presented, which has been conducted in 2010 in the Transonic Wind Tunnel in Göttingen (DNW-TWG), Germany. The main objective of the experiment was the investigation of the dynamic response problem of an elastic wing model concerning an encountering generic gust induced by a gust generator. Fluid–structure simulations, using a finite element structural model and a computational fluid dynamics model based on time-accurate, Reynolds-averaged Navier–Stokes equations, are compared to the experiment to validate the numerical methodology. Comparisons include steady and unsteady deflections of the elastic wing and pressure distributions. Finally, the results of simulated transfer functions of the gust generator to the elastic wing are presented in comparison to the test data.