带旋转自由度C^0类任意四边形板（壳）单元

基于Ｒｅｉｓｓｎｅｒ－Ｍｉｎｄｉｌｉｎ板弯曲理论和Ｖｏｎ－Ｋａｒｍａｎ大挠度理论,采用单元域内和边界位移插值一致性的概念,将四节点等参弯曲单元与Ａｌｌｍａｎ膜变形二次插值模式相结合,对层合板壳的大挠度分析提供了一种实用的带旋转自由度的四节点Ｃ＾０类板单元。大量算例表明：该单元对板壳结构的线性强度、稳定性和后屈曲分析都表现出良好的收敛性和足够的工程精度。     

非线性弹性基础上矩形板热后屈曲分析

给出非线性弹性基础上矩形板在均匀和非均匀（抛物型）热分布作用下的后屈曲分析。采用摄动——Ｇａｌｅｒｋｉｎ混合法给出完善和非完善矩形板热屈曲载荷和热后屈曲平衡路径。数值计算结果表明，非线性弹性基础上矩形板具有不稳定的热后屈曲平衡路径，且对初始几何缺陷是敏感的     

直接基于单元平衡的梁非线性分析方法

以单元力平衡条件为基础建立的梁单元在近年来受到广泛关注,与传统的以单元位移相容条件为基础建立的梁单元相比,基于平衡的梁单元可获得较高的精度。针对基于平衡的梁单元,给出了一种几何非线性分析方法,使用该方法进行几何非线性分析时不需进行单元内部位移场的构造,避免了使用现有位移场构造方法时可能引起的振荡问题,同时降低了有限元程...     

Isogeometric nonlocal strain gradient quasi-three-dimensional plate model for thermal postbuckling of porous functionally graded microplates with central cutout with different shapes

This study presents the size-dependent nonlinear thermal postbuckling characteristics of a porous functionally graded material(PFGM) microplate with a central cutout with various shapes using isogeometric numerical technique incorporating nonuniform rational B-splines. To construct the proposed non-classical plate model, the nonlocal strain gradient continuum elasticity is adopted on the basis of a hybrid quasithree-dimensional(3D) plate theory under through-thickness deformation conditions by only four variables. By taking a refined power-law function into account in conjunction with the Touloukian scheme, the temperature-porosity-dependent material properties are extracted. With the aid of the assembled isogeometric-based finite element formulations,nonlocal strain gradient thermal postbuckling curves are acquired for various boundary conditions as well as geometrical and material parameters. It is portrayed that for both size dependency types, by going deeper in the thermal postbuckling domain, gaps among equilibrium curves associated with various small scale parameter values get lower, which indicates that the pronounce of size effects reduces by going deeper in the thermal postbuckling regime. Moreover, we observe that the central cutout effect on the temperature rise associated with the thermal postbuckling behavior in the presence of the effect of strain gradient size and absence of nonlocality is stronger compared with the case including nonlocality in absence of the strain gradient small scale effect.     

Coupled buckling and postbuckling analysis of active laminated piezoelectric composite plates

A theoretical framework for analyzing the pre- and postbuckling response of composite laminates and plates with piezoactuators and sensors is presented. The mechanics include nonlinear effects due to large rotations and stress stiffening, and are incorporated into a coupled mixed-field piezoelectric laminate theory. Using the previous mechanics, a nonlinear finite element method and an incremental-iterative solution are formulated for the analysis of nonlinear adaptive plate structures subject to in-plane electromechanical loading. A novel eight-node nonlinear plate finite element is also developed. Evaluation cases predict the buckling and postbuckling response of adaptive composite beams and plates with piezoelectric actuators and sensors. The case of piezoelectric buckling and postbuckling induced by the actuators is addressed and quantified. Finally, the possibility to actively mitigate the mechanical buckling and postbuckling response of adaptive piezocomposite plates is illustrated.     

Axisymmetric static and dynamics snap-through phenomena in a thermally postbuckled temperature-dependent FGM circular plate

