Effective elastic moduli of composites reinforced by particle or fiber with an inhomogeneous interphase

The solution of the strain energy change of an infinite matrix due to the presence of one spherical particle or cylindrical fiber surrounded by an inhomogeneous interphase is the basis of solving effective elastic moduli of corresponding composites based on various micromechanics models. In order to find out the strain energy change, the composite sphere or cylinder, i.e., the spherical particle or cylindrical fiber together with its interphase, is replaced by an effective homogeneous particle or fiber. Independent governing differential equations for each modulus of the effective particle or fiber are derived by extending the replacement method [J. Mech. Phys. Solids 12 (1964) 199]. As far as the strain energy changes of the infinite matrix subjected to various far-field stress systems are concerned, the present model is simple. Meanwhile, FEM analysis is carried out for a verification, which shows that the model can lead to rather accurate results for most practical interphases. Besides, to check the validity of the model further when the interactions among composite cylinders exist, the two problems of an infinite matrix containing two composite cylinders and the effective moduli of composites with the equilateral triangular distribution of composite cylinders are analyzed using FEM. The FEM results show that the model is still rather accurate, especially for the case of interphase properties varying between those of fiber and matrix. Therefore, composite spheres or cylinders are assumed as the effective homogeneous particles or fibers and simple expressions of the effective moduli of composites containing the composite spheres or cylinders are obtained. Furthermore, the present model is compared with some existing models that are based on very complicated derivations.     

基于辛弹性力学解析本征函数的有限元应力磨平方法

在实际工程结构的结构强度与优化等力学数值分析中,应力计算结果的精度是非常重要的。有限元法是得到最广泛应用的一类数值方法,并形成了众多通用的有限元程序系统。这些程序系统采用的几乎都是基于最小总势能的位移法,虽然其分析给出的有限元位移场具有较高的精度,但所得到的有限元应力场的精度较位移场大大降低。基于极坐标辛对偶体系所提供的平面弹性力学的解析辛本征展开解,并借用有限元程序系统所得到的节点位移,本文提出了一个应力分析的改进方法。数值结果表明,本方法给出的应力分析精度得到大幅提高,并具有良好的数值稳定性,可用于有限元程序系统的后处理,以提高应力尤其是关键区域应力的分析精度。     

Prestressing effect of twisted fiber yarn in composites

By utilizing the behavior of twisted fiber yarn, prestressing composites can be made. In unloading, the matrix of the prestressing composite is applied in compression. If tensile load acts on the composite, the value of the tensile stress acting on the materials may decrease greatly. By twisted fiber yarn, therefore, the tensile strength of the composites can be improved. An analysis for extension of continuous fiber yarns is made here and residual stress after curing is taken into account. The prestressing intensity in the composite depends on the twist of fiber yarn. The photoelastic test and the analysis of electron micrograph are performed, and the theoretical method for calculating prestressing effect is presented in this paper.     

三维连续与非连续变形分析

将石根华博士所提出的二维非连续变形分析——Discontinuous Deformation Analysis（DDA）方法扩展到三维情况,并对三维不连续块体进行有限元网格剖分,即块体之间的接触采用DDA描述,块体内部的位移场和应力场则采用有限单元法描述,从而将三维DDA与有限元方法结合起来,增强了DDA方法与有限元方法解决实际工程问题的能力,实现了三维连续与非连续变形分析.给出了基本公式的推导过程和各子矩阵的形式.典型接触、碰撞算例证明了所提出方法的有效性和正确性.     

Experimental Investigation of Wheel/Rail Rolling Contact at Railhead Edge

Wheel-rail rolling contact at railhead edge, such as a gap in an insulated rail joint, is a complex problem; there are only limited analytical, numerical and experimental studies available on this problem in the academic literature. This paper describes experimental and numerical investigations of railhead strains in the vicinity of the edge under the contact of a loaded wheel. A full-scale test rig was developed to cyclically apply wheel/rail rolling contact load to the edge zone of the railhead. An image analysis technique was employed to determine the railhead vertical, lateral and shear strain components. The vertical strains determined using the image analysis method have been validated with the strain gauge measurements and used for the calibration of a 3D nonlinear Finite Element Model (FEM) that simulates the wheel/rail contact at the railhead edge and use suitable boundary conditions commensurate to the experimental setup. The FEM was then used to determine other states of strains.     

复合材料中加捻纤维束的预应力效应

