Yield anisotropy effects on buckling of circular tubes under bending

Relatively thin-walled tubes bent into the plastic range buckle by axial wrinkling. The wrinkles initially grow stably but eventually localize and cause catastrophic failure in the form of sharp local kinking. The onset of axial wrinkling was previously established by bifurcation analyses that use instantaneous deformation theory moduli. The curvatures at bifurcation were predicted accurately, but the wrinkle wavelengths were consistently longer than measured values. The subject is revisited with the aim of resolving this discrepancy. A set of new bending experiments is conducted on aluminum alloy tubes. The results are shown to be in line with previous ones. However, the tubes used were found to exhibit plastic anisotropy, which was measured and characterized through Hill’s quadratic anisotropic yield function. The anisotropy was incorporated in the flow theory used for prebuckling and postbuckling calculations as well as in the deformation theory used for bifurcation checks. With the anisotropy accounted for, calculated tube responses are found to be in excellent agreement with the measured ones while the predicted bifurcation curvatures and wrinkle wavelengths fall in line with the measurements also. The postbuckling response is established using a finite element model of a tube assigned an initial axisymmetric imperfection with the calculated wavelength. The response develops a limit moment that is followed by a sharp kink that grows while the overall moment drops. The curvature at the limit moment agrees well with the experimental onset of failure. From parametric studies of the various instabilities it is concluded that, for optimum predictions, anisotropy must be incorporated in both bifurcation buckling as well as in postbuckling calculations.     

Buckling and post-buckling of long pressurized elastic thin-walled tubes under in-plane bending

The present paper focuses on the structural stability of long uniformly pressurized thin elastic tubular shells subjected to in-plane bending. Using a special-purpose non-linear finite element technique, bifurcation on the pre-buckling ovalization equilibrium path is detected, and the post-buckling path is traced. Furthermore, the influence of pressure (internal and/or external) as well as the effects of radius-to-thickness ratio, initial curvature and initial ovality on the bifurcation moment, curvature and the corresponding wavelength, are examined. The local character of buckling in the circumferential direction is also demonstrated, especially for thin-walled tubes. This observation motivates the development of a simplified analytical formulation for tube bifurcation, which considers the presence of pressure, initial curvature and ovality, and results in closed-form expressions of very good accuracy, for tubes with relatively small initial curvature. Finally, aspects of tube bifurcation are illustrated using a simple mechanical model, which considers the ovalized pre-buckling state and the effects of pressure.     

热环境中功能梯度圆柱薄壳分岔屈曲的边界约束效应

基于前屈曲一致理论,研究了热环境中受轴压功能梯度圆柱薄壳分岔屈曲的边界约束效应问题.导出前屈曲变形的解析解,结合分离变量法与有限差分法求解分岔屈曲控制方程,由此导出确定临界轴压的非线性特征值问题.考虑材料热物性质与温度的相关性,分别就两端简支和两端固支边界条件,分析了温度梯度、初始几何缺陷、组分材料体积分数对分岔屈曲临...     

Long-term non-linear behaviour and buckling of shallow concrete-filled steel tubular arches

This paper presents a theoretical analysis for the long-term non-linear elastic in-plane behaviour and buckling of shallow concrete-filled steel tubular (CFST) arches. It is known that an elastic shallow arch does not buckle under a load that is lower than the critical loads for its bifurcation or limit point buckling because its buckling equilibrium configuration cannot be achieved, and the arch is in a stable equilibrium state although its structural response may be quite non-linear under the load. However, for a CFST arch under a sustained load, the visco-elastic effects of creep and shrinkage of the concrete core produce significant long-term increases in the deformations and bending moments and subsequently lead to a time-dependent change of its equilibrium configuration. Accordingly, the bifurcation point and limit point of the time-dependent equilibrium path and the corresponding buckling loads of CFST arches also change with time. When the changing time-dependent bifurcation or limit point buckling load of a CFST arch becomes equal to the sustained load, the arch may buckle in a bifurcation mode or in a limit point mode in the time domain. A virtual work method is used in the paper to investigate bifurcation and limit point buckling of shallow circular CFST arches that are subjected to a sustained uniform radial load. The algebraically tractable age-adjusted effective modulus method is used to model the time-dependent behaviour of the concrete core, based on which solutions for the prebuckling structural life time corresponding to non-linear bifurcation and limit point buckling are derived.     

