Wave propagation analysis in 2D nonlinear hexagonal periodic networks based on second order gradient nonlinear constitutive models

We analyze the propagation of nonlinear waves in homogenized periodic nonlinear hexagonal networks, considering successively 1D and 2D situations. Wave analysis is performed on the basis of the construction of the effective strain energy density of periodic hexagonal lattices in the nonlinear regime. The obtained second order gradient nonlinear continuum has two propagation modes: an evanescent subsonic mode that disappears after a certain wavenumber and a supersonic mode characterized by an increase of the frequency with the wavenumber. For a weak nonlinearity, a supersonic mode occurs and the dispersion curves lie above the linear dispersion curve (v_p =v_p0). For a higher nonlinearity, the wave changes from a supersonic to an evanescent subsonic mode at s=0.7 and the dispersion curves drops below the linear case and vanish for certain values of the wavenumber. An important decrease in the frequency occurs for both subsonic and supersonic modes when the lattice becomes auxetic, and the longitudinal and shear modes become very close to each other. The influence of the lattice geometrical parameters of the lattice on the dispersion relations is analyzed.     

空间网格结构动力分析的非线性模态方法

空间网格结构因自由度数多且无简化的力学模型,非线性动力分析通常要耗费大量时间．传统的非线性模态方法用于求解多高层结构的局部非线性问题已获得良好的效果,但对系统非线性问题的应用尚缺少研究．对比分析多高层结构和空间网格结构动力性能差异,指出网格结构动力非线性分析存在的问题．以主振型理论和切线刚度分离法为基础,将非线性模态方法用于几何非线性效应显著的空间网格结构动力分析．通过对运动方程的非线性恢复力进行拆分,形成线性表达形式,然后解耦到主振型所在的广义坐标系,以达到缩减自由度数量的目的．并通过实例验证非线性模态方法的高效性与适用性．     

Hysteretic linear and nonlinear acoustic responses from pressed interfaces

An analysis of the linear and nonlinear acoustic responses from an interface between rough surfaces in elastoplastic contact is presented as a model of the ultrasonic wave interactions with imperfect interfaces and closed cracks. A micromechanical elastoplastic contact model predicts the linear and second order interfacial stiffness from the topographic and mechanical properties of the contacting surfaces during a loading–unloading cycle. The effects of those surface properties on the linear and nonlinear reflection/transmission of elastic longitudinal waves are shown. The second order harmonic amplitudes of reflected/transmitted waves decrease by more than an order of magnitude during the transition from the elastic contact mode to the elastoplastic contact mode. It is observed that under specific loading histories the interface between smooth surfaces generates higher elastoplastic hysteresis in the interfacial stiffness and the acoustic nonlinearity than interfaces between rough surfaces. The results show that when plastic flow in the contacting asperities is significant, the acoustic nonlinearity is insensitive to the asperity peak distribution. A comparison with existing experimental data for the acoustic nonlinearity in the transmitted waves is also given with a discussion on its contact mechanical implication.     

Second-harmonic generation of two-dimensional elastic wave propagation in an infinite layered structure with nonlinear spring-type interfaces

The two-dimensional elastic wave propagation in an infinite layered structure with nonlinear interlayer interfaces is analyzed theoretically to investigate the second-harmonic generation due to interfacial nonlinearity. The structure consists of identical isotropic linear elastic layers that are bonded to each other by spring-type interfaces possessing identical linear normal and shear stiffnesses but different quadratic nonlinearity parameters. Explicit analytical expressions are derived for the second-harmonic amplitudes when a single monochromatic Bloch mode propagates in the structure in arbitrary directions by applying the transfer-matrix approach and the Bloch theorem to the governing equations linearized by a perturbation method. The second-harmonic generation by a single nonlinear interface and by multiple consecutive nonlinear interfaces are shown to be profoundly influenced by the band structure of the layered structure, the fundamental Bloch wave mode, and its propagation direction. In particular, the second harmonics generated at multiple consecutive interfaces are found to grow cumulatively with the propagation distance when the phase matching occurs between the Bloch modes at the fundamental and double frequencies.     

Nonlinear elastic wave propagation in a phononic material with periodic solid–solid contact interface

