高温环境下纤维复合材料蠕变损伤的细观机理研究

首先利用复合材料纤维断裂单胞模型，编制蠕变损伤子程序，对单胞模型进行蠕变损伤分析。分析了纤维／基体弹性模量比对蠕变变形、蠕变损伤以及应力场的影响。从计算结果发现，蠕变损伤首先在纤维断裂尖端起始，然后沿着一定的角度向基体外围延伸，直至完全损伤，而且纤维／基体模量比对高温环境下的复合材料蠕变损伤产生很大的影响；纤维与基体的模量相差越大，复合材料越容易变形，抵抗蠕变变形的能力就越小，蠕变损伤越严重。经过对不同韧性的基体材料进行研究，发现基体韧性低的复合材料蠕变损伤明显高于高韧性基体复合材料，表明低韧性基体复合材料抵抗蠕变破坏的能力较低。     

纤维增强金属基复合材料整体蠕变性能的细观力学分析

高温下金属基复合材料的蠕变主要由基体蠕变和界面扩散蠕变两部分构成,以往的研究中常常只考虑其中一种蠕变机理,从而导致得到的规律具有较大的局限性．本文提出了一种可预测金属基复合材料整体蠕变性能的细观力学方法,同时考虑了基体蠕变和界面扩散蠕变两种蠕变机理,导出了具有张量形式并满足不可压缩性的界面扩散蠕变应变表达式．采用Mori-Tanaka法和自洽法二者结果的平均以便更准确地计算纤维中的应力,揭示了两种蠕变机理相互影响的竞争关系．研究了恒定双轴荷载下的总体蠕变和固定位移约束下的应力松弛这两种常见蠕变问题,探究了基体蠕变与界面扩散蠕变两种蠕变机理在总蠕变中发挥的作用,考察了不同加载条件和不同纤维体积分数对复合材料整体蠕变行为的影响．     

界面特性对短纤维金属基复合材料蠕变行为的影响

基于短纤维增强金属基复合材料(MMC)的单纤维三维模型(三相)，利用粘弹性有限元分析方法对影响金属基复合材料的蠕变行为的因素进行了较为系统的分析。研究中主要讨论了界面特性和纤维取向角对金属基复合材料的蠕变性能的影响。研究结果发现，界面特性诸如厚度、模量和应力指数都对纤维最大轴应力和稳定蠕变率产生影响：稳态蠕变率随界面模量的增大而逐渐减小，当高于基体模量时基本保持不变；纤维轴应力的变化与蠕变率正好相反。稳态蠕变率随界面厚度、应力指数的增加而增大；而轴应力则随之减小。同时不同的纤维取向也影响金属基复合材料蠕变时的轴应力分布和稳态蠕变率。     

兼具基体蠕变和轻度界面蠕变的复合材料的本构特性

金属基复合材料的基体蠕变特性在一定的温度与应力下表现明显，当界面为非理想并具有粘弹性性质时，其对复合材料的整体蠕变亦有重要影响．本文利用具有非理想界面的复合材料的Ｍｏｒｉ－Ｔａｎａｋａ方法，研究了既具有基体蠕变又具有界面轻度蠕变的复合材料的本构关系，给出了各微结构参量对复合材料整体蠕变特性的影响     

兼具基体蠕变和轻度界面蠕变的复合材料的本构特性

金属基复合材料的基体蠕变特性在一定的温度与应力下表现明显，当界面为非理想并具有粘弹性性质时，其对复合材料的整体蠕变亦有重要影响．本文利用具有非理想界面的复合材料的Ｍｏｒｉ－Ｔａｎａｋａ方法，研究了既具有基体蠕变又具有界面轻度蠕变的复合材料的本构关系，给出了各微结构参量对复合材料整体蠕变特性的影响     

汽车玻璃高温蠕变特性的实验与理论研究

