The issues of the uniqueness and the stability of the homogeneous response in uniaxial tests with gradient damage models

We consider a wide class of gradient damage models which are characterized by two constitutive functions after a normalization of the scalar damage parameter. The evolution problem is formulated following a variational approach based on the principles of irreversibility, stability and energy balance. Applied to a monotonically increasing traction test of a one-dimensional bar, we consider the homogeneous response where both the strain and the damage fields are uniform in space. In the case of a softening behavior, we show that the homogeneous state of the bar at a given time is stable provided that the length of the bar is less than a state dependent critical value and unstable otherwise. However, we also show that bifurcations can appear even if the homogeneous state is stable. All these results are obtained in a closed form. Finally, we propose a practical method to identify the two constitutive functions. This method is based on the measure of the homogeneous response in a situation where this response is stable without possibility of bifurcation, and on a procedure which gives the opportunity to detect its loss of stability. All the theoretical analyses are illustrated by examples.     

On the influence of interfacial properties to the bending rigidity of layered structures

Layered structures are ubiquitous, from one-atom thick layers in two-dimensional materials, to nanoscale lipid bi-layers, and to micro and millimeter thick layers in composites. The mechanical behavior of layered structures heavily depends on the interfacial properties and is of great interest in engineering practice. In this work, we give an analytical solution of the bending rigidity of bilayered structures as a function of the interfacial shear strength. Our results show that while the critical bending stiffness when the interface starts to slide plastically is proportional to the interfacial shear strength, there is a strong nonlinearity between the rigidity and the applied bending after interfacial plastic shearing. We further give semi-analytical solutions to the bending of bilayers when both interfacial shearing and pre-existing crack are present in the interface of rectangular and circular bilayers. The analytical solutions are validated by using finite element simulations. Our analysis suggests that interfacial shearing resistance, interfacial stiffness and preexisting cracks dramatically influence the bending rigidity of bilayers. The results can be utilized to understand the significant stiffness difference in typical biostructures and novel materials, and may also be used for non-destructive detection of interfacial crack in composites when stiffness can be probed through vibration techniques.     

On the existence of the low-frequency surface waves in a porous mediumSur l'existence des ondes de surface basse fréquence en milieux poreux

The existence and propagation of the surface waves at a vacuum/porous medium interface are investigated in the low frequency range. Two types of surface waves are shown to be possible: the generalized Rayleigh wave, which always exists, and the Stoneley wave, which exists for a limited range of wave numbers. Moreover, within the k-domain of existence the Stoneley wave cannot appear for certain values of elastic parameters of the solid phase. The bifurcation behavior of both the Stoneley wave and the Biot (P2) bulk wave, depending on the wave number, is revealed. The asymptotic formulas for the phase velocities of the surface waves are derived. To cite this article: I. Edelman, C. R. Mecanique 332 (2004).     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

Positivité de la dissipation intrinsèque d'une classe de modèles d'endommagement anisotropes non standardsPositivity of intrinsic dissipation of a class of nonstandard anisotropic damage models

A nonstandard thermodynamics framework ensures the positivity of the dissipation due to degradation mechanisms for damage states represented by a symmetric second order tensor. The proof of the positivity of the intrinsic dissipation is given. An increasing damage, in terms of positive principal values of the damage rate tensor, guaranties this positivity for the considered class of models, extending then to induced anisotropy the isotropic case property of a positive damage rate. This result gives many possibilities of modeling for anisotropic damage. To cite this article: R. Desmorat, C. R. Mecanique 334 (2006).     

On localization in saturated porous continua

Localization in elastic-plastic saturated porous media is investigated here using a linear perturbation approach. The adopted localization criterion corresponds to unbounded growth of perturbations. The critical conditions are compared with those obtained by a classical band analysis. While for one phase materials these conditions coincide, in the present context the linear perturbation approach leads in the limit of unbounded growth to the singularity of the undrained acoustic tensor, while the band analysis leads to the singularity of the drained acoustic tensor. Some general results clarifying the hierarchy of these two conditions are established for a quite general set of constitutive equations.     

Inverse analysis determining interfacial properties between metal film and ceramic substrate with an adhesive layer

In the present study, peel tests and inverse analysis were performed to determine the interracial mechanical parameters for the metal film/ceramic system with an epoxy interface layer between film and ceramic. Al films with a series of thicknesses between 20 and 250 μm and three peel angles of 90°, 135° and 180° were considered. A finite element model with the cohesive zone elements was used to simulate the peeling process. The finite element results were taken as the training data of a neural network in the inverse analysis. The interracial cohesive energy and the separation strength can be determined based on the inverse analysis and peel experimental result     

