The mode III stress/electric intensity factors and singularities analysis for edge-cracked circular piezoelectric shafts

This paper introduces a new analytical method to determine the stress and electric intensity factors for edge-cracked circular piezoelectric shafts using Hamiltonian formalism. The singularities near the crack tip are represented in terms of exponential series that can show the boundary layer effects effectively. A symplectic system is established directly by introducing dual vectors in terms of the symplectic eigenfunctions. The coefficients of the series are determined from the lateral boundary conditions along the crack faces and the outer boundary conditions along the exterior geometric domain. The intensity factors are determined by the first two coefficients of non-zero eigenvalue solutions. Numerical examples for various boundary conditions are given. The influencing parameters on the intensity factors are investigated.     

含裂纹纳米板振动问题的哈密顿体系方法

基于裂纹处范德华力效应,采用非局部弹性理论构造纳米板模型,并通过导入哈密顿体系建立含裂纹纳米板振动问题的对偶正则控制方程组。在全状态向量表示的哈密顿体系下,将含裂纹纳米板的固有频率和振型问题归结为广义辛本征值和本征解问题。利用哈密顿体系具有的辛共轭正交关系,得到问题解的级数解析表达式。结合边界条件,得到固有频率与辛本征值的代数方程关系式,进而直接给出固有频率的表达式。数值结果表明,非局部尺寸参数和裂纹长度对纳米板振动的各阶固有频率有直接的影响。对比表明,辛方法是准确且可靠的,可为工程应用提供依据。     

Determination of stress intensity factors by the finite element discretized symplectic method

When rewriting the governing equations in Hamiltonian form, analytical solutions in the form of symplectic series can be obtained by the method of separation of variable satisfying the crack face conditions. In theory, there exists sufficient number of coefficients of the symplectic series to satisfy any outer boundary conditions. In practice, the matrix relating the coefficients to the outer boundary conditions is ill-conditioned unless the boundary is very simple, e.g., circular. In this paper, a new two-level finite element method using the symplectic series as global functions while using the conventional finite element shape functions as local functions is developed. With the available classical finite elements and symplectic series, the main unknowns are no longer the nodal displacements but are the coefficients of the symplectic series. Since the first few coefficients are the stress intensity factors, post-processing is not required. A number of numerical examples as well as convergence studies are given.     

An analytical symplectic approach to the vibration analysis of orthotropic graphene sheets

A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach.A Hamiltonian system is established by introduc-ing a total unknown vector consisting of the displacement amplitude,rotation angle,shear force,and bending moment. The high-order governing differential equation of the vibra-tion of SLGSs is transformed into a set of ordinary differential equations in symplectic space.Exact solutions for free vibra-tion are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary con-ditions. Vibration modes are expressed in terms of the symplectic eigenfunctions.In the numerical examples,com-parison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given.A para-metric study of the natural frequency is also included.     

STRESS SINGULARITY ANALYSIS AT CRACK TIP ON BI-MATERIAL INTERFACES BASED ON HAMILTONIAN PRINCIPLE

In this paper, the stress singularity analysis at the crack tip on elastic bi-material inter-faces is considered. The governing equations of plane elasticity in sectorial domain are derived to be inHamiltonian form via variable substitution and variational principle. The methods of separation ofvariables and conjugate symplectic eigen-function expansion are developed to solve the equations insectorial domain. The general formulae for the solution of stress singularities at the crack tip on bi-ma-terial interfaces are put forward, and a new solution technique for fracture problems is presented.     

Mixed-mode thermal stress intensity factors from the finite element discretized symplectic method

A finite element discretized symplectic method is introduced to find the thermal stress intensity factors (TSIFs) under steady-state thermal loading by symplectic expansion. The cracked body is modeled by the conventional finite elements and divided into two regions: near and far fields. In the near field, Hamiltonian systems are established for the heat conduction and thermoelasticity problems respectively. Closed form temperature and displacement functions are expressed by symplectic eigen-solutions in polar coordinates. Combined with the analytic symplectic series and the classical finite elements for arbitrary boundary conditions, the main unknowns are no longer the nodal temperature and displacements but are the coefficients of the symplectic series after matrix transformation. The TSIFs, temperatures, displacements and stresses at the singular region are obtained simultaneously without any post-processing. A number of numerical examples as well as convergence studies are given and are found to be in good agreement with the existing solutions.     

Hamiltonian principle based stress singularity analysis near crack corners of multi-material junctions

