Huimin Song¹, Dewey H. Hodges¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 239-278, 2012, DOI:10.3970/cmes.2012.085.239

Abstract The present paper presents a rigorous approach that can accurately and efficiently capture the linear, static and free-vibration behaviors of a beam-like structure by the rigorous combination of a one-dimensional beam model with a three-dimensional continuum model. This study focuses on coupling these disparate finite element types, putting them both into a single finite element model while making use of the asymptotically exact information available as part of the beam model, which itself is obtained by asymptotic dimensional reduction. The coupling is undertaken by use of appropriate transformation matrices at the interface together with stress and displacement recovery relations that… More >

Pengfei Yang¹, Yan Liu¹, Xiong Zhang^1,2, Xu Zhou³, Yuli Zhao³

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 207-238, 2012, DOI:10.3970/cmes.2012.085.207

Abstract The material point method is extended to the simulations of fragmentation driven by detonation. A crack modeling scheme based on contact algorithm with material failure process is developed to study the dynamic crack propagation in plastic media. When considering microscopic damage of material, the plastic behavior is described by Gurson model with randomly-distributed initial void of material points. Gurson model can degenerate to J2 plastic theory while the microscopic void is ignored, in which situation the Weibull random failure scheme will be used. Meanwhile, a background-grid-based searching method is proposed to capture the statistical feature of the fragmentation. The scaling… More >

H.W. Tang¹, Y.T. Yang¹, C.K. Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 193-206, 2012, DOI:10.3970/cmes.2012.085.193

Abstract In this article, a new method combining the Mathematical Programming and the Method of Weighted Residual called MP-MWR is presented. Under the validation of maximum principle, and up on the collocation method, the differential equation can be transferred into a bilateral inequality problem. Applying the genetic algorithms helps to find optimal solutions of upper and lower bounds which satisfy the inequalities. Here, the method is verified by analyzing the deflection of superelliptical clamped plate problem. By using this method, the good approximate solution and its error bounds can be obtained effectively and accurately. More >

Souli, M.¹, Gabrys, J.²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 177-192, 2012, DOI:10.3970/cmes.2012.085.177

Abstract Bird impacts on aircraft are very common and cause significant safety threats to commercial and military aircraft. According to FAA ( Federal American Aviation) regulations, aircraft should be able to land safely following specified types of bird impact on components such as radomes, windshields, engines leading edge structures and other exposed components. Thus exposed components are required to be certified for bird impact. In order to evaluate whether the aircraft is compliant to FAA requirements, several experimental tests and numerical simulations of bird impact on components need to be preformed. This paper presents an experimental and numerical investigation of bird… More >

V.G.Yakhno¹, D. Ozdek²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 157-176, 2012, DOI:10.3970/cmes.2012.085.157

Abstract The present paper describes computation of the time-dependent Green's function for the equations of longitudinal vibration of a multi-step rod with a piecewise constant varying cross-section. This computation is based on generalization of the Fourier series expansion method. The time-dependent Green's function is derived in the form of the Fourier series. The basic functions of this series are eigenfunctions of an ordinary differential equation with boundary and matching conditions. Constructing the eigenvalues and eigenfunctions of this differential equation and then derivation of the Fourier coefficients of the Green's function are main steps of the method. Computational experiments confirm the robustness… More >

Zhibo Yang¹, Xuefeng Chen², Bing Li¹, Zhengjia He¹, Huihui Miao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 129-156, 2012, DOI:10.3970/cmes.2012.085.129

Abstract The implementation of the B-spline Wavelet on the Interval (BSWI) for curved shell elements with rectangular planform is presented in this paper. By aid of the general shell theory, cylinder shells, doubly-curved shallow shells and hyperbolic paraboloidal shells BSWI elements are formulated. Instead of traditional polynomial interpolation, scaling functions at certain scale have been adopted to form the shape functions and construct wavelet-based elements. Because of the good character of BSWI scaling functions, the BSWI curved shell elements combine the accuracy of wavelet-based elements approximation and the character of B-spline functions for structural analysis. Different from the flat shell elements,… More >

Roman Chapko¹, B. Tomas Johansson²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 105-128, 2012, DOI:10.3970/cmes.2012.085.105

Abstract We consider a Cauchy problem for the Laplace equation in a 3-dimen -sional semi-infinite domain that contains a bounded inclusion. The canonical situation is the upper half-space in I\tmspace -.1667em R3 containing a bounded smooth domain. The function value of the solution is specified throughout the plane bounding the upper half-space, and the normal derivative is given only on a finite portion of this plane. The aim is to reconstruct the solution on the surface of the bounded inclusion. This is a generalisation of the situation in Chapko and Johansson (2008) to three-dimensions and with Cauchy data only partially given.… More >

K. Ibrahim¹, W.A. El-Askary^1,2, A. Balabel¹, I.M. Sakr¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 79-104, 2012, DOI:10.3970/cmes.2012.085.079

Abstract In the present paper, a numerical code has been developed with different turbulence models aiming at simulating turbulent bubbly flows in vertical circular pipes. The mass and momentum conservation equations are used to describe the motion of both phases (water/air). Because of the averaging process additional models are needed for the inter-phase momentum transfer and turbulence quantities for closure. The continuous phase (water) turbulence is represented using different turbulence models namely: two-equation k-ε, extended k-ε and shear-stress transport (SST) k-ω turbulence models which contains additional term to account for the effect of the dispersed phase (air) on the continuous phase… More >

Miaoxia Xie¹, Yueming Li¹, Hualing Chen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 65-78, 2012, DOI:10.3970/cmes.2012.085.065

Abstract Energy finite element method(EFEM) is a promising method to solve high-frequency vibro-acoustic problem. Energy boundary element method (EBEM) is an effective way to compute high-frequency sound radiation in the unbounded medium. Vibro-acoustic coupling of cavity structure in anechoic chamber includes both the interior acoustic field and unbounded exterior acoustic field. In order to predict this kind of high-frequency vibro-acoustic coupling problem in anechoic chamber, an approach combined EFEM and EBEM is developed in this paper. As a numerical example, the approach is applied to solve the high-frequency vibro-acoustic coupling response of a cubic cavity structure excited by a point sound… More >

Ahmad Shirzadi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 45-64, 2012, DOI:10.3970/cmes.2012.085.045

Abstract This paper presents a meshless method based on the meshless local integral equation (LIE) method for solving the two-dimensional diffusion and diffusion-convection equations subject to a non-local condition. Suitable finite difference scheme is used to eliminate the time dependence of the problem. A weak formulation on local subdomains with employing the fundamental solution of the Laplace equation as test function transforms the resultant elliptic type equations into local integral equations. Then, the Moving Least Squares (MLS) approximation is employed for discretizing spatial variables. Two illustrative examples with exact solutions being used as benchmark solutions are presented to show the efficiency… More >