Lower bounds to fundamental frequencies and buckling loads of columns and plates

A general derivation of expressions for lower bounds to fundamental frequencies and buckling loads is given for the class of structures governed by linear elastic theory in the prebuckling state. These expressions involve two Rayleigh quotients both of which are upper bounds for the fundamental frequency under a prescribed load. The displacement trial functions must satisfy force and kinematic continuity but no other conditions are required. Thus, if appropriate high order base functions are used, the finite element procedure can be used to systematically narrow the difference between the upper and lower bounds.The theory is illustrated with several column and plate problems. The finite element method is applied to uniform and nonuniform columns with a representative set of boundary conditions. Elementary trial functions are used to show that reasonable bounds can also be obtained for plates subjected to known states of stress. Since the lower bound is obtained with a variation of the classical technique of Rayleigh, these results indicate that the method may be suitable for conservatively estimating buckling loads and fundamental frequencies of engineering structures.     

A method for the estimation of plastic buckling loads

Certain elastic-plastic buckling problems require the solution of an appropriate incremental or “rate” boundary-value problem in order that physically meaningful results may be obtained. In this paper, it is shown that a recent general variational theorem by Neale[8] may be advantageous for the approximate solution of such problems. As an example, the buckling of elastic-plastic cylindrical shells under torsion is analyzed, wherein the material is assumed to obey the incremental theory of plasticity and the effects of initial imperfections in geometry are taken into account.     

有限长圆柱壳受径向冲击的塑性动力屈曲分析

讨论两端固定的圆柱薄壳在均布外部冲击下的塑性动力屈曲。依据实验现象对位移场作出假设，用扰动法求得屈曲的扰动控制方程组，降阶后用标准的龙格库塔法进行数值求解。其中物理关系采用Ｌｅｖｙ－Ｍｉｓｅｓ流动理论，材料为刚线性强化模型。计算结果同Ｆｌｏｒｅｎｃｅ等的Ｂｅｓｓｅｌ函数解作了比较，对边界对屈曲的影响作了讨论。本文的结果对潜艇抗水下爆炸的研究有重要意义。     

Bifurcation and post-buckling analysis of bimodal optimum columns

A mathematical formulation of column optimization problems allowing for bimodal optimum buckling loads is developed in this paper. The columns are continuous and linearly elastic, and assumed to have no geometrical imperfections. It is first shown that bimodal solutions exist for columns that rest on a linearly elastic (Winkler) foundation and have clamped-clamped and clamped-simply supported ends. The equilibrium equation for a non-extensible, geometrically nonlinear elastic column is then derived, and the initial post-buckling behaviour of a bimodal optimum column near the bifurcation point is studied using a perturbation method. It is shown that in the general case the post-buckling behaviour is governed by a fourth order polynomial equation, i.e., near the bifurcation point there may be up to four post-buckling equilibrium states emanating from the trivial equilibrium state. Each of these equilibrium states may be either supercritical or subcritical in the vicinity of the bifurcation point. The conditions for stability of these non-trivial post-buckling states are established based on verification of positive semi-definiteness of a two-by-two matrix whose coefficients are integrals of the buckling modes and their derivatives. In the end of the paper we present and discuss numerical results for the post-buckling behaviour of several columns with bimodal optimum buckling loads.     

黏弹性板的大挠度蠕变屈曲

研究了具有初始小挠度受轴向压载黏弹性板的蠕变屈曲问题,在建立控制方程时,利用了von Karman非线性应变-位移关系,并考虑了初始挠度,用标准线性固体模型描述材料的黏弹性特性,在求解非线性积分方程时,利用梯形公式计算记忆积分式,将非线性积分方程化为非线性代数方程进行数值求解,得到了结构的蠕变变形过程,又将问题退化到小挠度情况进行研究,得到了挠度随时间扩展的解析解,分析了瞬时失稳临界载荷、持久临界载荷的物理意义,讨论了考虑几何非线性对黏弹性板蠕变屈曲的影响。     

Predicting buckling loads from vibration data

There are analytical methods for predicting the buckling loads of columns with the boundaries ideally fixed, i.e., simply supported or built-in, or partially fixed. Vibration-test results may furnish a practical method of measuring the fixity. In this investigation a beam, that may or may not be loaded as a column, is assumed to have a torsional spring at each end such that a zero torsional stiffness corresponds to a simply supported end and an infinite torsional stiffness corresponds to a built-in end. From a Rayleigh-Ritz analysis, the buckling load and the fundamental frequency of the beam are each computed as a function of the torsional stiffness. This procedure leads to a one-to-one nondimensional relationship between the buckling load and the natural frequency. From these calculations, it is seen that regardless of the degree of clamping of one end relative to the other end, all that is needed to predict the buckling load within a 15-percent range is a knowledge of the theoretical buckling load of the simply supported column; the theoretical fundamental frequency of the simply supported beam; and the experimental fundamental frequency. Experimental results are presented to support the theory.     

