Analytical study of the load transfer in fibre-reinforced 2D composite materials

The phenomenon of the load diffusion from a fibre to a surrounding matrix is analysed for the 2D case. We use an approximate analytical approach based on the asymptotic reduction of the governing biharmonic problem into two harmonic problems. The comparison of the obtained solutions with known results of other authors shows an acceptable accuracy of the proposed asymptotic simplifications. All solutions are obtained in closed analytical form.The case of perfect bonding between fibre and matrix for a single fibre and for a periodic system of fibres is firstly considered. Then we study the influence of the interface stiffness. In the case when only a single fibre is loaded, the influence of all other fibres is predicted by means of a three-phase model. The proposed approach gives a possibility to solve the problems for a broken fibre and for a broken matrix, as well as for partly debonded fibres. The important problem of infinite matrix cracks is also solved in the present paper.The obtained results can be used for the calculation of pull-out and push-out tests, as well as for the investigation of the fracture of composite materials.     

Anti-plane shear waves in a fibre-reinforced composite with a non-linear imperfect interface

The propagation of non-linear elastic anti-plane shear waves in a unidirectional fibre-reinforced composite material is studied. A model of structural non-linearity is considered, for which the non-linear behaviour of the composite solid is caused by imperfect bonding at the “fibre–matrix” interface. A macroscopic wave equation accounting for the effects of non-linearity and dispersion is derived using the higher-order asymptotic homogenisation method. Explicit analytical solutions for stationary non-linear strain waves are obtained. This type of non-linearity has a crucial influence on the wave propagation mode: for soft non-linearity, localised shock (kink) waves are developed, while for hard non-linearity localised bell-shaped waves appear. Numerical results are presented and the areas of practical applicability of linear and non-linear, long- and short-wave approaches are discussed.     

The deformation of the front of a 3D interface crack propagating quasistatically in a medium with random fracture properties

One studies the evolution in time of the deformation of the front of a semi-infinite 3D interface crack propagating quasistatically in an infinite heterogeneous elastic body. The fracture properties are assumed to be lower on the interface than in the materials so that crack propagation is channelled along the interface, and to vary randomly within the crack plane. The work is based on earlier formulae which provide the first-order change of the stress intensity factors along the front of a semi-infinite interface crack arising from some small but otherwise arbitrary in-plane perturbation of this front. The main object of study is the long-time behavior of various statistical measures of the deformation of the crack front. Special attention is paid to the influences of the mismatch of elastic properties, the type of propagation law (fatigue or brittle fracture) and the stable or unstable character of 2D crack propagation (depending on the loading) upon the development of this deformation.     

Characterization of a liquid-metal microdroplet thermal interface material

This work presents the characterization of a thermal interface material consisting of an array of mercury microdroplets deposited on a silicon die. Three arrays were tested, a 40 × 40 array (1600 grid) and two 20 × 20 arrays (400 grid). All arrays were assembled on a 4 × 4 mm² silicon die. An experimental facility which measures the thermal resistance across the mercury array under steady state conditions is described. The thermal interface resistance of the arrays was characterized as a function of the applied load. A thermal interface resistance as low as 0.253 mm² K W⁻¹ was measured. A model to predict the thermal resistance of a liquid-metal microdroplet array was developed and compared to the experimental results. The contact resistance of the mercury arrays was estimated based on the experimental and model data. An average contact resistance was estimated to be 0.14 mm² K W⁻¹.     

Experimental investigation of high temperature thermal contact resistance with interface material

Thermal contact resistance plays a very important role in heat transfer efficiency and thermomechanical coupling response between two materials, and a common method to reduce the thermal contact resistance is to fill a soft interface material between these two materials. A testing system of high temperature thermal contact resistance based on INSTRON 8874 is established in the present paper, which can achieve 600°C at the interface. Based on this system, the thermal contact resistance between superalloy GH600 material and three-dimensional braid C/C composite material is experimentally investigated, under different interface pressures, interface roughnesses and temperatures, respectively. At the same time, the mechanism of reducing the thermal contact resistance with carbon fiber sheet as interface material is experimentally investigated. Results show that the present testing system is feasible in the experimental research of high temperature thermal contact resistance.     

Adaptive subdivision piecewise linear interface calculation (ASPLIC) for 2D multi‐material hydrodynamic simulation codes

