对称复合材料层合板弯曲的三维数值分析

本文采用三维各向异性有限单元模拟纯弯曲载荷下的复合材料层合板，给出了应力、应变沿厚度方向的变化规律。分析结果表明，在斜交对称铺层层合板的中心区域（远离加载端和自由边的区域），板的层间粘接界面附近，有根强的应力集中现象，可称之为层间效应。层间界面处力学性质的突变导致了层间应力的产生，并使层合板处于三向应力状态。三维数值模型给出的应变分布不同于基于Ｋｉｒｃｈｈｏｆｆ－Ｌｏｖｅ直法线假设的经典层合板理论给出的应变分布。     

Hamilton体系下层合板弱粘接模型的应用

层合板的层间粘接模型对整个层合板的结构有重要影响．运用Hamilton正则方程对层合板层间的不同类型的粘接模型进行了分析．结合弹性材料修正后的Hellinger Reissner变分原理和插值函数,构建了直角坐标系下8节点层合板的每一层的线性方程;考虑到脱层板的连接界面处应力和位移的关系,改进了现有的常用弱粘接模型,建立不同粘接模型的控制方程;最后通过求解整个板的控制方程,得到层合板的层间应力和位移．数值算例验证了该模型的正确性,并研究了层间界面为线性和非线性时的问题．结果表明:应用改进后的弱粘接模型,能够更好地模拟层合板的弱界面失效过程．     

杂交能变分原理及层合板三维理论的基础——分析复合材料层间应力的新方法

本文讨论了复合材料层合板层间界面处应力和应变的间断性问题,并提出了复合材料层合板的三维理论.本文还以在层合板中全局连续的变量为基础,提出了与传统的势能和余能不同的一种新的弹性能.并得到与复合材料三维理论相应的变分原理.本文得到的三维层合板理论和杂交能变分原理,可以作为校核层合板二维经典理论和确定自由边界附近层间应力的理论基础.     

带裂纹层合板复应力强度分析

本文用经典板理论求得了复合材料层合板的界面应力强度因子.利用裂纹尖端能量释放率和复应力强度因子间的关系,给出最一般的复合材料层合板,在相应外荷载和模型混合参数下的复应力强度因子的一个封闭形式的解.然后提出确定这一模型混合参数的步骤,给出某些层合板的数值结果.并给出在相应外荷载下的小范围接触条件.特别证明了界面韧度曲线的对称性质.最后讨论振荡指数消失后,预计的断裂荷载的精确性.还通过一个实例表明β=0法的有效性和局限性.     

带裂纹层合板能量释放率分析

为了进一步了解裂纹尖端应力场的特性,本文对复合材料层合板的界面裂纹作了分析.文中强调了能量释放率分量存在的条件,并给出能量释放率分量和应力强度因子间的关系式.结合经典板理论的分析结果,根据外荷作用及某些几何参数和材料多数,导出了一般复合材料层合板的应力强度因子的封闭形式解.为了得到在一般荷载条件下能量释放率的分量,必须分别确定模型混合参数Ω,文中讨论了确定参数Ω的方法,最后,应用本文方法于几种不同种类的复合材料层合板,证明其结果可应用于工程实践.     

基于辛对偶体系的层合板自由边缘效应的分析解

在原变量——位移和其对偶变量——应力组成的辛几何空间,建立了Pipes-Pagano模型的复合材料层合板问题的辛对偶求解体系．与传统的单类变量不同,辛对偶变量有利于同时描述层间位移连续性条件和应力平衡条件．进入辛对偶体系以后,就可以应用辛对偶体系的统一解析求解方法,如分离变量和辛本征展开的方法对层合板问题进行解析分析和求解．对层合板自由边缘效应的分析求解,验证了辛对偶体系的方法对层合板问题的分析求解是十分有效的．     

具界面损伤压电智能层合板的非线性自由振动分析

基于广义6自由度板理论、应变等效原理和Hamilton变分原理,通过引入三维弹性平衡方程和静电平衡方程的通解来构造满足界面间力电耦合关系和各类连续条件的位移、电势分布形函数,建立了具铺设层内和层间界面处损伤效应的压电智能层合板的非线性运动控制方程组,并运用Galerkin 方法进行求解.数值算例中,分别讨论了,不同损伤程度、压电层厚度、厚跨比及长宽比对四边简支非理想界面压电智能层合板线性自由振动频率和非线性幅频响应曲线的影响.     

