Variational principle of hybrid energy and the fundamentals of 3-D laminate theory—A new approach for the analysis of interlaminar stresses in composite laminates

This paper discusses the discontinuity of stresses and strains at interlaminar surfaces of the composue laminate and presents a 3-D laminate theory for composite materials. This paper also presents a new type of elastic energy based on the globally continuous variables in laminates, different from the traditional potential energy and complementary energy. Then, a variational principle corresponding to the 3-D laminate theory is developed. The theory and the principle could be a basis of verifying the 2-D laminate theory and determining the interlaminar stresses near the free edges.     

厚热塑性复合材料AS4/PEEK角铺设层合板的层间应力分析

将表征热塑性复合材料AS4 /PEEK非线性行为和应变速率相关行为的三维弹塑性模型通过程序加以实现。将程序计算结果和文献实验结果相比较可以发现,二者吻合较好,验证了所生成程序的有效性。计算了厚的AS4 /PEEK角铺设层合板[±25]s4在不同界面上的层间应力。由层间应力的三维分布图,分析了不同界面上层间应力的分布特征,并说明了可能引起层间分层的主要因素。     

Interlaminar hygrothermal stresses in laminated plates

Within the elasticity formulation the most general displacement field for hygrothermal problems of long laminated composite plates is presented. The equivalent single-layer theories are then employed to determine the global deformation parameters appearing in the displacement fields of general cross-ply, symmetric, and antisymmetric angle-ply laminates under thermal and hygroscopic loadings. Reddy’s layerwise theory is subsequently used to determine the local deformation parameters of various displacement fields. An elasticity solution is also developed in order to validate the efficiency and accuracy of the layerwise theory in predicting the interlaminar normal and shear stress distributions. Finally, various numerical results are presented for edge-effect problems of several cross-ply, symmetric, and antisymmetric angle-ply laminates subjected to uniform hygrothermal loads. All results indicate high stress gradients of interlaminar normal and shear stresses near the edges of laminates.     

材料非线性对复合材料层合板热自由边界效应的影响

本文用准三维有限元法研究了材料非线性对复合材料层合板热自山边界效应的影响,给出了修正型Hahn-Tsai非线性应力-应变关系的三维形式。由本文非线性分析方法得到的层间应力与以往由线性分析方法得到的层间应力做了比较,结果表明:材料非线性能显著降低层间剪应力的集中程度,但对层间正应力影响不太明显。     

THREE DIMENSIONAL STRESS ANALYSIS OF SYMMETRIC COMPOSITE LAMINATES UNDER UNIAXIAL EXTENSION AND IN－PLANE PURE SHEAR

ＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＳＴＲＥＳＳＡＮＡＬＹＳＩＳＯＦＳＹＭＭＥＴＲＩＣＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＳＵＮＤＥＲＵＮＩＡＸＩＡＬＥＸＴＥＮＳＩＯＮＡＮＤＩＮ－ＰＬＡＮＥＰＵＲＥＳＨＥＡＲＣｈｉｃｎｗｅｉ－ｚａｎｇ（钱伟长）,Ｈｕａｎｇ...     

具有中心圆孔的[0°/90°]_s复合材料层合板的边界元法分析

本文在文[1]的基础上,采用子结构法建立了多层复合板的边界元方法,对具有中心园孔[0°/90°]_s的层合板的层间应力作了计算,同有限元法的结果进行了比较,结果表明,应用边界元法处理这类问题,单元划分少,节约了计算机时,而且有较高的计算精度。     

基于树脂层强度预测分层萌生载荷

 轴向拉伸下的层合板在自由边缘附近存在奇异应力场，容易产生分层萌生进而导致结构破坏．而大多数基于材料力学方法的分层萌生研究对强度参数的确定有较强的主观性，缺乏合理解释．本文通过在层间插入一定厚度的树脂层，将分层萌生视为层间树脂层在三维应力状态下的强度破坏，根据有明确试验标准的树脂强度和Mohr 判据判断是否发生分层萌生．采用Pipes-Pagano 的均匀轴向拉伸模型计算应力场，对层间应力与树脂层面外应力进行了对比．取临界长度为4 个单层板厚度对T800/914 层合板的分层萌生进行了预测，结果表明预测值与实验值吻合良好．     

Steady-state aeroelasticity of fluid flow over a laminated composite plate

A stress-function-based variational approach is used to determine the interlaminar stresses in a multilayered strip of laminate subject to arbitrary combinations of axial extension, bending, and the steady-state aerodynamic loading of fluid flow over the upper surface of laminated composite plate. Symmetric four-layer, cross-ply and angle-ply laminates are considered in details. Some numerical solutions by using a personal computer are obtained. The present results for four-layer laminates show that the aerodynamic loading has significant effects on the interlaminar stresses near the free-edge regions.     

Layup optimization with GA for tapered laminates with internal plydrops

