A scaled boundary node method applied to two-dimensional crack problems

A boundary-type meshless method called the scaled boundary node method(SBNM) is developed to directly evaluate mixed mode stress intensity factors(SIFs) without extra post-processing.The SBNM combines the scaled boundary equations with the moving Kriging(MK) interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter.As a result,the SBNM requires only a set of scattered nodes on the boundary,and the displacement field is approximated by using the MK interpolation technique,which possesses the δ function property.This makes the developed method efficient and straightforward in imposing the essential boundary conditions,and no special treatment techniques are required.Besides,the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction.Therefore,the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip.Numerical examples using the SBNM for computing the SIFs are presented.Good agreements with available results in the literature are obtained.     

Influence of residual surface stress on the fracture of nanoscale piezoelectric materials with conducting cracks

In this paper,we analyze the stress and electric field intensity factors affected by residual surface stress for conducting cracks in piezoelectric nanomaterials.The problem is reduced to a system of non-linear singular integral equations,whose solution is determined by iteration technique.Numerical results indicate that the residual surface stress can significantly alter the crack tip fields at nanometer length scales.Due to the residual surface stress,281he electric field can produce stress around crack tip.This suggests a strong electromechanical coupling crack tip field for nanoscale piezoelectric materials.Such a finding is considerably different from the classical fracture mechanics results.A transit electric field to stress load ratio is identified,for which influences of residual surface stresses vanish.The research is useful for the applications of nanoscale piezoelectric devices.     

动态断裂力学的无网格流形方法

运用无网格流形方法求解动态断裂力学问题.该方法利用单位分解法和有限覆盖技术建立形函数，形函数的建立不受域内不连续的影响，可较好地求解裂纹问题.对于局部化问题，该方法的形函数构造较其他方法更为有效，避免了其他方法在建立试函数时没有考虑不连续尖端的缺点.由于采用有限覆盖技术建立试函数，该方法克服了不连续对试函数的影响，尤其当不连续变得复杂时，更能显示该方法在处理不连续方面的优点.在求解动态断裂力学问题时，弹性动力学积分弱形式的推导采用加权残数法，空间离散采用基于单位分解法的无网格流形方法，时间离散主要采用Newmark法.最后给出两个数值算例，将计算结果与解析解对比，说明该方法的正确性和可行性. 关键词： 有限覆盖 无网格流形方法 动态断裂力学 动态应力强度因子     

轴对称构件受力分析的插值粒子法

面对土木工程与机械工程中广泛存在的轴对称力学问题, 采用具有离散点插值特性的无网格方法形函数, 结合弹性力学空间轴对称问题的最小势能原理, 建立了轴对称构件力学分析的插值粒子法. 本文无网格法方法构造形函数不依赖网格, 也具有像有限元法一样可直接施加边界条件的优点. 本方法能直接获得全域连续应力场, 避免了有限元法应力后处理二次拟合带来的计算误差. 最后通过实例分析, 验证了所建立的无网格方法的有效性.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

磁致伸缩/压电层叠复合材料磁电效应分析

采用有限元分析软件COMSOL5.0建立了三维悬臂梁模型,分析了磁致伸缩/压电/磁致伸缩叠层复合材料的磁电系数α_(ME),并就几何参数对复合结构磁电系数的影响进行了优化分析.首先,利用稳态求解器研究了磁电层状复合结构内部的应力、应变、位移以及电势分布情况,利用瞬态求解分析了磁电复合结构各变量动态分布规律;其次,应用小信号频域分析研究了该结构的谐振频率以及在不同偏置磁场对输出电压的影响,结果表明,随着直流偏置磁场的增加,输出电压逐渐减小.改变复合材料不同层的厚度,分析了磁电层与压电层厚度比t_m/t_p对磁电系数的影响,结果表明,随着厚度比增加,α_(ME)逐渐增大,其增加速率逐渐减小;最后,分析了磁电系数α_(ME)随复合结构面积、长宽比的变化情况.分析表明,α_(ME)随磁电复合结构面积的增加逐渐增加,其增加速率逐渐减小;当磁电复合结构面积恒定时,其磁电系数随长宽比L/W增加表现出先增加后减小的趋势,存在最优值.     

