应用阵面追踪法对散心爆轰波传播的数值模拟

考虑到爆轰波阵面曲率及化学反应区宽度的作用 ,由修正Hugoniot关系式解析求解了波后状态 ,应用阵面追踪法 (FTM)数值模拟了曲面散心爆轰波的传播。计算给出的散心爆轰波走时及波阵面上物理量都达到了波阵面曲率的一阶精度。     

一点起爆半球形装药装置中散心爆轰波传播的数值模拟研究

运用自行研制的二维流体动力学程序TDY2D,对一点起爆半球壳装药装置中散心爆轰波的传播过程进行数值模拟计算.对冲击波到达炸药和飞层各界面的波形的计算值与实验值进行了对比,两者符合较好,从而验证了该程序的准确性和有效性;对爆轰波的传播过程进行了分析研究,获得了一些散心爆轰波传播的规律性认识,为进一步研究散心爆轰波特性提供了数值分析基础.     

爆轰波与金属的斜相互作用

应用二种截然不同的炸药PBX 95 0 1和T2 ,计算了爆轰波在不同入射角下与金属平板的斜相互作用。在正规反射区 ,计算结果与激波极线理论基本一致 ;应用燃烧模型 ,分别计算了Fe球壳装药JB90 0 3(HE)及JB90 14 (IHE)散心爆轰波的绕射传播 ,计算结果与实验很好地符合。在非正规反射区 ,二维拉氏程序计算结果明显地不同于经典理论结果 ,计算中没有出现Mach反射。计算结果显示 ,毗邻介质影响节点附近的爆轰波阵面形状及爆轰波速度 ;不同的反应率函数计算的节点图像不同。     

底面起爆的圆弧钝感炸药爆轰传播规律的数值模拟研究

用爆轰冲击波动力学的广义几何光学模型ＧＧＯ２Ｄ编码计算了底面起爆的圆弧钝感炸药的爆轰波传播过程,以及同一圆心角上不同半径的爆轰走时,计算结果与实验符合得很好。     

定常爆轰数值模拟中化学反应率与人为粘性的相关性

本文给出了用爆速燃烧、Cochran及Forest Fire三种不同反应率函数对正常爆轰波进行数值模拟时,分别满足的唯象条件。应用Cochran或ForestFire反应率时,适当选取人为粘性,其唯象条件直接给出了数值模拟正常爆轰应使用的空间步长,提供了数值计算结果。     

炸药爆轰产物驱动不可压缩刚粘塑性柱壳的运动

从质量和动量守恒关系出发 ,结合炸药爆轰产物压力剖面的分析和假设 ,导出了不可压缩刚粘塑性柱壳在柱面散心和聚心爆轰产物推动作用下的运动方程。应用该方程式计算爆轰驱动柱壳的运动与实验结果及动力有限元程序DYNA2D的数值模拟结果十分接近 ,说明散心爆轰波场合波后产物流场可用柱面Taylor波的近似解合理描述 ,有限初始半径聚心爆轰波后产物流场可用平面Taylor波近似描述。     

铝粉/空气二维黏性两相爆轰的数值模拟

为了深入研究爆轰波形成和传播的机理与特性,建立了管内铝粉/空气二维黏性两相爆轰过程的数学模型,采用守恒元与求解元方法进行数值计算,并对其物理参数的分布进行了分析。结果表明:管内燃烧转爆轰的初期,压力沿径向变化明显,与壁面碰撞有明显的反射波;燃烧转爆轰的中后期,压力沿轴向变化明显,但径向效应仍不能忽视,碰撞形成的反射波对最终稳定爆轰波的形成影响较大。研究结果同时表明:不仅铝粉颗粒初始半径对爆轰波的形成与传播有一定的影响;而且气体的黏性作用在研究爆轰管内近壁面处流场时不容忽视。研究结果有利于进一步揭示铝粉燃烧转爆轰的机理。     

气相爆轰物理的若干研究进展