利用加捻纤维束的力学性能,可以制成预应力复合材料.在不受外力的情况下,这种复合材料的基体受压应力的控制.当有拉伸力作用的时候,材料所受的拉应力值会大大减少,这样,利用这种现象可以提高纤维增强复合材料的抗拉强度.本文将分析连续的加捻纤维束受拉状况并考虑到材料固化后产生的残余应力的影响.这种方法所产生的预应力取决于纤维束的加捻程度.本文将给出计算这种预应力的方法,而且还提供加捻纤维束产生预应力的光弹测试结果以及这种纤维预应力复合材料的纤维破坏的特征.     

强度折减有限元在堤防稳定分析中的应用研究

采用强度折减有限元法,研究非均值堤防塑性区的开展特征和失稳破坏过程,并根据堤防应力场分布和临界滑动面的形成来分析堤岸整体失稳破坏的机理。结果表明:把强度折减有限元应用到非均质土层的堤防边坡稳定分析,在理论上、数值模拟实现上都是可行的; 有限元静力平衡和位移计算不收敛作为堤防边坡整体失稳的标志,同时考察滑移面上某些特征点的位移与荷载增量的关系、参考土体内坡脚向坡顶上方贯通的塑性滑动带的形成,可以准确判别堤防边坡失稳破坏状态; 此外,堤防失稳时的应力场分布特征表明,有必要应用张拉-剪切复合屈服准则,以真实反映堤防边坡的临界滑动面形成机理。     

Strength probability of unidirectional hybrid composites

A probability strength analysis of an unidirectional three-component hybrid composite (HC) is carried out for the cases of high modulus (HM) fibres and low modulus (LM) elastic fibers regularly embedded in a low elastic modulus matrix. Both single layer intraply and multilayer HC are considered. The fiber strength is assumed to be a random variable with a Weibull distribution. Breaking of the HM fibers are accumulated initially while probability of LM fiber failure is low. Failure modes tend to be covered by the two extreme cases of alternative failure of HM and LM fibers only. These modes can be categorized by using graph technique. Developed are the algorithm for finding the most probable pattern of fiber breaking and method for estimating the strength and fiber damage of a HC. The stress level at which the LM fibers are found to break represents a lower bound of the HC strength. Damage of HM fibers in a three-component HC is much higher than in a two-component HM fibre composite. Negative ‘hybrid effect’ for strength is obtained.     

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.     

基于神经网络与光纤传感阵列的结构状态监测方法

以探测机敏复合材料与结构内应变、应力及损伤等物理状态的位置信息为目的．提出了一种新颖的可内埋于材料与结构内作为结构状态监测用的强度调制型光纤传感阵列网络，并采用人工神经网络来处理传感阵列输出的并行分布式传感信号，阐述了适用的Ｋｏｈｏｎｅｎ网模型及其变化形式，给出了仿真实验结果     

Effects of non-uniform strains on tensile fracture of fiber-reinforced ceramic composites

Effects of non-uniform strains on tensile fracture of fiber-reinforced ceramic–matrix composites have not been satisfactorily explained by existing mechanics-based models. In this paper, we use an exact model of fiber fragmentation under global load sharing conditions to predict fracture in three model problems in which non-uniform strains occur: (i) an end-constrained plate subject to a linear transverse temperature gradient; (ii) an internally-pressurized cylindrical tube with a linear through-thickness temperature gradient; and (iii) a rectangular beam under combined bending and tension. Fracture is assumed to occur when the global load reaches a maximum value. Approximations to the exact fragmentation model are also assessed, with the goal of decoupling the effects of two important parts of the computed stress–strain response: the rate of post-peak strain softening and the magnitude of the plateau “flow” stress once fiber fragmentation is complete. We find that for cases in which the fiber Weibull modulus is low and hence its plateau strength is high relative to its peak and the loading yields a sufficiently high strain gradient, the failure strain lies in the plateau regime. Consequently, the results can be predicted with good accuracy using a perfectly-plastic representation of the post-peak response. In contrast, for cases in which the fiber Weibull modulus is high, the failure strain lies in the softening portion of the curve. Here a linear-softening model is found to yield accurate results. A preliminary assessment of the model has been made by comparing predicted and measured bending/tension strength and failure strain ratios for one specific composite. The correlations appear good, though additional experiments are required in order to critically assess the model predictions over a range of loading scenarios.     

含大量随机分布细小纤维的单向纤维增强复合材料横向弹性性能参数计算的内嵌区域模型法