Layup optimization of the maximum strength of laminated composites with internal ply-drops is performed by genetic algorithm (GA). Interlaminar stresses are considered in estimating the strength of laminates and calculated by the stress function based complementary virtual work principle. Out-of-plane stress functions are expanded in terms of harmonic series through the thickness direction and initially satisfied the traction free boundary conditions of laminates automatically. As the number of expansion terms is increased, stress concentration near the dropped plies is predicted with better accuracy. Since the proposed analysis is relatively simple and efficient in the prediction of interlaminar stress concentration near the ply-drops, the layup optimization of composite laminates with dropped plies considering interlaminar strength can be easily performed by GA. In the formulation of genetic algorithm, a repair strategy is adopted to satisfy given constraints and multiple elitism scheme is implemented to efficiently find multiple global optima or near-optima.     

Properties degradation induced by transverse cracks in general symmetric laminates

This paper presents the details of a methodology for predicting the thermoelastic properties degradation in general symmetric laminates with uniform ply cracks in some or all of the 90° layers. First, a stress transfer method is derived by using the concept of state space equation. The laminate can be subjected to any combination of in-plane biaxial and shear loading, and the uniform thermal loading is also taken into account. The method takes into account all independent material constants and guarantees continuous fields of all interlaminar stresses across interfaces between material layers. By this method, a laminate may be composed of an arbitrary number of monoclinic layers and each layer may have different material property and thickness. Second, the concept of the effective thermoelastic properties of a cracked laminate is introduced. Based on the numerical solutions of specially designed loading cases, the effective thermoelastic constants of a cracked laminate can be obtained. Finally, the applications of the methodology are shown by numerical examples and compared with numerical results from other models and experiment data in the literature. It is found that the theory provides good predictions of the thermoelastic properties degradation in general symmetric laminates.     

Numerical simulation of interlaminar damage propagation in CFRP cross-ply laminates under transverse loading

This paper proposes a numerical simulation of interlaminar damage propagation in FRP laminates under transverse loading, using the finite element method. First, we conducted drop-weight impact tests on CFRP cross-ply laminates. A ply crack was generated at the center of the lowermost ply, and then a butterfly-shaped interlaminar delamination was propagated at the 90/0 ply interface. Based on these experimental observations, we present a numerical simulation of interlaminar damage propagation, using a cohesive zone model to address the energy-based criterion for damage propagation. This simulation can address the interlaminar delamination with high accuracy by locating a fine mesh near the damage process zone, while maintaining computational efficiency with the use of automatic mesh generation. The simulated results of interlaminar delamination agreed well with the experiment results. Moreover, we demonstrated that the proposed method reduces the computational cost of the simulation.     

Thermal bending analysis of laminated orthotropic plates by the generalized differential quadrature method

The interlaminar stresses in a thin laminated rectangular orthotropic plate with four sides simply supported edges under bending was determined by using the generalized differential quadrature (GDQ) method involving the effects of thermal expansion strain and transverse load. The approximate stress and displacement solutions are obtained under the effects of thermal expansion force and uniform pressure load for eight-layer unidirectional laminates, symmetric cross-ply laminates. Numerical results on the dominant interlaminar stresses and displacement of bending analysis are compared to the Navier solution. The thermal induced forces have significant effect on the bending of plates.     

Free-edge interlaminar stress analysis of piezo-bonded composite laminates under symmetric electric excitation

A stress function-based approach is proposed to analyze the free-edge interlaminar stresses of piezo-bonded symmetric laminates. The proposed method satisfies the traction free boundary conditions, as well as surface free conditions. The symmetric laminated structure was excited under electric fields that can generate induced strain, resulting in pure extension in the laminated plate. The governing equations were obtained by taking the principle of complementary virtual work. To verify the proposed method, cross-ply, angle-ply and quasi-isotropic laminates were analyzed. The stress concentrations predicted by the present method were compared with those analyzed by the finite element method. The results show that the stress function-based analysis of piezo-bonded laminated composite structures is an efficient and accurate method for the initial design stage of piezo-bonded composite structures.     

Extended Kantorovich method for local stresses in composite laminates upon polynomial stress functions

