爆轰波透射孔栅形成的高速爆燃波的结构和行为

以实验为主,研究光滑直管中乙炔氧气爆轰波透射孔栅形成的高速爆燃波结构和行为。实验以基于纹影平台的高速转鼓摄影记录孔栅下游近场内初始爆燃波的结构和发展,并以压力传感器跟踪其后继走向。研究发现,高速爆燃波是前驱激波和火焰的组合结构,随着初始压力的提升,它分别表现为紧随于前驱激波的层流和湍流燃烧火焰。入射爆轰胞格尺度小于或与孔栅扰动尺度相当时,湍流燃烧在下游近场迅速形成;层流结构的爆燃波通常无法抵制背景稀疏波而走向衰弱,而湍流燃烧结构可发生加速和向爆轰的转捩;他们之间存在一个不稳定的临界状态,高速爆燃波得以以0.5~0.6倍CJ爆轰速度传播较长距离,这一状态对应于双间断Rankine-Hugoniot关系的等容燃烧解。     

可爆性气体爆炸极限和爆燃转变成爆轰的研究

我们成功地研制了长为6.6m、内径为100mm的柱形爆炸激波管。利用此激波管我们研究了氢气和氧气混合物的爆炸特性。研究表明:氢气和氧气混合物的可爆性极限为25％H2至84％H2（体积比）;混合物起爆的临界初始压力P0c近似与混合物的浓度无关。我们还用石英传感器技术测量了混合物从爆燃转变为爆轰的过程（DDT）,确定了爆炸速度、爆炸压力与混合物初始压力P0的关系。爆燃波速度D和压力P快速增加,并且DG0（声速）时,爆燃波能转变成爆轰。爆燃波速度D和压力P衰减时,爆燃波将熄灭。氢氧混合物浓度接近上限时,爆燃向爆轰转变时的爆轰呈过激励状态,然后逐渐趋于正常爆轰。C-J理论可近似预估气体爆轰参数。     

两种不同气体中的高速爆燃波及其向爆轰的转变

对2C2H2 5O2及2C2H2 5O2 80%Ar两种可燃混合气体中的高速爆燃波及其向爆轰的转变过程进行实验研究.高速爆燃波由孔栅干涉爆轰波的方法直接生成,观测手段则以高速转鼓摄影获取孔栅近场流场x-t纹影图,以传感器追踪波面的后继发展.研究发现,两种气体中的爆燃波具有迥异的特性.前者燃烧波面在较低初压条件下为层流结构,而较高初压下为湍流结构,向爆轰转变点可以延伸至下游较长距离;后者在不同初压条件下燃烧波面无明显差异,爆轰的再次形成只能在孔栅下游近场内建立.两种气体中高速爆燃波的维持和爆轰转变过程均非纯粹激波压缩所致,湍流输运在其中起着必不可少的作用.分析显示,激波压缩效应对纯氧炔气体的高速爆燃和DDT贡献较小,湍流输运占主导地位;而氩气稀释气体较为稳定,缺乏自行衍生剧烈湍流燃烧的能力,因而激波压缩和外界扰动对其高速爆燃传播和爆轰转变起十分重要的作用.     

爆轰波在突扩通道中传播的数值模拟研究

建立了描述甲烷 空气混合物爆轰波传播的单步化学反应爆轰模型 ,通过数值模拟研究了在二维突扩通道中爆轰波的强度变化和各种波行为。结果表明 :爆轰波在进入突扩通道初始阶段的衍射使爆轰波局部向爆燃转变 ;爆炸波在壁面发生马赫反射形成的高温高压区域将直接诱导自持爆轰波的重新形成。     

大型水平爆轰管中悬浮铝粉爆炸过程的实验研究

铝粉的燃烧与爆轰性能是粉尘爆炸领域研究的热点.利用长29.6m,内径199mm,配有40套喷粉扬尘系统的水平爆轰管,在40J电火花点火条件下,实现了悬浮铝粉-空气混和物火焰加速、爆燃、爆轰及其转捩过程,测得了爆炸波传播过程中的压力信号,并且观察到了爆轰波的稳定传播现象.实验结果表明,当铝粉浓度为300 g/m3时,在距离点火端10.15m(长径比L/D=51)处发生了DDT,测得的爆轰波传播过程中管内的最大爆速为1840m/s,最大峰值超压为10.5MPa.铝粉尘爆炸波在爆轰管内的传播过程可分为爆燃段、爆燃转爆轰(DDT)、爆轰增强以及稳态爆轰四个阶段.     

