Y. T. Wu¹, J. Z. Cui², Y. F. Nie³, Y. Zhang³

CMC-Computers, Materials & Continua, Vol.42, No.3, pp. 205-226, 2014, DOI:10.3970/cmc.2014.042.205

Abstract Core–shell particle–filled PA6/EPDM–g–MA/HDPE ternary blend has excellent mechanical properties. In this paper, effective elastic properties and tensile yield strength of the ternary blend are predicted by the second–order two– scale method, to investigate the relationship between morphology and mechanical properties. The method and the limit analysis for predicting mechanical properties of random heterogeneous materials are briefly introduced. Realistic morphology of the ternary blend including both core–shell particles and pure particles is simulated, and finite element mesh is generated. The unified strength theory is embedded in the method for the convenience of selecting a suitable yield criterion. The effective elastic moduli… More >

M. Negin¹, S. D. Akbarov^2,3, M. E. Erguven¹

CMC-Computers, Materials & Continua, Vol.42, No.1, pp. 25-62, 2014, DOI:10.3970/cmc.2014.042.025

Abstract Within the framework of the piecewise homogeneous body model the influence of the shear-spring type imperfect contact conditions on the dispersion relation of the generalized Rayleigh waves in the system consisting of the initially stressed covering layer and initially stressed half plane is investigated. The second version of the small initial deformation theory of the three-dimensional linearized theory of elastic waves in initially stressed bodies is applied and the elasticity relations of the materials of the constituents are described by the Murnaghan potential. The magnitude of the imperfectness of the contact conditions is estimated through the shear-spring type parameter. Consequently,… More >

M. Chatiri¹, A. Matzenmiller²

CMC-Computers, Materials & Continua, Vol.35, No.3, pp. 255-283, 2013, DOI:10.3970/cmc.2013.035.255

Abstract This article presents a three dimensional constitutive model for anisotropic damage to describe the elastic-brittle behavior of unidirectional fibrereinforced laminated composites. The primary objective of the article focuses on the three dimensional relationship between damage of the material and the effective elastic properties for the purpose of stress analysis of composite structures, in extension to the two dimensional model in Matzenmiller, Lubliner and Taylor (1995). A homogenized continuum is adopted for the constitutive theory of anisotropic damage and elasticity. Damage initiation criteria are based on Puck failure criterion for first ply failure and progressive micro crack propagation is based on… More >

F. Wang^1,2, D.M. Zhang^1,2,3, H.H. Zhu⁴, H.W. Huang^1,2, J.H. Yin⁵

CMC-Computers, Materials & Continua, Vol.34, No.1, pp. 63-81, 2013, DOI:10.3970/cmc.2013.034.063

Abstract This paper studies the impact of overhead excavation on an existing tunnel through both field monitoring and a full 3D numerical model. It is found that the excavation induced longitudinal heave of the tunnel is uneven with maximum heave occurring below the excavation center. Even at the same cross section, the excavation induced heave is not uniform with the most significant heave occurring at the tunnel crown. The bending moments of the tunnel lining is decreased due to the overhead excavation. The axial forces of the tunnel lining generally decrease except at the tunnel invert. The shear forces of the… More >

C. T. Wu¹, F.-L. Yang², D. L. Young³

CMC-Computers, Materials & Continua, Vol.23, No.2, pp. 135-154, 2011, DOI:10.3970/cmc.2011.023.135

Abstract This paper proposes a method that can calculate the hydrodynamic force of a non-circular object in an inviscid, irrotational, and incompressible flow with the presence of external flow singularities. In order to handle irregular object, the method of fundamental solutions (MFS) is employed to numerically construct the singularity system that describes the body and the flow motion and meets the boundary condition. The obtained singularity system is then integrated into the generalized Lagally theorem to compute the instantaneous hydrodynamic force via algebraic calculations and to describe the unsteady interaction of the object and its ambient flow. The proposed method is… More >

Liviu Marin¹

CMC-Computers, Materials & Continua, Vol.17, No.3, pp. 233-274, 2010, DOI:10.3970/cmc.2010.017.233

Abstract We investigate two algorithms involving the relaxation of either the given boundary temperatures (Dirichlet data) or the prescribed normal heat fluxes (Neumann data) on the over-specified boundary in the case of the iterative algorithm of Kozlov91 applied to Cauchy problems for two-dimensional steady-state anisotropic heat conduction (the Laplace-Beltrami equation). The two mixed, well-posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV)… More >

B. Tomas Johansson¹, Liviu Marin²

CMC-Computers, Materials & Continua, Vol.13, No.2, pp. 153-190, 2009, DOI:10.3970/cmc.2009.013.153

Abstract We propose two algorithms involving the relaxation of either the given Dirichlet data or the prescribed Neumann data on the over-specified boundary, in the case of the alternating iterative algorithm of Kozlov, Maz'ya and Fomin(1991) applied to Cauchy problems for the modified Helmholtz equation. A convergence proof of these relaxation methods is given, along with a stopping criterion. The numerical results obtained using these procedures, in conjunction with the boundary element method (BEM), show the numerical stability, convergence, consistency and computational efficiency of the proposed methods. More >

LiviuMarin ¹

CMC-Computers, Materials & Continua, Vol.12, No.1, pp. 71-100, 2009, DOI:10.3970/cmc.2009.012.071

Abstract In this paper, the alternating iterative algorithm originally proposed by Kozlov, Maz'ya and Fomin (1991) is numerically implemented for the Cauchy problem in anisotropic heat conduction using a meshless method. Every iteration of the numerical procedure consists of two mixed, well-posed and direct problems which are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise… 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 >

P. B. Abhale¹, N. N. Viswanathan¹, N. B. Ballal¹

CMC-Computers, Materials & Continua, Vol.10, No.2, pp. 195-216, 2009, DOI:10.3970/cmc.2009.010.195

Abstract Simulation of gas flow in a multilayered non-isothermal packed bed is useful for blast furnace operators in deciding appropriate charging strategy. While using an anisotropic form of Ergun equation to simulate gas flow through such systems, a new solution methodology for non-isothermal gas with varying density flowing through a layered burden has been proposed. This involves handling non-linearity due to gas density variation with pressure and temperature by solving for the square of pressure instead of pressure directly and handling the non-linearity due to |v| term in the Ergun equation by solving linearized form of Ergun equation and updating |v|… More >