Leiting Dong¹, Satya N. Atluri1¹

CMC-Computers, Materials & Continua, Vol.30, No.1, pp. 39-82, 2012, DOI:10.3970/cmc.2012.030.039

Abstract In this paper, as an extension to the authors's work in [Dong and Atluri (2011a,b, 2012a,b,c)], three-dimensional Trefftz Voronoi Cells (TVCs) with ellipsoidal voids/inclusions are developed for micromechanical modeling of heterogeneous materials. Several types of TVCs are developed, depending on the types of heterogeneity in each Voronoi Cell(VC). Each TVC can include alternatively an ellipsoidal void, an ellipsoidal elastic inclusion, or an ellipsoidal rigid inclusion. In all of these cases, an inter-VC compatible displacement field is assumed at each surface of the polyhedral VC, with Barycentric coordinates as nodal shape functions. The Trefftz trial displacement… More >

L. Dong¹, S. N. Atluri¹

CMC-Computers, Materials & Continua, Vol.29, No.2, pp. 169-212, 2012, DOI:10.3970/cmc.2012.029.169

Abstract In this paper, three-dimensionalT-Trefftz Voronoi Cell Finite Elements (VCFEM-TTs) are developed for micromechanical modeling of heterogeneous materials. Several types of VCFEMs are developed, depending on the types of heterogeneity in each element. Each VCFEM can include alternatively a spherical void, a spherical elastic inclusion, a spherical rigid inclusion, or no voids/inclusions at all.In all of these cases, an inter-element compatible displacement field is assumed at each surface of the polyhedral element, with Barycentric coordinates as nodal shape functions.The T-Trefftz trial displacement fields in each element are expressed in terms of the Papkovich-Neuber solution. Spherical harmonics… More >

K. Choukairy¹, R. Bennacer²

FDMP-Fluid Dynamics & Materials Processing, Vol.8, No.2, pp. 155-172, 2012, DOI:10.3970/fdmp.2012.008.155

Abstract This study carries the natural thermosolutal convection induced in heterogeneous porous media. The configuration considered is cartesian. The horizontal and vertical walls are submitted to different mass and heat transfer. The equations which govern this type of flow are solved numerically by using the finite volume method. The flow is considered two-dimensional and laminar. The model of Darcy and the approximation of the Boussinesq are taken into account. The parameters which control the problem are the thermal Darcy-Rayleigh number, Rt, the buoyancy ratio, N, the Lewis number, Le, the aspect ratio of the enclosure, A… More >

De-Shin Liu, Zhen-Wei Zhuang, Cho-Liang Chung

The International Conference on Computational & Experimental Engineering and Sciences, Vol.19, No.2, pp. 39-40, 2011, DOI:10.3970/icces.2011.019.039

Abstract A two-dimensional heterogeneous hybrid moisture element method (HHMEM) for modeling transient moisture diffusion in permeable fiber-reinforced polymer composites is proposed in this paper.
The HHMEM scheme is based on a heterogeneous hybrid moisture element (HHME), whose properties are determined by equivalent moisture capacitance and conductance matrixes calculated using the conventional finite element formulation with the similarity mass/stiffness property and matrix condensing operations. A coupled HHME-FE scheme is developed and implemented in computer codes MATLAB to analyze the transient moisture diffusion characteristics of polymeric composite materials containing multiple permeable fibers. The analysis commences by comparing the… More >

A.P. Selvadurai

The International Conference on Computational & Experimental Engineering and Sciences, Vol.16, No.3, pp. 67-68, 2011, DOI:10.3970/icces.2011.016.067

Abstract Theme Lecture More >

De-Shin Liu¹, Zhen-Wei Zhuang¹, Shaw-Ruey Lyu^2,3, Cho-Liang Chung⁴, Pai-Chen Lin¹

CMC-Computers, Materials & Continua, Vol.26, No.2, pp. 111-136, 2011, DOI:10.3970/cmc.2011.026.111

Abstract This study proposes a two-dimensional heterogeneous hybrid moisture element method (HHMEM) for modeling transient moisture diffusion in permeable fiber-reinforced polymer composites.
The HHMEM scheme is based on a heterogeneous hybrid moisture element(HHME), with properties determined through an equivalent hybrid moisture capacitance/conductance matrix. This matrix was calculated using the conventional finite element formulation in space discretization as well as the θ-method in time discretization with similar mass/stiffness properties and matrix condensing operations. A coupled HHME-FE scheme was developed and implemented in computer code MATLAB in order to analyze the transient moisture diffusion characteristics of composite materials containing multiple… More >

Marco Dondero¹, Adrián P. Cisilino^1,2, J. Pablo Tomba¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 225-246, 2011, DOI:10.3970/cmes.2011.078.225

Abstract This paper introduces a numerical methodology for the design of random micro-heterogeneous materials with functionally graded effective thermal conductivities (ETC). The optimization is carried out using representative volume elements (RVEs), a parallel Genetic Algorithm (GA) as optimization method, and a Fast Multipole Boundary Element Method (FMBEM) for the evaluation of the cost function. The methodology is applied for the design of foam-like microstructures consisting of random distributions of circular insulated holes. The temperature field along a material sample is used as objective function, while the spatial distribution of the holes is the design variable. There More >

A. Bakkali¹, L. Azrar^1,2, N. Fakri¹

CMC-Computers, Materials & Continua, Vol.23, No.3, pp. 201-232, 2011, DOI:10.3970/cmc.2011.023.201

Abstract In this paper an N-phase Incremental Self Consistent model is developed for magnetoelectroelastic composites as well as the N-phase Mori-Tanaka and classical Self Consistent. Our aim here is to circumvent the limitation of the Self Consistent predictions for some coupling effective properties at certain inclusion volume fractions. The anomalies of the SC estimates are more drastic when the void inclusions are considered. The mathematical modeling is based on the heterogeneous inclusion problem of Eshelby which leads to an expression for the strain-electric-magnetic field related by integral equations. The effective N-phase magnetoelectroelastic moduli are expressed as More >

Zhenhan Yao¹, Lingfei Gao¹

CMC-Computers, Materials & Continua, Vol.9, No.1, pp. 1-24, 2009, DOI:10.3970/cmc.2009.009.001

Abstract To improve the effectiveness of the lattice model for simulating the failure processes of heterogeneous brittle materials, each lattice element is refined as a subdomain with homogenous material, and is modeled by the boundary element method in this paper. For simplicity, each subdomain is modeled with constant boundary elements. To enhance the efficiency, a row of sub-domains is formed, and then a chain structure of such row-subdomain is constructed. The row-equation systems are solved one by one, and then back substituted, to obtain the final solution. Such a chain subdomain approach of the boundary element… More >

De-Shin Liu¹, Zhen-Wei Zhuang^1,2, Cho-LiangChung³, Ching-Yang Chen⁴

CMC-Computers, Materials & Continua, Vol.13, No.2, pp. 89-114, 2009, DOI:10.3970/cmc.2009.013.089

Abstract This paper employs a novel numerical technique, designated as the hybrid moisture element method (HMEM), to model and analyze moisture diffusion in a heterogeneous epoxy resin containing multiple randomly distributed particles. The HMEM scheme is based on a hybrid moisture element (HME), whose properties are determined by equivalent moisture capacitance and conductance matrixes calculated using the conventional finite element formulation. A coupled HME-FE scheme is developed and implemented using the commercial FEM software ABAQUS. The HME-FE scheme is then employed to analyze the moisture diffusion characteristics of a heterogeneous epoxy resin layer containing particle inclusions. The… More >