含分层损伤复合材料层合板分层扩展研究

采用基于Ｍｉｎｄｌｉｎ－阶剪切理论的四节点板单元,分析了含椭圆分层合板分层扩展行为。利用虚裂纹闭合技术计算分层前缘处的总能量释放率,并采用总能量释放率准则作为扩展准则,结合自适应网格移动技术,并考虑了分层前缘闭合接角效应,对 合材料层合板的分层扩展行为进行了模拟分析。结果表明,初始分层形状对其扩展有方式有限大影响。     

Interlaminar stresses and energy release rates in delaminated orthotropic composite plates

The classical laminated plate theory is applied to calculate the stresses and energy release rate function in symmetrically delaminated orthotropic plates. First, the equilibrium of classical plate forces, moments and interfacial shear stresses is formulated. Second, the displacement continuity between the interface plane of a double-plate model was considered. The governing equation system of the double-plate model consists of ten equations. As an example a delaminated, orthotropic, simply-supported plate subjected to a point force is analyzed. The distribution of the plate forces as well as the interlaminar shear stresses over the uncracked part were determined. Moreover, the mode-II and mode-III energy release rate distributions along the crack front were calculated by the J-integral. The 3D finite element model of the delaminated composite plate was also created. The results indicate a reasonably good agreement between analysis and numerical calculation.     

Solutions of beam-shaped-function for analysis of composite plates with embedded delaminations

Based on beam-shaped-function, the analytical solution for composite plates with arbitrary embedded delaminations is presented. The deflection function of the delaminated plate is composed by those of beams with the corresponding loading and support conditions, which can be easily and accurately derived from the beam analysis, and the deflection amplitude is derived by the minimum potential energy principle. The closed form solutions of displacements, stresses, and energy release rate of a composite plate containing an arbitrarily embedded rectangular delamination are obtained and compared with the three-dimensional finite element results to validate the accuracy of this present method. Furthermore, the influences of delamination depth, length, central position, and modulus mismatch ratio (E ₁/E ₂) of the upper and lower sublaminate on the energy release rate are discussed.     

板结构计算模型的新发展

现代复合材料层合板具有高强和轻型的突出优点,从而在军工和民用等诸多领域发挥着重要作用。这种板结构的特点是随着纤维走向的不同,层间材料的物理-力学特性发生剧烈变化。沿板厚方向变形的梯度比较陡峭,并在层间结合面处发生强不连续,呈现zig-zag （锯齿状）现象。这导致横向剪应变在板的静态和动态响应中发生重要作用,不计横向变形的经典组合板计算模型CLPT难以适应现代多层板计算分析的需要。考虑横向剪切变形影响的板的计算模型得到重视和发展。需要指出,现有各种考虑剪切变形影响的计算模型虽然有了很大的发展,但在全面和准确性上仍然存在一定的不足,难以适应现代多层组合板横向力和物理性能多变的情况。模型预测的沿板厚方向位移和应力的变化规律难以通过严格的检验。本文提出的以比例边界有限元为基础的正交各向异性板的数值计算模型,同时可适用于各种薄板与厚板的分析,对现代复合材料层合板的分析具有特殊的优越性。所得到的板的位移、正应力和剪应力沿板厚方向的变化,与三维弹性理论的标准解高度吻合。数值算例进一步表明,随着层间纤维走向的变化,板内位移场和应力场沿板厚方向剧烈变化所呈现的锯齿现象均可以精准地进行模拟。据此,本文建议方法对现代板分析的广泛适应性和高度准确性得到了充分论证。     

Numerical analysis of delamination growth for stiffened composite laminated plates

A study of postbuckling and delamination propagation behavior in delaminated stiffened composite plates was presented. A methodology was proposed for simulating the multi-failure responses, such as initial and postbuckling, delamination onset and propagation, etc. A finite element analysis was conducted on the basis of the Mindlin first order shear effect theory and the von-K6rm~n nonlinear deformation assumption. The total energy release rate used as the criteria of delamination growth was estimated with virtual crack closure technique (VCCT). A self-adaptive grid moving technology was adopted to model the delamination growth process. Moreover, the contact effect along delamination front was also considered during the numerical simulation process. By some numerical examples, the influence of distribution and location of stiffener, configuration and size of the delamination, boundary condition and contact effect upon the delamination growth behavior of the stiffened composite plates were investigated. The method and numerical conclusion provided should be of great value to engineers dealing with composite structures.     

