Singularities in 2D anisotropic potential problems in multi-material corners: Real variable approach

An analysis of singular solutions at corners consisting of several different homogeneous wedges is presented for anisotropic potential theory in plane. The concept of transfer matrix is applied for a singularity analysis first of single wedge problems and then of multi-material corner problems. Explicit forms of eigenequations for evaluation of singularity exponent in the case of multi-material corners are derived both for all combinations of homogeneous Neumann and Dirichlet boundary conditions at faces of open corners and for multi-material planes with singular interior points. Perfect transmission conditions at wedge interfaces are considered in both cases. It is proved that singularity exponents are real for open anisotropic multi-material corners, and a sufficient condition for the singularity exponents to be real for anisotropic multi-material planes is deduced. A case of a complex singularity exponent for an anisotropic multi-material plane is reported, apparently for the first time in potential theory. Simple expressions of eigenequations are presented first for open bi-material corners and bi-material planes and second for a crack terminating at a bi-material interface, as examples of application of the theory developed here. Analytical solutions of these eigenequations are presented for interface cracks with any combination of homogeneous boundary conditions along the interface crack faces, and also for a special case of a crack perpendicular to a bi-material interface. A numerical study of variation of the singularity exponent as a function of inclination of a crack terminating at a bi-material interface is presented.     

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.     

STRESS SINGULARITY ANALYSIS OF MULTI-MATERIAL WEDGES UNDER ANTIPLANE DEFORMATION

With the help of the coordinate transformation technique, the symplectic dual solving system is developed for multi-material wedges under antiplane deformation. A virtue of present method is that the compatibility conditions at interfaces of a multi-material wedge are expressed directly by the dual variables, therefore the governing equation of eigenvalue can be derived easily even with the increase of the material number. Then, stress singularity on multi-material wedges under antiplane deformation is investigated, and some solutions can be presented to show the validity of the method. Simultaneously, an interesting phenomenon is found and proved strictly that one of the singularities of a special five-material wedge is independent of the crack direction.     

多材料交接点裂纹无摩擦奇性应力场

在文献「１」基础上，针对工程中难于求解的多材料交接点裂纹尖应力奇性分析问题，在于哈密顿原理，通过分离变量与共轭辛本征函数地求解，利用材料间的界面连接条件与坐标变换关系，建立了应力奇性与本征函数求解的解析表达式。由于采取裂纹面接触区模型，因而不再发生振荡奇异性。     

基于哈密尔顿体系的裂纹尖端应力奇性分析及计量

对弹性平面扇形域问题,将径向坐标模拟成时间坐标,通过适当的变换,将扇形域问题导向哈密尔顿体系,利用分离变量法及本征函数向量展开等方法,推导出裂纹尖端的应力奇性解的计算公式,结合变分原理,提出一种解决应力奇性计算的断裂分析元,将此分析元与有限元法相结合,可以进行某些断裂力学或复合材料等应力奇性问题的计算及分析,数值计算结果表明,该方法具有精度高,使用十分方便,灵活等优点,是哈密尔顿体系和辛数学优越性的一次具体体现。     

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

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

各向异性复合材料尖劈和接头的奇性应力指数研究

提出了一个新的、基于位移的、求解三维尖劈端部奇性应力指数问题的非协调元特征分析法。该方法假定尖劈端部邻域内的位移场没有采用奇异变换技术，导出虚功方程的出发点不同于过去原有求解裂纹尖端近似场的有限元特征分析法，在有限元离散时采用的单元形式为非协调元。文中运用该方法给出了若干求解各向异性复合材料尖劈／接头端部奇性应力指数的算例。所有的计算结果表明，本文方法能够求解复杂尖劈／接头的全部奇性应力指数，使用的单元少而且精度高。     

基于哈密顿原理的两种材料界面裂纹奇性研究

研究了两种材料组成的弹性体在交界面上含裂纹时的裂纹尖端奇异场。通过变量代换及变分原理，将平面弹性扇形域的方程导向哈密体系，从而可通过分离变量及共轭辛本征函数展开法解析法求解扇形域方程，得到求解双材料界面裂纹尖点奇性的一般表达式，由此为该类问题的求解开辟了一条新途径。     

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.     

A clamped rectangular plate containing a crack

In this paper the general problem of a rectangular plate clamped along two parallel sides and containing a crack parallel to the clamps is considered. The problem is formulated in terms of a system of singular integral equations and the asymptotic behaviour of the stress state near the corners is investigated. Numerical examples are considered for a clamped plate without a crack and with a centrally located crack, and the stress intensity factors and the stresses along the clamps are calculated.     

Paradox solution on elastic wedge dissimilar materials

According to the Hellinger-Reissner variational principle and introducing proper transformation of variables , the problem on elastic wedge dissimilar materials can be led to Hamiltonian system, so the solution of the problem can be got by employing the separation of variables method and symplectic eigenfunction expansion under symplectic space, which consists of original variables and their dual variables . The eigenvalue - 1 is a special one of all symplectic eigenvalue for Hamiltonian system in polar coordinate . In general, the eigenvalue - 1 is a single eigenvalue, and the classical solution of an elastic wedge dissimilar materials subjected to a unit concentrated couple at the vertex is got directly by solving the eigenfunction vector for eigenvalue - 1. But the eigenvalue - 1 becomes a double eigenvalue when the vertex angles and modulus of the materials satisfy certain definite relationships and the classical solution for the stress distribution becomes infinite at this moment, that is, the para     

