梁杆结构二阶效应分析的一种新型梁单元

推导了一种计及梁杆二阶效应的新型两结点梁单元。首先依据插值理论构造了三结点Euler-Bernoulli梁单元的位移场:使用五次Hermite插值函数建立梁单元的侧向位移场,二次Lagrange插值函数建立梁单元的轴向位移场,进而由非线性有限元理论推导了单元的线性刚度矩阵和几何刚度矩阵,然后使用静力凝聚方法消除三结点梁单元中间结点的自由度,从而得到一种考虑轴力效应的新型两结点梁单元。实例分析表明,此新型梁单元具有很高的计算精度,使用此单元进行梁杆结构分析可获得相当准确的二阶位移和内力。     

Residual stress release characteristics of hole-drilling determined by in-plane three-directional optical interference moiré

By adopting a new stress calibration method, a grating rosette moiré interferometry and incremental hole-drilling combined system is developed to determine the magnitud of the residual stress in aluminium plate subjected to a uniform uniaxial tensile load. Performing in-plane three-directional fringe analysis, the optical data contained in the moiré interferograms are converted into values of the strains in three directions corresponding to the grating rosette. The evaluation is carried out through the measurement of the in-plane displacement field in three directions generated by the introduction of the small incremental hole. The in-plane three-directional displacement fields are determined from the calculation of the optical in-plane three-directional phase distribution by means of a phase-shifting method. A new finite element calibration analysis that is general axisymmetric elements instead of 3-D block element and harmonic axisymmetric element is adopted to relate the relieved displacement field to the magnitude of the residual stress. The magnitude of the principal stresses is finally evaluated through a least-squares calculation and is also compared with the stress value applied to the specimen measured with strain gauges.     

The curved beam/deep arch/finite ring element revisited

In this paper, an attempt is made to understand the errors arising in curved finite elements which undergo both flexural and membrane deformations. It is shown that with elements of finite size (i.e. a practical level of discretization at which reasonably accurate results can be expected), there can be errors of a special nature that arise because the membrane strain fields are not consistently interpolated with terms from the two independent field functions that characterize such a problem. These lead to errors, described here as of the ‘second kind’ and a physical phenomenon called ‘membrane locking’. The findings here emerge from recent research on the effect of reduced integration on shallow curved beam elements and on the use of coupled displacement fields in finite rings. The failures which have occurred in earlier attempts to use independent polynomial displacement fields for curved elements may not have been due to neglect of rigid body motions or failure to achieve constant strain states, but because of locking due to spurious constraints. These emerge in the penalty limits of extreme thinness (an inextensional regime), when exact integration of the energy functional of an element based on low order independent interpolations for the in-plane and normal displacements is used. It seems possible to determine optimal integration rules that will allow the extensional deformation of a curved beam/deep arch/finite ring element to be modelled by independently chosen low order polynomial functions and which will recover the inextensional case in the penalty limit of extreme thinness without spurious locking constraints. The much maligned ‘cubic in w–lincar in u’ curved beam element is now reworked to show its excellent behaviour in all situations. What is emphasized is that the choice of shape functions, or subsequent operations to determine the discretized functionals, must consistently model the physical requirements the problem imposes on the field variables. In this manner, we can restore an old element to respectability and thereby indicate clearly the underlying principles. These are: the importance of ‘field consistency’ so that arch and shell problems can be modelled consistently by independent polynomial displacement fields, and the role that reduced integration or some equivalent construction can play to achieve this.     

解析型几何非线性圆拱单元

为研究圆拱结构几何非线性内力和位移计算以及平面内稳定分析问题,基于圆拱几何非线性静力分析模型,构造了解析位移形函数,进而利用能量法和势能变分原理,构造了解析型几何非线性圆拱单元,给出了单元列式;同时将该文模型通过理论退化,得到了不考虑轴向变形的无压缩圆拱模型,进一步得到了解析型无压缩几何非线性圆拱单元,并将此单元应用于圆拱的平面内稳定分析。研究结果表明:采用本单元模型进行平面内分岔失稳分析得到的失稳临界荷载系数与经典Timoshenko模型、Dinnik模型及静力法模型计算结果一致,相对误差为0%;平面内分岔失稳的临界荷载与平面内极值点失稳的临界荷载一致。该文单元具有解析型、高精确性以及良好的适用性,可用于任意工况下圆拱结构几何非线性位移、内力和平面内稳定分析。     

