Non-linear dynamics of the sandwich double circular plate system

Multi-frequency vibrations of a system of two isotropic circular plates interconnected by a visco-elastic layer that has non-linear characteristics are considered. The considered physical system should be of interest to many researches from mechanical and civil engineering. The first asymptotic approximation of the solutions describing stationary and no stationary behavior, in the regions around the two coupled resonances, is the principal result of the authors. A series of the amplitude-frequency and phase-frequency curves of the two frequency like vibration regimes are presented. That curves present the evolution of the first asymptotic approximation of solutions for different non-linear harmonics obtained by changing external excitation frequencies through discrete as well as continuous values. System of the partial differential equations of the transversal oscillations of the sandwich double circular plate system with visco-non-linear elastic layer, excited by external, distributed, along plate surfaces, excitation are derived and approximately solved for various initial conditions and external excitation properties. System of differential equations of the first order with respect to the amplitudes and the corresponding number of the phases in the first asymptotic averaged approximation are derived for different corresponding multi-frequency non-linear vibration regimes. These equations are analytically and numerically considered in the light of the stationary and no stationary resonant regimes, as well as the multi-non-linear free and forced mode mutual interactions, number of the resonant jumps.     

Post-buckling of a hinged-fixed beam under uniformly distributed follower forces

Based on geometrically non-linear theory for extensible elastic beams, governing equations of statically post-buckling of a beam with one end hinged and the other fixed, subjected to a uniformly distributed, tangentially compressing follower forces are established. They consist of a boundary-value problem of ordinary differential equations with a strong non-linearity, in which seven unknown functions are contained and the arc length of the deformed axis is considered as one of the basic unknown functions. By using shooting method and in conjunction with analytical continuation, the non-linear governing equations are solved numerically and the equilibrium paths as well as the post-buckled configurations of the deformed beam are presented. A comparison between the results of conservative system and that of the non-conservative systems are given. The results show that the features of the equilibrium paths of the beams under follower loads are evidently different from that under conservative ones.     

Time-domain Green’s functions for unsaturated soils. Part II: Three-dimensional solution

The presented paper has been dedicated to complete the closed form three-dimensional fundamental solutions of the governing differential equations for an unsaturated deformable porous media with linear elastic behavior and a symmetric spherical domain in both Laplace transform and time domains. The governing differential equations consist of equilibrium, air and water transfer equations including the suction effect and dissolved air in water. The obtained Green’s functions have been derived exactly, for the first time, using the linear form of the governing differential equations and considering the effects of non-linearity of the governing equations and have been verified in both frequency and time domains.     

A robust hyper-elastic beam model under bi-axial normal-shear loadings

In this paper, a simple and robust constitutive model is proposed to simulate mechanical behaviors of hyper-elastic materials under bi-axial normal-shear loadings in the finite strain regime. The Mooney–Rivlin strain energy function is adopted to develop a two-dimensional (2D) normal-shear constitutive model within the framework of continuum mechanics. A motion field is first proposed for combined normal and shear deformations. The deformation gradient of the proposed field is calculated and then substituted into right Cauchy–Green deformation tensor. Constitutive equations are then derived for normal and shear deformations. They are two explicit coupled equations with high-level polynomial non-linearity. In order to examine capabilities of the developed hyper-elastic model, uniaxial tensile responses and non-linear stability behaviors of moderately thick straight and curved beams undergoing normal axial and transverse shear deformations are simulated and compared with experiments. Fused deposition modeling technique as a 3D printing technology is implemented to fabricate hyper-elastic beam structures from soft poly-lactic acid filaments. The printed specimens are tested under tensile/compressive in-plane and compressive out-of-plane forces. A finite element formulation along with the Newton–Raphson and Riks techniques is also developed to trace non-linear equilibrium path of beam structures in large defamation regimes. It is shown that the model is capable of predicting non-linear equilibrium characteristics of hyper-elastic straight and curved beams. It is found that the modeling of shear deformation and finite strain is essential toward an accurate prediction of the non-linear equilibrium responses of moderately thick hyper-elastic beams. Due to simplicity and accuracy, the model can serve in the future studies dealing with the analysis of hyper-elastic structures in which two normal and shear stress components are dominant.     

