A Refined Shear Deformation Theory of Bending of Isotropic Plates*

ABSTRACT

A nonlinear, in-plane displacement assumption is proposed, based on an undetermined variation df/dz of transverse shear strains through the plate thickness. A second-order ordinary differential equation for f(z) and two surface conditions, as well as a set of eighth-order partial differential equations and four associated boundary conditions, are derived from the principle of minimum potential energy. Coupling exists between the partial and ordinary differential equations. In the homogeneous solutions for the former, in addition to an interior solution contribution, there exist two edge-zone solution contributions, one of which induces self-equilibrated (in the thickness direction) boundary stresses. Three examples are calculated using the present theory. The last gives the stress couple and maximum-stress concentration factors at the free edge of a circular hole in a large bent plate. Numerical results for the examples are compared with those given by three-dimensional elasticity theory and several two-dimensional theories. It is found that the present theory can accurately predict nonlinear variations of in-plane stresses through the thickness of a plate.     

A new zig-zag theory and C<Superscript>0</Superscript> plate bending element for composite and sandwich plates

A higher-order zig-zag theory for laminated composite and sandwich structures is proposed. The proposed theory satisfies the interlaminar continuity conditions and free surface conditions of transverse shear stresses. Moreover, the number of unknown variables involved in present model is independent of the number of layers. Compared to the zig-zag theory available in literature, the merit of present theory is that the first derivatives of transverse displacement have been taken out from the in-plane displacement fields, so that the C⁰ interpolation functions is only required during its finite element implementation. To obtain accurately transverse shear stresses by integrating three-dimensional equilibrium equations within one element, a six-node triangular element is employed to model the present zig-zag theory. Numerical results show that the present zig-zag theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories. Moreover, it is convenient to obtain transverse shear stresses by integration of equilibrium equations, as the C⁰ shape functions is only used when implemented in a finite element.     

Efficient higher-order shear deformation theories for bending and free vibration analyses of functionally graded plates

In this paper, various efficient higher-order shear deformation theories are presented for bending and free vibration analyses of functionally graded plates. The displacement fields of the present theories are chosen based on cubic, sinusoidal, hyperbolic, and exponential variations in the in-plane displacements through the thickness of the plate. By dividing the transverse displacement into the bending and shear parts and making further assumptions, the number of unknowns and equations of motion of the present theories is reduced and hence makes them simple to use. Equations of motion are derived from Hamilton’s principle. Analytical solutions for deflections, stresses, and frequencies are obtained for simply supported rectangular plates. The accuracy of the present theories is verified by comparing the obtained results with the exact three-dimensional (3D) and quasi-3D solutions and those predicted by higher-order shear deformation theories. Numerical results show that all present theories can archive accuracy comparable to the existing higher-order shear deformation theories that contain more number of unknowns.     

A quadrilateral element based on refined global-local higher-order theory for coupling bending and extension thermo-elastic multilayered plates

In this paper a refined higher-order global-local theory is presented to analyze the laminated plates coupled bending and extension under thermo-mechanical loading. The in-plane displacement fields are composed of a third-order polynomial of global coordinate z in the thickness direction and 1,2–3 order power series of local coordinate ζ_k in the thickness direction of each layer, which is identical to the 1,2–3 global-local higher-order theory by Li and Liu [Li, X.Y., Liu, D., 1997. Generalized laminate theories based on double superposition hypothesis. Int. J. Numer. Methods Eng. 40, 1197–1212] Moreover, a second-order polynomial of global coordinate z in the thickness direction is chosen as transverse displacement field. The transverse shear stresses can satisfy continuity at interfaces, and the number of unknowns does not depend on the layer numbers of the laminate.Based on this theory, a quadrilateral laminated plate element satisfying the requirement of C¹ continuity is presented. By solving both bending and thermal expansion problems of laminates, it can be found that the present refined theory is very accurate and obviously superior to the existing 1,2–3 global-local higher-order theory. The most attractive feature of this theory is that the transverse shear stresses can be accurately predicted from direct use of constitutive equations without any post-processing method. It is also shown that the present quadrilateral element possesses higher accuracy.     

