复合材料层合板分层裂纹扩展的试验及数值模拟研究

由于复合材料层合板在制作加工过程中会产生初始的微裂纹,且复合材料具有各向异性的力学特征,因此在承载条件下,这些微裂纹可能在黏结面上及在层与层之间发生、扩展、发展成宏观层间裂纹,最终导致结构的破坏.本文对复合材料层合板试验件进行了拉伸试验,提出了一种复合材料层合板分层建模方法,利用这种方法对一组试验件进行了单向拉伸荷载作用下极限承载力的研究,得出其分层扩展的位置及层间应力分布规律,将试验结果与数值模拟结果进行了对比,论证了分层建模方法对分析复合材料层合板分层裂纹扩展问题的可行性与适用性.     

复合材料层板分层伴以横向裂纹扩展的三维有限元分析

本文采用三维有限元模型,对典型铺设的[0_2/±45_2/90_2]_s碳/环氧复合材料层板中分层伴以横向裂纹的产生和扩展导致的层间应力分布进行了分析。计算结果表明层间裂纹首先在90°层中部出现并开裂至相邻界面处而产生横向裂纹,横向裂纹的出现引起局部分层按三角形状扩展;并指出分层损伤过程是一个主导性的稳定扩展过程,这是导致刚度下降的主要因素。最后,数值计算结果与实验结果比较,两者是吻合的。     

微裂纹损伤对复合材料层板分层的影响

本文将双标量损伤模型推广,考虑层内及层间微裂纹损伤相互作用,建立了层合复合材料的多标量损伤本构关系,并将该本构关系应用于含贯穿分层层板的后屈曲分析,采用准三维有限元法研究了分层尖端场与极限损伤区演变的特点,结果表明：损伤导致裂尖区层间应力奇异性消失,引起层间应力重新分布;损伤区的发展方向取决于分层上下表面层的铺设角,能量释放率不能反映损伤和分层扩展特征,因而有必要建立基于损伤演化的分层扩展准则。     

层合圆柱壳体表层局部分层屈曲及其扩展

应用能量法研究了对称铺层的闭合圆柱形层合壳体的表层局部分层屈曲问题；表层分层形式有三角形，椭圆和双纽线；分析了母层壳，子层壳的几何参数，物理参数同分层屈曲临界应变值之间的关系；给出了子层分层为椭圆的分层扩展的可能发生方向。     

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

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

界面裂纹非对称扩展模型

提出了一种新的复合材料结构中层间脱层裂纹扩展模型，该模型能够揭示裂纹扩展对裂纹左右尖端能量释放率的相对大小及界面材料性能的依赖性，能够比较真实地揭示复合材料结构层间脱层裂纹扩展复杂性的特点．     

基于互补理论和扩展有限元法的摩擦接触裂纹应力强度因子计算

从虚功方程出发，结合扩展有限元离散技术与接触条件的非线性互补表述，建立了摩擦接触裂纹问题的扩展有限元非线性互补模型，将不等式接触条件转化为非线性互补类的非光滑方程组，并采用基于广义导数的非光滑阻尼牛顿法求解方程组，无需引入任何额外人工变量以及迭代求解。以含中心倾斜裂纹平板和边裂纹平板为例，运用相互作用积分法计算摩擦接触裂纹的应力强度因子，将其结果与理论解进行对比分析，该方法都能给出精确的计算结果；基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟，计算结果表明，受压裂纹数值结果与实验结果比较吻合，从而验证了本文方法的有效性与正确性。     

复合材料铺层拼接层间断裂韧性的试验研究

针对评价复合材料层合板层间断裂韧性的测量,提出了用拉伸试验法测定Ⅱ型层间断裂韧性,设计了内含铺层拼接区的分层破坏试验层合板,制备了拉伸试件.通过拉伸试验测得了拼接区开裂和分层裂纹稳态扩展过程中的载荷与变形规律:层间破坏具有Ⅱ型断裂特征,且裂纹扩展比较稳定.利用测试数据计算出断裂功,并以临界能量释放率表示层合板的Ⅱ型层间断裂韧性,结果表明用铺层拼接件拉伸法进行层间断裂韧性试验是可行的.     

界面超高速自相似动态分层分析：反平面情况

对材料界面超高速自相似动态分层的反平面问题进行了解析分析。分层模拟为界面裂纹由零长度自相似扩展，扩展速度为蹭音速或超音速。首先考虑运动集中载荷作用下界面动态分层的情况，利用界面裂纹自相似扩展的运动位错模型将问题归结为奇异积分方程，并求得解析解，分析了裂纹尖端的应力奇性，获得了动应力强度因子。最后，利用叠加原理给出了ｘ＾ｎ型载荷作用下界面动态分层的解。     

