Asymmetric bifurcations of thick-walled circular cylindrical elastic tubes under axial loading and external pressure

In this paper, we consider bifurcation from a circular cylindrical deformed configuration of a thick-walled circular cylindrical tube of incompressible isotropic elastic material subject to combined axial loading and external pressure. In particular, we examine both axisymmetric and asymmetric modes of bifurcation. The analysis is based on the three-dimensional incremental equilibrium equations, which are derived and then solved numerically for a specific material model using the Adams–Moulton method. We assess the effects of wall thickness and the ratio of length to (external) radius on the bifurcation behaviour.     

Localized bulging in an inflated cylindrical tube of arbitrary thickness – the effect of bending stiffness

We study localized bulging of a cylindrical hyperelastic tube of arbitrary thickness when it is subjected to the combined action of inflation and axial extension. It is shown that with the internal pressure P and resultant axial force F viewed as functions of the azimuthal stretch on the inner surface and the axial stretch, the bifurcation condition for the initiation of a localized bulge is that the Jacobian of the vector function $(P,F)$ should vanish. This is established using the dynamical systems theory by first computing the eigenvalues of a certain eigenvalue problem governing incremental deformations, and then deriving the bifurcation condition explicitly. The bifurcation condition is valid for all loading conditions, and in the special case of fixed resultant axial force it gives the expected result that the initiation pressure for localized bulging is precisely the maximum pressure in uniform inflation. It is shown that even if localized bulging cannot take place when the axial force is fixed, it is still possible if the axial stretch is fixed instead. The explicit bifurcation condition also provides a means to quantify precisely the effect of bending stiffness on the initiation pressure. It is shown that the (approximate) membrane theory gives good predictions for the initiation pressure, with a relative error less than 5%, for thickness/radius ratios up to 0.67. A two-term asymptotic bifurcation condition for localized bulging that incorporates the effect of bending stiffness is proposed, and is shown to be capable of giving extremely accurate predictions for the initiation pressure for thickness/radius ratios up to as large as 1.2.     

Flow pattern and stability of pseudoplastic axial Taylor–Couette flow

The effect of an axial flow on the stability of the Taylor–Couette flow is explored for pseudoplastic fluids. The fluid is assumed to follow the Carreau–Bird model and mixed boundary conditions are imposed while the axial flow can be independent of rotational flow. The four-dimensional low-order dynamical system, resulted from Galerkin projection of the conservation of mass and momentum equations, includes additional non-linear terms in the velocity components originated from the shear-dependent viscosity. In absence of axial flow the base flow loses its radial flow stability to the vortex structure at a lower critical Taylor number, as the pseudoplasticity effects increases. The emergence of the vortices corresponds to the onset of a supercritical bifurcation which is also seen in the flow of a linear fluid. However, unlike the Newtonian case, pseudoplastic Taylor vortices lose their stability as the Taylor number reaches a second critical number corresponding to the onset of a Hopf bifurcation. Existence of an axial flow, induced by a pressure gradient appears to further advance each critical point on the bifurcation diagram. Complete flow field together with viscosity maps are given for stability regions in the bifurcation diagram.     

Instability and failure of internally pressurized ductile metal cylinders

Forlong, ductile, thick-walled tubes under internal pressure instabilities and final failure modes are studied experimentally and theoretically. The test specimens are closed-end cylinders made of an aluminum alloy and of pure copper and the experiments have been carried out for a number of different initial external radius to internal radius ratios. The experiments show necking on one side of the tubes at a stage somewhat beyond the maximum internal pressure. All tubes, except for one aluminum alloy tube, failed by shear fracture under decreasing pressure. The aluminum alloy tubes exhibited localized shear deformations in the neck region prior to fracture and also occasionally surface wave instabilities. The numerical investigation is based on an elastic-plastic material model for a solid that develops a vertex on the yield surface, using representations of the uniaxial stress-strain curves found experimentally. In contrast to the simplest flow theory of plasticity this material model predicts shear band instabilities at a realistic level of strain. A rather sharp vertex is used in the material model for the aluminum alloy, while a more blunt vertex is used to characterize copper. The theoretically predicted bifurcation into a necking mode, the cross-sectional shape of the neck, and finally the initiation and growth of shear bands from the highly strained internal surface in the neck region are in good agreement with the experimental observations.     

Axisymmetric wrinkling of snug-fit lined pipe

