• Journal of Ocean Engineering and Technology
• /
• v.23 no.5
• /
• pp.1-8
• /
• 2009

This paper is a continuation of the recent development on Hermite-based divergence free element method and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square enclosure with temperature difference across the two sides. The basis functions for the velocity field consist of the Hermite function and its curl while the basis functions for the temperature field consists of the Hermite function and its gradients. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x and y derivatives. This paper presents numerical results for Ra = 105, and compares with those from a stabilized finite element method developed by Illinca et al. (2000). The comparison has been done on 32 by 32 uniform elements and the degree of approximation of elements used for the stabilized finite element are linear (Deg. 1) and quadratic (Deg. 2). The numerical results from both methods show well agreements with those of De vahl Davi (1983).

• Structural Engineering and Mechanics
• /
• v.7 no.1
• /
• pp.111-126
• /
• 1999

The work presented here focuses on the development of suitable discretised formulations, for large-displacement shape and non-shape design sensitivity analysis (DSA), which enable the straightforward incorporation of structural optimisation into established finite element analysis (FEA) codes. For the generalised displacement-based functional the design sensitivity vector has been expressed in terms of displacement sensitivity. The Total Lagrangian formulation is utilised for modelling of large deformation of truss structures. The variational formulation of the sensitivity analysis procedure is discretised by using "pseudo" - finite elements, Results are presented for the sensitivity analysis and optimisation of standard truss structures. For the purposes of this work, the analysis and optimisation procedures outlined below are incorporated into the FEA code ABAQUS.

• Journal of the Korean Society for Precision Engineering
• /
• v.19 no.3
• /
• pp.46-51
• /
• 2002

Recently, High capacity alternator are used fur some special equipments in industry. But, several serious problem are occured, especially, broken coil of rotor, caused by crash with stator on rotating the rotor. Although added coil cover thor protect coil of rotor, coil cover is broken. In this study purpose 2 step for corrected that problems. First, three dimensional finite element method far investigate what is most important point. For that purpose, performed stress analysis of coil and coil cover that modeling and finite analysis by ANSYS software. Second, Apply prestress when winding the coil on stator to modify direction of net force. Vector analysis is used for determine corrected prestress. Result of the analysis and prestress are reviewed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.11
• /
• pp.1730-1736
• /
• 2001

