• Structural Engineering and Mechanics
• /
• v.55 no.3
• /
• pp.473-494
• /
• 2015

A graded harmonic finite element formulation based on three-dimensional elasticity theory is developed for the structural analysis of 2D functionally graded axisymmetric structures. The mechanical properties of the axisymmetric solid structures composed of two different metals and ceramics are assumed to vary in radial and axial directions according to power law variations as a function of the volume fractions of the constituents. The material properties of the graded element are calculated at the integration points. Effects of material distribution profile on the static deformation, natural frequency and dynamic response analyses of particular axisymmetric solid structures are investigated by changing the power law exponents. It is observed that the displacements, stresses and natural frequencies are severely affected by the variation of axial and radial power law exponents. Good accuracy is obtained with fewer elements in the present study since Fourier series expansion eliminates the need of finite element mesh in circumferential direction and continuous material property distribution within the elements improves accuracy without refining the mesh size in axial and radial directions.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.38-42
• /
• 2000

In this study, the low velocity impact behavior of the composite laminates has been described by using 3 dimensional nonlinear finite elements. To describe the geometric nonlinearity due to large deformation, the dynamic contact problem is formulated using the exterior penalty finite element method on the base of Total Lagrangian formulation. The incremental decomposition is introduced, and the converged solution is attained by Newton-Raphson Method. The Newmark's constant-acceleration time integration algorithm is used. To make verification of the finite element program developed in this study, the solution of the nonlinear static problem with occurrence of large deformation is compared with ABAQUS, and the solution of the static contact problem with indentation is compared with the Hertz solution. And, the solution of low velocity impact problem for isotropic material is verificated by comparison with that of LS-DYNA3D. Finally the contact force of impact response from the nonlinear analysis are compared with those from the linear analysis.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.2410-2413
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAB, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. In the pre-processor of 3D FEM package, a new 3D mesh generating algorithm can make information on 3D Delaunay tetrahedral mesh elements for analyses of 3D FEM problems. The solver of the 3D FEM package offers three methods for solving the linear algebraic matrix equation, i.e., Gauss-Jordan elimination solver, Band solver, and Skyline solver. The post-processor visualizes the results for 3D FEM problems such as the deformed position and the stress. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.

• Proceedings of the KIEE Conference
• /
• 2001.04a
• /
• pp.167-170
• /
• 2001

A three dimensional adaptive finite element analysis algorithm is developed. In the method, the edge elements are used for field analysis, and the local error. In each element is estimated from the fact that the magnetic field should satisfy. The continuity condition at the interface of the two adjacent elements. Based on the estimated error, the elements which are considered to have big error are divided into several elements using the bisection method. The effectiveness of the developed algorithm is proved through numerical examples.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.700-705
• /
• 2004

Shell finite elements are widely used for the analysis of thin section objects such as sheet metal parts, automobile bodies and et al. due to their computational efficiency. Since many of input data for finite element analysis are given as solid models or triangulated surface models, one should extract midsurface information from these input data initially and then construct shell meshes on the extracted midsurfaces. In this paper, a method of generating shell elements on midsurfaces directly from input models have been proposed. In order to construct shell meshes, the input models should be triangulated on surfaces first, and then tetrahedral elements are generated by using an advancing front method, and finally mid shell surfaces are obtained from tetrahedral meshes. Some examples are given to demonstrate the efficiency of the proposed method.

• Smart Structures and Systems
• /
• v.15 no.1
• /
• pp.99-117
• /
• 2015