改进的虚裂纹闭合技术及其在复合材料脱层分析中的应用

提出一种改进的虚裂纹闭合技术（ＶＣＣＴ）对轮胎带束（一咱复合材料结构）不同脱层长度下的应变能释放率进行了计算。将强度裂纹闭合的位移作了一种可变约束施加在变形后的裂纹尖端,约束反力作了裂纹尖端力计算,应变能释放率由裂纹尖端力和约束位移的乘积得到。该方法的优点是简便易行,在物理理解上直观,实施起来方便。本文就用该技术对带束脱层进行了断裂力学分析,给出数值结果,并对带束的脱层寿命进行了评价。     

A crack tip tracking algorithm for cohesive interface element analysis of fatigue delamination propagation in composite materials

A novel approach is proposed for the use of cohesive elements in the analysis of delamination propagation in composite materials under high-cycle fatigue loading. The method is applicable to delamination propagation within the Paris-law regime and is suitable for the analysis of three-dimensional structures typical of aerospace applications. The major advantages of the proposed formulation are its complete independence of the cohesive zone length – which is a geometry-dependent parameter – and its relative insensitivity to mesh refinement. This is only possible via the introduction of three nonlocal algorithms, which provide (i) automated three-dimensional tracking of delamination fronts, (ii) an estimation of direction of crack propagation and (iii) accurate and mesh-insensitive integration of strain energy release rate. All calculations are updated at every increment of an explicit time-integration finite element solution, which models the envelopes of forces and displacements with an assumption of underlying constant cyclic loading. The method was implemented as a user-defined subroutine in the commercial finite element software LS-Dyna and supports the analysis of complex three-dimensional models. Results are presented for benchmark cases such as specimens with central cut plies and centrally-loaded circular plates. Accurate predictions of delamination growth rates are observed for different mesh topologies in agreement with the Paris-laws of the material.     

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.     

复合材料层压板低速冲击响应尺度效应数值模拟研究

为了研究尺度效应对于复合材料层压板在低速冲击作用下的动态响应和冲击损伤的影响,基于相似理论,建立了三种不同尺寸的层压板受冲击的三维有限元模型。在该模型中,针对层压板的面内损伤,采用改进的Chang-Chang准则进行预测;针对层压板内层间分层损伤,则使用Cohesive界面单元进行模拟。一旦复合材料层压板在低速冲击作用下产生损伤,则对出现损伤的区域进行材料参数退化。采用该模型对三种不同尺寸的层压板的冲击过程进行有限元分析,并将不同冲击速度下的冲击响应进行比较,得出了如下结论：在层压板内未发生冲击损伤时,冲击产生的挠度和冲击力与相似理论解十分吻合,一旦出现冲击损伤,则冲击力的变化与相似理论解有所差别;如果两个缩放模型的冲击速度之比等于缩放比例的平方根,则两个模型中的相对分层尺寸基本是相同的,这个结果与已有的实验结果吻合;而对冲击后面内损伤的分析表明,其损伤尺寸不符合这一相似规律。     

Bending solution of third-order orthotropic Reddy plates with asymmetric interfacial crack

In this paper Reddy’s third-order shear deformable plate theory is applied to asymmetrically delaminated orthotropic composite plates under antiplane–inplane shear fracture mode. A double-plate system is utilized to capture the mechanical behavior of the uncracked plate portion. An assumed displacement field is used and modified in order to satisfy the traction-free conditions at the top and bottom plate boundaries. Moreover, the system of exact kinematic conditions was also implemented into the novel plate model. An important improvement of this work compared to previous papers is the continuity condition of the shear strains at the interface of the double-plate system. Applying these conditions it is shown that the nineteen parameters of the third-order displacement field can be reduced to nine. Using the simplified displacement field the governing equations are derived, as well. The solution of a simply-supported delaminated plate is presented using the state-space model and the displacement, strain and stress fields are determined, respectively. The energy release rate and mode mixity distributions are calculated using the 3D J-integral. The analytical results are compared to those by finite element computations and it is concluded that the present model is the most accurate one among the previous plate theory-based approaches.     

