• Journal of Mechanical Science and Technology
• /
• v.14 no.4
• /
• pp.401-407
• /
• 2000

An elastic-plastic finite element analysis is performed to investigate detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The finite element analysis performed under plane stress using 4-node isoparametric elements can predict fatigue crack closure behavior. The mesh of constant element size along crack surface can not predict the opening level of fatigue crack. The crack opening level for the constant mesh size increases linearly from initial crack growth. The crack opening level for variable mesh size, is almost flat after crack tip has passed the monotonic plastic zone. The prediction of crack opening level using the variable mesh size proportioning the reversed plastic zone size with the opening stress intensity factors presents a good agreement with the experimental data regardless of stress ratios.

• Structural Engineering and Mechanics
• /
• v.3 no.1
• /
• pp.75-89
• /
• 1995

A transition membrane element designated as CLM which has variable mid-side nodes with drilling freedoms has been presented in this paper. The functional for the linear problem, in which the drilling rotations are introduced as independent variables, has been formulated. The transition elements with variable side nodes can be efficiently used in the local mesh refinement for the in-plane structures, which have stress concentrations. A modified Gaussian quadrature is needed to be adopted to evaluate the stiffness matrices of these transition elements mainly due to the slope discontinuity of displacement within the elements. Detailed numerical studies show the excellent performance of the new transition elements developed in this study.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.265-272
• /
• 2007

The finite element linear buckling analysis of folded plate structures using adaptive h-refinement methods is presented in this paper. The variable-node flat shell element used in this study possesses the drilling D.O.F. which, in addition to improvement of the element behavior, permits an easy connection to other elements with six degrees of freedom per node. The Box-typed structures can be analyzed using these developed flat shell elements. By introducing the variable-node elements some difficulties associated with connecting the different layer patterns, which are common in the adaptive h-refinement on quadrilateral mesh, can be overcome. To obtain better stress field for the error estimation, the super-convergent patch recovery is used. The convergent buckling modes and the critical loads associated with these modes can be obtained.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.463-470
• /
• 1998

The boundary/finite element analysis for the seismic wave amplifications due to nonhomogeneous alluvial deposits was performed in this study. For numerical analysis, the homogeneous linear elastic soil half-space was modeled by using the 3-node isoparametric boundary elements and the inhomogeneous alluvial deposit was modeled by using the 8-node isoparametric finite elements. The two elements at interface were coupled together by the equilibrium condition of the tractions and the compatibility condition of the displacements. As a prarmetric variable, the incident angle and the dimensionsless frequency of the SH, P and SV-waves and the shear wave velocity ratio and the mass density ratio between the half-space and the alluvial deposit were selected.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.1A
• /
• pp.21-33
• /
• 2006

The reliability analysis can be conducted more effectively by formulating the stochastic finite element method suitable for the reliability theory about the cable stayed bridge. After conducting the initial equilibrium analysis of the cable stayed bridge, the program which can conduct the linear and nonlinear stochastic finite element analysis using the perturbation method and the reliability analysis considered to the correlation of the random variable is developed. Using the results of this program about the cable stayed bridge, the characteristic of the node displacement, element force and cable tension according to the correlation of the random variable is investigated quantitatively. Also the reliability index and the failure probability are examined by the compounding the correlation of the random variable.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.53 no.2
• /
• pp.63-68
• /
• 2004

This paper presents an approach to optimum design of Moving Magnet Type Linear Oscillatory Actuator(MM-LOA). The Finite Element Method is applied to characteristic parameters for characteristic analysis and in order to reduce modeling time and efforts, the moving model node technique is used. In addition the neural network is used to reduce computational time of analysis according to changing design variable. To confirm the validity of this study, optimum design results are compared with results of analysis procedure that is verified by experiment.

• Steel and Composite Structures
• /
• v.43 no.1
• /
• pp.1-17
• /
• 2022

The free and forced nonlinear dynamic behaviors of Porous Functionally Graded Material (PFGM) plates are examined by means of a High-Order Implicit Algorithm (HOIA). The formulation is developed using the Third-order Shear Deformation Theory (TSDT). Unlike previous works, the formulation is written without resorting to any homogenization technique neither rule of mixture nor considering FGM as a laminated composite, and the distribution of the porosity is assumed to be gradually variable through the thickness of the PFGM plates. Using the Hamilton principle, we establish the governing equations of motion. The Finite Element Method (FEM) is used to compute approximations of the resulting equations; FEM is adopted using a four-node quadrilateral finite element with seven Degrees Of Freedom (DOF) per node. Nonlinear equations are solved by a HOIA. The accuracy and the performance of the proposed approach are verified by presenting comparisons with literature results for vibration natural frequencies and dynamic response of PFGM plates under external loading. The influences of porosity volume fraction, porosity distribution, slenderness ratio and other parameters on the vibrations of PFGM plate are explored. The results demonstrate the significant impact of different physical and geometrical parameters on the vibration behavior of the PFGM plate.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.412-417
• /
• 2008

In a traditional topology optimization method, material properties are usually distributed by finite element density and visualized by a gray level image. The distribution method based on element density is adequate for a great mass of 2-D topology optimization problems. However, when it is used for 3-D topology optimization, it is always difficult to obtain a smooth model representation, and easily appears a virtualconnect phenomenon especially in a low-density domain. The 3-D structural topology optimization method has been developed using the node density instead of the element density that is based on SIMP (solid isotropic microstructure with penalization) algorithm. A computer code based on Matlab was written to validate the proposed method. When it was compared to the element density as design variable, this method could get a more uniform density distribution. To show the usefulness of this method, several typical examples of structure topology optimization are presented.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.146-148
• /
• 1998

For the shape design of electromagnetic devices using the FEM, the choice of design parameters influence to the success of the optimization process. If the design parameter distribution has a one to one corespondence with finite element model, we can encounter not only serious accuracy problem but also obtain a zigzag shape along the interface. The nodes between those design parameters can be parameterized by interpolating using one among many interpolation methods. The conventional parameterization of design parameters has a limit of application for shape, because design parameters and movable nodes are linearly intepolated. In this paper, using the B-spline curve that use to present any interfaces in computer graphics, the curvilinear parameterization between design parameters and node points is compared with the linear parameterization.