P. H. Wen¹, Y. C. Hon²

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.3, pp. 177-192, 2007, DOI:10.3970/cmes.2007.021.177

Abstract In this paper, we perform a geometrically nonlinear analysis of Reissner-Mindlin plate by using a meshless collocation method. The use of the smooth radial basis functions (RBFs) gives an advantage to evaluate higher order derivatives of the solution at no cost on extra-interpolation. The computational cost is low and requires neither the connectivity of mesh in the domain/boundary nor integrations of fundamental/particular solutions. The coupled nonlinear terms in the equilibrium equations for both the plane stress and plate bending problems are treated as body forces. Two load increment schemes are developed to solve the nonlinear differential equations. Numerical verifications are… More >

L.N. Gergidis, D. Kourounis, S. Mavratzas, A. Charalambopoulos¹

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.2, pp. 157-176, 2007, DOI:10.3970/cmes.2007.021.157

Abstract A new complete set of scattering eigensolutions of Helmholtz equation in spheroidal geometry is constructed in this paper. It is based on the extension to exterior boundary value problems of the well known Vekua transformation pair, which connects the kernels of Laplace and Helmholtz operators. The derivation of this set is purely analytic. It avoids the implication of the spheroidal wave functions along with their accompanying numerical deficiencies. Using this novel set of eigensolutions, we solve the acoustic scattering problem from a soft acoustic spheroidal scatterer, by expanding the scattered field in terms of it. Two approaches concerning the determination… More >

Shunde Yin¹, Leo Rothenburg¹, Maurice B. Dusseault¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 111-120, 2007, DOI:10.3970/cmes.2007.019.111

Abstract Subsidence problem is of great importance in petroleum engineering and environmental engineering. In this paper, we firstly apply a hybrid Displacement Discontinuity-FEM modeling to this classic problem: the evaluation of subsidence over a compacting oil reservoir. We use displacement discontinuity method to account for the reservoir surrounding area, and finite element methods in the fully coupled simulation of the reservoir itself. This approach greatly reduces the number of degrees of freedom compared to an analyzing fully coupled problem using only a finite element or finite difference discretization. More >

T. Jarak¹, J. Sorić¹, J. Hoster¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 235-246, 2007, DOI:10.3970/cmes.2007.018.235

Abstract A meshless computational method based on the local Petrov-Galerkin approach for the analysis of shell structures is presented. A concept of a three dimensional solid, allowing the use of completely 3-D constitutive models, is applied. Discretization is carried out by using both a moving least square approximation and polynomial functions. The exact shell geometry can be described. Thickness locking is eliminated by using a hierarchical quadratic approximation over the thickness. The shear locking phenomena in case of thin structures and the sensitivity to rigid body motions are minimized by applying interpolation functions of sufficiently high order. The numerical efficiency of… More >

S.Y. Wang^1,2, K.M. Lim^2,3, B.C. Khoo^2,3, M.Y. Wang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 155-182, 2007, DOI:10.3970/cmes.2007.018.155

Abstract In this paper, a geometric deformation constrained level set method is presented as an effective approach for structural shape and topology optimization. A level set method is used to capture the motion of the free boundary of a structure. Furthermore, the geometric deformation of the free boundary is constrained to preserve the structural connectivity and/or topology during the level set evolution. An image-processing-based structural connectivity and topology preserving approach is proposed. A connected components labeling technique based on the 4-neighborhood connectivity measure and a binary image is used for the present region identification. The corresponding binary image after an exploratory… More >

M. Heshmatzadeh, G. E. Bridges¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.1, pp. 45-58, 2007, DOI:10.3970/cmes.2007.018.045

Abstract Cell Method, a Finite Formulation technique, is compared in detail with the Finite Element Method (FEM), a differential-based numerical technique. In the finite formulation technique, Poisson's equation is described starting from a topological foundation. The final set of algebraic equations resulting from the two approaches are compared in matrix form. The equivalence of the coefficient matrices is proven for a Voronoi dual mesh and linear shape functions in the FEM. The difference between the source (charge) vectors in the two approaches is described. It is shown that the use of linear shape functions in the FEM is equivalent to the… More >

M. Kröger¹, M. Müller², J. Nievergelt²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 559-570, 2003, DOI:10.3970/cmes.2003.004.559

Abstract We present a novel method for generating starting polymer structures for molecular simulations in the dense phase. The work describes the ingredients of an algorithm for the creation of large, dense or diluted amorphous polymeric systems close to equilibrium and provides measures for its quality. The model systems are made of semiflexible (wormlike) repulsive multibead chains. The key feature of the method is its efficiency, in particular for large systems, while approaching given local and global chain characteristics. Its output has been proven to serve as an excellent basis for subsequent off-lattice molecular dynamics computer simulation. By combining chain growing… More >

Wen Yi Lin¹, Kuo Mo Hsiao²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 259-272, 2003, DOI:10.3970/cmes.2003.004.259

Abstract The buckling and postbuckling behavior of spatial rods under different types of end torque and compressive axial force is investigated using finite element method. All coupling among bending, twisting, and stretching deformations for beam element is considered by consistent second-order linearization of the fully geometrically nonlinear beam theory. However, the third order term of the twist rate is also considered. An incremental-iterative method based on the Newton-Raphson method combined with constant arc length of incremental displacement vector is employed for the solution of nonlinear equilibrium equations. The zero value of the tangent stiffness matrix determinant of the structure is used… More >

K. Ijima¹, H. Obiya¹, S. Iguchi², S. Goto²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 239-248, 2003, DOI:10.3970/cmes.2003.004.239

Abstract Defining element coordinates in space frame, element end deformations become statically clear from the energy principle. Therefore, the deformations can be expressed by nodal displacement without any approximation. The paper indicates that the exact expressions of the deformations and the geometrical stiffness strictly based on the equations makes large displacement analysis of space frame possible with robustness on the computation. More >

T. Chekifi^1,2,*, B. Dennai², R. Khelfaoui²

FDMP-Fluid Dynamics & Materials Processing, Vol.14, No.3, pp. 201-212, 2018, DOI: 10.3970/fdmp.2018.00322

Abstract The fluidic oscillator is an interesting device developed for passive flow measurement. These microsystems can produce a high oscillating jet frequency with high flow velocity. The main advantages of fluidic oscillators are that no moving parts is included in the device. Commercial CFD code FLUENT was used to perform analysis of flows in fluidic oscillator. Numerical simulations were carried out for different flow conditions, where water and air were used as working fluids. The oscillation frequencies were identified by the discrete fast Fourier transform method (FFT). Furthermore a low-pressure vortex of fluid flow in the oscillating chamber was observed. The… More >