二维平纹编织复合材料压缩力学行为研究

通过对二维平纹编织复合材料压缩力学性能进行三维有限元计算,对复合材料在沿垂直于编织平面方向和纤维束方向压缩后的微观应力场进行了分析。对于沿垂直于编织平面方向的压缩(即横向压缩),发现在经纬向纤维束相互交替处基体材料的横向压应力最大,是其薄弱部位。分析表明,与0°/90°层合板相比,编织复合材料受到横向压缩载荷时,自由边界处层间剪应力较大,容易发生层间裂纹,与实验观察相吻合。通过改变纤维束截面尺寸和纤维束之间间隙大小研究了破坏面角度的变化。     

陶瓷基复合材料基体唇形裂纹失效规律研究

建立横向拉伸载荷下的唇形裂纹数学模型,采用复变函数的方法,通过保角映射,推导了唇形裂纹尖端应力场和位移场的解析解,建立了唇形裂纹的应力强度因子准则和最大能量释放率准则,结合算例分析陶瓷基复合材料基体唇形裂纹的几何参数、外载荷和纤维分布对失效准则的影响规律.结果 表明,(1)裂纹尖端应力场和位移场的解析解与有限元计算结果进行对比,验证了方法的有效性;(2)相较于Griffith裂纹和椭圆裂纹,基于唇形裂纹的失效准则对裂纹尖端的敏感性更高,适用于预测脆性陶瓷基体裂纹的扩展;(3)对于不同几何参数的唇形裂纹,采用最大能量释放率准则的基体裂纹的扩展速率要大于应力强度因子准则.     

基于弹性-塑性内聚力模型的纤维拔出界面宏观行为分析

用解析法分析了单纤维从聚合物基体中的拔出过程,采用弹性—塑性内聚力模型模拟裂纹的扩展和界面失效,确定了临界纤维埋入长度,该值区分两种不同长度的纤维拔出过程. 在纤维拔出过程,界面经历不同的阶段. 纤维埋长小于临界长度时,界面的脱粘载荷与纤维的埋长成正比;超过临界长度后,界面的脱粘载荷近似为常数. 分析了界面参数对脱粘载荷的影响:增加界面的剪切强度和界面的断裂韧性,或减小界面裂纹萌生位移,均能提高界面的脱粘载荷;界面脱粘后无界面摩擦应力时,拔出载荷—位移曲线的峰值载荷等于界面的脱粘载荷;界面摩擦应力存在时,使峰值载荷大于脱粘载荷,需要较长的纤维埋入长度和较大的界面摩擦应力.      

单向纤维增强复合材料的载荷集中因子和裂纹扩展统计理论

本文对单向纤维增强复合材料裂纹扩展统计理论的两个主要论据进行了研究.1.由于纤维断裂引起的应力重分布的计算:文中采用通常的剪滞模型,进行了弹塑性形变理论分析和弹性分析,弹性分析的结果与已有的J.M.Hedgepeth 的结果相符.2.裂纹扩展模型的探索:这里舍弃沿用Gücer-Gurland.Rosen 链式模型,提出了逐渐扩大断裂层的模型,参照A.S.Argon 等人的方法进行了统计分析,得到了单向纤维增强复合材料拉伸强度的计算公式.     

基于近场动力学理论的单层板损伤分析方法

基于键基近场动力学理论(PD)构建复合材料的单层板模型,在本构方程中引入连续变化的键刚度,以表示复合材料中与纤维夹角相关的力学性能的变化。推导了与复合材料主方向热膨胀系数相关的任意方向的热膨胀系数,为热载荷的施加提供理论支持,将热载荷和力学载荷统一于本构方程中。最后,对预先存在裂纹的单层板进行裂纹扩展模式的仿真,得到预期的仿真结果,并与已有实验结果对比验证模型的有效性。     

复合材料桥纤维拔出的动态裂纹模型

在无限大正交各向异性体弹性平面上对复合材料桥纤维平行自由表面的内部中央裂纹提出了桥纤维拔出的动态裂纹模型。通过复变函数将其转化为Reimann-Hilbert混合边界值问题。求得了裂纹在坐标原点受载荷Px/t、Px2/t作用的解析解。利用这一解析解可通过迭加原理求得任意复杂问题的解。     