考虑横法向热应变的C~0型Reddy板理论和三角形板单元

考虑横法向热变形,建议了C0型Reddy理论,并用于分析复合材料层合/夹层板热膨胀问题。虽然考虑了横法向热应变,但不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,构造有限元时仅需C0插值函数。基于这一模型,运用虚位移原理推导了复合材料板平衡方程以及构造了6节点三角形板单元,并分析了简支复合材料层合/夹层板的热膨胀问题。数值结果表明,建立的模型能准确分析复合材料层合/夹层板热膨胀问题,而忽略横法向热应变的理论分析热膨胀问题误差较大。     

Elastodynamic behavior of laminated orthotropic plates under initial stress

A finite element method of analysis of the vibrational and wave propagational characteristics is presented for a laminated orthotropic plate under initial stress. The plate may have an arbitrary number of bonded elastic orthotropic layers, each with distinct thickness, density and mechanical properties, and the analysis is capable of treating a completely arbitrary three-dimensional state of initial stress. Biot's theory for incremental elastic deformations of a stressed solid forms the basis for this study. A homogeneous, isotropic plate under two different states of initial stress was analyzed and their numerical results showed excellent correlation with those from an exact solution. Further examples of a three layer composite plate and a sandwich plate are offered to add some general insight to the physical behavior of such plates.     

On crack-tip strain energy release rate in non-principal directions of elasticity for simple layer plate of composite materials

In this paper,the fracture problem in non-principal directions of elasticity for a simple layer plate of linear-elastic orthotropic composite materials is studied.The formulae of transformation between characteristic roots,coefficients of elastic compliances in non-principal directions of elasticity and corresponding parameters in principal directions of elasticity are derived.Then,the computing formulae of strain energy release rate under skew-symmetric loading in terms of engineering parameters for principal directions of elasticity are obtained by substituting crack-tip stresses and displacements into the basic formula of the strain energy release rate.     

A hybrid stress approach for laminated composite plates within the First-order Shear Deformation Theory

This paper presents a hybrid stress approach for the analysis of laminated composite plates. The plate mechanical model is based on the so called First-order Shear Deformation Theory, rationally deduced from the parent three-dimensional theory. Within this framework, a new quadrilateral four-node finite element is developed from a hybrid stress formulation involving, as primary variables, compatible displacements and elementwise equilibrated stress resultants. The element is designed to be simple, stable and locking-free. The displacement interpolation is enhanced by linking the transverse displacement to the nodal rotations and a suitable approximation for stress resultants is selected, ruled by the minimum number of parameters. The transverse stresses through the laminate thickness are reconstructed a posteriori by simply using three-dimensional equilibrium. To improve the results, the stress resultants entering the reconstruction process are first recovered using a superconvergent patch-based procedure called Recovery by Compatibility in Patches, that is properly extended here for laminated plates. This preliminary recovery is very efficient from the computational point of view and generally useful either to accurately evaluate the stress resultants or to estimate the discretization error. Indeed, in the present context, it plays also a key role in effectively predicting the shear stress profiles, since it guarantees the global convergence of the whole reconstruction strategy, that does not need any correction to accommodate equilibrium defects. Actually, this strategy can be adopted together with any plate finite element. Numerical testing demonstrates the excellent performance of both the finite element and the reconstruction strategy.     

含有分层损伤的复合材料加筋层合板的屈曲性态研究

基于Ｍｉｎｄｉｎ假定推导了考虑剪切的复合材料加筋层板的有限元列式,并在此基础上计算出筋间基板含嵌入分层以及筋与基板连接处含穿透分层的加筋层合板在受压缩载荷情况下的屈曲模式和临界力。本文所给出的有限元方法及结论对从事复合材料结构设计的工程人员具有参考价值。     