利用有限元方法求取单向纤维增强复合材料的横向弹性性能参数的计算模型包括三维模型、两维平面应变模型、单胞模型等等.由于单胞模型仅仅适用于纤维规则排列情况.在纤维随机分布且纤维大小亦为随机时,单向纤维增强复合材料横向弹性性能参数必须通过对于复合材料块体的计算才能获得.同时在随机分布纤维的数量增大时,三维模型和二维平面应变模型的计算量急剧增加,模型的处理能力不强.该文提出一种利用内嵌区域模型来计算含大量随机大小、随机分布细小纤维的单向纤维增强复合材料块体的横向弹性性能参数的方法,有效降低了计算量.在较低的计算费用下,能够快速获得单向纤维增强复合材料的横向弹性性能参数.     

极端环境下连续碳纤维增韧的陶瓷基复合材料的力学行为

连续纤维增韧的碳化硅复合材料(以下简称C/SiC),作为超高速飞行器热结构使用时,有可能在高温环境下受到高速撞击的作用,因此,掌握其在极端环境(高温、高应变率)下的力学性能是进行结构安全设计的基础。本文采用具有高温实验能力的分离式Hopkinson杆,在293~1273K温度范围内进行了动态压缩力学性能测试,研究了环境温度和加载速率对材料力学性能的影响。结果表明:C/SiC复合材料的高温压缩力学性能主要受应力氧化损伤和残余应力的共同影响。实验温度低于873K时,应力氧化损伤的影响很小,而由于增强纤维和基体界面残余应力的释放使界面结合强度增大,复合材料的压缩强度随温度的升高而增大;当实验温度高于873K时,应力氧化损伤加剧,其对压缩强度的削弱超过残余应力释放对强度的贡献,材料的压缩强度随温度的升高逐渐降低。由于应力氧化损伤受应变率的影响很大,当温度由873K升高至1273K时,高应变率下压缩强度降低的程度要比应变率为0.0001/s时低得多。     

Constitutive modeling of fiber composites with a soft hyperelastic matrix

This paper presents an experimental and numerical study of unidirectional carbon fiber composites with a silicone matrix, loaded transversally to the fibers. The experiments show nonlinear behavior with significant strain softening under cyclic loading. The numerical study uses a plane-strain finite element continuum model of the composite material in which the fiber distribution is based on experimental observations and cohesive elements allow debonding to take place at the fiber/matrix interfaces. It is found that accurate estimates of the initial tangent stiffness measured in the experiments can be obtained without allowing for debonding, but this feature has to be included to capture the non-linear and strain-softening behavior.     

Damage simulation of repaired composite laminate with rectangular cut-out

Repaired composite laminate with rectangular cut-out is studied by tensile experiments and FEM simulation in this paper. A user damage subroutine is implemented to ABAQUS to simulate the damage of laminate. Good agreements between experimental data and numerical results are reached. From the experiment and FEM simulation results, the integrality of composite structure is destroyed due to the existence of cut-out. However, the asymmetry of repair gets this effect farther. From damage simulation, the initiation and development of three types of damages are predicted. The matrix cracking occurs at first, fiber–matrix shear and fiber break are subsequent. The effect of matrix cracking to the carrying capacity of structure is very small, but when fiber–matrix shear and fiber break occur, dangerous to structure, the carrying capacity of structure loses rapidly.     

Multiscale analysis of the transverse properties of Ti-based matrix composites reinforced by SiC fibres: from the grain scale to the macroscopic scale

It is well known that the presence of continuous fibres in SiC/Ti composites leads to a high mechanical anisotropy of the composite between the longitudinal and the two transverse directions. But it is also possible that the crystallographic texture of the matrix may lead to a non-negligible anisotropy of the mechanical properties of the composite. The crystallographic orientation of the matrix grains was determined using the Electron BackScattering Diffraction technique. A local texture of the matrix of the composite is thus evidenced. Finite Element calculations are used to determine the stress field in the matrix resulting from an applied transverse loading. The representative mechanical quantities thus determined are discussed in relation with the fundamental mechanisms of plastic deformation of the matrix. Finally, the crystallographic texture of the matrix of the composite is taken into account in the numerical modellings using a three-scale model that combines crystal plasticity, a polycrystalline aggregate model and a periodic homogenization through a Finite Element unit cell for the composite analysis.     

原位微拉曼测试技术在碳纳米管纤维和薄膜材料力学性能研究中的应用

