黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

On the Torsion of a Class of Nonlinear Fiber Composite Cylinders

This paper presents an analysis of the torsion of a solid or annular circular cylinder consisting of nonlinear material in the form of an elastic matrix with embedded unidirectional elastic fibers parallel to the cylinder axis. The specific class of composite considered is one for which nonlinear fiber-matrix interface slip is captured by uniform cohesive zones of vanishing thickness. Previous work on the effective antiplane shear response of this material leads to a stress–strain relation depending on the interface slip together with an integral equation governing its evolution. Here, we obtain an approximate single mode solution to the integral equation and utilize it to solve the torsion problem. Equations governing the radial distributions of shear stress and interface slip are obtained and formulae for torque–twist rate are presented. The existence of singular surfaces, i.e., surfaces across which the slip and the shear stress experience jump discontinuities are analyzed in detail. Specific results are presented for an interface force law that allows for interface failure in shear.     

Stress analysis of a layered elastic solid in contact with a rough surface exhibiting fractal behavior

A contact stress analysis is presented for a layered elastic half-space in contact with a rough surface exhibiting self-affine (fractal) behavior. Relationships for the mean contact pressure versus representative strain and the real half-contact width versus elastic properties of the layer and the substrate, asperity radius, layer thickness, and truncated half-contact width were derived from finite element simulations of a layered medium compressed elastically by a rigid cylindrical asperity. These relationships were incorporated in a numerical algorithm that was used to obtain the contact pressure distributions and stresses generated by the asperity contacts formed at the interface of the layered medium and the fractal surface. Analytical solutions illustrate the significance of the elastic material properties, layer thickness, and surface topography (roughness) on global parameters such as normal load and real contact area. Results for the contact pressure distribution and the surface and subsurface stresses provide insight into the initiation of yielding and the tendency for cracking in the layered medium. It is shown that cracking at the surface and the layer/substrate interface is more likely to occur in the case of a stiff layer, whereas surface cracking is more prominent for a relatively compliant layer.     

Interfacial slip and creep in rolling contact incorporating a cylinder with an elastic layer

A systematic approach for investigating the interfacial behaviour of tyred systems is proposed. A two-dimensional contact model of an elastic strip, shrink-fitted onto a wheel, and subjected to different rolling contact conditions, has been adopted to illustrate the method. The model combines existing techniques to explore individual elastic contact problems and it enables us to characterise the behaviour at the strip/substrate interface caused by loads induced by a quasi-static application of stationary and moving loads on the surface of the layer. The solution is compared to the stationary load case and regimes of local slip, full stick, separation and frictional creep are identified and collated for a variety of loading conditions, materials and geometries. Further, this article presents an investigation of frictional shakedown for layered systems subjected to periodic contact loading. The term shakedown is here referred to as the possibility of developing interfacial residual stresses at the layer/substrate interface such that frictional slip, originally activated by the applied external contact load, ceases after a few loading cycles. The possible applicability of the Melan’s theorem for elastic frictional system is investigated and preliminary results presented.     

Wave induced sliding at the interface between a layered elastic medium and a half-space

The sliding interface between an unrestrained elastic half-space and a grounded layered half-space excited by an incident harmonic wave is investigated. The problem is formulated considering various possible boundary conditions and boundary inequalities at the sliding interface. The Coulomb friction model without distinction between the static and kinetic coefficients of friction is considered to govern the sliding condition. Three possible bands at the interface, namely slip, stick, and separation, are considered. The interface is assumed to be preloaded under normal and shear stresses. The solution is developed by modifying the problem of welded interface, which then is reduced to a set of algebraic equations. The effects of the incident angle, layer thickness, friction coefficient and externally applied stresses on the drift velocity of the unrestrained half-space are studied numerically for a pair of materials. It is shown that the sliding interface, and hence the drift velocity of unrestrained half-space is noticeably influenced by the layered medium. These results are expected to be useful for the development of a new kind of ultrasonic drive in future.     