通过对汽车用挡风玻璃进行不同温度、不同应力水平下的单轴拉伸蠕变实验,揭示了其非线性蠕变变形特性.实验发现,不同温度和载荷下的蠕变曲线形状类似,高温(550℃～590℃)下该材料的减速蠕变变形不很明显,主要表现为稳态和加速蠕变变形,且其断裂数据符合Monkman-Grant关系,应力和温度对其影响不显著.基于实验结果,本文建立了非线性高温蠕变本构方程.理论预测与实验数据进行了比较,结果表明所提出的本构方程能较好地描述汽车玻璃高温蠕变变形的全过程.     

切应力对牛骨蠕变变形的影响研究

骨尤其是湿骨,在恒定载荷作用下会发生蠕变变形。为了确定切应力是否影响骨的蠕变变形,采用对骨薄板试样分别施加集中载荷和均布载荷的方式,测量试样挠度实时的变化曲线。结果显示,在载荷恒定时,骨试样的挠度随时间不断增加,体现了典型的蠕变特性。集中荷载下骨的蠕变变形远大于均布荷载下骨的蠕变变形,湿骨的蠕变位移比干骨高近7倍。分析认为,对试样粘弹性性质的影响不仅有正应力的作用,也有切应力的作用;切应力产生的蠕变变形约为正应力所产生蠕变变形的0.85倍。     

高应力高水压下砂岩三轴蠕变特性试验研究

对砂岩进行高围压高水压条件下的三轴压缩蠕变试验。试验表明,在整个加载过程中,孔隙水压力主要起到增强轴向变形和横向变形的作用;但在加载的初始阶段,孔隙水压力在一定程度上抑制了轴向变形。当应力水平高于屈服应力时,横向蠕变速率明显大于轴向蠕变速率,且横向蠕变率先进入加速蠕变阶段。本文提出一个新的非线性黏性元件,并引入一个能判断是否进入加速蠕变阶段的计时器元件,组建一个非线性黏塑性加速蠕变启动模型,将该黏塑性模型与Burgers模型串联,构建一个新的非线性黏弹塑性蠕变模型,推导了该模型在常规三轴应力状态下的本构方程。基于试验结果,通过对非线性优化算法(BFGS)的Matlab编程,实现对本文提出蠕变模型的参数识别,识别效果比较理想。对比试验曲线与拟合曲线,二者相当吻合,验证了新提出的非线性黏弹塑性蠕变模型的正确性。     

砂岩峰后蠕变失稳特性试验研究

为揭示岩石峰后蠕变失稳特性,利用MTS815电液伺服岩石力学试验系统,对红砂岩进行了峰后蠕变试验,采用线性回归的方法估算了蠕变失稳时间,探讨了损伤、初始围压以及围压增量对峰后蠕变特性的影响规律。结果表明:随着损伤的加大,蠕变变形速率增大,蠕变失稳时间缩短。同一围压下,各级应力水平下蠕变量随应力水平的增大而增大,且服从指数函数关系。增大初始围压,各级应力水平下的蠕变变形速率均有一定程度的降低,且对等速蠕变阶段蠕变速率的影响较加速蠕变阶段蠕变速率大,蠕变失稳时间延长。施加围压增量,大大减小了蠕变变形速率,促使了加速蠕变阶段向等速蠕变阶段转化,明显延长了蠕变失稳时间,减小了蠕变形变量。研究结果为工程围岩的长期稳定性分析提供了参考依据。     

考虑强度参数劣化的砂岩非定常蠕变模型研究

