Characteristics of elastic waves in FGM spherical shells,an analytical solution

In this paper, we investigate the propagation characteristics of elastic guided waves in FGM spherical shells with exponentially graded material in the radial direction. A new separation of variables technique to displacements is proposed to convert the governing equations of the wave motion to the second-order ordinary differential equations with variable coefficients. By further a variables transform technique, these equations are transformed to the Whittaker equations so that analytic solutions can be obtained. For the spherical shell case, by satisfying the traction-free boundary conditions on both the inner and outer surfaces of the shell, we obtain the dispersion equations, which show that both the SH and Lamb-type wave modes are generated in the structure. The calculated dispersion curves in the functionally graded shell demonstrate a clear influence of the gradient coefficient as compared to those of the homogeneous shell, with the Lamb-type waves more sensitive to the gradient coefficient. The mode shapes and the distributions of stresses in the shells for various gradient coefficients are also presented to illustrate their dependence on the gradient coefficient.     

One-dimensional analytical model for the response of elastic coatings to water blast

In previous studies, the response of sandwich structures to water blast is solved by envisaging the plate adjacent to water as a rigid one, while the effects of the elasticity in the longitudinal direction of the plate are rarely studied. In this paper, a monolithic elastic coating with varying stiffness and thickness is investigated by a one-dimensional analytical approach, based on linear wave motion theory, to reveal the elastic effect of the plate on the incident wave. One side of the coating is loaded by a planar shock wave; on the opposite side, rigid boundary or air-backed boundary is imposed. The fluid–structure interaction (FSI), cavitation phenomenon and large deformation of the coating are taken into account. In particular, the initiation and evolution of cavitation, including the propagations of breaking fronts and closing fronts, as well as the pressure histories of radiated wave by the closing front, are examined. The analytical solution has been compared with finite element (FE) predictions. The results are found to be in excellent agreement for the propagation of breaking front and closing front, as well as the pressure and particle velocity histories at the wet face before the cavitation reaches the wet face. However, when the wet face cavitates, the predictions provided by the analytical method are less accurate and the analytically-computed particle velocity can only be compared in an average sense with the FE predictions. For air-backed case, Taylor׳s model prior to cavitation becomes a trivial case of the analytical model and the comparison also indicates the validity of the analytical model. The validated analytical model is used to determine the dependence of the peak pressure at the wet face and the impulse transmitted to the coating on the coating properties, including the wave impedance and thickness.     

THE ANALYTICAL STUDY ON THE LASER INDUCED REVERSE-PLUGGINGEFFECT BY USING THE CLASSICAL ELASTIC PLATETHEORY (Ⅰ)-TEMPERATURE FIELDS

The temperature distributions in the metallic foils induced by spatially cylindrical long-pulsed laser is exammed in order to analyse the newly-discovered reverse-plugging effect (RPE). An exact solution for the temperature fields is derived by using the Hankel transform and Laplace transform. Numerical results are obtained for both spatial distributions with Gaussian and cylindrical types. The results show that the spatially cylindrical distribution of laser offers a formidable potential for the RPE. Project supported by the Field of Laser Technology of the National High-Tech. Research and Development Programme, i.e., the 863 plan.     

A model for piecewise linear stress near the concrete crack tip and the solution

A model for the piecewise linear stress near the tip of a crack in concrete is proposed. Based on this model and FCM, the analytical solutions of the stress and displacement fields for mode I crack in concrete are obtained. The functional relationship between the dimensions of fracture process zone and the initial crack length is established.     

Line-integral representations for the elastic displacements,stresses and interaction energy of arbitrary dislocation loops in transversely isotropic bimaterials

The elastic displacements, stresses and interaction energy of arbitrarily shaped dislocation loops with general Burgers vectors in transversely isotropic bimaterials (i.e. joined half-spaces) are expressed in terms of simple line integrals for the first time. These expressions are very similar to their isotropic full-space counterparts in the literature and can be easily incorporated into three-dimensional (3D) dislocation dynamics (DD) simulations for hexagonal crystals with interfaces/surfaces. All possible degenerate cases, e.g. isotropic bimaterials and isotropic half-space, are considered in detail. The singularities intrinsic to the classical continuum theory of dislocations are removed by spreading the Burgers vector anisotropically around every point on the dislocation line according to three particular spreading functions. This non-singular treatment guarantees the equivalence among different versions of the energy formulae and their consistency with the stress formula presented in this paper. Several numerical examples are provided as verification of the derived dislocation solutions, which further show significant influence of material anisotropy and bimaterial interface on the elastic fields and interaction energy of dislocation loops.     

