PBX材料宏细观断裂行为的数字散斑相关法实验研究

 采用半圆盘弯曲实验和数字散斑相关方法,对高聚物粘结炸药（PBX）的宏、细观断裂行为进行了实验研究。宏观上,带有预制裂纹的半圆盘试样发生拉伸破坏,利用数字散斑相关技术得到了试样的应变场和位移矢量场分布,定量分析了试样全场的变形特征,并测得了PBX材料的平面应变断裂韧性;细观上,用配有加载装置的扫描电子显微镜对含预制裂纹的半圆盘试样间接拉伸下的损伤演化和破坏过程进行了实时原位观察,借助于数字散斑相关方法,定量分析了试样损伤局部化特征。结果表明,将数字散斑相关方法用于研究PBX材料宏、细观尺度上的变形破坏问题是有效的。     

PBX炸药含各向异性损伤的黏弹性统计微裂纹本构模型初步研究

建立含损伤本构模型是研究炸药动态力学响应规律的基础。基于PBX炸药材料的宏观黏弹性特征和细观上微裂纹面的方向性,建立了含各向异性损伤的黏弹性统计微裂纹(Aniso-Visco SCRAM)本构模型,简化后得到单轴应力加载下的本构方程。利用数值计算程序,以PBX9501为例,分析了微裂纹扩展的各向异性、PBX炸药破坏强度及临界应变的拉压异性和应变率相关性,考察了微裂纹数密度、初始微裂纹尺寸、微裂纹面摩擦系数及断裂表面能4个主要参数的敏感性及影响规律。结果表明,它们对微裂纹的扩展演化有较大影响,进而导致材料表现出不同的力学响应。     

PBX炸药粘弹性统计微裂纹本构模型的改进

基于描述PBX炸药本构关系的粘弹性统计微裂纹模型(Visco-SCRAM),通过修正其微裂纹方向矢量的取值范围、初始缺陷尺寸和体量本构关系,得到模型的改进形式。将改进的模型嵌入有限元程序中,开展PBX炸药JO-9159的平板撞击实验数值模拟研究,讨论了微裂纹数密度参数对计算结果的影响。结合实测曲线,对模型改进前、后的计算结果进行对比和分析。结果表明,改进的模型能够更合理地描述低速冲击下PBX炸药的损伤和破坏。     

某PBX炸药的动态力学性能研究

 PBX炸药作为现代武器的主装药,它的力学行为决定着武器的生存能力。为了研究PBX炸药的动态力学特性,采用分离式霍普金森压杆（Split Hopkinson Pressure Bar,SHPB）作为加载手段,结合半导体应变片测试技术和压电晶体监测技术,保证了实验数据的有效性。利用SHPB加载波整形技术,实现了材料两端应力平衡和常应变率加载,得到了不同应变率（90~410 s^-1）下材料的应力-应变曲线。根据材料的模量、破坏强度和破坏应变随应变率的变化规律,采用粘性修正的Sargin模型,得到了该PBX炸药在单轴压缩下的唯象本构模型,模拟结果与实验曲线符合较好。     

压缩加载下三种含能材料细观破坏特征观察

 通过扫描电镜,对三种含能材料在不同应变率压缩条件下回收的试样断面进行了细观观察,得到了三种含能材料的细观破坏特征,并对其破坏机理进行了分析。同时结合高聚物粘结炸药（PBX炸药）和压装B炸药组分的比较,分析了两种炸药在高应变率条件下的不同破坏模式。     

JOB-9003炸药试件机械加载损伤超声检测

近来PBX材料损伤、损伤的非破坏过程检测已成为PBX材料研究和库存老化研究关注的热点和难点。为解决高聚物粘结炸药（PBXs）在加工或贮存过程中的损伤开裂的技术难题，开展了在典型力学加载条件下对炸药试样的超声波检测研究。依据材料的细观损伤的基本理论，对高于200kHz频率的线弹性各向同性体中损伤为一标量，泊松比不随损伤而变，且损伤因子D==1-p1C^2 L1/p0C^2 L0，p0，C L0为初始密度、声速，p1，CL1为有损伤时的密度、声速。     

三氨基三硝基苯基高聚物粘结炸药热力学性质的理论计算研究

高聚物粘结炸药(PBX)的热力学性质是用于炸药结构响应、安全性评估、数值模拟分析等的重要参数.由于PBX结构的多尺度特性,完全采取实验方法精细表征这些参数存在巨大的挑战.本文运用第一性原理和分子动力学计算的方法,系统研究了三氨基三硝基苯(TATB)基高聚物粘结炸药的热力学参数和界面热传导性质.利用散射失配模型研究了TATB与聚偏二氟乙烯(PVDF)界面的热传导过程,发现热导率随温度升高而上升,并且在高温情况下接近于定值.基于分子动力学获得的TATB热导率并结合界面热导率,分析了PBX炸药的热导与颗粒尺寸的关系,当颗粒尺寸大于100 nm时,界面热阻对于PBX热导率的影响有限.     

