Post-buckling analysis of FGM annular sector plates based on three dimensional elasticity graded finite elements

In this article, post-buckling and non-linear bending analysis of functionally graded annular sector plates based on three dimensional theory of elasticity in conjunction with non-linear Green strain tensor is considered. In-plane normal compressive loads have been applied to either radial, circumferential, or all edges of annular sector plates. Material properties are graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of constituents while Poisson׳s ratio is assumed to be constant. The governing equations are developed based on the principle of minimum total potential energy and solved based on graded finite element method. Non-linear equilibrium equations are solved based on iterative Newton–Raphson method. The effects of material gradient exponent, different sector angles, thickness ratio, loading condition and two different boundary conditions on the post-buckling behavior of FGM annular sector plates have been investigated. Results denote that due to the stretching–bending coupling effects of the FGMs, the post-buckling behavior of movable simply supported FGM plates is not of the bifurcation-type buckling. Moreover, FGM annular sector plates subjected to uniaxial compression at radial edges show a non-linear bending behavior with unique and stable equilibrium paths following a flattening feature.     

Mode interaction in triple-bay prestressed stayed columns

Prestressed stayed columns are an innovative type of structural system where the compressive load-carrying capacity can be enhanced through prestressed external cable stays. A nonlinear analytical model for prestressed stayed columns with multiple crossarm systems along the column length, based on the Rayleigh–Ritz method, is presented that captures modal interactions for perfect geometries explicitly for the first time. It is demonstrated that the theoretical compressive strength enhancements under certain configurations can only be obtained at the expense of triggering a sequence of destabilizing bifurcations. This can give rise to dangerously unstable interactive post-buckling behaviour including so-called ‘mode jumping’ and ‘snaking’ phenomena. Parametric spaces where the system is most susceptible to the modal interactions are identified and it is for those configurations that the system is likely to be highly sensitive to initial imperfections. The model is validated using a nonlinear finite element model formulated within the commercial code ABAQUS and excellent comparisons are obtained.     

Canonical dual finite element method for solving post-buckling problems of a large deformation elastic beam

This paper presents a canonical dual mixed finite element method for the post-buckling analysis of planar beams with large elastic deformations. The mathematical beam model employed in the present work was introduced by Gao in 1996, and is governed by a fourth-order non-linear differential equation. The total potential energy associated with this model is a non-convex functional and can be used to study both the pre- and the post-buckling responses of the beams. Using the so-called canonical duality theory, this non-convex primal variational problem is transformed into a dual problem. In a proper feasible space, the dual variational problem corresponds to a globally concave maximization problem. A mixed finite element method involving both the transverse displacement field and the stress field as approximate element functions is derived from the dual variational problem and used to compute global optimal solutions. Numerical applications are illustrated by several problems with different boundary conditions.     

腹板V形加劲冷弯薄壁卷边槽钢柱试验研究

为了解腹板V形加劲冷弯薄壁卷边槽钢构件的受压性能,对8根短柱和中长柱进行了轴心受压试验研究。在试验之前首先用有限条和有限元方法对柱试件进行了详细的设计,以确保发生畸变屈曲失效。柱受压试验结果发现：所有试件均发生畸变屈曲失效破坏,同时在峰值荷载时中长柱试件有绕弱轴的弯曲产生,畸变屈曲和局部屈曲一样也有一定的屈曲后强度。在总结分析现有国内外规范计算理论的基础上,提出了基于直接强度法设计腹板V形加劲冷弯薄壁卷边槽钢轴压构件畸变屈曲承载力的计算方法。采用该建议方法计算所得结果与试验值和参数分析值吻合较好,其设计过程简便且安全可靠。     

A non-linear 3-D finite element analysis of soil failure with tillage tools

A non-linear 3-D finite element model was developed to study the soil failure under a narrow tillage blade. The weighted residual method was applied to formulate the finite element model. The Duncan and Chang hyperbolic stress-strain model was used in the analysis. This finite element model also takes into account friction at the soil-tool interface, and progressive and continuous cutting. A FORTRAN program was written to carry out the finite element analysis. The results provided soil forces, a progressive developed failure zone, displacement field and stress distribution along the tool surface. Tillage were conducted in the laboratory soil bin to verify soil forces from the finite element analysis. The comparison between the results from the finite element model and those from the soil bin tests was reasonably good.     

复合材料加筋板极限压缩承载能力可靠性分析

基于Abaqus软件用户子程序,利用渐进失效分析方法对复合材料加筋板极限压缩承载能力进行预测。算例分析表明,对于加筋板1和2,本文方法给出的极限压缩强度与实验结果的误差分别为2.53%和1.68%。在结构可靠性分析过程中,为提高计算效率,利用屈曲载荷与极限压缩强度之间的关系建立功能函数,只需对极限压缩承载能力进行一次分析。对于加筋板1和2,本文方法相对经典可靠性方法的计算误差分别为-1.04%和-1.01%,计算时间仅为经典可靠性方法的1.18%和1.66%。     