Phononic materials enable enhanced dynamic properties, and offer the ability to engineer the material response. In this work we study the wave propagation in such a structure when introduced with nonlinearity. Our system is comprised of pre-compressed material with periodic solid–solid contacts, which exhibit a quadratic nonlinearity for small displacements. We suggest a new approach to modeling this system, where we discretize the unit cell in order to derive an approximate analytical solution using a perturbation method, which we are then able to easily validate numerically. With these methods, we study the band structure in the system and the second harmonic generation originating from the nonlinearity. We qualitatively analyze the second harmonic response of the system in terms of the single-crack response with linear band structure considerations. Significant band structure manipulation by changing system parameters is demonstrated, including possible in-situ tuning. The system also exhibits effective frequency doubling, i.e. the transmitted wave is primarily comprised of the second harmonic wave, for a certain range of frequencies. We demonstrate very high robustness to disorder in the system, in terms of band structure and second harmonic generation. These results have possible applications as frequency-converting devices, tunable engineered materials, and in non-destructive evaluation.     

Modeling non-collinear mixing by distributions of clapping microcracks

Nonlinear scattering by distributions of clapping cracks in a non-collinear wave mixing setting is modeled. Features of the nonlinear response discriminating distributions of clapping cracks from quadratic nonlinear damage are investigated for distributions of cracks that are parallel to each other or randomly oriented. The effective properties of these distributions are recovered extending an existing model that applies to open cracks. The equation of motion is solved using a perturbation approach, and its solutions are evaluated numerically. Their dependence on the amplitude of the incident field is found to be linear, in contrast with the quadratic dependence characterizing quadratic nonlinearity. The spectrum of the scattered field is shown to contain an infinite number of higher harmonics already at the first order of perturbation. Grating-like structures due to the opening and closing of cracks are responsible for adding diffraction peaks to the directivity functions of waves scattered by open cracks. The locations of the most prominent peaks of these functions do not satisfy the selection rules controlling nonlinear scattering by quadratic nonlinearity. Examples of these are given, together with others showing the possibility of using at least one of several discrimination modalities offered by non-collinear wave mixing.     

Rayleigh wave in a quadratic nonlinear elastic half-space (Murnaghan model)

The nonlinear equations that underlie the analysis of classical Rayleigh waves are derived for the two-dimensional case of nonlinear elastic deformation described by the Murnaghan model. In addition to the case of presence of both geometrical and physical nonlinearities, two special cases are considered where one only type of nonlinearity is taken into account. It is shown that unlike the one-dimensional problems for plane waves where only three types of nonlinear interaction should be allowed for, the two-dimensional problems should include 24 types of nonlinear interaction. In the case of geometrical nonlinearity alone, a preliminary analysis of the nonlinear equations is carried out. Second-order equations are derived. The second approximation includes the second harmonics of the wave itself and its attenuating amplitude and is nonlinearly dependent on the initial amplitude of the Rayleigh wave and linearly increasing with the distance traveled by the wave     

Quadratically Nonlinear Cylindrical Hyperelastic Waves: Primary Analysis of Evolution

The propagation and interaction of hyperelastic cylindrical waves are studied. Nonlinearity is introduced by means of the Murnaghan potential and corresponds to the quadratic nonlinearity of all basic relationships. To analyze wave propagation, an asymptotic representation of the Hankel function of the first order and first kind is used. The second-order analytical solution of the nonlinear wave equation is similar to that for a plane longitudinal wave and is the sum of the first and second harmonics, with the difference that the amplitudes of cylindrical harmonics decrease with the distance traveled by the wave. A primary computer analysis of the distortion of the initial wave profile is carried out for six classical hyperelastic materials. The transformation of the first harmonic of a cylindrical wave into the second one is demonstrated numerically. Three ways of allowing for nonlinearities are compared __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 7, pp. 73–82, July 2005.     

Internal waves in a two‐layer system using fully nonlinear internal‐wave equations

In order to understand the nonlinear effect in a two‐layer system, fully nonlinear strongly dispersive internal‐wave equations, based on a variational principle, were proposed in this study. A simple iteration method was used to solve the internal‐wave equations in order to solve the equations stably. The applicability of the proposed numerical computation scheme was confirmed to agree with linear dispersion relation theoretically obtained from variational principle. The proposed computational scheme was also shown to reproduce internal waves including higher‐order nonlinear effect from the analysis of internal solitary waves in a two‐layer system. Furthermore, for the second‐order numerical analysis, the balance of nonlinearity and dispersion was found to be similar to the balance assumed in the KdV theory and the Boussinesq‐type equations. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonlinear random coupled motions of structural elements with quadratic nonlinearities