为了更好地研究岩石在不同应力作用下的蠕变变形特性,采用TAW-2000岩石试验系统对围岩开展了单轴蠕变试验;将西原体模型蠕变参数转化为关于时间的函数,且将强度参数也转化为关于时间的函数,进一步将屈服函数非线性化,构建出一种新型的非定常分数阶蠕变模型,并通过试验数据来验证模型的正确性与合理性。结果表明:轴向蠕变变形总体上随应力水平的增加而增加,但蠕变变形占总变形的比例先增大再减小,反映了试样内部结构的不断损伤与劣化过程;该模型不仅准确地反映了衰减和稳定蠕变阶段的蠕变特性,也克服了传统西原体模型难以描述加速蠕变的缺点;同时考虑到蠕变参数随应力和时间变化发生的劣化,可以更好地描述不同应力状态下岩石的蠕变损伤过程,为建立非定常蠕变模型和确定模型参数提供一个新的思路。     

A variational method for unidirectional fiber-reinforced composites with matrix creep

A variational method is developed for analyzing the matrix creep induced time-dependent change in fiber stress profiles in unidirectional composites. A functional of admissible profiles of fiber stress rate is presented by supposing a fiber broken in matrix as well as a fiber pulled out from matrix. The functional is shown to have the stationary function satisfying an incremental differential equation based on the shear lag assumption. Then, the stationary function is approximately determined by assuming bilinear profiles of fiber stress and a power law of matrix creep, leading to analytical solutions for the time-dependent change in fiber stress profiles. The solutions are verified on the basis of an energy balance equation and a finite difference computation. Moreover, it is shown that the solution for the fiber pull-out model agrees well with an experiment on a single carbon fiber/acrylic model composite if the initial slip at fiber/matrix interface is taken into account. In addition, the solution for the fiber breakage model is used for evaluating the characteristic time in long-term creep rupture of unidirectional composite.     

A homogenization theory for time-dependentnonlinear composites with periodic internal structures

A homogenization theory for time-dependent deformation such as creep andviscoplasticity of nonlinear composites with periodic internal structures is developed. To beginwith, in the macroscopically uniform case, a rate-type macroscopic constitutive relation betweenstress and strain and an evolution equation of microscopic stress are derived by introducing twokinds of Y-periodic functions, which are determined by solving two unit cell problems.Then, the macroscopically nonuniform case is discussed in an incremental form using thetwo-scale asymptotic expansion of field variables. The resulting equations are shown to beeffective for computing incrementally the time-dependent deformation for which the history ofeither macroscopic stress or macroscopic strain is prescribed. As an application of the theory,transverse creep of metal matrix composites reinforced undirectionally with continuous fibers isanalyzed numerically to discuss the effect of fiber arrays on the anisotropy in such creep.     

On the effect of fiber creep-compliance in the high-temperature deformation of continuous fiber-reinforced ceramic matrix composites

Creep models for unidirectional ceramic matrix composites reinforced by long creeping fibers with weak interfaces are presented. These models extend the work of Du and McMeeking (1995) [Du, Z., McMeeking, R. 1995. Creep models for metal matrix composites with long brittle fibers. J. Mech. Phys. Solids 43, 701–726] to include the effect of fiber primary creep present in the required operational temperatures for ceramic matrix composites (CMCs). The effects of fiber breaks and the consequential stress relaxation around the breaks are incorporated in the models under the assumption of global load sharing and time-independent stochastics for fiber failure. From the set of problems analyzed, it is found that the high-temperature deformation of CMCs is sensitive to the creep-compliance of the fibers. High fiber creep-compliance drives the composite to creep faster, leading however to greater lifetimes and greater overall strains at rupture. This behavior is attributed to the fact that the greater the creep-compliance of the fibers, the higher the creep rate but the slower the matrix stress relaxation – since the matrix must deform with a rate compatible with the more creep-resistant fibers – and therefore the less the load carried by the main load-bearing phase, the fibers. As a result, fewer fibers fail and less damage is accumulated in the system. Moreover, the greater the creep-compliance of the fibers, the slower the matrix shear stress relaxation – and thus the lower the levels of applied stress for which this effect becomes important. The slower the shear stress relaxes, the slower the “slip” length increases. Due to the Weibull nature of the fibers, the fiber strengths at the smaller gauge length of the slip length are stronger; therefore fewer fibers undergo damage. Hence, high fiber creep-compliance is desirable (in the absence of any explicit creep-damage mechanism) in terms of composite lifetime but not in terms of overall strain. These results are considered of importance for composite design and optimization.     