This paper presents a new method for the stress singularity analysis near the crack corners of a multi-material junctions. The stress singularities near the crack corners of multi-dissimilar isotropic elastic material junctions are studied analytically in terms of the methods developed in Hamiltonian system. The governing equations of plane elasticity in a sectorial domain are derived in Hamiltonian form via variable substitution and variational principle respectively. Both of the methods of global state variable separation and symplectic eigenfunction expansion are used to find the analytical solution of the problem. The relationships among the state vectors in different material spaces are obtained by means of coordinate transformation and consistent conditions between the two adjacent domains. The expression of the original problem is thus changed into a new form where the solutions of symplectic generalized eigenvalues and eigenvectors are needed. The closed form of expressions is established for the stress singularity analysis near the corner with arbitrary vertex angles. Numerical results are presented with several chosen angles and multi-material constants. To show the potential of the new method proposed, a semi-analytical finite element is furthermore developed for the numerical analysis of crack problems.     

双层石墨烯系统稳态受迫振动问题中的辛方法

针对悬臂石墨烯系统提出一种求解其稳态受迫振动问题的辛解析方法。基于Eringen非局部理论,将石墨烯层板受迫振动问题导入哈密顿体系。采用边界条件分解技术,将问题化为三种边界条件的子问题。通过辛解析方法,得到由辛本征值和辛本征解表示的双层石墨烯系统受迫振动问题的解析解表达式。数值结果表明,辛本征解级数具有很好的收敛性和精度,并与文献结果吻合;在一定的外载激励下可发生同向振动模式和反向振动模式;在一定的参数下,得到一些新的现象和结论。     

哈密顿体系与弹性楔体问题

将哈密体系引入到级坐标下的弹性力学楔体问题,利用该体系辛空间的性质,将问题化为本征值和本征向量求解上,得到了完备的解空间,从而改变了弹性力学传统的拉格朗日体系以应力函数为特征的半逆法的讨论去解决该类问题的思路,给出了一条求解该类问题的直接法。     

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

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

Crack energy density of a piezoelectric material under general electromechanical loading

Crack energy density is considered and used as a possible fracture parameter in piezoelectricity under arbitrary electromechanical remote loads. The closed-form solution of a crack in a piezoelectric infinite plate subjected to general static electromechanical loading is obtained through a method alternative to the more common Stroh’s formalism. This analytical method, which is based on the spectral theorem of linear algebra, involves a transformation of similarity induced by the fundamental matrix in order to express the equations governing the problem in terms of complex potentials. The application of the mechanical boundary condition of stress-free crack and of one of the three considered electric boundary conditions (impermeable, permeable or semipermeable) leads then to the formulation of a Hilbert problem whose solution yields the stress and displacement fields. The crack energy density factors for mixed mode are then calculated under different mechanical and electrical loadings, as well as under different electric boundary conditions. The non-singular terms of the stress expressions are retained as well. The definition of the minimum energy density fracture criterion, as proposed by Sih, is given, and the influence of load biaxiality and positive or negative applied electric field on the criterion results is analyzed. The prediction of the incipient branching angle as from the energy density approach is also compared to that arising from the maximum circumferential stress theory for a mixed mode loading condition. Numerical results and graphs are presented and discussed for a PZT-4 piezoelectric ceramic.     

Closed form stress distribution in 2D elasticity for all boundary conditions

This paper applies a Hamiltonian method to study analytically the stress dis- tributions of orthotropic two-dimensional elasticity in(x,z)plane for arbitrary boundary conditions without beam assumptions.It is a method of separable variables for partial differential equations using displacements and their conjugate stresses as unknowns.Since coordinates(x,z)can not be easily separated,an alternative symplectic expansion is used. Similar to the Hamiltonian formulation in classical dynamics,we treat the x coordinate as time variable so that z becomes the only independent coordinate in the Hamiltonian ma- trix differential operator.The exponential of the Hamiltonian matrix is symplectic.There are homogenous solutions with constants to be determined by the boundary conditions and particular integrals satisfying the loading conditions.The homogenous solutions consist of the eigen-solutions of the derogatory zero eigenvalues(zero eigen-solutions) and that of the well-behaved nonzero eigenvalues(nonzero eigen-solutions).The Jordan chains at zero eigenvalues give the classical Saint-Venant solutions associated with aver- aged global behaviors such as rigid-body translation,rigid-body rotation or bending.On the other hand,the nonzero eigen-solutions describe the exponentially decaying localized solutions usually ignored by Saint-Venant's principle.Completed numerical examples are newly given to compare with established results.     

四边固支矩形薄板自由振动的哈密顿解析解

在哈密顿体系中利用辛几何方法求解了四边固支矩形薄板自由振动问题的解析解。首先,从基本方程出发,将问题表示成Hamilton正则方程,然后利用辛几何方法导出本征值问题,从而得到本征函数解,使之满足边界条件;再由方程组有非零解的条件,最终推导出四边固支矩形薄板的自振频率方程,得到频率的解析解。计算了不同长宽比情况下四边固支矩形薄板的频率,结果与已有文献完全一致。该解法有望推广至更多尚未得到解析解的矩形板的振动问题。     

