湿式发射对同心筒瞬态热冲击的影响机理

针对路基同心筒自力发射整体热环境恶劣的问题,依托弹性变形和域动分层相结合的动网格技术,基于均质多相流理论并耦合液态水专用汽化求解程序,建立在发射筒底部注水的三维气液两相流体动力学模型;以火箭发动机自由射流注水实验为基础,验证汽化程序三维计算的可靠性与有效性;通过瞬态数值计算,讨论筒底注水角度对导弹、内外筒热环境和导弹载荷特性的影响规律。分析表明:发射筒内发生了显著的汽化反应;导弹及发射系统总体热环境得到了显著改善,实现了发射系统持续降温的目的;在筒底注水后,弹底的附加推力及火箭发动机的推力有一定增加,随着注水量的减少,注水对导弹载荷的影响越来越弱。     

水下爆炸气泡运动的理论研究

在适当深度的无黏、无旋的流体中对水下爆炸气泡运动特性进行理论研究。综合运用势流理论、能量方程以及拉格朗日方程建立气泡在不可压缩流体中的运动方程。并以此为基础,考虑重力、浮力以及阻力等多种因素对气泡运动特性的影响,通过引入新的边界积分方程,结合分析力学中完整非保守系统的Hamilton原理建立气泡在可压缩流体中的运动微分方程,并对微分方程进行求解。将方程的数值解与MSC.DYTRAN非线性有限元软件的计算结果以及经验公式进行对比,方程数值解与二者都具有较好的一致性。结果表明,基于非保守系统可压缩流体建立的气泡运动方程正确、可行,相关的理论研究和计算具有一定参考价值。     

导弹水下发射燃气泡计算

采用一维非定常气流场模型和轴对称理想水流场模型，对水下发射导弹的尾部非定常燃气泡内外流场进行了耦合数值求解。考虑了高温燃气与水介质之间的传热、汽化等对泡内气体流动的影响，数值模拟了燃气泡的生长和脱落过程，揭示了燃气泡中气体流动的基本规律，计算了燃气泡发展对导弹所受水动力的影响。     

喷管对水下爆轰气泡形态和压力特征影响的实验研究

针对具有不同类型喷管的爆轰管在水下爆轰中形成的燃气射流问题,搭建了爆轰实验平台,研究了单次爆轰过程中尾部喷管对水下气泡形态与压力特征的影响。采用数字粒子图像测速技术对高速摄影机拍摄得到的气泡脉动图片进行流场可视化分析,得到各喷管工况下的气泡速度场。为了确认爆轰管内是否形成稳定爆轰波,并观察爆轰波在气液两相界面上的透反射特性,爆轰管尾部安装有2个动态压力传感器,同时在距离喷管一定距离处设置一个水下爆炸传感器,以监测水中传播的压力波。结果表明：扩张喷管工况下的气泡脉动过程与直喷管工况基本一致,但扩张喷管提高了燃气射流速度,气泡膨胀体积更大;因为燃气射流的持续性,收敛喷管工况下的气泡脉动过程具有明显差异,气泡膨胀体积较小,但气泡二次脉动时长相较于一次脉动时长衰减更小;扩张喷管提高了气泡脉动强度,扩张喷管工况下的气泡脉动压力与透射冲击波压力远大于直喷管工况下的气泡脉动压力与透射冲击波压力;收敛喷管工况下的气泡脉动压力与透射冲击波压力都较小,但收敛喷管燃气射流的持续性减缓了气泡脉动压力的衰减速度。相比于直喷管,扩张喷管工况下的气泡脉动时间、气泡脉动压力与透射冲击波压力都更大。收敛喷管工况下的气泡...     

新型路基同心筒热冲击机理与热环境影响因子

针对新型同心筒自力发射高速热冲击载荷下热环境评估与影响因子决策问题,结合弹性变形和域动分层结合的动网格技术,求解了二维轴对称Navier-Stokes方程,分析了新型路基同心筒流场机理与热冲击特性,并确定了热环境评价指标;通过建立以优化拉丁超立方试验设计和径向基神经网络为理论基础的近似数学模型,解决了CFD自动建模困难、计算量大的难点;结合径向基神经网络训练方法,对导弹热环境的影响因子进行了智能决策研究。分析表明:倒吸进入新型同心筒内筒的低温气体有力改善了同心筒热环境;建立的近似模型精度较高,满足工程需求;对导弹热环境的影响因子从大到小依次为筒底导流板直径、筒底导流板长度、导流器高度;为导弹热环境多学科优化设计提供参考。     

烟火药水下燃烧气泡的实验研究

