A. Tadeu¹, P. Stanak², J. Sladek², V. Sladek², J. Prata¹, N. Simões¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.6, pp. 489-516, 2013, DOI:10.3970/cmes.2013.093.489

Abstract This paper presents a technique that couples the boundary element method (BEM) with the meshless local Petrov-Galerkin (MLPG) method, formulated in the frequency domain. It is then used to study the transient heat diffusion through a two-dimensional unbounded medium containing confined subdomains where the material properties vary from point to point. To exploit the advantages of each method, the BEM is used for the homogeneous unbounded domain and the MLPG method is used for the non-homogeneous confined subdomains. The nodal points placed at the interface between the confined subdomains and the unbounded homogenous medium are used to couple the BEM… More >

W. Zhao^1,2, J.K. Liu³, X.Y. Li², Q.W. Yang⁴, Y.Y. Chen⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.6, pp. 469-488, 2013, DOI:10.3970/cmes.2013.093.469

Abstract In order to obtain reliable structural design, it is of extreme importance to evaluate the failure probability, safety levels of structure (reliability analysis) and the effect of a change in a variable parameter on structural safety (sensitivity analysis) when uncertainties are considered. With a computationally cheaper approximation of the limit state function, various response surface methods (RSMs) have emerged as a convenient tool to solve this especially for complex problems. However, the traditional RSMs may produce large errors in some conditions especially for those highly non-linear limit state functions. Instead of the traditional least squares approximation, in the present paper,… More >

D.S. Liu^1,2, Z.H. Fong¹, I.H. Lin¹, Z.W. Zhuang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.6, pp. 441-468, 2013, DOI:10.3970/cmes.2013.093.441

Abstract In this study, we proposed a novel numerical technique to simulate the transient moisture diffusion process and to apply it to heterogeneous composite resins. The method is based on a heterogeneous hybrid moisture element (HHME), with properties determined through an equivalent hybrid moisture capacitance/ conductance matrix that was calculated using the conventional finite element formulation in space discretization and the q-method in time discretization, with similar mass/stiffness properties and matrix condensing operations. A coupled HHME with finite element scheme was developed and implemented in the computer code by using the commercial software MATLAB to analyze the transient moisture diffusion process… More >

Marcin Kamiński^1,2, Bernd Lauke²

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.6, pp. 411-440, 2013, DOI:10.3970/cmes.2013.093.411

Abstract The main aim in this paper is a computational study devoted to the sensitivity gradients and probabilistic moments of the effective elastic parameters for the rubber-filled polymers. The methodology is based on least squares recovery of the polynomial functions relating the effective tensor components and the given input design/random parameters. All numerical experiments are provided with respect to Young’s moduli of the elastomer constituents. Computational analysis is possible thanks to the application of the Response Function Method, which is enriched in our approach with the weighting procedures implemented according to the Dirac-type distributions. The homogenized elasticity tensor components are derived… More >

Dan Xie¹, Min Xu²

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.5, pp. 377-409, 2013, DOI:10.3970/cmes.2013.093.377

Abstract We apply a simple proper orthogonal decomposition (POD) method to compute the nonlinear oscillations of a degenerate two-dimensional fluttering plate undergoing supersonic flow. First, the von Karman’s large deflection theory and quasi-steady aerodynamic theory are employed in constructing the governing equations of the simply supported plate. Then, the governing equations are solved by both the Galerkin method and the POD method. The Galerkin method is accurate but sometimes computationally expensive, since the number of degrees of freedom (dofs) required is relatively large provided that nonlinearity is strong. The POD method can be used to capture the complex dynamics of a… More >

K. Kakuda¹, T. Nagashima¹, Y. Hayashi¹, S. Obara¹, J. Toyotani¹, S. Miura², N. Katsurada³, S. Higuchi³, S. Matsuda³

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.5, pp. 363-376, 2013, DOI:10.3970/cmes.2013.093.363

Abstract The application of a GPU-based particle method to three-dimensional incompressible viscous fluid flow problems is presented. The particle approach is based on the MPS (Moving Particle Semi-implicit) scheme using logarithmic weighting function to stabilize the spurious oscillatory solutions for solving the Poisson equation with respect to the pressure fields by using GPU-based SCG (Scaled Conjugate Gradient) method. Numerical results demonstrate the workability and the validity of the present approach through the dam-breaking flow problem and flow behavior in a liquid ring pump with rotating impeller blades. More >

G. P. Nikishkov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.5, pp. 335-361, 2013, DOI:10.3970/cmes.2013.093.335

Abstract Accuracy of the quarter-point and transition elements is investigated on one- and two-dimensional problems with inverse square-root singularity. It is demonstrated that most coefficients of the stiffness matrix of the quarter-point element are unbounded. However, numerical integration produces finite values of these coefficients. Influence of several parameters on the error in determining the stress intensity factor is studied. Solution accuracy can be improved using special distribution of element sizes and increasing the element integration order in the radial direction. More >

Tai C. Lee¹, Huan J. Keh²

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.5, pp. 317-333, 2013, DOI:10.3970/cmes.2013.093.317

Abstract A theoretical study of the thermocapillary migration of a fluid sphere located at an arbitrary position inside a spherical cavity is presented in the quasisteady limit of small Reynolds and Marangoni numbers. The applied temperature gradient is perpendicular to the line through the drop and cavity centers. The general solutions to the energy and momentum equations governing the system are constructed from the superposition of their fundamental solutions in the spherical coordinates originating from the two centers, and the boundary conditions are satisfied by a multipole collocation method. Results for the thermocapillary migration velocity of the drop are obtained for… More >

Bing Li^1,2, Zhengjia He¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.4, pp. 293-316, 2013, DOI:10.3970/cmes.2013.093.293

Abstract The vibration-based crack identification problem insists of finding a measured vibration parameter from a complete crack-detection-database constructed by numerical simulation. It is one of the classical optimization problems. Many intelligence methods, such as neural network (NN), genetic algorithm (GA), determinant transformation (DT), and frequency contour (FC) etc., have been extensively employed as optimization tools to achieve this task. The aim of this paper is to propose a benchmark problem to compare these extensive-used optimization methods in terms of crack identification precision and computational time. The merit and demerits for each method are discussed. The results suggest that FC is a… More >

Xiaokui Yue¹, Yong Xu², Jianping Yuan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.4, pp. 277-291, 2013, DOI:10.3970/cmes.2013.093.277

Abstract In this paper, a stochastic averaging method is derived for a class of non-linear stochastic systems under parametrical Poisson white noise excitation, which may be used to model the beam-beam interaction models in particle accelerators. The averaged Generalized Fokker-Planck equation is derived and the approximate stationary solution of the averaged Generalized Fokker-Planck equation is solved by using perturbation method. The present method applied in this paper can reduce the dimensions of stochastic ODE from 2n to n, which simplify the complex stochastic ODE, and then the analytical stationary solutions can be obtained. An example is employed to demonstrate the procedure… More >