含分层损伤的变角度纤维层合板后屈曲力学行为研究

变角度（Variable angle tow,简称VAT）纤维复合材料层合板的纤维方向能够连续变化．相较于传统的直线纤维层合板,此类层合板通过刚度变化,整体的屈曲性能可以得到很大的提升．本文利用ABAQUS 自带的粘结单元(Cohesive Element)对预制圆形分层的变角度纤维复合材料层合板进行了后屈曲力学行为研究,得出载荷位移曲线,以及分层裂纹萌生和扩展的情况．然后本文分析了预制分层尺寸对板的刚度、前后屈曲阶段和裂纹萌生及扩展的影响．最后通过变角度纤维层合板和直线纤维层合板的后屈曲力学行为进行对比,深入探索了变角度复合材料层合板在抵抗分层裂纹萌生和扩展方面的优势．     

Study on failure process of delaminated stiffened composite plates under compression

Failure behavior of the delaminated stiffened composite plates under compression is studied by the finite element method, based on a Global-Local variational model. A virtual crack closure technique and a self-adaptive grid moving scheme are proposed to predict the delamination growth process. The contact effect along the delamination front is considered. The numerical results show that the influences of the distribution and location of the stiffeners, the configuration and size of the delamination, the boundary condition and the contact upon the failure behavior of the plates are significant. The project supported by the National Natural Science Foundation of China (59975014)     

A thermomechanical cohesive zone model for bridged delamination cracks

The coupled thermomechanical numerical analysis of composite laminates with bridged delamination cracks loaded by a temperature gradient is described. The numerical approach presented is based on the framework of a cohesive zone model. A traction-separation law is presented which accounts for breakdown of the micromechanisms responsible for load transfer across bridged delamination cracks. The load transfer behavior is coupled to heat conduction across the bridged delamination crack. The coupled crack-bridging model is implemented into a finite element framework as a thermomechanical cohesive zone model (CZM). The fundamental response of the thermomechanical CZM is described. Subsequently, bridged delamination cracks of fixed lengths are studied. Values of the crack tip energy release rate and of the crack heat flux are computed to characterize the loading of the structure. Specimen geometries are considered that lead to crack opening through bending deformation and buckling delamination. The influence of critical mechanical and thermal parameters of the bridging zone on the thermomechanical delamination behavior is discussed. Bridging fibers not only contribute to crack conductance, but by keeping the crack opening small they allow heat flux across the delamination crack to be sustained longer, and thereby contribute to reduced levels of thermal stresses. The micro-mechanism based cohesive zone model allows the assessment of the effectiveness of the individual mechanisms contributing to the thermomechanical crack bridging embedded into the structural analysis.     

FRP strengthened walls with delaminated regions – Geometrically nonlinear response

The paper studies the geometrically nonlinear behavior of walls that are strengthened with fiber reinforced polymer (FRP) composite materials but include pre-existing delaminated regions. The paper uses an analytical–numerical methodology. Three specially tailored finite elements that correspond to perfectly bonded regions, to delaminated regions where the debonded layers are in contact, and to delaminated regions where the debonded layers are not in contact are presented. All finite elements are based on a high order multi layered plate theory. The geometrical nonlinearity is introduced by means of the Von Karman nonlinear strains whereas the contact nonlinearity is handled iteratively. The validity and convergence of the finite element models is demonstrated for each type of element through comparison with closed form analytical solutions available for specific cases. The unified model that combines the three types of finite element is then used for studying the nonlinear behavior of a locally delaminated FRP strengthened wall under in-plane normal and in-plane shear loads. Finally, conclusions regarding the effect of the delamination on the response of the strengthening system, on the conditions that evolve in the bonded region that surrounds the delamination, and on the global response of the multi-layered structure are drawn. Additional conclusions regarding the application of the modeling approach to other delamination sensitive layered structural systems close the paper.     

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.     

具有分层损伤的不同加筋形式复合材料层合板的后屈曲性态研究

基于板的一阶剪切理论和V on-K arm an大挠度理论,分别推导了复合材料层合板和层合梁的几何非线性有限元列式,提出了含嵌入分层的复合材料加筋层合板在受压缩载荷作用下的后屈曲有限元分析方法,对在板厚方向具有不同分层位置的加筋板结构进行了有限元数值分析,研究了不同的加筋方式及筋的分布对具有分层损伤的复合材料加筋层合板的后屈曲性态的影响,所得结果对确定在压缩载荷作用下含损伤复合材料加筋层合板的剩余承载能力具有参考价值。     

