A stationary free boundary problem modeling electrostatic MEMS

We investigate a free boundary problem describing small deformations in a membrane based model of electrostatically actuated microelectromechanical systems (MEMS). The existence of stationary solutions is established for small voltage values and non-existence is obtained for high voltage values. We give a justification of the widely studied narrow-gap model by showing that steady state solutions of the free boundary problem converge toward stationary solutions of the narrow-gap model when the aspect ratio of the device tends to zero.     

压电弹性层合板静力机耦合特性的解析解

应用经典权壳理论和压电理论对层合板结构进行简化。对四边简支压电层合板在不同电学边界条件下,包括四边电学开路和电学短路。上下表面受外电压及无外加电压作用,进行分析,求得了电势和挠度的解析表达式,给出了压电层和基体的挠度、电势分布图。     

离子选择性表面电对流的多块LB模拟

离子选择性表面(如纳米通道、离子交换膜等)复杂的动力学现象为微纳流控技术的发展提供了新思路. 向带有离子选择性表面的电解质溶液施加电压, 通过液体的电流密度会经历复杂的非线性变化过程; 当电压超过某一临界值时会引发对流现象, 这种流动被称为第二类电渗或离子选择性表面的电对流, 关于此类问题的数值求解引发了大量的研究. 本文提出一种基于多块网格加密的格子玻尔兹曼方法(lattice Boltzmann method, LBM)的数值模型, 用于模拟第二类电渗流动. 结合该算法, 给出了求解流动、电势和离子浓度的网格信息交换方程, 较好地解决了此类问题中大浓度梯度边界对计算分辨率的要求. 利用该数值模型模拟获得的电流?电压特性曲线先随着电压升高而迅速增大, 随后达到饱和状态, 与理论解吻合良好. 此外, 模拟结果还表明, 当流动发生后, 相对低电压下的流动倾向于形成大涡且流动呈指数趋势增强; 而较大电压则会先激发多个小涡, 并逐渐合并为大涡流动, 且大涡流动有更高的离子输运效率. 此外, 除了模拟离子选择性表面的电对流现象, 本文提出的数值格式还可拓展到其他电流体动力学问题的模拟.      

功能梯度压电材料圆环在位移和反对称电势边值条件下的解析解

采用了在径向极化情况下横观各向同性的线性本构关系,考虑了材料性质沿径向的梯度分布,对功能梯度压电材料圆环在给定的位移和电势边界条件下,导出了问题的一般解．推导了外壁固定、接地,内壁沿垂向有一微小位移、电势为反对称分布问题的解析解,并计算了该问题在位移和电势作用情况下的位移、电势在不同梯度分布时的数值结果．     

Three-dimensional non-axisymmetric behavior of a penny shaped crack in a piezoelectric strip subjected to in-plane loads

The behavior of a penny shaped crack in a three-dimensional piezoelectric ceramic strip under non-axisymmetric in-plane normal mechanical and electrical loads is analyzed based on the continuous electric boundary conditions of the crack surface. The potential theory, Hankel transform and Fourier series are used to obtain the systems of dual integral equations, which are then expressed as Fredholm integral equations. The singular mechanical and electric fields and all mode-I field intensity factors are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic are shown graphically for an uniform load and a pair of concentrated load, respectively.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

电迁移诱发表面扩散下沿晶微裂纹的演化

随着微电子技术的快速发展,铜内连导线的失效问题日益受到关注。基于表面扩散和蒸发—凝结的经典理论及其弱解描述,建立描述电迁移下微结构演化的有限单元法,对铜材料内沿晶微裂纹在电迁移诱发表面扩散下的不稳定外形演化进行了有限元模拟。详细讨论了电场、形态比和晶界能与表面能比值对沿晶微裂纹演化的影响。结果表明：沿晶微裂纹在沿着晶界和电场方向发生漂移的过程中存在分节与不分节两种形貌演化分叉趋势,且裂腔分节存在临界形态比 和临界电场值 。当 或 时,沿晶微裂纹会沿着晶界分节成一大一小两个小沿晶微裂纹。沿晶微裂纹分节时间随形态比和电场的增大而减小,即形态比和电场的增大都使沿晶微裂纹加速分节。而临界电场值 随着形态比的增大而减小;临界形态比 随电场的增大而减小。也就是说形态比和电场的增大还将有助于沿晶微裂纹分节。此外,沿晶微裂纹分节时间要比晶内微裂纹的小,即晶界的存在有助于加速裂腔分节。     