Dynamic behaviour of a bubble near an elastic infinite interface

This paper presents a study to describe the behaviour of a non-equilibrium bubble in a fluid (Fluid 1) that is in contact with another fluid (Fluid 2). Fluid 2 is assumed to incorporate some elastic properties, which are modelled through a pressure term at the fluid–fluid interface. The Laplace equation is assumed to be valid in both fluids and the boundary integral method is employed to simulate the dynamics of the bubble and the fluid–fluid interface. Interesting characteristic phenomena concerning bubble oscillations and the deformation of the fluid–fluid interface are studied for a range of parameters (distance from the fluid–fluid interface, density ratios of the two fluids and elastic properties of Fluid 2). Some of the phenomena observed are jet formation in the bubble, bubble splitting, a ring bubble separating from the main bubble, mushroom-shaped bubbles and the dynamic elevation of the elastic interface. Most of these phenomena are only observed when Fluid 2 possesses some elastic properties (besides the usual formation of a high speed liquid jet). Comparisons with experimental observations confirm the validity of our simulations.     

V型缺口裂端的三维应力状态及约束分析

利用三维有限元方法,研究了有限厚度板中V型缺口根部穿透裂纹前沿的三维弹性应力场。对不同厚度、不同缺口张开角和裂纹长度对应力强度因子及裂尖附近的三维约束程度的影响进行了分析,同时还讨论了三维约束区的大小。研究结果显示：当缺口张开角小于60度时,不同缺口的应力强度因子和离面约束因子的分布基本一致,角度的影响不明显。应力强度因子是厚度的函数,板中面的应力强度因子随厚度的增加逐渐减小趋近干平面值,最大为1．08倍的平面值。当板厚超过15倍的缺口深度时,应力强度因子最大值将从中面转移至接近自由表面位置,距中面约0．4倍板厚。三维约束非常明显的区域在裂尖前沿约0.45倍厚度的范围内．二维到三维的过渡区在裂尖前沿1.5倍厚度的区域内;在中面上三维效应影响区最大,随着离中面距离的增加逐渐减小,在自由表面上降为0。     

A general analytical solution for elastic vibration of rectangular thin plates

A general solution of differential equation for lateral displacement lunction of rectangular elastic thin plates in free vibration is established in this paper. It can be used to solve the vibration problem of rectangular plate with arbitrary boundaries. As an example, the frequency and its vibration mode of a rectangular plate with four edges free have been solved.     

Some possible alternative solutions of near-tip fields for steady crack growth in elastic perfectly-plastic media

Some possible alternative solutions of near-tip fields are studied for plane-strain Mode—I quasi-static steady crack growth in incompressible (ν=1/2) elastic perfectly-plastic media. A group of four-sector solutions and a three-sector solution in which the elastic-unloading region vanishes are given. Stress functions, plastic flow factors and plastic strains in each region are also given. Project supported by the State Education Commission under a funding program for Excellent University Young Faculties and National Natural Science Foundation of China.     

Interfacial excess energy, excess stress and excess strain in elastic solids: Planar interfaces

In this paper, interfacial excess energy and interfacial excess stress for coherent interfaces in an elastic solid are reformulated within the framework of continuum mechanics. It is shown that the well-known Shuttleworth relationship between the interfacial excess energy and interfacial excess stress is valid only when the interface is free of transverse stresses. To account for the transverse stress, a new relationship is derived between the interfacial excess energy and interfacial excess stress. Dually, the concept of transverse interfacial excess strain is also introduced, and the complementary Shuttleworth equation is derived that relates the interfacial excess energy to the newly introduced transverse interfacial excess strain. This new formulation of interfacial excess stress and excess strain naturally leads to the definition of an in-plane interfacial stiffness tensor, a transverse interfacial compliance tensor and a coupling tensor that accounts for the Poisson's effect of the interface. These tensors fully describe the elastic behavior of a coherent interface upon deformation.     

