A contour integral method is developed for computation of stress intensity and electric intensity factors for cracks in continuously nonhomogeneous piezoelectric body under a transient dynamic load. It is shown that the asymptotic fields in the crack-tip vicinity in a continuously nonhomogeneos medium is the same as in a homogeneous one. A meshless method based on the local Petrov-Galerkin approach is applied for computation of physical fields occurring in the contour integral expressions of intensity factors. A unit step function is used as the test functions in the local weakform. This leads to local integral equations (LBIEs) involving only contour-integrals on the surfaces of subdomains. The moving least-squares (MLS) method is adopted for approximating the physical quantities in the LBIEs. The accuracy of the present method for computing the stress intensity factors (SIF) and electrical displacement intensity factors (EDIF) are discussed by comparison with available analytical or numerical solutions.

A meshless method based on the local Petrov-Galerkin approach is proposed for crack analysis in two-dimensional (2-D) and three-dimensional (3-D) axisymmetric magneto-electric-elastic solids with continuously varying material properties. Axial symmetry of geometry and boundary conditions reduces the original 3-D boundary value problem into a 2-D problem in axial cross section. Stationary and transient dynamic problems are considered in this paper. The local weak formulation is employed on circular subdomains where surrounding nodes randomly spread over the analyzed domain. The test functions are taken as unit step functions in derivation of the local integral equations (LIEs). The moving least-squares (MLS) method is adopted for the approximation of the physical quantities in the LIEs. The accuracy of the present method for computing the stress intensity factors (SIF), electrical displacement intensity factors (EDIF) and magnetic induction intensity factors (MIIF) are discussed by comparison with numerical solutions for homogeneous materials.

Keywords:Meshless local Petrov-Galerkin method (MLPG) Moving least-squares approximation FGM 2D problems Internal and edge cracks Houbolt method a b s t r a c t A meshless method based on the local Petrov-Galerkin approach is proposed, to solve initial-boundary value problems of magneto-electro-elastic solids with continuously varying material properties. Stationary and transient thermal problems are considered in this paper. The mechanical 2-D fields are described by the equations of motion with an inertial term. Nodal points are spread on the analyzed domain, and each node is surrounded by a small circle for simplicity. The spatial variation of displacements, electric and magnetic potentials is approximated by the moving least-squares (MLS) scheme. After performing the spatial integrations, one obtains a system of ordinary differential equations for certain nodal unknowns. That system is solved numerically by the Houbolt finite-difference scheme as a time stepping method.

a b s t r a c tThe meshless local Petrov-Galerkin (MLPG) method is used to analyze transient dynamic problems in 3D axisymmetric piezoelectric solids with continuously inhomogeneous material properties. Both mechanical and thermal loads are considered here. A 3D axisymmetric body is created by rotation of a cross section around an axis of symmetry. Axial symmetry of geometry and boundary conditions reduces the original 3D boundary value problem into a 2D problem. The cross section is covered by small circular sub-domains surrounding nodes randomly spread over the analyzed domain. A unit step function is chosen as test function, in order to derive local integral equations on the boundaries of the chosen sub-domains, called local boundary integral equations (LBIE). These integral formulations are either based on the Laplace transform technique or the time-difference approach. The local integral equations are non-singular and take a very simple form, despite of inhomogeneous and anisotropic material behaviour across the analyzed structure. Spatial variation of all physical fields (or of their Laplace transforms) at discrete time instants are approximated on the local boundary and in the interior of the sub-domain by means of the moving least-squares (MLS) method. The Stehfest algorithm is applied for the numerical Laplace inversion, in order to retrieve the time-dependent solutions.

No abstract