Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.21, No.2, pp. 87-106, 2011, DOI:10.3970/cmc.2011.021.087

Abstract We address a new numerical approach to deal with these multi-dimensional backward wave problems (BWPs) in this study. A fictitious time τ is utilized to transform the dependent variable u(x, y, z, t) into a new one by (1+τ)u(x, y, z, t)=: v(x, y, z, t, τ), such that the original wave equation is written as a new hyperbolic type partial differential equation in the space of (x, y, z, t, τ). Besides, a fictitious viscous damping coefficient can be employed to strengthen the stability of numerical integration of the discretized equations by using a group preserving scheme. Several numerical… More >

Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.19, No.3, pp. 285-314, 2010, DOI:10.3970/cmc.2010.019.285

Abstract In this article, we propose a new numerical approach for solving these multi-dimensional nonlinear and nonhomogeneous backward heat conduction problems (BHCPs). A fictitious time t is employed to transform the dependent variable u(x, y, z, t) into a new one by (1+t)u(x, y, z, t)=: v(x, y, z, t, t), such that the original nonlinear and nonhomogeneous heat conduction equation is written as a new parabolic type partial differential equation in the space of (x, y, z, t, t). In addition, a fictitious viscous damping coefficient can be used to strengthen the stability of numerical integration of the discretized equations… More >

A. Rama Ch,ra Murthy^1,2, G.S. Palani¹, Nagesh R. Iyer^1,3, M Srinivasa Pavan¹, Smitha Gopinath¹

CMC-Computers, Materials & Continua, Vol.19, No.2, pp. 135-154, 2010, DOI:10.3970/cmc.2010.019.135

Abstract This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SIF accounting for tension softening effect has been obtained as the difference of SIF obtained using linear elastic fracture mechanics (LEFM) principles and SIF due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SIF due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In… More >

Chih-Wen Chang¹, Chein-Shan Liu²

CMC-Computers, Materials & Continua, Vol.17, No.1, pp. 19-40, 2010, DOI:10.3970/cmc.2010.017.019

Abstract In this study, we employ a semi-analytical approach to solve a two-dimensional advection-dispersion equation (ADE) for identifying the contamination problems. First, the Fourier series expansion technique is used to calculate the concentration field C(x, y, t) at any time t < T. Then, we ponder a direct regularization by adding an extra termaC(x, y, 0) on the final time data C(x, y, T), to reach a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us obtaining a closed-form solution of the Fourier coefficients. A strategy to choose the regularization parameter is offered. The solver utilized in… 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 >

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 >

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 >

Amiya K. Samanta¹, Somnath Ghosh²

CMC-Computers, Materials & Continua, Vol.12, No.1, pp. 1-38, 2009, DOI:10.3970/cmc.2009.012.001

Abstract This paper presents a formulation of hypo-elasticity based RC beam model with bond-slip. Details of the constitutive model and analysis method used are provided. A procedure has been described to carry out three-dimensional analysis considering both geometrical as well as material nonlinearity for a simply supported RC beam employing finite element technique, which uses 8-noded isoparametric hexahedral element HCiS18. Enhanced assumed strain (EAS) formulation has been utilized to predict load-deformation and internal stresses both in the elastic as well as nonlinear regime. It models the composite behaviour of concrete and reinforcements in rigid /perfect bond situation and their mutual interaction… More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.11, No.1, pp. 15-32, 2009, DOI:10.3970/cmc.2009.011.015

Abstract Motivated by the evolutionary and dissipative properties of parabolic type partial differential equation (PDE), Liu (2008a) has proposed a natural and mathematically equivalent approach by transforming the quasilinear elliptic PDE into a parabolic one. However, the above paper only considered a rectangular domain in the plane, and did not treat the difficulty arisen from the quasilinear PDE defined in an arbitrary plane domain. In this paper we propose a new technique of internal and boundary residuals in a fictitious rectangular domain, which are driving forces for the ordinary differential equations based on the Fictitious Time Integration Method (FTIM). Several numerical… More >