Compound Bifurcations in the Buckling of a Delaminated Composite Strut

An analysis of the buckling and post-buckling of a delaminated composite strut is presented using a simple 4 degree of freedom nonlinear Rayleigh–Ritz formulation. Bifurcation analysis indicates that instability occurs in general at an asymmetric point of bifurcation. Depending on the depth of delamination both thin-film and overall buckling can occur in the post-buckling range, the transition being seen at a point of secondary bifurcation. For certain combinations of parameters this becomes a stellar bifurcation, associated with a double eigenvalue, where there are three possible subsequent routes for the post-buckling. The method used is fast and reliable and can be readily extended to modelling a composite with several layers.     

Thermal buckling and post-buckling of pinned–fixed Euler–Bernoulli beams on an elastic foundation

In this article, both thermal buckling and post-buckling of pinned–fixed beams resting on an elastic foundation are investigated. Based on the accurate geometrically non-linear theory for Euler–Bernoulli beams, considering both linear and non-linear elastic foundation effects, governing equations for large static deformations of the beam subjected to uniform temperature rise are derived. Due to the large deformation of the beam, the constraint forces of elastic foundation in both longitudinal and transverse directions are taken into account. The boundary value problem for the non-linear ordinary differential equations is solved effectively by using the shooting method. Characteristic curves of critical buckling temperature versus elastic foundation stiffness parameter corresponding to the first, the second, and the third buckling mode shapes are plotted. From the numerical results it can be found that the buckling load-elastic foundation stiffness curves have no intersection when the value of linear foundation stiffness parameter is less than 3000, which is different from the behaviors of symmetrically supported (pinned–pinned and fixed–fixed) beams. As we expect that the non-linear foundation stiffness parameter has no sharp influence on the critical buckling temperature and it has a slight effect on the post-buckling temperature compared with the linear one.     

Dynamic buckling of a 2-DOF imperfect system with symmetric imperfections

Non-linear dynamic buckling of a two-degree-of freedom (2-DOF) imperfect planar system with symmetric imperfections under a step load of infinite duration (autonomous system) is thoroughly discussed using energy and geometric considerations. This system under the same load applied statically exhibits (prior to limit point) an unstable symmetric bifurcation lying on the non-linear primary equilibrium path. With the aid of the total energy-balance equation of the system and the particular geometry (due to symmetric imperfections) of the plane curve corresponding to the zero level total potential energy “surface” exact dynamic buckling loads are obtained without solving the non-linear initial-value problem. The efficiency and the reliability of the technique proposed herein is demonstrated with the aid of various dynamic buckling analyses which are compared with numerical simulation using the Verner-Runge-Kutta scheme, the accuracy of which is checked via the energy-balance equation.     

Dynamics of a Flexible Cantilever Beam Carrying a Moving Mass

The motion of a flexible cantilever beam carrying a moving spring-mass system is investigated. The beam is assumed to be an Euler–Bernouli beam. The motion of the system is described by a set of two nonlinear coupled partial differential equations where the coupling terms have to be evaluated at the position of the mass. The nonlinearities arise due to the coupling between the mass and the beam. Due to the nonlinearities the system exhibits internal resonance which is investigated in this work. The equations of motion are solved numerically using the Rayleigh–Ritz method and an automatic ODE solver. An approximate solution using the perturbation method of multiple scales is also obtained.     

The Response of a Rayleigh–Liénard Oscillator to a Fundamental Resonance

A Rayleigh–Liénard oscillator excited by a fundamentalresonance is investigated by using an asymptotic perturbation method based on Fourier expansion and time rescaling. Two first-order nonlinear ordinarydifferential equations governing the modulation of the amplitude andthe phase of solutions are derived. These equations are used todetermine steady-state responses and their stability. Excitationamplitude-response and frequency-response curves are shown and checkedby numerical integration. Dulac's criterion, the Poincaré–Bendixsontheorem, and energy considerations are used in order to study the existenceand characteristics of limit cycles of the two modulation equations. Alimit cycle corresponds to a modulated motion for the Rayleigh–Liénardoscillator. For small excitation amplitude, the analytical results arein excellent agreement with the numerical solutions. In certain caseswhen the excitation amplitude is very low, an approximate analyticsolution corresponding to a modulated motion can be obtained andnumerically checked. Moreover, if the excitation amplitude is increased,an infinite-period bifurcation occurs because the modulation periodlengthens and becomes infinite, while the modulation amplitude remainsfinite and suddenly the attractor settles down into a periodic motion.     

