考虑两种晶界的各向异性双晶和三晶体晶界附近弹塑性应力场分析

采用率相关的晶体滑移有限元程序对具有不同晶体取向的双晶体晶界附近及三晶体三晶粒交汇处的弹塑性应力场进行了计算,考虑了几何晶界和物理晶界的影响.计算结果表明:双晶体及三晶体考虑几何晶界和物理晶界时,这两种晶界具有相同的应力分布趋势,只是物理晶界比几何晶界的应力集中程度小,双晶体晶界附近有较大的应力梯度,存在应力集中现象.三晶体三晶粒交汇处可能是应力集中之地也可能不造成应力集中,这主要取决于晶粒晶体取向及加载方向.由此可见,要准确理解金属材料的断裂过程,还需要从细观的角度对晶界的力学响应进行细致和深入的研究.     

63Sn-37Pb钎料合金在非比例循环加载下的非弹性流动特性研究

基于63Sn-37Pb钎料舍金在多种非比例应变循环加载下的实验结果,通过考察材料的非弹性应变率与偏应力之间的夹角随累积非弹性应变的变化规律,对63Sn-37Pb钎料合金的非弹性流动特性进行了定量分析。分析结果显示：在相同的非比例加载路径下,当加载等效应变幅值相同时,等效应变率越高,非弹性应变率与偏应力之间夹角平均水平越低,当等效应变率相同时,等效应变幅值越大,相应的夹角平均水平越低;在保持时间范围内,非弹性应变率方向与偏应力方向趋于一致;当非比例路径形状不同时,其非弹性应变率与偏应力之间的夹角随累积非弹性应变的变化趋势明显不同。研究表明,材料的非弹性流动特性强烈依赖于等效应变幅值、等效应变率、保持时间、非比例路径形状。     

RESEARCH ON DYNAMIC MECHANICAL RESPONSE AND CONSTITUTIVE MODEL OF PORCINE LIVER

在交通事故中,腹部器官常因冲击载荷作用而受到伤害,严重时甚至危及生命.肝损伤是腹部损伤中最为常见的一种,致死率很高,了解肝脏的动态力学性能对于事故中肝脏的损伤评估及防护设计有着重要的意义.从新鲜的猪肝组织中取肝实质部分制作试样,利用英斯特朗材料试验机对其进行两种加载率（0.004 s^-1,0.04 s^-1）和两种加载方向（垂直肝脏表面和平行于肝脏表面）的准静态压缩试验,并压缩至破坏.利用改进的分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）实验装置沿平行于肝脏表面方向进行三种高应变率（1 300 s^-1,2 400 s^-1,4 500 s^-1）的动态压缩试验.结果表明：所有应变率下的猪肝压缩应力应变曲线都呈非线性凹向上特征,初始阶段应力值很低,应变约30%后应力幅值显著增大;准静态压缩时,两种应变率（0.004 s^-1,0.04 s^-1）和两种加载方向下肝脏组织破坏应力和破坏应变等力学性能无显著不同,平均破坏应变为48%,平均破坏应力为0.45 MPa.高应变率下肝脏组织的流动应力明显高于准静态下的流动应力,表现出一定的率敏感性.采用Yeoh型超弹性本构模型描述猪肝组织准静态力学性能,基于黏超弹性模型理论,提出了一个能描述肝脏组织从低应变率到高应变率范围力学性能的率相关本构模型,该模型与实验结果有很好的一致性.     

猪肝动态力学性能及本构模型研究