爆轰现象的研究已经有一百多年的历史了,爆轰物理的研究取得了许多重要进展.本文从爆轰波的经典理论、胞格爆轰波的多波结构、气相爆轰波形成机理、气相爆轰波传播机制等方面综述了相关的若干研究进展,评述了这些进展的科学性与局限性,并探讨了将来可能的研究方向.这些研究进展主要包括:CJ(Chapman-Jouguet)理论和ZND(Zel'dovich,von Neumann,D?ring)模型、爆轰波多波结构、爆轰胞格特征、直接起爆和爆燃转爆轰过程、热点起爆机制、爆轰波稳定性、扰动爆轰波的传播等.爆轰波是以超声速传播的自持燃烧现象,涉及了激波相互作用、燃烧化学反应、湍流扩散和流动不稳定性等复杂的气动物理过程,相关研究具有重要的学科理论意义.另外,爆轰燃烧具有高效的热化学能释放特点,在先进的热力推进技术方面有着重要的应用背景,因此相关研究也具有重要的工程应用价值.      

提高高马赫数超燃冲压发动机推力的理论方法

斜爆轰发动机和激波诱导燃烧冲压发动机在高马赫数吸气式发动机中具有重要应用前景,但是斜爆轰发动机是否具有足够大的净推力,还是一个未知的问题,因此需要对高马赫数冲压发动机的推进性能以及提高推力的方法进行理论研究.本文主要分为3部分.第1部分理论研究了超燃冲压发动机中的爆燃波和爆轰波的传播特性.保证发动机稳定燃烧是提高推力的...     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

Application of the theory of branching of solutions of nonlinear equations in the estimation of pole perturbance of linear closed-loop systems

When the state and input matrices of a multivariable linear time-invariant system are perturbed, the problem of the estimation of pole perturbance of closed-loop system is considered by making use of the theory of branching of solutions for nonlinear equations.     

两自由度振动系统的斜碰撞分析

研究斜碰撞振动系统动力学的一个关键问题是对系统在碰撞前后的状态进行合理描述和正确计算．针对两弹性体斜碰撞问题，基于瞬间碰撞假设，提出了采用步进冲量来分析和求解斜碰撞前后的状态关系；并以弹簧摆和振子组成的两自由度斜碰撞振动系统为例，具体介绍了该算法如何实现．用解析方法讨论了该系统在斜碰撞过程中可能出现的各种力学现象，将冲量步进算法得到的数值解与解析结果进行对比，取得了完全一致的结果．该数值方法能适应多种斜碰撞问题的计算．     

The distribution of zeroes of solutions of neutral equations

1.IntroductionThestudyofoscillator}'behaviorofneutraldifferentialequationsisarelativealnewfieldandisveryinterestinginapplications.Inrecentyears.theoscillatorytheoryfortheseequationshasbeenextensivelydeveloped.WerefertotherecentbookbyGyoriandLadaSI'I,Baino…     

沙漠地貌演化过程的湍流大涡模拟研究

采用大涡模拟和混合粒径群体沙粒运动算法结合内置边界方法来模拟沙波纹与沙丘的演化和发展,在平坦沙区形成沙波纹的同时,在沙丘上同样也出现了沙波纹的结构,即沙波纹的发展出现了多尺度现象.此外,还分析了粒径对沙波纹与沙丘自身形态演化和发展的影响,以及沙丘对平坦沙区沙波纹形成发展的作用.     

综合管廊早龄期混凝土温度应力仿真模拟

以汉中市综合管廊为工程背景,通过有限元软件,首先模拟分析综合管廊早龄期温度场,随后将温度场结果附加到管廊结构中求解温度应力,从而分析综合管廊早龄期混凝土温度应力与浇筑长度的关系。模拟分析发现:综合管廊在混凝土浇筑早期,温度场沿长度方向变化不显著,且随着浇筑长度的增加,温度场不变;在墙体温度上升阶段,墙体内部沿长度方向的应力主要为压应力,当墙体混凝土内部温度达到最高,混凝土开始降温,结构内变形由膨胀逐渐变为收缩,内力由压应力逐渐变为拉应力,随着时间的变化,墙体内外温差趋于稳定,墙体温度应力也逐渐平稳,此时温度应力主要受外界环境变化影响;当综合管廊浇筑长度在60m范围内,浇筑长度越大,温度应力越大,但当浇筑长度超过60m,温度应力不再增加。     