The asymptotic solution of a dynamic buckling problem in elastic columns

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

A procedure of the method of reverberation ray matrix for the buckling analysis of a thin multi-span plate

A procedure of the method of reverberation ray matrix(MRRM) is developed to perform the buckling analysis of thin multi-span rectangular plates having internal line supports or stiffeners. A computation algorithm for the reverberation ray matrix in the MRRM is derived to determine the buckling loading. Specifically, the analytical solutions are presented for the buckling of the structure having two opposite simply-supported or clamped-supported edges with spans, while the constraint condition of two remaining edges may be in any combination of free, simply-supported, and clamped boundary conditions. Furthermore, based on the analysis of matrices relating to the unknown coefficients in the solution form for the deflection in terms of buckling modal functions, some recursive equations(REs) for the MRRM are introduced to generate a reduced reverberation ray matrix with unchanged dimension when the number of spans increases, which promotes the computation efficiency. Several numerical examples are given, and the present results are compared with the known solutions to illustrate the validity and accurateness of the MRRM for the buckling analysis.     

脱层屈曲的高阶摄动解

本文用摄动法求解了三维复合材料中有纤维搭桥的圆形脱层的屈曲问题。搭桥纤维被简化为连续分布于脱层面上的线性弹簧。用Taylor级数展开法确定了每一级摄动方程用Bessel函数表达的解析解．给出了脱层屈曲临界载荷、挠度曲线与脱层半径、搭桥因子等的封闭关系式。根据这一结果并考虑脱层只能向一个方向挠曲的约束条件，可以进一步探讨脱层屈曲的特征。     

Rayleigh–Ritz analysis for localized buckling of a strut on a softening foundation by Hermite functions

This paper proposes a Rayleigh–Ritz procedure for localized buckling of a strut on a non-linear elastic foundation. Firstly, the deflected shape of a strut is expanded into a series of Hermite orthogonal functions, which are proved energy-integrable in an infinite region. Secondly, the errors of the numerical integrations of Hermite functions on the infinite region are investigated and the suitable integral limit is proposed. Through the numerical investigation, it is demonstrated that the first thirty Hermite functions are usually enough to approximate the localized buckling pattern. The proposed method overcomes the disadvantages of the traditional methods, in which the trial functions in either Rayleigh–Ritz or Galerkin analysis are based on the perturbation analyses of the corresponding non-linear differential equation.     

Stochastic homogenization analysis on elastic properties of fiber reinforced composites using the equivalent inclusion method and perturbation method

This paper describes a methodology for evaluation of influence of microscopic uncertainty in material properties and geometry of a microstructure on a homogenized macroscopic elastic property of an inhomogeneous material. For the analysis of the stochastic characteristics of a homogenized elastic property, the first-order perturbation method is used. In order to analyze the influence of microscopic geometrical uncertainty, the perturbation-based equivalent inclusion method is formulated. In this paper, an analytical form of the perturbation term using the equivalent inclusion method is provided.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance of the homogenized elastic tensor of a unidirectional fiber reinforced plastic, which is caused by microscopic uncertainty in material properties or geometry of a microstructure, are estimated with computing the first order perturbation term of the homogenized elastic tensor. Compared the results of the proposed method with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated.     

复合材料层合板的二次屈曲和二次分枝点分析

为了研究复合材料层合板的二次分叉特性 ,利用能量变分原理和非线性几何方程建立了具有弹性约束的复合材料层合板在面内载荷作用下的非线性稳定性控制方程组。控制方程组用广义傅立叶级数法进行求解 ,并得到载荷 -挠度曲线。基于分叉理论中的 Lerray-Schaulder定理 ,采用小挠动法 ,直接导出了复合材料层合板的二次失稳方程。研究结果表明 ,非对称层板也可能存在分叉 ,弹性转动支持系数和铺层等因素对二次分叉有很重要的影响。随着弹性系数的增大 ,二次失稳载荷值与初次失稳载荷值之比下降     

Perturbation formulation of continuation method including limit and bifurcation points