Thermal postbuckling analysis and the axisymmetric static and dynamic snap-through phenomena due to static/sudden uniform lateral pressure in a thermally postbuckled functionally graded material circular plate are performed in this research. Plate is formulated using the first order shear deformation plate theory. Thermo-mechanical properties of the plate are assumed to be temperature dependent where dependency is described according to the higher order Touloukian representation. Two types of temperature loading are considered. Uniform temperature rise and heat conduction across the thickness direction. The one dimensional heat conduction equation in the thickness direction is obtained and discreted via the central finite difference method. The obtained system of equations is nonlinear since the thermal conductivity itself is a function of the unknown nodal temperatures. Using the von-Kármán assumptions, the governing equations of the plate are obtained in a matrix representation with the aid of the conventional Ritz method whose shape functions are developed using the Gram-Schmidt process. At first thermal postbuckling analysis is performed which is a nonlinear problem with respect to both temperature and displacements. Afterwards, response of the bulged thermally postbuckled plate is obtained under the static and dynamic uniform pressure. Snap-through phenomenon may be observed in both static and dynamic loading cases, due to the immovability of the edge of the plate and the initial deflection caused by postbuckling deflection. To capture the snapping phenomenon and trace the path beyond the limit loads, cylindrical arch-length technique is used. In dynamic snap-through analysis, the effect of structural damping is also included. Numerical results of this study reveal that the structure is sensitive to the initial deflection caused by thermal postbuckling load. Increasing the temperature prior to mechanical loads enhances the snap-through intensity and also increases both the upper and lower limit loads. As shown, dynamic snap-through loads are lower than the static ones, however dynamic snap-through intensity is more than the static snap-though intensity. Furthermore, structural damping enhances the dynamic buckling loads of the plate and decreases the dynamic postbuckling deflection of the plate.     

FRP strengthened walls with delaminated regions – Geometrically nonlinear response

The paper studies the geometrically nonlinear behavior of walls that are strengthened with fiber reinforced polymer (FRP) composite materials but include pre-existing delaminated regions. The paper uses an analytical–numerical methodology. Three specially tailored finite elements that correspond to perfectly bonded regions, to delaminated regions where the debonded layers are in contact, and to delaminated regions where the debonded layers are not in contact are presented. All finite elements are based on a high order multi layered plate theory. The geometrical nonlinearity is introduced by means of the Von Karman nonlinear strains whereas the contact nonlinearity is handled iteratively. The validity and convergence of the finite element models is demonstrated for each type of element through comparison with closed form analytical solutions available for specific cases. The unified model that combines the three types of finite element is then used for studying the nonlinear behavior of a locally delaminated FRP strengthened wall under in-plane normal and in-plane shear loads. Finally, conclusions regarding the effect of the delamination on the response of the strengthening system, on the conditions that evolve in the bonded region that surrounds the delamination, and on the global response of the multi-layered structure are drawn. Additional conclusions regarding the application of the modeling approach to other delamination sensitive layered structural systems close the paper.     

Thermal buckling and postbuckling of laminated composite beams with temperature-dependent properties

The thermal buckling and postbuckling analysis of laminated composite beams with temperature-dependent material properties is presented. The governing equations are based on the first-order shear deformation beam theory (FSDT) and the geometrical nonlinearity is modeled using Green's strain tensor in conjunction with the von Karman assumptions. The differential quadrature method (DQM) as an accurate, simple and computationally efficient numerical tool is adopted to discretize the governing equations and the related boundary conditions. A direct iterative method is employed to obtain the critical temperature (bifurcation point) as well as the nonlinear equilibrium path (the postbuckling behavior) of symmetrically laminated beams. The applicability, rapid rate of convergence and high accuracy of the method are established via different examples and by comparing the results with those of existing in literature. Then, the effects of temperature dependence of the material properties, boundary conditions, length-to-thickness ratios, number of layers and ply angle on the thermal buckling and postbuckling characteristic of symmetrically laminated beams are investigated.     

Nonlinear deformation and buckling of elastic inhomogeneous shells under thermomechanical loads

The paper outlines the fundamentals of the method of solving static problems of geometrically nonlinear deformation, buckling, and postbuckling behavior of thin thermoelastic inhomogeneous shells with complex-shaped midsurface, geometrical features throughout the thickness, or multilayer structure under complex thermomechanical loading. The method is based on the geometrically nonlinear equations of three-dimensional thermoelasticity and the moment finite-element scheme. The method is justified numerically. Results of practical importance are obtained in analyzing poorely studied classes of inhomogeneous shells. These results provide an insight into the nonlinear deformation and buckling of shells under various combinations of thermomechanical loads     

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

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

ANALYSIS OF ELASTO-PLASTIC POSTBUCKLING AND ENERGY RELEASE RATE FOR DELAMINATED FIBER METAL LAMINATED BEAMS IN THERMAL ENVIRONMENT

The elasto-plastic postbuckling of fiber metal laminated beams with delamination and the energy release rate along the delamination front are discussed in this paper. Considering geometrical nonlinearity, thermal environment and geometrical initial imperfection, the incremental nonlinear equilibrium equations of delaminated fiber metal laminated beams are established,which are solved using the differential quadrature method and iterative method. Based on these,according to the J-integral theory, the elasto-plastic energy release rate is studied. The effects of some important parameters on the elasto-plastic postbuckling behavior and energy release rate of the aramid reinforced aluminum laminated beams are discussed in details.     

任意外型薄板非线性后屈曲问题的边界元解