Low-velocity impact analysis of composite laminates under initial in-plane load

In this paper, low-velocity impact response and damage of composite laminates is analytically investigated. A modified displacement field of the plate considering initially loaded in-plane strain is proposed. From the displacement field, a finite element equation on the structural behavior of composite laminate is newly induced and a computational program is coded. Numerical results using the FEM code is compared with the numerical ones from reference. Additional numerical analysis is performed on another impact condition, and the effect of initial in-plane load is investigated. Potential delamination damage area in the first inter-ply surface from the bottom of the laminate is approximately estimated, and the effect of the initial in-plane load and the impact condition are also investigated. Consequently, it may be concluded that the initial in-plane load of the laminate does not affect so much on the impact damage area of the laminate.     

Stability Analysis of Laminated Soft Core Sandwich Plates Using Higher Order Zig-Zag Plate Theory

C₀ finite element model based on higher order zig-zag plate theory is used to study the stability analysis of laminated sandwich plates. The in-plane displacement field is obtained by superposing a global cubically varying displacement field on a zig-zag linearly varying displacement field with different slope in each layer. The transverse displacement assumes to have a quadratic variation within the core and constant in the faces. The conditions regarding transverse shear stress at layer interfaces and top and bottom are satisfied. Numerical examples covering different features of laminated sandwich plates are presented to illustrate the accuracy of the model.     

Static stability of plates using fully conforming element

In this paper the geometric stiffness matrices of a 12 degree-of-freedom rectangular element based on a fully compatible displacement field are developed for the first time and used to solve the elastic stability problem of rectangular plates under combined in-plane forces. The rapid convergence of buckling solutions based on the derived geometric stiffness matrices is demonstrated for plates with various boundary and loading conditions. Comparison of results shows that the present conforming element gives in all cases better results than the well known 12 degree-of-freedom non-conforming rectangular element.     

A two‐noded locking–free shear flexible curved beam element

A new two‐noded shear flexible curved beam element which is impervious to membrane and shear locking is proposed herein. The element with three degrees of freedom at each node is based on curvilinear deep shell theory. Starting with a cubic polynomial representation for radial displacement (w), the displacement field for tangential displacement (u) and section rotation (θ) are determined by employing force‐moment and moment‐shear equilibrium equations. This results in polynomial displacement field whose coefficients are coupled by generalized degrees of freedom and material and geometric properties of the element. The procedure facilitates quartic polynomial representation for both u and θ for curved element configurations, which reduces to linear and quadratic polynomials for u and θ, respectively, for straight element configuration. These coupled polynomial coefficients do not give rise to any spurious constraints even in the extreme thin regimes, in which case, the present element exhibits excellent convergence to the classical thin beam solutions. This simple C⁰ element is validated for beam having straight/curved geometries over a wide range of slenderness ratios. The results indicates that performance of the element is much superior to other elements of the same class. Copyright © 1999 John Wiley & Sons, Ltd.     

基于一阶剪切变形理论的新型复合材料层合板单元

基于一阶剪切变形理论(FSDT),本文构造一种新型的20自由度(每结点5个自由度),四边形复合材料层合板单元,适合于任意铺设情形的层合板的计算。它是按如下方式构造的:(1) 单元每边的转角和剪应变由Timoshenko层合厚梁理论来确定;(2) 对单元域内的转角场和剪应变场进行合理的插值;(3) 引入平面内双线性位移场来体现层合板面内与弯曲的耦合作用。本文单元,记为TMQ20,不存在剪切闭锁现象,在计算单层的各向同性板时可以退化为文[1]中优质的中厚板单元TMQ。在文[2]中将给出本文单元对于层合板问题的详细数值算例。     

引入泡状位移含旋转自由度的广义协调三角形膜元

