Experimental and numerical study of the fragmentation of expanding warhead casings by using different numerical codes and solution techniques

There has been increasing interest in numerical simulations of fragmentation of expanding warheads in 3D. Accordingly there is a pressure on developers of leading commercial codes, such as LS-DYNA, AUTODYN and IMPETUS Afea, to implement the reliable fracture models and the efficient solution techniques. The applicability of the Johnson–Cook strength and fracture model is evaluated by comparing the fracture behaviour of an expanding steel casing of a warhead with experiments. The numerical codes and different numerical solution techniques, such as Eulerian, Lagrangian, Smooth particle hydrodynamics (SPH), and the corpuscular models recently implemented in IMPETUS Afea are compared. For the same solution techniques and material models we find that the codes give similar results. The SPH technique and the corpuscular technique are superior to the Eulerian technique and the Lagrangian technique (with erosion) when it is applied to materials that have fluid like behaviour such as the explosive and the tracer. The Eulerian technique gives much larger calculation time and both the Lagrangian and Eulerian techniques seem to give less agreement with our measurements. To more correctly simulate the fracture behaviours of the expanding steel casing, we applied that ductility decreases with strain rate. The phenomena may be explained by the realization of adiabatic shear bands. An implemented node splitting algorithm in IMPETUS Afea seems very promising.     

SPH数值模拟中固壁边界的一种处理方法

在光滑粒子流体动力学 (SPH)数值模拟中尝试了一种处理固壁边界的边界力方法 ,给出了一种新的边界力形式。利用SPH方法及边界力方法对水坝坍塌和涌波进入静止水塘这两个自由表面流动问题作了数值模拟。数值模拟结果表明 ,在SPH计算中使用本文所给边界力处理固壁是行之有效的     

Fatigue crack propagation in a helicopter component subjected to impact damage

Damage tolerant methodology is increasingly used in aeronautical components,especially due the fact that the Aviation Regulation requires such an assessment in case an accidental damage occurs.At pre-sent,there is a strong and actual interest in applying such procedures to helicopter components that are subjected to high frequency cyclic loads.In this paper,an investigation on a damaged transmission shaft for a tail rotor transmission of an actual helicopter has been carried out focusing on the fatigue crack propagation.A complete sequence of experimental tests was performed in order to create an actual ballistic damage and to subsequently check the damage tolerant behaviour.The shaft was damaged by oblique ballistic impact and was subsequently subjected to torsional fatigue loading.During the fatigue cycles several cracks propagated from the ballistic damages.Both of these steps(impact and fatigue loading)were also simulated by a complex modelling approach based on Finite Element Models and fracture mechanics theory.The comparison between the experimental and numerical results shows a good agreement but it underlines the need for a very refined modelling technique capable to replicate all the features associated with the damage in order to reliably simulate the subsequent propagation phase.     

Adaptive phasefield modelling of crack propagation in orthotropic functionally graded materials

In this work, we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials (FGMs). A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement. The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency. The implementation is validated against experimental data and other nu-merical experiments on orthotropic materials with different material orientations. The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle. The study is then extended to the analysis of orthotropic FGMs. It is observed that, if the gradation in fracture properties is neglected, the material gradient plays a secondary role, with the fracture behaviour being dominated by the orthotropy of the material. However, when the toughness increases along the crack propagation path, a substantial gain in fracture resistance is observed.     

Characteristics structure analysis on debris cloud in the hypervelocity impact of disk projectile on thin plate

In this paper, the gauge points setting is introduced in the SPH simulation to analyze the debris cloud structure generated by the hypervelocity impact of disk projectile on thin plate. Compared with the experiments, more detailed information of the debris cloud structure can be classified from the numerical simulation. However, due to the solitary dispersion and overlap display of the particles in the SPH simulation, accurate comparison between numerical and experimental results is difficult to be performed. To track the velocity and spatial distribution of the particles in the debris cloud induced from disk and plate, gauge points are locally set in the single-layer profile in the SPH model. By analyzing the gauge points’ spatial coordinate and velocity, the location and velocity of characteristic points in the debris cloud are determined. The boundary of debris cloud is achieved, as well as the fragments distribution outside the main structure of debris cloud.     

Modeling mechanical behaviors of composites with various ratios of matrixeinclusion properties using movable cellular automaton method

Two classes of composite materials are considered: classical metaleceramic composites with reinforcing hard inclusions as well as hard ceramics matrix with soft gel inclusions. Movable cellular automaton method is used for modeling the mechanical behaviors of such different heterogeneous materials. The method is based on particle approach and may be considered as a kind of discrete element method. The main feature of the method is the use of many-body forces of inter-element interaction within the formalism of simply deformable element approximation. It was shown that the strength of reinforcing particles and the width of particle-binder interphase boundaries had determining influence on the service characteristics of metaleceramic composite. In particular, the increasing of strength of carbide inclusions may lead to significant increase in the strength and ultimate strain of composite material. On the example of porous zirconia ceramics it was shown that the change in the mechanical properties of pore surface leads to the corresponding change in effective elastic modulus and strength limit of the ceramic sample. The less is the pore size, the more is this effect. The increase in the elastic properties of pore surface of ceramics may reduce its fracture energy.     