In this work, we propose an adaptive subdivision piecewise linear interface calculation (PLIC) for 2D multimaterial hydrodynamic simulation codes. Classical volume‐of‐fluid PLIC technique uses one line segment and one given normal to separate two materials. Unfortunately, these paradigms are not sufficient when filaments occur, leading to the creation of flotsam and jetsam. We propose to detect such situations and to split the computational mixed cell into reconstruction subzones. Within these subzones, one computes a so‐called subgradient using an incomplete stencil of neighbors, and the material is distributed in these subzones. Given subzone volume fraction and the subgradient, one computes one line segment using classical PLIC method, leading to a modified PLIC method for subscale material entity. The subdivision procedure relies on a splitting point, which is chosen as a specific information about the relative location of the filament in the cell, leading to an adaptive subdivision for PLIC reconstructions. Numerical tests are carried out in a 2D Lagrange + Remap multimaterial hydrodynamics Eulerian code. Static and dynamic filaments and fragments are simulated in advection or stretched in vortex‐like motion. The full hydrodynamics equations are solved on a more realistic test (shock‐bubble impact). Results show that our approach supplements classical PLIC method for situations when filaments and fragments occur. Copyright © 2014 John Wiley & Sons, Ltd.     

A study ofJ-integral of the orthotropic composite material

The relation between J-integral near model I crack tip in the orthotropic plate and displacement derivative is derived in this paper. Meanwhile, the relation between stress intensity factor K _I and displacement is also given in this paper. With sticking film moire interferometry method, the three-point bending beam is tested, thus the values of J-integral and K _I can be obtained from the displacement field, and then the truth of relation formula between J-integral and K _I in the orthotropic composite materials is experimentally verified.     

On the composite Riemann problem for multi‐material fluid flows

We develop an Eulerian fixed grid numerical method for calculating multi‐material fluid flows. This approach relates to the class of interface capturing methods. The fluid is treated as a heterogeneous mixture of constituent materials, and the material interface is implicitly captured by a region of mixed cells that have arisen owing to numerical diffusion. To suppress this numerical diffusion, we propose a composite Riemann problem (CRP), which describes the decay of an initial discontinuity in the presence of a contact point between two different fluids, which is located off the initial discontinuity point. The solution to the CRP serves to calculate multi‐material no mixed numerical flux without introducing any material diffusion. We discuss the CRP solution and its implementation in the multi‐material fluid Godunov method. Numerical results show that a simple framework of the CRP greatly improves capturing material interfaces in the Godunov method and reproduces many of the advantages of more complicated interface tracking multi‐material treatments. Copyright © 2014 John Wiley & Sons, Ltd.     

ANALYSIS OF INTERPHASE MECHANICAL BEHAVIOR WITH INTERFACE ELEMENT IN THE COMPOSITE MATERIALS

ＡＮＡＬＹＳＩＳＯＦＩＮＴＥＲＰＨＡＳＥＭＥＣＨＡＮＩＣＡＬＢＥＨＡＶＩＯＲＷＩＴＨＩＮＴＥＲＦＡＣＥＥＬＥＭＥＮＴＩＮＴＨＥＣＯＭＰＯＳＩＴＥＭＡＴＥＲＩＡＬＳＡＮＡＬＹＳＩＳＯＦＩＮＴＥＲＰＨＡＳＥＭＥＣＨＡＮＩＣＡＬＢＥＨＡＶＩＯＲＷＩＴＨＩＮ...     

Dynamic stress analysis of the interface in a particle-reinforced composite

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｉｎｔｅｒｆａｃｉａｌｄｅｂｏｎｄｉｎｇｍａｙｃａｕｓｅｔｈｅｍｉｃｒｏｖｏｉｄｓ’ｎｕｃｌｅａｔｉｏｎｉｎａｐａｒｔｉｃｌｅ＿ｒｅｉｎｆｏｒｃｅｄｃｏｍｐｏｓｉｔｅ .Ｔｈｅｉｎｔｅｒｆａｃｉａｌｄｅｂｏｎｄｉｎｇｉｓｕｓｕａｌｌｙｇｏｖｅｒｎｅｄｂｙｔｈｅｔｅｎｓｉｏｎｓｔｒｅｓｓｃｒｉｔｅｒｉｏｎ[1,2 ]ｏｒｅｎｅｒｇｙｃｒｉｔｅｒｉｏｎ[3,4].Ｉｎｔｈｅｒｅｃｅｎｔｒｅｓｅａｒｃｈｗｏｒｋｓ,ｔｈｅｐｒｅｓｅｎｔａｕｔｈｏｒｓｐｏｉｎｔｅｄｏｕｔｔｈａｔｔｈｅａｎ…     

复合材料和粘弹性材料中的应力奇异性

在各种材料中,出现的应力奇异性已为工程技术人员的科学研究提供了丰富的素材,本文的主要目的在于对粘弹性材料的应力奇异性进行探讨,利用对应性原理从弹性材料导出粘生材料的应力奇异分布。     

压电复合材料粘接界面断裂有限元模拟

根据数字化FRMM(Fix-Ratio Mix-Mode)断裂试验,得到了压电复合材料试件的断裂韧性和位移及应变场。本文在试验的基础上,通过非线性有限元软件ABAQUS及用户子程序UMAT进行了模拟分析,采用基于损伤力学的粘聚区模型(CZM)对压电复合材料界面的起裂和脱胶扩展进行了分析,并与VCCT方法进行了比较。计算得到的荷载位移曲线更接近于试验结果,但在裂纹扩展路径上的吻合需要对粘聚区法则进一步修正。通过进一步对CZM参数进行分析,表明界面粘结强度和界面刚度对计算结果的影响很大。研究结果表明,粘聚区模型可以很好地表征压电复合材料弱粘接界面脱胶断裂问题。     

