高温空气下C/SiC复合材料断裂韧性实时测试和微观结构表征分析

为了研究高温空气下C/SiC复合材料断裂韧性和微观结构,采用单边切口梁三点弯曲法实时测试了C/SiC复合材料在高温空气下的断裂韧性,并采用电子扫描显微镜 (scanning electron microscope,SEM)和X 射线衍射分析仪 (X-ray diffraction, XRD)分析了复合材料在不同温度下的破坏断口和失效机制。研究结果表明随测试温度升高,C/SiC复合材料断裂韧性降低,材料的断裂形式由脆性断裂逐渐演变成塑性断裂。从室温升温到1 000 ℃测试温度条件下,C/SiC复合材料的断裂韧性由12.5 MPa·m1/2降低为10.96 MPa·m1/2,降幅仅为12%,C/SiC复合材料高温断裂韧性良好。不同温度下,材料呈现出不同形式的断裂形貌。常温下断口形貌主要可以看到纤维拔出的现象,随着温度的升高,该现象基本消失,断裂截面变得更平整,材料的强度主要取决于基体的强度。     

STUDY ON NONLINEAR CONSTITUTIVE RELATIONSHIP AND FRACTURE BEHAVIOR OF STITCHED C/SiC COMPOSITES$^{\bf 1)}$

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.     

三维编织CMC的弯曲断裂研究

从界面断裂的角度出发，对三维编织CMC的断裂作了理论研究和数值分析，对于三点弯曲试件，通过数值拟合修正了能量释放率G的理论表达式中的自由常数A，同时也研究了材料的各个参变量对于断裂韧性的影响，由此得出了一个基本完善的三点弯曲试件断裂韧性G的理论公式，该能量释放率方法可以应用于单试件的试验计算，与断裂韧性的柔度标定方法相比，该方法一方面可以减少试验件数量；另一方面，试验结果显示出在试件切口尺雨处于0．4≤a/W≤0.5时，可以获得比较稳定的断裂韧性值。     

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

Microvoid multistage nucleation model and its application in analysis of micro damage and fracture

In this paper, the microvoid multistage nucleation model [14,15] suggested by the authors of this paper has been studied on the micro ductile damage and fracture of metallic material under large elastic-plastic deformation.Using this model, the analyses of micro damage and fracture for various axisymmetric tensile specimens and for TPB and CCP cracked specimens have been carried out. And the results from these analyses on damage development and fracture are in good agreement with the experimental ones for axisymmetric specimens and reasonable for cracked specimens from the microscopic point of view.The project surported by National Science Foundations of China.     

Computer simulated analyses on deformation and fracture of non-cracked and cracked specimens

A unified damage and fracture model, the combinatory work density model, which is suitable for either non-cracked body or cracked body has been suggested^[t−7]. In the present paper, the deformation and fracture of the two kinds of tensile spceimen and TPB specimen made of 40Cr steel have been simulated by using the new model together with the large elastic-plastic deformation finite element method. The results give a good picture of the whole deformation and fracture processes of the specimens in experiments; especially, the results on the TPB specimen can be used to obtain the relationship between load and displacement at the loading pointP-Δ, and between crack extension and displacement at the loading point Δa-Δ, the resistance curveJ _R -Δa and the fracture toughnessJ _1C . All the results are in remarkable agreement with those obtained by experiments. Therefore the model suggested here can be used to simulate crack initiation and propagation in non-cracked body and fracture initiation and crack stable propagation in cracked body. The project supported by National Natural Science Foundation of China     

Toughening of ferroelectrics by the out-of-plane poling

Subjected to the prior out-of-plane poling, the ferroelectrics can be toughened considerably. The present paper describes the variation of the stress intensity factor (SIF) by 90° switching in ferroelectrics. The analysis is carried out for the combined mechanical and electrical loading, with simple relations obtained for the case of the purely electrical loading. The out-of-plane poling is found to raise the SIF for the crack initiation, but appreciably reduces the SIF for the crack growth in a steady state. More stable fracture resistance curves can be achieved by the out-of-plane poling. This prediction is supported quantitatively by the testing data of SENB specimens of PZT-5 samples, when the toughening effects of polings in three orthogonal directions are compared. The project supported by Sino-Germany Center for Science and Technology (G2050/8)     

