S.D. Akbarov^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.3, pp. 173-205, 2013, DOI:10.3970/cmes.2013.095.173

Abstract By employing the Three-dimensional Linearized Theory of Elastic Waves in Initially Stressed Bodies (TLTEWISB) the time-harmonic Lamb’s problem for a system comprising a finite pre-strained half-space and finite pre-strained covering layer made of incompressible materials is examined for the case where the material of the covering layer is stiffer than that of the half-space material. It is assumed that on the upper free face plane of the covering layer the point-located time-harmonic force acts. The elasticity relations of the materials are described through Treloar’s potential. The corresponding boundary-value problem is solved by employing the Hankel integral transformation. The corresponding inverse… More >

Zichun Yang^1,2,3, Kunfeng Li^1,4, Qi Cai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.2, pp. 143-171, 2013, DOI:10.3970/cmes.2013.095.143

Abstract The conventional probabilistic reliability model for structures is based on the “probability assumption” and “binary-state assumption”. These assumptions are often offset the reality of practical engineering and lead to a wrong conclusion. In fact, besides randomness, fuzziness which is different from randomness in nature is also a prevalent uncertainty factor and plays an important role in structural reliability assessment. In this paper, a novel structural reliability model with random variables and fuzzy variables is established by using the fuzzy set theory, possibility theory and probability measure for fuzzy events, based on the “mixed probability and possibility assumption” and “fuzzy state… More >

Sanderson L. Gonzaga de Oliveira¹, Mauricio Kischinhevsky², João Manuel R. S. Tavares³

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.2, pp. 119-141, 2013, DOI:10.3970/cmes.2013.095.119

Abstract A novel graph-based adaptive mesh refinement technique for triangular finite-volume discretizations in order to solve second-order partial differential equations is described. Adaptive refined meshes are built in order to solve timedependent problems aiming low computational costs. In the approach proposed, flexibility to link and traverse nodes among neighbors in different levels of refinement is admitted; and volumes are refined using an approach that allows straightforward and strictly local update of the data structure. In addition, linear equation system solvers based on the minimization of functionals can be easily used; specifically, the Conjugate Gradient Method. Numerical and analytical tests were carried… More >

R.A. Amaro Junior¹, L.Y. Cheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.2, pp. 87-118, 2013, DOI:10.3970/cmes.2013.095.087

Abstract In this paper simulations of brittle fracture and hydroelastic problems are carried out by using a numerical approach based on the Moving Particle Simulation (MPS) method. It is a meshless method used to model both fluid and elastic solid, and all the computational domain is discretized in Lagrangian particles. A higher order accuracy gradient operator is used herein by adopting a correction matrix. Also, in order to correctly simulate the collision of the fragments, a contact detection algorithm that takes into account the presence of the solid surfaces generated by brittle fracture is proposed. In case of fluid-structure interaction, a… More >

F. Wang¹, J.L. Li¹, B.X. Yang¹, N.A. Hill²

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.1, pp. 59-85, 2013, DOI:10.3970/cmes.2013.095.059

Abstract A gas-liquid two-phase model for the simulation of a power-law fluid mold filling process with the consideration of phase change is proposed, in which the governing equations for the melt and air in the cavity, including the mass conservation, momentum conservation and energy conservation equations, are unified into one system of equation. A revised Enthalpy method, which can be used for both the melt and air in the mold cavity, is proposed to describe the phase change during the mold filling. Finite volume method on non-staggered grid is used to solve the system. The level set method is used to… More >

K. Luo¹, Y. L. Pi¹, W. Gao¹, M. A. Bradford¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.1, pp. 31-58, 2013, DOI:10.3970/cmes.2013.095.031

Abstract This paper presents a theoretical analysis for the time-dependent nonlinear behaviour and buckling of shallow concrete-filled steel tubular (CFST) arches under a sustained central concentrated load. The virtual work method is used to establish the differential equations of equilibrium for the time-dependent behaviour and buckling analyses of shallow CFST arches, and the age-adjusted effective modulus method is adopted to model the creep behaviour of the concrete core. Analytical solutions of time-dependent displacements and internal forces of shallow CFST arches are derived. It has been found that under a sustained central concentrated load, the deformations and bending moments in a shallow… More >

Seyed Mahmoud Hosseini ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.1, pp. 1-29, 2013, DOI:10.3970/cmes.2013.095.001

Abstract In this article, the effects of carbon nanotubes distributions on natural frequency are studied for a functionally graded nanocomposite thick hollow cylinder reinforced by single-walled carbon nanotubes using a hybrid mesh-free method. The FG nanocomposite cylinder is excited by a shock loading, which is applied on the inner surface of cylinder. The first natural frequency is obtained for various nonlinear grading patterns of distributions of the aligned carbon nanotubes. The effects of various nonlinear grading patterns on natural frequency are obtained and discussed in details. The presented hybrid mesh-free method is based on the generalized finite difference (GFD) method for… More >

H. Yoshikawa¹, R. Matsuura², N. Nishimura¹

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.6, pp. 553-565, 2013, DOI:10.3970/cmes.2013.094.553

Abstract The convolution integrals with respect to time in the time domain boundary integral equation method (TD-BIEM) are calculated approximately using the Lubich convolution quadrature method (CQM). The influence matrices in the discretized boundary integral equation are computed with the Laplace transform of the fundamental solution in TD-BIEM with the Lubich CQM. These matrices, however, are dense, and both the computational cost and memory requirements are high. In this paper, we apply Adaptive Cross Approximation (ACA) to the influence matrices to achieve a fast solver of TD-BIEM with the Lubich CQM. Moreover, we reduce the computational time of TD-BIEM with the… More >

Hongtao Wang^1,2, Haitao Wang², Lie Jin², Zhenhan Yao³

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.6, pp. 529-552, 2013, DOI:10.3970/cmes.2013.094.529

Abstract Since only the boundary of the analyzed domain needs to be discretized, the boundary element method (BEM) inherently has the advantages of solving crack problems. In this paper, a micromechanics BEM scheme is applied to determine the effective elastic moduli of three-dimensional (3-D) solids containing a large number of parallel or randomly oriented microcracks. The 3-D analyses accelerated by the fast multipole method were carried out to investigate the relations between the effective elastic moduli and the microcrack density parameter. Numerical examples show that the results agree well with the available analytical solution and known micromechanics models. From the numerical… More >

Haitao Wang^1,2, Zhenhan Yao³

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.6, pp. 507-527, 2013, DOI:10.3970/cmes.2013.094.507

Abstract This paper presents the use of a three-dimensional time domain boundary element method (BEM) in conjunction with adaptive cross approximation method (ACA) for dynamic analyses of the HTR-PM nuclear island foundation. The advantage of this approach is that only foundation of the HTR-PM nuclear island and limited surfaces of the supporting half-space soil medium are modeled and analyzed in a direct time stepping scheme. In addition, the ACA can compress the BEM coefficient matrices at each time step efficiently, therefore allowing larger models to be analyzed compared with conventional BEMs. In order to discretize the boundary integral equation (BIE) we… More >