粘弹层合板的自由振动和横向应力

基于Reddy分层理论并在板厚方向取二次位移插值函数来推导出粘弹层合板的动力学方程,得到了简支粘弹夹层层合板的振动频率,其数值与已知文献数据吻合较好,且能够计算出协调的横向应力.在低频自由振动时,横向剪应力是造成粘弹层合板脱层的主要原因;高频时,横向正应力在脱层破坏中起主要作用.分析了粘弹材料模量对层合板横向应力的影响以及横向应力最大值与面内应力最大值的比值.结果表明所采用的算法、算式及所编写的程序是可靠的.     

关于层合板层间的连续条件

给出了保证满足复合材料层合板层间位移和应力连续条件的各种形式混合能量原理的泛涵,并通过算例证明了泛函的可靠性.     

复合材料层合板的弯曲与振动

本文应用应力杂交有限元方法分析了复合材料层合板的弯曲与振动.在本文中,首先根据修正的余能变分原理,构造了一个适合于复合材料层合板特点的矩形应力杂交板弯曲单元.在单元内,分层假设应力参数,在单元的边界上,根据YNS理论的假设确定边界位移场.这样使得构造出来的单元不仅能够考虑横向剪切变形的影响和局部扭曲效应,而且具有较少的自由度数.其次,用此单元求解了层合板的弯曲与振动问题,并将计算结果与精确解进行了比较,比较表明二者非常接近.这说明了在计算方面本文单元具有较高的精确度.     

A combined finite/discrete element algorithm for delamination analysis of composites

A combined finite/discrete element method is developed to model delamination behaviour in laminated composites. A penalty based algorithm is employed to evaluate the interlaminar stress state. The failure surface for delamination is defined by a Chang-Springer criterion, and the interlaminar crack propagation is achieved by a standard discrete element contact/release algorithm. The ability of the method for simulation of this behaviour is assessed by solving standard test cases available from the literature.     

Mode I,mode II,and mixed-mode I/II interlaminar fracture toughness of GFRP influenced by fiber surface treatment

The interlaminar fracture behavior of unidirectional glass fiber reinforced composites with fiber surface treatment has been investigated in modes I and II and for fixed mode I to mode II ratio of 1.33. The data obtained from these tests have been analyzed by using different analytical approaches. The present investigation is focused on the influence of the glass fiber surface treatment on the interlaminar fracture toughness of unidirectional laminates. Glass fibers with two different fiber surface treatments have been investigated. fiber surface treatment was carried out by using a polyethylene or silane coupling agent in combination with modifying agents. The glass fibers were embedded in the brittle epoxy matrix. Mode I, mode II, and mixed-mode I/II tests were performed in order to determine critical strain energy release rates. Double cantilever beam (DCB), end-notched flexure (ENF), and mixed-mode flexure (MMF) specimens were used. For both types of fiber surface treatment about the same values of mode I initiation fracture toughness G_IC ^init were obtained. It was observed that in mode I interlaminar crack growth in the DCB test for the composite sized by polyethylene, the crack propagation is accompanied by extensive fiber bridging. For both fiber surface treatments interlaminar fracture toughness increases considerably with increasing of crack length. For the fiber surface treatment with the silane coupling agent, the value of mode II initiation fracture toughness G_IIC ^init was about 2.5-times higher in comparison with that of a composite sized by polyethylene. For both types of fiber surface treatments the mixed-mode I/II test has shown a similar behavior to the mode I DCB test.     

Fracture toughness characteristics of laminated composites

Parameters characterizing the resistance of laminated composites to interlaminar fracture are discussed. The properties of the specific interlaminar fracture work, i.e., the amount of work spent on the formation of a unit of new surface of interlaminar crack, were examined. Taking account of the anisotropy of the material, this work may be characterized using a matrix. Upon change in the direction of crack growth, the matrix elements are transformed similarly to the components of a symmetrical second rank tensor. An interpretation is offered for the matrix elements. The proposed theoretical model was in accord with our experimental results.Moscow Engineering Institute. Technical University, Russia. Translated from Mekhanika Kompozitnykh Materialov, No. 1, pp. 21–31, January–February, 1996.     

Free-edge stress analysis of general composite laminates under extension, torsion and bending

