结晶特性与制造工艺对炸药件力学性能的影响

通过对塑料粘结炸药(PBX)在压制、热老化及贮存中炸药HMX、TATB和粘结剂F的性能变化规律的研究,揭示了加工和贮存条件对炸药件的力学性能的影响。得出如下结论:1)随着老化温度的提高,粘结剂F结晶度增加。TATB基PBX炸药经老化后力学性能没有明显变化,说明粘结剂结晶度和炸药颗粒度的变化对炸药总体性能影响不大。2)钢模压制的TATB基PBX药柱在经历多次温度循环后,TATB与粘结剂F界面的作用有所减弱,药柱内部产生由脱粘引起的缺陷,其力学强度下降。TATB基PBX药柱的力学强度与模量均随着环境温度的升高而呈下降趋势,而等静压成型能明显改善TATB基PBX的力学性能。     

甲烷及其氟氯代物在TATB晶体表面吸附的模

采用Discover/MaterialStudio中的分子动力学方法(COMPASS力场、NVT)及结构优化方法,模拟了甲烷及其氟氯代物(共15种化合物)在TATB超胞的(001)、(100)及(010)晶面的吸附作用。结果表明:(1)TATB晶面与大多数分子相互作用的强弱次序为:(010)(100)(001),这可能与晶面结构差异及氢键的形成有关;(2)三种晶面上,吸附能最大(为负值)的都是CH4、CF4,表明如果氟聚合物中氢或氟的含量过高均可能导致聚合物与TATB相互作用减弱。     

复合材料界面与其力学性能的关系

复合材料是一种具有多相结构的材料。它由基体(高聚物树脂、金属、陶瓷等)和增强剂(纤维或粉粒状填料)组合而成。这种材料的特点,在于它能根据使用的需要采用不同的宏观或微观的组合形式。换言之,可以进行材料设计。同时,复合材料不仅能综合各组成原材料的特色,而且还有可能具备原材料所没有的特性。因此,复合材料在材料科学领域中很有发展前途,目前已经获得广泛的应用。由于复合材料是一种多相结构的材料,因此必然包含着各组分之间构成的界面。这个界面对手复合材料的性质,尤其是力学性能,起了极其重要的作用。      

一种预测颗粒增强复合材料界面力学性能的新方法

界面在颗粒增强复合材料中起到传递载荷的关键作用, 界面性能对复合材料整体力学行为产生重要影响. 然而由于复合材料内部结构较为复杂, 颗粒与基体间的界面强度和界面断裂韧性难以确定, 尤其是法向与切向界面强度的分别预测缺乏有效方法. 本文以氧化锆颗粒增强聚二甲基硅氧烷(PDMS)复合材料为研究对象, 提出一种预测颗粒增强复合材料界面力学性能的新方法. 首先, 实验获得纯PDMS基体材料及单颗粒填充PDMS试样的单轴拉伸应力$\!-\!$应变曲线, 标定出PDMS基体材料的单轴拉伸超弹性本构关系; 其次, 建立与单颗粒填充试样一致的有限元模型, 选择特定的黏结区模型描述界面力学行为, 通过样品不同阶段拉伸力学响应的实验与数值结果对比, 分别给出颗粒与基体界面的法向强度、切向强度及界面断裂韧性; 进一步应用标定的界面力学参数, 开展不同尺寸及不同数目颗粒填充试样的实验与数值结果比较, 验证界面性能预测结果的合理性. 本文提出的界面力学性能预测方法简便、易操作、精度高, 对定量预测颗粒增强复合材料的力学性能具有一定帮助, 亦对定量预测纤维增强复合材料的界面性能具有一定参考意义.      

多种碳/环氧复合材料断口形貌及其断裂模型分析

对金属材料断口的形貌分析已有了较成熟的研究。复合材料是各向异性材料,对这种材料的微结构(纤维、基体及界面)的显微形貌,特别是从力学角度进行分析还缺少深入的研究。而且,复合材料本身是一个力学结构,纤维主要起承载作用,基体起连接纤维和传递载荷的作用。因此,研究复合材料的力学性能和破坏机理,仅 ...     

烤燃条件下HMX/TATB基混合炸药多步热分解反应计算

采用多步热分解反应动力学模型,描述单质炸药热分解反应,提出了多组分网格单元计算方法,对以HMX/TATB为基的多元混合炸药在烤燃条件下的热反应过程进行了计算。通过炸药烤燃实验测量了炸药内部温度,获得了炸药点火时间,验证了计算的准确性。分析了混合炸药组成比例的变化对炸药热反应性能的影响。在HMX/TATB混合炸药热反应阶段,主要是HMX发生分解反应释放热量,TATB的反应量很少,随着混合炸药中TATB含量的增多,炸药的点火时间逐渐增长,点火温度逐渐增高,炸药热安全性增强。     

共混/填充高分子复合材料界面力学行为实验研究进展

