Yunhua Li^1,2, Yunze Li³, Chieh-Li Chen⁴, Cha’o-Kuang Chen⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.1, pp. 1-14, 2010, DOI:10.3970/cmes.2010.058.001

Abstract Addressing the identification problem of the general Lagrange multiplier in the He's variational iteration method, this paper proposes a new kind of method based on Duhamel's principle for the dynamic system response analysis. In this method, we have constructed an analytical iteration formula in terms of the convolution for the residual error at the nth iteration, and have given a new interpretation to He's variational iteration method. The analysis illustrates that the computational result of this method is equal to that of He's variational iteration method on the assumption of considering the impulse response of the linear parts, or equal… More >

E. J. Sellountos¹, A. Sequeira¹, D. Polyzos²

CMES-Computer Modeling in Engineering & Sciences, Vol.57, No.2, pp. 109-136, 2010, DOI:10.3970/cmes.2010.057.109

Abstract A new Local Boundary Integral Equation (LBIE) method is proposed for the solution of plane elastostatic problems. Non-uniformly distributed points taken from a Finite Element Method (FEM) mesh cover the analyzed domain and form background cells with more than four points each. The FEM mesh determines the position of the points without imposing any connectivity requirement. The key-point of the proposed methodology is that the support domain of each point is divided into parts according to the background cells. An efficient Radial Basis Functions (RBF) interpolation scheme is exploited for the representation of displacements in each cell. Tractions in the… More >

B. Chandrasekhar¹

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.3, pp. 211-230, 2010, DOI:10.3970/cmes.2010.056.211

Abstract In this work, various basis functions used in the Method of Moments or Boundary Element (MoM/BEM) solution of acoustic scattering problems are compared with each other for their performance. Single layer formulation of the rigid bodies is considered in comparison of the solutions. Geometry of a scatterer is descritized using triangular patch modeling and basis functions are defined on triangular patches, edges and nodes for three different solutions. Far field scattering cross sections for different frequencies of incident acoustic wave are compared with the closed form solutions. Also, the errors of the solutions using these three types of basis functions… More >

Chia-Cheng Tsai¹, Chein-Shan Liu², Wei-Chung Yeih³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 131-152, 2010, DOI:10.3970/cmes.2010.056.131

Abstract The fictitious time integration method (FTIM) previously developed by Liu and Atluri (2008a) is combined with the method of fundamental solutions and the Chebyshev polynomials to solve Poisson-type nonlinear PDEs. The method of fundamental solutions with Chebyshev polynomials (MFS-CP) is an exponentially-convergent meshless numerical method which is able to solving nonhomogeneous partial differential equations if the fundamental solution and the analytical particular solutions of the considered operator are known. In this study, the MFS-CP is extended to solve Poisson-type nonlinear PDEs by using the FTIM. In the solution procedure, the FTIM is introduced to convert a Poisson-type nonlinear PDE into… More >

Sheng-Hu Ding^1,2, Xing Li², Yue-Ting Zhou^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.1, pp. 43-84, 2010, DOI:10.3970/cmes.2010.056.043

Abstract In this paper, the crack-tip fields in bonded functionally graded finite strips are studied. Different layers may have different nonhomogeneity properties in the structure. A bi-parameter exponential function was introduced to simulate the continuous variation of material properties. The problem was reduced as a system of Cauchy singular integral equations of the first kind by Laplace and Fourier integral transforms. Various internal cracks and edge crack and crack crossing the interface configurations are investigated, respectively. The asymptotic stress field near the tip of a crack crossing the interface is examined and it is shown that, unlike the corresponding stress field… More >

M. A. Kelmanson¹ and M. C. Tenwick¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 191-210, 2010, DOI:10.3970/cmes.2010.055.191

Abstract A method is presented for improving the accuracy of the widely used Gauss-Legendre Nyström method for determining approximate solutions of Fredholm integral equations of the second kind on finite intervals. The authors' recent continuous-kernel approach is generalised in order to accommodate kernels that are either singular or of limited continuous differentiability at a finite number of points within the interval of integration. This is achieved by developing a Gauss-Jacobi Nyström method that moreover includes a mean-value estimate of the truncation error of the Hermite interpolation on which the quadrature rule is based, making it particularly accurate at low orders. A… More >

A. Vukovic¹, P. Sewell¹, T. M. Benson¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 171-190, 2010, DOI:10.3970/cmes.2010.055.171

Abstract A fast and accurate method is developed for the analysis of a class of metal three-dimensional resonators with rotational symmetry. The analysis is formulated using the Body of Revolution approach and the Method of Analytical Regularization. This development is motivated by the need for three-dimensional analytical solvers that could enable fast and accurate analysis of photonic resonant structures which support very high Q whispering gallery modes and which are computationally challenging for numerical simulations. The paper outlines the formulation of the method and demonstrates the stability and the source of computation errors of the method. As a practical illustration, the… More >

J. Trevelyan¹and G. Coates¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 147-170, 2010, DOI:10.3970/cmes.2010.055.147

Abstract The terminology "wave boundary elements" relates to boundary elements enriched in the Partition of Unity sense by a multiple plane wave basis for the analysis of the propagation of short wavelength waves. This paper presents a variant of this approach in which the plane wave basis is selected adaptively according to an error indicator. The error indicator is residual based, and exhibits useful local and global properties. Model improvement in each adaptive iteration is carried out by the addition of new plane waves with no h-refinement. The convergence properties of the scheme are demonstrated. More >

Ting Du¹, J. P. Xu¹, Y.Q. Liu², Z. B. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.2, pp. 93-114, 2009, DOI:10.3970/cmes.2009.051.093

Abstract An eight-node hexahedral element with four-point quadrature in plane is developed using the assumed strain method, which can eliminate volumetric locking of incompressible material and absence of the portion of shear velocity strain related with hourglass mode to suppress hourglass mode and shear locking. In this approach, the radial return algorithm is adopted for more precise calculation of internal forces, stress and strain. In addition, a co-rotational coordinates system is established to make bending simulation much more effective, and the system is applicable to arbitrary 3D anisotropic yield criteria. A large elastic-plastic deformation of unconstrained thick plate bending example is… More >

Z.H.Yao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.1, pp. 21-46, 2009, DOI:10.3970/cmes.2009.050.021

Abstract The traditional time domain boundary integral equation (TDBIE) of elastodynamics is formulated based on the time dependent fundamental solution and the reciprocal theorem of elastodynamics. The time dependent fundamental solution of the elastodynamics is the response of the infinite elastic medium under a unit concentrate impulsive force subjected at a point and at an instant, including not only the pressure wave and shear wave, but also the Laplace wave with speed between that of P and S waves. In this paper, a new TDBIE is derived directly from the initial boundary value problem of the partial differential equation of elastodynamics,… More >