Ching Kong Chao^1,2, Chin Kun Chen¹, Fu Mo Chen³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.3, pp. 283-298, 2009, DOI:10.3970/cmes.2009.047.283

Abstract Analytical exact solutions of a fundamental heat conduction problem for a three-phase elliptical composite under a remote uniform heat flow are provided in this paper. The steady-state temperature and heat flux fields in each phase of an elliptical composite are analyzed in detail. Investigations on the present heat conduction problem are tedious due to the presence of material inhomogeneities and geometric discontinuities. Based on the technique of conformal mapping and the method of analytical continuation in conjunction with the alternating technique, the general expressions of the temperature and heat flux are derived explicitly in a More >

B.R.Kim¹ and H.K.Lee^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.3, pp. 253-282, 2009, DOI:10.3970/cmes.2009.047.253

Abstract Damage characteristic of particulated ductile composites is a complex evolutionary phenomenon that includes particle debonding and void evolution with the accumulation of the plastic straining of the ductile matrix. In this paper, a micromechanical elastoplastic damage model for ductile matrix composites considering gradually incremental damage (particle debonding and void evolution) is proposed to predict the overall elastoplastic behavior and damage evolution in the composites. The constitutive damage model proposed in an earlier work by the authors [Kim and Lee (2009)] considering particle debonding is extended to accommodate the gradually incremental damage and elastoplastic behavior of More >

Gang Lin¹, Yinghua Liu^1,2, Zhihai Xiang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.3, pp. 217-252, 2009, DOI:10.3970/cmes.2009.047.217

Abstract Deterioration of reinforced concrete (RC) structures due to chloride ingress followed by reinforcement corrosion is a serious problem all over the world, therefore prediction of chloride profiles is a key element in evaluating durability and integrity of RC structures exposed to chloride environment. In the present paper, an integrated finite element-based computational framework is developed for predicting service life of RC structures exposed to chloride environment, which takes environment temperature and humidity fluctuations, diffusion and convection, chloride binding, as well as the decay of durability of structures caused by coupled deterioration processes into account. The… More >

G. Kosec¹, B. Šarler²

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 191-216, 2009, DOI:10.3970/cmes.2009.047.191

Abstract This paper explores an application of a novel mesh-free Local Radial Basis Function Collocation Method (LRBFCM) [Sarler and Vertnik (2006)] in solution of coupled heat transfer and fluid flow problems with solid-liquid phase change. The melting/freezing of a pure substance is solved in primitive variables on a fixed grid with convection suppression, proportional to the amount of the solid fraction. The involved temperature, velocity and pressure fields are represented on overlapping sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective… More >

Chih-Chang Chi¹, Weichung Yeih^1,2, Chein-Shan Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 167-190, 2009, DOI:10.3970/cmes.2009.047.167

Abstract In this study, a novel method for solving the Cauchy problem of Laplace equation is developed. Through the fictitious time integration method (FTIM), the finding of the root of the resulting linear equations can be transformed into for finding the fixed point of a system of first order ordinary differential equations, in which a fictitious time variable is introduced. In such a sense, the inverse of ill-posed leading matrix is not necessary for the FTIM. This method uses the residual of each equation to control the evolution of unknowns in the fictitious time, and it More >

M. De Munck¹, D. Moens², W. Desmet³, D.Vandepitte³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 119-166, 2009, DOI:10.3970/cmes.2009.047.119

Abstract Deterministic simulation tools enable a very precise simulation of physical phenomena using numerical models. In many real life situations however, a deterministic analysis is not sufficient to assess the quality of a design. In a design stage, some physical properties of the model may not be determined yet. But even in a design ready for production, design tolerances and production inaccuracies introduce variability and uncertainty. In these cases, a non-deterministic analysis procedure is required, either using a probabilistic or a non-probabilistic approach. The authors developed an intelligent Kriging response surface based optimisation procedure that can More >

T.B Jönsthövel¹, M.B. van Gijzen², C.Vuik², C. Kasbergen¹, A. Scarpas¹

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 97-118, 2009, DOI:10.3970/cmes.2009.047.097

Abstract The introduction of computed x-ray tomography allows for the construction of high quality, material-per-element based 3D meshes in the field of structural mechanics. The use of these meshes enables a shift from meso to micro scale analysis of composite materials like cement concrete, rocks and asphalt concrete. Unfortunately, because of the extremely long execution time, memory and storage space demands, the majority of commercially available finite element packages are not capable of handling efficiently the most computationally demanding operation of the finite element solution process, that is, the inversion of the structural stiffness matrix. To More >

J. Sladek¹, V. Sladek¹, P. Solek², C.L. Tan³, Ch. Zhang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 61-96, 2009, DOI:10.3970/cmes.2009.047.061

Abstract The meshless local Petrov-Galerkin (MLPG) method for transient linear thermoelastic analysis is presented. Orthotropic material properties are considered here. In uncoupled thermoelasticity, the temperature field is not influenced by displacements. Therefore, in the first step, the heat conduction equation is solved for the temperature distribution in the domain. The equations of motion are then solved with the inertial term considered. A Heaviside step function as the test functions is applied in the weak-form to derive local integral equations for solving two- and three-dimensional problems. Local integral equations are written on small sub-domains with circular or More >

J. Noorzaei¹, K.H. Bayagoob², A.A. Abdulrazeg¹, M.S. Jaafar^1,1, T.A. Mohammed¹

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 43-60, 2009, DOI:10.3970/cmes.2009.047.043

Abstract This paper focuses on the development, verification and application of a three-dimensional finite element code for coupled thermal and structural analysis of roller compacted concrete dams. The Kinta RCC gravity dam, which is the first roller compacted concrete dam in Malaysia, has been taken for the purpose of verification of the finite element code. The actual climatic conditions and thermal properties of the materials were considered in the analysis. The structural stress analysis was performed using the elasto-plastic stress analysis. The Mohr yield criterion which is widely used for concrete plasticity modeling was adopted in More >

A. Djeukou¹, O. von Estorff²

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 23-42, 2009, DOI:10.3970/cmes.2009.047.023

Abstract The paper deals with the response of rectangular composite plates to low-velocity impact. A third-order shear deformation theory as well as the Newmark integration are used to determine the contact force history analytically. The interaction between the impactor and the plate is modeled with the help of a two degrees-of-freedom system, consisting of springs and masses. The Choi's linearized Hertzian contact model is used to determine the contact force. The maximum impact force is employed for a static damage analysis of the composite plate by means of the radial point interpolation method, while the Tsai-Wu More >