Non-linear in-plane buckling of rotationally restrained shallow arches under a central concentrated load

This paper investigates the non-linear in-plane buckling of pin-ended shallow circular arches with elastic end rotational restraints under a central concentrated load. A virtual work method is used to establish both the non-linear equilibrium equations and the buckling equilibrium equations. Analytical solutions for the non-linear in-plane symmetric snap-through and antisymmetric bifurcation buckling loads are obtained. It is found that the effects of the stiffness of the end rotational restraints on the buckling loads, and on the buckling and postbuckling behaviour of arches, are significant. The buckling loads increase with an increase of the stiffness of the rotational restraints. The values of the arch slenderness that delineate its snap-through and bifurcation buckling modes, and that define the conditions of buckling and of no buckling for the arch, increase with an increase of the stiffness of the rotational end restraints.     

Interactive buckling in long thin-walled rectangular hollow section struts

An analytical model describing the nonlinear interaction between global and local buckling modes in long thin-walled rectangular hollow section struts under pure compression founded on variational principles is presented. A system of nonlinear differential and integral equations subject to boundary conditions is formulated and solved using numerical continuation techniques. For the first time, the equilibrium behaviour of such struts with different cross-section joint rigidities is highlighted with characteristically unstable interactive buckling paths and a progressive change in the local buckling wavelength. With increasing joint rigidity within the cross-section, the severity of the unstable post-buckling behaviour is shown to be mollified. The results from the analytical model are validated using a nonlinear finite element model developed within the commercial package Abaqus and show excellent comparisons. A simplified method to calculate the local buckling load of the more compressed web undergoing global buckling and the corresponding global mode amplitude at the secondary bifurcation is also developed. Parametric studies on the effect of varying the length and cross-section aspect ratio are also presented that demonstrate the effectiveness of the currently developed models.     

Magnetoelastic buckling of a rectangular block in plane strain

Of interest here is the stability of a rectangular block subjected to a uniform magnetic field perpendicular to its longitudinal axis. The two ends of the block are frictionless and kept parallel to each other. This boundary value problem is motivated by the classical problem of magnetoelastic buckling in which a cantilever beam subjected to a transverse magnetic field buckles when the applied field reaches a critical value.This work presents a finite strain continuum mechanics formulation of the stability problem of a homogeneous, compressible, magnetoelastic rectangular block in plane strain subjected to a uniform transverse magnetic field. The applied variational approach employs an unconstrained energy minimization recently proposed by the authors.The analytical solution for the critical buckling fields for both the antisymmetric and symmetric modes are obtained for three different constitutive laws. The corresponding result for thin beams is extracted asymptotically for a special material and the solution is compared to previously published results. The critical magnetic field is shown to increase monotonically with the block's aspect ratio for each material and mode type. Antisymmetric modes are always the critical buckling modes for stress saturated and neo-Hookean materials, except for a narrow range of moderate aspect ratios (about 0.25) where symmetric modes become critical. For strain-saturated solids no buckling is possible above a maximum aspect ratio.     

Buckling analysis of pre-cracked beam-columns by Liapunov's second method

Interaction between reduced axial and lateral torsional buckling strengths of thin beam-columns of double symmetric sections with an edge crack at the middle is examined via Liapunov's second method. This is done in order to show how critical loads for pre-cracked thin-walled beam-columns can be obtained with out the rigour of meticulous solution of a set of differential equations. Having derived an asymptotic expression for the normalized energy release rates for the beam-column, a metric space and a Liapunov's Functional for a pre-cracked beam-column are proposed. By using Zubov's stability theorems and appropriate eigen value inequalities, expressions for reduced buckling strengths due to edge cracks are obtained. A crack reduction parameter k is proposed which is used to calculate reserve strengths for published lateral buckling strength equations. Results obtained using the proposed parameter are well in agreement with findings of earlier published works in the literature when the parameter is set to zero. The results also show that cracks cause loss of stiffness and can reduce the lateral buckling strengths of beam-columns considerably.     

Micro-buckling in the nanocomposite structure of biological materials

