R. Q. Rodríguez^1,2, C. L. Tan², P. Sollero¹, E. L. Albuquerque³

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.5, pp. 359-372, 2014, DOI:10.3970/cmes.2014.102.359

Abstract The efficient evaluation of the fundamental solution for 3D general anisotropic elasticity is a subject of great interest in the BEM community due to its mathematical complexity. Recently, Tan, Shiah, andWang (2013) have represented the algebraically explicit form of it developed by Ting and Lee (Ting and Lee, 1997; Lee, 2003) by a computational efficient double Fourier series. The Fourier coefficients are numerically evaluated only once for a specific material and are independent of the number of field points in the BEM analysis. This work deals with the application of hierarchical matrices and low rank More >

Y.C. Shiah¹, Chung-Lei Hsu¹, Chyanbin Hwu^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.4, pp. 257-270, 2014, DOI:10.3970/cmes.2014.102.257

Abstract As has been well documented for the boundary element method (BEM), a volume integral is present in the integral equation for thermoelastic analysis. Any attempt to directly integrate the integral shall inevitably involve internal discretization that will destroy the BEM’s distinctive notion as a true boundary solution technique. Among the schemes to overcome this difficulty, the exact transformation approach is the most elegant since neither further approximation nor internal treatments are involved. Such transformation for 2D anisotropic thermoelasticity has been achieved by Shiah and Tan (1999) with the aid of domain mapping. This paper revisits More >

L. Rodríguez-Tembleque¹, M.H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.3, pp. 183-210, 2014, DOI:10.3970/cmes.2014.102.183

Abstract This work presents new contact constitutive laws for friction and wear modelling in fiber-reinforced plastics (FRP). These laws are incorporated to a numerical methodology which allows us to solve the contact problem taking into account the anisotropic tribological properties on the interfaces. This formulation uses the Boundary Element Method for computing the elastic influence coefficients. Furthermore, the formulation considers micromechanical models for FRP that also makes it possible to take into account the fiber orientation relative to the sliding direction, the fiber volume fraction, the aspect ratio of fibers, or the fiber arrangement. The proposed More >

A.I. Dimitriadis¹, N.V. Kantartzis¹ and T.D. Tsiboukis¹

CMC-Computers, Materials & Continua, Vol.39, No.3, pp. 231-265, 2014, DOI:10.3970/cmc.2014.039.231

Abstract The most important surface susceptibility models for the electromagnetic characterization of periodic metafilms, based on the dipole approximation method, are systematically analyzed in this paper. Specifically, two well-known techniques, which lead to a set of local effective surface parameters, are investigated along with a new dynamic non-local modeling algorithm. The latter formulation is properly expanded, in order to be applicable for any arbitrary periodic metafilm, irrespective of its way of excitation. The featured schemes are then directly compared toward their ability to efficiently predict the reflection and transmission properties of several lossless and lossy metafilms. More >

Y.C. Shiah¹, Y.M. Lee², T.C. Huang²

CMC-Computers, Materials & Continua, Vol.39, No.1, pp. 49-71, 2014, DOI:10.3970/cmc.2014.039.049

Abstract A common difficulty arises in characterizing the anisotropic properties of a thin sheet of anisotropic material, especially in the transverse direction. This difficulty is even more phenomenal for measuring its mechanical properties on account of its thickness. As the prelude of such investigation, this paper proposes a novel approach to identify the thermal conductivities of an unknown thin layer of anisotropic material. For this purpose, the unknown layer is sandwiched in isotropic materials with known conductivities. Prescribing proper boundary conditions, one may easily measure temperature data on a few sample boundary points. Therefore, the anisotropic More >

Y.C. Shiah¹, Y.M. Lee², Chi-Chang Wang²

CMC-Computers, Materials & Continua, Vol.33, No.3, pp. 229-255, 2013, DOI:10.3970/cmc.2013.033.229

Abstract In this paper, the boundary integrals for treating 3D field problems are fully regularized for planar elements by the technique of integration by parts (IBP). As has been well documented in open literatures, these integrals appear to be strongly singular and hyper-singular for the associated fundamental solutions. In the past, the IBP approach has only been applied to regularize the integrals for 2D problems. The present work shows that the IBP can also be further extended to treat 3D problems, where two variables of the local coordinates are involved. The presented formulations are fully explicit 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 More >

L. Rodríguez-Tembleque¹, A. Sáez¹, F.C. Buroni¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.2, pp. 131-158, 2013, DOI:10.3970/cmes.2013.096.131

Abstract A boundary element based formulation is applied to study numerically the tribological behavior of fiber-reinforced plastics (FRP) under different frictional contact conditions, taking into account the micromechanics of FRP. Micromechanical models presented consider continuous and short fiber reinforced plastics configurations. The Boundary Element Method (BEM) with an explicit approach for fundamental solutions evaluation is considered for computing the elastic influence coefficients. Signorini’s contact conditions and an orthotropic law of friction on the potential contact zone are enforced by contact operators over the augmented Lagrangian. The proposed methodology is applied to study carbon FRP under frictional More >

Yuping Zhu^1,2, Tao Shi¹, Yuanbing Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.6, pp. 501-517, 2013, DOI:10.3970/cmes.2013.095.501

Abstract Single crystal ferromagnetic shape memory alloy is a kind of new intelligent materials, it shows obvious anisotropy. Micromechanics theory has been used to analyze the whole mechanical behaviors of this material. However, Eshelby’s tensor of this material which plays an important role has still not solved efficiently. Based on the existing micromechanics constitutive model, this paper analyzes the numerical calculation formula of Eshelby’s tensor of anisotropic ferromagnetic shape memory alloy. Adopting the way of Gauss integral, the optimal Gaussian integral points for different inclusion shapes and the corresponding numerical solution of Eshelby’s tensor are obtained.Furthermore, 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… More >