在交通事故中,腹部器官常因冲击载荷作用而受到伤害,严重时甚至危及生命.肝损伤是腹部损伤中最为常见的一种,致死率很高,了解肝脏的动态力学性能对于事故中肝脏的损伤评估及防护设计有着重要的意义.从新鲜的猪肝组织中取肝实质部分制作试样,利用英斯特朗材料试验机对其进行两种加载率(0.004 s~(-1),0.04 s~(-1))和两种加载方向(垂直肝脏表面和平行于肝脏表面)的准静态压缩试验,并压缩至破坏.利用改进的分离式霍普金森压杆(split Hopkinson pressure bar,SHPB)实验装置沿平行于肝脏表面方向进行三种高应变率(1 300 s~(-1),2 400 s~(-1),4 500 s~(-1))的动态压缩试验.结果表明:所有应变率下的猪肝压缩应力应变曲线都呈非线性凹向上特征,初始阶段应力值很低,应变约30%后应力幅值显著增大;准静态压缩时,两种应变率(0.004 s~(-1),0.04 s~(-1))和两种加载方向下肝脏组织破坏应力和破坏应变等力学性能无显著不同,平均破坏应变为48%,平均破坏应力为0.45 MPa.高应变率下肝脏组织的流动应力明显高于准静态下的流动应力,表现出一定的率敏感性.采用Yeoh型超弹性本构模型描述猪肝组织准静态力学性能,基于黏超弹性模型理论,提出了一个能描述肝脏组织从低应变率到高应变率范围力学性能的率相关本构模型,该模型与实验结果有很好的一致性.     

基于初始应变变化的自组装螺旋变形有限元分析

螺旋结构是自然界和工程应用中的基本形态,在不同几何和力学条件下可实现多种形态间的相互转变。本文以双层薄弹性材料结构模型为例(长度L>>宽度W>>厚度H),利用层间存在的不同初始应变,采用三维有限元模拟方法,分析了结构的形态转变情况。结果表明:当加载单层预应变(ε1=0.25)时,改变几何主轴和初始应变主轴方向间的错配角φ,可实现弹性带在圆环、螺旋形、不完整柱形间相互转变;当加载双层预应变(ε1=-ε2=0.25)时,同时改变几何错配角,可实现弹性带在"轮毂"形、纯扭转、一般扭转间的自组装螺旋相互转变。     

不同晶向金属纳米线拉伸力学性能分子动力学研究

在经典等温分子动力学框架下,采用位移加载方式,准静态条件下数值模拟常温条件下金属纳米线的单向拉伸,研究了面心立方晶格（FCC）单晶金属纳米线的弹塑性力学性能。研究发现 <100>,<110>,<111>三个不同晶向纳米线拉伸时呈现不同的拉伸变形力学性能。其中<111>晶向拉伸有最高的屈服强度,<110>晶向屈服屈服最小,特别的是<100>晶向拉伸时屈服应变最大。由于不同的晶向对应纳米线的不同表面,三个晶向的纳米线拉伸呈现不同的应力应变关系曲线,变形过程中的局部结构具有不同的演化方式。分析了纳米单晶铜线的三个晶向拉伸表现不同的等效弹性刚度和屈服强度,讨论了相关的局部位错结构演化过程和与位错发射分解剪应力相关的纳米线塑性变形机理。     

孪生诱发塑性对V-5Cr-5Ti合金应变硬化行为的影响

钒合金(V-Cr-Ti)作为潜在重要的聚变反应堆用结构材料, 近年来受到广泛的关注. 为了研究 V-5Cr-5Ti 合金不同应变率压缩下的应变硬化行为, 特别是孪生对塑性变形的影响, 以位错密度和孪晶演化为基础, 建立了该合金的应变硬化模型. 模型中考虑了孪晶中的位错滑移对材料塑性应变的贡献. 模拟结果表明, 由于孪生诱发塑性, 从而使动态压缩时的位错密度小于准静态加载时的, 这使得 V-5Cr-5Ti 合金在动态压缩时的应变硬化率比准静态加载时的小. 当孪晶形成后, 位错滑移引起的塑性应变率随应变增大而增大, 并逐渐接近加载应变率, 而孪生引起的塑性应变率则随应变增大而减小.      

含与不含晶界空穴的双晶体蠕变行为研究