Modelling of finite piezoelectric patches: Comparing an approximate power series expansion theory with exact theory

Plate equations for a plate consisting of one elastic layer and one piezoelectric layer with an applied electric voltage have previously been derived by use of power series expansions of the field variables in the thickness coordinate. These plate equations are here evaluated by the consideration of a time harmonic 2D vibration problem with finite layers. The boundary conditions at the sides of the layers then have to be considered. Numerical comparisons of the displacement field are made with solutions from two other theories; a solution with equivalent boundary conditions for a thin piezoelectric layer applied on an elastic plate and an exact solution based on Fourier series expansions. The two approximate theories are shown to be equally good and they both yield accurate results for low frequencies and thin plates.     

Asymptotic analysis of an impermeable crack in an electrostrictive material subjected to electric loading

The simple asymptotic problem of an impermeable crack in an electrostrictive ceramic under electric loading is analyzed. Closed form solutions of elastic fields are obtained by using the complex function theory. It is found that the K_I-dominant region is very small compared to the electric saturation zone. A fracture parameter for an electrostrictive material subjected to electric loading is discussed. In order to investigate the influence of the transverse electric displacement on fracture behavior under the small-scale conditions, we also consider the modified boundary layer problem of a crack in an electrostrictive material. Analytic solutions of electric displacement fields for the asymptotic problem are obtained based on the nonlinear dielectric theory from a modified boundary layer analysis. The shape of the electric displacement saturation zone is shown to depend on the transverse electric displacement. Stress intensity factors induced by the electrostrictive strains are evaluated using the nonlinear solution of the electric displacements. It is found that the transverse electric displacement affects strongly the variation of the mode mixity.     

Manipulation of particles using dielectrophoresis

A numerical scheme based on the distributed Lagrange multiplier method (DLM) is used to study the motion of particles of electrorheological suspensions subjected to non-uniform electric fields. At small Reynolds number, the time taken by the particles to collect at the minimums or maximums of the electric field is primarily determined by a parameter defined to be the ratio of the dielectrophoretic and viscous forces. Simulations show that in non-uniform electric fields the collection time is also influenced by a parameter defined by the ratio of the electrostatic particle–particle interaction and dielectrophoretic forces. The collection time decreases as this parameter decreases because when this parameter is less than one, particles move to the regions of high or low electric field regions individually. However, when this parameter is greater than one, particles regroup into chains which then move toward the electric field maximums or minimums without breaking. It is also shown that when the real part of the Clausius–Mosotti factor (β) is negative the positions of the local minimums of the electric field, and thus also the locations where particles collect, can be modified by changing the electric potential boundary conditions.     

Electrical nonlinear anti-plane shear crack in a functionally graded piezoelectric strip

The electrical nonlinear behavior of an anti-plane shear crack in a functionally graded piezoelectric strip is studied by using the strip saturation model within the framework of linear electroelasticity. The analysis is conducted on the electrically unified crack boundary condition with the introduction of the electric crack condition parameter that can describe all the electric crack boundary condition in accordance with the aspect ratio of an ellipsoidal crack and the permittivity inside the crack, in particular, including traditional permeable and impermeable crack boundary conditions. The resulting mixed boundary value problem is analysed and near tip field is obtained by using the integral transform techniques. Numerical results for the normalized five kinds of energy release rates under the small scale electric saturation condition are presented and compared to show the influences of the electric crack condition parameter with the variation of the ellipsoidal crack parameters, electric loads, functionally graded piezoelectric material gradation, crack length, electromechanical coupling coefficient, and crack location. It reveals that there are considerable differences between the results obtained from the traditional electric crack models and those obtained from the current unified crack model.     