Bending instabilities of elastic tubes

The present paper examines instabilities of long thin elastic tubes. Both initially straight and initially bent tubes are analyzed under in-plane bending. Tube response, a combination of ovalization instability and bifurcation instability (buckling), is investigated using a nonlinear finite element (FE) technique, which employs polynomial functions in the longitudinal tube direction and trigonometric functions to describe cross-sectional deformation. It is demonstrated that the interaction between the two instability modes depends on the value and the sign of the initial tube curvature. The ovalization of initially bent tubes is examined in detail and, in particular, the case of opening moments. Furthermore, the paper emphasizes on bifurcation instability. It is shown that buckling may occur prior to or beyond the ovalization limit point, depending on the value of the initial curvature. Using the nonlinear FE formulation, the location of bifurcation on the primary path is detected, post-buckling equilibrium paths are traced, and the corresponding wavelengths of the buckled configurations are calculated. Moreover, results over a wide range of initial curvature values are presented, extending the findings of previous works. Finally, several analytical approaches, introduced in previous research works, are also employed to estimate the moments causing ovalization and bifurcation instability. These approaches are based on nonlinear flexible shell theory or simplified ring analysis. The efficiency and accuracy of those analytical methods with respect to the nonlinear FE formulation are examined.     

Comparisons of test and theory for nonsymmetric elastic-plastic buckling of shells of revolution

Experimental and analytical buckling pressures are presented for very carefully fabricated thin cylindrical shells with 45, 60 and 75° conical heads and for cylindrical shells with torispherical heads pierced by axisymmetric cylindrical nozzles of various thicknesses and diameters. Nonsymmetric buckling occurs at pressures for which some of the material is loading plastically in the neighborhoods of stress concentrations caused by meridional slope discontinuities. The buckling pressures for the cone-cylinder vessels are predicted within 2.6 per cent and for the pierced torispherical vessels within 4.4 per cent with use of BOSOR5, a computer program based on the finite difference energy method in which axisymmetric large deflections, nonlinear material properties and nonsymmetric bifurcation buckling are accounted for. The predicted buckling pressures of the pierced torispherical specimens are rather sensitive to details of the analytical model in the neighborhood of the juncture between the nozzle and the head. The buckling pressures of the cone-cylinder vessels can be accurately predicted by treatment of the wall material as elastic, enforcement of the full compatibility conditions at the juncture in the prebuckling analysis, and release of the rotation compatibility condition in the bifurcation (eigenvalue) analysis.     

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).     

Wrinkling of tubes in bending from finite strain three-dimensional continuum theory

The problem of a tube under pure bending is first solved as a generalised plane strain problem. This then provides the prebifurcation solution, which is uniform along the length of the tube. The onset of wrinkling is then predicted by introducing buckling modes involving a sinusoidal variation of the displacements along the length of the tube. Both the prebuckling analysis and the bifurcation check require only a two-dimensional finite element discretisation of the cross-section with special elements. The formulation does not rely on any of the approximations of a shell theory, or small strains. The same elements can be used for pure bending and local buckling a prismatic beam of arbitrary cross-section. Here the flow theory of plasticity with isotropic hardening is used for the prebuckling solution, but the bifurcation check is based on the incremental moduli of a finite strain deformation theory of plasticity.For tubes under pure bending, the results for limit point collapse (due to ovalisation) and bifurcation buckling (wrinkling) are compared to existing analysis and test results, to see whether removing the approximations of a shell theory and small strains (used in the existing analyses) leads to a better prediction of the experimental results. The small strain analysis results depend on whether the true or nominal stress–strain curve is used. By comparing small and finite strain analysis results it is found that the small strain approximation is good if one uses (a) the nominal stress–strain curve in compression to predict bifurcation buckling (wrinkling), and (b) the true stress–strain curve to calculate the limit point collapse curvature.In regard to the shell theory approximations, it is found that the three-dimensional continuum theory predicts slightly shorter critical wrinkling wavelengths, especially for lower diameter-to-wall-thickness (D/t) ratios. However this difference is not sufficient to account for the significantly lower wavelengths observed in the tests.     

Numerical investigation of fluid–particle interactions for embolic stroke