为探讨烟火药水下燃烧的声辐射机理,基于气泡的瞬时性,采用高速摄影和录像的方法对烟火药水下燃烧气泡的整个物理过程进行了试验研究.发现烟火药水下燃烧过程可分为气泡形成、上浮、破裂成气泡云以及气泡云的水面破裂喷涌四个阶段.随着燃烧速度的增加,气泡演变更加剧烈,出现了合并现象,同时所辐射的噪声也随之增大.气泡上升的过程中并不是完整的球形,且总体上由于气泡体积的变化受力不均而呈曲线形上升.     

装药量及水深对水下爆炸气泡动态特性的影响

以气泡体积加速度模型为基础研究水下爆炸气泡运动的初始条件,采用MSC.DYTRAN 非线性 有限元软件,结合开发的定义流场初始条件与边界条件的子程序,研究水下爆炸气泡运动特性,包括气泡的脉 动、坍塌以及射流等运动特性,并将气泡脉动体积计算结果与实验及边界积分方法计算结果进行对比,验证了 有限元模型的正确性与有效性。以此为基础,得到初始水深、装药量与气泡的脉动体积、最大半径、周期以及 射流速度之间的关系,计算结果与经验公式具有较好的一致性。得到一些有规律性的曲线,可为相关水下爆 炸气泡动态特性研究提供参考。     

水下爆炸气泡脉动测量及分析

水下爆炸气泡破坏效应是水中兵器的重要毁伤模式之一。为研究水下爆炸气泡脉动现象,建立了小当量水下爆炸实验系统,并进行了爆炸当量分别为0.125g、1.0g、3.375g和8g TNT的水下爆炸实验。采用球形PETN装药并中心起爆,产生球形对称的气泡和冲击波载荷,并利用高速摄像系统记录水下爆炸气泡脉动过程,以及布置压力传感器测量水中冲击波压力。实验获得了清晰的水下爆炸气泡脉动过程图像,得到了冲击波和气泡脉动压力曲线。对数字化图像进行判读,得到气泡脉动直径和周期。另外根据冲击波曲线测量了气泡脉动周期,对比分析了气泡脉动相关参数。结果表明,高速摄像数据测量的气泡直径与经验公式较接近,高速摄像测量的气泡周期与冲击波曲线测量的气泡脉动周期以及经验公式结果具有较好的一致性。本文提出的实验技术安全、经济、可靠,气泡脉动参数判读精确,满足水下爆炸气泡脉动研究需求。     

Interaction of underwater explosion bubble with complex elastic-plastic structure

Considering the elastic-plasticity of the structure, the combination of boundary element method and finite element method （FEM） is employed to present the calculation method for solving the complex coupled dynamic problem of bubble, elastic-plastic structure and the free surface, and the complete three-dimensional calculation program is developed as well. The error between the calculated result and the experimental result is within 10%. Taking a surface ship for example, the three-dimensional calculation program is extended to engineering filed. By employing the program, the response of the ship under the bubble loading is analyzed. From the stress-time history curves of typical elements of the structure, it can be seen that the pressure reaches its maximum when the bubble collapses and this validates that the pressure generated by the bubble collapse and the jet can cause serious damage on the ship structure. From the dynamic process of the interaction between the three-dimensional bubble and the ship, the low order vertical mode of the ship is provoked and the ship presents whip-shaped motion. And the ship does elevation and subsidence movement with the expansion and shrinkage of the bubble. Some rules and conclusions which can be applied to the engineering problems are obtained from the analysis in this paper.     

水下爆炸气泡射流载荷计算的新方法研究

由于存在强冲击、多相界面复杂运动等强非线性效应,目前对于水下爆炸气泡运动的计算方法无法给出较为可信的水射流运动特征及其载荷形式.本文基于多相可压缩流体的five-equation计算模型,引入界面函数和界面密度压缩技术,采用5阶WENO重构与HLLC近似Riemann求解器进行空间离散,时间离散采用3阶TVD Rung...     

舱段结构在气泡射流作用下的毁伤效果

为研究水下爆炸气泡载荷对船舶局部结构的毁伤效果,设计建造了实尺度舱段模型,并进行水下爆炸实验,测量了水中压力、气泡脉动周期、舱段外板结构的动应变与塑性变形。分析了不同冲击因子作用下舱段结构的动态响应,并结合水中压力测量结果,探讨了气泡射流的成因。在此基础之上,分析了气泡射流载荷的影响范围以及外板在不同载荷作用下的毁伤模式。实验结果表明,相比冲击波载荷,在某些情况下射流冲击载荷引起的结构响应更剧烈,毁伤效果更严重,因此在中近场水下爆炸作用下,气泡射流冲击载荷对船体结构的影响不能忽视。     