Application of the extended traction boundary element-free method to the fracture of two-dimensional infinite magnetoelectroelastic solid

A novel extended traction boundary element-free method is proposed to analyze the crack problems of two-dimensional infinite magnetoelectroelastic solid.An extended traction boundary integral equation only involving Cauchy singularity is firstly derived.Then,the extended dislocation densities on the crack surface are expressed as the combination of a characteristic term and unknown weight functions,and the radial point interpolation method is adopted to approximate the unknown weight functions.The numerical...     

Analysis on vibration characteristics of large-size rectangular piezoelectric composite plate based on quasi-periodic phononic crystal structure

Based on the theory of composite materials and phononic crystals (PCs), a large-size rectangular piezoelectric composite plate with the quasi-periodic PC structure composed of PZT-4 and epoxy is proposed in this paper. This PC structure can suppress the transverse vibration of the piezoelectric composite plate so that the thickness mode is purer and the thickness vibration amplitude is more uniform. Firstly, the vibration of the model is analyzed theoretically, the electromechanical equivalent circuit diagram of three-dimensional coupled vibration is established, and the resonance frequency equation is derived. The effects of the length, width, and thickness of the piezoelectric composite plate at the resonant frequency are obtained by the analytical method and the finite element method, the effective electromechanical coupling coefficient is also analyzed. The results show that the resonant frequency can be changed regularly and the electromechanical conversion can be improved by adjusting the size of the rectangular piezoelectric plate. The effect of the volume fraction of the scatterer on the resonant frequency in the thickness direction is studied by the finite element method. The band gap in X and Y directions of large-size rectangular piezoelectric plate with quasi-periodic PC structures are calculated. The results show that the theoretical results are in good agreement with the simulation results. When the resonance frequency is in the band gap, the decoupling phenomenon occurs, and then the vibration mode in the thickness direction is purer.     

基于温度效应的无限长压电圆杆纵波分析

利用有限变形理论,以无限长压电圆杆为研究对象,考虑了在横向惯性、等效泊松比效应以及在热电弹耦合共同作用下,基于Hamilton原理,并引入Euler方程推导出压电圆杆的纵向波动方程.采用Jacobi椭圆函数展开法,求解压电圆杆的波动方程和对应的解.最后,通过Matlab软件得到不同波速比下的色散曲线,以及温度场对压电圆...     

A theoretical study of the propagation of Rayleigh waves in a functionally graded piezoelectric material (FGPM)

An exact approach is used to investigate Rayleigh waves in a functionally graded piezoelectric material (FGPM) layer bonded to a semi infinite homogenous solid. The piezoelectric material is polarized when the six fold symmetry axis is put along the propagation direction x₁. The FGPM character imposes that the material properties change gradually with the thickness of the layer. Contrary to the analytical approach, the adopted numerical methods, including the ordinary differential equation (ODE) and the stiffness matrix method (SMM), treat separately the electrical and mechanical gradients. The influences of graded variations applied to FGPM film coefficients on the dispersion curves of Rayleigh waves are discussed. The effects of gradient coefficients on electromechanical coupling factor, displacement fields, stress distributions and electrical potential, are reported. The obtained deviations in comparison with the ungraded homogenous film are plotted with respect to the dimensionless wavenumber. Opposite effects are observed on the coupling factor when graded variations are applied separately. A particular attention has been devoted to the maximum of the coupling factor and it dependence on the stratification rate and the gradient coefficient. This work provides with a theoretical foundation for the design and practical applications of SAW devices with high performance.     

基于虚拟仪器的压电陶瓷驱动技术研究

针对压电陶瓷驱动器的非线性应变,以压电陶瓷微位移驱动原理为基础,分析了扫频激光器腔长控制原理;在硬件时钟的定时下,分析了压电陶瓷的驱动特性,并通过调整驱动信号的步长电压来调节压电陶瓷的非线性。设计了基于虚拟仪器的控制系统,并对线性化方法进行了原理分析和实验。实验结果表明,该测量系统可以使压电陶瓷在其伸长范围内线性地伸长。     

基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.     

基于声粒子分布积分的无网格声场计算方法