Nanocomposite structure, consisting of hard mineral and soft protein, is the elementary building block of biological materials, where the mineral crystals are arranged in a staggered manner in protein matrix. This special alignment of mineral is supposed to be crucial to the structural stability of the biological materials under compressive load, but the underlying mechanism is not yet clear. In this study, we performed analytical analysis on the buckling strength of the nanocomposite structure by explicitly considering the staggered alignment of the mineral crystals, as well as the coordination among the minerals during the buckling deformation. Two local buckling modes of the nanostructure were identified, i.e., the symmetric mode and anti-symmetric mode. We showed that the symmetric mode often happens at large aspect ratio and large volume fraction of mineral, while the anti-symmetric happens at small aspect ratio and small volume fraction. In addition, we showed that because of the coordination of minerals with the help of their staggered alignment, the buckling strength of these two modes approached to that of the ideally continuous fiber reinforced composites at large aspect ratio given by Rosen's model, insensitive to the existing “gap”-like flaws between mineral tips. Furthermore, we identified a mechanism of buckling mode transition from local to global buckling with increase of aspect ratio, which was attributed to the biphasic dependence of the buckling strength on the aspect ratio. That is, for small aspect ratio, the local buckling strength is smaller than that of global buckling so that it dominates the buckling behavior of the nanocomposite; for comparatively larger aspect ratio, the local buckling strength is higher than that of global buckling so that the global buckling dominates the buckling behavior. We also found that the hierarchical structure can effectively enhance the buckling strength, particularly, this structural design enables biological nanocomposites to avoid local buckling so as to achieve global buckling at macroscopic scales through hierarchical design. These features are remarkably important for the mechanical functions of biological materials, such as bone, teeth and nacre, which often sustain large compressive load.     

Influence of the connecting condition on the dynamic buckling of longitudinal impact for an elastic rod

The stress wave propagation law and dynamic buckling critical velocity are formulated and solved by considering a general axial connecting boundary for a slender elastic straight rod impacted by a rigid body. The influence of connecting stiffness on the critical velocity is investigated with varied impactor mass and buckling time. The influences of rod length and rod mass on the critical velocity are also discussed. It is found that greater connecting stiffness leads to larger stress amplitude, and further results in lower critical velocity. It is particularly noteworthy that when the connecting stiffness is less than a certain value,dynamic buckling only occurs before stress wave reflects off the connecting end. It is also shown that longer rod with larger slenderness ratio is easier to buckle, and the critical velocity for a larger-mass rod is higher than that for a lighter rod with the same geometry.     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

弹性地基上多孔功能梯度材料矩形板的自由振动与临界屈曲载荷分析

多孔功能梯度材料(FGM)构件的特性与孔隙率和孔隙分布形式有密切关系。本文基于经典板理论,考虑不同孔隙分布形式时修正的混合率模型,研究Winkler弹性地基上四边受压多孔FGM矩形板的自由振动与临界屈曲载荷特性。首先利用Hamilton原理和物理中面的定义推导Winkler弹性地基上四边受压多孔FGM矩形板自由振动的控制微分方程并进行无量纲化,然后应用微分变换法(DTM)对无量纲控制微分方程和边界条件进行变换,得到计算无量纲固有频率和临界屈曲载荷的代数特征方程。将问题退化为孔隙率为零时的FGM矩形板并与已有文献进行对比以验证其有效性。最后计算并分析了梯度指数、孔隙率、地基刚度系数、长宽比、四边受压载荷及边界条件对多孔FGM矩形板无量纲固有频率的影响以及各参数对无量纲临界屈曲载荷的影响。     

Buckling pressure of “pumpkin” balloons