本文首先介绍了高聚物复合材料界面微观结构的特点和界面破坏的微观形式, 在此基础上,详细回顾了粒子填充高分子复合材料和共混高分子合金两类材料的界面力学性能实验研究的进展,总结了研究粒子/基体界面力学性能的实验方法和用这些方法获得的实验结果。介绍的实验方法从宏观拉伸实验到细观在位拉伸实验, 从定性的断口形貌分析到定量的界面强度和粘结能计算, 从不同层次和角度阐述了研究界面性能各种方法的优缺点。最后, 作者根据综合分析的结果, 指出了实验研究界面力学性能遇到的主要困难来自实验仪器的局限性和实验方法的匮乏性。为了获得可靠的界面性能, 一方面应该设计出利于表征的界面, 另一方面需要在实验技术上有所创新。相信该文对于深入了解和研究共混/填充高分子复合材料的界面问题, 以及其它类型材料的界面问题具有一定的借鉴意义。      

基于构型力内变量的界面损伤问题研究

本文基于材料构型力的基本理论和损伤力学中含有内变量的热力学框架,提出了新的损伤变量定义方式,为研究界面损伤问题提供了一种新思路。首先,基于双相弹性体的能量分析,给出界面材料构型力表达式,通过构型力的离散化方法,实现了其在有限元中的数值计算。其次,定义构型力为界面损伤内变量,进而提出一种新的损伤演化模型,并采用刚度劣化的方法,对该界面损伤模型进行数值实现。最后,通过对复合材料界面损伤问题（有裂纹或无裂纹）进行数值模拟,分析了其界面损伤发展趋势,探讨了此模型的合理性和优越性。基于构型力内变量的界面损伤模型,可为研究复合材料的界面损伤失效问题提供一种普适性的方法。     

界面强度对玻璃微珠填充聚丙烯力学性能的影响

本文针对玻璃微珠填充聚丙烯这一刚性粒子填充聚合物复合体系进行了实验研究。通过偶联剂改性对比,研究了该聚合物复合材料在不同界面粘结状态下的宏观拉伸、冲击力学性能。此外,根据冲击破坏断面的电镜观测结果,发现复合体系的断裂和增韧机制随界面粘结强度不同而发生改变,界面改性使得材料抗冲击破坏能力得到增强。本文还采用在位拉伸过程中的细观观测方法,观测到材料在一维应力作用下,刚性粒子和基体界面的脱粘、开裂过程,分析了该复合体系细观结构和宏观力学性能之间的关系,发现界面改性对于材料细观结构的界面脱粘和宏观屈服现象的重要影响,为发展新型复合材料提供了实验依据。     

用直径圆盘试验评价小样品塑料粘结炸药拉伸性能的初步研究

在圆盘试验的力学模型基础上 ,测定了JB90 14等七种TATB基PBX的间接拉伸性能 ,为进一步完善该方法做了初步研究。结果显示 :所有PBX均产生中心开裂现象 ,与该模型预测一致 ;同时该测试方法能有效地反映出JB90 14及其经不同改性后得到的三种PBX的整体力学性能变化趋势。这说明圆盘试验方法适用于PBX。通过比较这七种PBX拉伸性能 ,认为JB90 14及其改性的JB90 14 G1应值得注意。     

The effects of shear on delamination in layered materials

The effect of transverse shear on delamination in layered, isotropic, linear-elastic materials has been determined. In contrast to the effects of an axial load or a bending moment on the energy-release rate for delamination, the effects of shear depend on the details of the deformation in the crack-tip region. It therefore does not appear to be possible to deduce rigorous expressions for the shear component of the energy-release rate based on steady-state energy arguments or on any type of modified beam theory. The expressions for the shear component of the energy-release rate presented in this work have been obtained using finite-element approaches. By combining these results with earlier expressions for the bending-moment and axial-force components of the energy-release rates, the framework for analyzing delamination in this type of geometry has been extended to the completely general case of any arbitrary loading. The relationship between the effects of shear and other fracture phenomena such as crack-tip rotations, elastic foundations and cohesive zones are discussed in the final sections of this paper.     

Transmission of a screw dislocation across a coherent, non-slipping interface

Current research on nanocrystalline metals and nanoscale multilayer thin films suggests extraordinary plastic strength is due to confinement of slip to individual grains or layers. To assess the magnitude of confinement, a Peierls model of slip transmission of a screw dislocation across a coherent, non-slipping interface is presented. The results reflect that large interfacial barriers to transmission are generated by rapid fluctuations in dislocation line energy near the interface due to elastic modulus mismatch, stacking fault energy mismatch, and antiphase boundary energy for transmission into an ordered phase. Coherency stress is predicted to dramatically alter the dislocation core configuration and impart additional strength regardless of the sign. Contributions to strength are not additive due to nonlinear coupling via the dislocation core configuration. The predicted barrier strength for a coherent (0 0 1) Cu/Ni interface is comparable to atomistic (EAM) results but larger than estimates from hardness data.     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