Rayleigh-Taylor不稳定性的SPH模拟

采用SPH方法模拟了流体界面的Rayleigh Taylor不稳定性,给出了界面初始扰动的发展过程,得到了蘑菇状的扰动图像,并对两种不同初始分布的扰动过程进行了比较。分析表明,计算结果合理,SPH方法适合于界面不稳定性的模拟。     

内压作用下复合材料修复裂纹管道的失效分析

为准确评估在内压作用下纤维复合材料修复含裂纹管道的有效性及失效压力,建立复合材料修复后裂纹管道的失效数值模型。该数值模型通过扩展有限元法模拟管道裂纹的扩展,利用cohesive单元模拟胶层的脱粘失效,复合材料的失效通过最大应力失效准则进行判定。通过静水压爆裂试验对所提失效数值模型进行验证,实验结果与数值计算结果具有较好的一致性。失效数值分析结果表明:当内压增大至一定值后,未修复管道的初始裂纹沿轴向及壁厚方向逐渐扩展,进而使得管道内壁单元扩展成真实裂纹,此时真实裂纹贯穿整个壁厚方向,即认为裂纹管道发生爆裂失效,爆裂失效压力随初始裂纹半长呈指数形式下降。复合材料修复裂纹管道的不同修复工况呈现相同的失效模式:在单调递增的内压作用下,管道内表面首先出现黏结裂纹,而后其外表面裂纹张开趋势急剧上升,使得复合材料层应力急剧上升,达到极限强度而失效。且对于不同的初始裂纹尺寸,存在对应的复合材料缠绕层数临界值。     

二维损耗和色散介质光子晶体传输特性的研究

将传输矩阵法(TMM)用于光子晶体传输特性的研究,采用Mur近似吸收边界和周期边界来截断计算区域,计算了以TM模正入射时,二维方格子光子晶体在完整周期结构下的透过率谱;在微波波段制作了光子晶体模型,并设计了实验装置,实验与数值模拟计算结果相符合;另外还研究了有损介质光子晶体、色散和吸收介质光子晶体的传输特性,及其对光子禁带的影响。     

基于边界交叉法的高超声速滑翔飞行器多目标轨迹优化研究

针对高超声速滑翔飞行器复杂约束条件下多目标轨迹设计问题，基于边界交叉法和伪谱法提出了其多目标轨迹优化方法。首先，分析了高超声速滑翔飞行器复杂约束轨迹优化问题的特点，提出了多目标轨迹优化问题。然后，采用边界交叉法和伪谱法将多目标轨迹优化问题转化为一组单目标优化子问题，利用非线性规划算法分别求解。在优化过程中，将已求解子问题的解作为下一个子问题的初始值。利用上述方法求解了最大横程和最小峰值热流轨迹优化问题，仿真结果表明：本文方法能够有效搜索到优化轨迹的Pareto前沿，可以为高超声速滑翔飞行器轨迹设计提供参考。     

抗超高速撞击多层结构复合材料研究

本文介绍了抗超高速撞击多层结构复合材料的设计、制备和实验,并对结果进行了讨论,实验证明,以SiC陶瓷基复合材料为防护层、以空心微珠或多孔变密度复合材料为缓冲层、以编织物复合材料为结构层的多层结构复合材料,具有面密度低、抗撞击速度范围宽、抗撞击性能好的特点。     

Experimental and numerical study of the blast wave decrease using sandwich panel by granular materials core

Among the intrinsic properties of some materials, e.g., foams, porous materials, and granular materials, are their ability to mitigate shock waves. This paper investigated shock wave mitigation by a sandwich panel with a granular core. Numerical simulations and experimental tests were performed using Autodyn hydro-code software and a shock tube, respectively. The smoothed particle hydrodynamics (SPH) method was used to model granular materials. Sawdust and pumice, whose properties were determined by several compression tests, were used as granular materials in the sandwich panel core. These granular materials possess many mechanisms, including compacting (e.g., sawdust) and crushing (e.g., pumice) that mitigate shock/blast wave. The results indicated the ineffectiveness of using a core with low thickness, yet it was demonstrated to be effective with high thickness. Low-thickness pumice yielded better results for wave mitigation. The use of these materials with a core with appropriate core reduces up to 88% of the shock wave. The results of the experiments and numerical simulations were compared, suggesting a good agreement between the two. This indicates the accuracy of simulation and the ability of the SPH method to modeling granular material under shock loading. The effects of grain size and the coefficient of friction between grains have also been investigated using simulation, implying that increasing the grain size and coefficient of friction between grains both reduce overpressure.     

具有模糊系数约束的多目标线性规划

研究了一类具有模糊系数约束的多目标线性规划问题.根据各目标函数的梯度方向来量化目标之间的冲突程度,以此提出了一种确定目标权重的新方法,然后基于惩罚函数运用梯度上升算法求问题的有效解.最后给出了一个数值例子.     

改进的接触算法及其在光滑粒子流体动力学中的应用

