Y.C. Shiah¹, C. L. Tan², C.Y. Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 1-22, 2012, DOI:10.3970/cmes.2012.087.001

Abstract The fundamental solution, or Green's function, for 3D anisotropic elastostatics as derived by Ting and Lee (1997) [Q.J. Mech. Appl. Math.; 50: 407-426] is one that is fully explicit and algebraic in form. It has, however, only been utilized in boundary element method (BEM) formulations quite recently even though it is relatively straightforward and direct to implement. This Green's function and its derivatives are necessary items in this numerical analysis technique. By virtue of the periodic nature of the angles when it is expressed in the spherical coordinate system, the present authors have very recently represented the Green's function as… More >

Y.C. Shiah¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 95-108, 2011, DOI:10.3970/cmes.2011.078.095

Abstract By differentiating the Green function of Ting and Lee (1997) for 3D general anisotropic elastotatics in a spherical coordinate system as an intermediate step, and then using the chain rule, derivatives of up to the second order of this fundamental solution are obtained in exact, explicit, algebraic forms. No tensors of order higher than two are present in these derivatives, thereby allowing these quantities to be numerically evaluated quite expeditiously. These derivatives are required for the computation of the internal point displacements and stresses via Somigliana's identity in BEM analysis. Some examples are presented to demonstrate their successful implementation to… More >

K.P. Baxevanakis¹, A.E. Giannakopoulos²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.2, pp. 181-198, 2010, DOI:10.3970/cmes.2010.060.181

Abstract The present work gives a systematic and rigorous implementation of (edge type) circular Volterra dislocation loops in ordinary axisymmetric finite elements using the thermal analogue and the integral representation of dislocations through stresses. The accuracy of the proposed method is studied in problems where analytical solutions exist. The full fields are given for loop dislocations in isotropic and anisotropic crystals and the Peach-Koehler forces are calculated for loops approaching free surfaces and bimaterial interfaces. The results are expected to be very important in the analysis of plastic yield strength, giving quantitative results regarding the influence of grain boundaries, interstitial particles,… More >

C. L. Tan¹, Y.C. Shiah², C.W. Lin²

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.3, pp. 195-214, 2009, DOI:10.3970/cmes.2009.041.195

Abstract The explicit, closed-form expressions of the Green's functions for generally anisotropic elastic solids in three-dimensions that have been derived using Stroh's formalism are employed in a formulation of the boundary element method (BEM). Unlike several other existing schemes, the evaluation of these fundamental solutions does not require further numerical integration in the BEM algorithm; they have surprisingly not been implemented previously. Three numerical examples are presented to demonstrate the veracity of the implementation and the general applicability of the BEM for the 3D elastic stress analysis of generally anisotropic solids. The results are compared with known solutions in the literature… More >

B. Yang², V. K. Tewary³

CMC-Computers, Materials & Continua, Vol.8, No.1, pp. 23-32, 2008, DOI:10.3970/cmc.2008.008.023

Abstract A three-dimensional Green's function for a material system consisting of anisotropic and linearly elastic planar multilayers with interfacial membrane and flexural rigidities has been derived. The Stroh formalism and two-dimensional Fourier transforms are applied to derive the general solution for each homogeneous layer. The Green's function for the multilayers is then solved by imposing the surface boundary condition, the interfacial displacement continuity condition, and the interfacial traction discontinuity condition. The last condition is given by the membrane and bending equilibrium equations of the interphases modeled as Kirchhoff plates. Numerical results that demonstrate the validity and efficiency of the formulation are… More >