基于晶体滑移理论,建立了各向异性镍基合金双晶体的蠕变本构模型和蠕变寿命预测模型,通过MARC用户子程序CRPLAW将上述本构模型进行了有限元实现,并对双晶体蠕变行为进行了计算分析,考虑了:(1)晶体取向的影响;(2)垂直、倾斜和平行于外载方向的三种位向晶界情况;(3)晶界处引进空间空穴的影响。结果表明,双晶体上特别是微空穴和晶界附近区域的蠕变应力应变呈现不同的变化规律,对此晶粒晶体取向和晶界位向有较大的影响;微空穴的存在削弱了双晶体的承载能力,显著地影响了双晶体蠕变持久寿命;相同条件下,垂直晶界对双晶体模型的蠕变损伤影响最为强烈,倾斜晶界次之,平行晶界最小;微空穴的生长与晶界位向和晶体取向有强烈的依赖关系,其中垂直晶界更有利于晶体滑移和微空穴生长。     

低速冲击下混凝土少筋梁断裂性能

动态加载下,混凝土中钢筋的阻裂性能一直是冲击动力学研究领域的难点之一。利用落锤试验机对含缺口的混凝土少筋梁进行三点弯曲试验,分析了不同加载速率下梁的冲击力、跨中挠度、混凝土起裂应变率和钢筋应变。实验结果表明:在一定加载速率范围内（0.885～1.252 m/s）,混凝土预制裂缝尖端的裂纹起裂应变率、冲击力最大值、跨中挠度峰值与加载速率呈线性增长关系,当加载速率增至1.771 m/s时,增长趋势减弱;冲击力卸载时,钢筋部分弹性变形恢复导致裂纹产生闭合,裂纹嘴张开位移逐渐减小至恒定值,对裂纹嘴张开位移峰值前的部分曲线进行拟合后得到裂纹嘴张开位移率,结果表明裂纹嘴张开位移率随加载速率的提高而线性增大。     

闪长岩单轴加载过程中声发射和弹性波速度变化规律的研究

为研究闪长岩在单轴加载过程中的声发射和各向波速变化规律,在单轴阶段加载和循环阶段加载条件下,对闪长岩岩样破裂过程中的声发射累计数、不同应力水平不同方向的波速、切线模量、轴向应变速率进行了研究。实验结果表明:(1)随着应力水平的增高,声发射事件数不断增加,在高应力水平(约80%峰值强度)时,声发射累计数急剧增多,随后切线模量出现震荡变化。(2)在加载过程中,压密程度及裂纹扩展方向对波速产生了巨大的影响,导致不同方向波速在不同的应力水平呈现出不同的变化规律,由此可以推测破裂面位置和破裂模式。在较高应力水平下(约60%峰值强度),平行于加载方向的波速趋于稳定,而垂直于加载方向的波速则持续下降,故用垂直于加载方向传播的波速预测岩石的破坏更具可靠性。(3)随着应力的增加,应变速率有逐渐减小的趋势,但临近岩石破裂时无异常变化出现,说明利用变形观测难以预测此类岩石的破坏。以上研究表明,根据纵波波速、声发射累计数和切线模量的变化可以有效预测岩石的破坏。     

Non-slipping JKR model for transversely isotropic materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers.     

花岗岩的加载速率效应及能量机制研究

加载速率对岩石的力学性质以及变形破坏方式具有重要的影响。基于MTS810电液伺服材料试验系统与PCI-2声发射仪对岩样进行不同加载速率作用下的单轴压缩和声发射试验。研究结果表明：（1）在各级加载速率作用下,岩样单轴压缩应力-应变曲线大致经历了压密、弹性、屈服、破坏四个阶段。岩样峰后曲线在加载速率为0.001~0.01 mm/s时出现台阶型分段跌落状,在加载速率为0.01~0.1 mm/s时呈现光滑、陡峭的连续曲线。（2）岩样峰值强度、弹性模量随加载速率的增加而增大,与加载速率对数均呈现三次多项式拟合关系。峰值应变随加载速率的增加而减小,与加载速率对数呈现线性拟合关系。（3）随着加载速率由0.001mm/s增加至0.1mm/s,岩样吸收的总应变能 具有波动性,可释放的弹性应变能 增幅60.42%,耗散应变能 降幅 66.38%, 增幅43.33%, 降幅66.67%,岩样破裂模式由拉剪破坏逐渐向张拉劈裂破坏过渡,岩样破裂块数增多。（4）加载速率为0.001~0.1 mm/s时,岩样破坏方式有所不同,但破坏为同一类损伤过程。单轴压缩状态下,能量耗散使得岩样损伤致使强度丧失,而能量释放使得岩样宏观破裂面贯通,并向着能量释放的方向张裂或弹射破坏。     