On non-linear dynamics of a coupled electro-mechanical system

Electro-mechanical devices are an example of coupled multi-disciplinary weakly non-linear systems. Dynamics of such systems is described in this paper by means of two mutually coupled differential equations. The first one, describing an electrical system, is of the first order and the second one, for mechanical system, is of the second order. The governing equations are coupled via linear and weakly non-linear terms. A classical perturbation method, a method of multiple scales, is used to find a steady-state response of the electro-mechanical system exposed to a harmonic close-resonance mechanical excitation. The results are verified using a numerical model created in MATLAB Simulink environment. Effect of non-linear terms on dynamical response of the coupled system is investigated; the backbone and envelope curves are analyzed. The two phenomena, which exist in the electro-mechanical system: (a)?detuning (i.e. a natural frequency variation) and (b)?damping (i.e. a decay in the amplitude of vibration), are analyzed further. An applicability range of the mathematical model is assessed.     

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.     

Non-linear global dynamics of an axially moving plate

In the present study, the geometrically non-linear dynamics of an axially moving plate is examined by constructing the bifurcation diagrams of Poincaré maps for the system in the sub and supercritical regimes. The von Kármán plate theory is employed to model the system by retaining in-plane displacements and inertia. The governing equations of motion of this gyroscopic system are obtained based on an energy method by means of the Lagrange equations which yields a set of second-order non-linear ordinary differential equations with coupled terms. A change of variables is employed to transform this set into a set of first-order non-linear ordinary differential equations. The resulting equations are solved using direct time integration, yielding time-varying generalized coordinates for the in-plane and out-of-plane motions. From these time histories, the bifurcation diagrams of Poincaré maps, phase-plane portraits, and Poincaré sections are constructed at points of interest in the parameter space for both the axial speed regimes.     

Analysis of the von Karman equations by group methods

One of the systems of equations approximating the large deflection of plates consists of two coupled non-linear fourth order partial differential equations, known as the von Karman equations. The full symmetry group for the steady equations is a finitely generated Lie group with ten parameters. For the time-dependent system the full symmetry group is an infinite parameter Lie group. Several subgroups of the full group are used to generate exact solutions of the time-independent and the time-dependent systems. These include the dilatation group (similar solutions), rotation group, screw group and others. Physical implications and applications are discussed.     

A Newton's method scheme for solving free-surface flow problems

A Newton's method scheme is described for solving the system of non-linear algebraic equations arising when finite difference approximations are applied to the Navier–Stokes equations and their associated boundary conditions. The problem studied here is the steady, buoyancy-driven motion of a deformable bubble, assumed to consist of an inviscid, incompressible gas. The linear Newton system is solved using both direct and iterative equation solvers. The numerical results are in excellent agreement with previous work, and the method achieves quadratic convergence.     

Non-linear stability analysis of imperfect thin-walled composite beams

The static non-linear behavior of thin-walled composite beams is analyzed considering the effect of initial imperfections. A simple approach is used for determining the influence of imperfection on the buckling, prebuckling and postbuckling behavior of thin-walled composite beams. The fundamental and secondary equilibrium paths of perfect and imperfect systems corresponding to a major imperfection are analyzed for the case where the perfect system has a stable symmetric bifurcation point. A geometrically non-linear theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field. An initial displacement, either in vertical or horizontal plane, is considered in presence of initial geometric imperfection. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results are presented for a simply supported beam subjected to axial or lateral load. It is shown in the examples that a major imperfection reduces the load-carrying capacity of thin-walled beams. The influence of this effect is analyzed for different fiber orientation angle of a symmetric balanced lamination. In addition, the postbuckling response obtained with the present beam model is compared with the results obtained with a shell finite element model (Abaqus).     

Accurate analytical perturbation approach for large amplitude vibration of functionally graded beams