This paper presents a computational study of the critical buckling pressure of pumpkin balloons, which consist of a thin, compliant membrane constrained by stiff meridional tendons. The n-fold symmetric shape of a pumpkin balloon with n identical lobes is exploited by adopting a symmetry-adapted coordinate system, which leads to the tangent stiffness matrix in an efficient block-diagonal form; the smallest eigenvalue of a particular block leads to the buckling pressure for the balloon. Two different types of balloon design are considered. Extensive results are obtained for the buckling pressures of a set of 10 m diameter experimental balloons and also for an 80 m diameter flight balloon. The key findings are as follows: the same type of buckling mode, forming four circumferential waves is critical for most of the balloons that have been analysed; balloons with flatter lobes are more stable, and the buckling pressure varies with an inverse power-law of the number of lobes; increasing the Young’s modulus, the Poisson’s ratio of the membrane, or the diameter of the end fitting has the effect of increasing the buckling pressure; but increasing the axial stiffness of the tendons has the effect of decreasing the buckling pressure.     

Stability analysis of functionally graded rectangular plates under nonlinearly varying in-plane loading resting on elastic foundation

In this research work, an exact analytical solution for buckling of functionally graded rectangular plates subjected to non-uniformly distributed in-plane loading acting on two opposite simply supported edges is developed. It is assumed that the plate rests on two-parameter elastic foundation and its material properties vary through the thickness of the plate as a power function. The neutral surface position for such plate is determined, and the classical plate theory based on exact neutral surface position is employed to derive the governing stability equations. Considering Levy-type solution, the buckling equation reduces to an ordinary differential equation with variable coefficients. An exact analytical solution is obtained for this equation in the form of power series using the method of Frobenius. By considering sufficient terms in power series, the critical buckling load of functionally graded plate with different boundary conditions is determined. The accuracy of presented results is verified by appropriate convergence study, and the results are checked with those available in related literature. Furthermore, the effects of power of functionally graded material, aspect ratio, foundation stiffness coefficients and in-plane loading configuration together with different combinations of boundary conditions on the critical buckling load of functionally graded rectangular thin plate are studied.     

Symmetry breaking in an initially curved pre-stressed micro beam loaded by a distributed electrostatic force

The symmetric and asymmetric buckling of an initially curved micro beam subjected to an axial pre-stressing load and transversal distributed electrostatic force is studied. The analysis is based on a reduced order (RO) model resulting from the Galerkin decomposition with buckling modes of a straight beam used as the base functions. The criteria of symmetric limit point buckling and of non-symmetric bifurcation are derived in terms of the geometric parameters of the beam and the axial load. Two symmetry breaking conditions, defining the relations between the axial load and the geometric parameters of beams for which an asymmetric response bifurcates from the symmetric one, are obtained. The necessary criterion establishes the conditions for the appearance of bifurcation points on the unstable branch of the symmetric response curve; the sufficient criterion assures a realistic asymmetric buckling bifurcating from the stable branches of the symmetric response curve. A comparison between the RO model results and those obtained by direct numerical analysis shows good agreement between the two and indicates that the obtained criteria can be used to predict symmetric and non-symmetric buckling in electrostatically actuated curved pre-stressed micro beams. It is shown that while the symmetry breaking conditions are affected by the nonlinearity of the electrostatic force, its influence is less pronounced than in the case of the symmetric snap-through criterion. The nature of the latter and the relations between it and the symmetry breaking criteria are found to go through a prominent qualitative change as the initial distance between the beam and the electrode, characterizing the electrostatic force, changes.     

Exact solutions for buckling of non-uniform columns under axial concentrated and distributed loading

In this paper, the buckling problem of non-uniform columns subjected to axial concentrated and distributed loading is studied. The expression for describing the distribution of flexural stiffness of a non-uniform column is arbitrary, and the distribution of axial forces acting on the column is expressed as a functional relation with the distribution of flexural stiffness and vice versa. The governing equation for buckling of a non-uniform column with arbitrary distribution of flexural stiffness or axial forces is reduced to a second-order differential equation without the first-order derivative by means of functional transformations. Then, this kind of differential equation is reduced to Bessel equations and other solvable equations for 12 cases, several of which are important in engineering practice. The exact solutions that represent a class of exact functional solutions for the buckling problem of non-uniform columns subjected to axial concentrated and distributed loading are obtained. In order to illustrate the proposed method, a numerical example is given in the last part of this paper.     

THE TRANSVERSE NONLINEAR FREE VIBRATION AND THE BUCKLING OF A COMPRESSIVE BAR ON AN ELASTIC FOUNDATION

