Controlled buckling of thin film on elastomeric substrate in large deformation

Electronic systems with large stretchability have many applications. A precisely controlled buckling strategy to increase the stretchability has been demonstrated by combining lithographically patterned surface bonding chemistry and a buckling process. The buckled geometry was assumed to have a sinusoidal form, which may result in errors to determine the strains in the film. A theoretical model is presented in this letter to study the mechanics of this type of thin film/substrate system by discarding the assumption of sinusoidal buckling geometry. It is shown that the previous model overestimates the deflection and curvature in the thin film. The results from the model agree well with finite element simulations and therefore provide design guidelines in many applications ranging from stretchable electronics to micro/nano scale surface patterning and precision metrology.     

BUCKLING ANALYSIS OF STRETCHABLE FERROELECTRIC THIN FILM ON ELASTOMERIC SUBSTRATES

The thin stiff films on pre-stretched compliant substrates can form wrinkles, which can be controlled in micro and nanoscale systems to generate smart structures. Recently, buck- led piezoelectric/ferroelectrie nanoribbons have been reported to show an enhancement in the piezoelectric effect and stretchability, which can be applied in energy harvesting devices, sensors and memory devices instead of polymeric polyvinylidine fluoride （PVDF）. This paper studies the buckling and post-buckling process of ferroelectric thin films bonded to the pre-stretched soft layer, which in turn lies on a rigid support. Nonlinear electromechanical equations for the buckling of thin piezoelectric plates are deduced and employed to model the ferroelectric film poled in the thickness direction. Two buckling modes are analyzed and discussed： partially de-adhered buck- ling and fully adhered buckling. Transition from one buckling mode to the other is predicted and the effect of piezoelectricity on the critical buckling condition of piezoelectric film is examined.     

Buckling-induced dislocation emission in thin films on substrates

Atomistic simulations of the evolution of a strained thin film on a substrate has been reported and the formation of dislocations has been observed in the film/substrate interface after the film has buckled. In the framework of the linear elasticity theory, an analytical model has been developed to explain the buckle effect on the formation of the dislocations. A stability diagram with respect to the buckling and dislocation emission phenomena is finally presented for the film as a function of the uniaxial strain and the Burgers vector.     

水平圆管中受压扭作用管桩屈曲后的解析解

首先建立了考虑自重、受水平圆管约束的管柱在压担组合作用下的屈曲方程：此方程为四阶强非线性常微分方程,难以直接精确求解;通过管柱的正弦屈由分析,找到了方程中的小参数;利用小参数摄动法求解了强非线性屈曲方程,得到了管柱螺旋屈曲后的实际屈曲构形解析解,由此可进一步得到管柱发生螺旋屈曲的临界载荷;所得解析结果与数值计算结果及本文结果退化后与文献中的相关结果比较都有良好的一致性．     

基于辛Runge-Kutta方法的棋盘形褶皱二维薄膜-基底结构动力学特性研究

基于力学屈曲原理的褶皱薄膜-基底结构已成功应用于制备可延展无机电子器件。然而,该类电子器件在应用时需要服役于复杂动态环境中,针对棋盘形褶皱薄膜结构的动力学问题鲜有研究,此问题又是该类电子器件走向实际应用需要解决的关键问题之一。本文首先采用能量方法,分别计算了二维薄膜的弯曲能、膜弹性能和柔性基底中的弹性能以及薄膜动能;然后采用拉格朗日方程,推导出了该结构的振动控制方程;而该方程为非线性动力学方程,无法给出其解析解;因此,本文采用辛Runge-Kutta方法对其进行数值求解;数值结果表明,辛数值方法具有长期稳定的特性和系统结构特性,为高精度的可延展电子器件的动力学问题研究提供了优异的数值方法。     

Tensile and compressive buckling of plates weakened by cracks

Plates are susceptible to buckling under compression when the thickness dimension becomes sufficiently small. Such mode of structure failure can prevail even if the plates were extended in tension. Wribkling of stretched thin sheets is a commonly observed phenomenon that leads to complex deflection patterns, particularly in regions close to crack-like imperfections. Predictions of the buckled displacement modes for plates weakened by cracks will be made on the basis of a theory formulated by application of variational calculus. Finite element method is used such that defects of other shapes can also be analyzed. Various buckled displacement modes of a center-cracked plate are determined and displayed graphically. The critical buckling loads are found to decrease with increasing crack size. Moreover, local wrinkling of the plate surface becomes less pronounced for the higher buckling modes. The method of solution applies equally well to plates that are initially curved.     

