非理想爆源产生的爆炸场数值模拟

本文采用有限差分方法TVD(Total Variation Diminishing,A. Harten,1983)格式对非理想爆源在地面上形成的二维轴对称爆炸场进行了数值模拟,数值计算值与试验结果吻合得相当好。     

泄爆过程中二次爆炸的动力学机理研究

在容积为0.00814m3的柱形泄爆容器中，对泄爆现象进行实验研究. 容器内充满当量比为1的甲烷-空气预混气，采用底端中心点火，泄爆压力为230±15kPa. 基于k-ε湍流模型和EBU燃烧模型，利用同位网格的SIMPLE算法，对该现象进行了数值模拟. 实验和计算获得的外轴线上4个测压点的压力曲线和外流场的阴影和数值照片，形象地描述了高压泄爆时外部流场的变化. 数值结果与实验结果基本一致. 根据实验和数值结果，详细地讨论了泄爆过程中二次爆炸产生的动力学机理. 泄爆的初始阶段，在破膜激波的引导下，泄出的未燃气体因欠膨胀在外流场形成稀疏波低压区和悬激波高压区. 高压区可燃气体密度和温度上升，成为高密度的预热区域. 随后，火焰以射流形式从泄爆口泄出，点燃可燃气云. 受湍流等因素的影响，特别在高密度的预热区域，燃烧速率可能迅速增大，从而导致二次爆炸.     

水下爆炸作用下基于声学的爆源子结构输入方法

无限域吸收边界和爆炸波输入是水下爆炸载荷数值模拟的两个关键问题.本文借鉴基于内部子结构的地震波动输入方法和爆源子结构多尺度分析方法,考虑水下爆炸载荷与地震载荷同属于波动问题,提出一种水下爆炸作用下爆源子结构的爆炸波输入方法,该方法首先通过对爆源区域进行水下爆炸波自由场的波场分解,并利用地震波的内部子结构输入方法,将该自由波场运动转换为等效爆炸波输入载荷,从而实现了水下爆炸问题中冲击波的输入;本文采用圆形爆源子结构区域,并采用AUTODYN软件中一维自由场模型计算水下爆炸作用下该区域的自由波场压力时程.进一步,基于连分式近似方法提出了一种模拟无限域水体辐射效应的高精度时域人工边界条件,该吸收边界可设置于离结构和爆源子结构较近处.本文提出的分析方法通过圆形爆源子结构进行水下爆炸载荷的转化,并采用高精度吸收边界大大减少计算区域,既保证了计算精度,又降低了单元数量,具有较高的计算效率和较强的实用性.通过数值算例分析,验证了本文模型与方法的准确性,模拟了水下爆炸作用下的冲击波和气泡脉动阶段测点的压力时程曲线,并研究了圆形结构对水下近场爆炸波散射效应的影响规律.     

铁电体在当代高技术中的应用

回顾了铁电体问世以来，在军事，工业和科研应用领域的发展，着重介绍了目前国际，国内在高技术领域研究和应用趋势，特别在粒子加速器领域。铁电体作为高品质新型脉冲强流电子束源，展示的一系列独特优点，正受到高功率粒子加速器界和阴极电子学界的关注。近几年,随着高功率微波技术的发展，人们利用铁电体的高重复率激励能力，研制出一种触发可靠的高重复率，超低抖动（约70ps)的先进开关，这种开关将能满足超宽带（UWB）及高功率瞬变阵列等技术的精密要求。对铁电体这种新型功能材料的探索和开发应用将会给当代高技术的诸多方面注入新的活力，获得使用传统技术不可能达到的新进展。     

泄爆外流场的可视化

通过实验和数值方法,对泄爆外流场进行了可视化研究。实验中采用YA-16高速阴影系统,拍摄了泄爆外流场的时序阴影照片。基于K-湍流模型和漩涡破碎（eddy dissipation）燃烧模型,利用同位网格SIMPLE算法,对泄爆过程进行了数值模拟。根据计算结果,由计算光学获得泄爆外流场的时序计算阴影图。实验阴影图与计算阴影图较一致,都形象地揭示了泄爆后湍流火焰的发展及二次爆炸的产生和变化过程。     

多爆源云雾爆炸波相互作用的三维数值研究

从三维全N S方程出发,用二阶迎风TVD格式,针对贴地和近地三种不同的爆源位置,研究了两团云雾爆炸后爆炸场的相互作用。计算结果表明:本文的方法可模拟多爆源爆炸后冲击波的传播和爆炸场复杂的波系结构。压力分布可确定爆炸场任意位置的压力时间变化历程。     