Identification of viscoelastic characteristics of composite materials on the basis of results of an experimental and theoretical analysis of the dynamic behavior of hemispherical shells

A method is proposed for determining the stiffness and rheological characteristics of composite materials, which is based on minimizing the disagreement between experimental data and results of numerical simulations of deformation of hemispherical shells under explosive loading. The damping characteristics of randomly reinforced polymer materials are analyzed with the use of this method. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 126–133, May–June, 2006.     

On some problems of stability of systems with an infinite number of degrees of freedom

The problem of stability of equilibria of a physical pendulum with a nonstretchable thread attached to it is considered from the standpoint of the Lagrange theorem on stability and its inversion. Specific difficulties which one faces when studying an infinite dimensional mechanical system are discussed. A new approach to the study of stability with respect to two metrics is suggested. The influence of resonant phenomena on the motion of the shortened (linearized) system is considered.

---

Sommario Si considera il problema di stabilità di punti di equilibrio di un pendolo fisico con un filo inestensibile appeso ad esso dal punto di vista del teorema di Lagrange sulla stabilità e sulla sua inversione. Difficoltà specifiche relative allo studio di un sistema meccanico di dimensione infinita sono discusse. Si suggerisce un nuovo metodo per lo studio della stabilità rispetto a due metriche. L'influenza di fenomeni di risonanza sul moto del sistema ridotto (linearizzato) è considerata.

     

Evolution of viscoelastic properties of polyurethane in the course of curing

Simultaneous measurements of the components of dynamic modulus at the set of frequencies were made by Fourier-Transform Mechanical Spectroscopy (FTMS) to investigate the curing process of a segmented polyurethane.The whole process of curing is treated as consisting of three stages. In the first stage, the material can flow and this stage is complete at the gel-point. In the second stage, the curing proceeds in the network and ends with phase separation of aggregating segments (blocks) of polymeric chains. In the third stage, the curing continues in the two-phase system. The object of investigation was the transition from the second to the third stage. It was shown that the position of the maximum of loss tangent depends on frequency and does not coincide with different characteristic points marking the time dependencies of G and G. This is due to the formation of micro-particles of the newly formed phase. The FTMS method is used for separation of the modulus changes into parts related to formation of chemical and physical networks.Special means of treating calorimetric measurements identified a weak maximum heat output at the moment corresponding to the transition point.     

The Formation of Shock Waves of the Equations of Magnetohydrodynamics

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｎｅ＿ｄｉｍｅｎｓｉｏｎａｌＭＨＤｅｑｕａｔｉｏｎｓｕｎｄｅｒＬａｇｒａｎｇｅｃｏｏｒｄｉｎａｔｅａｒｅ τ ｔ- ｕ ｘ=0 , ｕ ｔ Ｐ ｘ 18π ｘ(Ｈ2 ｙ Ｈ2 ｚ) =0 , Ｈｙ ｔ Ｈｙτ ｕ ｘ =0 , Ｈｚ ｔ Ｈｚτ ｕ ｘ=0 , Ｓ ｔ=0 ,( 1 )ｗｈｅｒｅｕ ,τ,Ｈｙ,ＨｚａｎｄＳｔｈｅｖｅｌｏｃｉｔｙ ,ｓｐｅｃｉｆｉｃｖｏｌｕｍｅ,ｔｈｅｃｏｍｐｏｎｅｎｔｓｏｆｍａｇｎｅｔｉｃｆｉｅｌｄｉｎｙ＿ａｎｄｚ＿ｄｉｅｃｔｉｏｎｓａｎｄｔｈｅｅｎｔｒｏｐｙｒｅｓｐｅｃｔｉｖｅｌｙ .Ｔ…