Giacomo Landi¹, Surya R. Kalidindi²

CMC-Computers, Materials & Continua, Vol.16, No.3, pp. 273-294, 2010, DOI:10.3970/cmc.2010.016.273

Abstract In this paper, we present a computationally efficient multi-scale framework for predicting the local fields in the representative volume element of a multiphase material system subjected to thermo-mechanical loading conditions. This framework for localization relationships is a natural extension of our recent work on two-phase composites subjected to purely mechanical loading. In this novel approach, the localization relationships take on a simple structure expressed as a series sum, where each term in the series is a convolution product of local structure and the governing physics expressed in the form of influence coefficients. Another salient feature of this approach is its… More >

Xue-Qian Fang^1,2, Xiao-Hua Wang¹, Le-Le Zhang³

CMC-Computers, Materials & Continua, Vol.16, No.3, pp. 229-246, 2010, DOI:10.3970/cmc.2010.016.229

Abstract In the design of advanced micro- and nanosized materials and devices containing inclusions, the effects of surfaces/interfaces on the stress concentration become prominent. In this paper, based on the surface/interface elasticity theory, a two-dimensional problem of an elliptical nano-inhomogeneity under anti-plane shear waves is considered. The conformal mapping method is then applied to solve the formulated boundary value problem. The analytical solutions of displacement fields are expressed by employing wave function expansion method, the expanded mode coefficients are determined by satisfying the boundary conditions at the interfaces of the nano-inhomogeneity. Analyses show that the effect of the interfacial properties on… More >

W. F. Yuan^1,2, H. B. Zhang¹

CMC-Computers, Materials & Continua, Vol.16, No.2, pp. 117-126, 2010, DOI:10.3970/cmc.2010.016.117

Abstract Corrosion pits on a steel plate can reduce the strength of the plate. However, it is difficult to calculate the corrosion effect analytically since the pits are normally distributed on the plate's surface randomly. In this manuscript, a simple approach is proposed to convert the corroded plate into a perfect one. By this method, the corrosion pits are treated as inclusions embedded in the plate. Then the analytical mechanics model used for composite material can be adopted in the calculation of the steel plate's effective material properties. To verify the proposed approach, numerical simulation is conducted using finite element method. More >

Veturia Chiroiu¹, Ligia Munteanu¹, Pier Paolo Delsanto²

CMC-Computers, Materials & Continua, Vol.16, No.1, pp. 75-100, 2010, DOI:10.3970/cmc.2010.016.075

Abstract Conventional continuum theories are unable to capture the observed indentation size effects, due to the lack of intrinsic length scales that represent the measures of nanostructure in the constitutive relations. In order to overcome this deficiency, the Toupin-Mindlin strain gradient theory of nanoindentation is formulated in this paper and the size dependence of the hardness with respect to the depth and the radius of the indenter for multiple walled carbon nanotubes is investigated. Results show a peculiar size influence on the hardness, which is explained via the shear resistance between the neighboring walls during the buckling of the multiwalled nanotubes. More >

M. Vogler¹, G. Ernst¹, R. Rolfes¹

CMC-Computers, Materials & Continua, Vol.16, No.1, pp. 25-50, 2010, DOI:10.3970/cmc.2010.016.025

Abstract In this article, a constitutive formulation of a transversely-isotropic material and failure model for fiber-reinforced polymers is presented comprising pre-failure material nonlinearities, a novel invariant based quadratic failure criterion (IQC) as well as post failure material softening. The failure surface of the IQ criterion is assumed to take the influence of triaxiality on fracture into account. Further, a distinction between fiber failure and inter-fiber failure is conducted. Material softening is governed by a fracture energy formulation and the introduction of an internal length. The constitutive model is implemented into a programming user interface of the commercial finite element program Abaqus.… More >

Chen Y J^1,2, Huber N^2,3

CMC-Computers, Materials & Continua, Vol.16, No.1, pp. 1-24, 2010, DOI:10.3970/cmc.2010.016.001

Abstract An efficient wear integration algorithm is crucial for the simulation of wear in complex transient contact situations. By rewriting Archard's wear law for two dimensional problems, the wear integration can be replaced by the total contact force. This avoids highly resolved simulations in time and space, so that the proposed method allows a significant acceleration of wear simulations. All quantities, including the average contact velocity, slip rate and total contact force, which are required for the pressure-force transformation, can be determined from geometric and motion analysis, or alternatively, from Finite Element simulations. The proposed CForce method has been implemented into… More >

Luming Shen¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 241-260, 2010, DOI:10.3970/cmc.2010.015.241

Abstract In this study, molecular dynamics (MD) simulations are performed to form glassy silica from meltedb-cristobalite using cooling rates of 2, 20 and 200 K/ps. The resulting glassy silica samples are then shocked at particle velocities ranging from 0.3 to 11 km/s in the MD simulations. The effect of the cooling rate on the shock wave velocity is observed for particle velocities below 2 km/s. Moreover, the simulated pressure and density of the shocked glassy silica increase as the cooling rate increases. As compared with the experimental data, the MD simulation can approximately identify the initiation of densification and predict the… More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 221-240, 2010, DOI:10.3970/cmc.2010.015.221

Abstract We propose a simple numerical scheme for solving the space- and time-fractional derivative Burgers equations: D_t^αu + εuu_x = vu_xx + ηD_x^βu, 0 < α, β ≤ 1, and u_t + D_*^β(D_*^1-βu)²/2 = vu_xx, 0 < β ≤ 1. The time-fractional derivative D_t^αu and space-fractional derivative D_x^βu are defined in the Caputo sense, while D_*^βu is the Riemann-Liouville space-fractional derivative. A fictitious time τ is used to transform the dependent variable u(x,t) into a new one by (1+τ)^γu(x,t) =: v(x,t,τ), where 0 < γ ≤ 1 is a parameter, such that the original equation is written as a new functional-differential… More >

F.Y.Li^1,2, G.Y.Li¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 199-220, 2010, DOI:10.3970/cmc.2010.015.199

Abstract In this paper, an uncertain multi-objective optimization method is suggested to deal with crashworthiness design problem of vehicle, in which the uncertainties of the parameters are described by intervals. Considering both lightweight and safety performance, structural weight and peak acceleration are selected as objectives. The occupant distance is treated as constraint. Based on interval number programming method, the uncertain optimization problem is transformed into a deterministic optimization problem. The approximation models are constructed for objective functions and constraint based on Latin Hypercube Design (LHD). Thus, the interval number programming method is combined with the approximation model to solve the uncertain… More >

Hrvoje Gotovac¹, Vedrana Kozulić², Blaž Gotovac¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 173-198, 2010, DOI:10.3970/cmc.2010.015.173

Abstract The space-time Adaptive Fup Collocation Method (AFCM) for solving boundary-initial value problems is presented. To solve the one-dimensional initial boundary value problem, we convert the problem into a two-dimensional boundary value problem. This quasi-boundary value problem is then solved simultaneously in the space-time domain with a collocation technique and by using atomic Fup basis functions. The proposed method is a generally meshless methodology because it requires only the addition of collocation points and basis functions over the domain, instead of the classical domain discretization and numerical integration. The grid is adapted progressively by setting the threshold as a direct measure… More >