An exact symplectic solution for the dynamic analysis of shear deformable antisymmetric angle-ply laminated plates

From the mixed variational principle, by the selection of the state variables and its dual variables, the Hamiltonian canonical equation for the dynamic analysis of shear deformable antisymmetric angle-ply laminated plates is derived, leading to the mathematical frame of symplectic geometry and algorithms, and the exact solution for the arbitrary boundary conditions is also derived by the adjoint orthonormalized symplectic expansion method. Numerical results are presented with the emphasis on the effects of length/thickness ratio, arbitrary boundary conditions, degrees of anisotropy, number of layers, ply-angles and the corrected coefficients of transverse shear.     

紧凑拉伸试件应力强度因子的解析-变分解法以及有关系统计算结果

本文推导了在材料断裂性能测试中常见的受钉传载荷含边缘裂纹试件应力与位移的函数项级数表达式。该级数逐项满足弹性力学所有基本方程、裂纹表面边界条件与绕钉孔的合力平衡条件以及位移单值条件。通过以最小势能原理为基础的变分方程满足其余的静力边界条件,从而求解级数中的待定系数并确定应力强度因子。计算结果表明,级数收敛迅速、正确,计算节省机时,简化数据准备工作。本文还通过计算指出了目前通用的有关矩形紧凑拉伸试件应力强度因子计算公式与曲线的不准确性并且给出了正确、系统的计算曲线,同时还提供了圆形紧凑拉伸试件系统的计算结果。     

反对称铺设层合板动力问题的Hamilton体系及辛几何解法

基于ＹＮＳ层合格理论，建立反对称铺设层合板动力问题的Ｈａｍｉｌｔｏｎ正则方程，并采用共轭辛正交归一关系给出一对边简支，另一对边为任意支承层合板自振频率的精确解，数值算例讨论了长宽比，铺设角，层数及剪切修正系数的影响。     

Numerical methods for solving singular integral equations obtained by fracture mechanical analysis of cracked wedge简

The numerical solutions to the singular integral equations obtained by the fracture mechanical analyses of a cracked wedge under three different conditions are considered. The three considered conditions are: (i) a radial crack on a wedge with a non-finite radius under the traction-traction boundary condition, (ii) a radial crack on a wedge with a finite radius under the traction-traction boundary condition, and (iii) a radial crack on a finite radius wedge under the traction-displacement boundary condition. According to the boundary conditions, the extracted singular integral equations have different forms. Numerical methods are used to solve the obtained coupled singular integral equations, where the Gauss-Legendre and the Gauss-Chebyshev polynomials are used to approximate the responses of the singular integral equations. The results are presented in figures and compared with those obtained by the analytical response. The results show that the obtained Gauss-Chebyshev polynomial response is closer to the analytical response.     

弹性矩形板问题的Hamilton正则方程

为了采用辛算法求出弹性矩形板问题的解析解，中直接从弹性矩形板的控制方程出发推导了弹性矩形板，其中包括弹性矩形薄板和厚板问题以及弹性地基上矩形薄板和厚板问题的Hamilton正则方程，为利用辛几何方法求出任意边界条件下这类问题的理论解奠定了基础．     

Interfacial fracture analysis of bonded dissimilar strips with a functionally graded interlayer under antiplane deformation

This paper provides the solution to the problem of dissimilar, homogeneous semi-infinite strips bonded through a functionally graded interlayer and weakened by an embedded or edge interfacial crack. The bonded system is assumed to be under antiplane deformation, subjected to either traction-free or clamped boundary conditions along its bounding planes. Based on the Fourier integral transform, the problem is formulated in terms of a singular integral equation which has a simple Cauchy kernel for the embedded crack and a generalized Cauchy kernel for the edge crack. In the numerical results, the effects of geometric and material parameters of the bonded system on the crack-tip stress intensity factors are presented in order to quantify the interfacial fracture behavior in the presence of the graded interlayer.     

Non-local theory solution for a Mode I crack in piezoelectric materials

In this paper, the non-local theory of elasticity is applied to obtain the behavior of a Griffith crack in the piezoelectric materials subjected to a uniform tension loading. The permittivity of the air in the crack is considered. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, in which the unknown variables are the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surfaces are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the crack length, the materials constants, the electric boundary conditions and the lattice parameter on the stress and the electric displacement fields near the crack tips. It can be obtained that the effects of the electric boundary conditions on the electric displacement fields are large. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularities are present at the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips, thus allowing us to use the maximum stress as a fracture criterion.