Bending and stability analysis of gradient elastic beams

The problems of bending and stability of Bernoulli–Euler beams are solved analytically on the basis of a simple linear theory of gradient elasticity with surface energy. The governing equations of equilibrium are obtained by both a combination of the basic equations and a variational statement. The additional boundary conditions are obtained by both variational and weighted residual approaches. Two boundary value problems (one for bending and one for stability) are solved and the gradient elasticity effect on the beam bending response and its critical (buckling) load is assessed for both cases. It is found that beam deflections decrease and buckling load increases for increasing values of the gradient coefficient, while the surface energy effect is small and insignificant for bending and buckling, respectively.     

Size-dependent axial buckling analysis of functionally graded circular cylindrical microshells based on the modified strain gradient elasticity theory

In this paper, a size-dependent first-order shear deformable shell model is developed based upon the modified strain gradient theory (MSGT) for the axial buckling analysis of functionally graded (FG) circular cylindrical microshells. It is assumed that the material properties of FG materials, which obey a simple power-law distribution, vary through the thickness direction. The principle of virtual work is utilized to formulate the governing equations and corresponding boundary conditions. Numerical results are presented for the axial buckling of FG circular cylindrical microshells subject to simply-supported end conditions and the effects of material length scale parameter, material property gradient index, length-to-radius ratio and circumferential mode number on the size-dependent critical buckling load are extensively studied. For comparison purpose, the critical buckling loads predicted by modified couple stress theory (MCST) and classical theory (CT) are also presented. Results show that the size effect plays an important role for lower values of dimensionless length scale parameter. Moreover, it is observed that the critical buckling loads obtained based on MSGT are greater than those obtained based on MCST and CT.     

CLASSICAL AND HOMOGENIZED EXPRESSIONS FOR BUCKLING SOLUTIONS OF FUNCTIONALLY GRADED MATERIAL LEVINSON BEAMS

The relationship between the critical buckling loads of functionally graded material(FGM) Levinson beams(LBs) and those of the corresponding homogeneous Euler-Bernoulli beams(HEBBs) is investigated. Properties of the beam are assumed to vary continuously in the depth direction. The governing equations of the FGM beam are derived based on the Levinson beam theory, in which a quadratic variation of the transverse shear strain through the depth is included.By eliminating the axial displacement as well as the rotational angle in the governing equations,an ordinary differential equation in terms of the deflection of the FGM LBs is derived, the form of which is the same as that of HEBBs except for the definition of the load parameter. By solving the eigenvalue problem of ordinary differential equations under different boundary conditions clamped(C), simply-supported(S), roller(R) and free(F) edges combined, a uniform analytical formulation of buckling loads of FGM LBs with S-S, C-C, C-F, C-R and S-R edges is presented for those of HEBBs with the same boundary conditions. For the C-S beam the above-mentioned equation does not hold. Instead, a transcendental equation is derived to find the critical buckling load for the FGM LB which is similar to that for HEBB with the same ends. The significance of this work lies in that the solution of the critical buckling load of a FGM LB can be reduced to that of the HEBB and calculation of three constants whose values only depend upon the throughthe-depth gradient of the material properties and the geometry of the beam. So, a homogeneous and classical expression for the buckling solution of FGM LBs is accomplished.     

Analytical solution for buckling of asymmetrically delaminated Reissner’s elastic columns including transverse shear

The exact analytical solution of buckling in delaminated columns is presented. In order to investigate analytically the influence of axial and shear strains on buckling loads the geometrically exact beam theory is employed with no simplification of the governing equations. The critical forces are then obtained by the linearized stability theory. In the paper, we limit the studies to linear elastic columns with a single delamination, but with arbitrary longitudinal and vertical asymmetry of delamination and arbitrary boundary conditions. The studies of quantitative and qualitative influence of transverse shear are shown in detail and extensive results for buckling loads with respect to delamination length, thickness and longitudinal position are presented.     

复合材料脱层梁的屈曲分析