基于Hamilton 原理,运用假设时间模态法,得到了弹性基础上压杆的横向非线性自由振动与屈曲的位移型常微分控制方程. 考虑一端固定另一端可移简支边界条件,采用打靶法得到了结构第一至第三阶结构频率与一阶屈曲载荷的数值结果. 结果表明：随轴心压力增加,结构频率减小;随弹性基础刚度增加,结构频率与屈曲载荷均增加;弹性基础刚度对结构频率的影响随振型阶数增加在减小;在小振幅的情形下,不同振型对一阶屈曲载荷的影响很小.     

Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

The purpose of the present work is to study the buckling problem with plate/shell topology optimiza-tion of orthotropic material.A model of buckling topology optimization is established based on the independent,con-tinuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function(PF)as filter functions are introduced to recognize the element mass,the element stiffness matrix,and the ele-ment geometric stiffness matrix.The filter functions of the orthotropic material stiffness are deduced. Then these fil-ter functions are put into buckling topology optimization of a differential equation to analyze the design sensitiv-ity.Furthermore,the buckling constraints are approximately expressed as explicit functions with respect to the design vari-ables based on the first-order Taylor expansion.The objective function is standardized based on the second-order Taylor expansion. Therefore,the optimization model is translated into a quadratic program.Finally,the dual sequence quadratic programming(DSQP)algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions(CEF and PF)are applied to solve the opti-mal model.Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.     

Postbuckling of long unsymmetrically laminated composite plates under axial compression

An approximate solution is given for the postbuckling of infinitely long and unsymmetrically laminated composite plates. This solution is obtained by superposing a polynomial transverse displacement given by bending due to unsymmetric laminate configurations and a simple functional representation for the buckling mode in conjunction with the Galerkin method. Nondimensional parameters are used to express the approximate solution in a very simple and clear formulation. The results given by this solution for axial compression in the longitudinal direction are compared with the results given by the nonlinear finite element method (FEM) for finite length rectangular long plates. The influence of the boundary conditions on postbuckling response is also studied. For the FEM analysis, two different simply supported boundary conditions on the long edges of the plate are considered. It is found that these two sets of boundary conditions give different results for the buckling and postbuckling finite element analysis. In most cases the FEM analysis overestimate and, respectively, underestimate the approximate closed form solution, depending on the type of simply supported boundary condition considered. Thus, the approximate solution appears useful for design purposes as an averaged quantity between the two FEM analyses. Also, it is found that the reduced bending stiffness method can be successfully used for determining the approximate solution.     

Finite element analysis of post-buckling dynamics in plates—Part I: An asymptotic approach

Various static and dynamic aspects of post-buckled thin plates, including the transition of buckled patterns, post-buckling dynamics, secondary bifurcation, and dynamic snapping (mode jumping phenomenon), are investigated systematically using asymptotical and non-stationary finite element methods. In part I, the secondary dynamic instability and the local post-secondary buckling behavior of thin rectangular plates under generalized (mechanical and thermal) loading is investigated using an asymptotic numerical method which combines Koiter’s nonlinear instability theory with the finite element technique. A dynamic multi-mode reduction method—similar to its static single-mode counterpart: Liapunov–Schmidt reduction—is developed in this perturbation approach. Post-secondary buckling equilibrium branches are obtained by solving the reduced low-dimensional parametric equations and their stability properties are determined directly by checking the eigenvalues of the resulting Jacobian matrix. Typical post-secondary buckling forms—transcritical, supercritical and subcritical bifurcations are observed according to different combinations of boundary conditions and load types. Geometric imperfection analysis shows that not only the secondary bifurcation load but also changes in the fundamental post-secondary buckling behavior are affected. The post-buckling dynamics and the global analysis of mode jumping of the plates are addressed in part II.     

STRESS-STRAIN STATE OF ELASTOPLASTIC BODIES WITH CRACK