Crack in a 2D beam lattice: Analytical solutions for two bending modes

We consider an infinite square-cell lattice of elastic beams with a semi-infinite crack. Symmetric and antisymmetric bending modes of fracture under remote loads are examined. The related long-wave asymptotes corresponding to a continuous anisotropic bending plate are also considered. In the latter model, the symmetric mode is characterized by the square-root type singularity, whereas the antisymmetric mode results in a hyper-singular field. A solution for the continuous plate with a finite crack is also presented. These closed-form continuous solutions describe the fields in the whole plane. The main goal is to establish analytical connections between the ‘macrolevel’ state, defined by the continuous asymptote of the lattice solution, and the maximal bending moment in the crack-front beam, that is, to determine the resistance of the lattice with an initial crack to the crack advance. The solutions are obtained in the same way as for mass-spring lattices. Considering the static problems we use the discrete Fourier transform and the Wiener-Hopf technique. Monotonically distributed bending moments ahead of the crack are determined for the symmetric mode, and a self-equilibrated transverse force distribution is found for the antisymmetric mode. It is shown that in the latter case only the crack-front beam resists to the fracture development, whereas the forces in the other beams facilitate the fracture. In this way, the macrolevel fracture energy is determined in terms of the material strength. The macrolevel energy release is found to be much greater than the critical strain energy of the beam, especially in the hyper-singular mode. In both problems, it is found that among the beams surrounding the crack the crack-front beam is maximally stressed, and hence its strength defines the strength of the structure.     

In-plane displacement measurement of composite material by objective speckle

The objective speckle was used to measure the in-plane displacement around the hole in a coupon specimen of composite material under tensile loading. Mirror transplantation method was adopted to ensure the high reflectivity of specimen surface and to obtain high quality double exposure specklegram. The adjustable spatial frequency of the whole-field pattern of the displacement made it possible to measure a broad range of strain, from elastic to plastic. The results obtained by objective speckle agree well with those by Moire method.     

复合材料切口应力奇性指数计算

提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段.     

Comparative study of an interface crack for different wedge-interface models

Summary  An interface crack problem is investigated under various assumptions on an interface between two elastic materials. The interface is modeled by an additional third structure (thin elastic wedge of differing elastic properties) matching the bonded materials, or by introducing special boundary conditions on the crack line ahead. The main emphasis of the paper is placed on a comparison of the asymptotic expansion of the elastic solutions near the crack tip obtained for the different models. In particular, the behaviour of the stress singularity exponent and the generalized SIF are discussed. Numerical examples are presented. Received 16 August 2000; accepted for publication 26 May 2001     

Structural optimization and experimental analysis of composite material panels for naval use

The purpose of the work is the realization of a composite material with long glass fibers having better characteristic than a fiber random composite, to permit the reduction of weight and costs to shipyards for pleasure craft. Structural optimization is performed by ANSYS for the choice of the layers disposition to obtain the maximum stiffness with minimum material employment, saving weight. The study is centered on the research of the better configurations of plies packing in relation of pure shear stress for four different plies. Unidirectional plies, both symmetric orthotropic and symmetric non-orthotropic ones, are realized successively by the vacuum bag technique. Experimental tests of traction, bending, inter laminar shear and pure shear are executed to characterize the three different type of material. Experimental results are compared to ones obtained numerically to validate the procedures; the comparison with the analytical results permitted to attribute an adequate value to shape factor of the fibers. In all the cases the optimization permitted the construction of much more resistant plies than random ones, with a lower thickness.     

Case studies on the elastic-plastic creep in the metallic matrix of fibre-reinforced composite materials

Summary FEM-microanalyses for cyclic thermomechanical loads are performed on fibre-reinforced metallic matrix composites in order to study the damage history. First, material modelling is carried out in which, in addition to the elastic deformation, inelastic deformations (plasticity and creep) are taken into account. After the implementation of the model in the FEM, numerical shakedown simulations are performed in case studies.     

An experimental study of spallation criterion in phenolic-resin based woven roving glass fibre reinforced composite material

An experimental study of spallation in phenolic-resin based woven roving glass fibre reinforced composite material under uniaxial strain condition is reported in this paper. Experiments were carried out with a 101 mm bore single-stage light gas gun. The stress waves induced by plate impact propagate along the normal to glass fibre mats. The experiments show that the threshold of incipient spallation in this material depends strongly on the duration of stress wave and that this incipient spallation appears essentially in one layer. The criterion of spallation in this material given in this paper is in good agreement with experimental results. Besides, the spallation in phenolic-resin based obliquely laminated woven roving glass fibre reinforced composite material is also discussed in this paper.