Axisymmetric Deformation of a Transversely Isotropic Cylindrical Body: A Hamiltonian State-Space Approach

A state-space approach for exact analysis of axisymmetric deformation and stress distribution in a circular cylindrical body of transversely isotropic material is developed. By means of Hamiltonian variational formulation via Legendre’s transformation, the basic equations in cylindrical coordinates are formulated into a state-space framework in which the unknown state vector comprises the displacements and associated stress components as the dual variables and the system matrix possesses the symplectic characteristics of a Hamiltonian system. Upon delineating the symplecticity of the formulation, a viable solution approach using eigenfunction expansion is developed. For illustration, an exact analysis of a finite thick-walled circular cylinder under internal and external pressures is presented, with emphasis on the end effects.     

Anti-plane shear crack normal to and terminating at the interface of two bonded piezoelectric ceramics

The electroelastic analysis of two bonded dissimilar piezoelectric ceramics with a crack perpendicular to and terminating at the interface is made. By using Fourier integral transform, the associated boundary value problem is reduced to a singular integral equation with generalized Cauchy kernel, the solution of which is given in closed form. Results are presented for a permeable crack under anti-plane shear loading and in-plane electric loading. Obtained results indicate that the electroelastic field near the crack tip in the homogeneous piezoelectric ceramic is dominated by a traditional inverse square-root singularity, while the electroelastic field near the crack tip at the interface exhibits the singularity of power law r^−α, r being distance from the interface crack tip and α depending on the material constants of a bi-piezoceramic. In particular, electric field has no singularity at the crack tip in a homogeneous solid, whereas it is singular around the interface crack tip. Numerical results are given graphically to show the effects of the material properties on the singularity order and field intensity factors.     

Symplectic solution for three dimensional couple stress problem and its variational principle

A new state vector is presented for symplectic solution to three dimensional couple stress problem. Without relying on the analogy relationship, the dual PDEs of couple stress problem are derived by a new state vector. The duality solution methodology in a new form is thus extended to three dimensional couple stress. A new symplectic orthonormality relationship is proved. The symplectic solution to couple stress theory based a new state vector is more accordant with the custom of classical elasticity and is more convenient to process boundary conditions. A Hamilton mixed energy variational principle is derived by the integral method.     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

有限尺寸下双材料界面附近裂纹位置对裂纹尖端的影响

随着复合材料的应用和发展,不同材料组成的界面结构越来越受到人们的重视。界面层两侧材料的性能相异会引起材料界面端奇异性,同时界面和界面附近存在裂纹会引起裂尖处的应力奇异性。因此双材料界面附近的力学分析是比较复杂的。本文建立双材料直角界面模型,在材料界面附近预设初始裂纹,计算了有限材料尺寸对界面应力场及其附近裂纹应力强度因子的影响。运用弹性力学中的 Goursat 公式求得直角界面端在有限尺寸下的应力场以及其应力强度系数。通过叠加原理和格林函数法进一步得到在直角界面端附近的裂纹尖端应力强度因子。计算结果表明,在适当范围内改变材料内裂纹与界面之间的距离,界面附近裂纹尖端的应力强度因子随着裂纹与界面距离的增加而减少,并且逐渐趋于稳定。分析结果可以为预测双材料结构复合材料界面失效位置提供参考。     

Fracture analysis of mode III crack problems for the piezoelectric bimorph

In this paper, a symplectic method based on the Hamiltonian system is proposed to analyze the interfacial fracture in the piezoelectric bimorph under anti-plane deformation. A set of Hamiltonian governing equations is derived from the Hamiltonian function by introducing dual variables of generalized displacements and stresses which can be expanded in series in terms of the symplectic eigensolutions. With the aid of the adjoint symplectic orthogonality, coefficients of the series are determined by the boundary conditions along the crack faces and along the external geometry. The stress\electric displacement intensity factors and energy release rates (G) directly relate to the first few terms of the nonzero eigenvalue solutions. The two ideal crack boundary conditions, namely the electrically impermeable and permeable crack assumptions, are considered. Numerical examples including the complex mixed boundary conditions are considered to show fracture behaviors of the interface crack and discuss the influencing factors.     

A new variable-order singular boundary element for calculating stress intensity factors in three-dimensional elasticity problems

A new three-dimensional variable-order singular boundary element has been constructed for stress analysis of three-dimensional interface cracks and internal material junctions. The singular fields in the vicinity of crack front or junction have been accurately represented by the singular elements by taking account the variable order of singularities and the angular profiles of field variables. Both the singular stress fields and displacement fields are independently formulated by the element’s shape functions. Different kinds of displacement formulations are investigated. The formulation combining singular and linear terms is found to be the most accurate one. The mixed-mode stress intensity factors are treated as nodal unknowns. The variation of stress intensity factors along the line of singularity can be obtained directly from the final system of equations and thus no post processing, such as three-dimensional J-integral or domain integral, is necessary. Numerical examples involving stress singularity, such as penny-shaped cracks in homogeneous and dissimilar material interface, plates with through-thickness cracks, and a dissimilar inclusion, are investigated. The analysis results are in good agreement with those reported in the literature.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.