On variational approaches in NRT continua

In the paper some features of the theory of Not Resisting Tension (NRT) material are deepened. In details, one first introduces the basic NRT model, which is proved to simply and effectively interpreting the behaviour of mechanical bodies made by not-cohesive materials; thereafter one analyses energetic approaches and limit analysis tools for problems relevant to NRT continua. Afterward, on the basis of the fundamental variational theorems, the main rules governing the NRT behaviour are demonstrated, by imposing Kuhn–Tucker stationarity conditions for the stated constrained optimisation procedures. Finally an application is operated of the presented theory to an elastic NRT semi-plane subject to a distributed load, reproducing the stress situation induced in the soil by a foundation structure.     

Mo-ZrC梯度金属陶瓷的冲击响应行为

分层梯度材料特定的梯度变化能有效增强材料性能。为研究梯度结构、冲击方向对分层梯度材料冲击响应的影响,利用分离式霍普金森压杆结合高速摄影技术对Mo-ZrC分层梯度金属陶瓷进行了动态压缩实验,基于数字图像相关技术讨论了梯度结构、冲击方向对金属陶瓷材料破坏模式的影响,利用Mori-Tanaka理论计算得到金属陶瓷等效性质,结合应力波理论研究波在分层梯度复合材料中的传播规律。结果表明：(1)相同加载条件下,梯度结构对材料的强度、韧性和破坏产物的完整性具有重要影响,在冲击过程中,样品响应可以分为压紧阶段、裂纹成核发展阶段和贯穿阶段,对于不同梯度结构和冲击方向,样品在加载过程中呈现出不同的破坏时序和失效模式;(2)利用数字图像相关方法跟踪分层梯度陶瓷的局部变形发展,分析发现局部增量达到临界状态后,局部变形发展转化为微裂纹的形成和累积,最终导致整体性破碎失效;(3)通过分层梯度材料一维应力波传播理论推导得到,改变冲击梯度方向对应力波透反射系数存在一定影响,不同梯度结构设计对改变冲击梯度方向敏感性不同,且存在极值情况。     

On the non-conservativeness of a class of anisotropic damage models with unilateral effectsSur le caractère non-conservatif de quelques modèles d'endommagement anisotrope avec conditions unilatérales

The aim of this Note is to show that a class of anisotropic elastic-damage models including unilateral effects can be considered, for constant damage values, as non-linear and non-conservative elastic. The conservative character of corresponding constitutive models is related to the symmetry of the Hessian tensor. For the models under consideration, it is shown that the condition of conservativeness (existence of the elastic potential energy function) is obtained only when there is coaxiality of the strain and damage tensors. To cite this article: N. Challamel et al., C. R. Mecanique 334 (2006).     

Interaction of an anti-plane singularity with interfacial anti-cracks in cylindrically anisotropic composites

Anti-plane problem for a singularity interacting with interfacial anti-cracks (rigid lines) under uniform shear stress at infinity in cylindrically anisotropic composites is investigated by utilizing a complex potential technique in this paper. After obtaining the general solution for this problem, the closed solution for the interface containing one anti-crack is presented analytically. In addition, the complex potentials for a screw dislocation dipole inside matrix are obtained by the superimposing method. Expressions of stress singularities around the anti-crack tips, image forces and torques acting on the dislocation or the center of dipole are given explicitly. The results indicate that the anisotropy properties of materials may weaken the stress singularity near the anti-crack tip for the singularity being a concentrated force but enhance the one for the singularity being a screw dislocation and change the equilibrium position of screw dislocation. The presented solutions are valid for anisotropic, orthotropic or isotropic composites and can be reduced to some new or previously known results.     

Strain energy release rate for interfacial cracks in hybrid beams

The finite element modeling and fracture mechanics concept were used to study the interfacial fracture of a FRP-concrete hybrid structure. The strain energy release rate of the interfacial crack was calculated by the virtual crack extension method. It is shown that the crack growth has three phases, namely, cracking initiation, stable crack growth and unstable crack propagation. The effects of geometric and physical parameters of the hybrid beam on the energy release rate were considered. These parameters include Young’s moduli of the FRP, the concrete and the adhesive, thickness of the FRP plate and adhesive, and the distance of FRP plate end from the beam end. The numerical results show that the energy release rate of the interfacial crack is influenced considerably by these parameters. The present investigation can contribute to the mechanism understanding and engineering design of the hybrid structures.     

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and micro-devices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either ‘microhard’ (impenetrable to dislocations) or ‘microfree’ (an infinite dislocation sink).     