复合材料桥连的断裂动力学模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂．在数学方法的处理上有很大困难，至今人们研究大多是桥连的静力学问题．而对其动力学问题研究得很少。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时．其纤维也连续地断裂。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。通过复变函数论的方法，将所讨论的问题转化为Riemann—Hilbert问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力P及阶跃载荷作用下位移、应力和动态应力强度因子的解析解，并通过叠加原理，最终求得了该模型的解。     

复合材料的裂纹动力学问题的模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂，在数学方法的处理上有很大困难，至今人们研究的大多是桥连的静力学问题，而对其动力学问题研究得很少。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时，其纤维也连续地断裂。通过复变函数论的方法，将所讨论的问题转化为Riemann-Hilbet问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力Px/t及均布载荷作用下位移、应力和动态应力强度因子的解析解，并通过迭加原理，最终求得了该模型的解。     

纤维增强聚合物布加固裂纹木梁的稳定性

梁中横向裂纹等效为无质量内部转动弹簧，假定纤维增强聚合物(FRP)布与梁表面紧密粘贴，建立了考虑轴向压力二阶效应FRP 布加固裂纹梁线性弯曲的控制方程，并得到其显式解析通解．在此基础上，研究了FRP加固简支裂纹木梁的稳定性，通过数值求解方程，分析了纤维增强聚合物(CFRP)布含量、裂纹深度和位置以及数量等因素对CFRP 布加固简支裂纹杉木梁临界载荷的影响，结果表明：CFRP 加固可明显减小裂纹深度和数量等对裂纹杉木梁临界载荷的影响，且裂纹处弯矩较大或裂纹较深时加固效应愈加显著；CFRP 加固裂纹木梁临界载荷随CFRP 布加固层含量的增加而增加，但当CFRP 布含量达到一定值后，进一步增加CFRP 含量对CFRP加固裂纹梁临界载荷提高并不明显．     

倾斜内锁型三维机织陶瓷基复合材料的细观力学模型

利用平均化方法提出了倾斜内锁型三维机织陶瓷基复合材料弹性性能分析的三维细观力学模型,对材料的弹性性能进行了预测。这个力学模型考虑了倾斜内锁型三维机织陶瓷基复合材料经向纤维束的弯曲和纬向纤维束的平直,纤维束的横截面形状尺寸和相邻纤维束之间的孔洞以及材料制造过程中碳纤维性能下降对弹性性能的影响。基于层合板理论,提出两种单胞应变状态假设分别对材料的九个弹性常数进行了推导计算,结果表明两种方法理论的预测值非常接近。计算结果与实验值比较吻合,表明所提出的细观力学模型是合理的,可以为纺织陶瓷基复合材料的优化设计提供有价值的参考。     

Two-parameter failure criterion for elastoplastic bodies with mode I cracks

The generalized Dugdale crack model is used to formulate two-parameter failure criteria for the cases of quasibrittle state and developed plastic zones at a mode I crack tip. The failure criteria relate the fracture strength characteristics and the stress mode at the crack tip through the plastic constraint factor. The critical state of bodies with cracks under uni-and biaxial loading is analyzed in the cases of plane stress and plane strain using the Tresca and von Mises yield criteria. A small-scale yield criterion, which is an analytic relation between the critical stress intensity factor and T-stresses, is established __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 7, pp. 47–57, July 2007.     

大开口复合材料层合板强度破坏研究

复合材料层合板的各向异性及非均质,使得复合材料层合板内部的破坏形式非常复杂.在复合材料结构的设计中,为满足制造及使用功能上的需求,在复合材料层合板承力结构件上不可避免地需要设计各种开口.然而,含大开口复合材料层合板的强度破坏问题变得更为复杂,使得现有的强度理论面临新的挑战.针对碳纤维增强复合材料大开口层合板受单向拉伸载荷作用下的强度破坏问题进行了数值分析和实验研究.首先,根据Hashin准则和刚度退化模型,对含不同圆形开口尺寸的[0]_(10)单向铺层、[0/90]_5和[±45]_5正交铺层的层合板,进行了单向拉伸载荷作用下渐进失效的数值模拟分析,获得了对应结构的极限载荷和破坏模式.在此基础上,采用数字图像相关方法,进行复合材料大开口层合板强度破坏的实验研究.研究结果表明,大开口复合材料层合板在单向拉伸加载下主要呈现脆性破坏形式,破坏起始位置处于应力集中区.此外,破坏强度和失效模式与复合材料铺层方式和开口尺寸大小密切相关.其中[±45]_5铺层的开口层合板承载能力最弱,分层破坏最严重.开口尺寸越大,结构的极限载荷值越低.同实验测试结果相比,数值模拟对复合材料层合板的损伤失效分析略显不足,往往很难全面分析复合材料层合板破坏失效过程中的各种因素的影响.     