An analytical solution for underground structure-country rock dynamic interaction

In this paper according to the results of a great quantity of tests and numerical calculations, it is pointed out that the surrounding rock with thickness of 1/3 span has mechanical characteristics of a thick flexural member in underground rock cave subjected to transverse blast loading. However, it approaches the stress state of free field outside country rock and the equation of thick plate theory under the loading of free field pressure may be applied to the solution of this problem. Therefore, the underground structure-country rock dynamic interaction may be described by dynamic equations of flexural members of thin plate and thick plate, which express linear and country rock respectively. The interaction force between the liner and country rock is expressed as contacting pressure function q(x, l). By solving system of simultaneous equations, the analytical solution for the dynamic analysis of arch-straight liner including elastic half space interaction effect is given and the analytical expression of function q(x, t) is obtained. This analytic solution will contribute to the study of some substantive problems of underground structure-medium dynamic interaction theoretically. (Instutute of Engineering Mechanics, Academia Sinica)     

An analytical expression of reliability solution for Druker-Prager criterion

Strength parameters of rock mass are taken as random variables. Based on the reliability theory, an analytical expression of reliability solution for Drucker-Prager criterion is given. The Monte-Carlo method and FORM （first order reliability method） are used to test the correctness of the solution. Some influences of variation of strength parameters on the judging results of Drucker-Prager criterion are discussed by the use of the solution. Conclusions of discussion show that variations of strength parameters have different influences on the reliability probability of Drucker-Prager criterion. When the coefficients of variation of strength parameters are great, their influences on the reliability probability of Drucker-Prager criterion are obvious and can not be neglected. The conclusions of this paper provide a new way for the reliability judgment of yield on rock m＇ass.     

The lowest order solution of the crack–inhomogeneity interaction for mode I crack

The lowest order solution to crack–inhomogeneity interaction is derived for mode I crack. The basic solution evaluates the variation of near-tip stress intensity factors induced by an inhomogeneity of arbitrary shape. A set of simplified forms of the basic solution is also obtained for several special inhomogeneity shapes. As validated by numerical examples, the approximate solutions have good accuracy.     

Acoustic radiation force on an elastic cylinder in a Gaussian beam near an impedance boundary

Acoustic radiation force (ARF) is studied by considering an infinite elastic cylinder near an impedance boundary when the cylinder is illuminated by a Gaussian beam. The surrounding fluid is an ideal fluid. Using the method of images and the translation-addition theorem for the cylindrical Bessel function, the resulting sound field including the incident wave, its reflection from the boundary, the scattered wave from the elastic cylinder, and its image are expressed in terms of the cylindrical wave function. Then, we deduce the exact equations of the axial and transverse ARFs. The solutions depend on the cylinder position, cylinder material, beam waist, reflection coefficient, distance from the impedance boundary, and absorption in the cylinder. To analyze the effects of the various factors intuitively, we simulate the radiation force for non-absorbing elastic cylinders made of stainless steel, gold, and beryllium as well as for an absorbing elastic cylinder made of polyethylene, which is a well-known biomedical polymer. The results show that the impedance boundary, cylinder material, absorption in the cylinder, and cylinder position in the Gaussian beam significantly affect the magnitude and direction of the force. Both stable and unstable equilibrium regions are found. Moreover, a larger beam waist broadens the beam domain, corresponding to non-zero axial and transverse ARFs. More importantly, negative ARFs are produced depending on the choice of the various factors. These results are particularly important for designing acoustic manipulation devices operating with Gaussian beams.     

A constrained theory of a Cosserat point for the numerical solution of dynamic problems of non-linear elastic rods with rigid cross-sections