Roughly one-third of all strokes are caused by an embolus traveling to a cerebral artery and blocking blood flow in the brain. The objective of this study is to gain a detailed understanding of the dynamics of embolic particles within arteries. Patient computed tomography image is used to construct a three-dimensional model of the carotid bifurcation. An idealized carotid bifurcation model of same vessel diameters was also constructed for comparison. Blood flow velocities and embolic particle trajectories are resolved using a coupled Euler–Lagrange approach. Blood is modeled as a Newtonian fluid, discretized using the finite volume method, with physiologically appropriate inflow and outflow boundary conditions. The embolus trajectory is modeled using Lagrangian particle equations accounting for embolus interaction with blood as well as vessel wall. Both one- and two-way fluid–particle coupling are considered, the latter being implemented using momentum sources augmented to the discretized flow equations. It was observed that for small-to-moderate particle sizes (relative to vessel diameters), the estimated particle distribution ratio—with and without the inclusion of two-way fluid–particle momentum exchange—were found to be similar. The maximum observed differences in distribution ratio with and without the coupling were found to be higher for the idealized bifurcation model. Additionally, the distribution was found to be reasonably matching the volumetric flow distribution for the idealized model, while a notable deviation from volumetric flow was observed in the anatomical model. It was also observed from an analysis of particle path lines that particle interaction with helical flow, characteristic of anatomical vasculature models, could play a prominent role in transport of embolic particle. The results indicate therefore that flow helicity could be an important hemodynamic indicator for analysis of embolus particle transport. Additionally, in the presence of helical flow, and vessel curvature, inclusion of two-way momentum exchange was found to have a secondary effect for transporting small to moderate embolus particles—and one-way coupling could be used as a reasonable approximation, thereby causing substantial savings in computational resources.     

Buckling of elastic-plastic shells of revolution with discrete elastic-plastic ring stiffeners

The theory is summarized for axisymmetric prebuckling and nonsymmetric bifurcation buckling of ring-stiffened shells of revolution. The analysis is based on finite difference energy minimization in which moderately large meridional rotations, elastic-plastic effects, and primary or secondary creep are included. This theory is implemented in a computer program called BOSOR5, for the analysis of segmented and branched ring-stiffened shells of revolution of multi-material construction.Comparisons between test and theory are given for axisymmetric collapse and nonsymmetric bifurcation buckling of 69 machined ring-stiffened aluminum cylinders submitted to external hydrostatic pressure. Because most of the cylinders fail at an average stress which corresponds to the knee of the stress-strain curve, the analytical predictions are not very sensitive to modeling particulars such as nodal point density or boundary conditions. Agreement between test and theory is improved if the analytical model reflects the fact that the shell and rings intersect over finite axial lenths.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

Bifurcation buckling of shells of revolution including large deflections,plasticity and creep

A summary is first presented of the conceptual difficulties and paradoxes surrounding plastic bifurcation buckling analysis. Briefly discussed are nonconservativeness, loading rate during buckling, and the discrepancy of buckling predictions with use of J₂ flow theory vs J₂ deformation theory. The axisymmetric prebuckling analysis, including large deflections, elastic-plastic material behavior and creep is summarized. Details are given on the analysis of nonsymmetric bifurcation from the deformed axisymmetric state. Both J₂ flow theory and J₂ deformation theory are described. The treatment, based on the finite-difference energy method, applies to layered segmented and branched shells of arbitrary meridional shape composed of a number of different elastic-plastic materials. Numerical results generated with a computer program based on the analysis are presented for an externally pressurized cylinder with conical heads.     

Mechanical behavior and wrinkling of lined pipes

The paper focuses on wrinkling of lined pipes (sometimes referred to as clad pipes) under bending loading, where a corrosion-resistant thin-walled liner is fitted inside a carbon–steel outer pipe. The problem is solved numerically, using nonlinear finite elements to simulate liner pipe deformation and its interaction with the outer pipe. Stresses and strains are monitored throughout the deformation stage, detecting possible detachment of the liner from the outer pipe and the formation of wrinkles. The wrinkling behavior of elastic and elastic–plastic (steel) lined pipes under bending is examined. The results indicate that the lateral confinement of the liner pipe due to the deformable outer pipe and its interaction with the outer pipe has a decisive influence on the wrinkling behavior of the lined pipe. It is also shown that the behavior is characterized by a first bifurcation in a uniform wrinkling pattern, followed by a secondary bifurcation. The values of corresponding buckling curvature are determined and comparison with available experimental results is conducted in terms of wrinkle height development and the corresponding buckling wavelength. The results of the present research can be used for safer design of lined pipes in pipeline applications.     