Mircopolar Model of Fracture for Composite Material

Failure of a composite is a complex process accompanied by irreversible changes in the microstructure of the material. Microscopic mechanisms are known of the accumulation of damage and failure of the type of localized and multiple ruptures of the fibers delamination along interphase boundaries, and also mechanisms associated with fracture of fibers. In this work, we propose a mathematical model of the local mechanism of failure of a composite material randomly reinforced with a system of short fibers. We implement the Cosserat moment model of crack tip for filament material, reinforced with whiskers or in fiber- reinforced polycrystalline materials. It is assumed that the angular distribution of the fibers is isotropic and the elastic characteristics of the fibers are considerably higher than the elastic constants of the matrix. We implement the homogenization procedure for the effective Cosserat constants similarly to the effective elastic constants. The singular solution in the vicinity of the crack tip in the Cosserat moment model is found. Using this solution, we examine the bending stresses in the filaments due to effective moment stresses in the material. The constructed model describes the phenomenon of fracture of the fibers occurring during crack propagation in those composites. The following assumptions are used as the main hypotheses for the micromechanical model. The matrix contains a nucleation crack. When the load is increased the crack grows and its boundary comes into contact with the reinforcing fibers. A further increase of the stress causes bending of the fiber. When~the fiber curvature reaches a specific critical value, the fiber ruptures. If the stress at infinity is given, the fibers no longer delay the development of failure during crack propagation The degree of bending distortion of the fiber in the vicinity of the boundary of the crack is determined by the moment model of the material. The necessity to take into account the moment stresses in the failure theory of the reinforced material was stressed in [Muki and Sternberg (1965) Zeitschrift f angew Math und Phys 16:611–615; Garajeu and Soos (2003) Math Mech Solids 8(2):189–218; Ostoja-Starzewski et al (1999) Mech Res Commun 26:387–396]. The moment Cosserat stresses were accounted also for inhomogeneous biomechanical materials by Buechner and Lakes (2003) Bio Mech Model Mechanobiol 1: 295–301. We should also mention the important methodological studies [Sternberg and Muki (1967) J Solids Struct 1:69–95; Atkinson and Leppington (1977) Int J Solids Struct 13: 1103–1122] concerned with the moment stresses in homogeneous fracture mechanics.     

Quasi-brittle fracture diagram of structured bodies in the presence of edge cracks

The Neuber-Novozhilov approach is used to obtain necessary and sufficient fracture criteria. Using a modified Leonov-Panasyuk-Dugdale model, simple relations for the critical fracture parameters are derived for opening mode edge cracks for the case where the diameter of the prefracture zone coincides with the diameter of the plasticity zone. These relations are suitable for studying fracture where the crack length is negligibly small. A fracture diagram using critical stresses under both criteria is proposed for a wide range of crack length. At a certain level of loading, three regions are identified, in the first of which the crack is stable, in the second, the crack extends but remains stable, and in the third, the crack is unstable. Experimental data on the fracture of specimens with edge cracks are obtained. It is established that the theoretical critical fracture curves are in good agreement with the obtained critical parameters for flat tensile specimens with two collinear edge cracks.     

Slope stability analysis using fracture mechanics approach

The analysis of slope failure is complicated due to its mechanism as well as the geological history of the slope. In classical slope stability analysis, the slope failure is assessed using the basic continuum mechanics or the limit equilibrium approach. This analysis, however, must be slightly modified when tension cracks exist at the upper edge of the slope. From post-mortem analyses of slope failures in the past, it was found that tension cracks had some considerable importance affecting the failure mechanism of the slope.Some attempts had also been directed in the past towards simple vertical cut slopes with tension cracks. Terzaghi [40]. Using an elasto-plastic model and a few analytical assumptions, Terzaghi [40] obtained an expression for a critical height of a vertical cut slope. However, since the stress distribution is changing during crack propagation, the failure mechanism cannot in principle conjecture from a static (undisturbed) stress field alone. In this report, a unified numerical approach combining Finite Element Method (FEM), fracture mechanics and remeshing technique is used to model the failure analysis of vertical cut slopes with tension cracks. Using this approach, the crack can be extended incrementally under a certain energy based failure criterion (Strain Energy Density criterion - SED) while taking into account the existing stress singularity field in the vicinity of the crack-tip.The use of fracture mechanics in soils, however, poses some difficulties in obtaining its relevant parameters due to the granular structure of soil. Because of its granular structure, the shear strength of soil depends on cohesion and the applied confining pressure making the determination of fracture parameters in soil difficult. These relevant parameters are not yet available in literature and hence, certain assumptions and interpolations had to be made in obtaining these parameters. Future research direction in this area is badly needed. Using a certain set of soil properties and crack geometry, a series of curves relating parameters S (Strain Energy Density Factor) and a non-dimensional variable N ( = H/C - slope height vs. crack distance from edge ratio) for vertical cut slopes can be constructed through numerous parametric studies. These curves can serve as a quick guide to obtain the critical height for a given vertical cut slope geometry with tension cracks.     