Crack-tip damage and fracture toughness of borsic/titanium composite

The fracture behavior of center-notched unidirectional borsic/titanium composites under uniaxial tensile loading has been investigated. Load-crack-opening displacement (COD) curves were obtained for crack length-to-width ratios ranging from 0.05 to 0.5. The COD was measured across the crack surfaces by means of laser interferometric displacement gage (gage length of 500 μm). As expected, the interferometric displacement gage (IDG) was very sensitive to the appearance of damage at the crack tip, yielding nonlinear load-COD curves. The crack-tip damage in borsic/titanium was found to grow coplanarly with the original crack and consist of fiber breaks and matrix plastic deformation, the former preceding the latter. The failure mode and the crack-tip-damage growth were found to be different from those observed in unidirectional boron/aluminum and these differences are discussed. The examination techniques include radiography, SEM, photomicrographs, interference microscope and optical-interferometry technique. It was found that the fracture strengths and crack-tip-damage growth can be analyzed by the analytical models proposed in the literature. The resistance-curve method was also tried and found to be applicable. A correlation has been established among the fracture strength, crack-opening displacement and failure mode.     

不同撞击速度下穿燃弹侵彻陶瓷/铝合金复合靶板时弹芯破碎失效特性研究

为了研究12.7 mm穿燃弹以不同速度撞击陶瓷/铝合金复合靶板时弹芯的破碎及失效特性,开展了12.7 mm穿燃弹以434.5～844.6 m/s速度撞击SiC陶瓷/6061T6铝合金复合靶板的弹道试验,分析了弹靶的失效模式。弹芯在侵彻靶板后会产生不同尺寸的碎片,使用回收箱收集弹芯碎片并用不同孔径筛网对其进行筛分、称重,得到了不同撞击速度下弹芯碎片的质量分布,并对不同部位的弹芯碎片断口形貌进行了宏观和微观观测分析。研究结果表明:背板失效模式为碟形变形-剪切穿孔-花瓣形失效,试验后的弹芯碎片累积质量分布符合Rosin-Rammler幂率分布规律,且随着着靶速度的增大,小质量碎片质量增加;弹芯在冲击过程中等效直径较大碎片（大于8 mm）失效模式为拉伸脆性断裂,而等效直径小于2 mm的碎片上存在局部塑性剪切断裂。     

纤维─基体界面剪切强度与断裂韧性的表征和测定

根据单向复合材料纤维分布、剪切变形与裂纹扩展的细观分析,提出纤维一基体界面剪切强度和断裂韧性的简易细观表征与测定方法．     

纤维─基体界面剪切强度与断裂韧性的表征和测定

根据单向复合材料纤维分布、剪切变形与裂纹扩展的细观分析,提出纤维一基体界面剪切强度和断裂韧性的简易细观表征与测定方法．     

Brittle fracture: Variation of fracture toughness with constraint and crack curving under mode I conditions

The effect of constraint on brittle fracture of solids under predominantly elastic deformation and mode I loading conditions is studied. Using different cracked specimen geometry, the variation of constraint is achieved in this work. Fracture tests of polymethyl methacrylate were performed using single edge notch, compact tension and double cantilever beam specimens to cover a bread range of constraint. The test data demonstrate that the apparent fracture toughness of the material varies with the specimen geometry or the constraint level. Theory is developed using the critical stress (strain) as the fracture criterion to show that this variation can be interpreted using the critical stress intensity factorK _Cand a second parameterT orA ₃,whereT andA ₃are the amplitudes of the second and the third term in the Williams series solution, respectively. The implication of this constraint effect to the ASTM fracture toughness value, crack tip opening displacement fracture criterion and energy release rateG _Cis discussed. Using the same critical stress (strain) as the fracture criterion, the theory further predicts crack curving or instability under mode I loading conditions. Experimental data are presented and compared with the theory.     

Effect of specimen size and crack depth on 3D crack-front constraint for SENB specimens

