短切纤维增韧泡沫夹芯复合材料梁界面断裂过程的物质点方法模拟

对短切纤维增韧泡沫夹芯复合材料梁界面韧性试验结果进行了分析讨论并基于物质点方法(material pointmethod, MPM)对试验过程进行了数值模拟。在MPM方法中,通过可视准则引入不连续性来处理裂纹问题,发展了包含裂纹的MPM 算法,模拟了泡沫夹芯复合材料梁界面断裂试验过程,数值分析结果和试验数据取得了良好的一致性。研究结果表明短切纤维增韧工艺能够显著提高泡沫夹芯复合材料结构的界面韧性和承载能力,同时表明该文推导的包含裂纹的MPM方法处理断裂问题的精确性和有效性。     

温度对缝合气凝胶夹芯复合材料面内压缩性能的影响

分别开展缝合气凝胶夹芯复合材料在不同温度下的面内压缩试验,研究材料在室温、300℃、600℃和800℃下的面内压缩力学性能,并采用微焦点工业CT扫描的方法对试样内部结构进行分析,结合有限元分析方法,探究其结构破坏机制。结果表明:在面内压缩载荷作用下,材料存在极限载荷,面板的局部屈曲、芯层的剪切破坏以及缝线柱的断裂是材料破坏的主要方式。随着温度的升高,材料的面内压缩模量和极限载荷也逐渐升高,面板破坏处的断口逐渐呈现出类似脆性的断裂。300℃、600℃和800℃下材料的面内压缩模量分别为室温的1.05倍、1.57倍和1.65倍;极限载荷分别为室温的1.14倍、1.46倍和1.67倍。室温下有限元分析结果和试验结果的对比,验证了缝合气凝胶夹芯复合材料面内压缩破坏模式的合理性。     

结构参数对树脂基纤维编织复合材料折叠夹芯结构力学性能的影响

折叠夹芯结构是一种新型的复合材料夹芯结构,其结构参数对力学性能有重要的影响。文中以碳纤维和Kevlar平纹编织预浸料为芯材原料,采用热压工艺,制备了复合材料折叠夹芯结构试样。通过压缩试验得到不同条件下折叠夹芯结构在静态压缩载荷作用下的力-位移变化曲线。构建了复合材料折叠夹芯结构有限元模型,对不同结构参数复合材料折叠夹芯的力学性能进行了数值模拟分析,并将模拟结果与实验结果进行对比验证了模型的可靠性。实验及数值模拟的分析结果表明,随着芯材厚度的增加,折叠夹芯层的压缩强度呈线性增加,其破坏形式由假塑性变形逐渐向脆性破坏转化;面板对夹芯层的约束作用能够极大地提高压缩模量和强度,而且上下面板对压缩性能曲线有着不同的影响;折叠夹芯单元的高度、长度、折叠夹角等参数对其力学性能具有不同程度的影响。     

纤维复合材料的损伤韧度设计

文中以裂纹连续扩展的断裂力学方法, 研究单向纤维复合材料的细观损伤演化与韧度计算。用基于有限元的数值方法模拟复合材料的多模式破坏过程, 并计算复合材料在裂纹不同扩展路径下的有效韧度。研究表明界面的两个重要断裂参数(界面断裂能和混合度) 是控制复合材料破坏模式与增韧机理的关键参数。提出了复合材料韧度设计的基本原理。      

多层抗爆结构冲击响应无网格MPM法分析

在两层钢板中间夹衬泡沫铝等多孔材料构成多层复合结构被应用于抗爆、抗冲击的结构设计中,能够有效地降低冲击载荷对结构的破坏作用。为了研究多层复合结构的抗爆机理和变形破坏过程,使用材料非线性本构模型和无网格物质点法对在高速冲击载荷作用下各层材料的弹塑性大变形进行数值模拟。MPM法利用了欧拉法和拉格朗日法两者的优点,不仅与网格无关,也避免了有限元法中网格畸变,而且在对涉及多物质分界面的问题计算时,因MPM法的耦合条件自动满足,不需要考虑材料界面的变形和破坏,为计算多层抗爆结构的冲击响应建立了一个有效的无网格法数值模拟平台。     

