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 >

Hayder A. Abdulbari ^1,2, Hassan D. Mahammed¹, Z. Hassan, Wafaa K. Mahmood³

FDMP-Fluid Dynamics & Materials Processing, Vol.12, No.2, pp. 86-101, 2016, DOI:10.3970/fdmp.2016.012.086

Abstract An experimental investigation was carried out to study the response of a turbulent pressure drop fluctuations to longitudinal groove riblets, involved two configurations being triangular and spaced triangular grooves with height 600, 800, 1000 μm and peak to peak spacing 1000 μm and 2000 μm respectively. Experiments were therefore performed at free stream velocity up to 0.44 m/sec, which were corresponding to Reynolds number (Re) 53000. The development of the obtained turbulent layer downstream of the grooves was then compared with the results from the corresponding smooth-wall case. To conclude, the effect of the spaced triangular riblets on the turbulent… More >

A. Filali¹, L. Khezzar^1,2

FDMP-Fluid Dynamics & Materials Processing, Vol.9, No.3, pp. 307-329, 2013, DOI:10.3970/fdmp.2013.009.307

Abstract Numerical investigations of purely-elastic instabilities occurring in creeping flows are reported in planar cross-slot geometries with both sharp and round corners. The fluid is described by the upper-convected Maxwell model, and the governing equations are solved using the finite element technique based on a steady (non-iterative) direct solver implemented in the POLYFLOWcommercial software (version 14.0). Specifically, extensive simulations were carried out on different meshes, with and without the use of flow perturbations, for a wide range of rheological parameters. Such simulations show the onset of flow asymmetries above a critical Deborah number (De). The effect of rounding the corners is… More >

Priyankar Datta¹, Manas C. Ray¹

CMC-Computers, Materials & Continua, Vol.49-50, No.1, pp. 47-80, 2015, DOI:10.3970/cmc.2015.049.047

Abstract This paper deals with the implementation of the one dimensional form of the fractional order derivative constitutive relation for three dimensional analysis of active constrained layer damping (ACLD) of geometrically nonlinear laminated composite plates. The constraining layer of the ACLD treatment is composed of the vertically/obliquely reinforced 1–3 piezoelectric composites (PZCs). The von Kármán type nonlinear strain displacement relations are used to account for the geometric nonlinearity of the plates. A nonlinear smart finite element model (FEM) has been developed. Thin laminated substrate composite plates with various boundary conditions and stacking sequences are analyzed to verify the effectiveness of the… More >