This paper gives the perturbation formulation of continuation method for nonlinear equations. Emphasis is laid on the discussion of searching for the singular points on the equilibrium path and of tracing the paths over the limit or bifurcation points. The method is applied to buckling analysis of thin shells. The pre-and post-buckling equilibrium paths and deflections can be obtained, which are illustrated in examples of buckling analysis of cylindrical and toroidal shells.     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式,运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后,利用非 正交多项式展开式和正交多项式展开式的关系矩阵,可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明,当随机量涨落较大时,递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果,显示了该方法对几何非线性 随机问题求解的有效性.     

薄宽带材局部纵向屈曲的辛弹性力学求解方法

局部纵向屈曲是普遍存在于薄宽带材生产过程的板形缺陷,是屈曲研究的难点,精确的解析求解方法对局部纵向屈曲形成机理的研究和板形质量的提高具有重要意义。本文将任意位置的局部纵向屈曲分为带材边部和内部两类,采用辛弹性力学方法直接推导得到了局部纵向屈曲区域承受不同边界约束条件时的临界屈曲应力和屈曲挠度函数,并将求解结果与有限元和相关文献结果进行了对比。结果表明：辛弹性力学方法与有限元方法相比具有相同计算精度和更高的计算效率,计算精度高于传统能量法;带材边界的约束条件对临界屈曲应力、屈曲区域几何形状和屈曲挠度函数均存在显著影响,验证了传统能量法求解的不足,有利于提高局部屈曲计算精度。     

A boundary layer theory for the buckling of thin cylindrical shells under external pressure

Based on the boundary layer theory for the buckling of thin elastic shells suggested in ref. [14]. the buckling and postbuckling behavior of clamped circular cylindrical shells under lateral or hydrostatic pressure is studied applying singular perturbation method by taking deflection as perturbation parameter. The effects of initial geometric imperfection are also considered. Some numerical results for perfect and imperfect cylindrical shells are given. The analytical results obtained are compared with some experimental data in detail, which shows that both are rather coincident.     

Finite strip buckling analysis of curved plate assemblies under biaxial loading

A finite strip method is presented for calculating the linear buckling stresses of structural assemblies of long, thin plate components which, in general, are curved and which are rigidly joined together at their longitudinal edges. It is assumed that on buckling under the action of a biaxial direct stress field the perturbation forces and displacements all vary sinusoidally in the longitudinal direction. A stiffness matrix relating the amplitudes of the perturbation forces and displacements is developed for the curved strip on the further assumption of relatively high-order polynomial variations of the displacement components around the plate width. Numerical results are presented of the application of the curved strip in calculating the buckling stresses of plates, cylinders, panels and formed sections.     

On dynamic buckling phenomena in axially loaded elastic-plastic cylindrical shells

Some characteristic features of the dynamic inelastic buckling behaviour of cylindrical shells subjected to axial impact loads are discussed. It is shown that the material properties and their approximations in the plastic range influence the initial instability pattern and the final buckling shape of a shell having a given geometry. The phenomena of dynamic plastic buckling (when the entire length of a cylindrical shell wrinkles before the development of large radial displacements) and dynamic progressive buckling (when the folds in a cylindrical shell form sequentially) are analysed from the viewpoint of stress wave propagation resulting from an axial impact. It is shown that a high velocity impact causes an instantaneously applied load, with a maximum value at t=0 and whether or not this load causes an inelastic collapse depends on the magnitude of the initial kinetic energy.     

Thermal buckling and postbuckling behavior of FG-GRC laminated cylindrical shells with temperature-dependent material properties

Thermal postbuckling analysis is presented for graphene-reinforced composite (GRC) laminated cylindrical shells under a uniform temperature field. The GRC layers are arranged in a functionally graded (FG) graphene reinforcement pattern by varying the graphene volume fraction in each GRC layer. The GRCs possess temperature dependent and anisotropic material properties and the extended Halpin–Tsai model is employed to evaluate the GRC material properties. The governing equations are based on a higher order shear deformation shell theory and include the von Kármán-type kinematic nonlinearity and the thermal effects. A singular perturbation method in conjunction with a two-step perturbation approach is applied to determine the thermal postbuckling equilibrium path for a GRC shell with or without geometric imperfection. An iterative scheme is developed to obtain numerical thermal buckling temperatures and thermal postbuckling load–deflection curves for the shells. The results reveal that the FG-X piece-wise FG graphene distribution can enhance the thermal postbuckling capacity of the shells when the shells are subjected to a uniform temperature loading.