小尺度封闭空间内部声场的数值计算是声学设计、噪声控制等领域的关键技术。由于波动声学及几何声学方法计算频率上的限制,中频段声场计算问题一直是个难点。本文以声学无网格法为基础,提出了一种基于声粒子分布积分的无网格声场数值计算方法。文中利用声线跟踪理论计算声场中的声粒子分布,并以某个时间点上的声粒子作为蒙特卡罗法中的积分点,将其应用于无网格法中,从而获得声场中的节点声压。利用该方法对一个矩形封闭空间的中低频声场进行了计算,并与模态叠加法、商用声场计算软件、经典无网格法的结果进行了对比,证明基于声粒子分布积分的无网格声场数值计算方法在中低频段相较于传统基于网格的方法具有更高的精度。     

Identification of the piezoelectric material coefficients using the finite element method with an asymptotic waveform evaluation

A rapid identification of the piezoelectric material constants for a piezoelectric transducer is proposed. The validity of a three-dimensional finite element routine was confirmed experimentally. The asymptotic waveform evaluation (AWE) was adopted for a fast frequency sweep of the finite element analysis. The three-dimensional finite element method with an AWE and a design sensitivity method was used for a material inversion scheme of piezoelectric transducers. In order to confirm the inversion routine of the material constants, the mechanical displacements, which mean the mode shape, were calculated along the vertical and lateral position of the sample transducer.     

Element-free Galerkin (EFG) method for analysis of the time-fractional partial differential equations

The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α (0<α ≤1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.     

外加磁场压电悬臂梁能量采集系统的磁化电流法磁力研究

以外加磁场压电悬臂梁能量采集系统结构为研究对象, 根据磁化电流方法探讨了具有悬臂梁特征的系统结构的磁场作用力及其计算方法, 给出了相应的磁力计算模型, 并将计算结果与实验数据进行了对比. 研究表明, 磁化电流方法导出的磁力计算模型存在偏差, 其磁力计算误差随着磁铁间距缩小而增大. 通过引入悬臂梁末端磁铁的偏转角度, 对磁化电流法计算模型进行改进, 得到合理的外加磁场压电悬臂梁能量采集系统的磁力计算模型, 为该能量采集系统的进一步研究提供了可靠的磁力计算理论依据.     

复合声子结构的声场模拟及噪声-压电转化设计

城市噪声作为一种污染源,给人们的生活带来不便;同时噪声中蕴含的能量没有被合理利用,造成浪费。为了解决此问题,该文设计引入了点缺陷的声子晶体、纤维层、压电材料复合结构,对入射噪声进行有效吸收并在点缺陷处对入射声压有较明显的声压加强。采用有限元和边界元的方法,对该复合结构的降噪发电效果进行数值模拟。结果显示在输入声压为2 Pa时,声子晶体结构带隙1.1 kHz处,吸声系数达到0.6,压电片有最高的输出电压0.5 V,输出功率密度达到308.49μW/cm~3。与传统降噪结构相比能更好地实现对声能的吸收与利用。     

Optical method of caustics applied in viscoelastic fracture analysis

The optical method of caustics is developed here to study the fracture of viscoelastic materials. By adopting a distribution of viscoelastic stress fields near the crack tip, the method of caustics is used to determine the viscoelastic fracture parameters from the caustic patterns near the crack tip. Two viscoelastic materials are studied. These are PMMA and ternary composites of HDPE/POE-g-MA/CaCO₃. The transmitted and reflective methods of caustics are performed separately to investigate viscoelastic fracture behaviors. The stress intensity factors (SIFs) versus time is determined by a series of shadow spot patterns combined with viscoelastic parameters evaluated by creep tests. In order to understand the viscoelastic fracture mechanisms of HDPE/POE-g-MA/CaCO₃ composites, their fracture surfaces are observed by a Scanning Electron Microscope (SEM). The results indicate that the method of caustics can be used to characterize the fracture behaviors of viscoelastic materials and further to optimize the design of polymer composites.     

基于分布源边界点法的局部近场声全息技术

为了克服基于分布源边界点法的近场声全息技术在小全息孔径条件下造成的重建误差问题,提出了基于分布源边界点法的局部近场声全息技术.该技术运用分布源边界点法,采用测得的较小全息面上的声压数据来外推较大全息面上的声压数据,然后用外推的数据进行全息重建.仿真和实验结果验证了采用该技术在小全息孔径条件下进行声源局部重建的有效性.