The second-order closure method is used to analyze the nonlinear response of two-degree-of-freedom systems with quadratic nonlinearities. The excitation is assumed to be the sum of a deterministic harmonic component and a random component. The case of primary resonance of the second mode in the presence of a two-to-one internal (autoparametric) resonance is investigated. The method of multiple scales is used to obtain four first-order ordinary-differential equations that describe the modulation of the amplitudes and phases of the two modes. Applying the second-order closure method to the modulation equations, we determine the stationary mean and mean-square responses. For the case of a narrow-band random excitation, the results show that the presence of the nonlinearity causes multi-valued mean-square responses. The multi-valuedness is responsible for a jump phenomenon. Contrary to the results of the linear analysis, the nonlinear analysis reveals that the directly excited second mode takes a small amount of the input energy (saturates) and spills over the rest of the input energy into the first mode, which is indirectly excited through the autoparametric resonance.     

闭环硅微加速度计非线性补偿（英文）

为了降低闭环硅微加速度计的非线性,分析了其主要误差源并提出了相应的补偿方法。首先,分析了闭环状态下检测质量块偏离几何中心位置所造成的非线性问题,并确定了电路零位是主要误差源;其次,利用闭环反馈控制进行了非线性的优化分析;最后,提出了非线性补偿的工程调试方法。离心试验结果表明,采用该调试方法可将加速度计的非线性减小一个数量级以上。该结果验证了非线性误差分析和补偿方法的有效性,且适用于同批次加工的其它加速度计。     

Gravity-Capillary Waves in Layered Fluid

In this work, we analyze the stability of a gravity wave generated on the separation surface of two immiscible liquids inside a moving container and perturbed by a capillary wave. Such a phenomenon is experimentally observed when the amplitude and the frequency of the motion imposed to the container attain certain values. The evolution of the system is described by the variational principle. We assume that the motion of the system is decomposed into two modes: the gravity mode and the capillary mode. With suitable scaling assumptions, it is possible to show that the evolution of the gravity mode is determined by the forcing motion, while the capillary mode is excited by the nonlinear interactions between the capillary and gravity modes. Finally, an analytic dispersion relation is obtained for the pulsation of the capillary mode. This relation is a function of several quantities, all depending on the capillary wave number and the characteristics of the exciting motion.     

Do we observe Gerstner waves in wave tank experiments?

We investigate theoretically the effects of viscosity and surface films on small-amplitude Gerstner waves in deep water. The analysis is performed by using a Lagrangian formulation of fluid motion. For inviscid fluids with a free surface Gerstner waves of arbitrary amplitude are exact solutions of the nonlinear Lagrangian equations. These waves have a trochoidal surface shape. They possess vorticity, but have no mean wave momentum, i.e. induce no net drift in the fluid. By expanding the wave motion after the wave steepness as a small parameter, we demonstrate how Gerstner waves to second order in wave steepness change due to viscosity, leading to a mean drift near the surface and a backward drift beneath the surface layer, so that they conserve total (zero) mean wave momentum. In addition, if the surface is covered by a freely floating inextensible film, the mean drift at the surface (the film speed) increases dramatically. A comparison with experimental data for the drift of thin plastic sheets in wave tanks is made, showing that the presence of viscosity-modified Gerstner waves cannot be ruled out on the basis of these observations.     

An interfacial Gerstner-type trapped wave

Coastally trapped rotational interfacial waves are studied theoretically by using a Lagrangian formulation of fluid motion in a rotating ocean. The waves propagate along the interface between two immiscible inviscid incompressible fluid layers of finite depths and different densities, and are trapped at a straight wall due to the Coriolis force. For layers of finite depth, solutions are sought as series expansions after a small parameter. Comparison is made with the irrotational interfacial Kelvin wave. Both types of waves are identical to first order, having zero vorticity. The second order solution yields a relation between the vorticity and the velocity shear in the wave motion. Requiring that the mean motion in both layers is irrotational, then follows the well-known Stokes drift for interfacial Kelvin waves. On the other hand, if the mean forward drift is identically zero, we obtain the second order vorticity in the Gerstner-type wave. The solutions in both layers for the Gerstner-type interfacial wave are given analytically to second order. It is shown that small density differences and thin upper layers both act to yield a shape of the material interfacial with broader crests and sharper troughs. These effects also tend to make the particle trajectories at the interface in both layers become distorted ellipses which are flatter on the upper side. It is concluded that the effect of air excludes the possibility of observing the exact Gerstner wave in deep water.     

