I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Theoretical investigation of turbulent boundary layer over a wavy surface

The important features of the two-dimensional incompressible turbulent flow over a wavy surface of wavelength comparable with the boundary layer thickness are analyzed.

A turbulent field method using model equation for turbulent shear stress similar to the scheme of Bradshaw, Ferriss and Atwell (1967) is employed with suitable modification to cover the viscous sublayer. The governing differential equations are linearized based on the small but finite amplitude to wavelength ratio. An orthogonal wavy coordinate system, accurate to the second order in the amplitude ratio, is adopted to avoid the severe restriction to the validity of linearization due to the large mean velocity gradient near the wall. Analytic solution up to the second order is obtained by using the method of matched-asymptotic-expansion based on the large Reynolds number and hence the small skin friction coefficient.

In the outer part of the layer, the perturbed flow is practically "inviscid." Solutions for the velocity, Reynolds stress and also the wall pressure distributions agree well with the experimental measurement. In the wall region where the perturbed Reynolds stress plays an important role in the process of momentum transport, only a qualitative agreement is obtained. The results also show that the nonlinear second-order effect is negligible for amplitude ratio of 0.03. The discrepancies in the detailed structure of the velocity, shear stress, and skin friction distributions near the wall suggest modifications to the model are required to describe the present problem.

The effect of surface hardening on the elastic shakedown of elliptical contacts

The effect of varying both the aspect ratio and the coefficient of friction of contacts with elliptical geometry on their elastic shakedown performance has been examined theoretically for surfaces with two types of subsurface hardness or strength profiles. In stepwise hardening the hard layer is of uniform strength while in linear hardening its strength reduces from a maximum at the surface to that of the core at the base of the hardened layer. The shakedown load is expressed as the ratio of the maximum Hertzian pressure to the strength of the core material. As the depth of hardening, expressed as a multiple of the elliptical semi-axis, is increased so the potential shakedown load increases from a level that is appropriate to a uniform half-space of unhardened material to a value reflecting the hardness of the surface and near-surface material. In a step-hardened material, the shakedown limit for a surface 'pummelled' by the passage of a sequence of such loads reaches a cut-off or plateau value, which cannot be exceeded by further increases in hardening depth irrespective of the value of the friction coefficient. For a linear-hardened material the corresponding plateau is approached asymptotically. The work confirms earlier results on the upper bounds on shakedown of both point and line contacts and provides numerical values of shakedown loads for intermediate geometries. In general, the case depth required to achieve a given shakedown limit reduces in moving from a transversely moving nominal line load to an axisymmetric point load.

Rolling of thin strip and foil: Application of a tribological model for "mixed" lubrication

A mechanical model of cold rolling of foil is coupled with a sophisticated tribological model. The tribological model treats the "mixed" lubrication regime of practical interest, in which there is "real" contact between the roll and strip as well as pressurized oil between the surfaces. The variation of oil film thickness and contact ratio in the bite is found by considering flattening of asperities on the foil and the build-up of hydrodynamic pressure through the bite. The boundary friction coefficient for the contact areas is taken from strip drawing tests under similar tribological conditions. Theoretical results confirm that roll load and forward slip decrease with increasing rolling speed due to the decrease in contact ratio and friction. The predictions of the model are verified using mill trials under industrial conditions. For both thin strip and foil, the load predicted by the model has reasonable agreement with the measurements. For rolling of foil, forward slip is overestimated. This is greatly improved if a variation of friction through the bite is considered.

Electrotactile touch surface by using transparent graphene

In this work we present a flexible Electrostatic Tactile (ET) surface/display realized by using new emerging material graphene. The graphene is transparent conductor which successfully replaces previous solution based on indium-thin oxide (ITO) and delivers more reliable solution for flexible and bendable displays. The electrostatic tactile surface is capable of delivering programmable, location specific tactile textures. The ET device has an area of 25 cm 2, and consists of 130 μm thin optically transparent (>76%) and mechanically flexible structure overlaid unobtrusively on top of a display. The ET system exploits electro vibration phenomena to enable on-demand control of the frictional force between the user's fingertip and the device surface. The ET device is integrated through a controller on a mobile display platform to generate fully programmable range of stimulating signals. The ET haptic feedback is formed in accordance with the visual information displayed underneath, with the magnitude and pattern of the frictional force correlated with both the images and the coordinates of the actual touch in real time forming virtual textures on the display surface (haptic virtual silhouette). To quantify rate of change in friction force we performed a dynamic friction coefficient measurement with a system involving an artificial finger mimicking the actual touch. During operation, the dynamic friction between the ET surface and an artificial finger stimulation increases by 26% when the load is 0.8 N and by 24% when the load is 1 N. © 2012 ACM.

Flow physics of a normal-hole bled supersonic turbulent boundary layer

An experimental study on normal hole bleed in a supersonic turbulent boundary layer has been conducted. A combination of LDV, Schlieren imagery and oil flow visualization were used to provide a better understanding of the three-dimensional flow field surrounding a supersonic bleed array. Experiments were performed at Mach numbers of 1.8 and 2.5, while previously published results at Mach numbers of 1.3 and 1.5 were also incorporated. The bleed system was capable of removing up to approximately 10% of the incoming boundary layer through a tunnel-spanning array of discrete holes with diameters the same order of magnitude of boundary layer displacement thickness. Inspection of boundary layer profiles downstream of the bleed region indicates that vorticity generated by the discrete holes can have a substantial influence on changes to the boundary layer shape factor and skin friction coefficient, through modification of the lower 20% of the boundary layer. This vorticity was visualized through oil-flow visualization, and LDV measurements, showing the development of two vortices off each bleed hole, and corresponding upwash and downwash regions with far-reaching three dimensional effects. © 2013 by J. M. Oorebeek and H. Babinsky.

