韧性材料在强动载作用下的损伤演化

采用细观动力学分析方法在高加载率条件下得到了孔洞的动态增长模型。模型考虑了应变率效率、惯性效应和热效应对孔洞增长的影响，给出了孔洞增长的应力临界值表达式。模型的数值分析表明：初始温度Ｔ０对应力临界值Ｐｃｒｉｔ有较大影响；在高应变率条件下惯性和热效应对韧性孔洞动态增长有明显的影响，惯性表现为对孔洞增长的阻碍作用，而热效应应则表现为对孔洞增长的促进作用，数值分析还表明，在低应变率条件下惯性效应可以忽略     

纯铜一维动态损伤研究

在Ｊｏｈｎｓｏｎ方法的基础上，本文采用细观方法提出了一个延性材料动态损伤模型。模型中考虑了材料的硬化，应变率效应和惯性效应。同时进行了纯铜一维层裂实验，对软回收试件进行了显微金相和定量金相分析，观察到了一些有意义的力学现象，利用一维Ｌａｇｒａｎｇｅ有限差分动力学程序。采用本文提出的损伤模型对纯铜的层裂过程进行了数值模拟，得到了比较满意的结果。     

延性介质动态损伤模型

本文借鉴Johnson的基本方法,采用高应变率下的应力应变关系,以一个球形实体中含有球形空穴为基本力学单元,通过分析基元的受力状况,采用奇异摄动方法,推导出了空穴膨胀比的演化方程。方程中考虑了材料的硬化和应变率效应,同时又计及了惯性效应,在考虑空穴成校对损伤的影响时,采用了Chu和Needleman的结果,利用本文提出的模型对两个一维高速碰撞实验进行了数值模拟,计算结果与实验吻合较好。     

纯铜—维动态损伤研究

在Johnson方法的基础上,本文采用细观方法提出了一个延性材料动态损伤模型.模型中考虑了材料的硬化、应变率效应和惯性效应.同时进行了纯铜一维层裂实验,对软回收试件进行了显微金相和定量金相分析,观察到了一些有意义的力学现象.利用一维Lagrange 有限差分动力学程序,采用本文提出的损伤模型对纯铜的层裂过程进行了数值模拟,得到了比较满意的结果.     

混凝土类材料动态压缩强度在多维应力状态下的应变率效应

对混凝土类材料动态压缩应变率效应研究的发展及问题进行了概述,对比不同应力状态下混凝土类材料动态压缩应变率效应的表现特征,揭示了不同加载路径下实测动态强度提高系数的显著差异。研究表明,在高应变率下,基于初始一维应力加载路径的试件将因横向惯性效应导致的侧向围压而演化至多维应力状态,传统霍普金森杆技术无法获得高应变率下基于真实一维应力路径的动态强度提高系数,在强度模型中直接应用实测数据将过高估计材料的动态强度。鉴于应变率效应的加载路径依赖性,将仅包含应变率的强度提高系数模型扩展至同时计及应变率和应力状态的多维应力状态模型,并结合Drucker-Prager准则在强度模型中给予了实现。针对具有自由和约束边界试件开展的数值霍普金森杆实验表明,多维应力状态下的应变率效应模型可以考虑应变率效应随应力状态改变的特点,从而准确预测该类材料的动态压缩强度。研究结果可为正确应用霍普金森杆技术确定脆性材料的动态压缩强度提供参考。

     

含氦泡辐照老化材料层裂损伤计算方法分析

辐照条件下,一些材料内部产生大量的氦泡等微缺陷,氦泡的大小和数密度随着辐照年限的增长而增长。氦泡分布特征的变化不仅影响材料本身的物理、力学性质,而且直接影响材料层裂损伤演化后期材料破坏颗粒度的分布特征。延性材料的层裂损伤演化过程一般包括孔洞的成核、增长和汇合,但因已有孔洞对新成核孔洞存在抑制作用,当初始孔洞数密度达到一定临界值时,材料内部没有新的孔洞成核,因此,层裂损伤的计算可以不考虑新孔洞成核的影响。本文中基于损伤早期演化的特征,给出了这一临界值的计算方法,并进一步探讨了含氦泡辐照老化钚材料层裂损伤的计算方法。同时,在完善孔洞增长(void growth, VG)层裂损伤模型中参数的确定方法的基础上,借助含氦泡常规铝材料的层裂实验结果,对此问题进行了定性的分析：在氦泡尺寸变化不大的情况下,当氦泡浓度低于临界氦泡浓度时,需要考虑初始氦泡以及新增孔洞的综合影响;反之,可以采用简单的层裂损伤模型,不需要计算孔洞成核,但由于增长孔洞之间的相互影响,损伤模型的初始损伤参数需要重新确定。     