A constrained theory of a Cosserat point has been developed for the numerical solution of non-linear elastic rods. The cross-sections of the rod element are constrained to remain rigid but tangential shear deformations and axial extension are admitted. As opposed to the more general theory with deformable cross-sections, the kinetic coupling equations in the numerical formulation of the constrained theory are expressed in terms of the simple physical quantities of force and mechanical moment applied to the common ends of neighboring elements. Also, in contrast with standard finite element methods, the Cosserat element uses a direct approach to the development of constitutive equations. Specifically the kinetic quantities are determined by algebraic expressions which are obtained by derivatives of a strain energy function. Most importantly, no integration is needed over the element region. A number of example problems have been considered which indicate that the constrained Cosserat element can be used to model large deformation dynamic response of non-linear elastic rods.     

An analytical solution for VOCs sorption on dry building materials

Based on the diffusion within the material and the mass transfer through the air boundary layer, a VOCs sorption model with an arbitrary inlet concentration is proposed. With the use of Laplace transform and convolution theorem, a method to get the analytical solution has been developed, which helps to understand the mechanism of the sorption theoretically and makes the computation convenient. Experimental results on VOCs sorption in an environmental chamber are used to validate the present model and a good agreement has been obtained. Two cases of different inlet concentrations are discussed in detail. The patterns of the concentration distribution in the air vary with the inlet concentration. The influences of model parameters have been investigated. With the increase of the diffusion coefficient and the partition coefficient, the rate of the adsorption increases. Relatively, the partition coefficient plays a more important role in the sorption than the diffusion coefficient.     

An analytical solution for the population balance equation using a moment method

Brownian coagulation is the most important inter-particle mechanism affecting the size distribution of aerosols. Analytical solutions to the governing population balance equation (PBE) remain a challenging issue. In this work, we develop an analytical model to solve the PBE under Brownian coagulation based on the Taylor-expansion method of moments. The proposed model has a clear advantage over conventional asymptotic models in both precision and efficiency. We first analyze the geometric standard deviation (GSD) of aerosol size distribution. The new model is then implemented to determine two analytic solutions, one with a varying GSD and the other with a constant GSD. The varying solution traces the evolution of the size distribution, whereas the constant case admits a decoupled solution for the zero and second moments. Both solutions are confirmed to have the same precision as the highly reliable numerical model, implemented by the fourth-order Runge–Kutta algorithm, and the analytic model requires significantly less computational time than the numerical approach. Our results suggest that the proposed model has great potential to replace the existing numerical model, and is thus recommended for the study of physical aerosol characteristics, especially for rapid predictions of haze formation and evolution.     

On the effect of interactions of inhomogeneities on the overall elastic and conductive properties

A simple method of estimating the effect of inhomogeneity interactions on the overall properties (elastic and conductive) is developed. It is formulated in terms of property contribution tensors that give the contribution of an inhomogeneity to the overall properties. The method can be viewed as further development of the approach of Rodin and Hwang (1991) and Rodin (1993) that generalized the method of analysis of crack interactions (Kachanov, 1987) to inhomogeneities. We also extend the method to the conductive properties. Considering the effect of interactions on the property contribution tensors on the example of pores we find that this effect is generally moderate, at most (even when pores touch one another) – in contrast with the effect on local fields. On example of two spheres, we compare the interaction effects on the elastic and the conductive properties, and discuss the impact of interactions on the cross-property connections.     

An analytical solution for the population balance equation using a moment method

Brownian coagulation is the most important inter-particle mechanism affecting the size distribution of aerosols.Analytical solutions to the governing population balance equation(PBE) remain a challenging issue.In this work,we develop an analytical model to solve the PBE under Brownian coagulation based on the Taylor-expansion method of moments.The proposed model has a clear advantage over conventional asymptotic models in both precision and efficiency.We first analyze the geometric standard deviation(GSD) of aerosol size distribution.The new model is then implemented to determine two analytic solutions,one with a varying GSD and the other with a constant GSD.The varying solution traces the evolution of the size distribution,whereas the constant case admits a decoupled solution for the zero and second moments.Both solutions are confirmed to have the same precision as the highly reliable numerical model,implemented by the fourth-order Runge-Kutta algorithm,and the analytic model requires significantly less computational time than the numerical approach.Our results suggest that the proposed model has great potential to replace the existing numerical model,and is thus recommended for the study of physical aerosol characteristics,especially for rapid predictions of haze formation and evolution.