本文就薄板后屈曲问题建立一组新型的边界元计算公式，用这组公式求解能方便处理各种边界问题，另外文中将面内应力分解成基本部份和附加部份，并利用微分算子分解理论导得了挠度的一个不同形式的基本解，由于计算公式中，实现了面内位移和挠度的解耦，从而使迭代过程得到简化，文末还对圆板后屈曲路径进行了计算，得到了满意的结果。     

钢管混凝土拱稳定分析的有限元法

在 2 0结点三维块体等参元及 1 6结点相对位移板壳元 [1 ] 的基础上、引入梁的基本假定 ,采用等效数值积分法 ,构造出 1 2 - 2 0结点三维退化平面和柱面层合曲梁单元 ,同时给出用于梁或拱线弹性稳定性分析的有限元列式 ,最后 ,以绍兴轻纺大桥为工程背景 ,计算出轻纺大桥钢管混凝土拱面内屈曲及面外侧倾屈曲的临界荷载。     

Nonlinear stability of functionally graded material (FGM) sandwich cylindrical shells reinforced by FGM stiffeners in thermal environment

In this paper, Donnell's shell theory and smeared stiffeners technique are improved to analyze the postbuckling and buckling behaviors of circular cylindrical shells of stiffened thin functionally graded material(FGM) sandwich under an axial loading on elastic foundations, and the shells are considered in a thermal environment. The shells are stiffened by FGM rings and stringers. A general sigmoid law and a general power law are proposed. Thermal elements of the shells and reinforcement stiffeners are considered. Explicit expressions to find critical loads and postbuckling load-deflection curves are obtained by applying the Galerkin method and choosing the three-term approximate solution of deflection. Numerical results show various effects of temperature, elastic foundation, stiffeners, material and geometrical properties, and the ratio between face sheet thickness and total thickness on the nonlinear behavior of shells.     

Buckling and postbuckling of magnetoelastic flat plates carrying an electric current

The basic equations of a fully nonlinear theory of electromagnetically conducting flat plates carrying an electric current and exposed to a magnetic field of arbitrary orientation are derived. The relevant equations have been obtained by considering that both the elastic and electromagnetic media are homogeneous and isotropic. The geometrical nonlinearities are considered in the von-Kármán sense, and the soft ferromagnetic material of the plate is assumed to feature negligible hysteretic losses. Based on the electromagnetic and elastokinetic field equations, by using the standard averaging methods, the 3-D coupled problem is reduced to an equivalent 2-D one, appropriate to the theory of plates. Having in view that the elastic structures carrying an electric current are prone to buckling, by using the presently developed theory, the associated problems of buckling and postbuckling are investigated. In this context, the problem of the electrical current inducing the buckling instability of the plate, and its influence on the postbuckling behavior are analyzed. In the same context, the problem of the natural frequency–electrical current interaction of flat plates, as influenced by a magnetic field is also addressed.     

Thermomechanical Postbuckling of Imperfect Moderately Thick Plates on Nonlinear Elastic Foundations*

ABSTRACT

A postbuckling analysis is presented for a moderately thick rectangular plate subjected to combined axial compression and uniform temperature loading, and resting on a softening nonlinear elastic foundation. The cases of (1) thermal postbuckling of initially compressed plates and (2) compressive post-buckling of initially heated plates are considered. The initial geometrical imperfections of the plates are taken into account. Formulations are based on Reissner-Mindlin plate theory, considering first-order shear deformation effects, and including plate-foundation interaction and thermal effects. The analysis uses a deflection-type perturbation technique to determine buckling loads and postbuckling equilibrium paths. Numerical examples include the performance of perfect and imperfect, moderately thick plates resting on softening nonlinear elastic foundations. Typical results are presented in dimensionless graphical form.     

Patterned bilayer plate microstructures subjected to thermal loading: Deformation and stresses

We studied the deformation of a series of gold/polysilicon patterned plate microstructures fabricated by surface micromachining. The patterned plate microstructures were subjected to a uniform temperature change from 100 °C to room temperature that was intended to induce linear and geometrically nonlinear deformation. We used interferometry to measure full-field deformed shapes of the microstructures. From these measurements we determined the spatially-averaged curvature of the deformed microstructures within individual lines and across the entire plate. The deformation response of the patterned plates can be broadly characterized in terms of the average curvature as a function of temperature change and exhibits linear and geometrically nonlinear behavior. We modeled the deformation response of the patterned plates using geometrically nonlinear plate theory with the finite element method. Good agreement was obtained between predictions and measurements for both local curvature variations across lines and for the evolution of curvature of the entire plate with temperature change. Using a generalized plane strain approach with the finite element method we also modeled the spatial dependence of the stress distribution in the lines and substrate. For thick plates, our results agree with those of previous studies, showing a decrease in the von Mises stress in the metal lines with decreasing linewidth. For thinner substrates, though, we find the behavior with linewidth is opposite and there is a critical substrate thickness (about 10 μm for the system in our study) where the behavior with linewidth changes. These results have important implications in the design of patterned structures for micro-electro-mechanical systems (MEMS) applications where films are of comparable thickness to the underlying substrate.     