STATISTIC MODELING OF THE CREEP BEHAVIOR OF METAL MATRIX COMPOSITES BASED ON FINITE ELEMENT ANALYSIS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｒｅｅｐｂｅｈａｖｉｏｒｏｆｓｈｏｒｔｆｉｂｅｒｒｅｉｎｆｏｒｃｅＭｅｔａｌＭａｔｒｉｘＣｏｍｐｏｓｉｔｅｓ (ＭＭＣｓ)ｄｅｐｅｎｄｓｏｎｔｈｅｆｏｌｌｏｗｉｎｇｆａｃｔｏｒｓ,ｓｕｃｈａｓｔｈｅｃｒｅｅｐｐｒｏｐｅｒｔｙｏｆｔｈｅｍａｔｒｉｘ ,ｅｌａｓｔｉｃａｎｄｆｒａｃｔｕｒｅｓｐｒｏｐｅｒｔｉｅｓｏｆｔｈｅｆｉｂｅｒ,ｇｅｏｍｅｔｒｉｃｐａｒａｍｅｔｅｒｓｏｆｔｈｅｆｉｂｅｒｓ,ａｒｒａｎｇｅｍｅｎｔｏｆｔｈｅｆｉｂｅｒｓａｎｄｔｈｅｐｒｏｐｅｒｔｙｏｆｔｈｅｆ…     

纤维增强韧性基体界面力学行为

分析了纤维增强韧性基体的界面力学行为及其失效机理,按剪滞理论和应变理化规律研究微复合材料的弹塑性变形和应力状态,讨论了幂硬化和线性硬化基体的弹塑性变形和界面应力分布,并给出纤维应力和位移的表达式。按最大剪应力强度理论建立了纤维／基体界面失效准则,推导出弹塑性界面失效的平均剪应力随纤维埋入长度的变化关系。     

Creep-rupture lifetime simulation of unidirectional metal matrix composites with and without time-dependent fiber breakage

A numerical simulation for predicting the axial creep-rupture lifetime of continuous fiber-reinforced metal matrix composites is proposed, based on the finite element method. The simulation model is composed of line elements representing the fibers and four-node isoparametric plane elements representing the matrix. While the fibers behave as an elastic body at all times, the matrix behaves as an elasto-plastic body at the loading process and an elasto-plastic creep body at the creep process. It is further assumed in the simulation that the fibers are fractured not only in stress criterion but time-dependently with random nature. Simulation results were compared with the creep-rupture lifetime data of a boron-aluminum composite with 10% fiber volume fraction experimentally obtained. The simulated creep-rupture lifetimes agreed well with the averages of the experimental data. The proposed simulation is further carried out to predict a possibility of creep-rupture for the composite without time-dependent fiber breakage. It is finally concluded that the creep-rupture of a boron-aluminum composite is closely related with the shear stress relaxation occurring in the matrix as well as time-dependent fiber breakage.     

A time dependent micro-mechanical fiber composite model for inelastic zone growth in viscoelastic matrices

In this work, a fiber composite model is developed to predict the time dependent stress transfer behavior due to fiber fractures, as driven by the viscoelastic behavior of the polymer matrix, and the initiation and propagation of inelastic zones. We validate this model using in situ, room temperature, micro-Raman spectroscopy fiber strain measurements. Multifiber composites were placed under constant load creep tests and the fiber strains were evaluated with time after one fiber break occurred. These composite specimens ranged in fiber volume fraction and strain level. Comparison between prediction and MRS measurements allows us to characterize key in situ material parameters, the critical matrix shear strain for inelastic zones and interfacial frictional slip shear stress. We find that the inelastic zone is predominately either shear yielding or interfacial slipping, and the type depends on the local fiber spacing.