Dynamic thermal shock in a layered cylinder with initial interface pressure

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｙｎａｍｉｃｐｒｏｂｌｅｍｏｆｓｔｒｕｃｔｕｒｅｓｔｒｅｎｇｔｈｉｓａｎｏｕｔｓｔａｎｄｉｎｇｐｒｏｂｌｅｍ．Ａｆｔｅｒｃｏｎｓｕｌｔｉｎｇｓｏｍｅｐａｐｅｒｓ,ｗｅｃａｎｋｎｏｗｔｈａｔｔｈｉｓｓｕｂｊｅｃｔｈａｓｂｅｅｎｓｔｕｄｉｅｄｅｘｔｅｎｓｉｖｅｌｙｂｙｓｅｖｅｒａｌａｕｔｈｏｒｓ［１～７］．Ｔｈｅｉｒｉｎｖｅｓｔｉｇａｔｉｏｎｓｏｎｌｙｂａｓｅｄｏｎｃａｌｃｕｌａｔｉｎｇｖｉｂｒａｔｉｏｎｉｎａｌａｙｅｒｅｄｓｔｒｕｃｔｕｒｅａｎｄｐｒｅｓｅｎｔｉｎｇ…     

Discrete dislocation dynamics modelling of mechanical deformation of nickel-based single crystal superalloys

Discrete dislocation dynamics (DDD) has been used to model the deformation of nickel-based single crystal superalloys with a high volume fraction of precipitates at high temperature. A representative volume cell (RVC), comprising of both the precipitate and the matrix phase, was employed in the simulation where a periodic boundary condition was applied. The dislocation Frank-Read sources were randomly assigned with an initial density on the 12 octahedral slip systems in the matrix channel. Precipitate shearing by superdislocations was modelled using a back force model, and the coherency stress was considered by applying an initial internal stress field. Strain-controlled loading was applied to the RVC in the [0 0 1] direction. In addition to dislocation structure and density evolution, global stress-strain responses were also modelled considering the influence of precipitate shearing, precipitate morphology, internal microstructure scale (channel width and precipitate size) and coherency stress. A three-stage stress-strain response observed in the experiments was modelled when precipitate shearing by superdislocations was considered. The polarised dislocation structure deposited on the precipitate/matrix interface was found to be the dominant strain hardening mechanism. Internal microstructure size, precipitate shape and arrangement can significantly affect the deformation of the single crystal superalloy by changing the constraint effect and dislocation mobility. The coherency stress field has a negligible influence on the stress-strain response, at least for cuboidal precipitates considered in the simulation. Preliminary work was also carried out to simulate the cyclic deformation in a single crystal Ni-based superalloy using the DDD model, although no cyclic hardening or softening was captured due to the lack of precipitate shearing and dislocation cross slip for the applied strain considered.     

Dynamics of Phase Transitions and Hysteresis in a Viscoelastic Ericksen's Bar on an Elastic Foundation

This work is a follow-up on the study [32] of interface dynamics and hysteresis in materials undergoing solid-solid phase transitions. We consider the dynamics of a viscoelastic bar with a nonmonotone stress-strain relation and viscous stress linearly proportional to the strain rate. The bar is placed on an elastic foundation with stiffness β mimicking the interaction of phases in higher dimensions. Time-dependent displacement-controlled loading of the bar results in a tilted and serrated hysteresis loop, in qualitative agreement with some experimental observations in shape-memory alloys. The model exhibits three phase transition processes: phase nucleation, interface slip and phase annihilation. Between these dynamic processes the system gets stuck in local minimizers of the potential energy. As β increases from zero, a slip-dominated solution behavior transforms to the one where slip and annihilation events are preceded by a step-by-step nucleation process. We show that this transition is caused by an interplay between the slip-favoring inertia term and the nucleation-favoring elastic foundation terms. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Force on a Dislocation Near a Weakly Bonded Interface