Effect of counterpart on the tribological behavior and tribo-induced phase transformation of Si

The tribological behaviors and phase transformation of single crystal silicon against Si3N4, Ruby and steel were investigated in this study. It was found that the strong chemical action between silicon and Fe was the key factor to the tribological behavior of silicon as slid against steel. SEM and Raman spectroscopy indicated that phase transformation of single crystal silicon occurred during the running-in period at low sliding velocity as slid against Si3N4 and Ruby. and gave birth to single or a mixture phase of Si-III, Si-XII and amorphous silicon. The high hardness of counterpart and the absence of chemical action between silicon and counterpart facilitated the phase transformation of single crystal silicon. (C) 2008 Elsevier Ltd. All rights reserved.

Study on mechanical properties of GaN epitaxy films grown on sapphire by MOCVD

GaN epitaxy films were grown on (0001) oriented sapphire substrate by metal-organic vapor deposition(MOCVD). AFM and SEM were used to analyze the surface morphology of GaN films. Hardness and critical load of GaN films were measured by an nano-indentation tester, friction coefficient by reciprocating UMT-2MT tribometer. It is found that the surface of GaN film is smooth and the epitaxial growth mechanism is in two-dimension mode, GaN epitaxy films also belong to ultra-hardness materials, whose hardness is 22.1 MPa and elastic modulus is 299.5 GPa. Adhesion strength of epitaxial GaN to sapphire is high, and critical load reaches 1.6 N. Friction coefficient against GCr15 ball is steadily close to 0.13, while GaN films turns to be broken rapidly by using Si3N4 ceramic ball as counterpart.

潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

大型铰接并联六维测力平台摩擦建模

在已有制造工艺及标定技术基础上,为进一步改善大型铰接并联六维测力平台的测量精度,本文基于螺旋理论和影响系数原理,引入符号函数建立了Stewart结构大型铰接六维测力平台的摩擦模型。文中提出了关节摩擦对铰接并联六维测力平台测量精度的影响矩阵及I、H类误差表达式,绘制了在不同外载和关节摩擦系数条件下,六维测力平台的I、II类误差曲线,并总结丁关节摩擦和平台自重对测力平台测量精度的影响规律。为具有普通球形铰链人型Stewart平台六维测力下台精度的提高和改善提供了理论基础。

蛇形机器人行波运动的研究

对蛇形机器人的行波运动机理进行了研究,基于Serpenoid曲线建立了蛇形机器人行波运动的形状控制模型、运动学模型、机构动力学模型及环境动力学模型,给出了求解过程。仿真结果表明:在行波运动中,各关节输入力矩大小呈周期变化,且随其与蛇体重心距离的增加而减小。蛇体中心关节的最大输入力矩为蛇形机器人所需最大输入力矩。对称关节输入力矩幅值变化规律相同,但相差不同;运动初始弯角保持不变,关节输入力矩随着环境摩擦因数的增加而增加。环境摩擦因数保持不变,关节输入力矩随着初始弯角的增加而减小。

Fully developed laminar flow and heat transfer in smooth trapezoidal microchannel

Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal micro-channels has been studied with uniform wall heat flux boundary conditions. Through coordinate transformation, the governing equations are transformed from physical plane to computational domain, and the resulting equations are solved by a finite-difference scheme. The influences of velocity slip and temperature jump on friction coefficient and Nusselt number are investigated in detail. The calculation also shows that the aspect ratio and base angle have significant effect on flow and heat transfer in trapezoidal micro-channel. (c) 2005 Elsevier Ltd. All rights reserved.

Fission dynamics for capture reactions in 58,64ni + 208pb systems: New results in terms of thermal energy and neutron multiplicity correlated distributions

The neutron multidetector DéMoN has been used to investigate the symmetric splitting dynamics in the reactions 58.64Ni + 208Pb with excitation energies ranging from 65 to 186 MeV for the composite system. An analysis based on the new backtracing technique has been applied on the neutron data to determine the two-dimensional correlations between the parent composite system initial thermal energy (EthCN) and the total neutron multiplicity (νtot), and between pre- and post-scission neutron multiplicities (νpre and νpost, respectively). The νpre distribution shape indicates the possible coexistence of fast-fission and fusion-fission for the system 58Ni + 208Pb (Ebeam = 8.86 A MeV). The analysis of the neutron multiplicities in the framework of the combined dynamical statistical model (CDSM) gives a reduced friction coefficient β = 23 ± 2512 × 1021 s-1, above the one-body dissipation limit. The corresponding fission time is τf = 40 ± 4620 × 10-21 s. © 1999 Elsevier Science B.V. All rights reserved.

An isothermal-isobaric Langevin thermostat for simulating nanoparticles under pressure: Application to Au clusters

We present a method for simulating clusters or, molecules subjected to an external pressure, which is exerted by a pressure-transmitting medium. It is based on the canoninical Langevin thermostat, but extended in such a way that the Brownian forces are allowed to operate only from the region exterior to the cluster. We show that the frictional force of the Langevin thermostat is linked to the pressure of the reservoir in a unique way, and that this property manifests itself when the particle it acts upon is not pointlike but has finite dimensions. By choosing appropriately the strength of the random forces and the friction coefficient, both temperature and pressure can be controlled independently. We illustrate the capabilities of this new method by calculating the compressibility of small gold clusters under pressure.

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.