基体微尺度效应对弹塑性多孔洞材料本构势及孔洞长大的影响

基于尺度相关的应变梯度塑性SG理论,对含孔洞的球形体胞模型进行了解析分析,得到了在基体梯度塑性环境下球形孔洞的演化规律,给出了弹塑性孔洞材料的宏观屈服面方程.与现有的基于尺度无关塑性理论的Gurson模型相比,该模型考虑了基体材料的特征长度l与孔洞半径a之比λ(λ=l/a)对多孔材料宏观屈服面和孔洞演化规律的影响.当不计基体材料的塑性梯度效应和硬化效应时,该模型能退化到经典的Gurson模型.     

延性材料动态损伤演化模型研究

从微损伤系统的统计描述出发,应用Ｍｏｒｉ－Ｔａｎａｋａ理论的物理思想,考虑了微孔洞间的相互作用,建立了延性材料动态损伤的统计演化模型。结果能够描述材料在外载作用下的损伤弱化效应。还分析了材料变形中的两种不同的耗散机制,考虑了含损伤的本构关系。将上述模型用于几种材料的层裂过程数值模拟,并与实验结果作了初步比较。     

基于对数应力率大应变本构关系的参数标定研究

在所有率型弹塑性本构模型中,只有对数应力率对应的本构模型能够满足自适应准则.基于对数应力率,采用实心圆轴扭转实验,对大应变弹塑性本构模型中的参数标定问题进行了讨论.推导出了考虑Swift效应时端部自由实心圆轴扭转变形的变形率、对数旋率、Kirchhoff应力及Kirchhoff应力的对数应力率.对于等向强化大应变弹塑性本构关系,给出了由实心圆轴扭转实验标定的、基于Kirchhhoff应力对数应力率的本构关系中塑性刚度函数的表达式.分析了扭转圆轴的Swift效应对塑性刚度函数的影响.结果表明,实心圆轴扭转的轴向伸长变形和径向变形对基于对数应力率大应变本构关系中的塑性刚度函数都有影响.当不考虑Swift效应时,所得塑性刚度函数表达式与不考虑Swift效应时基于Jaumann应力率的塑性刚度函数表达式相同.     

十次对称二维准晶中的热应力分析

推导了当考虑热效应时十次对称二维准晶体平面应变问题的通解表示．作为应用，采用所获得的通解首先得到了十次对称二维准晶体中的一个点热源所引起的声子场和相位子场，给出了点热源所引起的声子场和相位子场应力分量的解析表达式；接着获得了在均匀热流作用下十次对称二维准晶体中-绝缘椭圆孔洞所引起的热应力问题的弹性解答，给出了沿椭圆边界环向应力分布的解析表达式；当椭圆的短轴趋于零时，则获得了裂纹问题的解答，给出了应力强度因子、裂纹表面张开位移及能量释放率的解析表达式；推导了在任意热载荷作用下裂尖附近的渐近场．     

EVOLUTION OF VOIDS IN DUCTILE POROUS MATERIALS AT HIGH STRAIN RATE

In the present work,a dynamic damage model in ductile materials underthe application of dynamic general stresses loading is presented.The evolution equationof ductile voids has the closed form,in which work-hardening,the change of surfaceenergy of voids,rate-dependent,inertial effects are taken into account.Theexpressions of critical stresses for the growth and compaction of voids are directlyobtained from the evolution equations of voids.From the expressions,the resultobtained by Carroll and Holt,as a special example,is given.Numerical analysis ofthe model indicates that the growth of voids is sensitive to the strain rates.The voidsgrow quickly as the increase of strain rates.It is also shown that the influence of theinertial effects on the void growth is great at high loading rates.It appears to resist thegrowth of voids.In addition,a dynamic collapse model of ductile voids is alsoproposed,which can be applied to study the problems of compaction in powder andother materials.     

延性金属层裂自由面速度曲线特征多尺度模拟研究