An analytical bifurcation solution is presented for axisymmetric wrinkling on a lined pipe under axial compression without internal pressure. The internal liner consists of corrosion-resistant alloy (CRA), it is not metallurgically bonded to the carbon steel backing pipe, and it is assumed to be in a snug fit condition: i.e. there is no gap between the liner and the backing pipe, but also no prestress that would lead to a positive contact or gripping pressure between the liner and the backing. The backing is assumed to be much thicker than the liner, so that wrinkling-related deformations of the backing pipe can be neglected.The solution indicates that the incipient wrinkling strain for the snug-fit pipe without any imperfections is the same as the incipient wrinkling strain for a single pipe with (5/3) times the wall thickness of the liner, and the same midsurface diameter, as determined by the solution of Batterman (1965) for the case of small strains, or Peek (2000a) for the case of finite strains. For the case when the liner-pipe friction is included the factor (5/3) increases slightly.A positive contact pressure due to prestress or internal pressure raises the wrinkling strain, whereas imperfections (e.g. at seam or girth welds) reduces it. The snug-fit solution accounts for neither, but nevertheless provides a useful reference wrinkling strain, and can be used to validate numerical solutions, and it gives a bifurcation modeshape and wrinkle length that can be used in numerical models to investigate post-bifurcation behaviour.     

An analytic study of the pressure distribution in operation of multistratum petroleum deposits

A study is made of the pressure distributions in the two strata of a circular deposit (radius R₊) extracted by different well systems. It is assumed that the two strata are separated by a stratum with much inferior collector properties.We are indebted to V. N. Shchelkachev for proposing this topic and for valuable comments on the work.     

The Effect of Elastic Surface Coating on the Finite Deformation and Bifurcation of a Pressurized Circular Annulus

A plane-strain theory of an elastic solid coated with a thin elastic film on part or all of its boundary was developed recently by Steigmann and Ogden (1997a). In this paper the theory is applied to the (plane-strain) problem of a thick-walled circular cylindrical tube which is subject to both internal and external pressure and which has an elastic coating on one or both of its circular cylindrical boundaries. The effect of the coating on the symmetrical response of the annular cross-section of the tube is determined first. It is noted, in particular, that while the pressure may exhibit a maximum followed by a minimum during inflation for an uncoated tube it may be a monotonic increasing function of the radius for a coated tube with coating elastic modulus sufficiently large. Next, the possibility of bifurcation from a symmetrical configuration is examined and again the influence of the coating is analysed. The effect of a coating on the outer boundary is compared with that on the inner boundary. Specifically, during compression, coating on the outer boundary delays bifurcation compared with the uncoated case. On the other hand, when the coating is on the inner boundary, bifurcation is either delayed or advanced relative to the uncoated situation depending on the values of the bending stiffness and tube thickness parameters. Generally, bifurcation is delayed by an increase in the magnitude of the bending stiffness of the coating at fixed values of the other parameters. This revised version was published online in August 2006 with corrections to the Cover Date.     

Effect of time-dependent speed on frictional heat generation and wear in transient axisymmetrical contact of sliding

Summary A transient contact problem with frictional heating and wear for two nonuniform sliding half-spaces is considered. One of the two half-spaces is assumed to be slightly curved to give a Hertzian initial pressure distribution: the other is a rigid nonconductor. Under the assumption that the contact pressure distribution could be described by Hertz formulas during all the process of interaction, the problem is formulated in terms of one integral equation of Volterra type with unknown radius of contact area. A numerical solution of this equation is obtained using a piecewise-constant presentation of an unknown function. The influence of operating parameters on the contact temperature and the radius of the contact area is studied. Accepted for publication 3 November 1996     

Instability analysis in pressurized transversely isotropic Saint-Venant–Kirchhoff and neo-Hookean cylindrical thick shells