黄土中爆炸挤密实验与数值模拟

为研究爆炸挤密加固技术在黄土中的应用,规避在既有公路上进行大规模爆炸挤密现场实验的风险,验证借助计算机软件进行数值模拟的可行性和可靠性,先设计实施了小型爆炸挤密室外实验,再利用室外实验的材料参数和几何尺寸建立与各实验工况相对应的有限元模型,利用ANSYS/LS-DYNA进行数值模拟,通过将爆腔体积、爆后土壤密度和作用于土壤的峰值压应力3个方面对数值模拟结果和实测结果进行对比,验证了将ANSYS/LS-DYNA用于数值模拟爆炸挤密技术加固黄土的可行性和可靠性,并得出上述三个方面的变化规律,可为根据现场路基状况和几何尺寸进行数值模拟提供借鉴和参考。     

一维轴对称杆件爆源模型及其在台阶爆破模拟中的应用

为了解决在实体炮孔建模时需要加密网格、计算量大等问题,提出了一种一维轴对称爆源模型:利用一维线状杆件表述炮孔及炸药,实体单元表述周围岩体,通过杆件节点与实体单元的拓扑关系,将杆件节点上的爆生气体压力施加至周围实体单元上,并根据实体单元的体应变计算出杆件节点处的横截面变化情况,从而实现模拟炸药与周围岩体的相互作用。通过与实体炮孔模型的数值对比分析,发现压力衰减指数为1.25时,一维轴对称爆源模型获得的径向质点峰值振动速度（peak particle velocity, PPV）衰减规律及振动速度时程曲线与实体炮孔模型基本一致,证明了该模型在模拟爆破问题中的精确性。针对混凝土块动态爆破破坏特性的研究,通过与文献对比分析,进一步验证了该模型的正确性。为验证该一维轴对称爆源模型在台阶爆破模拟中的应用,以鞍钢露天铁矿台阶爆破开采为研究对象,建立了5排50炮孔三维台阶爆破概化模型,模拟了爆区内露天边坡的损伤破坏状态。数值计算结果表明,爆区内拉伸破坏为主,并且除了离爆源较近的第一个测点外,其余测点处的振动速度峰值大小及其随距离的变化规律与现场实测数据基本一致,证明了所提爆源模型在三维台阶爆破远场模拟的可行性。     

FAE爆炸波对地面目标作用的三维数值研究

利用二阶迎风TVD格式 ,对多组分全N S方程进行全耦合数值求解 ,研究了近地空中环氧丙烷 (C3H6 O)形成的圆柱形云雾爆炸波和地面作用问题 ,得到了对目标作用前、后冲击波的传播和波系结构的演化以及不同时刻爆源组分的分布。对目标表面计算单元的压力进行面积加权平均 ,得到目标各表面平均超压随时间变化的曲线。压力值的动画显示表明了目标表面所受载荷的大小和分布随时间变化的规律。修改初、边值条件 ,本文的方法和开发的软件可对任意形状云雾爆炸波与目标作用问题进行综合数值分析。     

二级高压驱动阵列弹珠同步弹射微型爆源的研制

针对当前大当量地下爆炸真空室模拟试验中爆源起爆方式高度依赖火药制品等问题，基于地下爆炸相似理论和二级气炮原理，自主研制了二级高压驱动阵列弹珠同步弹射微型爆源装置。装置利用二级高压气体驱动弹珠同步击碎玻璃球壳，释放球内高压气体，以模拟真实爆炸气体生成物的推出。整套爆源装置的发射参数:高压气室充气压力4 MPa，玻璃球壳内残余稳态气体压力约为3 kPa，能够用于0~20 kt当量地下爆炸成坑效应的真空室模拟。爆源适用性试验验证表明，该爆源装置的爆破机制和爆破效果满足大当量地下抛掷爆炸真空室模拟试验的功能需求，且具有较高的安全性、可控性和可操作性，为开展相关模拟试验提供了新的技术方法。     

Viscoelastic properties of polymerlike reverse micelles

A dramatic increase in the viscosity of reverse micellar solutions of lecithin in a variety of organic solvents of up to a factor of 10⁶ upon the addition of a small amount of water can be observed. The formation of viscoelastic solutions can be explained by a water-induced aggregation of lecithin molecules into flexible cylindrical reverse micelles and the subsequent formation of a transient network of entangled micelles. The viscoelastic properties of these solutions are characterized as a function of water content and temperature for different organic solvents by means of dynamic shear viscosity measurements. The results are interpreted by making analogies to the behavior of semidilute polymer solutions and living polymers.Dedicated to Prof. Dr. J. Meissner on the occasion of his 60th birthday.     