以延性金属钽为研究对象,对钽在平板撞击下的层裂行为进行了多尺度下的数值模拟研究,从微观视角对自由面速度曲线上的典型特征进行了新的解读。在宏观尺度,对比分析了光滑粒子流体动力学法（smootfied particle hydrodynamics, SPH）与Lagrange网格法以及几种本构模型的模拟结果及其适用性。通过与实验数据的对比表明,Steinberg-Cochran-Guinan本构模型在层裂模拟中与实验数据吻合较好,通过改变加载条件获得了不同应变率下的自由面速度曲线,分析了不同应变率下的自由面速度曲线中的典型特征。在微观尺度,采用分子动力学方法获得层裂区域内损伤演化情况,揭示了宏观尺度自由面速度曲线典型特征所蕴含的物理内涵。分析表明,层裂表现为材料内部微孔洞形核、长大和聚集的损伤演化过程,自由面速度曲线上的典型特征与层裂区域的损伤演化过程存在密切关联。Pullback信号是层裂区域内微孔洞形核的宏观表征;自由面速度曲线的下降幅值在一定程度上反映了微孔洞的形核条件,由此计算得到的层裂强度实际上是微孔洞的形核强度。此外,Pullback信号后的速度回跳速率反映了微损伤演化的速率。     

Combined grain size, strain rate and loading condition effects on mechanical behavior of nanocrystalline Cu under high strain rates

Molecular dynamics simulations of nanocrys-talline Cu with average grain sizes of 3.1 nm, 6.2 nm, 12.4 nm and 18.6 nm under uniaxial strain and stress tension at strain rates of 10 8 s 1 , 10 9 s 1 and 10 10 s 1 are performed to study the combined grain size, strain rate and loading condition effects on mechanical properties. It is found that the strength of nanocrystalline Cu increases as grain size increases regardless of loading condition. Both the strength and ductility of nanocrystalline Cu increase with strain rate except that there is no monotonic relation between the strength and strain rate for specimens under uniaxial strain loading. Moreover, the strength and ductility of specimens under uniaxial strain loading are lower than those under uniaxial stress loading. The nucleation of voids at grain boundaries and their subsequent growth characterize the failure of specimens under uniaxial strain loading, while grain boundary sliding and necking dominate the failure of specimens under uniaxial stress loading. The rate dependent strength is mainly caused by the dynamic wave effect that limits dislocation motion, while combined twinning and slipping mechanism makes the material more ductile at higher strain rates.     

Ductile fracture by cavity nucleation between larger voids

A mechanism of ductile fracture involving the interaction of relatively large voids with small-scale voids is studied by a computational model. The larger voids are described as circular cylindrical holes arranged in a doubly periodic array in the initial state. In the matrix material between these voids the nucleation and growth of much smaller voids is accounted for by using approximate constitutive equations for a ductile, porous medium. The computations show bands of highly localized straining and void growth, initiating at the surfaces of larger voids and growing into the matrix material, until the bands connect two neighbouring voids. The materials are analysed both under plane strain conditions and under conditions approximating those in a round tensile bar. The failure strains obtained under different principal stress ratios show rather good agreement when plotted against a measure of the stress-triaxiality.     

损伤演化对韧性金属碎裂过程的影响

固体在冲击拉伸载荷作用下会断裂成多个碎片,基于线性内聚力断裂假设的Mott-Grady模型能较好地预测碎裂过程所产生的平均碎片尺度的下限。然而实际上,韧性金属的损伤演化是多元化的,为此通过数值模拟方法研究了不同损伤演化规律对韧性碎裂过程的影响。利用ABAQUS/Explicit动态有限元软件数值再现了韧性金属杆（45钢）在高应变率下拉伸碎裂的过程,分析了线性和非线性损伤演化对韧性碎裂过程的影响规律。结果表明：损伤演化规律对韧性金属的碎裂过程具有显著影响,非线性指标α越大,碎裂过程产生的碎片数越少;Grady-Kipp碎裂公式仍能在一定范围内预测韧性碎裂过程中产生的碎片尺寸;当非线性指标α远大于零时,在较低冲击拉伸载荷作用下,数值模拟结果和Grady-Kipp模型预测值偏差较大,随着应变率增大,数值模拟结果与Grady-Kipp模型预测值吻合较好。

     

The dynamic growth of a single void in a viscoplastic material under transient hydrostatic loading