本文通过在单元位移场上增加含内参的广义泡状位移方法,构造了一个具有旋转自由度的三角形膜元。由于直接从完全多项式中应用广义协调方法确定附加位移,有效地扩大了单元位移场的解空间,所构造的单元具有计算精度高、对计算网格畸变不敏感的优良特性。本文为改善位移元世态提供了一个简单有效的方法。     

Modeling of functionally graded sandwich shells accounting for variation in transverse displacement

In this study, a simple C₀ isoparametric finite element formulation based on higher-order shear deformation theory is presented for static analysis of functionally graded material sandwich shells (FGMSS). To characterize the membrane-flexure behavior observed in a functionally graded shell, a displacement field involving higher-order terms in in-plane and transverse fields is considered. The proposed kinematics field incorporates for transverse normal deformation, transverse shear deformation, and nonlinear variation of the in-plane displacement field through the thickness to predict the overall response of the shell in an accurate sense. To develop the efficient C₀ formulation, the derivatives of transverse displacement are treated as independent field variables (nodal unknowns). Voigt's rule of mixture is employed to ascertain the mechanical properties of each layer's constituents along the thickness direction. A wide range of numerical problems are solved assuming various parameters: side-thickness ratio, curvature-side ratio, and gradation parameter, and their interactions with regard to static analysis of FGMSS are discussed in brief. Deflection and stresses incorporating different thickness schemes of sandwich shells are presented in the form of figures. To validate the results, a functionally graded shell without sandwich arrangement is considered. Since no results are available on static analysis of FGMSS, the present 2D model based on the finite element method might be helpful in assessing the applicability of other analytical and numerical models in this area in the future.     

Analytical and experimental studies on the low-velocity impact response and damage of composite laminates under in-plane loads with structural damping effects

In this paper, low-velocity impact response and damage of composite laminates under in-plane loads are analytically and experimentally investigated. The authors recently proposed a modified displacement field of plate theory, considering the effect of initially loaded in-plane strain, and used a finite element program to analyze the structural behavior of the composite laminate. In this study, the program is upgraded to account for the structural damping effect of the laminate. A pendulum type impact test system and an in-plane loading fixture are constructed for the experimental study. The analytical and experimental impact behaviors are compared at different impact energy levels for cases with an initial in-plane tensile load and a compressive load, as well as cases without the initial in-plane load. The results show good correspondence between the analytical and experimental impact force histories. The effect of the initial in-plane load reduces for higher impact energies. The numerical estimation of the damaged area is in good agreement with the results from C-scanning experiments.     

Fringe formation in multiaperture speckle shear interferometry

Multiaperture speckle shear interferometry exhibits sensitivity to in-plane displacement components. The response of the interferometer to in-plane displacement, however, depends on the location of the shear element on an aperture in front of the imaging lens. Two configurations with three apertures, with the shear element positioned at two different locations, are examined theoretically and experimentally.     

Consistent multiscale analysis of heterogeneous thin plates with smoothed quadratic Hermite triangular elements

A consistent multiscale formulation is presented for the bending analysis of heterogeneous thin plate structures containing three dimensional reinforcements with in-plane periodicity. A multiscale asymptotic expansion of the displacement field is proposed to represent the in-plane periodicity, in which the microscopic and macroscopic thickness coordinates are set to be identical. This multiscale displacement expansion yields a local three dimensional unit cell problem and a global homogenized thin plate problem. The local unit cell problem is discretized with the tri-linear hexahedral elements to extract the homogenized material properties. The characteristic macroscopic deformation modes corresponding to the in-plane membrane deformations and out of plane bending deformations are discussed in detail. Thereafter the homogenized material properties are employed for the analysis of global homogenized thin plate with a smoothed quadratic Hermite triangular element formulation. The quadratic Hermite triangular element provides a complete C¹ approximation that is very desirable for thin plate modeling. Meanwhile, it corresponds to the constant strain triangle element and is able to reproduce a simple piecewise constant curvature field. Thus a unified numerical implementation for thin plate analysis can be conveniently realized using the triangular elements with discretization flexibility. The curvature smoothing operation is further introduced to improve the accuracy of the quadratic Hermite triangular element. The effectiveness of the proposed methodology is demonstrated through numerical examples.     

