• Journal of Navigation and Port Research
• /
• v.34 no.5
• /
• pp.367-373
• /
• 2010

When the numerical analysis is carried out, it is necessary to set proper elements as a feature of analysis domains for more accurate simulations. In this study, Distinct Element Method(DEM) is applied, only considering repulsive force and tensile force except for frictional force and resisting force of particle. When the filled particles with initial Quad-tree type is relocated by DEM, a blank space existing among the particles can be minimized because the shape of particle is circular. Finally, it is the effective feature that the centroidal disposion of the particles is similar to an equilateral triangle. Triangular mesh are formed by using the Delaunay triangular technique on these relocated particles, the quality of triangular mesh is more improved by carrying out Laplace interpolations. The compared result of Aspect Ratio before and after the Laplace interpolation is shown that although the quality of triangular mesh made by DEM is good, the later triangular mesh are higher quality than the formers. In this study, although the developed technique takes a longer calculational time than the previous technique to generate triangular mesh, it is considered that the applicable possibility is very high in the generation of finite element mesh about wave analysis and various numerical simulation to need a complex or reappearance of exact topography.

• Structural Engineering and Mechanics
• /
• v.14 no.5
• /
• pp.575-593
• /
• 2002

Based on the Mindlin plate theory, a refined discrete 15-DOF triangular laminated composite plate finite element RDTMLC with the re-constitution of the shear strain is proposed. For constituting the element displacement function, the exact displacement function of the Timoshenko's laminated composite beam as the displacement on the element boundary is used to derive the element displacements. The proposed element can be used for the analysis of both moderately thick and thin laminated composite plate, and the convergence for the very thin situation can be ensured theoretically. Numerical examples presented show that the present model indeed possesses the properties of higher accuracy for anisotropic laminated composite plates and is free of locking even for extremely thin laminated plates.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.10
• /
• pp.1475-1482
• /
• 2004

In this paper, the rigid-viscoplastic finite element solutions obtained by MINI-elements based triangular elements and tetrahedral elements are compared with those obtained from numerically well-behaved rectangular and hexahedral elements. The theoretical background of the MINI-elements is introduced in detail and the rigid-viscoplastic finite element formulation is also given. Discussion on the results of the MINI-elements is made with emphasis on the effect of a stabilizer simplifying velocity-bubble coupled terms.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.297-320
• /
• 2005

In recent years the pure displacement formulation for plate elements has not been as popular as other formulations. We revisit the pure displacement formulation for shear-deformable plate elements and propose a family of N-node, displacement-compatible, fully-integrated, pure-displacement, triangular, Mindlin plate elements, MIN-N. The development has been motivated by the relative simplicity of the pure displacement formulation and by the success of the existing 3-node plate element, MIN3. The formulation of MIN3 is generalized to obtain the MIN-N family, which possesses complete, fully compatible kinematic fields, in which the interpolation functions for transverse displacement are one degree higher than those for rotations. General element-level formulas for the thin-limit Kirchhoff constraints are developed. The 6-node, 18 degree-of-freedom element MIN6, with cubic displacement and quadratic rotations, is implemented and tested extensively. Numerical results show that MIN6 exhibits good performance for both static and dynamic analyses in the linear, elastic regime. The results illustrate that the fully-integrated MIN6 element has excellent performance in the thin limit, even for coarse meshes, and that it does not require shear relaxation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.39-44
• /
• 2002

A three node triangular element with drilling rotations incorporating Improved Layerwise Zig-zag Theory(HZZT) is developed to analyze the vibration of spinning pretwisted composite blades with embedded damping layer. Matching conditions at the interfaces between the damping material and the border material are enforced by setting the shear forces matched and different shear strains along the interfaces. The natural frequencies and modal loss factors of cantilevered pretwisted composite blade with damping core are calculated with the present triangular element enforcing the matching conditions and compared to experimental results and MSC/NASTRAN results using a layered combination of plate and solid elements.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.4 s.10
• /
• pp.85-92
• /
• 2003

In this work, a finite element model is presented for geometrically non-linear analysis of shell structures. Finite element by using a three-node flat triangular shell element is formulated. The non-linear incremental equilibrium equations are formulated by using an updated Lagrangian formulation and the solutions are obtained with the incremental/iterative Newton-Raphson method and arc length method. Some of results are presented for shell structures. The obtained results are in good agreement with the results available in existing literature.

• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

• Structural Engineering and Mechanics
• /
• v.43 no.2
• /
• pp.253-271
• /
• 2012

To analyze the bending and transverse shear effects of laminated composite plates, a thirteen nodes triangular element will be presented. The suggested formulations consider a parabolic variation of the transverse shear strains through the thickness. As a result, there is no need to use shear correction coefficients in computing the shear stresses. The proposed element can model both thin and thick plates without any problems, such as shear locking and spurious modes. Moreover, the effectiveness of $w_{,n}$, as an independent degree of freedom, is concluded by the present study. To perform the accuracy tests, several examples will be solved. Numerical results for the orthotropic materials with different boundary conditions, shapes, number of layers, thickness ratios and fiber orientations will be presented. The suggested element calculates the deflections and stresses more accurate than those available in the literature.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.229-236
• /
• 2004

A new three-node triangular shell element based on higher order zig-zag theory is developed for laminated composite shells with multiple delaminations. The present higher order zig-zag shell theory is described in a general curvilinear coordinate system and in general tensor notation. All the complicated curvatures of surface including twisting curvatures can be described in an exact manner in the present shell element because this element is based on geometrically exact surface representation. The displacement field of the proposed finite element includes slope of deflection. which requires continuity between element interfaces. Thus the nonconforming shape function of Specht's three-node triangular plate bending element is employed to interpolate out-of-plane displacement. The present element passes the bending and twisting patch tests in flat surface configurations. The developed element is evaluated through the buckling problems of composite cylindrical shells with multiple delaminations. Through the numerical examples it is demonstrated that the proposed shell element is efficient because it has minimal degrees of freedom per node. The accuracy of the present element is demonstrated in the prediction of buckling loads and buckling modes of shells with multiple delaminations. The present shell element should serve as a powerful tool in the prediction of buckling loads and modes of multi-layered thick laminated shell structures with arbitrary-shaped multiple delaminations.