F.C. Araújo¹, L.J. Gray²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 103-122, 2008, DOI:10.3970/cmes.2008.024.103

Abstract In recent years, carbon nanotubes (CNTs) have been widely employed to build advanced composites. In this work, a Boundary Element Method (BEM) is applied to 3D representative volume elements (RVEs) to estimate mechanical properties of CNT-based composites. To model the thin-walled nanotubes, special integration procedures for calculating nearly-strongly-singular integrals have been developed. The generic BE substructuring algorithm allows modeling complex CNT-reinforced polymers, containing any number of nanotubes of any shape (straight or curved). The subregion-by-subregion strategy, based on Krylov solvers, makes the independent generation, assembly, and storage of the many parts of the complete BE model possible. Thus, significant memory… More >

Wenzhong Tang¹, Suresh G. Advani^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.1, pp. 31-40, 2007, DOI:10.3970/cmes.2007.022.031

Abstract In this paper, non-equilibrium molecular dynamics (MD) simulations were performed to investigate water flow around a single-walled carbon nanotube. In the simulation, the nanotube was modeled as a rigid cylinder of carbon atoms. Water molecules were described with the extended simple point charge (SPC/E) model. The nanotube-water interactions were calculated with a Lennard-Jones potential between carbon-oxygen pairs. The water-water interactions comprised a Lennard-Jones potential between the oxygen-oxygen pairs and a Coulomb potential between all charge sites on interactive water molecules. It was shown that classical continuum mechanics does not hold when the drag forces on the nanotube are considered. In… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.3, pp. 157-176, 2007, DOI:10.3970/cmes.2007.020.157

Abstract In many engineering applications the Lie group calculation is very important. With this in mind, the subject of this paper is for an in-depth investigation of time-varying linear systems, and its accompanied Lie group calculations. In terms of system matrix A in Eq. (11) and a one-order lower fundamental solution matrix associated with the sub-state matrix function A_s^s, we propose two methods to nilpotentlize the time-varying linear systems. As a consequence, we obtain two different calculations of the general linear group. Then, the nilpotent systems are further transformed to a unique new system Ż(t) = B(t)Z(t), which having a… More >

H. Okada ¹, T. Kamibeppu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.3, pp. 229-238, 2005, DOI:10.3970/cmes.2005.010.229

Abstract In this paper, a three-dimensional VCCM (Virtual Crack Closure-Integral Method) for evaluating the energy release rate and the stress intensity factor is presented. Many engineers and researchers believe that hexahedral finite elements should be used to perform three-dimensional fracture analyses. Previous VCCM formulations assume the use of hexahedral finite elements. In present study, the authors have been developing a VCCM that works with tetrahedral finite elements. In the field of large-scale computation, the use of tetrahedral finite elements has becoming very popular as high performance mesh generation programs became available. Therefore, building a large and complex analysis model with hexahedral… More >

M.Ganesh¹, D. Sheen²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 183-194, 2001, DOI:10.3970/cmes.2001.002.183

Abstract A parallel numerical algorithm is introduced and analyzed for solving inhomogeneous initial-boundary value parabolic problems. The scheme is based on the method recently introduced in Sheen, Sloan, and Thomée (2000) for homogeneous problems. We give a method based on a suitable choice of multiple parameters. Our scheme allows one to compute solutions in a wide range of time. Instead of using a standard time-marching method, which is not easily parallelizable, we take the Laplace transform in time of the parabolic problems. The resulting elliptic problems can be solved in parallel. Solutions are then computed by a discrete inverse Laplace transformation.… More >

R. Mustata¹, S. D. Harris², L. Elliott¹, D. Lesnic¹, D. B. Ingham¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 107-116, 2000, DOI:10.3970/cmes.2000.001.409

Abstract An inverse boundary element method is developed to characterise the components of the hydraulic conductivity tensor K of anisotropic materials. Surface measurements at exposed boundaries serve as additional input to a Genetic Algorithm (GA) using a modified least squares functional that minimises the difference between observed and BEM-predicted boundary pressure and/or hydraulic flux measurements under current hydraulic conductivity tensor component estimates. More >

Qiping Chen¹, Yu Liu¹, Liping Zeng¹, Qiang Xiao¹, Conghui Zhou¹, Sheng Kang¹

FDMP-Fluid Dynamics & Materials Processing, Vol.14, No.4, pp. 281-291, 2018, DOI:10.32604/fdmp.2018.03721

Abstract The brake-by-wire system requirement is promoted owing to the development of green energy vehicle, and the brake pressure control method is needed. A control method for the master cylinder hydraulic pressure based on the sliding mode control approach is proposed to provide the Electro-Hydraulic Brake system (EHB) of electric vehicles with superior system performances. An assessment is carried out about the complex nonlinear characteristics and sensitivity to the external environment of these systems, which include illustrating the working principle of the EHB system, establishing the dynamic models of the key components of the EHB system. The sliding mode control method… More >

Hayder A. Abdulbari ^1,2, Hassan D. Mahammed¹, Z. Hassan, Wafaa K. Mahmood³

FDMP-Fluid Dynamics & Materials Processing, Vol.12, No.2, pp. 86-101, 2016, DOI:10.3970/fdmp.2016.012.086

Abstract An experimental investigation was carried out to study the response of a turbulent pressure drop fluctuations to longitudinal groove riblets, involved two configurations being triangular and spaced triangular grooves with height 600, 800, 1000 μm and peak to peak spacing 1000 μm and 2000 μm respectively. Experiments were therefore performed at free stream velocity up to 0.44 m/sec, which were corresponding to Reynolds number (Re) 53000. The development of the obtained turbulent layer downstream of the grooves was then compared with the results from the corresponding smooth-wall case. To conclude, the effect of the spaced triangular riblets on the turbulent… More >

H. T. Liu¹, Z. D. Han¹, A. M. Rajendran², S. N. Atluri³

CMC-Computers, Materials & Continua, Vol.4, No.1, pp. 43-54, 2006, DOI:10.3970/cmc.2006.004.043

Abstract The Rajendran-Grove (RG) ceramic damage model is a three-dimensional internal variable based constitutive model for ceramic materials, with the considerations of micro-crack extension and void collapse. In the present paper, the RG ceramic model is implemented into the newly developed computational framework based on the Meshless Local Petrov-Galerkin (MLPG) method, for solving high-speed impact and penetration problems. The ability of the RG model to describe the internal damage evolution and the effective material response is investigated. Several numerical examples are presented, including the rod-on-rod impact, plate-on-plate impact, and ballistic penetration. The computational results are compared with available experiments, as well… More >

Qifeng Fan¹, Yaping Zhang², Leiting Dong^1,3, Shu Li¹, Satya N. Atluri⁴

CMC-Computers, Materials & Continua, Vol.47, No.2, pp. 65-88, 2015, DOI:10.3970/cmc.2015.047.065

Abstract A very simple displacement-based hexahedral 32-node element (denoted as DPH32), with over-integration in the thickness direction, is developed in this paper for static and dynamic analyses of laminated composite plates and shells. In contrast to higher-order or layer-wise higher-order plate and shell theories which are widely popularized in the current literature, the proposed method does not develop specific theories of plates and shells with postulated kinematic assumptions, but simply uses the theory of 3-D solid mechanics and the widely-available solid elements. Over-integration is used to evaluate the element stiffness matrices of laminated structures with an arbitrary number of laminae, while… More >