弹性基础上的杆在轴压下的二次屈曲

本文研究弹性基础上受轴向加载的两端铰支杆当其最低两屈曲荷载很近时的后屈曲行为。首先,使用Liapunov-Schmidt约化并借助稳定性分析,揭示了杆的二次屈曲现象;基于分叉方程给出了原始后屈曲分支及二次分支的渐近展开。其次,我们使用作者建立的二次分叉的计算方法对杆的二次屈曲做了数值计算,数值结果与渐近展开符合得很好     

The plastica: The large elastic-plastic deflection of a strut

An extension of the theory of the Elastica is developed to determine the shape of a strut (or a cantilever) which undergoes large plastic deflection. The equations governing such a behaviour and known as the Plastica equations are set up and then solved by a perturbation method and by numerical integration. The post-buckling elastic-plastic deformation of a strut after initial elastic buckling is analysed, and some numerical results are given to show the variation of its load-carrying capacity with the development of the plastic region.     

Stability of Orthotropic Corrugated Cylindrical Shells

A technique for stability analysis of cylindrical shells with a corrugated midsurface is proposed. The wave crests are directed along the generatrix. The relations of shell theory include terms of higher order of smallness than those in the Mushtari–Donnell–Vlasov theory. The problem is solved using a variational equation. The Lamé parameter and curvature radius are variable and approximated by a discrete Fourier transform. The critical load and buckling mode are determined in solving an infinite system of equations for the coefficients of expansion of the resolving functions into trigonometric series. The solution accuracy increases owing to the presence of an aggregate of independent subsystems. Singularities in the buckling modes of corrugated shells corresponding to the minimum critical loads are determined. The basic, practically important conclusion is that both isotropic and orthotropic shells with sinusoidal corrugation are efficient only when their length, which depends on the waveformation parameters and the geometric and mechanical characteristics, is small     

Bifurcation Analysis of Parametrically Excited Rayleigh–Liénard Oscillators

A parametrically excited Rayleigh–Liénard oscillator is investigatedby an asymptotic perturbation method based on Fourier expansion and timerescaling. Two coupled equations for the amplitude and the phase ofsolutions are derived and the stability of steady-state periodic solutionsas well as parametric excitation-response and frequency-response curvesare determined. Comparison with the parametrically excited Liénardoscillator is performed and analytic approximate solutions are checkedusing numerical integration. Dulac's criterion, thePoincaré–Bendixson theorem, and energy considerations are used in order to study the existence and characteristics of limit cycles of the twocoupled equations. A limit cycle corresponds to a modulated motion forthe Rayleigh–Liénard oscillator. Modulated motion can be also obtainedfor very low values of the parametric excitation, and in this case, anapproximate analytic solution is easily constructed. If the parametricexcitation is increased, an infinite-period bifurcation is observed because the modulation period lengthens and becomes infinite, while themodulation amplitude remains finite and suddenly the attractor settlesdown into a periodic motion. Floquet's theory is used to evaluatethe stability of the periodic solutions, and in certain cases,symmetry-breaking bifurcations are predicted. Numerical simulationsconfirm this scenario and detect chaos and unbounded motions in theinstability regions of the periodic solutions.     

Influence of heat transfer on the dynamic response of a spherical gas/vapour bubble