Micromechanical analysis of kinematic hardening in natural clay

This paper presents a micromechanical analysis of the macroscopic behaviour of natural clay. A microstructural stress–strain model for clayey material has been developed which considers clay as a collection of clusters. The deformation of a representative volume of the material is generated by mobilizing and compressing all the clusters along their contact planes. Numerical simulations of multistage drained triaxial stress paths on Otaniemi clay have been performed and compared the numerical results to the experimental ones in order to validate the modelling approach. Then, the numerical results obtained at the microscopic level were analysed in order to explain the induced anisotropy observed in the clay behaviour at the macroscopic level. The evolution of the state variables at each contact plane during loading can explain the changes in shape and position in the stress space of the yield surface at the macroscopic level, as well as the rotation of the axes of anisotropy of the material.     

Finite element modeling of elasto-plastic contact between rough surfaces

This paper presents a finite element calculation of frictionless, non-adhesive, contact between a rigid plane and an elasto-plastic solid with a self-affine fractal surface. The calculations are conducted within an explicit dynamic Lagrangian framework. The elasto-plastic response of the material is described by a J₂ isotropic plasticity law. Parametric studies are used to establish general relations between contact properties and key material parameters. In all cases, the contact area A rises linearly with the applied load. The rate of increase grows as the yield stress σ_y decreases, scaling as a power of σ_y over the range typical of real materials. Results for A from different plasticity laws and surface morphologies can all be described by a simple scaling formula. Plasticity produces qualitative changes in the distributions of local pressures in the contact and of the size of connected contact regions. The probability of large local pressures is decreased, while large clusters become more likely. Loading-unloading cycles are considered and the total plastic work is found to be nearly constant over a wide range of yield stresses.     

Iterative dipole moment method for calculating dielectrophoretic forces of particle-particle electric field interactions

Dielectrophoresis (DEP) is one of the most popular techniques for bio-particle manipulation in microfluidic systems. Traditional calculation of dielectrophoretic forces of single particle based on the approximation of equivalent dipole moment (EDM) cannot be directly applied on the dense particle interactions in an electrical field. The Maxwell stress tensor (MST) method is strictly accurate in the theory for dielectrophoretic forces of particle interaction, but the cumbersome and complicated numerical computation greatly limits its practical applications. A novel iterative dipole moment (IDM) method is presented in this work for calculating the dielectrophoretic forces of particle-particle interactions. The accuracy, convergence, and simplicity of the IDM are confirmed by a series of examples of two-particle interaction in a DC/AC electrical field. The results indicate that the IDM is able to calculate the DEP particle interaction forces in good agreement with the MST method. The IDM is a purely analytical operation and does not require complicated numerical computation for solving the differential equations of an electrical field while the particle is moving.     

Influence of anisotropic elasticity on the mechanical properties of fivefold twinned nanowires

Previous atomistic simulations and experiments have shown an increased Young's modulus and yield strength of fivefold twinned (FT) face-centered cubic metal nanowires (NWs) when compared to single crystalline (SC) NWs of the same orientation. Here we report the results of atomistic simulations of SC and FT Ag, Al, Au, Cu and Ni NWs with diameters between 2 and 50 nm under tension and compression. The simulations show that the differences in Young's modulus between SC and FT NWs are correlated with the elastic anisotropy of the metal, with Al showing a decreased Young's modulus. We develop a simple analytical model based on disclination theory and constraint anisotropic elasticity to explain the trend in the difference of Young's modulus between SC and FT NWs. Taking into account the role of surface stresses and the elastic properties of twin boundaries allows to account for the observed size effect in Young's modulus. The model furthermore explains the different relative yield strengths in tension and compression as well as the material and loading dependent failure mechanisms in FTNWs.     

Mean field modelling of the plastic behaviour of co-continuous dual-phase alloys with strong morphological anisotropy

This work addresses the plastic flow properties of a composite material in which the reinforcing phase is continuous and cannot be suitably represented by isolated ellipsoidal inclusions. The dual-phase metal under consideration is composed of a network of Inconel-601 fibres infiltrated by pure aluminium. Hence, both phases exhibit elastic–plastic behaviour and are continuous in the three dimensions of space. The fibre network presents a large morphological anisotropy that is reflected in the mechanical response of the composite. The modelling is based on Eshelby’s equivalent inclusion theory. Strain partitioning between the phases is computed incrementally based on tangent operators derived from the isotropic response of individual phases. Assessment of the model relies on extensive experimental data. Uniaxial tensile tests, involving measurement of the Lankford coefficient, have been performed at various temperatures on samples containing different volume fractions of fibres. Measurement of the phase stresses by neutron diffraction supplements the information provided by the macroscopic stress–strain curves. It is demonstrated that predictions are valid only when the micro–macro averaging scheme accounts for the co-continuous character of the constitutive phases.