利用宏观应力联合原位微拉曼测试技术对双壁碳纳米管(Carbon Nanotube,CNT)纤维和薄膜材料的力学性能进行了实验分析,探讨了拉伸加载期间纤维和薄膜内CNT的载荷响应及其与宏观力学性能的关联,揭示了两种材料力学性能差异性的微观机理。实验分析表明,CNT纤维和薄膜的拉伸变形呈现弹性、强化和损伤断裂三个阶段,但其内的CNT只发生弹性变形,没有塑性形变,且没有明显的损伤或键的断裂,纤维和薄膜呈现阶段性拉伸变形的原因可归结为滑移。纤维的弹性模量显著高于薄膜,是薄膜的4.7倍,原因是弹性阶段纤维中CNT的轴向伸长对宏观应变的贡献较大。纤维和薄膜的拉伸强度相差较小,原因是强化阶段薄膜内不断有大量CNT进入承载队伍,这也使得薄膜具有比纤维更高的韧性。     

The effect of brittle interfacial compounds on deformation and fracture of molybdenum-aluminum fiber composites

The effect of brittle intermetallic compounds at the fiber-matrix interface on the deformation characteristics of molybdenum-aluminum fiber composites was investigated. If the filament is ductile and notch-insensitive, then composite strength degradation is relatively minor and composite strength can be predicted by a modified mixture-rule which neglects the strength contribution of the brittle compound. For the case of notch-sensitive filaments, severe filament degradation occurs upon compound formation. The degradation was shown to result from cracks formed during deformation at the roots of compound nodules. The presence of 10 per cent compound by volume results in a 50 per cent decrease in tensile strength, but larger amounts of compound cause little additional strength reduction. At filament volume fractions of 25 and 34 per cent and compound volume fractions less than 10 per cent, composite fracture occurs by the statistical accumulation of fiber necks or fractures depending on the notch sensitivity of the fiber. At high fiber or compound volume fractions, composite failure occurs upon the first or the second filament fracture.     

Prediction of progressive failure in multidirectional composite laminated panels

A mechanism-based progressive failure analyses (PFA) approach is developed for fiber reinforced composite laminates. Each ply of the laminate is modeled as a nonlinear elastic degrading lamina in a state of plane stress according to Schapery theory (ST). In this theory, each lamina degrades as characterized through laboratory scale experiments. In the fiber direction, elastic behavior prevails, however, in the present work, the phenomenon of fiber microbuckling, which is responsible for the sudden degradation of the axial lamina properties under compression, is explicitly accounted for by allowing the fiber rotation at a material point to be a variable in the problem. The latter is motivated by experimental and numerical simulations that show that local fiber rotations in conjunction with a continuously degrading matrix are responsible for the onset of fiber microbuckling leading to kink banding. These features are built into a user defined material subroutine that is implemented through the commercial finite element (FE) software ABAQUS in conjunction with classical lamination theory (CLT) that considers a laminate as a collection of perfectly bonded lamina (Herakovich, C.T., 1998. Mechanics of Fibrous Composites. Wiley, New York). The present model, thus, disbands the notion of a fixed compressive strength, and instead uses the mechanics of the failure process to provide the in situ compression strength of a material point in a lamina, the latter being dictated strongly by the current local stress state, the current state of the lamina transverse material properties and the local fiber rotation. The inputs to the present work are laboratory scale, coupon level test data that provide information on the lamina transverse property degradation (i.e. appropriate, measured, strain–stress relations of the lamina transverse properties), the elastic lamina orthotropic properties, the ultimate tensile strength of the lamina in the fiber direction, the stacking sequence of the laminate and the geometry of the structural panel. The validity of the approach advocated is demonstrated through numerical simulations of the response of two composite structural panels that are loaded to complete failure. A flat, 24-ply unstiffened panel with a cutout subjected to in-plane shear loading, and a double notched 70-ply unstiffened stitched panel subjected to axial compression are selected for study. The predictions of the simulations are compared against experimental data. Good agreement between the present PFA and the experimental data are reported.     

非对称铺层复合材料层合板质量优化设计

复合材料广泛应用于航空航天等领域,追求轻量化设计已经成为研究重点.对复合材料层合板质量优化设计,可以减少层合板的纤维用量,减小层合板的质量,降低成本.首先研究复合材料层合板在承受轴向载荷时,产生的形变量、应力示意图,分析容易发生失效部位;以层合板铺层厚度为设计变量,最大应变、铺层比例等为约束条件,最小化层合板质量为优化...