The standard approach to analyse the bubble motion is the well known Rayleigh–Plesset equation. When applying the toolbox of nonlinear dynamical systems to this problem several aspects of physical modelling are usually sacrificed. Particularly in vapour bubbles the heat transfer in the liquid domain has a significant effect on the bubble motion; therefore the nonlinear energy equation coupled with the Rayleigh–Plesset equation must be solved. The main aim of this paper is to find an efficient numerical method to transform the energy equation into an ODE system, which, after coupling with the Rayleigh–Plesset equation can be analysed with the help of bifurcation theory. Due to the strong nonlinearity and violent bubble motions the computational effort can be high, thus it is essential to reduce the size of the problem as much as possible. In the first part of the paper finite difference, Galerkin and spectral collocation methods are examined and compared in terms of efficiency. In the second part free and forced oscillations are analysed with an emphasis on the influence of heat transfer. In the case of forced oscillations the unstable branches of the amplification diagrams are also computed.     

Large thermal deflections of Timoshenko beams under transversely non-uniform temperature rise

Geometrically non-linear deformation of axially extensional Timoshenko beams subjected mechanical as well thermal loadings were characterized by a system of 7 coupled and highly non-linear ordinary differential equations, which results in a complicated two-point boundary-value problem. By using shooting method this kind of problem can be numerically solved efficiently. Based on the above-mentioned mathematical formulation and numerical procedure, analysis of large thermal deflections for Timoshenko beams, subjected transversely non-uniform temperature rise and with immovably pinned–pinned as well as fixed–fixed ends, is presented. Characteristic curves showing the relationships between the beam deformation and temperature rise are illustrated. Especially, the effects of shear deformation on the bending and buckling response are quantitatively investigated. The numerical results show, as we know, that shear deformation effects become significant with the decrease of the slenderness and with the increase of the shear flexibility.     

A geometric approach for establishing dynamic buckling loads of autonomous potential N-degree-of-freedom systems

Non-linear dynamic buckling of autonomous non-dissipative N-degree-of-freedom systems whose static instability is governed either by a limit point or by an unstable symmetric bifurcation is thoroughly discussed using energy and geometric considerations. Characteristic distances associated with the geometry of the zero level total potential energy “hypersurface” in connection with total energy-balance equation lead to dynamic (global) instability criteria. These criteria allow the determination of “exact” dynamic buckling loads without solving the non-linear initial-value problem. The reliability and efficiency of the proposed geometric approach is demonstrated via several dynamic buckling analyses of 3-degree-of-freedom systems which subsequently are compared with corresponding numerical analyses based on the Verner–Runge–Kutta scheme.     

A new numerical procedure for free vibrations of pretwisted plates

Summary  A numerical procedure is proposed for the analysis of free vibrations of pretwisted thin plates. An accurate strain–displacement relationship based on the thin-shell theory combined with the finite element method using triangular plate elements with three nodes and nine degrees of freedom for each node is utilized. The vibration characteristics of pretwisted thin plates with different twist rates and aspect ratios are studied. The numerical results are compared with the previous results obtained by various types of finite elements and by the Rayleigh–Ritz method. The effect of the twist rate on the vibration characteristics is studied briefly. Received 28 February 2001; accepted for publication 18 December 2001     

Free vibration of an elastic bottom plate of a partially fluid-filled cylindrical container with an internal body

An analytical method is developed to consider the free vibration of an elastic bottom plate of a partially fluid-filled cylindrical rigid container with an internal body. The internal body is a rigid cylindrical block that is concentrically and partially submerged inside the container. The developed method captured the analytical features of the velocity potential in a non-convex, continuous, and simply connected fluid domain including the interaction between the fluid and the structure. The interaction between the fluid and the bottom plate is included. The Galerkin method is used for matching the velocity potentials appropriate to two distinct fluid regions across the common horizontal boundary (artificial horizontal boundary). Then, the Rayleigh–Ritz method is also used to calculate the natural frequencies and modes of the bottom plate of the container. The results obtained for the problem without internal body are in close agreement with both experimental and numerical results available in the articles. A finite element analysis is also used to check the validity of the present method in the presence of the internal body. Furthermore, the influences of various variables such as fluid level, internal body radius, internal body length, and the number of nodal diameters and circles on the dynamic behaviour of the coupled system are investigated.     