甲烷-空气爆轰波传播的数值模拟

研究了甲烷-空气混合物爆轰波在二维突扩通道中传播的各种复杂行为。结果表明：爆轰波在进入突扩通道的初始阶段有局部向爆燃转变的现象；爆炸波在壁面发生马赫反射形成的高温高压区域将直接驱动自持爆轰波的重新形成。     

气体爆轰波在声学吸收壁下游的再加强过程

实验研究了爆轰波在声学吸收段下游发生的再恢复过程。利用声学吸收段首先将一个稳定的爆轰波衰减,于是产生一个非稳定爆轰波,速度降为 45% ~60?J。该非稳定爆轰波传出声学吸收段后再次加强,并在距声学吸收段出口端一定距离处突然发展成超驱动爆轰。实验利用烟迹技术捕捉了爆轰波从发生衰减到再次加强的全过程,获得了爆轰波再加强距离随初压变化的实验规律。由此,本文建立了再加强段内加速传播火焰的理论模型,称为激波 -爆燃波复合态,并解释了火焰加速和超驱动爆轰形成的机理。     

气相爆轰波绕射流场显示研究

采用基于红宝石激光器（波长694.3 nm）的纹影系统,对气相爆轰波绕射进行了初步的流场显示研究。采用单色激光和合适半带宽（15 nm）的滤光片,有效地消除了爆轰波自发光对流场显示的影响。合理设置激光器同步控制系统的触发延时,得到了序列的爆轰波阵面纹影照片。结果表明:图像清晰地显示了爆轰波阵面的诱导激波、横波及化学反应区。当爆轰波在左尖点处绕射,受稀疏波作用,诱导激波与化学反应区明显分离,导致爆轰波衰减为爆燃。分离的诱导激波和折皱的化学反应区在纹影图上清晰可见。诱导激波在垂直支管右壁面反射,诱导二次起爆。畸变爆轰波在水平和垂直支管中均发生马赫反射。提高初压,爆轰波受分叉口几何属性的影响减小,畸变爆轰波在水平和垂直支管下游较易恢复为自持爆轰波。     

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

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

压装TNT的SDDT特性研究

本文报告厂将水箱技术同改型隔板试验结合起来研究凝聚炸药由冲击波经爆燃向爆轰转变(SDDT)的方法。实测了压装TNT(装药密度1.615g/cm~3)的SDDT曲线,并确定了其燃烧阈值和爆轰阈值(分别为2.1GPa和2.5-3.0GPa)。按照流体动力学特征线法计算了不同强度起爆冲击波作用下在炸药中形成的一系列冲击波、爆燃波和爆轰波的剖面,这些剖面反映了炸药SDDT的基本特征。所取得的信息对于炸药的实际应用和对冲击波起爆的深入研究均有价值。     

骨陷窝-骨细胞形状和方向对骨单元内液体流动行为的影响

骨组织内的流体流动不仅为骨细胞的生存提供了充足营养供应及代谢物排放途径,也在骨重建过程中起到关键作用. 为了更精确地阐明骨内液体流动的具体形式,这项研究利用骨陷窝-骨细胞的密度,形态和方向等参数来计算骨单元内液体的流动行为. 首先,计算出不同形状和方向的骨陷窝周围骨小管的数量及分布情况,其次利用算出的参数以及骨组织其他微结构数据来估计骨组织的渗透率和孔隙率等参数,最后根据计算所得的参数建立骨单元的多孔弹性力学有限元模型,并分析了在轴向位移载荷作用下骨陷窝形状和方向对骨单元内液体渗流行为的影响. 结果表明,在所研究的参数范围内不同骨单元模型的相同区域上,骨陷窝形状影响下的骨单元最大压力和流速比最小的分别增加了86%和18%;骨陷窝方向影响下的最大压力和流速比最小的分别增加了125%和56%. 伸长形骨陷窝对单个骨单元局部压力的影响远大于扁平形和圆形骨陷窝. 骨陷窝从0°绕$x$轴旋转到90°过程中压力是逐渐降低的,且30°,45°和60°的模型对骨单元内局部流速有显著影响. 该模型表示骨陷窝的形状和方向以及骨小管的三维分布对骨单元内液体压力和流速幅值及沿不同方向的流动差异有显著的影响. 这项研究将有助于精确量化描述骨内液体的流体行为.      

