G. Q. Xie^1,2, S. Y. Long^1,3

CMC-Computers, Materials & Continua, Vol.4, No.1, pp. 11-20, 2006, DOI:10.3970/cmc.2006.004.011

Abstract The concept of the micropolar theory is employed to investigate vibration behaviors of carbon nanotubes. The constitutive relation has been deduced from the two-dimensional analysis of the microstructure of the carbon nanotube. Van der Waals interactions are simulated by a weak spring model. Hamilton's principle is employed to obtain dynamics equations of the multi-walled carbon nanotube. Numerical examples for both single-walled and double-walled carbon nanotubes are presented and the significant difference in vibration behaviors between them has been distinguished. Numerical results show that fundamental frequencies for the cantilever single-walled carbon nanotube decreases with increase of More >

D.L. Young^1,2, C.C. Tsai³, Y.C. Lin¹, C.S. Chen⁴

CMC-Computers, Materials & Continua, Vol.4, No.1, pp. 1-10, 2006, DOI:10.3970/cmc.2006.004.001

Abstract This paper describes the method of fundamental solutions (MFS) to solve eigenfrequencies of plate vibrations by utilizing the direct determinant search method. The complex-valued kernels are used in the MFS in order to avoid the spurious eigenvalues. The benchmark problems of a circular plate with clamped, simply supported and free boundary conditions are studied analytically as well as numerically using the discrete and continuous versions of the MFS schemes to demonstrate the major results of the present paper. Namely only true eigenvalues are contained and no spurious eigenvalues are included in the range of direct More >

S. Straface¹, S. Troisi, V. Gagliardi

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 155-166, 2006, DOI:10.3970/cmc.2006.003.155

Abstract The present work shows an alternative to the classical methods to solve the Richards' Equation (RE), used to model flow in unsaturated porous media. This alternative is named Cell Method (CM). The CM is based on a preliminary reformulation of the mathematical model in a partially discrete form, which preserves as much as possible the physical and geometrical content of the original problem, and is made possible by the existence and properties of a common mathematical structure of field theories. The goal is to maintain the focus, both in the modelling and discretization steps, on… More >

T. Nishioka¹, S. Tchouikov¹, T. Fujimoto¹

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 147-154, 2006, DOI:10.3970/cmc.2006.003.147

Abstract In this study, phenomena of multiple branching of dynamically propagating crack are investigated numerically. The complicated paths of cracks propagating in a material are simulated by moving finite element method based on Delaunay automatic triangulation (MFEM BODAT), which was extended for such problems. For evaluation of fracture parameters for propagating and branching cracks switching method of the path independent dynamic J integral was used. Using these techniques the generation phase simulation of multiple dynamic crack branching was performed. Various dynamic fracture parameters, which are almost impossible to obtain by experimental technique alone, were accurately evaluated. More >

G. Haasemann¹, M. Kästner¹ and V. Ulbricht¹

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 131-146, 2006, DOI:10.3970/cmc.2006.003.131

Abstract Novel textile reinforced composites provide an extremely high adaptability and allow for the development of materials whose features can be adjusted precisely to certain applications. A successful structural and material design process requires an integrated simulation of the material behavior, the estimation of the effective properties which need to be assigned to the macroscopic model and the resulting features of the component. In this context two efficient modelling strategies - the Binary Model (Carter, Cox, and Fleck (1994)) and the Extended Finite Element Method (X-FEM) (Moës, Cloirec, Cartraud, and Remacle (2003)) - are used to… More >

X.G. Yuan¹, R.J. Zhang²

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 119-130, 2006, DOI:10.3970/cmc.2006.003.119

Abstract In this paper, the problem of cavity formation and motion in an incompressible transversely isotropic nonlinearly elastic solid sphere, which is subjected to a uniform radial tensile dead load on its surface, is examined in the context of nonlinear elastodynamics. The strain energy density associated with the nonlinearly elastic material may be viewed as the generalized forms of some known material models. It is proved that some determinate conditions must be imposed on the form of the strain energy density such that the surface tensile dead load has a finite critical value. Correspondingly, as the More >

B.B. Karki¹, D. Bhattarai¹, L. Stixrude²

CMC-Computers, Materials & Continua, Vol.3, No.3, pp. 107-118, 2006, DOI:10.3970/cmc.2006.003.107

Abstract Computer modeling of liquid phase poses tremendous challenge: It requires a relatively large simulation size, long simulation time and accurate interatomic interaction and as such, it produces massive amounts of data. Recent advances in hardware and software have made it possible to accurately simulate the liquid phase. This paper reports the details of methodology used in the context of liquid simulations and subsequent analysis of the output data. For illustration purpose, we consider the results for the liquid phases of two geophysically relevant materials, namely MgO and MgSiO₃. The simulations are performed using the parallel first-principles More >

Kunigal N. Shivakumar¹, Anil Bhargava¹, Sameer Hamoush²

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 97-106, 2006, DOI:10.3970/cmc.2006.003.097

Abstract A detailed three-dimensional finite element stress analysis was conducted on straight-shank and countersunk rivet holes in a plate subjected to tension loading. The study included a wide range of plate width to radius, thickness to radius, countersunk depth to thickness ratios and countersunk angles(θ_c). The stress concentration is maximum at or near the countersunk edge. The stress concentration depends on countersunk depth, plate thickness and width and it is nearly independent of the countersunk angle for 80° ≤ θ_c ≤ 120°. Using the finite element results and limiting conditions, an equation for stress concentration factor is More >

Hamid Sharifi¹ and Augustin Gakwaya¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 77-96, 2006, DOI:10.3970/cmc.2006.003.077

Abstract In this paper, an object-oriented modeling of solid material constitutive behavior using the UML notation is presented. Material properties are first classified into large and small deformation kinematical models. In the small deformation package, we keep classes such as Elastic, ElastoPlastic, ViscoElastic and ViscoPlastic. In the large deformation package, we store classes such as ElastoPlastic, HyperElastic, HyperPlastic, HyperViscoElastic, HyperViscoPlastic and so on. The hierarchical structure, the association relationships as well as key attributes and methods of these classes are presented. We used a C++ implementation of the above model for developing HyperElastic, HyperElastoPlastic and Contact More >

P. B. Wang¹, Z. H. Yao^1,2, T. Lei¹

CMC-Computers, Materials & Continua, Vol.3, No.2, pp. 65-76, 2006, DOI:10.3970/cmc.2006.003.065

Abstract The fast multipole method (FMM) is applied to the dual boundary element method (DBEM) for the analysis of finite solids with large numbers of microcracks. The application of FMM significantly enhances the run-time and memory storage efficiency. Combining multipole expansions with local expansions, computational complexity and memory requirement are both reduced to O(N), where N is the number of DOFs (degrees of freedom). This numerical scheme is used to compute the effective in-plane bulk modulus of 2D solids with thousands of randomly distributed microcracks. The results prove that the IDD method, the differential method, and the More >