Shells with random geometric imperfections simulation — based approach

The paper concerns the stochastic analysis of the static response of materially and geometrically non-linear models of imperfect shell structures. Geometric imperfections of shells are described using random variables or random fields. The Monte Carlo method combined with a finite element program analysis is employed. The proposed approach is used to analyse the non-linear post-buckling behaviour of a shallow cylindrical shell. The numerically obtained critical load histograms allow to estimate the structure reliability. The accuracy of simulation-based approach is discussed.     

基于响应面法的复合材料连接结构可靠性分析

复合材料的力学性能参数比金属材料呈现出更大的分散性。本文提出采用加权线性响应面法对复合材料单钉连接结构进行可靠性分析。采用三维有限元模型进行复合材料单钉连接的确定性渐进失效分析。单向铺层的力学性能参数和强度参数均考虑为输入随机变量,用于可靠性分析。详细研究了不同失效准则和不同刚度退化方法的影响。结果表明,Olmedo-Santiuste失效准则与Camanho-Matthews刚度退化方法的组合能够给出与试验均值最接近的强度分析结果。另外,还仔细研究了响应面参数f的影响。结果表明,取f=1的线性加权响应面能够为复合材料单钉连接可靠性分析提供较好的统计特性逼近值。     

深海保温输油管后屈曲及屈曲传播的有限元分析

本文运用ABAQUS有限元软件模拟深海输油管道的后屈曲及屈曲传播现象．将深海输油管道视为内、外层为合金钢材,夹心为聚氨酯泡沫的输油管道．采用线弹性-线性硬化的本构关系,运用Riks分析法获得后屈曲平衡路径,以此模拟管道中的屈曲传播过程,并获得屈曲传播压力．通过数值算例,综合讨论了不同初始缺陷、不同径厚比、不同夹层与夹心厚度比对、不同的材料弹性模量对屈曲传播过程和压力的影响．最后将输油管退化为单层管,将传播压力的有限元结果与实验结果、 Palmer理论解比较．结果表明有限元值、实验结果、理论值三者比较吻合．验证了有限元分析的正确性．     

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.     

Stretch-induced stress patterns and wrinkles in hyperelastic thin sheets

Wrinkles are commonly observed in stretched thin sheets and membranes. This paper presents a numerical study on stretch-induced wrinkling of hyperelastic thin sheets based on nonlinear finite element analyses. The model problem is set up for uniaxial stretching of a rectangular sheet with two clamped ends and two free edges. A two-dimensional stress analysis is performed first under the plane-stress condition to determine stretch-induced stress distribution patterns in the elastic sheets, assuming no wrinkles. As a prerequisite for wrinkling, development of compressive stresses in the transverse direction is found to depend on both the length-to-width aspect ratio of the sheet and the applied tensile strain in the longitudinal direction. A phase diagram is constructed with four different distribution patterns of the stretch-induced compressive stresses, spanning a wide range of aspect ratio and tensile strain. Next, an eigenvalue analysis is performed to find the potential buckling modes of the elastic sheet under the prescribed boundary conditions. Finally, a nonlinear post-buckling analysis is performed to show evolution of stretch-induced wrinkles. In addition to the aspect ratio and tensile strain, it is found that the critical condition for wrinkling and the post-buckling behavior both depend sensitively on the sheet thickness. In general, wrinkles form only when both the magnitude and the distribution area of the compressive stresses are sufficiently large. The wrinkle wavelength decreases with increasing strain, in good agreement with the prediction by a scaling analysis. However, as the tensile strain increases, the wrinkle amplitude first increases and then decreases, eventually flattened beyond a moderately large critical strain, in contrast to the scaling analysis.     

APPLICATION OF STOCHASTIC FINITE ELEMENT METHOD TO STRENGTH AND STABILITY ANALYSIS OF EARTH DAMS

This paper applies the stochastic finite element method to analyse thestatistics of stresses in earth dams and assess the safety and reliability of the dams.Formulations of the stochastic finite element method are briefly reviewed and theprocedure for assessing dam's strength and stability is described.As an example,adetailed analysis for an actual dam-Nululin dam is performed.A practical method forstudying built-dams based on the prototype observation data is described.     

Compressive strength of fibre composites with random fibre waviness

