The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Waverider configurations derived from general conical flowfields

A method based on the computational fluid dynamics (CFD) is presented for a flexible waverider's design. The generating bodies of this method could be any cones. In addition, either the leading edge or the profile of the scramjet's inlet is used as the waverider's definition curve, parameterized by the quadric function, the sigmoid function or the B-spline function. Furthermore, several numerical examples are carried out to validate the method and the relevant codes. The CFD results of the configurations show that all the designs are successful. Moreover, primary suggestions are proposed for practical design by comparing the geometrical and aerodynamic performances of the cone-derived waveriders at Mach 6.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

Electric characteristics of pseudo-spark during breakdown phase

An analysis of the time-dependent resistive voltage and power deposition during the breakdown phase of pseudo-spark is presented. The voltage and current were measured by specially designed low-inductance capacitive voltage divider and current measuring resistor. The measured waveforms of voltage and current are digitized and processed by a computer program to remove the inductive component, so as to obtain resistive voltage and power deposition. The influence of pressure, cathode geometry and charging voltage of storage capacitors on the electrical properties in the breakdown phase are investigated. The results suggest that the breakdown phase of pseudo-spark consists of three stages. The first stage is mainly hollow cathode discharge. In the second stage, field-enhanced thermionic emission takes place, resulting in a fast voltage drop and sharp rise of discharge current. The third stage of discharge depends simply on the parameters of the discharge circuit.

A general theory on media with randomly distributed inclusions .1. The average field behaviors

In this paper, a theory is developed to calculate the average strain field in the materials with randomly distributed inclusions. Many previous researches investigating the average field behaviors were based upon Mori and Tanaka's idea. Since they were restricted to studying those materials with uniform distributions of inclusions they did not need detailed statistical information of random microstructures, and could use the volume average to replace the ensemble average. To study more general materials with randomly distributed inclusions, the number density function is introduced in formulating the average field equation in this research. Both uniform and nonuniform distributions of inclusions are taken into account in detail.

A ground hydrologic model with inclusion of a layer of vegetation canopy that can interact with general-circulation model

In this paper, a ground hydrologic model(GHM) is presented in which the vapor, heat and momentum exchanges between ground surface covers (including vegetation canopy) and atmosphere is described more realistically. The model is used to simulate three sets of field data and results from the numerical simulation agree with the field data well. GHM has been tested using input data generated by general circulation model (GCM) runs for both the North American regions and the Chinese regions, The results from GHM are quite different from those of GHMs in GCMs. It shows that a more active concerted effort on the land surface process study to provide a physically realistic GHM for predicting the exchange between land and atmosphere is important and necessary.

Explorations of electric-current system in solar active regions .1. Empirical inferences of the current flows

In this paper we explore techniques to identify sources of electric current systems and their channels of flow in solar active regions. Measured photospheric vector magnetic fields (VMF) together with high-resolution white-light and H filtergrams provide the data base to derive the current systems in the photosphere and chromosphere. Simple mathematical constructions of fields and currents are also adopted to understand these data. As an example, the techniques are then applied to infer current systems in AR 2372 in early April 1980. The main results are: (i) In unipolar sunspots the current density may reach values of 103 CGSE, and the Lorentz force on it can accelerate the Evershed flow, (ii) Spots exhibiting significant spiral pattrn in the penumbral filaments are the sources of vertical major currents at the photospheric surface, (iii) Magnetic neutral lines where the transverse field was strongly sheared were channels along which strong current system flows, (iv) The inferred current systems produced oppositely-flowing currents in the area of the delta configuration that was the site of flaring in AR 2372.

Evolution and propagation of density waves on galactic disk .1. General evolution equation of wavelet and discussion of its solution

This paper presents a general self-consistent theory of evolution and propagation of wavelets on the galactic disk. A simplified model for this theory, i. e. the thin transition-layer approximation is proposed.There are three types of solutions to the basic equation governing the evolution of wavelets on the disk: (ⅰ) normal propagating type; (ⅱ) swing type; (ⅲ) general evolving type. The results show that the first two types are applicable to a certain domain on the galactic disk and a certain region of the wave number of wavelets. The third is needed to join the other two types and to yield a coherent total picture of the wave motion. From the present theory, it can be seen that the well-known "swing theory" of the G-L sheet model holds only for a certain class of basic states of galaxies.

The influence of fluctuation fields on the force-free field .2. The local expansion

In Paper I (Hu, 1982), we discussed the the influence of fluctuation fields on the force-free field for the case of conventional turbulence and demonstrated the general relationships. In the present paper, by using the approach of local expansion, the equation of average force-free field is obtained as (1+b)×B 0=(#x002B;a)B 0#x002B;a (1)×B 0#x002B;K. The average coefficientsa,a (1),b, andK show the influence of the fluctuation fields in small scale on the configurations of magnetic field in large scale. As the average magnetic field is no longer parallel to the average electric current, the average configurations of force-free fields are more general and complex than the usual ones. From the view point of physics, the energy and momentum of the turbulent structures should have influence on the equilibrium of the average fields. Several examples are discussed, and they show the basic features of the fluctuation fields and the influence of fluctuation fields on the average configurations of magnetic fields. The astrophysical environments are often in the turbulent state, the results of the present paper may be applied to the turbulent plasma where the magnetic field is strong.

国际理论和应用力学联合会全代会(IUTAM General Assembly)会议简介

<正> 1国际理论和应用力学联合会全代会的介绍 国际理论和应用力学联合会(IUTAM)全代会(Gen-eral Assembly)是全世界力学界的最高组织IUTAM的最高权力机构.全代会按照IUTAM诸参加国的等级(分5等)来划分代表名额.只有美国、日本等国由于交纳最高等级的会费而拥有5个代表名额,而我国在代表大会中有4位代表(郑哲敏,何友声,白以龙和杨卫).此外,在决策IUTAM