Anti-electroviscous effect of near-surface 5CB liquid crystal and its boundary lubrication property

The rheological properties of the near-surface 5CB (4-pentyl-4′-cyanobiphenyl) liquid crystal (LC) under an external electric field (EEF) are investigated using a dynamic quartz crystal microbalance (QCM). The real-time film-forming process and shearing response of the EEF-induced absorbed LC liquid are studied. The results indicate that the EEF-induced adsorbed LC is composed of a bulk liquid layer and a near-surface boundary layer. Within the boundary layer, the nematic 5CB LC presents an anti-electroviscous effect which may be due to the EEF-induced ordered structure. According to the estimation from QCM measurements, this near-surface layer is about 100 nm thick under certain electric strength and is irreversible even after the EEF is removed. Based on a QCM model, the near-surface LC layer presents a decreased and irreversible viscosity as the EEF voltage increases against the reversible electro-viscous effect of the bulk 5CB liquid crystal measured by rheometer and Raman spectrum measurements. The anti-electroviscous effect of the near-surface 5CB layer is also proved by an improved boundary lubrication property tested on a tribometer. The coefficient of friction of 5CB LC after a preliminary induction of EEF is the lowest one compared to those without EEF and during the application of EEF. The unique anti-electroviscous property of near-surface 5CB LC revealed in this article suggests a potential method to actively reduce shear resistance in boundary lubrication and in microfluidics.     

Effect of the electric field on the wave flow regimes of a thin film of a viscous dielectric fluid

Waves on the surface of a thin film of a viscous dielectric fluid flowing down the inner surface of one plate of a plane capacitor with alternating voltage applied is considered. It is shown that the volume forces acting from the inhomogeneous electric field are negligibly small in the case of long waves, and the influence of the electric field reduces to the influence of additional pressure onto the film surface. A model equation for determining the deviation of the film thickness from the undisturbed value is derived in the long-wave approximation. Some numerical solutions of this equation are given.     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

Buckling and post-buckling analyses of size-dependent piezoelectric nanoplates

This paper attempts to investigate the buckling and post-buckling behaviors of piezoelectric nanoplate based on the nonlocal Mindlin plate model and von Karman geometric nonlinearity. An external electric voltage and a uniform temperature rise are applied on the piezoelectric nanoplate. Both the uniaxial and biaxial mechanical compression forces will be considered in the buckling and post-buckling analysis. By substituting the energy functions into the equation of the minimum total potential energy principle,the governing equations are derived directly, and then discretized through the differential quadrature(DQ) method. The buckling and post-buckling responses of piezoelectric nanoplates are calculated by employing a direct iterative method under different boundary conditions. The numerical results are presented to show the influences of different factors including the nonlocal parameter, electric voltage,and temperature rise on the buckling and post-buckling responses.     

电活性聚合物圆柱壳静态与动态电压下的响应及稳定性

摘要：在电活性聚合物圆柱壳内外表面施加电压,圆柱壳会变薄并且伸长,因此相同的电压会在圆柱壳内产生更大的电场。这个正反馈可能使圆柱壳厚度不断变薄,最终导致其失稳破坏。本文研究了电活性聚合物圆柱壳在静态和周期电压作用下的响应及稳定性问题。采用neo-Hookean材料模型得到描述圆柱壳表面运动的非线性常微分方程。给出了圆柱壳在不同厚度和边界条件下外加电压随圆柱壳变形的变化曲线,结果表明存在一个临界电压,当外加电压大于这一临界值时,圆柱壳将被破坏。同时,也讨论了厚度和边界条件对临界电压的影响。圆柱壳在正弦周期电压作用下,其运动随时间的变化是周期性的或拟周期性的非线性振动。给出了圆柱壳振动固有频率的计算结果,采用打靶法得到圆柱壳振动的周期解,并且用数值法研究了周期解的稳定性。采用数值仿真得到圆柱壳振动振幅随外加动态电压激励频率的变化曲线,结果表明圆柱壳会发生多频共振,共振时圆柱壳振幅发生跳跃,导致圆柱壳失稳破坏。最后给出共振点临近点的振动曲线和相图,并对其振动特性进行讨论。     