Crack initiation direction from interface of bonded dissimilar media

Debonded region of an interface between two dissimilar materials are modeled as a line crack that tends to enhance the initiation of failure by fracture. Depending on the load that interacts with dissimilar materials, no a priori knowledge of how failure would initiate from an existing interface crack is assumed. By application of the strain energy density criterion, potential crack initiation sites are obtained for different biaxial loading states and materials with dissimilar properties.Numerical results are obtained for an epoxy/aluminum medium. In each case, a finite line segment of debonding is assumed. Uniform stresses are applied normal and parallel to the interface so that a biaxial load factor k determines the relative magnitude of biaxiality. Positive and negative k correspond, respectively, to applied tension and compression parallel to the interface. For a fixed ratio of the elastic moduli, crack initiation angles measured from the interface would increase with positive k and decrease with an increase of negative k. These findings are presented for different values of k. The direction of maximum yield initiation could also be determined from the stationary values of the strain density function. These locations are identified with elements that undergo excessive distortion while the possible fracture sites are assumed to coincide with regions where dilatational effects would dominate.     

Study on fracture behavior of ferroelectric ceramics under combined electromechanical loading by using a moiré interferometry technique

This paper discusses an in situ observation of fracture behavior around a crack tip in ferroelectric ceramics under combined electromechanical loading by use of a moiré interferometry technique. The deformation field induced by the electric field and the stress concentration near the crack tip in three-points bending experiments was measured. By analysis of the moiré images it is found that under a constant mechanical load, the electric field almost has no effect on the crack extension in the case that the directions of the poling, electric field and crack extension are perpendicular to each other. When the poling direction is parallel to the crack extension direction and perpendicular to the electric field, the strain decreases faster than that calculated by FEM with and without electrical loading as one goes away from the crack tip. In addition, as the electric field intensity increases, the strain near the crack tip increases, and the strain concentration becomes more significant. The project supported by the National Natural Science Foundation of China (10132010, 10025209, 10232023)     

Buckling of an orthotropic graded coating with an embedded crack bonded to a homogeneous substrate

To simulate buckling of nonuniform coatings, we consider the problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subjected to a compressive loading. The coating is graded in the thickness direction and the material gradient is orthogonal to the crack direction which is parallel with the free surface. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The loading consists of a uniform compressive strain applied away from the crack region. The graded coating is modeled as a nonhomogeneous medium with an orthotropic stress–strain law. Using a nonlinear continuum theory and a suitable perturbation technique, the plane strain problem is reduced to an eigenvalue problem describing the onset of buckling. Using integral transforms, the resulting plane elasticity equations are converted analytically into singular integral equations which are solved numerically to yield the critical buckling strain. The Finite Element Method was additionally used to model the crack problem. The main objective of the paper is to study the influence of material nonhomogeneity on the buckling resistance of the graded layer for various crack positions, coating thicknesses and different orthotropic FGMs.     

Thermal stresses in a two-dimensional infinite medium containing a rigid inclusion embedded in a line crack

This paper examines the problem of finding thermal stresses, caused by a symmetric indentation of a line crack by an inclusion in an infinite isotropic elastic heat conducting solid. The thermal and elastic problems are reduced to a system of triple integral equations. In each case the solution of the triple integral equations is obtained in a closed form. The expressions for the stress intensity factor at the edge of the line crack, the strain energy density function and the resultant pressure applied to the inclusion are obtained. The expression for the displacement component is also obtained. Finally the results for the physical quantities are displayed graphically.     

A unified expression of elastic–plastic notch stress–strain calculation in bodies subjected to multiaxial cyclic loading