Non-linear static-dynamic finite element formulation for composite shells

Three non-linear finite element formulations for a composite shell are discussed. They are the simplified large rotation (SLR), the large displacement large rotation (LDLR), and the Jaumann analysis of general shells (JAGS). The SLR and the LDLR theories are based on total Lagrangian approach, and the JAGS is based on a co-rotational approach. Both the SLR and LDLR theories represent the in-plane strains exactly the same as Green's strain-displacement relations, whereas, only linear displacement terms are used to represent the transverse shear strain. However, a higher order kinematic through the thickness assumption is used in the SLR theory, which leads to parabolic transverse shear stress distribution compared to a first order kinematic through the thickness relationship used in the LDLR theory that leads to linear transverse shear stress distribution. Furthermore, the LDLR theory uses an Euler-like angle in the kinematics to account for the large displacement and rotation. The JAGS theory decomposes the deformation into stretches and rigid body rotations, where an orthogonal coordinate system translates and rotates with the deformed infinitesimal volume element. The Jaumann stresses and strains are used. Layer-wise stretching and shear warping through the thickness functions are used to model the three-dimensional behavior of the shell, where displacement and stress continuities are enforced along the ply interfaces. The kinematic behavior is related to the original undeformed coordinate system using the global displacements and their derivatives. Numerical analyses of composite shells are performed to compare the three theories. The commercial code ABAQUS is also used in this investigation as a comparison.     

On consistent plate theories

Summary  Applying the uniform-approximation technique, consistent plate theories of different orders are derived from the basic equations of the three-dimensional linear theory of elasticity. The zeroth-order approximation allows only for rigid-body motions of the plate. The first-order approximation is identical to the classical Poisson-Kirchhoff plate theory, whereas the second-order approximation leads to a Reissner-type theory. The proposed analysis does not require any a priori assumptions regarding the distribution of either displacements or stresses in thickness direction. Received 10 January 2002; accepted for publication 16 April 2002     

一种预测层合板层间剪应力的新后处理方法

层合板是航空航天领域典型的承力构件,过大的层间应力是导致其分层失效的主要原因.准确的层间应力预测往往依赖于三维平衡方程后处理方法(TPM).然而,该方法需要计算面内应力的一阶导,使得基于C0型板理论构造的线性单元无法使用TPM计算横向剪应力.本文在三维平衡方程后处理方法的基础上,提出了一种新后处理方法(NPM).新后处理方法通过虚功等效法消除了三维平衡方程后处理方法中产生的位移参数的高阶导.基于提出的新后处理方法和C0型板理论,仅需使用线性单元就可以预测层合板的横向剪应力.为了验证所提方法的有效性,本文基于修正锯齿理论(RZT)和所提方法构造了一种C0连续的三节点三角形线性板单元.数值算例表明,所提方法和三维平衡方程后处理方法具有相同的计算精度,提出的板单元能够准确高效地预测层合板的横向剪应力.此外,所提方法便于结合现有的有限元商用软件使用,基于商用软件中板壳单元获得的节点位移,使用新后处理方法极易获得准确的层间剪应力.     

An extension of generalized plane stress for problems with out-of-plane restraint

The standard concept of generalized plane stress is extended to obtain a new mathematical model for studying the effect of local out-of-plane displacement restraint on the in-plane stresses and displacements in thin plates. It is pointed out how this model may be used by the photoelastician, whose otherwise plane-stress experiment introduces an unavoidable out-of-plane restraint condition in the model, to obtain some estimate of the deviation to be expected between the results of his experiment and the actual plane-stress solution of the problem. In this way, the model may be applied to aid in the interpretation of a large class of two-dimensional photoelastic analyses involving the determination of stresses near rigid inclusions and rigid boundaries. The extended model is then applied to the problem of an annular disk subjected to thermal shrinkage and completely restrained at its outer boundary. In view of the simplicity of the model, the predicted radial and circumferential stress distributions agree remarkably well with existing photoelastic data. In contrast, results obtained from standard generalized plane-stress theory, which cannot account for the out-of-plane displacement restraint at the outer boundary, show substantial deviation from experimental values, especially near the restrained boundary.     