Loss of adhesion at the tip of an interface crack

A model is constructed to analyze adhesive bond failure at the tip of an interface crack. The model is based on the assumption that there are zones of bounded cohesive tensile and shear stresses near a crack tip. Within the context of certain broad a-priori assumptions on the distributions of certain stress and displacement components in the cohesive zones, the requirement thatall stresses in the two materials remain bounded provides a method to compute the specific details for these zones. It is assumed that bond failure occurs when the extension of the bond fiber at the crack tip exceeds a critical value. For an interface crack in a uniform tension field computations for two alternate formulations suggest that this failure criterion is independent of the precise distribution of the cohesive stresses, but rather depends only upon their averaged values. Combined loading with a dominant tensile component has also been analyzed. If the critical extension of bond fibers and the maximum value of the cohesive tensile stress are known, the model provides the maximum allowable interface stresses for given crack dimension and material parameters.     

Construction of quasi-brittle and quasi-ductile fracture diagrams based on necessary and sufficient criteria

Materials with a regular structure characterized by quasi-brittle and quasi-ductile fractures are considered in the case where the characteristic linear dimension of the structural element is known. Necessary and sufficient fracture criteria are constructed using the Neuber-Novozhilov approach. A modified Leonov-Panasyuk-Dugdale model for an opening mode crack is proposed where the width of the prefracture zone coincides with the width of the plasticity zone. For the critical parameters of quasi-brittle fracture (tensile stress, length of prefracture zones, stress intensity factors), relations are obtained that allow material fracture to be considered in the case where the crack length is negligible compared to the characteristic linear dimension of the structural element. A fracture diagram obtained using the critical stresses calculated from the necessary and sufficient criteria is considered in a wide range of crack lengths. The elastoplastic problem of extension of a plate with a central crack is solved using the finite-element method. The dimensions and shape of the plastic zone near the crack tip are determined for different levels of loads corresponding to quasi-brittle and quasi-ductile fracture. The obtained results are analyzed to estimate the width of the prefracture zone and the critical crack opening.     

近场动力学在断裂力学领域的研究进展

近场动力学采用非局部积分计算节点内力, 利用统一数学框架描述空间连续与非连续, 避免了非连续区局部空间导数引起的应力奇异, 数值上具有无网格属性, 可自然模拟材料结构的断裂问题. 本文概述了近场动力学的弹性本构力模型, 系统介绍了近场动力学临界伸长率、临界能量密度以及材料强度相关的键失效准则. 详细介绍了近场动力学在断裂力学领域的研究进展, 包括断裂参数能量释放率与应力强度因子的求解、J积分、混合型裂纹、弹塑性断裂、黏聚力模型、动态断裂、材料界面断裂以及疲劳裂纹扩展等. 最后讨论了断裂问题近场动力学研究的发展方向.      

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

An extended strain energy density failure criterion by differentiating volumetric and distortional deformation

We extend Sih’s strain energy density criterion (Sih, 1974) for crack kinks and material failure by weighting differently the volumetric and distortional parts in the extended strain energy density factor. The work is inspired by the factor that failure by microscopic shearing governed by distortion and microscopic separation controlled by hydrostatic tension represent distinct deformation processes, and should be treated differently as we count their influences to material failure. With the weight parameter introduced to the extended strain energy density factor criterion, we explain satisfactorily several critical experiments which reported crack kink in samples subjected to mixed-mode loading. The extended strain energy density idea is also used to derive a generalized pressure-dependent yielding criterion, which supplies a theoretical basis for those novel strength criteria for materials like bulk metallic glasses. Corresponding methods to determine the two material parameters, the critical strain energy density factor and the weight parameter quantifying the relative contribution by distortion over volumetric deformation, are discussed.