Three-dimensional (3D) elastic–plastic finite element analyses (FEA) are performed to study constraint effect on the crack-front stress fields for single-edge notched bend (SENB) specimens. Both rectangular and square cross-section of the specimens with a deep crack of a/W=0.5 are considered to investigate the effect of specimen size. A square-cross-section specimen with a shallow crack of a/W=0.15 is also considered to examine the effect of crack depth. Stresses from FEA at the crack front on different planes of the specimen are compared with those determined by the J–A₂ three-term solution. Results show that in-plane stress fields can be characterized by the three-term solution throughout the thickness even in the region near the free surface. Cleavage fracture toughness data is compared to predict the effects of specimen size and crack depth on fracture behavior. It is found that the distributions of crack opening stress are nearly the same for the SENB specimens at the critical J which is consistent with the RKR model. Furthermore our results indicate that there is a distinct relationship between the crack-front constraint and the cleavage fracture toughness. By introducing the failure curves, the minimum fracture toughness and scatter band can be well captured using the J–A₂ approach.     

Experimental investigation of dynamic punch tests on isotropic and composite materials

An experimental investigation is conducted on the two-dimensional punch problem for isotropic materials and unidirectional fiber-reinforced composite materials under quasi-static and impact loading. Singular stresses are generated in the specimen near the punch corners, and the stress intensity factorK _Iis introduced to describe the singular stress field. Laser interferometry is used to measure in-plane stresses (transmission mode) and out-of-plane displacements (reflection mode) and then estimate the stress intensity factor. In the dynamic case, a high-speed photography technique was employed to capture the transient response of the specimen and measureK(t) just after the impact. In all the cases, a good agreement between the measurements ofK and theoretical predictions was found.     

RHEOLOGICAL FRACTURE ANALYSIS OF CABLE BY CONSIDERING THE SECOND ORDER EFFECTS OF RANDOM SLACK

In this paper,we consider the cable as a bundle consisting of n sub-bundles,with mparallel tension members per sub-bundle,and the tension members themselves are polymeric yarns im-pregnated with a resin matrix.The nonfailed members at any instant must share an applied system loadaccording to some rule,since there is a clearly expressed dependence of the fracture on the durationand character of the loading.So then,the fracture of cable is a process of nonlinear dynamic evolu-tion,which accommodates to the non-equilibrium thermodynamics of irreversible processes by itself.Let us assume that the polymeric yarns are as viscoelastic solid,under certain probabilistic assump-tions,according to the principles of theology of bodies with defects,the relationship between the singlemember loading and failure and the bundle loading are investigated.It can be shown that the bundlefailure time is asymptotically normally distributed as the number of members grows large.After astudy of the second order effects of random slack,it is known that the asymptotic mean and varianceare functions of the parameters of loading and single member rhcological behavior.Hence the loss inthe asymptotic bundle strength mean brought about by random member and sub-bundle slack,L.,andthe loss in the asymptotic bundle strength variance caused by random member slack,Δ_(?),are deter-mined.And finally,it is known that the asymptotic time of failure can make up a considerable part ofthe fracture of cable,and the fracture of cable is a time-dependent process of rheological fracture.     

A study of the in situ transverse tensile lamina strength and the classical solution of the opening mode crack in a 0/90/0 laminate

In this paper, the classical solution of the opening mode crack in the 90° layer of 0/90/0 laminates has been determined by means of Fourier transformations and the procedure of Copson for a pair of dual integral equations. The fracture behavior and the in situ transverse strength of the 90° layer have been quantitatively studied in graphite/epoxy laminates, based on the solution obtained above. The results show that the stress intensity factor of this kind of laminates, which is different from that of a single unidirectional 90° layer, decreases with the increase in thicknessb, or modulusE _L orG _LT of the 0° layer and also decreases with the decrease in the thickness of the 90° layer. So the lamination effect manifests itself and thein situ transverse strength of the 90° layer is thereby enchanced. The theoretical calculations agree with the experimental data presented by D.L. Flaggs. Supported by National Natural Science Fundation of China.     

金属试件变形断裂过程的计算模拟分析——组合功密度模型的应用