The effect of dissolved polymers on the rheological properties of coating colours

Coating colours used for the coating of paper and board consist mainly of a mineral pigment, which is very often clay, a synthetic binder such as a styrenebutadiene latex, dispersion agents and water retention aids. The latter are often water soluble polymers. These polymers have a very strong influence on the rheological properties of the coating colours, both on the strain rate dependence of the apparent viscosity and on the viscoelasticity. The effects of two different grades of carboxymethylcellulose (CMC) and one grade of hydroxyethylcellulose (HEC), on the rheological properties at room temperature of a clay-based coating colour at pH 8, have been investigated. It is concluded that the high values of the dynamic modulus of the colours are due to interactions between the cellulose derivatives and the solid particles, i.e. mainly the clay particles. For HEC this interaction is associated with adsorption of the polymeric molecules on the clay particles. In the case of CMC, the adsorption is strongly retarded by the presence of the dispersant (a polyacrylate salt). It is suggested that the marked elasticity of the CMC-containing colour in addition to a possible polymer adsorption may be due to charge interactions and/or depletion flocculation. The effect of CMC and HEC on the water-retention properties of the colour is also discussed.     

Nonlinear study of the dynamic behavior of a string-beam coupled system under combined excitations

In this paper,the nonlinear dynamic behavior of a string-beam coupled system subjected to external,parametric and tuned excitations is presented.The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system which are described by a set of ordinary differential equations with two degrees of freedom.The case of 1:1 internal resonance between the modes of the beam and string,and the primary and combined resonance for the beam is considered.The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system and obtain approximate solutions up to and including the second-order approximations.All resonance cases are extracted and investigated.Stability of the system is studied using frequency response equations and the phase-plane method.Numerical solutions are carried out and the results are presented graphically and discussed.The effects of the different parameters on both response and stability of the system are investigated.The reported results are compared to the available published work.     

Ultrasonic devulcanization of waste rubbers: Experimentation and modeling

Continuous ultrasonic devulcanization of ground tire rubber (GRT) and styrene-butadiene rubber (SBR) is considered. Experiments are performed under various processing conditions. Two recipes of SBR with different amounts of polysulfidic linkages are utilized. Gel fraction and crosslink density of devulcanized rubbers are measured and a unique relationship between them is established. Die characteristics with and without imposition of ultrasonic waves are determined. Devulcanized samples are revulcanized and mechanical properties are measured. In some cases, properties of revulcanized SBR samples exceeded those of virgin vulcanizates. This is explained based on the presence of a double network in the revulcanized rubber. A modification of acoustic cavitation and flow modeling of ultrasonic devulcanization of SBR and GRT is proposed using a concept of effective viscosity characterizing the flow of vulcanized particles before devulcanization combined with a shear rate, temperature and gel fraction-dependent viscosity of devulcanized rubber. Velocity, shear rate, pressure, and temperature field along with gel fraction, crosslink density and number of bonds broken are simulated. Predicted data on gel fraction, crosslink density, and pressure using the present modification of the model are found to be closer to experimental data then previously reported.Dedicated to the memory of Professor Tasos C. Papanastasiou     

Impact non-pénétrant sur le thorax : influence de la courbure sur la propagation des ondesNonpenetrating impact on the thorax: influence of curvature on wave propagation

Theoretical studies of the propagation of impact waves through the thorax are needed to improve the design of bulletproof jackets and blast protections (Fung in ‘Biomechanics Motions, Flow, Stress, and Growth’, Springer-Verlag, 1990; Cooper et al., J. Trauma 40 (1996) S38–S41). The influence of the weak acoustic coupling at the interface between the thoracic wall and the lung were described in (Grimal et al., C. R. Acad. Sci. IIB 329 (2001) 655–662); in this work, we study, within the frame of elastodynamics and with an approximate analytical method, the effects of the curvature of this interface. Results are given in terms of strain energy for the pressure wave, transmitted or converted. Focalisation of energy in the medium representing the lung is important for curvatures measured in humans. To cite this article: Q. Grimal et al., C. R. Mecanique 330 (2002) 569–574.     

Identification of multi-modal random variables through mixtures of polynomial chaos expansions