缝合泡沫夹芯复合材料低速冲击的多尺度数值方法

为了发展缝合泡沫夹芯复合材料低速冲击损伤的多尺度分析方法, 建立了缝合泡沫简化力学模型, 将缝合泡沫等效为缝线树脂柱增强的正交各向异性芯材, 其材料参数由各组分性能及所占体积分数根据均一化理论计算得出; 同时, 建立冲击试验有限元模型, 通过界面元模拟面板与芯材之间的层间分层。采用GENOA渐进损伤分析模块对缝合结构冲击动态响应过程进行数值模拟, 并将计算结果与试验记录进行对比分析。结果表明: 缝合可以减小面板破坏面积, 抑制面板与泡沫分层的扩展; 但缝纫会对结构造成初始损伤, 较高的缝合密度使芯材刚度增加, 不利于泡沫结构的缓冲吸能。数值模拟结果与试验记录吻合良好, 验证了多尺度分析方法的正确性。     

等腰梯形蜂窝芯玻璃钢夹芯板的面内压缩性能

为研究等腰梯形蜂窝芯玻璃钢夹芯板面内压缩破坏机制, 利用材料试验机对夹芯板面内压缩性能进行了试验测试, 并开展了模拟研究。结果表明: 夹芯板的面内压缩破坏方式主要有面板折断、夹芯板屈曲失稳和夹芯板中面板与蜂窝芯脱粘3种类型。面板为夹芯板面内压缩的主要承载构件, 蜂窝芯对面板起到固支作用。面板结构参数与材料参数为影响夹芯板面内压缩抗压强度与抗压刚度主要因素, 多数蜂窝芯的结构参数与材料参数对夹芯板面内压缩抗压强度的影响微弱, 而个别蜂窝芯的结构参数对夹芯板面内压缩抗压刚度的影响比较显著。夹芯板体积一定时, 随着蜂窝芯胞体单元数量的增加, 夹芯板面内压缩的抗压强度与抗压刚度逐渐增大。      

复合材料泡沫夹芯板局部连接拉脱破坏试验与数值仿真

对复合材料泡沫夹芯板局部连接拉脱破坏进行了试验研究,分析了接头的破坏模式、失效载荷和面板对接头的影响。采用ABAQUS有限元软件进行了数值分析,通过与实验结果对比验证其模型的可靠性,预测分析内部的破坏模式以及结构参数对接头破坏的影响,研究了泡沫芯体内部的渐进破坏以及面板和泡沫芯体之间的胶层脱粘破坏。结果表明:泡沫夹芯板预埋螺栓连接结构灌封胶边缘的泡沫先产生裂纹后向中间扩展,中间区域全部开裂时两端裂纹沿着45°方向向上扩展。胶层开裂的区域呈弧形条状,分布在螺栓紧固件的两侧,在面板宽度方向,开裂的区域贯穿两侧。随着预埋件深度的增加最大破坏载荷也在增加,随着预埋件直径的增加亚临界破坏载荷和最大破坏载荷没有比较明显的变化,但最大破坏位移在减小。      

粘接界面泡沫铝夹芯板的三点弯曲失效数值模拟

对粘接界面泡沫铝夹芯板三点弯曲载荷下的变形特性进行了实验和数值模拟方面的研究。基于有限元软件ABAQUS建立了泡沫铝夹芯板的三维有限元模型,应用内聚力模型对三点弯曲过程中典型的破坏模式——面板与芯层的界面脱粘给予了合理的模拟,模拟所得的结果与实验结果比较吻合。并在此基础上分析了面板和芯层厚度对夹芯板承载能力和吸收能量能力的影响。结果表明,增加芯层的厚度能够更大程度上提高泡沫铝夹芯板的承载能力和吸收能量的能力。     