Energy Scaling of Compressed Elastic Films –Three-Dimensional Elasticity¶and Reduced Theories

We derive an optimal scaling law for the energy of thin elastic films under isotropic compression, starting from three-dimensional nonlinear elasticity. As a consequence we show that any deformation with optimal energy scaling must exhibit fine-scale oscillations along the boundary, which coarsen in the interior. This agrees with experimental observations of folds which refine as they approach the boundary. We show that both for three-dimensional elasticity and for the geometrically nonlinear Föppl-von Kármán plate theory the energy of a compressed film scales quadratically in the film thickness. This is intermediate between the linear scaling of membrane theories which describe film stretching, and the cubic scaling of bending theories which describe unstretched plates, and indicates that the regime we are probing is characterized by the interplay of stretching and bending energies. Blistering of compressed thin films has previously been analyzed using the Föppl-von Kármán theory of plates linearized in the in-plane displacements, or with the scalar eikonal functional where in-plane displacements are completely neglected. The predictions of the linearized plate theory agree with our result, but the scalar approximation yields a different scaling.     

A two-dimensional continuum theory for angle-ply laminates

A two-dimensional continuum theory of microstructure is developed for stress analysis of angle-ply laminates under in-plane loading. An example problem is used to evaluate the results of the theory against a reference solution obtained by the finite element method. The results are in satisfactory agreement; they also show that the in-plane stresses reach somewhat higher peak values than reported in previous literature.The theory is also presented in a simplified version, which is found to be adequate for predicting interlaminar stresses and in-plane stress resultants, but does not give acceptable results for the variation of in-plane stresses through the thickness of the laminations.     

Higher-order shear deformable theories for flexure of sandwich plates—Finite element evaluations

A simple isoparametric finite element formulation based on a higher-order displacement model for flexure analysis of multilayer symmetric sandwich plates is presented. The assumed displacement model accounts for non-linear variation of inplane displacements and constant variation of transverse displacement through the plate thickness. Further, the present formulation does not require the fictitious shear correction coefficient(s) generally associated with the first-order shear deformable theories. Two sandwich plate theories are developed: one in which the free shear stress conditions on the top and bottom bounding planes are imposed and another, in which such conditions are not imposed. The validity of the present development(s) is established through, numerical evaluations for deflections/stresses/stress-resultants and their comparisons with the available three-dimensional analyses/closed-form/other finite element solutions. Comparison of results from thin plate. Mindlin and present analyses with the exact three-dimensional analyses yields some important conclusions regarding the effects of the assumptions made in the CPT and Mindlin type theories. The comparative study further establishes the necessity of a higher-order shear deformable theory incorporating warping of the cross-section particularly for sandwich plates.     

AN ACCURATE TWO-DIMENSIONAL THEORY OF VIBRATIONS OF ISOTROPIC,ELASTIC PLATES

An infinite system of two-dimensional equations of motion of isotropic elastic plates with edge and corner conditions are deduced from the three-dimensional equations of elasticity by expansion of displacements in a series of trigonometrical functions and a linear function of the thickness coordinate of the plate. The linear term in the expansion is to accommodate the in-plane displacements induced by the rotation of the plate normal in low-frequency flexural motions. A system of first-order equations of flexural motions and accompanying boundary conditions are extracted from the infinite system. It is shown that the present system of equations is equivalent to the Mindlin’s first-order equations, and the dispersion relation of straight-crested waves of the present theory is identical to that of the Mindlin’s without introducing any corrections. Reduction of present equations and boundary conditions to those of classical plate theories of flexural motions is also presented.     

A layerwise C0-type higher order shear deformation theory for laminated composite and sandwich plates