Electro-mechanical coupling bifurcation and bulging propagation in a cylindrical dielectric elastomer tube

This paper explores the critical and post-bulging bifurcation of a cylindrical dielectric elastomer (DE) tube undergoing finite deformation under electro-mechanical coupling loading. Explicit expressions for the critical conditions of electro-mechanical bifurcation are derived by using a simplified mathematical method. The post-bifurcation path is comprehensively investigated by specifying the material model as ideal dielectric elastomer. In the post-bifurcation analysis, we analytically establish conditions for the phase coexistence of steady propagation and analyze the physical implications. We demonstrate a global instability under force or voltage control and a localized instability under volume or charge control. Cylindrical tube experiments have been carried out under electro-mechanical coupling loading to verify the theoretical predictions. Good agreements on the critical conditions as well as the post-bifurcation path are obtained. This work characterizes the bifurcation mechanism of rubber-like materials under complex coupling loading.     

Postbuckling behavior of circular rings with two or four concentrated loads

Exact solutions are obtained for the deformations of circular rings which have buckled under two or four equal and equally spaced normal or centrally directed concentrated loads. The ring is assumed to be inextensible and the local moment to be proportional to the local change of curvature. The results indicate that the ring with four loads is neutrally stable at its bifurcation load while the ring with two loads is unstable.     

Plane strain dynamics of elastic solids with intrinsic boundary elasticity, with application to surface wave propagation

In this paper, in a development of the static theory derived by Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853), we establish the equations of motion for a non-linearly elastic body in plane strain with an elastic surface coating on part or all of its boundary. The equations of (linearized) incremental motions superposed on a finite static deformation are then obtained and applied to the problem of (time-harmonic) surface wave propagation on a pre-stressed incompressible isotropic elastic half-space with a thin coating on its plane boundary. The secular equation for (dispersive) wave speeds is then obtained in respect of a general form of incompressible isotropic elastic strain-energy function for the bulk material and a general energy function for the coating material. Specialization of the form of strain-energy function enables the secular equation to be cast as a quartic equation and we therefore focus on this for illustrative purposes. An explicit form for the secular equation is thereby obtained. This involves a number of material parameters, including residual stress and moment in the properties of the coating. It is shown how this equation relates to previous work on waves in a half-space with an overlying thin layer set in the classical theory of isotropic elasticity and, in particular, the significant effect of omission of the rotatory inertia term, even at small wave numbers, is emphasized. Corresponding results for a membrane-type coating, for which the bending moment, inertia and residual moment terms are absent, are also obtained. Asymptotic formulas for the wave speed at large wave number (high frequency) are derived and it is shown how these results influence the character of the wave speed throughout the range of wave number values. A bifurcation criterion is obtained from the secular equation by setting the wave speed to zero, thereby generalizing the bifurcation results of Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853) to the situation in which residual stress and moment are present in the coating. Numerical results which show the dependence of the wave speed on the various material parameters and the finite deformation are then described graphically. In particular, features which differ from those arising in the classical theory are highlighted.     

Hysteretic behavior of rectangular tube (box) sections based on Preisach model

In the present paper, the Preisach model of hysteresis is applied to model cyclic behavior of elasto-plastic material. Rate of loading and viscous effects will not be considered. The problem of axial loading of rectangular cross-section and cyclic bending of rectangular tube (box) will be studied in details. Hysteretic stress–strain loop for prescribed history of stress change is plotted for material modeled by series connection of three unite element. Also moment–curvature hysteretic loop is obtained for a prescribed curvature change of rectangular tube (box). All obtained results clearly show advantages of the Preisach model for describing cyclic behavior of elasto-plastic material.     

鳗鲡模式游动的鱼体变形曲率波的传播特性研究