The Stonecutters Bridge (SCB) in Hong Kong is the third-longest cable-stayed bridge in the world with a main span stretching 1,018 m between two 298 m high single-leg tapering composite towers. A Wind and Structural Health Monitoring System (WASHMS) is being implemented on SCB by the Highways Department of The Hong Kong SAR Government, and the SCB-WASHMS is composed of more than 1,300 sensors in 15 types. In order to establish a linkage between structural health monitoring and maintenance management, a Structural Health Rating System (SHRS) with relevant rating tools and indices is devised. On the basis of a 3D space frame finite element model (FEM) of SCB and model updating, this paper presents the development of an SHR-oriented 3D multi-scale FEM for the purpose of load-resistance analysis and damage evaluation in structural element level, including modeling, refinement and validation of the multi-scale FEM. The refined 3D structural segments at deck and towers are established in critical segment positions corresponding to maximum cable forces. The components in the critical segment region are modeled as a full 3D FEM and fitted into the 3D space frame FEM. The boundary conditions between beam and shell elements are performed conforming to equivalent stiffness, effective mass and compatibility of deformation. The 3D multi-scale FEM is verified by the in-situ measured dynamic characteristics and static response. A good agreement between the FEM and measurement results indicates that the 3D multi-scale FEM is precise and efficient for WASHMS and SHRS of SCB. In addition, stress distribution and concentration of the critical segments in the 3D multi-scale FEM under temperature loads, static wind loads and equivalent seismic loads are investigated. Stress concentration elements under equivalent seismic loads exist in the anchor zone in steel/concrete beam and the anchor plate edge in steel anchor box of the towers.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.1067-1089
• /
• 2014

In this paper, three dimensional static and dynamic analyses of two dimensional functionally graded annular sector plates have been investigated. The material properties vary through both the radial and axial directions continuously. Graded finite element and Newmark direct integration methods have been used to solve the 3D-elasticity equations in time and space domains. The effects of power law exponents and different boundary conditions on the behavior of FGM annular sector plate have been investigated. Results show that using 2D-FGMs and graded elements have superiority over the homogenous elements and 1D-FGMs. The model has been compared with the result of a 1D-FGM annular sector plate and it shows good agreement.

• Proceedings of the KWS Conference
• /
• 2003.05a
• /
• pp.269-271
• /
• 2003

A simplified analysis method is proposed in this paper, through which the deformation due to spot welding is obtained. In order to analyze the spot welding process, the finite elements are axisymmetrically modeled, and the results from them are used for analyzing the deformation of spot welded structure. It takes more little time in this simplified analysis than the analysis by 3-D continuum elements, and this analysis produces the results that nearly agree with the measured data.

• Journal of the Chungcheong Mathematical Society
• /
• v.35 no.3
• /
• pp.227-234
• /
• 2022

Let 𝔽_q be a finite field with odd q elements. In this article, we prove that if E ⊆ 𝔽^d_q, d ≥ 2, and |E| ≥ q, then there exists a set Y ⊆ 𝔽^d_q with |Y| ~ q^d such that for all y ∈ Y, the number of distances between the point y and the set E is ~ q. As a corollary, we obtain that for each set E ⊆ 𝔽^d_q with |E| ≥ q, there exists a set Y ⊆ 𝔽^d_q with |Y| ~ q^d so that any set E ∪ {y} with y ∈ Y determines a positive proportion of all possible distances. The averaging argument and the pigeonhole principle play a crucial role in proving our results.

• Geomechanics and Engineering
• /
• v.12 no.6
• /
• pp.1021-1046
• /
• 2017

The assessment of slope stability is an essential task in geotechnical engineering. In this paper, a three-dimensional (3D) finite element analysis (FEA) was employed to investigate the performance of different shear pin arrangements to increase the stability of a soil block resting on an inclined plane with a low-interface friction plane. In the numerical models, the soil block was modeled by volume elements with linear elastic perfectly plastic material in a drained condition, while the shear pins were modeled by volume elements with linear elastic material. Interface elements were used along the bedding plane (bedding interface element) and around the shear pins (shear pin interface element) to simulate the soil-structure interaction. Bedding interface elements were used to capture the shear sliding of the soil on the low-interface friction plane while shear pin interface elements were used to model the shear bonding of the soil around the pins. A failure analysis was performed by means of the gravity loading method. The results of the 3D FEA with the numerical models were compared to those with the physical models for all cases. The effects of the number of shear pins, the shear pin locations, the different shear pin arrangements, the thickness and the width of the soil block and the associated failure mechanisms were discussed.