本文研究了具有任意位置透型脱层的复合材料梁的屈曲问题。基于弹性理论建立了复合材料脱层梁的基本方程式。对脱层梁进行了分区处理,利用B样条函数作为位移型函数的基函数,方便地描述了脱层长度、脱层位置。考虑边界条件、区间位移连续性条件和弯矩剪力的平衡条件以及纵向内力的附加条件,对基本方程式进行了求解。得出了脱层位置不同,脱层长度不同的屈曲荷载的变化规律,并与轴对称脱层时的屈曲荷载进行了比较,认为层合梁考虑脱层对屈曲的影响是非常必要的。     

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

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

Elastica of sandwich panels with a transversely flexible core—A high-order theory approach

The elastica behavior of an extensional sandwich panel with a “soft” core when subjected to in-plane compressive loads is presented and it is compared with the response of its extensional equivalent single layer (ESL) with shear deformations model. The field equations along with the appropriate boundary conditions for the sandwich and the ESL panels have been derived through a variational approach following the High-order SAndwich Panel Theory (HSAPT) approach that takes into account the vertical flexibility of the core. The governing equations include the effects of the extension of the mid-surfaces of the face sheets of the sandwich panel or the mid-plane of the ESL model which the classical elastica approach misses. The results of the elastica response of a clamped-simply-supported sandwich panel and its ESL counterpart are presented and compared. They include the response along the panel, deformed shapes and equilibrium curves of in-plane loads versus structural quantities such as displacements and internal stress resultants and stresses. These results reveal that the predicted buckling load of the ESL panel is larger than that of the sandwich panel and that deep in the non-linear range the upper face sheet wrinkles with increasing overall and edge displacements and a release of the load. Hence, the use of an equivalent single layer panel especially when a sandwich panel with a compliant core is considered may lead to unsafe and unreliable predictions when large displacements and large rotations are considered.     

Initial thermo-mechanical post-buckling of beams with temperature-dependent physical properties

The objective of this work is to analyze the elastic buckling and initial post-buckling behavior of slender beams subjected to uniform heating. The beams are assumed to be double-hinged with fixed ends, preventing thermal expansion. Consequently, destabilizing compressive forces arise that may lead to beam buckling. When the temperature is further increased, the beam experiences finite displacements, with the result that the analysis is geometrically non-linear. The modulus of elasticity and the thermal induced strain, key material properties for this problem, are temperature-dependent. Thus, the coefficients of the governing equations are not constant. This suggests the physical non-linearity of the mathematical model. Hence, the analysis is geometrically and physically non-linear. The analysis is sensitive to the beam initial temperature, as the thermal strain is a function of the initial and final temperatures. The material is considered to be linear elastic, and consequently viscoelastic and plastic effects are not taken into account. Furthermore, the beam cross-section properties are assumed to be constant, which is consistent with the small strain formulation. A perturbation method is applied to the governing non-linear differential equations so that the initial post-buckling behavior may be analytically determined when temperatures above the critical temperature are applied to the beam. To illustrate the application of the formulation we present a case study for the aluminum 7075-T6 alloy, a material commonly used in aerospace and naval industries. Nonetheless, it is expected similar behavior for other metallic materials. The curves that define the variation of the modulus of elasticity, the thermal strain and the yield stress with temperature are considered in our analysis. The change in length, reaction forces at the supports and geometric configurations are obtained as a function of temperature and the beam slenderness ratio. The critical buckling loads and temperatures and the initial post-buckling analysis are also calculated in the context of the temperature-independent physical properties. Our results emphasize the importance of modeling the material's non-linearity if accuracy is required. However, from a practical application point of view results are acceptable if temperature-independent physical properties are employed, especially for large slenderness ratios.     

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

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

Dynamic buckling mode of an elastic bar

The conception of buckling relative initial imperfection is presented in this paper. According to Boulli-Euler beam equation, the dynamic buckling mode of an elastic bar under the homogeneous boundary conditions can be derived by applying the preferred mode analytical method. As an example, the dynamic buckling mode of an elastic bar clamped at both ends is discussed.     