We have examined the problem of the dynamic growth of a single spherical void in an elastic-viscoplastic medium, with a view towards addressing a number of problems that arise during the dynamic failure of metals. Particular attention is paid to inertial, thermal and rate-dependent effects, which have not previously been thoroughly studied in a combined setting. It is shown that the critical stress for unstable growth of the void in the quasistatic case is strongly affected by the thermal softening of the material (in adiabatic calculations). Thermal softening has the effect of lowering the critical stress, and has a stronger influence at high strain hardening exponents. It is shown that the thermally diffusive case for quasistatic void growth in rate-dependent materials is strongly affected by the initial void size, because of the length scale introduced by the thermal diffusion. The effects of inertia are quantified, and it is demonstrated that inertial effects are small in the early stages of void growth and are strongly dependent on the initial size of the void and the rate of loading. Under supercritical loading for the inertial problem, voids of all sizes achieve a constant absolute void growth rate in the long term. Inertia first impedes, but finally promotes dynamic void growth under a subcritical loading. For dynamic void growth, the effect of rate-hardening is to reduce the rate of void growth in comparison to the rate-independent case, and to reduce the final relative void growth achieved.     

聚氯乙烯弹性体静动态力学性能及本构模型

为揭示聚氯乙烯弹性体在静、动态载荷下的力学性能,采用万能材料试验机和改进的分离式霍普金森压杆实验装置获得了材料在应变率为0.001、0.01、0.1、1 510、2 260和3 000 s⁻¹下的应力应变曲线,并以屈服强度为整形器优选参数,对比了紫铜、铜版纸和铅等3种整形器材料的整形效果。使用修正的ZWT非线性黏弹性本构模型描述聚氯乙烯弹性体在静、动态载荷下的力学性能。结果表明：聚氯乙烯弹性体在静态载荷下具有应变率效应和显著的超弹性特性,动态载荷下表现出较明显的应变率效应和较强的抗变形能力,且静动态载荷下的力学行为受应变历史影响较大。3种整形器材料中铜版纸的整形效果最好。修正后的ZWT非线性黏弹性本构模型能够得到统一参数的本构表达式,且各应变率下的拟合结果与实验结果具有较好的一致性。

     

Modelling of compressive behaviour of concrete-like materials at high strain rate

Uniaxial compression tests are the most common tests for characterizing the strength of concrete-like materials. The dynamic compression strength of concrete-like material is typically obtained by Split Hopkinson Pressure Bar (SHPB) tests. The increase in material strength under dynamic loading is usually attributed to the strain rate effect and modelled with a dynamic increase factor (DIF). However, it was observed by some researchers that the radial inertial confinement caused apparent increase of dynamic strength of concrete-like specimen in SHPB tests. They attributed the material strength increase to this inertial effect, instead of the strain rate effect. In the present study, numerical analyses are performed to investigate the compressive behaviour of concrete-like material at high strain rates. A homogeneous macroscale model and a heterogeneous mesoscale model are developed in the study. In the macroscale model, the material is assumed to be homogeneous and isotropic. In the mesoscale model, the test sample is modelled as a three-phase composite consisting of aggregate, mortar matrix and interfacial transaction zone (ITZ) between the aggregate and the mortar matrix. The aggregate is assumed to be circular and the ITZ is modelled as a thin boundary around the aggregate. In the both models, the materials are assumed to be insensitive to the strain rate first. Therefore, the obtained strength enhancement is only due to the inertial confinement. Strain rate sensitive material properties are then used in the two models in the calculations. Numerical simulations of the concrete samples under compression at different strain rates are carried out. The relative contribution of the inertial effect and the strain rate effect on the compressive strength DIF is examined based on the numerical results. The failure process of concrete specimen is also studied.     

Effect of stress-state and spacing on voids in a shear-field

Unit cell model analyses are carried out for a material with a periodic array of voids, subject to shear loading. Thus the focus is on ductile fracture in conditions of low stress triaxiality. It has been shown recently that voids in shear are flattened out to micro-cracks, which rotate and elongate until interaction with neighboring micro-cracks gives coalescence, so that the failure mechanism is very different from that under tensile loading. In the present studies the plane strain unit cell has fully periodic boundary conditions, so that any combination of the stress components in the overall average stress state can be prescribed. This also allows for studies of the effect of different initial void spacing in the two in-plane coordinate directions. The stress states considered are essentially simple shear, with various levels of tensile stresses or compressive stresses superposed, i.e. low positive stress triaxiality or even negative stress triaxiality. For high aspect ratio unit cells a clear localization band is found inside the cell, which actually represents several parallel bands, due to periodicity. In the materials represented by a low aspect ratio unit cell localization would also occur after that the maximum shear stress has been passed, but this is not shown when periodicity is enforced. The effect of superposed tensile or compressive stresses is found to be bigger for high aspect ratio unit cells than for low aspect ratios.