近边界三维水下爆炸气泡动态特性研究

模拟了近壁、近自由面的水下爆炸气泡的非线性动态特性,假定水下爆炸气泡脉动阶段的流场是无旋、不可压缩的,采用高阶曲面三角形单元离散三维气泡表面,用边界积分法求解气泡的运动,在计算奇异积分时通过重新构造双层位势的主值积分消除双层奇异积分的奇异性,得到更精确的结果,并通过合理的加权方法精确的求解边界面上各节点的真实速度,结合弹性网格技术（elastic mesh technique, EMT）得到优化速度,在整个模拟过程中不需要采用数值光顺。将本文的三维模型与轴对称模型进行的对比分析表明,两种模型计算结果吻合很好,并用三维模型模拟了气泡与自由表面及圆筒的相互作用,水下爆炸气泡在自由表面及圆筒的联合作用下呈现出强非线性。     

Numerical study of bubble and particle motion in a turbulent boundary layer using proper orthogonal decomposition

We have studied the motion of bubbles and particles in the near-wall region of a turbulent boundary layer, to investigate the influence of the unsteady turbulent structure. The velocity field was computed using Proper Orthogonal Decomposition (POD), and the trajectories of bubbles and particle have been computed by integrating their equation of motion. We have used this to investigate the roles, and the relative importance, of the different forces acting on bubbles and particles, We find that the unsteady turbulent structure plays an important role in the preferential accumulation of bubbles and particles. The accumulation of bubbles depends on a rather complicated interaction between the pressure gradient and the lift force; neither is sufficient, acting on its own, to explain the strong accumulation observed when they act together.     

Advances in the simulation of multi‐fluid flows with the particle finite element method. Application to bubble dynamics

In this work, we extend the Particle Finite Element Method (PFEM) to multi‐fluid flow problems with the aim of exploiting the fact that Lagrangian methods are specially well suited for tracking interfaces. We develop a numerical scheme able to deal with large jumps in the physical properties, included surface tension, and able to accurately represent all types of discontinuities in the flow variables. The scheme is based on decoupling the velocity and pressure variables through a pressure segregation method that takes into account the interface conditions. The interface is defined to be aligned with the moving mesh, so that it remains sharp along time, and pressure degrees of freedom are duplicated at the interface nodes to represent the discontinuity of this variable due to surface tension and variable viscosity. Furthermore, the mesh is refined in the vicinity of the interface to improve the accuracy and the efficiency of the computations. We apply the resulting scheme to the benchmark problem of a two‐dimensional bubble rising in a liquid column presented in Hysing et al. (International Journal for Numerical Methods in Fluids 2009; 60 : 1259–1288), and propose two breakup and coalescence problems to assess the ability of a multi‐fluid code to model topology changes. Copyright © 2010 John Wiley & Sons, Ltd.     

On the sensitive dynamical system and the transition from the apparently deterministic process to the completely random process

The detailed analysis of the dynamical process of coin tossing is made. Through calculations, it is illustrated how and why the result is extremely sensitive to the initial conditions. It is also shown that, as the initial height of the mass center of the coin increases, the final configuration, i.e. head or tail, becomes more and more sensitive to the initial parameters (the initial velocity angular velocity, and the initial orientation), the coefficient of the air drag, and the energy absorption factor of the surface on which the coin bounces. If we keep the head upward initially but allow a small range for the change of some other initial parameters, the frequency that the final configuration is head, would be 1 if the initial height h of the mass center is sufficiently small, and would be clo to 1/2 if h is sufficiently large. An interesting question is how this frequency changes continuously from 1 to 1/2 as h increases. Detailed calculations show that such a transition is very similar to the transition from laminar to turbulent flows. A basic difference between the transition stage and the completely random stage is indicated: In the completely random stage, the deterministic process of the individual case is extremely sensitive to the initial conditions and the dynamical parameters, out the statistical properties of the ensemble are insensitive to the small changes of the initial conditions and the dynamical parameters. On the contrary, in the transition stage, both the deterministic process of the individual case and the statistical properties of the ensemble are sensitive to the initial conditions and the dynamical parameters. The mechanism for this feature of the transition stage is the existence of the long-train structure in the parameter space. The illuminations of this analysis on some other random phenomena are discussed.     

Understanding the effect of boundary conditions on damage identification process when using non-linear elastic wave spectroscopy methods