Analytical solutions for buckling of size-dependent Timoshenko beams

The inconsistences of the higher-order shear resultant expressed in terms of displacement(s) and the complete boundary value problems of structures modeled by the nonlocal strain gradient theory have not been well addressed. This paper develops a size-dependent Timoshenko beam model that considers both the nonlocal effect and strain gradient effect. The variationally consistent boundary conditions corresponding to the equations of motion of Timoshenko beams are reformulated with the aid of the weighted residual method. The complete boundary value problems of nonlocal strain gradient Timoshenko beams undergoing buckling are solved in closed forms. All the possible higher-order boundary conditions induced by the strain gradient are selectively suggested based on the fact that the buckling loads increase with the increasing aspect ratios of beams from the conventional mechanics point of view. Then, motivated by the expression for beams with simply-supported(SS) boundary conditions, some semiempirical formulae are obtained by curve fitting procedures.     

Instability of functionally graded micro-beams via micro-structure-dependent beam theory

This paper focuses on the buckling behaviors of a micro-scaled bi-directional functionally graded (FG) beam with a rectangular cross-section, which is now widely used in fabricating components of micro-nano-electro-mechanical systems (MEMS/NEMS) with a wide range of aspect ratios. Based on the modified couple stress theory and the principle of minimum potential energy, the governing equations and boundary conditions for a micro-structure-dependent beam theory are derived. The present beam theory incorporates different kinds of higher-order shear assumptions as well as the two familiar beam theories, namely, the Euler-Bernoulli and Timoshenko beam theories. A numerical solution procedure, based on a generalized differential quadrature method (GDQM), is used to calculate the results of the bi-directional FG beams. The effects of the two exponential FG indexes, the higher-order shear deformations, the length scale parameter, the geometric dimensions, and the different boundary conditions on the critical buckling loads are studied in detail, by assuming that Young’s modulus obeys an exponential distribution function in both length and thickness directions. To reach the desired critical buckling load, the appropriate exponential FG indexes and geometric shape of micro-beams can be designed according to the proposed theory.     

铁磁梁式板磁弹性初始后屈曲及缺陷敏感性分析

从铁磁板的磁弹性广义变分原理出发，通过对铁磁板内外的磁场和板变形 场的摄动技术以及采用经典Koiter理论，对横向磁场中铁磁悬臂、简支、固支梁式薄 板磁弹性稳定性的初始后屈曲行为及缺陷敏感性进行了定性研究. 解析地给出了悬臂板对初 始缺陷敏感，简支和固支板对初始缺陷不敏感等结果，对悬臂梁式板理论预测的临界磁场值 往往大于实验观测结果的现象从定性上给予了合理解释.     

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress,the modified couple stress theory(MCST),and the nonlocal elasticity theories using the differential quadrature method(DQM)is presented.Main advantages of the MCST over the classical theory(CT)are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter.Based on the nonlinear von K′arm′an assumption,the governing equations of equilibrium for the micro-classical plate considering midplane displacements are derived based on the minimum principle of potential energy.Using the DQM,the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained.Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature.A parametric study is conducted to show the effects of the aspect ratio,the side-to-thickness ratio,Eringen’s nonlocal parameter,the material length scale parameter,Young’s modulus of the surface layer,the surface residual stress,the polymer matrix coefficients,and various boundary conditions on the dimensionless uniaxial,biaxial,and shear critical buckling loads.The results indicate that the critical buckling loads are strongly sensitive to Eringen’s nonlocal parameter,the material length scale parameter,and the surface residual stress effects,while the effect of Young’s modulus of the surface layer on the critical buckling load is negligible.Also,considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate.The results show that the critical biaxial buckling load increases with an increase in G12/E2and vice versa for E1/E2.It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude.Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios,it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.     

Thermal buckling analysis of functionally graded beam with longitudinal crack