The physical system being studied here is an elastic dislocation in proximity to a material interface which has limited strength in shear. The configuration is relevant to elastic strain relaxation in bonded thin-film semiconductor material structures. A plane strain boundary value problem is formulated and an exact solution for the force on the dislocation versus distance from the interface is obtained. If the dislocation is close enough to the interface, its stress field causes irreversible slip across the interface. This slip, in turn, induces an attractive force on the dislocation. For the case of an edge dislocation with Burgers vector normal to the interface, it is found that the dislocation has a stable equilibrium position at a small distance from the interface. This revised version was published online in August 2006 with corrections to the Cover Date.     

Viscoelastic layer under the action of a moving load

Recently, problems concerning the dynamic behavior of imperfect continuous media under various types of actions have been widely investigated. The method of Laplace transformation is very convenient for describing physical processes concerning unsteady phenomena. In viscoelastic media two complications are added: the representation of the properties of a medium depending on time, and the inversion of the obtained solutions containing this additional complication. Certain approximate methods of inversion in the analysis of viscoelastic stresses are discussed in [1]. In [2, 3] a discussion is given for an effective method of constructing the solution of unsteady problems for finite and for infinite imperfect media using auxiliary functions, and a solution is presented for a half-space. Making use of the idea of the inversion of transforms, discussed in [4], in [5] a solution is obtained and a complete picture is presented for the dynamics of the variation of the stress field in a viscoelastic half-space. In the present study we consider the action of a normal moving load that is suddenly applied to the free surface of a viscoelastic layer. By Laplace and Fourier integral transformations we obtain a solution in the form of a uniformly converging series based on longitudinal and transverse waves reflected in the layer. By means of inverting the transforms by the method discussed in [4, 5], we obtain an exact solution for the stress field in the medium under investigation. We consider the special case of a viscoelastic medium of Boltzmann type, for which we obtain a numerical realization of the solution on a digital computer.     

土遗址楠竹锚固界面力学行为与传力机理研究

楠竹锚固技术是处理纵向裂隙引起的大型土遗址稳定问题的重要手段,但此类锚固系统传力机理尚不明确,成为制约其科学化、规模化应用的瓶颈。通过现场拉拔试验研究锚固界面力学行为,将识别出的完全脱粘现象通过摩擦段后的零剪应力段进行表征,建立适用于此类锚固界面的修正三线型粘结-滑移模型。在此基础上,提出基于ANSYS中非线性弹簧单元的锚固界面力学行为数值模拟方法,系统研究锚固系统传力机理。最后,通过与试验结果的对比验证了该模拟方法的可靠性。结果表明：锚杆/浆体界面的滑移失效是锚固系统破坏的主要原因,此类界面微段的受力过程可分为弹性、软化、摩擦和完全脱粘四个阶段,随荷载增加锚固界面高应力区逐渐由加载端向锚固端转移,界面进入完全脱粘阶段后剪应力趋近于零,极限锚固力和有效锚固长度均存在临界值;第一锚固界面传递至第二锚固界面的剪应力非常有限,遗址土体应力处于较低水平。研究结果能够为基于“安全第一,最小干预”原则的土遗址文物锚固优化设计提供参考。     

Bending and buckling of a class of nonlinear fiber composite rods