Flow patterns and interfacial velocities near a moving contact line

Experimental data on velocity fields and flow patterns near a moving contact line is shown to be at variance with existing hydrodynamic theories. The discrepancy points to a new hydrodynamic paradox and suggests that the hydrodynamic approach may be incomplete and further parameters or forces affecting the surfaces may have to be included. A contact line is the line of intersection of three phases: (1) a solid, (2) a liquid, and (3) a fluid (liquid or gas) phase. A moving contact line develops when the contact line moves along the solid surface. A flat plate moved up and down, inside and out of a liquid pool defines a simple, reliable experimental model to characterize dynamic contact lines. Highlighted are three important conclusions from the experimental results that should be prominent in the development of new theoretical models for this flow. First, the velocity along the streamline configuring the liquid–fluid interface is remarkably constant within a distance of a couple of millimeters from the contact line. Second, the relative velocity of the liquid–fluid interface, defined as the ratio of the velocity along the interface to the velocity of the solid surface, is independent of the solid surface velocity. Third, the relative interface velocity is a function of the dynamic contact angle.     

On moving continua with contacts and sliding friction: Modeling,general properties and examples

Many technical applications feature contacts between moving continua, where the contact zone is more or less spatially fixed. This contribution describes the modeling of such systems by means of weak formulations derived from Hamilton’s Principle in a material as well as a spatial formulation. Furthermore, linearizations about a stationary solution are carried out, the arising differential operators are stated and discussed with respect to symmetry. Particular attention is paid to the structure and physical interpretation of the contributions arising from the contacts as well as those due to the transport motion. The article closes with several examples.     

Subsonic interfacial fracture using strain gages in isotropic–orthotropic bimaterial

An experimental study has been conducted in which strain fields were used to investigate the behavior of subsonic crack propagation along the interface of an isotropic–orthotropic bimaterial system. Strain field equations were developed from available field equations and critically evaluated in a parametric study to identify optimum strain gage location and orientation. Bimaterial specimens were prepared with PSM-1 polycarbonate and Scotchply^® 1002 unidirectional, glass-fiber-reinforced, epoxy composite. Dynamic experiments were conducted using these specimens with strain gages mounted on the composite half to obtain values of the dynamic complex stress intensity factor, K=K₁+iK₂, in the region of the crack tip while photoelasticity was used on the PSM-1 half. Results show that the trend and magnitude of K obtained using strain gages compare favorably with those obtained using photoelasticity.     

Approche micromécanique du couplage perméabilité–endommagementMicromechanical approach to the coupling between permeability and damage

A self-consistent scheme is used in order to determine the permeability of a cracked porous medium. For weak values of the permeability of the uncracked porous matrix, it is found that the order of magnitude of the permeability increases beyond a critical threshold of the crack density parameter. In the framework of a micromechanical model, it is shown that both the evolution of crack opening and the crack propagation are controlled by Terzaghi's effective stress which therefore captures the coupling between permeability and mechanical loading. To cite this article: L. Dormieux, D. Kondo, C. R. Mecanique 332 (2004).     

Time and spatially resolved measurements of interfacial waves in vertical annular flow

Film thickness measurements have been performed in a vertical air/water annular flow in a pipe of 0.05 m diameter. A sensor has been built which allows to measure the film thickness evolution in time at 320 positions, such that the interface of the vertical annular flow can be reconstructed. The large-scale structures moving on the interface are described statistically, with a special attention to the disturbance waves. Probability density functions and mean statistics are given for the height, length, velocity, frequency and spatial distribution of the disturbance waves. In particular, it is shown that the disturbance waves are three-dimensional structures with large height fluctuations in the circumferential and axial direction, giving a meandering path between the maximum height around the circumference. It is also shown that the disturbance waves can flow with a slight inclination with respect to the axial direction. Finally, the disturbance waves are shown to be located randomly in space, within a Gamma distribution whose order only depends on the liquid superficial velocity. Due to the nature of the Gamma distribution, it could indicate that the spatial distribution of the disturbance waves results from a cascade of coalescence processes between the original disturbance waves on the film.