This paper presents an experimental campaign aimed at understanding the limitations and capabilities of non-linear elastic wave spectroscopy (NEWS) non-destructive technique (NDT) methods in the presence of variable boundary conditions. In particular, the objective was to understand if the contact between the structures under investigation and the clamps used to hold the structures could generate non-classical non-linear effects that could affect the damage detection process by producing false-positive indications of defects presence.Two different techniques were analysed with varying degree clamping torque. The first approach evaluates the resonance frequency shift as a function of the external load amplitude, and it is called non-linear resonant ultrasound spectroscopy (NRUS). The second method used, called non-linear wave modulation spectroscopy (NWMS), monitors the generation of sidebands and harmonics when the structure is excited by a double tone external load.The results showed that the non-classical hysteretical non-linear effects were dependent on the boundary conditions, highlighted by the presence of resonance shift and harmonics and sidebands in an undamaged sample. This research shows that more work is needed to demonstrate the effectiveness of the methods and the ease of implementation in a structural health monitoring system and further research studies and methodology development are needed to discern non-classical non-linear effects generated by contacts between mating parts (clamps and sample) from that generated due to the presence of damage.     

基于扩展有限元法的混凝土细观断裂破坏过程模拟

扩展有限元法（XFEM）是分析不连续力学问题（特别是断裂问题）的一种有效的数值方法。在常规的有限元位移模式中,基于单位分解的思想加入一个跳跃函数和渐进缝尖位移场来对不连续体附近的节点自由度进行局部加强,从而反映了位移的不连续性。介绍了扩展有限元的基本原理,给出了扩展有限元进行混凝土开裂及裂纹扩展的分析方法,最后采用扩展有限元法模拟了湿筛混凝土单轴拉伸作用下及WinklerL-型混凝土板的细观断裂破坏过程。分析了混凝土裂纹萌生、扩展的过程及破坏形态,数值结果与实验结果吻合良好。研究表明：扩展有限元法通过特定的位移模式,使裂纹两侧不连续位移场的表达独立于网格划分,能有效地模拟混凝土材料细观断裂破坏过程。     

On the homogenization of metal matrix composites using strain gradient plasticity

The homogenized response of metal matrix composites（MMC） is studied using strain gradient plasticity.The material model employed is a rate independent formulation of energetic strain gradient plasticity at the micro scale and conventional rate independent plasticity at the macro scale. Free energy inside the micro structure is included due to the elastic strains and plastic strain gradients. A unit cell containing a circular elastic fiber is analyzed under macroscopic simple shear in addition to transverse and longitudinal loading. The analyses are carried out under generalized plane strain condition. Micro-macro homogenization is performed observing the Hill-Mandel energy condition,and overall loading is considered such that the homogenized higher order terms vanish. The results highlight the intrinsic size-effects as well as the effect of fiber volume fraction on the overall response curves, plastic strain distributions and homogenized yield surfaces under different loading conditions. It is concluded that composites with smaller reinforcement size have larger initial yield surfaces and furthermore,they exhibit more kinematic hardening.     

On the transverse compression response of Kevlar KM2 using fiber-level finite element model

Flexible textile composites like woven Kevlar fabrics are widely used in high velocity impact (HVI) applications. Upon HVI they are subjected to both longitudinal tensile and transverse compressive loads. To understand the role of transverse properties, the single fiber and tow transverse compression response (SFTCR and TTCR) of Kevlar KM2 fibers are numerically analyzed using plane strain finite element (FE) models. A finite strain formulation with a minimum number of 84 finite elements is determined to be required for the fiber cross section to capture the finite strain SFTCR through a mesh convergence study. Comparison of converged numerical solution to the experimental results indicates the dominant role of geometric stiffening at finite strains due to growth in contact width. The TTCR is studied using a fiber length scale FE model of a single tow comprised of 400 fibers transversely loaded between rigid platens. This study along with micrographs of yarn after mechanical compaction illustrates fiber spreading and fiber–fiber contact friction interactions are important deformation mechanisms at finite strains. The TTCR is also studied using homogenized yarn level models with properties from the literature. Comparison of TTCR between fiber length scale and homogenized yarn length scale models indicate the need for a nonlinear material model for homogenized approaches to accurately predict the transverse compression response of the fabrics.     

应力波作用下岩石声发射能量特征

通过霍普金森(SHPB)实验系统,对应力波作用下岩石破坏进行声发射实验,获得了该加载条件花岗岩声发射能量的变化规律.实验结果表明,声发射能量变化呈现出两种明显不同的特征:I型,声发射峰值能量之后,能量迅速衰减,到加载的末期,能量出现一定的回升,产生“拐点”;Ⅱ型,声发射峰值能量之后,能量衰减相对Ⅰ型较慢,且不出现“拐点”.通过对岩石破碎块度的分析,得到了声发射峰值能量与岩石破碎分维之间的关系.本文的研究有助于重新寻求应力波下岩石破坏的声发射前兆规律.