Characterisation of non-linear viscoelastic foods by the indentation technique

A method to determine the non-linear viscoelastic constitutive constants from indentation force–displacement data corresponding to different indentation speeds has been developed. The method consists of two parts. In the first part, the force–displacement data is expressed as two functions which represent the strain and the time-dependent responses, respectively. From these functions, the time-dependent constants and the instantaneous force–displacement response are obtained. In the second part, the strain-dependent variables are determined from the instantaneous force–displacement response through an inverse analysis based on the Levenberg–Marquardt method. The method was verified by numerical experiments using the properties of cheese as examples.     

Numerical results for the conservation of energy of Maxwell media for the Rayleigh–Stokes problem

This publication continues our studies of analytical solutions of the Rayleigh–Stokes problem for Maxwell fluids [J. Zierep, C. Fetecau, Energetic balance for the Rayleigh–Stokes problem of a Maxwell fluid, Int. J. Eng. Sci. 45 (2007) 617–627]. We start from the Fourier sine transform. The numerical result is given and discussed for the velocity u, the power of the wall shear stresses L, the dissipation Φand the boundary layer thickness δ. These new results are important for nature and technology.     

Second-order approximation of primary resonance of a disk-type piezoelectric stator for traveling wave vibration

Considering the quadratic nonlinear constitutive relations of piezoelectric materials, a traveling wave dynamic model for a lead zirconate titanate stator of a traveling wave ultrasonic motor is established using Hamilton’s principle and the Rayleigh–Ritz method. Applying the method of multiple scales, the second-order approximation of the primary resonance for traveling wave vibration of the stator is investigated. The second harmonic component is found in the primary response of the stator, which arises from the quadratic stiffness in the condition of weak excitation. In the region of the resonance, the two coupled modals are split and the lower-order peak bends to the left, hence a jump and delay exist in the response. In this way numerical results are given to verify the feasibility of the analytical approach. The results provide a theoretical foundation for further nonlinear dynamic analysis and design of the traveling wave ultrasonic motor.     

Stability Analysis of a Porous Layer Heated from Below and Subjected to Low Frequency Vibration: Frozen Time Analysis

We investigate the convection amplitude in an infinite porous layer subjected to a vibration body force that is collinear with the gravitational acceleration and heated from below. The analysis focuses on the specific case of low frequency vibration where the frozen time approximation is used. The results reveal that for moderate Vadasz numbers, increasing the magnitude of the acceleration stabilizes the convection. The results of the large Vadasz number analysis reveals that the acceleration plays a passive role in the stability of convection and the classical stability criteria for Rayleigh–Benard convection applies.     

Hydroelastic vibration of two annular plates coupled with a bounded compressible fluid

A theoretical study on a linear hydroelastic vibration of two annular plates coupled with a bounded fluid is presented. The proposed method, based on the Rayleigh–Ritz method and the finite Hankel transform, is verified through a finite element analysis by using a commercial computer code, with an excellent accuracy. It is assumed that plates with an unequal thickness and with an unequal inner radius are clamped along their edges and an inviscid compressible fluid fills the space between the annular plates and the outer rigid vessel. When the two annular plates are identical, distinct in-phase and out-of-phase modes are observed. By increasing the difference in the plate thickness, the symmetric in-phase and out-of-phase modes with respect to the middle plane of the system are gradually shifted to pseudo in-phase and out-of-phase modes, and eventually they are changed to mixed modes. It is found that the natural frequencies decrease with an increase of the fluid compressibility, and additional modes due to a fluid concentration are observed when the plates are coupled with a compressible fluid. The fluid compressibility effect on the natural frequency is dominant in the out-of-phase modes and the higher modes. Also, the effects of the fluid thickness or the distance between the plates and the inner radius of the plates on the natural frequencies of the wet modes are investigated.