The thermal buckling problem of functionally graded beam with longitudinal crack is presented in the paper. The whole beam is divided into four sub-beams and each one is modeled as a Timoshenko beam. The buckling governing equation of each sub-beam in thermal environment is established by using Hamilton Principle. Combining with the boundary conditions, the continuous conditions of the displacements and the forces, the buckling governing equations are solved by both the analytical and numerical methods. The buckling modes and critical buckling temperatures are obtained, and the effects of the functionally graded index, crack length, crack depth, and crack longitudinal location on the buckling characteristics of beams are discussed in numerical examples.     

Post-buckling analysis of slender elastic vertical rods subjected to terminal forces and self-weight

This article presents the behavior of slender elastic rods subjected to axial terminal forces and self-weight. The mathematical formulation is presented, a solution is sought for a double-hinged boundary condition and the analysis is carried out for different values of non-dimensional weight. The formulation derives from geometrical compatibility, equilibrium of forces and moments and constitutive relations yielding a set of six first order non-linear ordinary differential equations with boundary conditions specified at both ends, which characterizes a complex two-point boundary value problem. Furthermore, a perturbation method is used to find the critical buckling loads and initial post-buckling solutions. A numerical integration scheme based on a three parameter shooting method is employed in the post-buckling solutions.     

Continuum-based stability of anisotropic and auxetic beams

The critical buckling loads of pinned-pinned and cantilever beams are computed using the equations of three-dimensional elasticity rather than typical beam theories. These loads are influenced both by the nature of the assumed displacement field over the beam cross-section and by the inclusion of the terms from the full constitutive tensor. Of special interest are beams that are either anisotropic or auxetic. For anisotropic beams, an increased ratio of longitudinal to shear modulus for cantilevered beams increases the generation of shear buckling rather than flexural buckling. For isotropic auxetic beams, the values of Poisson ratio that define the limit between buckling loads that approach the classical buckling load from above or below are discussed.     

Two-dimensional elasticity solution for bending of functionally graded beams with variable thickness

This paper studies the stress and displacement distributions of functionally graded beam with continuously varying thickness, which is simply supported at two ends. The Young’s modulus is graded through the thickness following the exponential-law and the Poisson’s ratio keeps constant. On the basis of two-dimensional elasticity theory, the general expressions for the displacements and stresses of the beam under static loads, which exactly satisfy the governing differential equations and the simply supported boundary conditions at two ends, are analytically derived out. The unknown coefficients in the solutions are approximately determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beams. The effect of Young’s modulus varying rules on the displacements and stresses of functionally graded beams is investigated in detail. The two-dimensional elasticity solution obtained can be used to assess the validity of various approximate solutions and numerical methods for the aforementioned functionally graded beams.     

Investigation of the Nonlinear Behavior of Beams*

ABSTRACT A system of equations that describes the nonlinear behavior of teams is presented. The differential equations are solved by using the Galerkin method, and a system of nonlinear algebric equations is obtained. A method to deal with discontinuities in structural properties and load distribution along a beam is presented. The derivation includes general expressions for gravity-type loads. Two methods of calculating the resultant moment along the beam are described. One method incorporates differentiation of the expressions for displacements, while the second method is based on integration of loads along the beam. The numerical model is used in order to investigate displacements and moment resultants along a cantilevered beam. Different problems that are associated with nonlinear behavior are presented and discussed. Good agreement is obtained between theoretical and experimental results.     

脱层梁屈曲的高阶剪切理论

脱层的存在将会大大降低层合结构的屈曲载荷。该文将含任意位置脱层的两端固支梁分成多段子层,用厚度的三次多项式模拟脱层梁屈曲时子层的轴向位移,利用变分原理和欧拉方程导出了脱层梁的屈曲方程和定解条件,并用状态空间方法进行求解。通过与一阶剪切理论和经典理论的比较,指出了它们各自的适用范围;考虑了脱层梁三种不同的屈曲模态。分析了脱层长度、深度、位置和材料的铺层方向对脱层梁屈曲载荷的影响;最后给出了多处简单脱层的屈曲分析。