Finite width effect of thin-films buckling on compliant substrate: Experimental and theoretical studies

Buckling of stiff thin films on compliant substrates has many important applications ranging from stretchable electronics to precision metrology and sensors. Mechanics plays an indispensable role in the fundamental understanding of such systems. Some existing mechanics models assume plane-strain deformation, which do not agree with experimental observations for narrow thin films. Systematic experimental and analytical studies are presented in this paper for finite-width stiff thin films buckling on compliant substrates. Both experiments and analytical solution show that the buckling amplitude and wavelength increase with the film width. The analytical solution agrees very well with experiments and therefore provides valuable guide to the precise design and control of the buckling profile in many applications. The effect of film spacing is studied via the analytical solutions for two thin films and for periodic thin films.     

Postbuckling behavior of composite-stiffened-curved panels loaded in compression

In this paper, a methodology is developed for the design of a weight-efficient, composite, curved-stiffened panel, loaded in compression well beyond the initial buckling load. A stiffened fuselage panel is designed to satisfy typical design load criteria for the moderately loaded sections of a typical fighter aircraft. Several stiffened panels are fabricated. Some panels are tested to determine experimentally the static strength and the remaining panels are subjected first to severe fatigue loading and then tested statically to determine the effect of fatigue loading on the postbuckling strengh. The experimentally observed behavior is compared with analytical predictions. The weight efficiencies of buckled and unbuckled construction are also compared.     

Energy method solution for the postbuckling response of an axially loaded bilaterally constrained beam

Bistable and multistable structures have shown great usefulness in many applications such as MEMS actuation and energy harvesting. Bistability of structures can be achieved through buckling. Confining a buckled beam between two lateral constraints allows it to buckle into higher modes as the axial load increases. This paper presents a theoretical study of the postbuckling response of a bilaterally constrained elastica subjected to gradually increased axial load. Equilibrium states are determined using an energy method. Under small deformation assumptions, the total potential energy is minimized under the defined constraints. The presented model allows for an accurate representation of the flatting behavior and the increase in the length of contact areas with the lateral constraints before the sudden snapping between equilibrium states. Mode transitions are manifested by jumps in the response curves. Previously developed models based on geometry and symmetries overestimate the required forces for higher equilibrium modes and do not match experimental observations. Results are validated with experimental force–displacement measurements under both force- and displacement-control. The kinetic energy released during buckling mode transitions is determined by a dynamic analysis.     

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.     

Inhomogeneous large deformation study of temperature-sensitive hydrogel

In this paper, inhomogeneous deformation of a temperature-sensitive hydrogel has been studied and analyzed under arbitrary geometric and boundary conditions. We present the governing equations and equilibrium conditions of an isothermal process based on the monophase gel field theory of hydrogel via a variational approach. We have adopted and implemented an explicit form of energy for temperature-sensitive hydrogel in a three-dimensional finite element method (FEM) using a user-supply subroutine in ABAQUS. For verification purpose, a few numerical results obtained by the proposed approach are compared with existing experimental data and analytical solutions. They are all in good agreement. We also provide several examples to show the possible applications of the proposed method to explain various complex phenomena, including the bifurcation, buckling of membrane, buckling of thin film on compliant substrate and the opening and closure of flowers.     

Nonlinear buckled states of rectangular plates

A constructive method is developed to establish the existence of buckled states of a thin, flat elastic plate that is rectangular in shape, simply supported along its edges, and subjected to a constant compressive thrust applied normal to its two short edges. Under the assumption that the stress function and the deformation of the plate are described by the nonlinear von Kármán equations, the approach used yields information regarding not only the number of buckled states near an eigenvalue of the linearized problem, but also the continuous dependence of such states on the load parameter and the possible selection of that buckled state “preferred” by the plate. In particular, the methods used provide a rigorous approach to studying the existence of buckled states near the first eigenvalue of the linearized problem (that is, near the “buckling load”) even when the first eigenvalue is not simple.     

Exact solution and stability of postbuckling configurations of beams