A physical cut model is used to describe the changes in the stress-strain state(SSS)in elastoplastic bodies weakened by cracks. The distance between the crack edges is considered to be finite in contrast to the mathematical cut. The interactive layer with a thickness limited by the possibility of using the hypothesis of continuity is distinguished on the physical cut extension.Distribution of stresses and strains over the layer thickness is constant and does not depend on the geometry of the boundary between the cut and the interactive layer. The relationship between stresses and strains is determined by the deformation plasticity theory. The problem of plane strain or plane stress state of an arbitrary finite body weakened by a physical cut is reduced to solving a system of two variational equations for displacement fields in the body parts adjacent to the interactive layer. The proposed approach eliminates the singularity in stress distribution in contrast to the mathematical cut model. Use of local strength criteria allows us to determine the time, point and direction of the fracture initiation. Possibilities of the proposed model are illustrated by solving the problems of determining the SSS of a rectangular body weakened by a physical cut under symmetric and antisymmetric loadings.     

Buckling load optimization of beams

In this paper, shape optimization is used to optimize the buckling load of a Euler–Bernoulli beam having constant volume. This is achieved by varying appropriately the beam cross section so that the beam buckles with the maximum or a prescribed buckling load. The problem is reduced to a nonlinear optimization problem under equality and inequality constraints as well as specified lower and upper bounds. The evaluation of the objective function requires the solution of the buckling problem of a beam with variable stiffness subjected to an axial force. This problem is solved using the analog equation method for the fourth-order ordinary differential equation with variable coefficients. Besides its accuracy, this method overcomes the shortcomings of a possible FEM solution, which would require resizing of the elements and recomputation of their stiffness properties during the optimization process. Several example problems are presented that illustrate the method and demonstrate its efficiency.     

ANALYSIS ON THE MAGNETO-ELASTIC-PLASTIC BUCKLING/SNAPPING OF CANTILEVER RECTANGULAR FERROMAGNETIC PLATES

An analysis of buckling/snapping and bending behaviors of magneto-elastic-plastic interaction and coupling for cantilever rectangular soft ferromagnetic plates is presented.Based on the expression of magnetic force from the variational principle of ferromagnetic plates,the buckling and bending theory of thin plates,the Mises yield criterion and the increment theory for plastic deformation,we establish a numerical code to quantitatively simulate the behaviors of the nonlinearly multi-fields coupling problems by the finite element method.Along with the phenom- ena of buckling/snapping and bending,or the characteristic curve of deflection versus magnitude of applied magnetic fields being numerically displayed,the critical loads of buckling/snapping, and the influences of plastic deformation and the width of plate on these critical loads,the plastic regions expanding with the magnitude of applied magnetic field,as well as the evolvement of deflection configuration of the plate are numerically obtained in a case study.     

U型波纹管及相关结构环向屈曲的有限元分析(Ⅰ)--基本方程及环板的屈曲

U型波纹管是现代管道系统中最常见的一种位移补偿器 ,它由环板和具有正、负Gauss曲率的半圆环壳组成 ,在管道所传输的介质的压力作用下会发生屈曲。其中环向屈曲最为复杂 ,精确的理论分析非常困难 ,有限元分析也不多见。作者在分析前人工作的基础上 ,以圆环壳段为单元 (特定的旋转壳段单元 ,能自动退化成环板单元 ) ,限于弹性范围和线性化特征值问题 ,对介质压力作用下U型波纹管及其相关结构 (圆环板、圆环壳、半圆环壳 )的环向屈曲问题进行了分析。考虑了结构屈曲前的弯曲 ,计及压力的二次势能 ,导出的应力刚度矩阵和载荷刚度矩阵是非对称的。全部工作分为三部分 :(Ⅰ )基本方程 ,环板的屈曲 ;(Ⅱ )圆环壳、半圆环壳的屈曲 ;(Ⅲ )波纹管平面失稳的机理。本文为第一部分 ,除推导公式外 ,对不同边界和不同内外径之比的环板在径向均匀压力作用下的环向屈曲进行了计算 (轴对称的径向屈曲作为特例得到 ) ,给出了前屈曲应力分布、临界载荷及相应的屈曲模态 ,并将临界压力的值与前人基于vonK偄ｒｍ偄ｎ大挠度板的精确解进行了比较 ,吻合良好。