A novel layerwise C⁰-type higher order shear deformation theory (layerwise C⁰-type HSDT) for the analysis of laminated composite and sandwich plates is proposed. A C⁰-type HSDT is used in each lamina layer and the continuity of in-plane displacements and transverse shear stresses at inner-laminar layer is consolidated. The present layerwise theory retains only seven variables without increasing the number of variables when the number of lamina layers are intensified. The shear stresses through the plate thickness derived from the constitutive equation of the present theory have the same shape as those calculated from the equilibrium equation. In addition, the artificial constraints are added in the principle of virtual displacements (PVD) and are certainly fulfilled through a penalty approach. In this paper, two C⁰-continuity numerical methods, such as the Finite Element Method (FEM) and Bézier isogeometric element (BIEM) are utilized to solve a discrete system of equations derived from the PVD. Several numerical examples with various geometries, aspect ratios, stiffness ratios, and boundary conditions are investigated and compared with the 3D elasticity solution, the analytical, as well as, numerical solutions based on various plate theories.     

An exact solution for the steady-state thermoelastic response of functionally graded orthotropic cylindrical shells

We analyze the steady-state response of a functionally graded thick cylindrical shell subjected to thermal and mechanical loads. The functionally graded shell is simply supported at the edges and it is assumed to have an arbitrary variation of material properties in the radial direction. The three-dimensional steady-state heat conduction and thermoelasticity equations, simplified to the case of generalized plane strain deformations in the axial direction, are solved analytically. Suitable temperature and displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the thermoelastic equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which are then solved by the power series method. In the present formulation, the cylindrical shell is assumed to be made of an orthotropic material, although the analytical solution is also valid for isotropic materials. Results are presented for two-constituent isotropic and fiber-reinforced functionally graded shells that have a smooth variation of material volume fractions, and/or in-plane fiber orientations, through the radial direction. The cylindrical shells are also analyzed using the Flügge and the Donnell shell theories. Displacements and stresses from the shell theories are compared with the three-dimensional exact solution to delineate the effects of transverse shear deformation, shell thickness and angular span.     

Fatigue-crack-tip plastic strains by the stereoimaging technique

The stereoimaging technique is an accurate, high-resolution means of measuring the in-plane displacements resulting from the deformation of a specimen so that the corresponding components of the strain tensor can be computed independently of the stresses. The example used in this paper is a fatigue-cracked specimen of a microscopically homogeneous experimental powder-metallurgy aluminum alloy, analyzed to determine the displacement and strain fields accompanying the opening of the fatigue crack. The displacement measurements are processed by a computer program which compensates for measurement fluctuations in the displacement data by smoothing, and derives the strain magnitudes. The principal strains and the maximum shear strain are determined using Mohr's circle, and the latter strain is then used to estimate the plastic-zone size. The crack-opening mode may be inferred from the displacement map, and the state of stress (plane stress or plane strain) inferred by applying the in-plane compatibility equation.     

Photoelastic stress analysis of an elliptical hole in a thick plate subjected to uniform in-plane compressive Loading

A three-dimensional photoelastic analysis was conducted to determine the stress distribution and concentration around the periphery of a centrally located elliptical hole in a plate of finite thickness. The edge of the plate was subjected to a uniformly distributed compressive uniaxial in-plane load. The principle of superposition was employed to study the effect of uniform biaxial loading.Elliptical holes with five different major/minor axis ratios () ranging from 1.0 to 2.64 were investigated. Among the results of this study, it was established that the variation of the principal stresses at the edge of the hole is not linear across the plate thickness. It was also found that in loading the plate in a direction parallel to the major axis of the ellipse, the value of the maximum tangential principal stress () occurs in a plane other than the middle plane of the plate. However, in loading the plate in a direction either parallel or perpendicular to the major axis, the maximum transverse stress ( _z ) occurs at the middle plane. In addition, the maximum value of ( _z ) was about 20 percent of the maximum value of the tangential stress for all models tested. Furthermore, the effect of the bixial loading has reduced the value of the maximum tangential stress at the periphery of the hole as compared with uniaxial loading.As a three-dimensional theoretical solution does not exist for this problem, the present findings were correlated with the well established two-dimensional solutions.     

Plane-stress fracture testing of finite sheets under biaxial loads