Under the effect of an inner pressure, a thick hyperelastic shell cylinder is susceptible to developing mechanical instabilities, leading to a bifurcated shape that no longer has the initial cylindrical symmetry. The considered problem has a strong resonance in a biomedical context, considering pathologies encountered for arteries, such as aneurysm. In this contribution, the mechanical behavior of thick elastic shells has been analyzed, considering a thick-walled cylindrical hyperelastic model material obeying a transversely isotropic behavior, first in a large-displacement situation and then in a large-deformation case. The response of the material is assumed to be instantaneous, so that time-dependent effects shall not be considered in this paper. The case of a Saint-Venant–Kirchhoff material is considered with special emphasis to exemplify a large-displacement small-strain situation; the neo-Hookean behavior is next considered to enlarge the constitutive law toward consideration of large strains. The stability conditions of the shell are studied and bifurcation conditions formulated in terms of the applied pressure and of the geometrical and mechanical parameters that characterize the shell. Analytical solutions of some bifurcation points are evidenced and calculated when the direction of the fibers coincide with the cylinder axis.     

基于气压鼓泡法研究石墨烯力学性能

通常假设二维材料为连续介质薄膜,然后采用连续介质薄膜的研究方法进行二维材料力学性能研究,其中气压鼓泡法是一种主要测试方法．但实验观测发现,悬空石墨烯并非处于气压鼓泡测试分析模型中假设的固支边界条件,而是处于一种粘附边界条件：靠近孔壁边界处,有小部分材料通过范德华吸引粘附在基底柱形孔的侧壁上,而且粘附部分可以在极小载荷作用下剥离．这导致石墨烯悬空部分的实际半径大于基底孔半径,即鼓泡实验中的石墨烯是一种松弛薄膜,而非通常认为的预拉伸薄膜．通过有限元数值模拟研究发现,可基于含有名义松弛应变的鼓泡分析模型获得处于粘附边界条件下的石墨烯弹性模量．     

Pressure dependence of the instability of multiwalled carbon nanotubes conveying fluids

Based on an elastic beam model, the instability of multiwalled carbon nanotubes (MWCNTs) induced by the moving fluid inside is investigated. At critical flow velocities, the MWCNTs become unstable and undergo pitchfork bifurcation and subsequently Hopf bifurcation. These critical velocities are found to increase very quickly with respect to decreasing inner radius and are inversely proportional to the length-to-outer-radius ratio. The effect of the van der Waals (vdW) interaction between tubes is investigated and it is found that the vdW interaction can enhance the stability of MWCNTs in general, but the vdW interaction reduces the stability capacity of MWCNTs with very small inner radius.     

Pulsatile blood flow in a resistance vessel

The quasi-one-dimensional Newtonian fluid flow in an active vessel in which the inlet pressure oscillates about a mean value is considered. The vascular wall is assumed to display the Bayliss effect and a reaction to changes in blood flow. The first is characterized by the dependence of the contractility of the vascular smooth-muscle cells on the intravascular pressure and the second by the sensitivity of the endothelial cells to the wall shear stress. The dependence of the period-average radius and blood flow rate on the mean value and pulsation amplitude and frequency of the inlet pressure is investigated numerically. The characteristics are compared for the passive vessel, the vessel with the Bayliss effect, and the vessel with both reactions.     

Chaotic behavior of gas bubble in non-Newtonian fluid: a numerical study

In the present paper, the nonlinear behavior of bubble growth under the excitation of an acoustic pressure pulse in non-Newtonian fluid domain has been investigated. Due to the importance of the bubble in the medical applications such as drug, protein or gene delivery, blood is assumed to be the reference fluid. Effects of viscoelasticity term, Deborah number, amplitude and frequency of the acoustic pulse are studied. We have studied the dynamic behavior of the radial response of bubble using Lyapunov exponent spectra, bifurcation diagrams, time series and phase diagram. A period-doubling bifurcation structure is predicted to occur for certain values of the effects of parameters. The results show that by increasing the elasticity of the fluid, the growth phenomenon will be unstable. On the other hand, when the frequency of the external pulse increases the bubble growth experiences more stable condition. It is shown that the results are in good agreement with the previous studies.     

On the stability of an elastic column positioned on an elastic half space

Summary Stability of a heavy elastic column loaded by a concentrated force at the top is analysed. It is assumed that the column is fixed to a rigid circular plate that is positioned on a homogeneous, isotropic, linearly elastic half-space. The constitutive equations for the column are assumed in the form that allows axial compressibility and takes into account the influence of shear stresses. It is shown that eigenvalues of the linearized equations determine the bifurcation points of the full non-linear system of equilibrium equations. Also, the type of bifurcation at the lowest eigenvalue is examined and shown that it could be both super-and sub-critical. The post-critical shape of the column is determined by numerical integration of the equilibrium equations. Received 13 June 1998; accepted for publication 12 November 1998     