A finite element analysis for a rotating cantilever beam is presented in this study. Based on a dynamic modeling method using the stretch deformation instead of the conventional axial deformation, three linear partial differential equations are (derived from Hamilton's principle. Two of the linear differential equations show the coupling effect between stretch and chordwise deformations. The other equation is an uncoupled one for the flapwise deformation. From these partial differential equations and the associated boundary conditions, two weak forms are derived: one is for the chordwise motion and the other is fur the flptwise motion. The weak farms are spatially discretized with newly defined two-node beam elements. With the discretized equations or the matrix-vector equations, the behaviors of the natural frequencies are investigated for the variation of the rotating speed.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.147-161
• /
• 2023

Based on the ideas of variational differential quadrature (VDQ) and finite element method (FEM), a numerical approach named as VDQFEM is applied herein to study the large deformations of plate-type structures under static loading with arbitrary shape hole made of functionally graded graphene platelet-reinforced composite (FG-GPLRC) in the context of higher-order shear deformation theory (HSDT). The material properties of composite are approximated based upon the modified Halpin-Tsai model and rule of mixture. Furthermore, various FG distribution patterns are considered along the thickness direction of plate for GPLs. Using novel vector/matrix relations, the governing equations are derived through a variational approach. The matricized formulation can be efficiently employed in the coding process of numerical methods. In VDQFEM, the space domain of structure is first transformed into a number of finite elements. Then, the VDQ discretization technique is implemented within each element. As the last step, the assemblage procedure is performed to derive the set of governing equations which is solved via the pseudo arc-length continuation algorithm. Also, since HSDT is used herein, the mixed formulation approach is proposed to accommodate the continuity of first-order derivatives on the common boundaries of elements. Rectangular and circular plates under various boundary conditions with circular/rectangular/elliptical cutout are selected to generate the numerical results. In the numerical examples, the effects of geometrical properties and reinforcement with GPL on the nonlinear maximum deflection-transverse load amplitude curve are studied.

• Computers and Concrete
• /
• v.2 no.5
• /
• pp.345-366
• /
• 2005

This work presents an assessment of the computational performance of a vector-parallel implementation of probabilistic model for concrete cracking in 3D. This paper shows the continuing efforts towards code optimization as reported in earlier works Paz, et al. (2002a,b and 2003). The probabilistic crack approach is based on the direct Monte Carlo method. Cracking is accounted by means of 3D interface elements. This approach considers that all nonlinearities are restricted to interface elements modeling cracks. The heterogeneity governs the overall cracking behavior and related size effects on concrete fracture. Computational kernels in the implementation are the inexact Newton iterative driver to solve the non-linear problem and a preconditioned conjugate gradient (PCG) driver to solve linearized equations, using an element by element (EBE) strategy to compute matrix-vector products. In particular the paper analyzes code behavior using OpenMP directives in parallel vector processors (PVP), such as the CRAY SV1 and CRAY T94. The impact of the memory architecture on code performance, and also some strategies devised to circumvent this issue are addressed by numerical experiment.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.275-284
• /
• 2007

For precise interpretation of magnetotelluric (MT) data distorted by irregular surface terrain, topography effects are investigated by computing apparent resistivities, phases, tippers and induction vectors for a three-dimensional (3D) hill-and-valley model. To compute MT responses for the 3D surface topography model, we use a 3D MT modeling algorithm based on an edge finite-element method which is free from vector parasites. Distortions on the apparent resistivity and phase are mainly caused by distorted currents that flow along surface topography. The distribution of tipper amplitudes over both hill and valley are the same, while the tipper points toward the center of hill and the base of the valley. The real part of induction vector also points in the same direction as that of tipper, while the imaginary part in the opposite direction.

• Journal of Environmental Science International
• /
• v.16 no.10
• /
• pp.1169-1177
• /
• 2007

The main objective of this research was to evaluate the effects of process variables which were forming ability, flow displacement, effective stress, effective strain, fluid vector and products defects on manufactured artificial lightweight aggregate made of both incinerated sewage sludge ash and clay by means of the numerical analysis of a rigid-plastic finite element method. CATIA (3D CAD program) was used for an extrusion metal mold design that was widely used in designing aircraft, automobile and metallic molds. A metal forming analysis program (ATES Co.) had a function of a rigid-plastic finite element method was used to analyze the program. The result of extrusion forming analysis indicated clearly that a shape retention of the manufactured artificial light-weight aggregate could be maintained by increasing the extrusion ratio (increasing compressive strength inside of extrusion die) and decreasing the die angle. The stress concentration of metal mold was increased by increasing an extrusion ratio, and it was higher in a junction of punch and materials, friction parts between a bottom of the punch and inside of a container, a place of die angle and a place of die of metal mold. Therefore, a heat treatment as well as a rounding treatment for stress distribution in the higher stress concentration regions were necessary to extend a lifetime of the metallic mold. A deformity of the products could have made from several factors which were a surface crack, a lack of the shape retention and a crack of inside of the products. Specially, the surface crack in the products was the most notably affected by the extrusion ratio.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• Proceedings of the Korean Society of Computer Information Conference
• /
• 2011.01a
• /
• pp.185-188
• /
• 2011

In order to estimate spatial charging process in the corona charging which has been used to make polymer electret, the electrical properties of polypropylene film were obtained from Thermally Stimulated Current(TSC) measurements after corona charging between knife electrode and cylinder electrode with the voltages of -5, -6, -7 and -8[kV], respectively. And then the electrostatic contour and the electric field vector were also simulated by using Finite Element Method(FEM). The edge effect around edge of knife electrode affected the electrostatic contour on surface of specimen and the electric field concentration inside specimen. The uneven charging state in the electret due to the mistake on design could be calculated and so the optimal design of corona charging device which is appropriate to various materials is come to be practicable.