Full-field assessment of the damage process of laminated composite open-hole tensile specimens. Part I: Methodology

The present paper is the first part of a study dedicated to the assessment of the damage taking place in composite open-hole tensile specimens using full-field strain measurements. It concentrates on the thorough validation of the methodology used to process the full-field data. First, the grid method used to provide the two components of the in-plane displacement field at the surface of the test specimens is described. Then, the data processing procedure to obtain the displacement maps during the test is discussed. The main part of the paper is dedicated to the derivation of strains from displacements, investigating two procedures: local differentiation and polynomial fitting. The effect of the grid is also investigated. Finally, successful experimental results are presented.     

Higher-order MITC general shell elements

Two mixed-interpolated general shell finite elements for non-linear analysis-a 9-node element and a 16-node element-are presented. The elements are based on the Mixed Interpolation of Tensorial Components (MITC) approach in which the covariant strain component fields for the in-plane and shear actions are interpolated and tied to the also interpolated displacement field. Both the 9-node element, referred to as the MITC9 element, and the 16-node element, referred to as the MITC16 element, are tested numerically and found to have high predictive capabilities.     

Nonlinear oscillations of laminated plates using an accurate four-node rectangular shear flexible material finite element

The objective of the present paper is to investigate the large amplitude vibratory behaviour of unsymmetrically laminated plates. For this purpose, an efficient and accurate four-node shear flexible rectangular material finite element (MFE) with six degrees of freedom per node (three displacements (u, v, w) along thex, y andz axes, two rotations (θ _x and θ _y ) abouty andx axes and twist (θ _xy )) is developed. The element assumes bi-cubic polynomial distribution with sixteen generalized undetermined coefficients for the transverse displacement. The fields for section rotations θ _x and θ _y , and in-plane displacementsu andv are derived using moment-shear equilibrium and in-plane equilibrium equations of composite strips along thex- andy-axes. The displacement field so derived not only depends on the element coordinates but is a function of extensional, bending-extensional coupling, bending and transverse shear stiffness as well. The element stiffness and mass matrices are computed numerically by employing 3×3 Gauss-Legendre product rules. The element is found to be free ofshear locking and does not exhibit any spurious modes. In order to compute the nonlinear frequencies, linear mode shape corresponding to the fundamental frequency is assumed as the spatial distribution and nonlinear finite element equations are reduced to a single nonlinear second-order differential equation. This equation is solved by employing thedirect numerical integration method. A series of numerical examples are solved to demonstrate the efficacy of the proposed element.     

A generalized higher-order theory for buckling of thick multi-layered composite plates with normal and transverse shear strains

The onset of buckling in square laminated multi-layered composite plates, subject to unidirectional in-plane loads, is investigated within the framework of a generalized higher-order shear deformation theory suitable to capture significant transverse shear and thickness-wise deformation effects. The displacement field is expanded in a Taylor series of the thickness coordinate with arbitrary polynomial degree; in turn, the series coefficients, expressed as a superposition of admissible functions, are determined according to the Rayleigh–Ritz method. Truly higher-order polynomial terms, along with a sufficient number of in-plane admissible functions, are shown to be necessary for convergence towards the fundamental buckling load multiplier. As a by-product, reduced-order models are identified for various plate geometries and lamination schemes. The sensitivity of the lowest buckling load with respect to the nondimensional parameters (the thickness ratio, the ratio between the elastic moduli, the ply angle) is investigated. In particular, the attention is focused on the cross-over phenomenon between the lowest two buckling eigenvalues in multi-layered composite square plates with different lamination schemes. The presented results shed light onto the buckling behavior of thick shear-deformable multi-layered plates.     

Simulation of tension fields with in-plane rotational degrees of freedom

This paper introduces a novel multigrid approach for the geometric non-linear simulation of tension fields on the basis of a three-node membrane finite element. The element possesses, in addition to the nodal displacement degrees of freedom, an in-plane rotational degree of freedom inside the element domain that controls the direction of the tension field. This rotational degree of freedom allows the enforcement of continuity and tension field boundary conditions on the basis of a coarser mesh with varying size.