半圆柱壳挠曲问题解析解精确性的研究

壳体力学已于上世纪由多位专家发展成熟,其中简支柱壳挠曲问题采用改进莱维解法的三角级数法解出,但是其解法复杂,手算难以完成．为讨论其结果的精确性,通过编写运行基于MATLAB的运算程序导出实例化解析解,与基于力学基本理论的推想假设对比,再引入有限元计算结果进行比较研究．最终发现,理论解析解应力和位移具有分布形式大致准确性,但仍存在不容忽视的细节与局部性问题．研究表明,理论解法工程意义有限,结果尚需改进．     

Viscoplasticity of agglomerated suspensions

This work is a theoretical study on the effects of agglomeration on the fluidity and plasticity of a suspension of neutrally buoyant particles in a Newtonian fluid. The dynamics of a cluster of permanently attached spherical particles in a simple shear field is analyzed. The viscous and plastic components of the drag force acting on each of the agglomerated particles is then calculated and found to depend on the size of the individual particle unit, its location being relative to the center of the cluster and the material properties of the engulfing fluid. This information in conjunction with the knowledge of the interparticle cohesive forces is used to establish criteria for the agglomerate size reduction during dispersive mixing. From the kinematics of the cluster movement and the forces acting on its particulate components the rate of energy dissipation is calculated and utilized to evaluate the viscosity and yield stress of the suspension. These rheological parameters depend on the volume fraction and architecture of the agglomerate, the number of fused particles per cluster, and the viscosity of the suspending fluid. The analysis is also extended to include the case of polydispersity of agglomerate sizes.     

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.     

复合柔性结构全局模态函数提取与状态空间模型构建

随着航空航天等领域中实际工程结构的大型化和柔性化,结构的非线性振动和主动振动控制问题越来越凸显.分析和处理此类结构出现的复杂振动问题的关键在于建立系统的非线性动力学模型与状态空间模型.对于由柔性部件、刚体、连接部件构成的复合柔性结构,由于各部件之间的振动耦合效应,单个柔性部件在悬臂、简支和自由等静定边界下的模态与结构的真实模态有较大差异.为此,本文提出复合柔性结构全局模态的解析提取方法,通过全局模态离散得到系统非线性动力学模型,从而构建状态空间模型.该方法采用笛卡尔坐标描述系统的运动,建立系统的运动方程;结合描述柔性部件的偏微分方程、刚体的常微分运动方程、连接界面处力、力矩、位移和转角的匹配条件以及系统的边界条件,利用分离变量法给出统一形式的频率方程,获取系统的固有频率和解析函数表征的全局模态.这里提出的全局模态提取方法不仅便于复合柔性结构固有频率和全局模态的参数化分析,而且为建立复合柔性结构低维非线性动力学模型和状态空间模型提供了有效的途径,对于推进这类结构的非线性动力学分析与主动振动控制研究具有重要意义.      

STEADY-STATE RESPONSE OF A TIMOSHENKO BEAM ON AN ELASTIC HALF-SPACE UNDER A MOVING LOAD

By introducing the equivalent stiffness of an elastic half-space interacting with a Timoshenko beam, the displacement solution of the beam resting on an elastic half-space subjected to a moving load is presented. Based on the relative relation of wave velocities of the half-space and the beam, four cases with the combination of different parameters of the half-space and the beam, the system of soft beam and hard half-space, the system of sub-soft beam and hard half-space, the system of sub-hard beam and soft half-space, and the system of hard beam and soft half-space are considered. The critical velocities of the moving load are studied using dispersion curves. It is found that critical velocities of the moving load on the Timoshenko beam depend on the relative relation of wave velocities of the half-space and the beam. The Rayleigh wave velocity in the half-space is always a critical velocity and the response of the system will be infinite when the load velocity reaches it. For the system of soft beam and hard half-space, wave velocities of the beam are also critical velocities. Besides the shear wave velocity of the beam, there is an additional minimum critical velocity for the system of sub-soft beam and hard half-space. While for systems of (sub-) hard beams and soft half-space, wave velocities of the beam are no longer critical ones. Comparison with the Euler-Bernoulli beam shows that the critical velocities and response of the two types of beams are much different for the system of (sub-) soft beam and hard half-space but are similar to each other for the system of (sub-) hard beam and soft half space. The largest displacement of the beam is almost at the location of the load and the displacement along the beam is almost symmetrical if the load velocity is smaller than the minimum critical velocity (the shear wave velocity of the beam for the system of soft beam and hard half-space). The largest displacement of the beam shifts behind the load and the asymmetry of the displacement along the beam increases with the increase of the load velocity due to the damping and wave radiation. The displacement of the beam at the front of the load is very small if the load velocity is larger than the largest wave velocity of the beam and the half space. The results of the present study provide attractive theoretical and practical references for the analysis of ground vibration induced by the high-speed train.     