弹粘塑性双球壳塌缩热点反应模型

 基于Kim的弹粘塑性单球壳塌缩模型,考虑PBX炸药中的粘结剂效应,假设炸药和粘结剂均为弹粘塑性材料,建立了弹粘塑性双球壳塌缩热点反应模型,给出了炸药球壳在冲击压力作用下的速度、应变、温度和化学反应速率的时空分布,以及新的热点反应速率理论表达式。把新的热点反应项与Kim的低压下慢反应项和张震宇提出的高压反应速率方程相结合,得到了新的冲击起爆三项式细观反应速率模型。把该模型加入DYNA2D中,模拟了PBX-9501炸药的一维冲击起爆过程,结果表明：该模型除了可以解释炸药颗粒度和孔隙度的影响外,还可以较好地描述粘结剂强度和含量对PBX炸药冲击起爆感度的影响。     

一种高聚物粘结炸药和B炸药的本构关系研究

 利用万能材料试验机对一种高聚物粘结炸药（简称PBX炸药）和压装B炸药进行了准静态压缩实验,得到了两种炸药在不同应变率下的应力应变曲线。通过两种炸药在不同应变率下力学性能的研究,建立了两种炸药含应变率效应的本构方程。结果表明,建立的本构方程能较好地描述两种炸药在弹性阶段和强化阶段的力学性能。     

炸药落锤实验及样品厚度效应的三维数值模拟

为了研究落锤实验中炸药的概率点火行为以及试样厚度对点火的影响,对脆性炸药PBX-2的落锤实验开展了三维数值模拟。采用有限元与离散元相结合的方法,并考虑了炸药材料的非均匀性,获得了炸药在不同落高下的点火概率分布,落高计算结果与文献报道的实验结果相符。此外,研究了落锤实验中试样厚度对炸药温升过程及点火阈值的影响,拟合得到了压力峰值和点火阈值随试样厚度变化的估算公式,可以为实验设计提供参考。     

Mechanical characterization of polytetrafluoroethylene polymer using full-field displacement method