An investigation of the mechanics of bending and buckling is carried out for a class of nonlinear fiber composite rods composed of embedded unidirectional fibers parallel to the rod axis. The specific class of composite considered is one in which the fibers interact with the matrix through a nonlinear Needleman-type cohesive zone [Needleman, A., 1987. A continuum model for void nucleation by inclusion debonding. ASME J. Appl. Mech. 54, 525-531; Needleman, A., 1992. Micromechanical modelling of interfacial decohesion. Ultramicroscopy 40, 203-214]. The primary decohesive mechanism active in bending and buckling of these composite rods is shear slip along the fiber-matrix interfaces allowing the use of a previously developed constitutive relation for antiplane shear response [Levy, A.J., 2000b. The fiber composite with nonlinear interface—part II: antiplane shear. ASME J. Appl. Mech. 67, 733-739]. The formulation requires the specification of a potential interface force-slip law that is assumed to permit interface failure in shear.Four cases of the bending and shearing of beams (concentrated or uniform load on a cantilever or a simply supported beam) are analyzed, each of which exhibits qualitatively distinct response. For certain values of interface parameters, the beam deflection or its gradient at a fixed location can change discontinuously with load. Furthermore, for interface parameter values within a certain range, singular surfaces will exist in uniformly loaded beams where there is a non-uniform distribution of shear stress along the beam length. These singular surfaces divide the beam into regions of maximal and minimal fiber slip and propagate with a rate that varies inversely as the square of the applied load. For other parameter values, singular surfaces will not exist and fiber slip will be diffuse.For the class of nonlinear composite considered, bifurcation and imperfection buckling of pinned-pinned columns is analyzed. For bifurcation buckling, a nonlinear eigenvalue problem is derived and the solution is obtained by Galerkin's method. It is demonstrated that critical loads are influenced by the initial slope, and hence the linear portion, of the interface force-slip relation but the post-buckling response, which in some sense resembles that of plastic buckling, is affected by the entire interface constitutive relation. Imperfection buckling is analyzed in a similar manner by assuming a slight initial curvature of the rod. Sensitivity of the response to imperfection magnitude is discussed as well.     

APPROXIMATE THEORY AND ANALYTICAL SOLUTION FOR DYNAMIC CALCULATION OF VISCOELASTIC LAYERED PERIODIC PLATE

由于周期性隔振结构动力计算中较少考虑轨道交通载荷及材料黏弹性,因此,本文以黏弹性层状周期板为研究对象,提出了垂向移动简谐载荷下,可以考虑材料黏弹性及板内横向剪切变形的黏弹性层状周期板动力计算近似理论并给出解析解答.设板中性面的横向剪切变形为横截面的整体剪切变形,利用Reissner-Mindlin假设及提出的剪切变形补充计算条件,得到了中性面法线转角与中性面剪应力的关系.基于平衡方程和应力连续条件,建立了黏弹性层状周期板振动控制方程,推导了对边简支对边自由条件下,板垂向位移的简化Fourier级数形式解.与经典层合板模型和有限元计算结果进行了比较,验证了本文解答的有效性.结果表明：（1）黏弹性层状周期板可以显著降低单一材料板在自振频率处的振动响应,但会引起局部低频频段的振动放大;（2）板的垂向位移随着载荷速度的增大而增大,当载荷速度超过300 km/h后,其对板振动响应的影响减弱;（3）黏弹性层剪切模量存在最佳设计值,可使结构的隔振性能最佳;（4）黏弹性层的阻尼特性在低频范围内对结构振动影响较小;（5）可在满足工程实际的情况下适当增加板长,以提高结构的隔振性能.     

Debonding of graded coatings under in-plane compression

Graded materials are multiphase composites with continuously varying thermophysical properties. The concept provides material scientists and engineers with an important tool to develop new materials tailored for some specific applications. One such application of this new class of materials is as top coats or interfacial regions in thermal barrier systems. A widely observed failure mode in these layered materials is known to be interfacial cracking that leads to spallation. In many cases it is the buckling instability of coating under mechanically or thermally induced compressive stresses that triggers spallation. Under in-plane loading since the linear elastic small deformation theory gives only a trivial solution, in this study the plane strain interface crack problem for a graded coating bonded to a homogeneous substrate is formulated by using a kinematically nonlinear continuum theory. Both the instability and the postbuckling problems are considered. The main objective of the study is the investigation of the influence of material nonhomogeneity, kinematic nonlinearity and plate approximation on the critical instability load and on such fracture mechanics parameters as strain energy release rate, stress intensity factors and crack opening displacements.     

A numerical shock profile in a slowly relaxing medium

A numerical solution of the equation which governs the shock solution in a viscoelastic medium displaying stress relaxation is obtained. This equation has been found to be the one which defines the lowest-order approximation to the inner solution of a “generalized Burgers' equation”. The case when the stress relaxation function is algebraically decaying is considered and the results are shown to be consistent with the expected behaviour of the solution.     