Effective characteristics of corrugated plates

Corrugated plates are widely used in modern constructions and structures, because they, in contrast to plane plates, possess greater rigidity. In many cases, such a plate can be modeled by a homogeneous anisotropic plate with certain effective flexural and tensional rigidities. Depending on the geometry of corrugations and their location, the equivalent homogeneous plate can also have rigidities of mutual influence. These rigidities allow one to take into account the influence of bending moments on the strain in the midplane and, conversely, the influence of longitudinal strains on the plate bending [1]. The behavior of the corrugated plate under the action of a load normal to the midsurface is described by equations of the theory of flexible plates with initial deflection. These equations form a coupled system of nonlinear partial differential equations with variable coefficients [2]. The dependence of the coefficients on the coordinates is determined by the corrugation geometry. In the case of a plate with periodic corrugation, the coefficients significantly vary within one typical element and depend on the values of local variables determined in each of the typical elements. There is a connection between the local and global variables, and therefore, the functions of local coordinates are simultaneously functions of global coordinates, which are sometimes called rapidly oscillating functions [3].One of the methods for solving the equations with rapidly oscillating coefficients is the asymptotic method of small geometric parameter. The standard procedure of this method usually includes preparatory stages. At the first stage, as a rule, a rectangular periodicity cell is distinguished to be a typical element. At the second stage, the scale of global coordinates is changed so that the rectangular structure periodicity cells became square cells of size l × l. The third stage consists in passing to the dimensionless global coordinates relative to the plate characteristic dimension L. As a result, the dependence between the new local variables and the new scaled dimensionless variables is such that the factor 1/α, where α=l/L ? 1 is a small geometric parameter, appears in differentiating any function of the local coordinate with respect to the global coordinate. After this, the solution of the problem in new coordinates is sought as an asymptotic expansion in the small geometric parameter [1], [4–10].We note that, in the small geometric parameter method, the asymptotic series simultaneously have the form of expansions in the gradients of functions depending only on the global coordinates. This averaging procedure can be applied to linear and nonlinear boundary value problems for differential equations with variable periodic coefficients for which the periodicity cell can be affinely transformed into the periodicity cube. In the case of an arbitrary dependence of the coefficients on the coordinates (including periodic dependence), another averaging technique can be used in linear problems. This technique is based on the possibility of the integral representation of the solution of the original problem for the linear equation with variables coefficients in terms of the solution of the same problem for an equation of the same type but with constant coefficients [11–13]. The integral representation implies that the solution of the original problem can be represented in the form of the series in the gradients of the solution of the problem for the equation with constant coefficients [13].The aim of the present paper is to develop methods for calculating effective characteristics of corrugated plates. To this end, we first write out the equilibrium equations for a flexible anisotropic plate, which is inhomogeneous in the thickness direction and in the horizontal projection, with an initial deflection. We write these equations in matrix form, which allows one to significantly reduce the length of the expressions and to simplify further calculations. After this, we average the initial matrix equations with variable coefficients. The averaging procedure implies the statement of problems such that, after solving them, we can calculate the desired effective characteristics. By way of example, we consider the case of a corrugated plate made of a homogeneous isotropic material whose corrugations are hexagonal in the horizontal projection. In this case, we obtain approximate expressions for the components of the effective tensors of flexural rigidity and longitudinal compliance and expressions for the effective plate thickness.     

具有分层损伤的不同加筋形式复合材料层合板的后屈曲性态研究

基于板的一阶剪切理论和V on-K arm an大挠度理论,分别推导了复合材料层合板和层合梁的几何非线性有限元列式,提出了含嵌入分层的复合材料加筋层合板在受压缩载荷作用下的后屈曲有限元分析方法,对在板厚方向具有不同分层位置的加筋板结构进行了有限元数值分析,研究了不同的加筋方式及筋的分布对具有分层损伤的复合材料加筋层合板的后屈曲性态的影响,所得结果对确定在压缩载荷作用下含损伤复合材料加筋层合板的剩余承载能力具有参考价值。     

考虑变形耦合的几何非线性空间梁单元

以"精确几何模型梁单元"为代表的很多几何非线性梁单元,在构造过程中分别对描述截面转动的转角和描述截面形心位置的位移进行了独立插值,由此引起了诸如运动学描述冗余和剪切闭锁等困难。其根本原因在于单元形函数没能体现细长梁中的变形耦合关系。本文对这类传统单元进行了改造,通过深入研究单元变形之间的内在联系,提出了一种变形场完全满足Bernoulli梁变形耦合关系的新单元,避免了构造过程中对转动矢量的插值,并通过数值算例检验了单元的有效性。