In this study, based on the reduced form of elasticity displacement field for a long laminate, an analytical method is established to exactly obtain the interlaminar stresses near the free edges of generally laminated composite plates subjects to extension, torsion, and bending. The constant parameters being in the displacement field, which describe the global deformation of a laminate, are appropriately calculated by using the improved first-order shear deformation theory. Reddy’s layerwise theory is subsequently employed for analytical and numerical examinations of the boundary layer stresses within arbitrary laminated composite plates. Various numerical results are developed for the interlaminar normal and shear stresses along the interfaces and through the thickness of laminates near the free edges. Finally the effects of end conditions of laminates and geometric parameters on the boundary-layer stress are studied.     

Computation of the Three-Dimensional Stress State in Composite Shell Structures with Mixed Finite Elements

Being able to compute the complete three-dimensional stress state in layered composite shell structures is essential in order to examine complicated interlaminar failure modes such as delamination. We lay out a mixed finite element formulation with independent displacements, rotations, stress resultants and shell strains. A mixed hybrid shell element with 4 nodes and 5 or 6 nodal degrees of freedom is developed, so that the element formulation can also be used for problems with shell intersections. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Experimental-theoretical analysis of the state of stress and strain of stiffened glass-reinforced plastic shells subjected to local loads

The state of stress and strain of stiffened glass-reinforced plastic shells subjected to local loads applied to the stiffeners is investigated with allowance for interlaminar shear strains over the thickness of the glass-reinforced plastic. The solution is based on Vlasov's semimembrane theory. A comparative experimental-theoretical investigation has been made using models subjected to radial and axial loads and twisting moments.     

Effect of certain structural parameters of reinforced plastics on strength and stiffness in interlaminar shear

This study is concerned with the manner in which strength and stiffness in interlaminar shear of orthogonally reinforced laminates depend on the mechanical characteristics of their components, on the volume fraction of reinforcing material, on the coefficients of reinforcement distribution in the stacking direction, and on the microporosity of the polymer matrix, as well as other microdiscontinuities in the various structural components of the plastic. For the purpose of assessing the strength and the stiffness in interlaminar shear, with the volume fraction of reinforcing material and the interlaminar microporosity of the composite material varied, approximate formulas have been derived and their practical applicability subsequently confirmed in tests on glass-resin textolite specimens. It is shown that the indicated mechanical properties of a plastic are appreciably affected by these two structural parameters, especially so by its interlaminar microporosity.     

Numerical elastic-plastic simulation of interlaminar stresses in a notched angle-ply thermoplastic composite laminate

A thermoplastic angle-ply AS4/PEEK laminate with a hole is considered. The interlaminar stresses along the hole edge at different interfaces under uniaxial extension are investigated. According to the symmetries of the structure and loading, a suitable finite-element model is developed. Utilizing a three-dimensional elastic-plastic finite-element procedure elaborated previously, a finite-element modeling of the interlaminar stresses in a thick angle-ply composite laminate is carried out. Based on the interlaminar stresses obtained, the dangerous locations of delamination initiation are predicted. The results obtained indicate that there is some relationship between the dangerous locations and fiber orientations in the adjacent layers, and it maybe inferred that the critical locations are near the regions where the hole edge is tangent to the fiber orientation in the layers adjacent to the interface. The interlaminar stresses at the same interfaces are not sensible to distances from the midplane of the laminate. Very high interlaminar tensile stresses are found to exist on the hole edge at the +25°/+25° or –25°/–25° interfaces, and delaminations can initiate there first. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 45, No. 3, pp. 427-440, May-June, 2009.     

Asymptotic Analysis of Free-Edge and Free-Corner Effects in Laminate Structures

Stress fields in the vicinity of free edges and corners of composite laminates exhibit singular characteristics and may lead to premature interlaminar failure modes like delamination fracture. It is of practical interest to investigate the nature of the arising free-edge and free-corner stress singularities - i.e. the singularity orders and modes - closely. The present investigations are performed using the Boundary Finite Element Method (BFEM) which in essence is a fundamental-solution-less boundary element method employing standard finite element formulations. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Free-Edge Stresses around Holes in Laminates by the Boundary Finite-Element Method

In the present paper the boundary finite-element method is presented as a highly efficient technique for the numerical investigation of the free-edge stresses around a circular hole in laminates. In this method, as in the boundary element method, only the boundary needs to be discretized, whereas the element formulation in essence is finite-element based. The surface discretization provides a high numerical efficiency and requires less computation time compared to finite-element analyses. Numerical results for the concentration of interlaminar stresses at holes in composite laminates show a very good agreement with comparative finite-element calculations.