The compressive strength of unidirectional long fibre composites is predicted for plastic microbuckling from a random two-dimensional distribution of fibre waviness. The effect of the physical size of waviness is addressed by using couple stress theory, with the fibre bending resistance scaling with the fibre diameter d. The predicted statistical distribution of compressive strength is found using a Monte Carlo method. An ensemble of fibre waviness profiles is generated from an assumed spectral density of waviness and the compressive strength for each such realisation is calculated directly by the finite element method. The average predicted strength agrees reasonably with practical values, confirming the hypothesis that microbuckles can be initiated by fibre misalignment. It is found that the probability distribution of strength is well matched by a Weibull fit, and the dependence of the Weibull parameters upon the spectral density of waviness is determined. For the practical range of fibre distributions considered, it is concluded that the strength depends mainly upon the root mean square amplitude of fibre misalignment, with the shape of the power spectral density function playing only a minor role. An engineering model for predicting the compressive strength is proposed, akin to weakest link theory for materials containing flaws. A specimen containing randomly distributed waviness is examined to locate regions of high-fibre misalignment. The strength of each of these weak regions is estimated from a look-up table derived from calculations with idealised circular or elliptical patches of waviness. The strength of the composite is given by the failure stress associated with the weakest such patch. For random distributions of waviness, the predictions using this engineering approach are in good agreement with the direct calculations of strength using the finite element method.     

受弯脱层层板后屈曲有限元分析

脱层是复合材料层板结构中主要的缺陷形式之一。当脱层层板受到压力载荷的作用会造成脱层的局部屈曲和扩展,从而使结构的强度和刚度大为降低。含脱层层板的弯曲问题包含了脱层的压缩问题,却比压缩问题更加复杂。本文对含穿透脱层层板在纯弯载荷作用下的后屈曲问题进行了基于一阶剪切层板理论的几何非线性有限元分析,运用虚裂纹闭合技术求解了纯弯载荷作用下的脱层尖端的能量释放率各型分量,并用脱层扩展判据求解了脱层起始扩展载荷。     

含分层损伤复合材料层合板的压缩强度研究

给出了基于一阶剪切变形理论的含分层损伤层合板有限元分析模型,将含分层损伤层合板在压缩载荷作用下的强度破坏分析和屈曲破坏分析统一起来。先区分其破坏形式,然后再进行具体破坏分析,在屈曲特性分析中考虑了铺层强度破坏引起的刚度折减的影响,数值结果表明,该文给出的方法和结论对含分层损伤复合材料层合板的设计更具参考价值。     

IMPROVEMENTS OF GEOMETRICALLY NONLINEAR ANALYSIS ALGORITHMS FOR SPATIAL FRAME STRUCTURES

以几何精确梁理论为基础,分别采用高阶拉格朗日插值和埃米特插值构造高精度空间梁单元。提出基于单元层次平衡迭代的自由度凝聚方法,以保证单元的通用性。实现了基于载荷控制或柱面弧长控制的结构几何非线性分析算法。算例研究结果表明,提出的改进方法不但提高了计算效率,而且还具有较高的数值稳定性;特别是基于三次埃米特插值构造的单元表现出较好的性态,适用于结构屈曲后分析。     

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.     

结构可靠度分析的改进BP神经网络响应面法

对功能函数不能明确表达即具有隐式功能函数的问题进行可靠度分析,常采用响应面法。其中二次多项式响应面法应用较为广泛,采用与此方法相同的思路,人们提出了BP神经网络响应面法。在此基础之上,笔者提出了改进的BP神经网络响应面法。对以上方法计算效果,通过算例进行了对比分析,其中改进BP神经网络响应面法计算精度较好,进行有限元分析次数较少。该方法用于大型复杂结构的可靠度分析,可相应提高工作效率和解题质量,具有一定实际应用价值。     

Numerical study of strength properties for a composite material with short reinforcing fibers

A numerical approach to the determination of strength properties for concrete with short reinforcing fibers on the basis of the finite element method is proposed. The mathematical model takes into account various processes of nonlinear deformation of the concrete matrix under compressive and tensile loads, the possibility of developing the inelastic strains in the concrete matrix and reinforcing fibers and the nonlinear interaction between them. The effect of fiber concentration, various loading surfaces for the material matrix, and the bonding type on the deformation of a composite material is analyzed. Numerical examples of strength analysis are given.     

Secondary buckling of a flat plate under uniaxial compression : Part 1: theoretical analysis of simply supported flat plate

While studies of post-buckling behavior and load-carrying capacities of thin plates subjected to uniaxial compression have been limited to stable conditions, further post-buckling loading generates an unstable condition. The secondary buckling which occurs with snap-through to higher-order deflections under such unstable conditions has not been analyzed in detail as yet. In the first part of this paper, a thin square plate under uniaxial compression, which is simply supported along four edges, is considered. A method based on the second variation of the total potential energy is then proposed for evaluating the stability of the post-buckling equilibrium state and inevitable secondary buckling is derived analytically. The effects of various factors, such as initial imperfections, assumed virtual displacement pattern, post-buckling deflection pattern and in-plane boundary conditions, on the secondary buckling values are discussed. In part 2, secondary buckling of clamped plates is analyzed by use of the finite element method and the resultant numerical results are compared with experimental results.