Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

Numerical Analysis of Theoretical Model of the Rf MEMs Switches

An improved electromechanical model of the RF MEMS (radio frequency microelectromechanical systems) switches is introduced, in which the effects of intrinsic residual stress from fabrication processes, axial stress due to stretching of beam, and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model. A semi-analytical method is developed to calculate the behavior of the RF MEMS switches. Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared, and the corresponding analysis with the dimensionless numbers is conducted too. The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

A Numerical Model of Rapid Solidification Processing of Ni-Al Alloy in Planar Flow Casting

A numerical model has been developed for simulating the rapid solidification processing (RSP) of Ni-Al alloy in order to predict the resultant phase composition semi-quantitatively during RSP. The present model couples the initial nucleation temperature evaluating method based on the time dependent nucleation theory, and solidified volume fraction calculation model based on the kinetics model of dendrite growth in undercooled melt. This model has been applied to predict the cooling curve and the volume fraction of solidified phases of Ni-Al alloy in planar flow casting. The numerical results agree with the experimental results semi-quantitatively.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

Electromechanical Model of RF MEMS Switches

With the recent rapid growth of Radio Frequency Micro-Electro-Mechanical Systems (RF MEMS) switches, there has developed an emergent requirement for more accurate theoretical models to predict their electromechanical behaviors. Many parameters exist in the analysis of the behavior of the switch, and it is inconvenient for further study. In this paper, an improved model is introduced, considering simultaneously axial stress, residual stress, and fringing-field effect of the fixed-fixed bridge structure. To avoid any unnecessary repetitive model tests and numerical simulation for RF MEMS switches, some dimensionless numbers are derived by making governing equation dimensionless. The electromechanical behavior of the fixed-fixed bridge structure of RF MEMS switches is totally determined by these dimensionless numbers.

Mathematical Modeling of Near-Wall Flows of Two-Phase Mixture with Evaporating Droplets

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream. 在双连续介质理论框架下,采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程,给出了控制汽雾流的相似判据。对于沿曲面的流动,边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体验点附近蒸汽。液滴边界层的数值计算结果。它们表明：由于蒸发,在边界层内近壁处形成了一个无液滴区域；在该区上边界处,液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加,甚至在自由来流中液滴质量浓度很低时此效应依然存在。

Unsteady Model of Drop Marangoni Migration in Microgravity

The experiments of drop Marangoni migration have been performed by the drop shift facility of short period of 4.5 s, and the drop accelerates gradually to an asymptotic velocity during the free fall. The unsteady and axisymmetric model is developed to study the drop migration for the case of moderate Reynolds number Re = O(1), and the results are compared with the experimental ones in the present paper. Both numerical and experimental results show that the migration velocity for moderate Reynolds number is several times smaller than that given by the linear YGB theory.

A simulation model of a floating half zone

A simulation model of a floating half zone was suggested by steady numerical simulation and experiment respectively, in the previous papers [Q.S. Chen, W.R. Hu, Int. J. Mass Heat Transfer 40 (1997) 757; J.H. Han, Y. Ar, R. Zhou, W.R. Hue, Int. J. Mass Heat Transfer 40 (1997) 2671]. In the present paper, the simulation model is studied by using the method of unsteady and three-dimensional numerical simulation, and the transient process from steady convection to oscillatory convection is especially analyzed. Comparison of onsets of oscillation for both simulation model and the usual model were obtained, and the results show that the critical Marangoni number of the simulation model is obviously smaller than that of the usual model for the same slender liquid bridge. This implies that the usual model of a floating half zone gives a lower estimation on the onset of oscillation for floating zone convection.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

A new model of thermocapillary convection in floating half zone

A new model of thermocapillary convection in floating half zone was suggested in the present paper. The liquid bridge floats between two co-axis rods, the lower rod consists of metal with constant temperature and the upper rod consists of thermal insulating materials, where the normal gradient of temperature is nearly zero. In this case, the new model is relatively closer to simulate a half part of floating full zone in microgravity environment in comparison with the usual model of floating half zone. Basic features of the new model were studied by both numerical simulation and experiments, and the comparisons with the usual model were also discussed.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.