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 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 type partial differential equation 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 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 More >

Hamad F. Al-Harbi¹, Marko Knezevic^1,2, Surya R. Kalidindi^1,3

CMC-Computers, Materials & Continua, Vol.15, No.2, pp. 153-172, 2010, DOI:10.3970/cmc.2010.015.153

Abstract In recent work, we have demonstrated the viability and computational advantages of DFT-based spectral databases for facilitating crystal plasticity solutions in face-centered cubic (fcc) metals subjected to arbitrary deformation paths. In this paper, we extend and validate the application of these novel ideas to body-centered cubic (bcc) metals that exhibit a much larger number of potential slip systems. It was observed that the databases for the bcc metals with a larger number of slip systems were more compact compared to those obtained previously for fcc metals with a smaller number of slip systems. Furthermore, we More >

Dustin D. Gerrard¹, David T. Fullwood¹, Denise M. Halverson², Stephen R. Niezgoda³

CMC-Computers, Materials & Continua, Vol.15, No.2, pp. 129-152, 2010, DOI:10.3970/cmc.2010.015.129

Abstract The effective properties of composite materials are often strongly related to the connectivity of the material components. Many structure metrics, and related homogenization theories, do not effectively account for this connectivity. In this paper, relationships between the topology, represented via homology theory, and the effective elastic response of composite plates is investigated. The study is presented in the context of popular structure metrics such as percolation theory and correlation functions. More >

Lei Bao^1,2, Christophe Schuman¹, Jean-sébastien Lecomte¹, Marie-Jeanne Philippe¹, Xiang Zhao², Liang Zuo², Claude Esling¹

CMC-Computers, Materials & Continua, Vol.15, No.2, pp. 113-128, 2010, DOI:10.3970/cmc.2010.015.113

Abstract The mechanisms of small plastic deformation of titanium (T40) during cold rolling and channel die compression by means of "interrupted in situ" EBSD orientation measurements were studied. These interrupted EBSD orientation measurements allow to determine the rotation flow field which leads to the development of the crystallographic texture during the plastic deformation. Results show that during rolling, tension twins and compression twins occur and various glide systems are activated, the number of grains being larger with twins than with slip traces. In channel die compression, only tension twins are observed in some grains, whereas slip More >

Oliver K. Johnson¹, Calvin J. Gardner¹, David T. Fullwood¹, Brent L.Adams¹, Nathan Hansen², George Hansen²

CMC-Computers, Materials & Continua, Vol.15, No.2, pp. 87-112, 2010, DOI:10.3970/cmc.2010.015.087

Abstract A novel nickel nanostrand-silicone composite material at an optimized 15 vol% filler concentration demonstrates a dramatic piezoresistive effect with a negative gauge factor (ratio of percent change in resistivity to strain). The composite volume resistivity decreases in excess of three orders of magnitude at a 60% strain. The piezoresistivity does decrease slightly as a function of cycles, but not significantly as a function of time. The material's resistivity is also temperature dependent, once again with a negative dependence.
The evidence indicates that nickel strands are physically separated by matrix material even at high volume fractions, and More >

Chia-Ming Fan^1,2, Chein-Shan Liu³, Weichung Yeih¹, Hsin-Fang Chan¹

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 67-86, 2010, DOI:10.3970/cmc.2010.015.067

Abstract In this study, the nonlinear obstacle problems, which are also known as the nonlinear free boundary problems, are analyzed by the scalar homotopy method (SHM) and the finite difference method. The one- and two-dimensional nonlinear obstacle problems, formulated as the nonlinear complementarity problems (NCPs), are discretized by the finite difference method and form a system of nonlinear algebraic equations (NAEs) with the aid of Fischer-Burmeister NCP-function. Additionally, the system of NAEs is solved by the SHM, which is globally convergent and can get rid of calculating the inverse of Jacobian matrix. In SHM, by introducing More >

Chih-Wen Chang¹, Chein-Shan Liu², Jiang-Ren Chang³

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 45-66, 2010, DOI:10.3970/cmc.2010.015.045

Abstract In this article, we propose a semi-analytical method to tackle the two-dimensional backward heat conduction problem (BHCP) by using a quasi-boundary idea. First, the Fourier series expansion technique is employed to calculate the temperature field u(x, y, t) at any time t < T. Second, we consider a direct regularization by adding an extra termau(x, y, 0) to reach a second-kind Fredholm integral equation for u(x, y, 0). The termwise separable property of the kernel function permits us to obtain a closed-form regularized solution. Besides, a strategy to choose the regularization parameter is suggested. When several numerical examples were More >