复合材料夹芯梁屈曲破坏模式及极限承载

为研究复合材料夹芯梁在轴压作用下的屈曲、后屈曲特性及承载能力,进行了试验研究与有限元仿真。首先,开展了系列复合材料夹芯梁屈曲特性试验,研究了铺层比例、梁长度、表层厚度及芯层厚度等因素对其屈曲、后屈曲破坏模式及极限承载的影响;然后,基于非线性屈曲理论,采用三维内聚力界面单元模拟面芯脱粘,并引入初始预变形及材料损伤准则对复合材料夹芯梁在轴压下的屈曲特性及极限承载进行仿真研究。结果显示:界面脱粘是屈曲破坏的重要模式;仿真计算的极限承载与试验结果相比,误差控制在10%以内。所得结论表明该方法可有效预报复合材料夹芯梁的后屈曲路径、破坏模式及极限承载。      

Debonding and crack kinking in foam core sandwich beams—II. Experimental investigation

Debonding and crack kinking in sandwich beams was experimentally examined, and also analyzed using the finite element method. Double cantilever beam (DCB) and shear fracture specimens employing aluminum facings bonded to a wide range of PVC and PMI foam cores using two types of adhesives were considered. It was found that the Young modulus of the core has a profound effect on the tendency of the facing/core interfacial crack to deflect (kink) into the core in DCB testing. In shear testing, crack kinking occurred for all core materials considered. The type of adhesive strongly influences the debond fracture resistance, but not the kink resistance and kink angle. The critical load for onset of kinking increased with increased core density. Finite element analysis of the fracture specimens enabled determination of mixed mode interfacial fracture toughness for the specimens that failed by debonding. For specimens that failed by kinking, interfacial stress intensity factors at the onset of kinking were determined. Measured kink angles compared favorably with kink angles calculated based on the interfacial stress intensity factors prior to kinking.     

梯度复合材料裂纹扩展路径和起裂载荷的有限元分析

为了模拟功能梯度材料(FGM)在工程应用中可能会出现的断裂问题并计算相应的开裂载荷,通过编写用户自定义UEL子程序将梯度扩展单元嵌入到ABAQUS软件中模拟功能梯度材料的物理场,并编写交互能量积分后处理子程序计算裂纹尖端的混合模式应力强度因子(SIF),采用最大周向应力准则编写子程序计算裂纹的偏转角,并模拟了裂纹扩展路径,计算了裂纹的起裂载荷。讨论了材料梯度参数对裂纹扩展路径以及起裂载荷的影响规律。通过与均匀材料的对比,验证了功能梯度材料断裂性能的优越性。研究表明:外载平行于梯度方向时,垂直梯度方向的初始裂纹朝着等效弹性模量小的方向扩展,且偏转角在梯度指数线性时出现峰值,并随着组分弹性模量比的增加而变大;当外载和初始裂纹均平行于梯度方向时,材料等效弹性模量和断裂韧性的增加或者梯度指数的减小都导致起裂载荷变大。     

Crack kinking of a delamination at an inclined core junction interface in a sandwich beam

The paper considers a general interface delamination and crack kinking from an inclined core junction in a sandwich beam. This particular problem is relevant for a newly developed peel-stopper component for sandwich structures.A finite element model (FEM) was developed and calibrated against a known model by He and Hutchinson. The numerically and analytically determined solution coefficients were in a perfect agreement with each other, so the necessary generalisation of results can be obtained through the application of FEM-analyses. The FE-model was used to determine solution coefficients for a number of interface compositions of practical interest. As expected some of the coefficients were quite sensitive to the specific material combination, which confirms that accurate solution strategies are important.The solution coefficients obtained were further applied to the analysis of the crack propagation and kinking process in three different sandwich beam configurations, each of which contained an inclined junction of 20°, 30°, or 40°. The objective was to examine how the core junction angle and the fracture mechanical properties of the sandwich components influenced the crack kinking tendency. The latter is vital for the design and functionality of a newly developed peel-stopper. It was shown that smaller core junction angles will lead to longer crack propagation (delamination) along the core-core interface prior to a possible kinking. The physical insight obtained is essential for optimal design of peel-stoppers.     

