BEHAVIOR OF OBSTRUCTED SQUARE BUOYANT VERTICAL JETS IN STATIC AMBIENT （Ⅰ）—— VERIFICATION OF MATHEMATICAL MODEL AND NUMERICAL METHOD

Some experiments were made for the buoyant jet from a square orifice with a square disc placed on it in static ambient and concentration along the axis in self-similar area behind disc was measured. And at the same time a three-dimensional mathematical model was established to simulate the whole flowing under different conditions. All the results predicted by the numerical calculation were substantiated by the experiments. The results were compared with experiential formula for obstructed round buoyant vertical jets in static ambient and it was found that the two concentration distributions had good accordance. Star shape of temperature isolines on cross-sections in the near areas from the disc was found and it was a very special figure for obstructed square buoyant vertical jets with a square disc. The shape will transform to concentric circles gradually alike to the round buoyant vertical jet in self-similar area with increasing of the distance from the disc.     

STRUCTURAL DAMAGE MODEL OF UNSATURATED EXPANSIVE SOIL AND ITS APPLICATION IN MULTI-FIELD COUPLE ANALYSIS ON EXPANSIVE SOIL SLOPE

Focused on the sensitivity to climate change and the special mechanical characteristics of undisturbed expansive soil, an elastc-plastic damage constitutive model was proposed based on the mechanics of unsaturated soil and the mechanics of damage. Undisturbed expansive soil was considered as a compound of non-damaged part and damaged part. The behavior of the non-damaged part was described using non-linear constitutive model of unsaturated soil. The property of the damaged part was described using a damage evolution equation and two yield surfaces, i.e., loading yield (LY) and shear yield (SY). Furthermore, a consolidation model for unsaturated undisturbed expansive soil was established and a FEM program named UESEPDC was designed. Numerical analysis on solid-liquid-gas tri-phases and multi-field couple problem was conducted for four stages and fields of stress, displacement, pore water pressure, pore air pressure, water content, suction, and the damage region as well as plastic region in an expansive soil slope were obtained.     

CFD SIMULATION OF FLUID CATALYTIC CRACKING IN DOWNER REACTORS

A mathematical model has been developed for the simulation of gas-particle flow and fluid catalytic cracking in downer reactors. The model takes into account both cracking reaction and flow behavior through a four-lump reaction kinetics coupled with two-phase turbulent flow. The prediction results show that the relatively large change of gas velocity affects directly the axial distribution of solids velocity and void fraction, which significantly interact with the chemical reaction. Furthermore, model simulations are carried out to determine the effects of such parameters on product yields, as bed diameter, reaction temperature and the ratio of catalyst to oil, which are helpful for optimizing the yields of desired products. The model equations are coded and solved on CFX4.4.     

Continuum modeling for two-lane traffic flow

In this paper, we study the continuum modeling of traffic dynamics for two-lane freeways. A new dynamics model is proposed, which contains the speed gradient-based momentum equations derived from a car-following theory suited to two-lane traffic flow. The conditions for securing the linear stability of the new model are presented. Numerical tests are carried out and some nonequilibrium phenomena are observed, such as small disturbance instability, stop-and-go waves, local clusters and phase transition. The project supported by the National Natural Science Foundation of China (70521001) The English text was polished by Yunming Chen.     

The indenter tip radius effect in micro- and nanoindentation hardness experiments.

Nix and Gao established an important relation between the microindentation hardness and indentation depth. Such a relation has been verified by many microindentation experiments (indentation depths in the micrometer range), but it does not always hold in nanoindentation experiments (indentation depths approaching the nanometer range). Indenter tip radius effect has been proposed by Qu et al. and others as possibly the main factor that causes the deviation from Nix and Gao's relationship. We have developed an indentation model for micro- and nanoindentation, which accounts for two indenter shapes, a sharp, conical indenter and a conical indenter with a spherical tip. The analysis is based on the conventional theory of mechanism-based strain gradient plasticity established from the Taylor dislocation model to account for the effect of geometrically necessary dislocations. The comparison between numerical result and Feng and Nix's experimental data shows that the indenter tip radius effect indeed causes the deviation from Nix-Gao relation, but it seems not be the main factor. The project supported by the National Natural Science Foundation of China (10121202) and the Ministry of Education of China (20020003023)     

基于迭代更新近似模型的车内噪声优化

基于声固耦合分析与迭代更新近似模型对车身板件进行了优化。首先通过板件声学贡献量分析,确定影响车内噪声的关键板件,将其厚度作为设计变量,在减少计算量的同时去除了无关变量的干扰。然后建立不同激励下驾驶员右耳处最大声压与关键板件厚度之间的组合近似模型,对其进行多目标优化并采用序列二次规划和K均值聚类算法对近似模型进行更新,使得采样样本逐渐靠近Pareto前沿,从而提高近似模型在Pareto前沿上的精度。优化结果表明,迭代更新近似模型能够有效优化车身结构,在保证轻量化的同时,降低了车内噪声。     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

磁流体特性对磁流变火炮后坐阻尼器性能的影响

针对火炮后坐磁流变阻尼器的特点,以某单管25 mm火炮实验用磁流变阻尼器为研究对象,基于Herschel-Bulkley本构模型,建立了该炮磁流变后坐阻尼器的轴对称一维层流模型,获得了不同磁场作用下阻尼力随活塞速度的变化规律。运用火炮的后坐运动方程,计算了不同磁流体特性指数下,火炮的后坐位移和后坐速度。计算结果表明,磁流体特性指数的变化对磁流变后坐阻尼器的性能影响显著。     

Hahn-Tsai复合材料的非线性杂交应力有限元方法

成功建立了Hahn-Tsai复合材料模型的非线性杂交应力有限元方程,采用Newton-Raphson迭代法求解结构的非线性位移方程。在迭代过程中,为了提高计算效率可采用简单迭代法由节点位移求解单元应力场。但是,当载荷增加到一定程度以后,非线性应力场由于循环迭代而无法收敛,显然,一般的加速方法不能解决这种循环迭代的发散问题。因此,本文发展了一种确实有效的非线性应力场迭代新方法,在不增加计算工作量的情况下,不仅极大地提高了收敛速度,而且对于较大载荷也能够很好地收敛,从而解决了大载荷下非线性杂交元方法失败的关键问题。数值算例表明该方法是确实可行的。     

混凝土单轴受拉的非局部本构模型

混凝土受拉本构行为存在很强的局部软化现象,使得单轴受拉试验无法给出应力-应变关系,而只能给出应力-位移关系。本文根据内变量理论和等效应变假设建立了基于真实应变的混凝土单轴受力本构方程,并根据Weibull分布可以描述混凝土等脆性材料断裂过程的试验现象,建立了关于弹性应变的损伤演化规律。然后,通过假设平均应变与真实弹性应变的函数关系,在应力-平均应变的本构关系中采用平均弹性应变以描述其非局部行为,而在材料的损伤演化规律中采用真实弹性应变以描述其局部行为,由此建立了单轴受拉荷载条件下的非局部本构模型。最后,对一个单调受拉试验和一个反复受拉试验的仿真结果表明所提出的非局部本构模型可以准确地模拟试验结果。