M. R. Hematiyan^1,*, B. Jamshidi¹, M. Mohammadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.130, No.3, pp. 1349-1369, 2022, DOI:10.32604/cmes.2022.018235

Abstract The method of fundamental solutions (MFS) is a boundary-type and truly meshfree method, which is recognized as an efficient numerical tool for solving boundary value problems. The geometrical shape, boundary conditions, and applied loads can be easily modeled in the MFS. This capability makes the MFS particularly suitable for shape optimization, moving load, and inverse problems. However, it is observed that the standard MFS lead to inaccurate solutions for some elastostatic problems with stress concentration and/or highly anisotropic materials. In this work, by a numerical study, the important parameters, which have significant influence on the accuracy of the MFS for… More >

Yui-Chuin Shiah^{1, *}, Sheng-Chi Huang¹, M. R. Hematiyan²

CMC-Computers, Materials & Continua, Vol.64, No.2, pp. 701-727, 2020, DOI:10.32604/cmc.2020.010417

Abstract In engineering practice, analysis of interfacial thermal stresses in composites is a crucial task for assuring structural integrity when sever environmental temperature changes under operations. In this article, the directly transformed boundary integrals presented previously for treating generally anisotropic thermoelasticity in two-dimension are fully regularized by a semi-analytical approach for modeling thin multi-layers of anisotropic/isotropic composites, subjected to general thermal loads with boundary conditions prescribed. In this process, an additional difficulty, not reported in the literature, arises due to rapid fluctuation of an integrand in the directly transformed boundary integral equation. In conventional analysis, thin adhesives are usually neglected due… More >

P.D. Shah¹, Ch. Song², C.L. Tan¹, X. Wang¹

Structural Durability & Health Monitoring, Vol.2, No.4, pp. 225-238, 2006, DOI:10.3970/sdhm.2006.002.225

Abstract Numerical T-stress solutions in two dimensional anisotropic cracked bodies are very scarce in the literature. Schemes to evaluate this fracture parameter in anisotropy have been reported only fairly recently. Among them are those developed in conjunction with two different computational techniques, namely, the Boundary Element Method (BEM) and the Scaled Boundary Finite-Element Method (SBFEM). This paper provides a review of the respective schemes using these techniques and demonstrates their efficacy with three examples. These examples, which are of engineering importance, involve cracks lying in a homogeneous medium as well as at the interface between dissimilar media. The numerical T-stress solutions… More >

J. Sladek¹, V. Sladek, Ch.Zhang²

Structural Durability & Health Monitoring, Vol.1, No.2, pp. 131-144, 2005, DOI:10.3970/sdhm.2005.001.131

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for crack analysis in two-dimensional (2-d), anisotropic and linear elastic solids with continuously varying material properties. Both quasi-static and transient elastodynamic problems are considered. For time-dependent problems, the Laplace-transform technique is utilized. A unit step function is used as the test function in the local weak-form. It is leading to local boundary integral equations (LBIEs) involving only a domain-integral in the case of transient dynamic problems. The analyzed domain is divided into small subdomains with a circular shape. The moving least-squares (MLS) method is adopted for approximating the physical… More >

J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.5, pp. 477-490, 2004, DOI:10.3970/cmes.2004.006.477

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for solution of static and elastodynamic problems in a homogeneous anisotropic medium. The Heaviside step function is used as the test functions in the local weak form. It is leading to derive local boundary integral equations (LBIEs). For transient elastodynamic problems the Laplace transfor technique is applied and the LBIEs are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry such as circles in 2-d problems. The final form of local integral equations has a pure contour character only in… More >

Y.C. Shiah¹, C.L. Tan², Y.H. Chen³

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.4, pp. 243-257, 2013, DOI:10.3970/cmes.2013.096.243

Abstract The present authors have recently proposed an efficient, alternative approach to numerically evaluate the fundamental solution and its derivatives for 3D general anisotropic elasticity. It is based on a double Fourier series representation of the exact, explicit form of the Green’s function derived by Ting and Lee (1997). This paper reports on the successful implementation of the fundamental solution and its derivatives based on this Fourier series scheme in the boundary element method (BEM) for 3D general anisotropic elastostatics. Some numerical examples of stress concentration problems and a crack problem are presented to demonstrate the veracity of the implementation. The… More >

Tao Zhang^1,2, Leiting Dong^2,3, Abdullah Alotaibi⁴, Satya N. Atluri^2,5

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.1, pp. 1-28, 2013, DOI:10.3970/cmes.2013.094.001

Abstract In this paper, a novel Meshless Local Petrov-Galerkin (MLPG) Mixed Collocation Method is developed for solving the inverse Cauchy problem of linear elasticity, wherein both the tractions as well as displacements are prescribed/measured at a small portion of the boundary of an elastic body. The elastic body may be isotropic/anisotropic and simply connected or multiply-connected. In the MLPG mixed collocation method, the same meshless basis function is used to interpolate both the displacement as well as the stress fields. The nodal stresses are expressed in terms of nodal displacements by enforcing the constitutive relation between stress and the displacement gradient… More >

Y.C. Shiah¹, C. L. Tan², R.F. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.2, pp. 167-198, 2010, DOI:10.3970/cmes.2010.069.167

Abstract In this paper, fully explicit, algebraic expressions are derived for the first and second derivatives of the Green's function for the displacements in a three dimensional anisotropic, linear elastic body. These quantities are required in the direct formulation of the boundary element method (BEM) for determining the stresses at internal points in the body. To the authors' knowledge, similar quantities have never previously been presented in the literature because of their mathematical complexity. Although the BEM is a boundary solution numerical technique, solutions for the displacements and stresses at internal points are sometimes required for some engineering applications. To this… More >

Y. C. Shiah¹, C. L. Tan², V.G. Lee³

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.3, pp. 205-226, 2008, DOI:10.3970/cmes.2008.034.205

Abstract The main impediment to the development of efficient algorithms for the stress analysis of 3D generally anisotropic elastic solids using the boundary element method (BEM) and the local boundary integral equation (LBIE) meshless method over the years is the complexity of the fundamental solutions and the computational burden to evaluate them. The ability to analytically simplify and reduce them into as explicit a form as possible so that they can be directly computed will offer significant cost savings. In addition, they facilitate easy implementation using existing numerical algorithms with the above-mentioned methods that have been developed for isotropy. In this… More >

P.D. Shah¹, C.L. Tan^1,2, X. Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 185-198, 2006, DOI:10.3970/cmes.2006.013.185

Abstract In this paper, the path independent mutual or M-integral for the computation of the T-stress for interface cracks between dissimilar anisotropic, linear elastic solids, is developed. The required auxiliary field solution is derived from the solution of the problem of an anisotropic composite wedge subjected to a point force at its apex. The Boundary Element Method (BEM) is employed for the numerical stress analysis in which special crack-tip elements with the proper oscillatory traction singularity are used. The successful implementation of the procedure for evaluating the T-stress in a bi-material interface crack and its application are demonstrated by numerical examples. More >