The present work derives the accurate analytical solutions for large amplitude vibration of thin functionally graded beams. In accordance with the Euler–Bernoulli beam theory and the von Kármán type geometric non-linearity, the second-order ordinary differential equation having odd and even non-linearities can be formulated through Hamilton's principle and Galerkin's procedure. This ordinary differential equation governs the non-linear vibration of functionally graded beams with different boundary constraints. Building on the original non-linear equation, two new non-linear equations with odd non-linearity are to be constructed. Employing a generalised Senator–Bapat perturbation technique as an ingenious tool, two newly formulated non-linear equations can be solved analytically. By selecting the appropriate piecewise approximate solutions from such two new non-linear equations, the analytical approximate solutions of the original non-linear problem are established. The present solutions are directly compared to the exact solutions and the available results in the open literature. Besides, some examples are selected to confirm the accuracy and correctness of the current approach. The effects of boundary conditions and vibration amplitudes on the non-linear frequencies are also discussed.     

Solitary wavepackets from oscillatory non-linear equations with damping

In this paper, an infinite family of solutions describing solitary wave packets with a finite number of nodes is presented. These structures arise from the study of damping in the framework of non-linear ordinary differential equations with oscillatory behaviour. Usually one expects to find effects of this kind in physical systems described by a set of partial differential equations. The standard argument is that the non-linear term acts against the dispersive flux and this balance explains the appearance of solitary waves. Here we show that the non-linear oscillatory behaviour can also balance the effect of damping in special cases. The theory used to discriminate among the various possibilities is plain Painlevé analysis. Several physical applications are briefly discussed.     

Dynamic response of various von-Kármán non-linear plate models and their 3-D counterparts

Dynamic von-Kármán plate models consist of three coupled non-linear, time-dependent partial differential equations. These equations have been recently solved numerically [Kirby, R., Yosibash, Z., 2004. Solution of von-Kármán dynamic non-linear plate equations using a pseudo-spectral method. Comp. Meth. Appl. Mech. Eng. 193 (6–8) 575–599 and Yosibash, Z., Kirby, R., Gottlieb, D., 2004. Pseudo-spectral methods for the solution of the von-Kármán dynamic non-linear plate system. J. Comp. Phys. 200, 432–461] by the Legendre-collocation method in space and the implicit Newmark-β scheme in time, where highly accurate approximations were realized.Due to their complexity, these equations are often reduced by discarding some of the terms associated with time derivatives which are multiplied by the plate thickness squared (being a small parameter). Because of the non-linearities in the system of equations we herein quantitatively investigate the influence of these a-priori assumption on the solution for different plate thicknesses. As shown, the dynamic solutions of the so called “simplified von-Kármán” system do not differ much from the complete von-Kármán system for thin plates, but may have differences of few percent for plates with thicknesses to length ratio of about 1/20. Nevertheless, when investigating the modeling errors, i.e. the difference between the various von-Kármán models and the fully three-dimensional non-linear elastic plate solution, one realizes that for relatively thin plates (thickness is 1/20 of other typical dimensions), this difference is much larger. This implies that the simplified von-Kármán plate model used frequently in the literature is as good as an approximation as the complete (and more complicated) model. As a side note, it is shown that the dynamic response of any of the von-Kármán plate models, is completely different compared to the linearized plate model of Kirchhoff–Love for deflections of an order of magnitude as the plate thickness.     

Absorbing Balls for Equations Modeling Nonuniform Deformable Bodies with Heavy Rigid Attachments

We study degenerate nonlinear partial differential equations with dynamical boundary conditions describing the forced motions of nonuniform deformable bodies with heavy rigid attachments. We prove that the dynamical system generated by a discretization of these equations has an absorbing ball whose size is independent of the order of the discretization. This result implies the existence of an absorbing ball for the infinite-dimensional dynamical system corresponding to the original degenerate partial differential equation and thereby serves as a critical step for establishing the existence of global attractors for this system. Our results also address the interesting mechanical question of how nonuniformity complicates the longterm dynamics of the coupled systems we consider.     

Non-linear flexural–torsional dynamic analysis of beams of arbitrary cross section by BEM