本文作者已建立了一个既适用于裂纹体亦适用于无裂纹体的统一损伤断裂模型——组合功密度模型。在本文中,作者应用该模型并结合大应变弹塑性有限元方法对40Cr制的圆棒光滑拉伸试件、圆棒切口拉伸试件和三点弯曲裂纹试件的变形和断裂过程进行了计算模拟。模拟结果很好地再现了两种圆棒拉伸试件的实验过程;而对于三点弯曲试件,模拟得到的载荷——加载位移关系曲线、裂纹扩展量——加载位移关系曲线、J_R阻力曲线和断裂韧性J_(IC)值等均与实测结果相当吻合。证实了该模型既能用于模拟预测无裂纹体中的裂纹萌生和扩展,也能模拟预测裂纹体中的裂纹起裂和稳定扩展过程。     

Stress–strain and fracture behaviour of 0°/90° and plain weave ceramic matrix composites from tow multi-axial properties

A computationally economic finite-element-based multi-linear elastic orthotropic materials approach has been developed to predict the stress–strain and fracture behaviour of ceramic matrix composites with strain-induced damage. The finite element analysis utilises a solid element to represent a homogenised orthotropic medium of a heterogeneous uni-directional tow. The non-linear multi-axial stress–strain behaviour has been discretised to multi-linear elastic curves, which have been implemented by a user defined subroutine or UMAT in the commercial finite element package, ABAQUS. The model has been used to study the performance of two CMC composites: a SiC (Nicalon) fibre/calcium aluminosilicate (CAS) matrix 0°/90° cross-ply laminate Nicalon/CAS; and, a carbon fibre/carbon matrix–SiC matrix (C/C–SiC) plain weave laminate DLR-XT. The global stress–strain curves with catastrophic fracture behaviour and effects of fibre waviness have been predicted. Comparisons have been made between the predictions and experimental data for both materials. The predicted results when fibre waviness is taken into account compare well with the experimental data.     

Fracture analysis of U-notched disc-type graphite specimens under mixed mode loading

Fracture phenomenon was investigated both experimentally and theoretically for a type of coarse-grained polycrystalline graphite weakened by a U-shaped notch under mixed mode loading. First, 36 disc-type graphite specimens containing a central U-notch, so called in literature as the U-notched Brazilian disc (UNBD), were prepared for four different notch tip radii and the fracture tests were performed under mode I and mixed mode I/II loading conditions. Then, the experimentally obtained fracture loads and the fracture initiation angles were predicted by using the U-notched maximum tangential stress (UMTS) and the newly formulated U-notched mean stress (UMS) fracture criteria. Both the criteria were developed in the form of the fracture curves and the curves of fracture initiation angle, in terms of the notch stress intensity factors (NSIFs). The results showed that while the criteria could predict successfully the experimental notch fracture toughness values, the UMS criterion provides slightly better predictions than the UMTS criterion, particularly for shear-dominant deformations. Also, found in this research was that the curves of fracture initiation angle were almost identical for the two criteria which both could predict well the experimental results.     

半结晶聚合物损伤演化的实验表征与数值模拟

损伤本构模型对研究材料的断裂失效行为有重要意义, 但聚合物材料损伤演化的定量表征实验研究相对匮乏. 通过4种高密度聚乙烯(high density polythylene, HDPE)缺口圆棒试样的单轴拉伸实验获得了各类试样的载荷-位移曲线和真应力-应变曲线, 采用实验和有限元模拟相结合的方法确定了HDPE材料不同应力状态下的本构关系, 并建立了缺口半径与应力三轴度之间的关系;采用两阶段实验法定量描述了4种HDPE试样单轴拉伸过程中的弹性模量变化, 并建立了基于弹性模量衰减的损伤演化方程, 结合中断实验和扫描电子显微镜分析了应力状态对HDPE材料微观结构演化的影响. 结果表明缺口半径越小, 应力三轴度越大, 损伤起始越早、演化越快; 微观表现为: 高应力三轴度促进孔洞的萌生和发展, 但抑制纤维状结构的产生;基于实验和有限元模拟获得的断裂应变、应力三轴度、损伤演化方程等信息提出了一种适用于聚合物的损伤模型参数确定方法, 最后将本文获得的本构关系和损伤模型用于HDPE平板的冲压成形模拟, 模拟结果与实验结果吻合良好.