湿热环境对挖补修理后复合材料层合板振动特性的影响

通过对湿热环境下阶梯挖补胶接修理后碳纤维/双马树脂(T300/QY8911)层合板振动特性的分析,探究了挖补修理后碳纤维/双马树脂层合板在不同湿热效应下的振动特性。基于Mindlin-Reissner板理论与最小势能变分原理,在考虑湿热效应前提下,推导了正交各向异性复合材料层合板的本构关系及四边简支条件下的自由振动控制方程;其次,利用ABAQUS软件建立了三维黏弹性有限元挖补修理模型,应用Fortran语言将湿热耦合下的应力-应变本构关系编入UMAT二次开发文件中,采用三维八节点五自由度等参单元进行了离散。通过与参考文献对照,验证了胶接修理后T300/QY8911层合板有限元模型和编译的湿热环境下UMAT子程序的合理性与有效性。算例结果表明:湿热环境对挖补修理后层合板固有频率的影响比对完整层合板严重;温度对材料性能的影响大于湿度;相比于材料性能变化对振动特性的影响,湿热残余应力对胶接修理后层合板振动特性的影响占主导;参考环境下,胶接修理后层合板的固有频率增大;湿应力对层合板固有频率的影响大于热应力。     

Dynamic analysis of composite plate with multiple delaminations based on higher-order zigzag theory

In a recent paper, Cho and Kim [Journal of Applied Mechanics] proposed a higher-order cubic zigzag theory of laminated composites with multiple delaminations. The proposed theory is not only accurate but also efficient because it work with a minimal number of degrees of freedom with the application of interface continuity conditions as well as bounding surface conditions of transverse shear stresses including delaminated interfaces. In this work, we investigate the dynamic behavior of laminated composite plates with multiple delaminations. A four-node finite element based on the efficient higher-order zigzag plate theory of laminated composite plates with multiple delaminations is developed to refine the prediction of frequencies, mode shape, and time response. Through the dynamic version of the variational approach, the dynamic equilibrium equations and variationally consistent boundary conditions are obtained. Natural frequency prediction and time response analysis of a composite plate with multiple delaminations demonstrate the accuracy and efficiency of the present finite element method. To prevent penetration violation at the delamination interfaces, unilateral contact constraints by Lagrange multiplier method are applied in the time response analysis. The present finite element is suitable for the prediction of dynamic response of thick composite plates with multiple and arbitrary shaped delaminations.     

Analytical model for delamination growth during small mass impact on plates

An analytical model is presented for delamination initiation and growth and the resulting response during small mass impact on orthotropic laminated composite plates, which typically is caused by runway debris and other small objects. The solution is obtained by a fast stepwise numerical solution of a single integral equation. Delamination size, load and deflection history are predicted by extension of an earlier elastic impact model by the author. Good agreement is demonstrated in comparisons with finite element simulations and experiments.     

材料沿厚度方向的非线性对复合材料层板分层的影响

基于修正的Ｈａｈｎ－Ｔｓａｉ非线性本构关系,用准三维有限单元分析了由层间应力引起的材料沿厚度方向的非线性对复合材料层板分层后屈及能量释放率的影响,结果表明材料非线性对能量释放率的影响与后屈曲形态密切相关,其中对Ⅱ型分层能量释放率的影响较大,而分层尖端应力场降低,局部变形增大。     

用高阶位移模式分析复合材料层合板的自由振动

建立一种新的高阶位移模式：分层假设复合材料层合板的位移场，利用各层间应力及位移的连续性条件，导出了作为整体的层合板控制微分方程。最后，采用这一高阶剪切理论来分析复合材料矩形层合板的自由振动。数值计算采用有限差分法，并编制了计算程序，计算得出了六种不同边界条件、不同铺层、不同宽厚比工况下的矩形层合板的自振频率。同经典理论解进行比较，可以看出本文方法简便易行，且精度较高，可以在微机上实用。