The stress field of a sliding inclusion

This paper discusses the stress fields when a spheroidal inclusion, free to slip along its interface, is subjected to a constant nonshear eigenstrain, and when an elastic body containing the inhomogeneity is under all-around tension or uniaxial tension at infinity. In each case the stress field in the inclusion or the inhomogeneity is not constant, contrary to Eshelby's solution. When sliding takes place, the stress increases locally compared with the perfect bonding case, but the elastic energy decreases due to the relaxation. The relative displacement (slip) along the interface is also evaluated.     

Homogenization of nonlinear problems in the mechanics of composites

A further development of the homogenization method is proposed to solve the physically nonlinear equilibrium problems for the laminated plates or the plates made of functionally graded materials. In the linear case, according to this method, the corresponding solution is a superposition of the solution to the global problem in the entire domain and the solution to the local problem in a representative domain, e.g., in a periodicity cell. In the nonlinear case, such a superposition is not valid, which complicates the application of the homogenization method. In order to eliminate this difficulty, it is possible to combine the homogenization method and the linearization method when solving a boundary value problem or a variational problem. In the mechanics of deformable solids, the constitutive relations can be considered as equations with respect to velocities or the stress and strain differentials in time or in the loading parameter. When these equations are linear with respect to velocities, it is possible to use the homogenization method. In this paper such an approach is illustrated by the example of a symmetric laminated plate bent under a uniformly distributed time-dependent load.     

线黏弹性球面发散应力波的频率响应特性

基于线黏弹性球面波Laplace域的理论解, 得到了不同传播距离处粒子速度、粒子位移、应力、应变等力学量的传递函数。以标准线性固体模型为例, 重点讨论了粒子速度频率响应函数的传播特征, 指出随着传播距离的增加, 粒子速度幅频响应函数的高频响应会低于低频响应, 而在理想弹性条件下, 粒子速度幅频响应函数的高频响应一直高于低频响应。以弹性半径为0.025 m的空腔爆炸为例, 采用Laplace数值逆变换方法对粒子速度波形的演化进行了分析, 给出了粒子速度强间断幅值及粒子速度峰值随传播距离变化的衰减规律曲线, 指出黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。     

Dynamic effects of moving load on a poroelastic soil medium by an approximate method

The problem of the dynamic response of a fully saturated poroelastic soil stratum on bedrock subjected to a moving load is studied by using the theory of Mei and Foda under conditions of plane strain. The applied load is considered to be the sum of a large number of harmonics with varying frequency in the form of a Fourier expansion. The method of solution considers the total field to be approximated by the superposition of an elastodynamic problem with modified elastic constants and mass density for the whole domain and a diffusion problem for the pore fluid pressure confined to a boundary layer near the free surface of the medium. Both problems are solved analytically in the frequency domain. The effects of the shear modulus, permeability and porosity of the soil medium and the velocity of the moving load on the dynamic response of the soil layer are numerically evaluated and compared with those obtained by the exact solution of the problem. It is concluded that for fine poroelastic materials, the accuracy of the present method against the exact one is excellent.     

A semi-analytic analysis of shape memory alloy thick-walled cylinders under internal pressure

In this paper, a new method for analysis of the pseudoelastic response of shape memory alloy thick-walled cylinders subjected to internal pressure is proposed. Two cases of short and long cylinders are considered by assuming the plane stress and plane strain conditions. In each case, a three-dimensional phenomenological SMA constitutive model is simplified to obtain the corresponding two-dimensional constitutive relations. The cylinder is partitioned into a finite number of narrow annular regions, and appropriate assumptions are made in order to find a closed-form solution for the equilibrium equations in each annular region. The global solution is obtained by enforcing the stress continuity condition at the interface of the annular regions and imposing the boundary conditions. Several numerical examples are presented to demonstrate the efficiency of the proposed method, and the results are compared with three-dimensional finite element simulations.