采用以黎曼解描述粒子间相互作用的接触算法,模拟一维情况下包含间断的流场。在弱间断的条件下,采用弱波近似的黎曼解来描述粒子间的相互作用;在强间断的条件下,则引入非迭代黎曼解来描述粒子间的相互作用。并引入Taylor展开思想,提高了接触算法在自由边界的计算精度。该算法无需引入人为粘性,程序结构简洁。数值计算表明,改进后的接触算法能很好地描述强间断,并有效改善了原始接触算法在自由面的计算缺陷。     

有初始间隙弹性摩擦接触问题的罚有限元法

针对三维有初始间隙带摩擦的弹性接触问题,提出了一种新的迭代求解方法——罚有限元法。根据接触问题局部非线性的特点,将接触区域的法向非嵌入条件及切向连续性条件作为罚因子引入系统的总势能泛函中,由最小势能原理,得出了求解有初始间隙带摩擦弹性接触问题的罚有限元公式系统。最后通过两个数值算例验证方法的正确性和有效性。     

界面热阻对金刚石/银复合材料导热率的影响

研究了金刚石／银复合材料的导热率与金刚石的含量、粒度等之间的关系,着重讨论了界面热阻对复合材料导热率的影响。并阐述了金刚石表面镀钛,能有效改善金刚石对银的润湿性,降低界面热阻,从而提高金刚石／银复合材料的导热率。     

Effect of nano powder (Al2O3-CaO) addition on the mechanical properties of the polymer blend matrix composite

This work describes the preparation and study of the properties of composite nanoparticles prepared by the sol-gel method which consists of two materials (Al₂O₃-CaO), and study the effect of these nanoparticles on the mechanical behavior of a polymer blend (EP 4% + 96% UPE). The powder was evaluated by X-ray diffraction analysis, scanning electron microscopy analysis (SEM), particle size analysis, and energy dispersive X-ray analysis (EDX). The mechanical behavior of the composite material was assessed by tensile test, bending test and hardness test. The evaluation results of the composite nanoparticles showed good distribution of the chemical composition between aluminum oxide and calcium oxide, smoothness in particles' size at calcination in high and low temperatures, formation of different shapes of nanoparticles and different (kappa and gamma) phases of the Al₂O₃ particles. The results of mechanical behavior tests showed marked improvement in the mechanical properties of the resulted composite material, especially at 1.5%, compared with polymer blend material without nano powder addition. The tensile properties improved about (24 and 14.9) % and bending resistance about (23.5 and 16.8) % and hardness by (25 and 22) % when adding particles of size (63.8 and 68.6) respectively. Therefore, this reflects the efficiency of the proposed method to manufacture the nanocomposite powder and the possibility of using this powder as a strengthening material for the composite materials and using these composite materials in bio applications, especially in the fabrication of artificial limbs.     

用能量法确定层状复合材料疲劳裂纹的扩展方向

含裂纹体疲劳寿命的预测,必须首先了解裂纹扩展路径的全过程,而在层状复合材料中,裂纹扩展的路径非常复杂,不仅可以拐弯,而且可以分叉.因此,首先提出一个简单、适用、统一的裂纹扩展方向的相对强度准则,即应变能释放率比值判据;然后,利用有效解决高梯度问题的数值方法--任意线法,对实测试件的关键时刻(例如,拐弯或进入界面),进行了数值分析;并对典型试件进行了疲劳断裂实验研究.3方面的结果,得到相互的验证.这种一致性,不仅初步证明了这一方法的正确性和适用性,而且为疲劳裂纹扩展全过程的解决,提供了必要的基础.     

A stochastic programming approach for chemotherapy appointment scheduling

Chemotherapy appointment scheduling is a challenging problem due to the uncertainty in premedication and infusion durations. In this paper, we formulate a two‐stage stochastic mixed integer programming model for the chemotherapy appointment scheduling problem under limited availability of nurses and infusion chairs. The objective is to minimize the expected weighted sum of nurse overtime, chair idle time, and patient waiting time. The computational burden to solve real‐life instances of this problem to optimality is significantly high, even in the deterministic case. To overcome this burden, we incorporate valid bounds and symmetry breaking constraints. Progressive hedging algorithm is implemented in order to solve the improved formulation heuristically. We enhance the algorithm through a penalty update method, cycle detection and variable fixing mechanisms, and a linear approximation of the objective function. Using numerical experiments based on real data from a major oncology hospital, we compare our solution approach with several scheduling heuristics from the relevant literature, generate managerial insights related to the impact of the number of nurses and chairs on appointment schedules, and estimate the value of stochastic solution to assess the significance of considering uncertainty.     

多铁性磁电智能复合材料的力学研究进展

介绍了多铁性磁电智能复合材料的概念、分类及发展现状;总结了两类典型多铁性复合材料的力学问题（包括波动与振动问题、夹杂问题、断裂问题等）的研究进展;分析了层状多铁性复合材料力学行为的特殊性;基于应变介导的磁电耦合机制,阐述了界面断裂力学研究对于层状多铁性复合材料研制的重要性和必要性;最后指出了该领域下一步的研究方向。