Crack propagation in structures subjected to periodic excitation

In the present paper, a simple mechanical model is developed to predict the dynamic response of a cracked structure subjected to periodic excitation, which has been used to identify the physical mechanisms in leading the growth or arrest of cracking. The structure under consideration consists of a beam with a crack along the axis, and thus, the crack may open in Mode I and in the axial direction propagate when the beam vibrates. In this paper, the system is modeled as a cantilever beam lying on a partial elastic foundation, where the portion of the beam on the foundation represents the intact portion of the beam. Modal analysis is employed to obtain a closed form solution for the structural response. Crack propagation is studied by allowing the elastic foundation to shorten (mimicking crack growth) if a displacement criterion, based on the material toughness, is met. As the crack propagates, the structural model is updated using the new foundation length and the response continues. From this work, two mechanisms for crack arrest are identified. It is also shown that the crack propagation is strongly influenced by the transient response of the structure.

A moving-coil designed micro-mechanics tester with application on MEMS

A moving-coil designed micro-mechanics tester, named as MicroUTM (universal testing machine), is in-house developed in this paper for micro-mechanics tests. The main component is a moving coil suspended in a uniform magnetic field through a set of springs. When a current passes through the coil, the electromagnetic force is proportional to the magnitude of the current, so the load can easily be measured by the current. The displacement is measured using a capacitive sensor. The load is calibrated using a Sartorius BP211D analytical balance, with a resolution/range of 0.01 mg/80 g or 0.1 mg/210 g. The displacement is calibrated using a HEIDENHAIN CT-6002 length gauge with an accuracy of +/- 0.1 mu m. The calibration results show that the load range is +/- 1 N and the displacement range is +/- 300 mu m. The noise levels of the load and displacement are 50 mu N and 150 nm, respectively. The nonlinearity of the load is only 0.2%. Several in-plane load tests of the MEMS micro-cantilever are performed using this tester. Experimental results, with excellent repeatability, demonstrate the reliability of the load measurement as well as the flexible function of this tester.

Effect of a laser pre-quenched steel substrate surface on the crack driving force in a coating-steel substrate system

A mechanical model of a coating/laser pre-quenched steel substrate specimen with a crack oriented perpendicular to the interface between the coating and the hardened layer is developed to quantify the effects of the residual stress and hardness gradient on the crack driving force in terms of the J-integral. It is assumed that the crack tip is in the middle of the hardened layer of the pre-quenched steel substrate. Using a composite double cantilever beam model, analytical solutions can be derived, and these can be used to quantify the effects of the residual stress and the hardness gradient resulting from the pre-quenched steel substrate surface on the crack driving force. A numerical example is presented to investigate how the residual compressive stress, the coefficient linking microhardness and yield strength and the Young's modulus ratio of the hardened layer to the coating influence the crack driving force for a given crack length. (C) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Deflections and curvatures of a film-substrate structure with the presence of gradient stress in MEMS applications

The close form solutions of deflections and curvatures for a film-substrate composite structure with the presence of gradient stress are derived. With the definition of more precise kinematic assumption, the effect of axial loading due to residual gradient stress is incorporated in the governing equation. The curvature of film-substrate with the presence of gradient stress is shown to be nonuniform when the axial loading is nonzero. When the axial loading is zero, the curvature expressions of some structures derived in this paper recover the previous ones which assume the uniform curvature. Because residual gradient stress results in both moment and axial loading inside the film-substrate composite structure, measuring both the deflection and curvature is proposed as a safe way to uniquely determine the residual stress state inside a film-substrate composite structure with the presence of gradient stress.

The interlaminar fracture-behavior and toughening mechanisms of new carbon-fiber-reinforced bismaleimide composites

The interlaminar fracture behaviour of carbon fibre-reinforced bismaleimide (BMI) composites prepared by using a new modified BMI matrix has been investigated by various methods. Laminates of three typical stacking sequences were evaluated. Double cantilever beam, end-notch flexure and edge-delamination tension tests were conducted under conventional conditions and in a scanning electron microscope. The strain energy release rates in Mode I and Mode III G(lc) and G(llc), as well as the total strain energy release rate, G(mc), have been determined and found to be higher than those for laminates with an epoxy matrix. Dynamic delamination propagation was also studied. The toughening mechanisms are discussed.

不确定边界条件下悬臂梁裂纹参数识别方法研究

基于Bernoulli-Euler梁振动理论,以等效弹簧来模拟裂纹引起的局部软化效应和由非完全固支边界条件引起的转角效应.推导了悬臂梁在不确定边界条件下确定其振动频率的特征方程,直接利用该特征方程,提出一种有效估计裂纹参数的优化方法,通过计算测量频率和理论频率之间的误差目标函数最小化即可识别裂纹参数-裂纹位置和深度.最后,应用两个实例-理想固支边界条件下和非完全固支边界条件下的悬臂梁实验来说明本文方法的有效性.实验结果表明:只需梁结构前三阶频率即可识别裂纹位置和深度.对于理想边界条件下的裂纹参数识别,在测量频率存在小误差情况下,该方法仍能给出比较满意的结果,对于非完全固支边界条件下的裂纹参数识别,利用本文方法能得到比Narkis的方法更精确的裂纹位置识别结果.同时本文方法还能给出比较满意的裂纹深度识别结果.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

