Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

Anisotropic thermopiezoelectric solids with an elliptic inclusion or a hole under uniform heat flow

The two-dimensional problem of a thermopiezoelectric material containing an elliptic inclusion or a hole subjected to a remote uniform heat flow is studied. Based on the extended Lekhnitskii formulation for thermopiezoelectricity, conformal mapping and Laurent series expansion, the explicit and closed-form solutions are obtained both inside and outside the inclusion (or hole). For a hole problem, the exact electric boundary conditions on the hole surface are used. The results show that the electroelastic fields inside the inclusion or the electric field inside the hole are linear functions of the coordinates. When the elliptic hole degenerates into a slit crack, the electroelastic fields and the intensity factors are obtained. The effect of the heat how direction and the dielectric constant of air inside the crack on the thermal electroelastic fields are discussed. Comparison is made with two special cases of which the closed solutions exist and it is shown that our results are valid.

Elastic Waves in a Hybrid Multilayered Piezoelectric Plate

An analytical-numerical method is presented for analyzing dispersion and characteristic surface of waves in a hybrid multilayered piezoelectric plate. In this method, the multilayered piezoelectric plate is divided into a number of layered elements with three-nodal-lines in the wall thickness, the coupling between the elastic field and the electric field is considered in each element. The associated frequency dispersion equation is developed and the phase velocity and slowness, as well as the group velocity and slowness are established in terms of the Rayleigh quotient. Six characteristic wave surfaces are introduced to visualize the effects of anisotropy and piezoelectricity on wave propagation. Examples provide a full understanding for the complex phenomena of elastic waves in hybrid multilayered piezoelectric media.

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacement D-2 induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with real H, the normal electric displacement D-2 is constant along the crack faces and electric field E-2 has the singularity ahead of the crack tip and a jump across the interface.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

An Electrowetting Model For Rough Surfaces Under Low Voltage

Electrowetting is one of the most effective methods to enhance wettability. A significant change of contact angle for the liquid droplet can result from the surface microstructures and the external electric field, without altering the chemical composition of the system. During the electrowetting process on a rough surface, the droplet exhibits a sharp transition from the Cassie-Baxter to the Wenzel regime at a low critical voltage. In this paper, a theoretical model for electrowetting is put forth to describe the dynamic electrical control of the wetting behavior at the low voltage, considering the surface topography. The theoretical results are found to be in good agreement with the existing experimental results. (c) Koninklijke Brill NV, Leiden, 2008.

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.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total 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. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for 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. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

A new model for the floating potential of fine particles in plasma

In the plasma processing of ultrafine particles of material, the heat transfer and force are considerably affected by particle charging. In this communication a new model, including thermal electron emission and incorporating the effect of electric field near the particle surface, is developed for metallic spherical particles under the condition of a thin plasma sheath. Based on this model, the particle floating potential, and thus the heat transfer and force, can be detemined more accurately and more realistically than previously.

Experimental Characterization Of Electrical Current Leakage In Poly(Dimethylsiloxane) Microfluidic Devices

Poly(dimethylsiloxane) (PDMS) is usually considered as a dielectric material and the PDMS microchannel wall can be treated as an electrically insulated boundary in an applied electric field. However, in certain layouts of microfluidic networks, electrical leakage through the PDMS microfluidic channel walls may not be negligible, which must be carefully considered in the microfluidic circuit design. In this paper, we report on the experimental characterization of the electrical leakage current through PDMS microfluidic channel walls of different configurations. Our numerical and experimental studies indicate that for tens of microns thick PDMS channel walls, electrical leakage through the PDMS wall could significantly alter the electrical field in the main channel. We further show that we can use the electrical leakage through the PDMS microfluidic channel wall to control the electrolyte flow inside the microfluidic channel and manipulate the particle motion inside the microfluidic channel. More specifically, we can trap individual particles at different locations inside the microfluidic channel by balancing the electroosmotic flow and the electrophoretic migration of the particle.

变形镜对千赫兹掺钛蓝宝石飞秒激光频谱相位补偿

尝试用光谱展宽的方法从频谱相位(而非时域相位)的角度利用变形反射镜来补偿1kHz飞秒激光系统输出光路的频谱相位畸变，从而提高飞秒激光脉冲的时域强度衬比度，改善其光束质量。频谱相位补偿实验是在一台1kHz掺钛蓝宝石飞秒激光系统输出光路中，针对超短脉冲光束通过传输介质后的频谱相位畸变，引入变形反射镜进行补偿。应用频谱相位干涉直接电场重构(SPIDER)方法和仪器作为测量手段，建立了一套相位测量补偿系统。实验结果表明用变形反射镜可使激光脉冲的相位畸变得到较好的补偿，脉冲的光束质量得到改善。这种方法的主要思想就是

双色激光场中激光强度和脉宽对线性多原子分子离子电离的影响

利用强激光场电离和离解分子来研究分子激发态的波包结构是强场物理的重要研究方向。利用短时指数传播子对称分割法和快速傅里叶变换技术。数值求解了一维含时Schr（oe）dinger方程，探讨了双色激光场中激光的基波和谐波强度之间的不同配比以及脉宽对线性多原子分子离子电离的影响。理论计算结果表明：基波和谐波的相对相位为π时，尽管随着激光的基波和谐波强度之间配比的变化，电离几率随原子间距变化的趋势基本保持不变，但在一定的激光基波强度下(1．2×10^13～1．2×10^15W／cm^2)，激光基波强度的变化可以明显

外加高压电场下空气中激光等离子体通道寿命研究

通过对飞秒激光在空气中产生的等离子体通道两端外加高压,来研究通道的寿命变化情况。实验得到,当在等离子体通道两端外加高压时(350 kV/m),等离子体通道寿命延长了近3倍。理论模拟和分析结果表明在外加电场条件下,碰撞电离得到增强,吸附作用相对减弱,解离复合系数随着电子平均能量的增加而下降的趋势更为剧烈,这进一步引起了等离子体通道寿命的延长。实验结果与理论分析共同表明了利用外加电场对空气中激光等离子体通道寿命进行延长的可行性。

Fabrication of microelectrodes deeply embedded in LiNbO3 using a femtosecond laser

We present a novel technique to fabricate deeply embedded microelectrodes in LiNbO3 using femtosecond pulsed laser ablation and selective electroless plating. The fabrication process mainly consists of four steps, which are (1) micromachining of microgrooves on the surface of LiNbO3 by femtosecond laser ablation; (2) formation of AgNO3 films on substrates; (3) scanning the femtosecond laser beam in the fabricated microgrooves for modi. cation of the inner surfaces; and (4) electroless copper plating. The void-free electroless copper plating is obtained with appropriate cross section of microgrooves and uniform initiation of the autocatalytic deposition on the inner surface of grooves. The dimension and shape of the microelectrodes could be accurately controlled by changing the conditions of femtosecond laser ablation, which in turn can control the distribution of electric field inside LiNbO3 crystal for various applications, opening up a new approach to fabricate three-dimensional integrated electro-optic devices. (C) 2008 Elsevier B. V. All rights reserved.