In this paper, a simplified thermodynamics analysis of cyclic plastic deformation is performed in order to establish an energy transition relation for describing the elastic–plastic stress and strain behavior of the notch-tip material element in bodies subjected to multiaxial cyclic loads. Based on the thermodynamics analysis, it is deduced that in the case of elastic–plastic deformation, Neuber’s rule inevitably overestimates the actual stress and strain at the notch tip, while the equivalent strain energy density (ESED) method tends to underestimate the actual notch-tip stress and strain. According to the actual energy conversion occurring in the notch-tip material element during cyclic plastic deformation, a unified expression for estimating the elastic–plastic notch stress–strain responses in bodies subjected to multiaxial cyclic loads is developed, of which Neuber’s rule and the ESED method become two particular cases, i.e. upper and lower bound limits of the notch stress and strain estimations. This expression is verified experimentally under both proportional and non-proportional multiaxial cyclic loads and a good agreement between the calculated and the measured notch strains has been achieved. It is also shown that in the case of multiaxial cyclic loading, the unified expression distinctly improves the accuracy of the notch-tip stress–strain estimations in comparison with Neuber’s rule and the ESED method. The unified expression of the notch stress–strain calculation developed in this paper can thus provide a more logical approximate approach for estimating the elastic–plastic notch-tip stress and strain responses of components subjected to lengthy multiaxial cyclic loading histories for local strain approach-based fatigue-crack-initiation life prediction.     

An anisotropic micromechanics model for predicting the rafting direction in Ni-based single crystal superalloys

An anisotropic micromechanics model based on the equivalent inclusion method is developed to investigate the rafting direction of Ni-based single crystal superalloys. The micromechanical model considers actual cubic structure and orthogonal anisotropy properties. The von Mises stress, elastic strain energy density, and hydrostatic pressure in different inclusions of micromechanical model are calculated when applying a tensile or compressive loading along the [001] direction. The calculated results can successfully predict the rafting direction for alloys exhibiting a positive or a negative mismatch, which are in agreement with pervious experimental and theoretical studies. Moreover, the elastic constant differences and mismatch degree of the matrix and precipitate phases and their influences on the rafting direction are carefully discussed.     

Channel cracking in inelastic film/substrate systems

Studies on channel cracking are generally limited to elastic films on elastic or inelastic substrates. There are important applications were the cracking process involves extensive plasticity in both the film and substrate, however. In this work steady-state channel cracking in inelastic thin-film bilayers undergoing large-scale yielding from thermal or mechanical loading is studied with the aid of a plane-strain FEA. The plasticity of the film and substrate, represented by a Ramberg–Osgood constitutive law, each increases the energy release rate (ERR) relative to the linearly-elastic case. This effect is more pronounced under mechanical loading where the entire bilayer undergoes large-scale yielding. To help assess the analytic approach some fragmentation tests are performed using a well-bonding epoxy/aluminum system. The analysis reproduced well the observed dependence of crack initiation strain on film thickness.Ultra-thin films may be well represented by an elastic-perfectly plastic response. For such films on a flexible support the ERR remains fixed as the applied strain exceeds the yield strain of the film. Accordingly, a critical coating thickness exists below which no channel cracking is possible. The explicit relations and graphical data presented may be used for optimal design of such structures against premature failure as well as for determining fracture energy of ductile thin films.     

Splitting of coherent precipitates caused by elastic non-equilibrium

During the process of nickel-based superalloys, gamma prime particles often split into doublets, quartets, or octets. It is an interesting instability phenomenon as it defies conventional surface thermodynamics. Particle splitting is examined through the discrete atom method, which is based on the combination of statistical mechanics and linear elasticity. Splitting phenomena may be classified into two: commensurate elastic instability and incommensurate elastic instability. In the former, the elastic anisotropy of a coherent particle is commensurate with that of the matrix phase, and a non-equilibrium elastic state may cause particle splitting provided that relaxation of the elastic strain energy can exceed the accompanying interfacial energy increase. Incommensurate elastic instability arises when the anisotropy ratios of the precipitate and matrix phase have opposite signs, for example, when the elastically soft direction of the matrix is parallel to the hard direction of the particle phase. The strain energy tends to relax along the elastically soft directions of both phases: thus the particle stretches along its own soft directions while yielding its hard directions to the matrix phase. Such a splitting can end up with a symmetry in the arrangement of particles, but requires a condition that the bulk moduli of both matrix and precipitate phase are comparable to each other.