BEM for Static and Dynamic Fracture Analysis in Thin Piezoelectric Structures

A 2D time-domain boundary element method (BEM) is developed to study the fracture problems in thin piezoelectric structure. The nearly singular integrals arisen in thin structures are calculated in two ways. One is based on a nonlinear coordinate transformation for curve-quadratic element, and the other one is an analytical integration method for straight quadratic element. Numerical examples are presented to verify the effectiveness and stability of the present BEM in thin piezoelectric structure.

A Two-Level Preconditioned Conjugate-Gradient Method in Distorted and Structured Grids

In this paper, we propose a new two-level preconditioned C-G method which uses the quadratic smoothing and the linear correction in distorted but topo-logically structured grid. The CPU time of this method is less than that of the multigrid preconditioned C-G method (MGCG) using the quadratic element, but their accuracy is almost the same. Numerical experiments and eigenvalue analysis are given and the results show that the proposed two-level preconditioned method is efficient.

A NONSINGULAR BOUNDARY ELEMENT METHOD FOR THE TORSION PROBLEM OF THE ANISOTROPIC UNIFORM BAR

The presentation is mainly devoted to the research on the regularized boundary integral equations (BIEs) with indirect unknowns for torsion problem of the anisotropic uniform bar. Based on a new view and idea, a novel regularization technique is pursued, in which the nonsingular indirect BIE (IBIE) excluding the CPV and HFP integrals is established. Such torsion problems can be solved directly by using the presented technique without transforming them into isotropic ones, for this reason, no inverse transform is required. Moreover, a unique feature of the shear stress BIEs expressed by density functions is that they are independent of the warp BIEs and, as such, can be collocated at the same locations as the warp BIEs. This provides additional and concurrently useable equations for various purposes. Besides, in the numerical implementation, the boundary geometric is depicted by exact elements, while the distribution of the boundary quantity on each element is approximated by a discontinuous quadratic element. Some numerical examples will be applied to validate the current scheme. It is shown that a better precision and high-computational efficiency can be achieved by the presentation.