Tienfuan Kerh^1,2, J.S. Lai¹, D. Gunaratnam², R. Saunders²

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.1, pp. 1-12, 2008, DOI:10.3970/cmes.2008.026.001

Abstract Taiwan frequently suffers from strong ground motion, and the current building code is essentially based on two seismic zones, A and B. The design value of horizontal acceleration for zone A is 0.33g, and the value for zone B is 0.23g. To check the suitability of these values, a series of actual earthquake records are considered for evaluating peak ground acceleration (PGA) for each of the zones by using neural network models. The input parameters are magnitude, epicenter distance, and focal depth for each of the checking stations, and the peak ground acceleration is calculated as the output with the… More >

P. Le¹, N. Mai-Duy², T. Tran-Cong³, G. Baker⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 43-68, 2008, DOI:10.3970/cmes.2008.025.043

Abstract This paper describes an integrated radial basis function network (IRBFN) method for the numerical modelling of the dynamics of strain localization due to strain softening in quasi-brittle materials. The IRBFN method is a truly meshless method that is based on an unstructured point collocation procedure. We introduce a new and effective regularization method to enhance the performance of the IRBFN method and alleviate the numerical oscillations associated with weak discontinuity at the elastic wave front. The dynamic response of a one dimensional bar is investigated using both local and non-local continuum models. Numerical results, which compare favourably with those obtained… More >

Shardool U. Chirputkar¹, Dong Qian²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 185-202, 2008, DOI:10.3970/cmes.2008.024.185

Abstract A multiscale method based on the extended space-time finite element method is developed for the coupled atomistic/continuum simulation of nanoscale material systems. Existing single scale approach such as the finite element method has limited capability of representing the fine scale physics in both the spatial and temporal domains. This is a major disadvantage for directly incorporating FEM in coupled atomistic/continuum simulations as it results in errors such as spurious wave reflections at the atomistic/continuum interface. While numerous efforts have been devoted to eliminating the interfacial mismatch effects, less attention has been paid to developing fine scale, atomistic level representations within… More >

G. Fitzgerald¹, G. Goldbeck-Wood², P. Kung¹, M. Petersen¹, L. Subramanian¹, J. Wescott²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 169-184, 2008, DOI:10.3970/cmes.2008.024.169

Abstract Molecular modeling has established itself as an important component of applied research in areas such as drug discovery, catalysis, and polymers. Algorithmic improvements to these methods coupled with the increasing speed of computational hardware are making it possible to perform predictive modeling on ever larger systems. Methods are now available that are capable of modeling hundreds of thousands of atoms, and the results can have a significant impact on real-world engineering problems. The article reviews some of the modeling methods currently in use; provides illustrative examples of applications to challenges in sensors, fuel cells, and nanocomposites; and finally discusses prospects… More >

Shueei-Muh Lin¹, Sen-Yung Lee², Yu-Sheng Lin³

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 175-186, 2008, DOI:10.3970/cmes.2008.023.175

Abstract The blade of a horizontal-axis wind power turbine is modeled as a rotating beam with pre-cone angles and setting angles. Based on the Bernoulli-Euler beam theory, without considering the axial extension deformation and the Coriolis forces effect, the governing differential equations for the bending vibration of the beam are derived. It is pointed out that if the geometric and the material properties of the beam are in polynomial forms, then the exact solution for the system can be obtained. Based on the frequency relations as revealed, without tedious numerical analysis, one can reach many general qualitative conclusions between the natural… More >

A.P.S. Selvadurai¹, H. Ghiabi²

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.1, pp. 53-74, 2008, DOI:10.3970/cmes.2008.023.053

Abstract This paper deals with the problem of the consolidation of a composite consisting of alternate layers of soft clay and a granular material. A series of experiments were conducted on components to develop the constitutive models that can be implemented in a computational approach. The constitutive response of the soft clay is represented by a poro-elasto-plastic Cam clay-based model and the granular medium by an elasto-plastic model with a Drucker-Prager type failure criterion and a non-associated flow rule. The computational poro-elasto-plastic model is used to calibrate the experimental results derived from the one-dimensional tests and to establish the influence of… More >

Ya Yan Lu¹, Jianxin Zhu²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 235-248, 2007, DOI:10.3970/cmes.2007.022.235

Abstract The perfectly matched layer (PML) is a widely used technique for truncating unbounded domains in numerical simulations of wave propagation problems. In this paper, the PML technique is used with a standard one-way model to solve a benchmark problem for underwater acoustics modeling. Accurate solutions are obtained with a PML layer with a thickness of only a quarter of the wavelength. The effect of a PML is analyzed in a perturbation analysis for waveguides. More >

C. Jiang¹, X. Han^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 97-118, 2007, DOI:10.3970/cmes.2007.022.097

Abstract A new uncertain optimization method is developed based on intervals and an approximation management model. A general uncertain optimization problem is considered in which the objective function and constraints are both nonlinear and uncertain, and intervals are used to model the uncertainty existing in the system. Based on a possibility degree of interval,anonlinear interval number programming (NINP) method is proposed. A deterministic objective function is constructed to maximize the possibility degree of the uncertain objective function, and the uncertain constraints are changed into deterministic ones by introducing some possibility degree levels. If the optimal possibility degree of the objective function… More >

Shamsh Tabrez, Mira Mitra, S. Gopalakrishnan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.1, pp. 77-90, 2007, DOI:10.3970/cmes.2007.022.077

Abstract In this paper, wave propagation is studied in degraded composite beam due to moisture absorption. The obtained wave responses are then used for diagnosis of the degraded zone. Moisture absorption causes an irreversible hygrothermal deterioration of the material. The change in temperature and moisture absorption changes the mechanical properties. Thus this affects the structure in dimensional stability as well as material degradation due to reduction in mechanical properties. Here, the composite beam is modeled as Timoshenko beam using wavelet based spectral finite element (WSFE) method. The WSFE technique is especially tailored for simulation of wave propagation. It involves Daubechies scaling… More >

L.N. Gergidis, D. Kourounis, S. Mavratzas, A. Charalambopoulos¹

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.2, pp. 157-176, 2007, DOI:10.3970/cmes.2007.021.157

Abstract A new complete set of scattering eigensolutions of Helmholtz equation in spheroidal geometry is constructed in this paper. It is based on the extension to exterior boundary value problems of the well known Vekua transformation pair, which connects the kernels of Laplace and Helmholtz operators. The derivation of this set is purely analytic. It avoids the implication of the spheroidal wave functions along with their accompanying numerical deficiencies. Using this novel set of eigensolutions, we solve the acoustic scattering problem from a soft acoustic spheroidal scatterer, by expanding the scattered field in terms of it. Two approaches concerning the determination… More >