We present an exact solution for the postbuckling configurations of beams with fixed–fixed, fixed–hinged, and hinged–hinged boundary conditions. We take into account the geometric nonlinearity arising from midplane stretching, and as a result, the governing equation exhibits a cubic nonlinearity. We solve the nonlinear buckling problem and obtain a closed-form solution for the postbuckling configurations in terms of the applied axial load. The critical buckling loads and their associated mode shapes, which are the only outcome of solving the linear buckling problem, are obtained as a byproduct. We investigate the dynamic stability of the obtained postbuckling configurations and find out that the first buckled shape is a stable equilibrium position for all boundary conditions. However, we find out that buckled configurations beyond the first buckling mode are unstable equilibrium positions. We present the natural frequencies of the lowest vibration modes around each of the first three buckled configurations. The results show that many internal resonances might be activated among the vibration modes around the same as well as different buckled configurations. We present preliminary results of the dynamic response of a fixed–fixed beam in the case of a one-to-one internal resonance between the first vibration mode around the first buckled configuration and the first vibration mode around the second buckled configuration.     

Buckling of delaminated bi-layer beam-columns

An improved analytical model is presented to analyze the delamination buckling of a bi-layer beam-column with a through-the-width delamination. Both the transverse shear deformation and local delamination tip deformations are taken into consideration, and two delaminated sub-layers as well as two substrates in the intact (un-delaminated) regions are modeled as individual Timoshenko beams. A deformable interface is introduced to establish the continuity condition between the two substrates in the intact regions. Consequently, a flexible joint is formed at the delamination tip, and it is different from the conventional rigid joint given in most of studies in the literature, in which the local delamination tip deformations are completely ignored. In contrast to the local delamination buckling in our previous study (Qiao et al., 2010), the present model accounts for the global deformations of the intact region in the delaminated composite beam-column, thus capable of capturing the buckling mode shape transitions from the global, to global–local coexistent, and to local buckling for asymmetric delamination as the interface delamination increases. Good agreement of the present analytical solutions with the full 2-D elastic finite element analysis demonstrates the local deformation effects around the delamination tip and verifies the accuracy of the present model. Parametric studies are conducted to investigate the effects of loading eccentricity, delaminated sub-layer thickness ratio, and interface compliance on the critical buckling load for the delaminated composite beam-column. Transitions of buckling modes from the global to local delamination buckling are also disclosed as the thickness of one sub-layer reduces from the thick sub-layer to a thin film. The developed delamination buckling solution facilitates the design analysis and optimization of laminated composite structures, and it can be used with confidence in buckling analysis of delaminated composite structures.     

Folding of Thin Walled Tubes as a Free Gradient Discontinuity Problem

The aim of this paper is to describe buckling deformations of hollow cylinders whose buckled configurations consist of inextensional deformations of their original middle surface, through minimization of two competing energy forms: a bulk elastic energy and an interface bending energy. These minimal energy configurations are obtained through descent minimization within the class of folding deformations. The non-standard problem of minimization over variable triangulations is considered.     

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.     

Static instability analysis for travelling membranes and plates interacting with axially moving ideal fluid

The out-of-plane instability of a moving plate, travelling between two rollers with constant velocity, is studied, taking into account the mutual interaction between the buckled plate and the surrounding, axially flowing ideal fluid. Transverse displacement of the buckled plate (assumed cylindrical) is described by an integro-differential equation that includes the centrifugal force, the aerodynamic reaction of the external medium, the vertical projection of membrane tension, and the bending force. The aerodynamic reaction is found analytically as a functional of the displacement. To find the critical divergence velocity of the moving plate and its corresponding buckling mode, an eigenvalue problem and variational principle are derived. Plate divergence, both within a vacuum and when submerged in an external medium, is investigated with the application of analytical and numerical techniques.     

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.     

Buckling and Post-buckling Behaviour of Prismatic Aluminium Columns Submitted to a Series of Compressive Loads

The buckling and post-buckling behaviour of prismatic aluminium columns from stocky to very slender shapes is investigated. The unconventional, in terms of buckling tests, displacement control of compressive load and a series of loadings provided an enhanced insight into the buckling process. A phenomenon of buckling load drop has been detected in columns of intermediate slenderness, reaching over 20% of the load early critical value. This newly observed occurrence resembles finite disturbance instability, which until recently was commonly believed to only appear in cases of thin walled cylindrical shells, but not columns. The observation is in contradiction to predicted results from the elasto-plastic buckling models of Engesser or Shanley, with constant or growing values of load during the post-buckling process. Further tests on columns of intermediate slenderness, with strain gauges glued at node and anti-node locations of the buckled profiles, revealed that even minute buckling results in fields of highly non-symmetric residual microplastic strain. The results of the present study indicate that running column buckling tests under displacement control is worthy of being adopted as common practice. The envelope of column post-buckling states can be conveniently determined. This information will in turn allow for the quick and reliable estimation of the safety of a column, which has undergone accidental or deliberate damage in the form of limited buckling when under operational load.     

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

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