Ankush Aggarwal¹, Jiun-Shyan Chen², William S. Klug³

CMES-Computer Modeling in Engineering & Sciences, Vol.98, No.1, pp. 69-99, 2014, DOI:10.3970/cmes.2014.098.069

Abstract Mechanical properties of proteins play an important role in their biological function. For example, microtubules carry large loads to transport organelles inside the cell, and virus shells undergo changes in shape and mechanical properties during maturation which affect their infectivity. Various theoretical models including continuum elasticity have been applied to study these structural properties, and a significant success has been achieved. But, the previous frameworks lack a connection between the atomic and continuum descriptions. Here this is accomplished through the development of a meshfree framework based on reproducing kernel shape functions for the large deformation mechanics of protein structures. The… More >

Chu Wang, Lucy T. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.98, No.1, pp. 41-67, 2014, DOI:10.3970/cmes.2014.098.041

Abstract Bloodstain pattern analysis (BPA) in forensic science is an important tool to solve crime scenes. The complex dynamic behavior of blood drops poses great challenges for accurate fluid dynamic simulations. In this paper, we specifically focus on simulations of blood drop spreading and impact, which may involve contact line hysteresis and spattering of drops as they interact with solid surfaces. Here, we set up a numerical framework that combines (1) the connectivity-free front tracking (CFFT) method for modeling multiphase (air and liquid) flows and (2) a dynamic contact line model for modeling fluid-solid contact line. Both components are necessary in… More >

R.A. Regueiro^1,2, B. Zhang², S.L. Wozniak³

CMES-Computer Modeling in Engineering & Sciences, Vol.98, No.1, pp. 1-39, 2014, DOI:10.3970/cmes.2014.098.001

Abstract The paper presents three-dimensional, large deformation, coupled finite element analysis (FEA) of dynamic loading on soft biological tissues treated as biphasic (solid-fluid) porous media. An overview is presented of the biphasic solidfluid mixture theory at finite strain, including inertia terms. The solid skeleton is modeled as an isotropic, compressible, hyperelastic material. FEA simulations include: (1) compressive uniaxial strain loading on a column of lung parenchyma with either pore air or water fluid, (2) out-of-plane pressure loading on a thin slab of lung parenchyma with either pore air or water fluid, and (3) pressure loading on a 1/8th symmetry vertebral disc… More >

K.P. Baxevanakis¹, A.E. Giannakopoulos²

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.6, pp. 535-544, 2014, DOI:10.3970/cmes.2014.097.535

Abstract This article has no abstract. More >

Tao Zhang^1,2, Yiqian He³, Leiting Dong⁴, Shu Li¹, Abdullah Alotaibi⁵, Satya N. Atluri^2,5

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.6, pp. 509-533, 2014, DOI:10.3970/cmes.2014.097.509

Abstract In this article, the Meshless Local Petrov-Galerkin (MLPG) Mixed Collocation Method is developed to solve the Cauchy inverse problems of Steady- State Heat Transfer In the MLPG mixed collocation method, the mixed scheme is applied to independently interpolate temperature as well as heat flux using the same meshless basis functions The balance and compatibility equations are satisfied at each node in a strong sense using the collocation method. The boundary conditions are also enforced using the collocation method, allowing temperature and heat flux to be over-specified at the same portion of the boundary. For the inverse problems where noise is… More >

Xiaolin Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.6, pp. 483-507, 2014, DOI:10.3970/cmes.2014.097.483

Abstract Combining moving least square (MLS) approximations and boundary integral equations, a symmetric and boundary-only meshless method, the Galerkin boundary node method (GBNM), is developed in this paper for two- and threedimensional elasticity problems with mixed boundary conditions. Unlike other MLS-based meshless methods, boundary conditions in this meshless method can be applied directly and easily. In the GBNM, the stiffness matrices so obtained are symmetric. The property of symmetry is an added advantage in coupling the GBNM with the finite element method (FEM). Thus, a symmetric coupling of the GBNM and the FEM is also discussed for elasticity problems. Error analysis… More >

Jin Li^1,2, Hongxing Ru¹, Dehao Yu^3,4

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.6, pp. 463-482, 2014, DOI:10.3970/cmes.2014.097.463

Abstract The computation with Simpson’s rule for the supersingular integrals on circle is discussed. When the singular point coincides with some priori known point, the convergence rate of the Simpson rule is higher than the globally one which is considered as the superconvergence phenomenon. Then the error functional of density function is derived and the superconvergence phenomenon of composite Simpson rule occurs at certain local coordinate of each subinterval. Based on the error functional, a modify quadrature is presented. At last, numerical examples are provided to validate the theoretical analysis and show the efficiency of the algorithms. More >

Ji-Chuan Liu¹, Jun-Gang Wang²

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.5, pp. 437-462, 2014, DOI:10.3970/cmes.2014.097.437

Abstract We consider the determination of heat flux within a body from the Cauchy data. The aim of this paper is to seek an approach to solve the onedimensional heat equation in a bounded domain without initial value. This problem is severely ill-posed and there are few theoretic results. A quasi-reversibility regularization method is used to obtain a regularized solution and convergence estimates are given. For numerical implementation, we apply a method of lines to solve the regularized problem. From numerical results, we can see that the proposed method is reasonable and feasible. More >

Mao Liu^1,2, Jiawei Xiang¹, Haifeng Gao¹, Yongying Jiang¹, Yuqing Zhou¹, Fengping Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.5, pp. 425-436, 2014, DOI:10.3970/cmes.2014.097.425

Abstract A wavelet finite element method (WFEM) is developed to analyze the dispersion relation for one-dimensional phononic crystals (1DPCs). In order to calculate the band gaps (BGs) of 1DPCs, the wavelet finite element model is constructed using a slender beam element based on B-spline wavelet on the interval (BSWI). Combining with the Bloch-Floquet theorem and ω(k) technique, the model will be simplified as a simple eigenproblem. The performance of the proposed method has been numerically verified by one numerical example. More >

S.D. Akbarov ^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.97, No.5, pp. 407-424, 2014, DOI:10.3970/cmes.2014.097.407

Abstract Flexural wave dispersion in a bi-material solid and hollow circular cylinders is investigated with the use of the three-dimensional linear theory of elastodynamics. It is assumed that on the interface surface of the cylinders the complete contact conditions satisfy. The analytical solution of the corresponding field equations is presented and, using these solutions, the dispersion equations for the cases under consideration are obtained. The dispersion equations are solved numerically and based on these solutions, dispersion curves are constructed for the concrete selected pairs of materials such as Tungsten (inner cylinder material) + Aluminum (outer cylinder material) and Steel (inner cylinder… More >