A methodology is introduced for the identification of a multi-modal real-valued random variable from a collection of samples. The random variable is seen as a finite mixture of uni-modal random variables. A functional representation of the random variable is used, which can be interpreted as a mixture of polynomial chaos expansions. After a suitable separation of samples into sets of uni-modal samples, the coefficients of the expansion are identified by using an empirical projection technique. This identification procedure allows for a generic representation of a large class of multi-modal random variables with low-order generalized polynomial chaos representations.     

Response of blood flow through an artery under stenotic conditions

This paper presents an analytical study on the behavoiur of blood flow in an artery having a stenosis. This is basically formulated through the use of a suitable mathematical model. The arterial segment under consideration is simulated by an anisotropically elastic cylindrical tube filled with a viscous incompressible fluid representing blood. The analysis is carried out for an artery with mild local narrowing in its lumen forming a stenosis. Particular emphasis has been paid to the effect of the surrounding connective tissues on the motion of the arterial wall. Blood is treated as a Newtonian fluid. The analysis is restricted to propagation of small amplitude harmonic waves, generated due to the flow of blood whose wave length is large compared to the radius of the arterial segment. The effect of the shape of stenosis on the resistance to blood flow has been well illustrated quantitatively through numerical computations of the resulting expressions. A quantitative analysis is also made for the variation of the phase velocity, as well as the velocity of wave propagation and the flow rate, in order to illustrate the applicability of the model.     

The motion equations of rotors,gyrostats and gyrodesics: application to powertrain-vehicle mechanical systems

The kinematics of wheels and rotors is described using a new auxiliary frame called gyrodesic frame or simply gyrodesic. By this, the absolute motion of the wheel becomes a serial composite of two motions: (1) the gyrodesic motion and (2) the wheels eigenmotion (or spin), i.e., the motion relative to the gyrodesic. The eigenmotion is described by an equation called rotor-equation. Gyrodesic coordinates turn out to be a particular useful tool in powertrain- and vehicle-dynamics as well as for general multibody systems. They allow a proper separation of the rotor- from the vehicle-equations and provide a rigorous method of coupling the powertrain-model into full spatial multibody-systems vehicle model. Some common misconceptions regarding this subject are identified and dispelled. The method is generalized to be applicable to the study of motion of general systems of rigid bodies with gyrostats or rotors as subsystems. The usefulness of the formalism is demonstrated by means of an illustrative example of non-trivial nature: the gyrostatic chain. Gyrodesic coordinates lead to a better grasp and deeper understanding of the structure of the dynamic equations of spatial vehicles in particular and of the motion of multibody-systems with rotors in general. The investigation reveals an interesting analogy to concept of parallel transport of vector fields in the sense of Levi-Civita. Dedicated to Prof. J. Wittenburg at the occasion of his 70th birthday.     

Effect of Elevated Velocity of Particles in Groundwater Flow and Its Role in Colloid-facilitated Transport of Radionuclides in Underground Medium

Colloid-facilitated transport of radionuclides by groundwater can increase the level of ecological hazard from radioactive contaminant migration in geological medium. The reason for this is that the migration velocity of a radioactive colloid can be higher than that of radionuclides carried by the groundwater as a solute. On the basis of their field and laboratory studies, a few researchers have concluded that the velocity of the colloid can even exceed the interstitial velocity of the groundwater by a few times. A theoretical analysis of this effect is carried out in this article. The analysis is based on the assumption that the only mechanism responsible for the effect is caused by a redistribution of the colloid in a cross section of the groundwater flow in a representative volume to such domains of the cross section where the local velocity of the flow is higher than the average velocity over the whole cross section. This redistribution can be caused by drift forces arising as a result of the Magnus effect. The influence of these drift forces on a particle’s movement is considered in two extreme cases; that of relatively large and relatively small colloid particles. Particles are considered relatively small if the thermal motion of water molecules exerts a decisive influence on the particles’ movement. Otherwise the particles are considered as being relatively large. It is shown that in the case of relatively large particles this redistribution can be caused by the instability of their movement. The redistribution in the case of relatively small particles can be caused by an influence of the drift forces on characteristics of Brownian motion. It follows from the results of the theoretical analysis that an influence of the drift forces in both cases does not lead to an increase in the particles’ migration velocity at near-horizontal direction of the groundwater flow. Data from experimental studies of the elevated velocity of colloids in porous medium are analyzed. It is shown that some findings of colloidal migration velocity exceeding the interstitial velocity of the groundwater are a result of misinterpretation of experimental results.