Wave propagation in nonlinear and hysteretic media––a numerical study

This paper considers the problem of one dimensional wave propagation in nonlinear, hysteretic media. The constitutive law in the media is prescribed by an integral relationship based on the Duhem model of hysteresis. It is found that the well known nonlinear elastic stress–strain relationship is a special case of this integral relationship. It is also shown that the stress–strain relationship from the McCall and Guyer model of hyesteretic materials can also be derived from this integral stress–strain relationship. The first part of this paper focuses on a material with a quadratic stress–strain relationship, where the initial value problem is formulated into a system of conservation laws. Analytical solutions to the Riemann problem are obtained by solving the corresponding eigenvalue problem and serve as reference for the verification and illustration of the accuracy obtained using the applied numerical scheme proposed by Kurganov and Tadmor. The second part of this research is devoted to wave propagation in hysteretic media. Several types of initial excitations are presented in order to determine special characteristics of the wave propagation due to material nonlinearity and hysteresis. The results of this paper demonstrate the accuracy and the robustness of this numerical scheme to analyze wave propagation in nonlinear materials.     

Nonlinear dynamics of a pipe conveying pulsating fluid with parametric and internal resonances

In this paper, the nonlinear planar vibration of a pipe conveying pulsatile fluid subjected to principal parametric resonance in the presence of internal resonance is investigated. The pipe is hinged to two immovable supports at both ends and conveys fluid at a velocity with a harmonically varying component over a constant mean velocity. The geometric cubic nonlinearity in the equation of motion is due to stretching effect of the pipe. The natural frequency of the second mode is approximately three times the natural frequency of the first mode for a range of mean flow velocity, resulting in a three-to-one internal resonance. The analysis is done using the method of multiple scales (MMS) by directly attacking the governing nonlinear integral-partial-differential equations and the associated boundary conditions. The resulting set of first-order ordinary differential equations governing the modulation of amplitude and phase is analyzed numerically for principal parametric resonance of first mode. Stability, bifurcation, and response behavior of the pipe are investigated. The results show new zones of instability due to the presence of internal resonance. A wide array of dynamical behavior is observed, illustrating the influence of internal resonance.     

Nonlinear dynamic instability analysis of laminated composite thin plates subjected to periodic in-plane loads

In this paper, the dynamic instability of thin laminated composite plates subjected to harmonic in-plane loading is studied based on nonlinear analysis. The equations of motion of the plate are developed using von Karman-type of plate equation including geometric nonlinearity. The nonlinear large deflection plate equations of motion are solved by using Galerkin’s technique that leads to a system of nonlinear Mathieu-Hill equations. Dynamically unstable regions, and both stable- and unstable-solution amplitudes of the steady-state vibrations are obtained by applying the Bolotin’s method. The nonlinear dynamic stability characteristics of both antisymmetric and symmetric cross-ply laminates with different lamination schemes are examined. A detailed parametric study is conducted to examine and compare the effects of the orthotropy, magnitude of both tensile and compressive longitudinal loads, aspect ratios of the plate including length-to-width and length-to-thickness ratios, and in-plane transverse wave number on the parametric resonance particularly the steady-state vibrations amplitude. The present results show good agreement with that available in the literature.     

Pseudo-Spectral simulation of 1D nonlinear propagation in heterogeneous elastic media

In this work, a 1D Pseudo-Spectral Time Domain (PSTD) algorithm has been developed for solving elastic wave equation in nonlinear heterogeneous solids using FFTs for calculation of the spatial differential operator on staggered grid. The solver uses a staggered fourth order Adams–Bashforth method, by which stress and particle velocity are updated at alternating half time steps, to integrate forward in time. To circumvent wraparound inherent to FFT-based pseudo-spectral simulation, Convolution Perfectly Matched Layer (CPML) boundary condition has been used to eliminate implementation problems linked in classical PML to the introduction in nonlinear elasticity of a time dependent bulk modulus. Different kinds of nonlinear elastic models (quadratic and cubic nonlinearity, Nazarov hysteretic nonlinearity, bi-modular nonlinearity, PM-Space nonlinearity) have been implemented. The present study will focus on the comparison of nonlinear signature (harmonics generation, shock, frequency shift and attenuation) of these different kinds of nonlinearity for rod resonance, shock wave generation. These results are expected to be useful in helping to determine the predominant nonlinear mechanism in a specific experiment.