Spiral vortices between concentric cylinders

Spiral vortices appearing in Couette-Taylor flows are studied by means of numerical simulation. Transition curves from Couette to spiral vortices for different radius ratios and wavenumbers have been calculated in order to test our technique. Critical Reynolds numbers, angular velocities and slopes of the spirals at the onset of the instability agree with previous results [1]. Non-linear solutions obtained by a pseudospectral collocation method are studied, and they show a weak net axial flow. In order to counteract this effect, which is absent in the usual experimental set-up, an axial pressure gradient has been included. This procedure has proved to be sufficient to make the axial flow negligible. The onset of a quasiperiodic flow for larger Reynolds numbers, corresponding to a secondary bifurcation is also presented.     

Dynamic bifurcations analysis of a micro rotating shaft considering non-classical theory and internal damping

In this study, the dynamic bifurcation of a viscoelastic micro rotating shaft is investigated. The non-classical theory (the modified couple stress theory) and the Kelvin Voigt model are used for modeling the viscoelastic micro shaft. The transverse equations of motion are derived using the variational approach. The reduced order model of the system is obtained by the Galerkin method. Using the Routh–Hurwitz criteria the stability regions of the system are extracted in which the effect of the length scale parameter is significant. Using the center manifold theory and the normal form method the double zero eigenvalue bifurcation is analyzed. The results show that the internal and external damping coefficients, the rotational speed and the material length scale parameter influence the critical speed, amplitude, and phase of a non-trivial solution, and radius of limit cycle (periodic solution). Also, it is seen that by increasing the dimensionless length scale parameter (material length scale per radius of the shaft) the radius of the limit cycle is decreased, whereas the critical rotational speed and the rate of the phase are increased. However, the radius of the limit cycle concerning the classical theory is higher than that of regarding the modified couple stress theory. Furthermore, with an increase of the external damping coefficient the radius of the limit cycle is linearly decreased; however, the critical speed of the system is increased. Additionally, by decreasing length scale parameter the results of the modified couple stress theory approach the classical theory ones.     

Bifurcation of a running crack in antiplane strain

An elastodynamic explanation of running crack bifurcation is explored. The geometry is a semi-infinite body in a state of antiplane strain, which contains a two-dimensional edge crack. It is assumed that a quasi-static increase of the external loads gives rise to rapid crack propagation at time t = 0, with an arbitrary and time-varying speed, but in the plane of the crack. A short time later the crack is assumed to bifurcate at angles −κπ and +gkπ, and with velocities v. The elastodynamic intensity factors are computed, and the balance of rates of energies is employed to discuss the conditions for bifurcation.     

Periodic Delay Effects on Cutting Dynamics

We consider a delay equation whose delay is perturbed by a small periodic fluctuation. In particular, it is assumed that the delay equation exhibits a Hopf bifurcation when its delay is unperturbed. The periodically perturbed system exhibits more delicate bifurcations than a Hopf bifurcation. We show that these bifurcations are well explained by the Bogdanov-Takens bifurcation when the ratio between the frequencies of the periodic solution of the unperturbed system (Hopf bifurcation) and the external periodic perturbation is 1:2. Our analysis is based on center manifold reduction theory.     

DERIVATION OF NONLINEAR WAVE EQUATION FOR FLEXURAL MOTIONS OF AN ELASTIC BEAM TRAVELLING IN AN AIR-FILLED TUBE

Asymptotic analysis is carried out to derive a nonlinear wave equation for flexural motions of an elastic beam of circular cross-section travelling along the centre-axis of an air-filled, circular tube placed coaxially. Both the beam and tube are assumed to be long enough for end-effects to be ignored and the aerodynamic loading on the lateral surface of the beam is considered. Assuming a compressible inviscid fluid, the velocity potential of the air is sought systematically in the form of power series in terms of the ratios of the tube radius to a wavelength and of a typical deflection to the radius. Evaluating the pressure force acting on the lateral surface of the beam, the aerodynamic loading including the effects of finite deflection as well as of air's compressibility and axial curvature of the beam are obtained. Although the nonlinearity arises from the kinematical condition on the beam surface, it may be attributed to the presence of the tube wall. With the aerodynamic loading thus obtained, a nonlinear wave equation is derived, whereas linear theory is assumed for the flexural motions of the beam. Some discussions are given on the results.