An Experimental Study on Flame Propagation in Cornstarch Dust Clouds

Following the quantitative determination of dust cloud parameters, this study investigates the flame propagation through cornstarch dust clouds in a vertical duct of 780 mm height and 160 x 160 mm square cross section, and gives particular attention to the effect of small scale turbulence and small turbulence intensity on flame characteristics. Dust suspensions in air were produced using an improved apparatus ensuring more uniform distribution and repeatable dust concentrations in the testing duct. The dispersion-induced turbulence was measured by means of a particle image velocimetry (PIV) system, and dust concentrations were estimated by direct weighing method. This quantitative assessment made it possible to correlate observed flame behaviors with the parameters of the dust cloud. Upward propagating dust flames, from both closed/open bottom end to open/closed top end of the duct, were visualized by direct light and shadow photography. From the observation of propagation regimes and the measurements of flame velocity, a critical value of the turbulence intensity can be specified below which laminar flame propagation would be established. This transition condition was determined to be 10 cm/s. Laminar flames propagated with oscillations from the closed bottom end to the open top end of the testing duct, while the turbulent flames accelerated continuously. Both laminar and turbulent flames propagated with steady velocity from the open bottom end to the closed top end of the duct. The measured propagation velocity of laminar flames appeared to be in the range of 0.45-0.56 m/s, and it was consistent with the measurements reported in the literature. In the present experimental study, the influence of dust concentration on flame propagation was also examined, and the flame propagation velocity was found weakly sensitive to the variations in dust concentration. Some information on the flame structure was revealed from the shadow records, showing the typical heterogeneous feature of the dust combustion process.

Rayleigh-Marangoni-Bénard Convective Instability

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) and flow patterns in the two-layer system of silicon oil 10cSt and Fluorinert FC70 liquids are studied theoretically and experimentally. Both linear instability analysis and 2D numerical simulation (A=L/H=10) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. Time-dependent oscillations arising at the onset of convection were investigated in a larger various range of two-layer depth ratios (Hr=H1/H2) from 0.2 to 5.0 for different total depth less than 12mm. Our results are different from the previous study on the Rayleig-B閚ard instability and show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. Primary experimental results of the critical instability parameters and the convective structure in the R-M-B convection have been obtained by using the digital particle image velocimetry (DPIV) system, and a good agreement in comparison with the results of numerical simulation was obtained.

Oscillatory Instabilities of Two-Layer Rayleigh-Marangoni-Benard Convection

The Rayleigh-Marangoni-Benard convective instability (R-M-B instability) in the two-layer systems such as Silicone oil (10cSt)/Fluorinert (FC70) and Silicone oil (2cSt)/water liquids are studied. Both linear instability analysis and nonlinear instability analysis (2D numerical simulation) were performed to study the influence of thermocapillary force on the convective instability of the two-layer system. The results show the strong effects of thermocapillary force at the interface on the time-dependent oscillations at the onset of instability convection. The secondary instability phenomenon found in the real two-layer system of Silicone oil over water could explain the difference in the comparison of the Degen's experimental observation with the previous linear stability analysis results of Renardy et al.

Experimental Study on Flow Induced Vibration of a Cylinder with Two Degrees of Freedom Near a Rigid Wall

The flow-induced vibration of a cylinder with two degrees of freedom near a rigid wall under the action of steady flow is investigated experimentally. The vibration amplitude and frequency of the cylinder and the vortex shedding frequency at the wake flow region of the cylinder are measured. The influence of gap-to-diameter ratio upon the amplitude response is analyzed. The experimental results indicate that when the reduced velocity (Vr) is in the range of 1.2 < Vr < 2.6, only streamwise vibration with small amplitude occurs, whose frequency is quite close to its natural frequency in the still water. When the reduced velocity Vr > 3.4, both the streamwise and transverse vibrations of the cylinder occur. In this range, the amplitudes of transverse vibration are much larger than those of streamwise vibrations, and the amplitudes of the streamwise vibration also get larger than those at the range of 1.2 < Vr < 2.6. At the range of Vr > 3.4, the frequency of streamwise vibration undergoes a jump at certain values of Vr, at which the streamwise vibrating frequency is twice as much as the transverse one. However, when the streamwise vibration does not experience a jump, its frequency is the same as that of the transverse vibration. The maximum values of second streamwise and transverse vibration amplitudes increase with increasing gap-to-diameter ratios.

Experimental investigation of Very-High-Cycle Fatigue behavior for a structural steel

Very-High-Cycle Fatigue (VHCF) test for a medium carbon structural steel (40Cr) has been performed and a stepwise S-N curve was obtained by employing cantilever-type rotary bending fatigue machine with hourglass shape specimen. The S-N curve was well explained as a combination of curves for surface-induced fracture and interior inclusion-induced fracture with fish-eye patterns. The morphology of the fish-eye pattern was illustrated in order to clarify subsurface crack initiation and propagation behavior.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.