Element-failure concepts for dynamic fracture and delamination in low-velocity impact of composites

An element-failure algorithm is proposed and incorporated into a finite element code for simulating dynamic crack propagation and impact damage in laminated composite materials. In this algorithm, when a crack is propagating within a finite element, the element is deemed to have partially failed, but not removed from the computations. Consequently, only a fraction of the stresses that were computed before the crack tip entered the element contribute to the nodal forces of the element. When the crack has propagated through the element, the element is completely failed and therefore can only resist volumetric compression. This treatment of crack propagation in isotropic solids allows fracture paths within individual elements and is able to accommodate crack growth in any arbitrary direction without the need for remeshing. However, this concept is especially powerful when extended to the modeling of damage and delamination in fibre-reinforced composite laminates. This is because the nature of damage in composite laminates is generally diffused, characterized by multiple matrix cracks, fibre pullout, fibre breakage and delaminations. It is usually not possible to define or identify crack tips in the tradition of fracture mechanics. Since parts of a damaged composite structure are often able to partially transmit load despite the presence of some damage, it is advantageous to model the damaged portions with partially failed elements. The damage may be efficiently modeled and tracked using element-failure concepts, with the application of appropriate failure criteria and damage evolution laws. The idea is to embody the effects of damage into the effective nodal forces of the finite element. In this paper, we report the novel use of element-failure concepts in the analysis of low-velocity impact damage of composite laminates. The initiation and propagation of delaminations arising from the impact are predicted and the results show qualitative agreement with experimental observation of the formation of multiple delaminations in impact-damaged specimens. While such delaminations do not permit transmission of tensile stress waves across the cracked surfaces, transmission of compressive stress waves are allowed in the simulation. It is further shown that, when elements are allowed to fail, the dynamic stress wave distributions are altered significantly. In the element-failure algorithm, the issue of interpenetration of delamination surfaces in the model does not arise. This is a significant advantage over the conventional method of explicitly modeling the delamination surfaces and crack front, where generally, much computational time must be spent in employing contact algorithms to ensure physically admissible solutions. Finally, we also demonstrate the simulation of crack propagation of pre-notched specimens of an isotropic material under initial conditions of mode II loading using the element-failure algorithm. The numerical results showed that the cracks propagated at an angle of about 70° with respect to the notches, in agreement with the experimental results of Kalthoff.     

Effects of voids on postbuckling delamination growth in unidirectional composites

This work examines the effects of manufacturing induced voids on the postbuckling behavior of delaminated unidirectional composites. In the finite element model developed, a through-width delamination is introduced close to one surface of a flat panel, and a void is placed in the delamination plane ahead of each delamination front. The panel is subjected to compression in the fiber direction. The postbuckling delamination growth is studied by calculating the strain energy release rate (SERR) using the virtual crack closure technique. Local stress analyses of the region near the delamination front are also performed to further investigate the void effects. It is found that although the presence of void does not significantly alter the postbuckling transverse displacement of the delaminated panel, the induced stress perturbation by the void affects the SERR. The Mode II SERR as well as the total SERR increase depending on the size of the void and its distance from the delamination front. Since the Mode I SERR shows non-monotonic behavior with the applied load, the effects of voids are studied on its maximum value.     

Mode III fracture behavior of laminated composite with edge crack in torsion

A three-dimensional (3-D) finite element analysis was performed on a [90,(+45/−45)_n,(−45/+45)_n,90]_s class of laminated composites under the edge crack torsion (ECT) test configuration. Finite element delamination models were established and formulas for calculating the Mode III fracture toughness from 3-D finite element models were developed. The relations between the interlaminar fracture behavior and various configuration parameters were investigated and the effects of point loads, ends, geometry, Mode II interference, and friction were evaluated. Results showed that with proper selection of ECT specimen configuration and layup, the delamination could grow in pure Mode III in the middle region of the specimen. Specimen end effect played an important role in the ECT test. A Mode II component occurred in the end regions but it did not interfere significantly with the Mode III delamination state. Specimen dimension ratio, layup, and crack length exhibited significant effect on the interlaminar fracture behavior and the calculated strain energy release rates. However, friction between crackfaces was found to have negligible effect on the interlaminar properties.