A procedure which combines the Williams series-type stress- and displacement-field expressions at the crack-tip neighborhood with a suitable numerical scheme away from the crack-tip was employed in the determination of the plane-stress fracture properties of four finite 7076-T6 aluminum sheets containing cracks emanating from a circular hole under four biaxial loads. The compatibility of the analytical and numerical displacements at the nodal points along the boundary of the crack-tip neighborhood was utilized in formulating displacement-continuity expressions containing some undetermined constants which solution depends on the nature of the boundary loading conditions. By linear superposition of the displacement due to remote uniaxial load and the displacements due to remotely applied transverse load in the neighborhood of the crack-tip, biaxial-displacement-continuity expressions containing these important fracture properties—namely, the opening Mode I stress-intensity factorK, the nonsingular stress term associated with the stresses in the direction parallel to the plane of cracksA and the integration termB associated with the displacement in this direction—were evaluated. Because no known biaxial testing of this geometry had been reported prior to this research, the analytical procedure was used to select the optimum geometry required in a biaxial fracture test of a finite-sheet specimen containing cracks emanating from a circular hole. This geometric optimization of the specimen guaranteed uniformity of stress all over the volume of specimen and also made the alteration of the existing MTS test fixtures unnecessary. Four square sheets of 7075-T6 aluminum alloy containing a central hole with two collinear cracks emanating radially at the edge of the hole were then fabricated in accordance with the analytically determined geometric requirements. The biaxial fracture test was then conducted under four biaxial load factors (λ) of 0.0, 0.5, 1.0 and 1.5. The fracture toughness obtained in this research was compared with those reported for uniaxial loading of large panels. It was found that there is a good correlation between the reported fracture toughness and this work.     

An advanced higher-order theory for laminated composite plates with general lamination angles

This paper proposes a higher-order shear deformation theory to predict the bending response of the laminated composite and sandwich plates with general lamination configurations.The proposed theory a priori satisfies the continuity conditions of transverse shear stresses at interfaces.Moreover,the number of unknown variables is independent of the number of layers.The first derivatives of transverse displacements have been taken out from the inplane displacement fields,so that the C 0 shape functions are only required during its finite element implementation.Due to C 0 continuity requirements,the proposed model can be conveniently extended for implementation in commercial finite element codes.To verify the proposed theory,the fournode C 0 quadrilateral element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate.Numerical results show that following the proposed theory,simple C 0 finite elements could accurately predict the interlaminar stresses of laminated composite and sandwich plates directly from a constitutive equation,which has caused difficulty for the other global higher order theories.     

厚度效应对梁冲击响应的影响

用一种半解析法——间接模态叠加法,研究了质点与弹性力学梁的冲击问题,这种方法避免了具有未知奇异载荷项的平衡微分方程求解问题。由于可以用解析方法得到简支弹性力学梁的模态函数,并且能够以显式形式给出其频率方程,因此以质点与简支弹性力学梁的冲击问题为例,来考察厚度效应对瞬态响应的影响,并将所得结果与用Timoshenko梁理论所得结果进行了比较,说明了厚度效应在梁冲击问题中的重要影响。讨论了纵波和剪切波对撞击力等动力响应的影响。     

Micron-Scale Residual Stress Measurement by Micro-Hole Drilling and Digital Image Correlation

This paper reports a new technique, namely the incremental micro-hole-drilling method (IμHD) for mapping in-plane residual or applied stresses incrementally as a function of depth at the micron-scale laterally and the sub-micron scale depth-wise. Analogous to its macroscale counterpart, it is applicable either to crystalline or amorphous materials, but at the sub-micron scale. Our method involves micro-hole milling using the focused ion beam (FIB) of a dual beam FEGSEM/FIB microscope. The resulting surface displacements are recorded by digital image correlation of SEM images recorded during milling. The displacement fields recorded around the hole are used to reconstruct the stress profile as a function of depth. In this way residual stresses have been characterized around a drilled hole of 1.8microns. diameter, enabling the profiling of the stress variation at the sub-micron scale to a depth of 1.8 microns. The new method is used to determine the near surface stresses in a (peened) surface-severe-plastically-deformed (S²PD) Zr₅₀Cu₄₀Al₁₀ (in atomic percent, at.%) bulk metallic glass bar. In plane principal stresses of -800 MPa ± 90 MPa and −600 MPa ± 90 MPa were measured, the maximum compressive stress being oriented 15° to the axis of the bar.