In this paper, a boundary element method is developed for the non-linear flexural–torsional dynamic analysis of beams of arbitrary, simply or multiply connected, constant cross section, undergoing moderately large deflections and twisting rotations under general boundary conditions, taking into account the effects of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting and/or axial loading. Four boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to the angle of twist and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique leads to a system of non-linear coupled Differential–Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold–Gear Backward Differentiation Formula (BDF), a linear multistep method for differential equations coupled to algebraic equations. The geometric, inertia, torsion and warping constants are evaluated employing the Boundary Element Method. The proposed model takes into account, both the Wagner's coefficients and the shortening effect. Numerical examples are worked out to illustrate the efficiency, wherever possible the accuracy, the range of applications of the developed method as well as the influence of the non-linear effects to the response of the beam.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

Non-linear non-planar vibrations of geometrically imperfect inextensional beams. Part II—Bifurcation analysis under base excitations

The non-linear non-planar steady-state responses of a near-square cantilevered beam (a special case of inextensional beams) with general imperfection under harmonic base excitation is investigated. By applying the combination of the multiple scales method and the Galerkin procedure to two non-linear integro-differential equations derived in part I, two modulation non-linear coupled first-order differential equations are obtained for the case of a primary resonance with a one-to-one internal resonance. The modulation equations contain linear imperfection-induced terms in addition to cubic geometric and inertial terms. Variations of the steady-state response amplitude curves with different parameters are presented. Bifurcation analyses of fixed points show that the influence of geometric imperfection on the steady-state responses can be significant to a great extent although the imperfection is small. The phenomenon of frequency island generation is also observed.     

Dispersion relations and mode shapes for waves in laminated viscoelastic composites by variational methods

The propagation of oscillatory waves through periodic elastic composites has been analysed on the basis of the Floquet theory. This leads to self-adjoint differential equation systems which it was proved convenient to solve by variational methods. Many composites, such as the light-weight high-strength boron-epoxy material, consist of strong reinforcing components in a plastic matrix. The latter can exhibit viscoelastic properties which can have a significant influence on wave propagation characteristics. Replacement of the elastic constant by the viscoelastic complex modulus changes the mathematical structure so that the differential equation system is no longer self-adjoint. However, a modification of the variational principles is suggested which retains formal self-adjointness, and yields variational principles which contain additional boundary terms. These are applied to the determination of wave speeds and mode shapes for a laminated composite made of homogeneous elastic reinforcing plates in a homogeneous viscoelastic matrix for plane waves propagating normally to the reinforcing plates. These results agree well with the exact solution which can be evaluated in this simple case. The variational principles permit solutions for periodic, but otherwise arbitrary variation of material properties.     

Non-linear buckling and postbuckling behavior of thin-walled beams considering shear deformation

The static stability of thin-walled composite beams, considering shear deformation and geometrical non-linear coupling, subjected to transverse external force has been investigated in this paper. The theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field (accounting for bending and warping shear) considering moderate bending rotations and large twist. This non-linear formulation is used for analyzing the prebuckling and postbuckling behavior of simply supported, cantilever and fixed-end beams subjected to different load condition. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results show that the beam loses its stability through a stable symmetric bifurcation point and the postbuckling strength is in relation with the buckling load value. Classical predictions of lateral buckling are conservative when the prebuckling displacements are not negligible and the non-linear buckling analysis is required for reliable solutions. The analysis is supplemented by investigating the effects of the variation of load height parameter. In addition, the critical load values and postbuckling response obtained with the present beam model are compared with the results obtained with a shell finite element model (Abaqus).     

Unsteady heat and mass transfer in MHD flow over an oscillatory stretching surface with Soret and Dufour effects

In this paper, we study the unsteady coupled heat and mass transfer of two-dimensional MHD fluid over a moving oscillatory stretching surface with Soret and Dufour effects. Viscous dissipation effects are adopted in the energy equation. A uniform magnetic field is applied vertically to the flow direction. The governing equations are reduced to non-linear coupled partial differential equations and solved by means of homotopy analysis method （HAM）. The effects of some physical parameters such as magnetic parameter, Dufour number, Soret number, the Prandtl num- ber and the ratio of the oscillation frequency of the sheet to its stretching rate on the flow and heat transfer characteristics are illustrated and analyzed.