J. Zidek¹ , J. Kucera¹, J. Jancar¹

CMC-Computers, Materials & Continua, Vol.16, No.1, pp. 51-74, 2010, DOI:10.3970/cmc.2010.016.051

Abstract It has been shown that in particulate filled composites, a cross-property relationship exists between various transport properties (e.g., electrical conductivity, mechanical reinforcement, gas permeation) of a macroscale composite. Thus, knowledge of the effective mechanical properties of a composite immediately places bounds on its electrical conductivity or gas permeation behavior. Using these bounds allows us to predict the phase dispersion state that optimizes one or multiple properties of the composite and, thus, the knowledge of how spatial arrangement of filler particles at their given content affects physical properties of the composite can be valuable. In this paper, a new numerical model… 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 >

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 >

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 a scalar homotopy function and… 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 tested,… More >

Rasa Kazakevičiũtè-Makovska¹

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 27-44, 2010, DOI:10.3970/cmc.2010.015.027

Abstract The phenomenon of stress softening observed in the cyclic inflation of spherical balloons or membranes is quantitatively and qualitatively examined. A new measure of the stress softening extent is proposed which correctly captures the main feature of this phenomenon. This measure of the stress softening is related to the relevant response functions in the constitutive models proposed in the literature to describe this effect. Using these relationships, the predictive capability of the theoretical models is examined. It is shown that only those theoretical models which admit a non-monotone character of the stress softening can properly describe this phenomenon. More >

J. Sladek¹, V. Sladek¹, P. Stanak¹, Ch. Zhang²

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 1-26, 2010, DOI:10.3970/cmc.2010.015.001

Abstract A meshless local Petrov-Galerkin (MLPG) method is applied to solve laminate plate problems described by the Reissner-Mindlin theory. Both stationary and transient dynamic loads are analyzed here. The bending moment and the shear force expressions are obtained by integration through the laminated plate for the considered constitutive equations in each lamina. The Reissner-Mindlin theory reduces the original three-dimensional (3-D) thick plate problem to a two-dimensional (2-D) problem. Nodal points are randomly distributed over the mean surface of the considered plate. Each node is the center of a circle surrounding this node. The weak-form on small subdomains with a Heaviside step… More >

Y.B. Zhang¹, A. Godfrey², D. Juul Jensen¹

CMC-Computers, Materials & Continua, Vol.14, No.3, pp. 197-208, 2009, DOI:10.3970/cmc.2009.014.197

Abstract Two methods to quantify protrusions/retrusions and to estimate local boundary curvature from sample plane sections are proposed. The methods are used to evaluate the driving force due to curvature of the protrusions/retrusions for partially recrystallized pure nickel cold rolled to 96% reduction in thickness. The results reveal that the values calculated by both these methods are reasonable when compared with the stored energy measured by differential scanning calorimetry. The relationship between protrusions and the average stored energy density in the deformed matrix is also investigated for partially recrystallized pure aluminum cold rolled to 50%. The results show that the local… More >