The aim of this work is to estimate two important material properties of the polytetrafluoroethylene (PTFE) polymer by means of a single experimental test. The displacement fields around a crack tip are used for estimating the modulus of elasticity (or, Young's modulus) and Poisson's ratio. These parameters are evaluated by fitting linear fracture mechanic expression of displacement fields in the vicinity of the crack, for mode I, to the experimental data. Measurements of these displacements are carried out using digital image correlation (DIC) method. In this way, the experimental procedure is conducted by loading a double-edge-cracked plate specimen. In order to validate the results, two available experimental tests have been performed. The modulus of elasticity is determined by means of the tensile test, using a standard test machine. Moreover, the Poisson's ratio is obtained by measuring lateral compressive and longitudinal extensional strain using DIC method.     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses.In the finite element model,the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method,and the dynamical interaction between both crack faces is modeled using a contact-impact theory.In the simulations,the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis,especially,the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail.Meanwhile,the related ultrasonic infrared images are also obtained experimentally,and the theoretical simulation results are in agreement with that of the experiments.The results show that,by using the theoretical method,a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

In situ observations of compressive behaviour of aluminium foams by local tomography using high-resolution X-rays

We have investigated the compressive behaviour of closed-cell aluminium foams using a high-resolution X-ray CT. The microstructures of cell walls or Plateau borders in the foams were visualized in 3D using the local tomography technique which is a high-resolution CT method to reconstruct a region of interest within a large sample. The shapes and the 3D distribution of micropores, particles, and regions of solute segregation in the foams are evaluated, comparing the cell walls with the Plateau borders. Under compressive loads, the damage behaviour of such microstructures has been observed using an in situ test rig. It is found that the microcracks were mainly initiated from the cell walls and the micropores with large diameters were also damaged. The crack initiation sites are classified from the results. In addition, a method for non-destructive characterization of elastic and plastic deformation in the foams, which is called a 3D microstructure gauge (MG) method, is presented. Thousands of micropores as markers on each load were automatically matched by the information of those volumes and surface areas. The local strain mapping by the MG indicates that the edges of the micropores with large diameters have large strain under compression and this is consistent with the crack analyses.     

Laser cutting silicon-glass double layer wafer with laser induced thermal-crack propagation

This study was aimed at introducing the laser induced thermal-crack propagation (LITP) technology to solve the silicon-glass double layer wafer dicing problems in the packaging procedure of silicon-glass device packaged by WLCSP technology, investigating the feasibility of this idea, and studying the crack propagation process of LITP cutting double layer wafer. In this paper, the physical process of the 1064 nm laser beam interact with the double layer wafer during the cutting process was studied theoretically. A mathematical model consists the volumetric heating source and the surface heating source has been established. The temperature and stress distribution was simulated by using finite element method (FEM) analysis software ABAQUS. The extended finite element method (XFEM) was added to the simulation as the supplementary features to simulate the crack propagation process and the crack propagation profile. The silicon-glass double layer wafer cutting verification experiment under typical parameters was conducted by using the 1064 nm semiconductor laser. The crack propagation profile on the fracture surface was examined by optical microscope and explained from the stress distribution and XFEM status. It was concluded that the quality of the finished fracture surface has been greatly improved, and the experiment results were well supported by the numerical simulation results.     

数字图像相关中基于位移场局部最小二乘拟合的全场应变测量

为了从含噪声的位移场中计算得到可靠的应变场,提出一种基于位移场局部最小二乘拟合的全场应变求解方法。介绍了数字图像相关方法的原理,阐述了基于位移场局部最小二乘拟合的全场应变求解方法,并讨论了计算区域边界、孔洞及裂纹附近区域等情况下的应变计算。对均匀变形和中心带圆孔的薄铝板拉伸实验的计算结果表明,该方法能有效地从原始位移场数据中提取全场应变信息。在均匀变形情况下应选择大的应变计算窗口,计算结果更逼近真值;在非均匀变形情况下,如果位移场中包含较强的噪声,则应选择较大的应变计算窗口,而位移场精度很高时可选择更小的应变计算窗口。     

A simple discrete-element-model of Brazilian test

We present a statistical model which is able to capture some interesting featuresexhibited in the Brazilian test of rock samples. The model is based on elements whichbreak irreversibly when the force experienced by the elements exceed their own loadcapacity. If an element breaks the load capacity of the neighboring elements are decreasedby a certain amount, assuming weakening effect around the defected zone. From the model wenumerically investigate the stress-strain behavior, the strength of the system, how itscales with the system size and also its fluctuation for both uniform and Weibulldistribution of breaking thresholds in the system. To check the validity of ourstatistical model we perform few Brazilian tests on Sandstone and Chalk samples. Thestress-strain curve from model results agree qualitatively well with the lab-test data.Also, the damage profile right at the point when the stress-strain curve reaches itsmaximum is seen to mimic the crack patterns observed in our Brazilian testexperiments.     

Inspection of cracks using laser-induced ultrasound with shadow method: Modeling and validation

A numerical model based on the finite element method (FEM) is developed for crack inspection using laser-induced ultrasonic waves with the shadow method. Before the simulation, the numerical code is verified by comparing the directivity patterns of ultrasonic waves with the experiment. The use of thermo-elastically generated ultrasonic waves for crack detection and evaluation with the shadow method is investigated by simulation. The simulation is validated by comparing with the experiment.     

Crack level estimation approach for planetary gearbox based on simulation signal and GRA

The planetary gearbox is a critical mechanism in helicopter transmission systems. A crack level estimation methodology for planetary gearbox has been devised by integrating a physical model for simulation signal generation, a statistic algorithm for feature selection and a grey relational analysis (GRA) algorithm for damage level estimation. The physical model is used to generate simulation datasets for developing and evaluating the diagnostic scheme that will be further calibrated with real world test data during application. The proposed method was calibrated with historical test data and then validated with real-time test data. The estimation results coincide with the actual test records, showing the effectiveness and accuracy of this method in providing a novel way for more accurate health monitoring and condition prediction.     

Finite strain stress fields near the tip of an interface crack between a soft incompressible elastic material and a rigid substrate

We present a numerical study of finite strain stress fields near the tip of an interface crack between a rigid substrate and an incompressible hyperelastic solid using the finite element method (FEM). The finite element (FE) simulations make use of a remeshing scheme to overcome mesh distortion. Analyses are carried out by assuming that the crack tip is either pinned, i.e., the elastic material is perfectly bonded (no slip) to the rigid substrate, or the crack lies on a frictionless interface. We focus on a material which hardens exponentially. To explore the effect of geometric constraint on the near tip stress fields, simulations are carried out under plane stress and plane strain conditions. For both the frictionless interface and the pinned crack under plane stress deformation, we found that the true stress field directly ahead of the crack tip is dominated by the normal opening stress and the crack face opens up smoothly. This is also true for an interface crack along a frictionless boundary in plane strain deformation. However, for a pinned interface crack under plane strain deformation, the true opening normal stress is found to be lower than the shear stress and the transverse normal stress. Also, the crack opening profile for a pinned crack under plane strain deformation is completely different from those seen in plane stress and in plane strain (frictionless interface). The crack face flips over and the tip angle is almost tangential to the interface. Our results suggest that interface friction can play a very important role in interfacial fracture of soft materials on hard substrates.