Fracture criterion for kinking cracks in a tri-material adhesively bonded joint under mixed mode loading

The fracture behavior of a composite/adhesive/steel bonded joint was investigated by using double cantilever beam specimens. A starter crack is embedded at the steel/adhesive interface by inserting Teflon tape. The composite adherend is a random carbon fiber reinforced vinyl ester resin composite while the other adherend is cold rolled steel. The adhesive is a one-part epoxy that is heat cured. The Fernlund-Spelt mixed mode loading fixture was employed to generate five different mode mixities. Due to the dissimilar adherends, crack turning into the adhesive (or crack kinking) associated with joint failure, was observed. The bulk fracture toughness of the adhesive was measured separately by using standard compact tension specimens. The strain energy release rates for kinking cracks at the critical loads were calculated by a commercial finite element analysis software ABAQUS in conjunction with the virtual crack closure technique. Two fracture criteria related to strain energy release rates were examined. These are (1) maximum energy release rate criterion (G_max) and, (2) mode I facture criterion (G_II = 0). They are shown to be equivalent in this study. That is, crack kinking takes place at the angle close to maximum G or G_I (also minimum G_II, with a value that is approximately zero). The average value of G_IC obtained from bulk adhesive tests using compact tension specimens is shown to be an accurate indicator of the mode I fracture toughness of the kinking cracks within the adhesive layer. It is concluded that the crack in tri-material adhesively bonded joint tends to initiate into the adhesive along a path that promotes failure in pure mode I, locally.     

Numerical simulation of cracking processes in dissimilar media

This paper presents a novel numerical technique to simulate the crack propagation process for inhomogeneous solids. The criteria of crack growth and path selection are expressed in terms of some simple functions of fracture energy and loading phase angle. The energy release rate and loading phase angle for mixed mode cracks are evaluated by the finite element method with a mesh-adaptive technique. Examples are considered to show that the method is effective for simulating and predicting crack paths in inhomogeneous solids. They include crack growth or kinking in bi-materials and damage processes of fiber-reinforced composites.     

7050铝合金Ⅰ-Ⅱ复合型疲劳裂纹研究

为研究7050铝合金在Ⅰ-Ⅱ型复合加载下疲劳裂纹扩展规律,在Amsler HFP5000高频试验机上利用Richard加载装置,完成紧凑拉剪(CTS)试样疲劳裂纹扩展试验,利用有限元对Ⅰ-Ⅱ复合型裂纹进行数值模拟,采用APDL命令流计算不同裂纹长度的应力强度因子,并引入最大周向应力准则计算裂纹扩展角,用有限元计算等效应...     

Cohesive layer models for predicting delamination growth and crack kinking in sandwich structures

There are potentially two types of fracture that sandwich structures with strong and stiff facing sheets and lightweight cores are liable to suffer. These are the delamination growth at the face-sheet core interface and crack kinking into the sandwich core, respectively. The paper proposes computational models to simulate these failure mechanisms. The models employ the cohesive layer concept and are so constructed as to ensure that the crack advance is controlled by the critical value of strain energy release rate in mode I fracture. Of these, the first model can treat only delamination along a predetermined plane and is designated as CLD (cohesive layer delamination model). The performance of this model is thoroughly investigated in the light of experimental results. The influence of the key parameters of the model, viz. the thickness of the cohesive layer and the strength and stiffness of the cohesive layer material, have been studied. It is found that the model, as developed in this study, is fairly robust and is not sensitive to changes in parameters other than the critical strain energy release rate. The second model can track crack growth which is not predetermined in its direction. This it does by identifying the element in which the maximum principal tensile stress exceeds a critical value; once a crack is nucleated, the stress across the crack is relieved so that the right amount of energy is released when the crack is fully developed - much in the same manner as in a cohesive layer model. This model is designated as CLDK (Cohesive Layer Delamination and Kinking) model as it deals with interfacial delamination and crack kinking- whichever is the preferred mode of fracture. Experimental results of three sandwich specimens, viz. bottom restrained beams with 0° and –10° tilt angle, respectively, and a compressed beam, were used for comparison. The results indicate that the both the models are able to capture the initiation and track the growth of the interfacial delamination. The CLDK model is capable in addition to track the crack kinking into the core, and its subsequent return to the face sheet-core interface.     