Modeling and dynamic analysis of a magnetically actuated butterfly valve

In this work, a magnetically actuated butterfly valve is considered and a novel and accurate mathematical model is derived. The equilibrium of the system is investigated and the effects of the inlet velocity and direct current voltage (DC) on the stable rotation angle of the valve are presented. Considering a time periodic perturbation arising from electric circuit, the effects of the operating angle, inlet velocity, and driving parameters on the periodic and chaotic dynamics of the system are investigated. It is observed that, for an opening angle less than the cut-off angle, there exists a unique DC voltage for a stable equilibrium. The stability of this equilibrium depends nonlinearly on the inlet velocity and the seating torque. An expression is derived for the threshold value for the stability of the valve. Under periodic voltage, the inlet velocity and stable angle induce a backward shift on the resonant frequency, and jump phenomena and subharmonics are observed for some values of the driving amplitude. The highest amplitudes of vibration are detected for a fully open valve, for an almost closed valve, and for a valve with large inlet velocity. Using bifurcation diagrams and Lyapunov exponents, it is shown that the system exhibits a route to chaos with windows of period doubling and unbounded motion. Some guidance for design of magnetically actuated butterfly valves is proposed as well as recommendations for future work.     

Exact solutions for the plane problem in piezoelectric materials with an elliptic or a crack

Based on the complex potential approach, the two-dimensional problems in a piezoelectric material containing an elliptic hole subjected to uniform remote loads are studied. The explicit, closed-form solutions satisfying the exact electric boundary condition on the hole surface are given both inside and outside the hole. When the elliptic hole degenerates into a crack, the field intensity factors are obtained. It is shown that the stress intensity factors are the same as that of isotropic material, while the electric displacement intensity factor depends on both the material properties and the mechanical loads, but not on the electric loads. In other words, the uniform electric loads have no influence on the field singularities. It is also shown that the impermeable crack assumption used previously to simply the electric condition is not valid to crack problems in piezoelectric materials.     

A new two-dimensional model for electro-mechanical response of thick laminated piezoelectric actuator

This paper investigates the electro-mechanical behaviour of a thick, laminated actuator with piezoelectric and isotropic lamina under externally applied electric loading using a new two-dimensional computational model. The elastic core is relatively thick and thus it is modelled by Timoshenko thick-beam theory. Although the piezoelectric lamina is a beam-like layer, it is formulated via a two-dimensional model because of not only the strong electro-mechanical coupling, but also of the presence of a two-dimensional electric field. It is shown in this paper that a one-dimensional model for the piezoelectric beam-like layer is inadequate. The piezoelectric model is constructed within the scope of linear piezoelectricity. The actuation response is induced through the application of external electric voltage. Under the strong coupling of elasticity and electricity, the strain energy and work of electric potential are presented. The electro-mechanical response of the laminated Timoshenko beam is formulated and determined via a variational energy principle. Numerical examples presented illustrate convincing comparison with finite element solutions and existing published data. New numerical solutions are also presented to investigate the geometric effect on the electro-mechanical bending behaviour.     

Synchrony and lag synchrony on a neuron model coupling with time delay

A neuron model of the Morris and Lecar form is investigated, which is composed of two individuals and is considered to be functioned by the gap junction coupling. When the level of the reversal potential in the calcium ion channel is small, neurons adjust their activities to the common asymptotic states. However, if we increase the level of the reversal voltage in the calcium ion channel, the exact synchrony firing of neurons is produced. Patterns of synchrony activity and the stability are observed to vary with the choice of time delay, which also enhances the multi-variety of the spike bursting firing rhythm. The lag synchrony of time trajectories of the voltage is illustrated near the boundary of the synchrony regime.