The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work, the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is estab-lished to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method (FEM). The convergent stresses have good agreements with those results obtained by three dimensional (3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kan-torovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.     

Application of microline/grid methods to temperature dependent microscopic deformation and damage in CFRP laminates

Tensile tests of CFRP symmetric cross-ply laminates are carried out in a scanning electron microscope (SEM), and microscopic interlaminar deformation and damage near the transverse crack tip are visualized by microlines or microgrids printed on the specimen edge surface. The local deformation around the transverse crack tip is observed at 20°C, 80°C, 120°C and 160°C to evaluate the effect of thermal residual stress on the microscopic deformation and damage in the interlaminar region near the transverse crack tip. Temperature dependence of the axial crack opening displacement (COD) is also measured. The displacement field of the specimen edge surface obtained from these experimental results is compared with the theoretical model proposed by Lee, Allen and Harris. The analysis is modified to consider the temperature effect, and a good agreement is obtained between the modified theoretical predictions and experimental results.     

基于整体局部1,3位移模式的层合板层间应力研究

为提高层合板层间应力计算的准确性,对Reddy型高阶剪切理论的基本位移模式进行改进,提出整体-局部1,3高阶位移模式．在满足层间位移连续,层间剪切应力连续,以及上下表面自由的条件下,与前人提出的整体-局部1,2-3位移模式相比,层合板板结构每个节点的独立变量由13缩减到11,并且不随层数的增加而变化．将整体-局部1,3高阶位移模式位移和应力的数值解与解析解进行对比,验证了整体-局部1,3位移模式的准确性,可应用于复合材料层合板的位移和应力分析．     

On the description of anisotropic damage in composite laminates

A general anisotropic damage theory of cracked laminates is formulated here. The deformation of composite laminates is composed of matrix elastic strains, pseudo-elastic damage strains due to cracking and permanent damage strains due to interlaminar slip. The surface of damage initiation is constructed according to the concept of linear elastic fracture mechanics for the virgin material. After the initial damage, a pesudo-elastic damage can be used to describe the damage behaviour if interlaminar slip is negligible. Damage evolution, load induced anisotropy and interlaminar intralaminar interaction for composite laminates are examined; the latter can perturb the normality structure of damage strain rate. Explicit expressions are given for pseudo-elastic (or secant) moduli of the damaging composite laminates, under a non-interacting assumption imposed on the cracks between different families.     

Effects of material and stacking sequence on behavior of composite plates with holes

Strain distributions to failure, tensile and compressive strain-concentration factors, and strength-reduction factors were determined for glass-, boron-, and graphite-epoxy plates with holes loaded in tension. Strain gages, photoelastic coatings and moiré techniques were used. Ten variations of layup and stacking sequence were studied. The boron-epoxy composite was found to be the stiffest and strongest of the three. The graphite laminate with the highest stress concentration and the most linear strain response exhibited the highest strength-reduction factor. In all cases, the maximum strain at failure on the hole boundary was higher than the ultimate tensile-coupon strain. In general, it was found that, the higher the stress-concentration factor, the higher the strength-reduction factor. Thus, the [0/90/0/90]_s layup with a stress-concentration factor of 4.82 had a strength-reduction factor of 3.18. At the other extreme, the most flexible layup [±45/±45]_s with the lowest stress-concentration factor of 2.06 had the lowest strength-reduction factor of 1.10. Stacking sequences associated with the tensile interlaminar normal stress or high interlaminar shear stress near the boundary, resulted in laminates 10 to 20 percent weaker than corresponding alternate stacking sequences. Furthermore, it was found that stacking-sequence variations can alter the mode of failure from catastrophic to noncatastrophic.     

A STRESS VECTOR-BASED CONSTITUTIVE MODEL FOR COHESIONLESS SOIL (Ⅰ)-THEORY

On the basis of the sufficient consideration of vectorial characteristics of stress, a new nonlinear constitutive model for cohesionless soil under plane strain and 3-D conditions was presented in a way that the action effects of stress vector are decomposed into the action effect of mean effective stress and that of the stress ratio vector (ratio of deviatoric stress vector to mean effective stress). The constitutive model can take account of the influence of both numerical and directional changes of stress vector on deformation of soil simultaneously, and is applicable of both static and dynamic loading. Biography: Shi Hong-yan (1964-)     

A general nonlinear global-local theory for bending and buckling analyses of imperfect cylindrical laminated and sandwich shells under thermomechanical loads

The accurate shell theories proposed so far have been calibrated based on linear kinematic relations. Many of them have ignored either the interlaminar stress continuity conditions at the interfaces or the transverse flexibility of the layers. Therefore, the available shell theories may encounter accuracy problems when analyzing the nonlinear behaviors, especially for sandwich shells with soft cores. Moreover, almost all of the available shell theories have been proposed employing the Love-Timoshenko assumption. Ideas of the previous global-local plate theory of the author are extended to develop the present nonlinear high-order global-local shell theory. The present theory has the advantages of: (1) suitability for non-linear analyses, (2) higher accuracy due to satisfying the complete interlaminar kinematic and transverse stress continuity conditions at the layer interfaces under thermo-mechanical loads, employing the exact Green’s strain tensor of the curvilinear coordinates, considering the transverse flexibility, and releasing the Love-Timoshenko assumption, (3) less required computational time due to using the global-local technique and matrix formulations, and (4) capability of investigating the local phenomena. To enhance the accuracy of the results, compatible Hermitian elements are employed. Various comparative examples are included in the present paper to validate the theory and to examine its accuracy and efficiency.