Crack propagation from a pre-existing flaw at a notch root – II: Detailed form of the stress intensity factors at the initial crack tip and conclusion.

This paper pursues the study of crack kinking from a pre-existing crack emanating from some notch root. It was shown in Part I that the stress intensity factors at the tip of the small initial crack are given by universal (that is, applicable in all situations, whatever the geometry of the body and the loading) formulae; they depend only on the `stress intensity factor of the notch' (the multiplicative coefficient of the singular stress field near the apex of the notch in the absence of the crack), the length of the crack, the aperture angle of the notch and the angle between its bisecting line and the direction of the crack. Here we identify the universal functions of the two angles just mentioned which appear in these formulae, by considering the model problem of an infinite body endowed with a notch with straight boundaries and a straight crack of unit length. The treatment uses Muskhelishvili's complex potentials formalism combined with some conformal mapping. The solution is expressed in the form of an infinite series involving an integral operator, which is evaluated numerically. Application of Goldstein and Salganik's principle of local symmetry then leads to prediction of the kink angle of the crack extension. It is found that although the direction of the crack is closer to that of the bisecting line of the notch after kinking than before it, the kink angle is not large enough for the crack tip to get closer to this line after kinking, except perhaps in some special situations.     

Experimental study on fracture behaviour of three‐flawed sandstone specimens after high‐temperature treatments

The crack coalescence of rocks significantly affects the stability of rock engineering, and extensive studies have been performed on preflawed rock specimens without thermal treatment. However, the fracturing behaviour of preflawed specimens after thermal treatment has not been investigated comprehensively. In this study, three‐flawed sandstone specimens with different flaw inclinations after high‐temperature treatments were tested under uniaxial compression. Photographic, acoustic emission and digital image correlation techniques were used to investigate the crack initiation, propagation and coalescence behaviour. Experimental results show that the peak strength, elastic modulus and peak strain of the three‐flawed specimens were lower than those of intact specimens and that they gradually recovered with increasing flaw angle. The peak strength and elastic modulus first increased and then decreased, whereas the peak strain increased with temperature. Noncoalescence, indirect coalescence and direct coalescence were three patterns observed between the two adjacent pre‐existing flaws. Finally, the mechanism of high temperature in alteration of the mechanical properties of sandstone was revealed through microobservations.     

Reprint of “Cracks in inhomogeneous materials: Comprehensive assessment using the configurational forces concept”

The concept of configurational forces is applied to demonstrate the application of the concept of configurational forces in the numerical simulation of crack growth and fracture processes. It is shown, how material property variations at an interface affect the crack driving force and how the criterion of maximum dissipation is used to evaluate the direction of crack propagation. Fatigue crack growth experiments were conducted on diffusion welded bimaterial specimens consisting of a high-strength steel and soft ARMCO iron. Two cases are considered: (1) specimens with an interface perpendicular to the initial crack orientation, and (2) specimens with an inclined interface. The numerical simulation with the concept of configurational forces show that not only variations of the elastic modulus and/or the yield stress have a tremendous influence on the crack driving force, the crack growth rate, and the curvature of the crack path, but also the thermal residual stresses that resulted from a rather small difference of the coefficient of thermal expansion.