THE HISTORY OF DEPARTMENT OF WATER CONSERVANCY OF PEIYANG UNIVERSITY 1)

本文基于部分历史史实,对北洋大学(今天津大学)水利系的发展进行了简要的梳理,并扼要地回顾了那些年,北洋大学水利那些事儿、那些人。重点从北洋大学-天津大学水利学科、水利系的历史延续及其在中国水利高等教育中的地位、中国第一水工试验所及 其发展过程、天津大学水利馆及知名教授和杰出校友的相关故事出发,以期较为全面地展示北洋大学-天津大学水利系125年的历 史传承和现代化发展成就,激励新一代天大水利人,续写北洋雄风,并深刻缅怀余元君等已故校友。     

分数导数粘弹性模型描述的土中管桩竖向振动的频域解

在分数导数粘弹性本构模型的基础上综合考虑桩周土和桩芯土的平衡方程和几何方程建立了桩周土和桩芯土的竖向运动的控制方程．在频率域内利用分离变量法和分数导数的性质求解了桩周土和桩芯土竖向振动控制方程．考虑管桩与桩周土、管桩与桩芯土的边界连续性条件以及三角函数的正交性得到了分数导数粘弹性模型描述的土中管桩的竖向振动,通过数值分析研究了管桩和土体模型参数和几何参数对管桩的桩顶复刚度的影响规律．结果显示：桩芯土本构模型的分数导数的阶数对管桩竖向振动的影响较桩周土本构模型的阶数要小,且与频率有一定关系;桩芯土与桩周土的模型参数比τ1 和τ2 对等效阻尼的影响较对刚度因子的影响要大;管桩桩周和桩芯的直径比d 越小,管桩复刚度的实部和虚部就越大;土体力学性能对管桩竖向振动的影响要比管桩桩身力学性能的影响小．     

受压油管屈曲变形与内外层杆管接触分析

针对机械采油工程中杆管偏磨问题, 考虑直井内受压油管屈曲变形与套管、 抽油杆的双层接触, 以及油管和抽油杆的变截面结构和环空间隙, 采用有限元结构稳定性分 析方法对受压油管进行失稳计算, 将求得的多阶失稳波形叠加值作为油管柱的初始位移, 再 采用大位移和接触非线性分析方法模拟受压油管屈曲变形, 以及与套管和抽油杆间的接触摩 擦, 为套管内受压油管屈曲变形分析提供了一种计算方法.     

基于数字图像的土、岩和混凝土内部结构定量分析和力学数值计算的研究进展

土、岩及混凝土是一种内部结构极其复杂的材料,内部细观组成和结构决定了它们在外力作用下内部应力和应变等物理场的分布状态,在很大程度上控制了它们的宏观力学响应和破坏机理和过程。自1995年以来,数字图像处理技术作为一种材料细观空间结构及几何形态的精确量测和数字表述手段快速地被应用于沥青混凝土、水泥混凝土、土和岩体材料细观结构定量分析中来。本文在研究相关文献的基础上,分析和研究了数字图像技术在土、岩及混凝土内部结构定量分析和力学数值计算中的发展过程和研究成果,内容包括:单个骨料颗粒形态特征定量分析;土、岩和混凝土内部结构定量分析;岩体结构定量分析中的应用;基于数字图像处理的土、岩和混凝土细观结构力学计算。在此基础上,本文进一步对数字图像技术在岩土领域中的应用进展及未来的发展趋势进行了探讨,认为数字图像作为一种空间分布测量手段和数字表述方法使得它具有巨大的潜力,特别是在岩土领域的纵深发展过程中,它可能实现考虑细微观结构性的岩土力学和工程分析和预测的、新的方法和理论。