Wei-Ming Lee¹, Jeng-Tzong Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 141-160, 2009, DOI:10.3970/cmes.2009.050.141

Abstract This paper presents the multipole Trefftz method to derive an analytical model describing the free vibration of a circular plate with multiple circular holes. Based on the addition theorem, the solution of multipoles centered at each circle can be expressed in terms of multipoles centered at one circle, where boundary conditions are specified. In this way, a coupled infinite system of simultaneous linear algebraic equations is derived for the circular plate with multiple holes. The direct searching approach is employed in the truncated finite system to determine the natural frequencies by using the singular value decomposition (SVD). After determining the… More >

D. Soares Jr.¹, J. Sladek², V. Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 115-140, 2009, DOI:10.3970/cmes.2009.050.115

Abstract In recent years the idea of using mesh-free or mesh-less methods for numerical solution of partial differential equations has received much attention, due to their potential advantage in eliminating the costly effort of mesh generation and re-meshing. A variety of meshless methods has been proposed so far. Many of them are derived from a weak-form formulation on global domain or a set of local subdomains. In the global formulation background cells are required for the integration of the weak form. In methods based on local weak-form formulation no cells are required and therefore they are often referred to as truly… More >

Delfim Soares Jr. ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 97-114, 2009, DOI:10.3970/cmes.2009.050.097

Abstract In this work, meshless methods based on the local Petrov-Galerkin (MLPG) approach are presented to analyse electromagnetic wave propagation problems. Formulations adopting the Heaviside step function and the Gaussian weight function as the test functions in the local weak form are considered. The moving least square (MLS) method is used to approximate the physical quantities in the local integral equations. After spatial discretization is carried out, a system of ordinary differential equations of second order is obtained. This system is solved in the time-domain by the Houbolt's method, allowing the computation of the so-called primary fields (either the electric or… More >

Sudib K. Mishra¹, J. K. Paik², S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.1, pp. 67-96, 2009, DOI:10.3970/cmes.2009.050.067

Abstract We present a simulation based study, by combining several models involving multiple time scales and physical processes, which govern the interfacial stress corrosion cracking (SCC) in grain boundaries, layered composites or bi-materials, and the mechanisms of inhibition using `smart' agents. The inhibiting agents described herein, automatically sense the initiation of damage, migrate to the sites and delay the corrosion kinetics involved in the process. The phenomenon of SCC is simulated using the lattice spring model (for the mechanical stresses), coupled with a finite difference model of diffusing species, causing the dissolution of the interfacial bonds. The dissolution is expressed through… More >

Tai-Ming Chang¹, Chien-Chou Weng¹, Mei-Jiau Huang^1,2, Chun-KaiLiu², Chih-Kuang Yu²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.1, pp. 47-66, 2009, DOI:10.3970/cmes.2009.050.047

Abstract We employ the non-equilibrium molecular dynamics (NEMD) simulation to calculate the in-plane thermal conductivity of silicon thin films of thickness 2.2nm and 11nm. To eliminate the finite-size effect, samples of various lengths are simulated and an extrapolation technique is applied. To perform the quantum correction which is necessary as the MD simulation temperature is lower than Debye temperature, the confined phonon spectra are obtained in advance via the EMD simulations. The investigation shows the thermal conductivities corrected based on the bulk and thin-film phonon densities of states are very close and they agree excellently with the theoretical predictions of a… 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 >

Yaoming Zhang¹, Yan Gu¹, Jeng-Tzong Chen²

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

Abstract There exist nearly singular integrals for thin walled structures in the boundary element method (BEM). In this paper, an efficient analytical method is developed to deal with the nearly singular integrals in the boundary integral equations (BIEs) for 2-D thin walled structures. The developed method is possible for problems defined in high-order geometry elements when the nearly singular integrals need to be calculated. For the analysis of nearly singular integrals with high-order geometry elements, much fewer boundary elements can be used to achieve higher accuracy. More importantly, computational models of thin walled structures or thin shapes in structures demand a… More >

S. Tangaramvong, F. Tin-Loi

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.3, pp. 269-296, 2009, DOI:10.3970/cmes.2009.049.269

Abstract We consider the optimal synthesis, namely minimum weight design, of rigid perfectly-plastic structures for which some supports involve unilateral frictional contact. This problem is of interest as it is not only encountered in practice but it also involves, in the general friction case, a nonassociative complementarity condition that makes it theoretically and numerically challenging. For simplicity of exposition, we focus on the class of bar structures for which yielding is governed by either pure bending or by combined axial and flexural forces. In view of possible multiplicity of solutions due to nonassociativity, a direct optimization formulation may lead to an… More >

Haiping Ye¹, Yongsheng Ding²

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.3, pp. 255-268, 2009, DOI:10.3970/cmes.2009.049.255

Abstract A fractional-order model for the immunological and therapeutic control of HIV is studied qualitatively. The equilibria are found and their local stability are investigated. Also the global stability of the infection-free equilibrium is established. The optimal efficacy level of anti-retroviral therapy needed to eradicate HIV from the body of an HIV-infected individual is obtained. More >

Gerardo Anguiano-Orozco^1,2, Rubén Avila³

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.3, pp. 217-254, 2009, DOI:10.3970/cmes.2009.049.217

Abstract A numerical investigation of the transient, three dimensional, laminar natural convection of a fluid confined in a spherical container is carried out. Initially the fluid is quiescent with a uniform temperature T_i equal to the temperature of the wall of the container. At time t=0, the temperature of the wall is suddenly lowered to a uniform temperature T_w=0. The natural convection, that conducts to a vortex ring formation within the sphere, is driven by a terrestrial gravity force (laboratory gravity) and by the step change in the temperature of the wall. A scaling analysis of a simplified transient, two dimensional… More >