将鳗鲡模式游动的七鳃鳗简化成材料性质均匀的变截面黏 弹性梁,通过数值方法求解鱼体在主动弯矩波(作为激励的驱动波)的驱动下匀速游动时身体变形曲率波的传播特性. 结果表明,当主动弯矩 的驱动频率高于鱼体结构基频时,可以观察到曲率波相对于驱动波存在相位滞后,且越靠近尾部滞后现象越明显,这意味着曲率波的波速 小于驱动波的波速,也间接地验证了前人的实验结果. 通过参数研究发现,鱼体变形曲率波与驱动波的波速比与表征流体黏性作用的雷诺数无关,而与表征驱动波和鱼体材料属性的 无量纲激励频率、激励波长及鱼体黏性系数有关. 对于鳗鲡模式游动的鱼类,曲率波与驱动波的波速比随着无量纲激励频率和波长的增大而降低,随着鱼体黏性系数的增大而增大. 进一步研究发现,通过小扰动分析得到的组合相似性参数$\varPi$可以统一描述波速比与激励参数、材料参数之间的关系.      

Mircopolar Model of Fracture for Composite Material

Failure of a composite is a complex process accompanied by irreversible changes in the microstructure of the material. Microscopic mechanisms are known of the accumulation of damage and failure of the type of localized and multiple ruptures of the fibers delamination along interphase boundaries, and also mechanisms associated with fracture of fibers. In this work, we propose a mathematical model of the local mechanism of failure of a composite material randomly reinforced with a system of short fibers. We implement the Cosserat moment model of crack tip for filament material, reinforced with whiskers or in fiber- reinforced polycrystalline materials. It is assumed that the angular distribution of the fibers is isotropic and the elastic characteristics of the fibers are considerably higher than the elastic constants of the matrix. We implement the homogenization procedure for the effective Cosserat constants similarly to the effective elastic constants. The singular solution in the vicinity of the crack tip in the Cosserat moment model is found. Using this solution, we examine the bending stresses in the filaments due to effective moment stresses in the material. The constructed model describes the phenomenon of fracture of the fibers occurring during crack propagation in those composites. The following assumptions are used as the main hypotheses for the micromechanical model. The matrix contains a nucleation crack. When the load is increased the crack grows and its boundary comes into contact with the reinforcing fibers. A further increase of the stress causes bending of the fiber. When~the fiber curvature reaches a specific critical value, the fiber ruptures. If the stress at infinity is given, the fibers no longer delay the development of failure during crack propagation The degree of bending distortion of the fiber in the vicinity of the boundary of the crack is determined by the moment model of the material. The necessity to take into account the moment stresses in the failure theory of the reinforced material was stressed in [Muki and Sternberg (1965) Zeitschrift f angew Math und Phys 16:611–615; Garajeu and Soos (2003) Math Mech Solids 8(2):189–218; Ostoja-Starzewski et al (1999) Mech Res Commun 26:387–396]. The moment Cosserat stresses were accounted also for inhomogeneous biomechanical materials by Buechner and Lakes (2003) Bio Mech Model Mechanobiol 1: 295–301. We should also mention the important methodological studies [Sternberg and Muki (1967) J Solids Struct 1:69–95; Atkinson and Leppington (1977) Int J Solids Struct 13: 1103–1122] concerned with the moment stresses in homogeneous fracture mechanics.     

CAVITATED BIFURCATION FOR INCOMPRESSIBLE HYPERELASTIC MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＣａｖｉｔａｔｉｏｎｐｈｅｎｏｍｅｎｏｎ ,ｔｈｅｓｕｄｄｅｎｆｏｒｍａｔｉｏｎａｎｄｇｒｏｗｔｈｏｆｖｏｉｄｓｉｎｓｏｌｉｄｍａｔｅｒｉａｌｓ,ｈａｓｌｏｎｇａｔｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎｂｅｃａｕｓｅｏｆｔｈｅｆｕｎｄａｍｅｎｔａｌｒｏｌｅｉｔｐｌａｙｓｏｎｔｈｅｌｏｃａｌｆａｉｌｕｒｅａｎｄｆｒａｃｔｕｒｅｏｆｍａｔｅｒｉａｌｓ.Ｔｈｅａｐｐｅａｒａｎｃｅｏｆｖｏｉｄｉｎｔｈｅｉｎｔｅｒｉｏｒｏｆｒｕｂｂｅｒｏｒｒｕｂｂｅｒ＿ｌｉｋｅｍａｔｅｒｉａｌｓｗａ…