• Journal of Welding and Joining
• /
• 제21권6호
• /
• pp.71-76
• /
• 2003

For the purpose of clarifying the influence of welding residual stress to the fatigue crack propagations behaviour, an analytical investigation based on finite element method is performed to examine the opening behaviour of tip-closed crack in the compressive residual stress. A finite element model comprised of contact elements for the crack plane and plane stress elements for the base material is used to evaluate crack opening stress of the crack existing in the residual stress field. Also an analytical method based on the superposition principle to estimate the length of opened part of tip closed crack and the stress distribution adjacent to the crack during uploading is applied to the finite element model. The software for the analysis is ABAQUS, which is a general purpose finite element package. The results show that stresses distributed on the crack surfaces are reduced and approached to zero as the applied stresses are increased up to crack tip opening stress and no mechanical discontinuity is found at the boundary of contact elements and plane stress elements. It is verified that the opening behavior of the fatigue crack in the residual stress can be predicted by finite element method with the proposed analytical method.

• Structural Engineering and Mechanics
• /
• 제30권1호
• /
• pp.119-129
• /
• 2008

A small strain and elastoplastic formulation of Polygonal Element Method (PEM) is developed for efficient analysis of elastoplastic solids. In this work, the polygonal elements are constructed based on traditional triangular finite meshes. The construction method of polygonal mesh can directly utilize the sophisticated triangularization algorithm and reduce the difficulty in generating polygonal elements. The Wachspress rational finite element basis function is used to construct the approximations of polygonal elements. The incremental variational form and a von Mises type model are used for non-linear elastoplastic analysis. Several small strain elastoplastic numerical examples are presented to verify the advantages and the accuracy of the numerical formulation.

• Proceeding of EDISON Challenge
• /
• 제5회(2016년)
• /
• pp.195-202
• /
• 2016

In this paper, crack detection method using eigen value analysis is presented. Three methods are used: theoretical analysis, finite element method with the cracked beam elements and finite element method with three dimensional continuum elements. Finite element formulation of the cracked beam element is introduced. Additional term about stress intensity factor based on fracture mechanics theory is added to flexibility matrix of original beam to model the crack. As using calculated stiffness matrix of cracked beam element and mass matrix, natural frequencies are calculated by eigen value analysis. In the case of using continuum elements, the natural frequencies could be calculated by using EDISON CASAD solver. Several cases of crack are simulated to obtain natural frequencies corresponding the crack. The surface of natural frequency is plotted as changing with crack location and depth. Inverse analysis method is used to find crack location and depth from the natural frequencies of experimental data, which are referred by another papers. Predicted results are similar with the true crack location and depth.

• Journal of KSNVE
• /
• 제10권1호
• /
• pp.97-106
• /
• 2000

The dynamic behavior of crankshaft-bearing system in scroll compressor has been investigated using the combined methodologies of finite elements and transfer matrices. The finite element formulation is proposed including the field element for a shaft section and the point element at balancer weight locations, bearing locations, etc., whereas the conventional method is used with the elements. The Houbolt method is used to consider the time march for the integration of the system equations. The linear stiffness and damping coefficients are calculated for a finite cylindrical fluid-film bearing by solving the Reynolds equation, using finite difference method. The orbital response of crankshaft supported on the linear bearing model is obtained, considering balancer weights of motor rotor. And, the steady state displacement of crankshaft are compared with a variation in balancer weight. The loci of crankshaft at bearing locations are composed of the synchronous whirl component and the non-synchronous whirl component.

• Journal of Mechanical Science and Technology
• /
• 제15권11호
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

• Korean Journal of Mathematics
• /
• 제13권2호
• /
• pp.127-137
• /
• 2005

The objective of numerical analysis is to devise and analyze efficient algorithms or numerical methods for equations arising in mathematical modeling for science and engineering. In this article, we present some recent topics in computational mathematics, specially in the finite element method and overview the development of the mixed finite element method in the context of second order elliptic and parabolic problems. Multiscale methods such as MsFEM, HMM, and VMsM are included.

• Coupled systems mechanics
• /
• 제2권4호
• /
• pp.389-410
• /
• 2013

Three-dimensional Lagrangian fluid finite element is applied to seismic response analysis of an oil storage tank with a floating roof. The fluid element utilized in the present analysis is formulated based on the displacement finite element method considering only volumetric elasticity and its element stiffness matrix is derived by using one-point integration method in order to avoid volumetric locking. The method usually adds a rotational penalty stiffness to satisfy the irrotational condition for fluid motion and modifies element mass matrices through the projected mass method to suppress spurious hourglass-mode appeared in compensation for one-point integration. In the fluid element utilized in the present paper, a small hourglass stiffness is employed. The fluid and structure domains for the objective oil storage tank are modeled by eight-node solid elements and four-node shell elements, respectively, and the transient response of the floating roof structure or the free surface are evaluated by implicit direct time integration method. The results of seismic response analyses are compared with those by other method and the validation of the present analysis using three-dimensional Lagrangian fluid finite elements is shown.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• 제28권10호
• /
• pp.1517-1525
• /
• 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 has been proposed, in which midsurface generating process can be omitted. 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.

• Steel and Composite Structures
• /
• 제6권3호
• /
• pp.183-202
• /
• 2006

The behavior of steel-concrete composite beams is strongly influenced by the type of shear connection between the steel beam and the concrete slab. For accurate analytical predictions, the structural model must account for the interlayer slip between these two components. This paper focuses on a procedure for response sensitivity analysis using state-of-the-art finite elements for composite beams with deformable shear connection. Monotonic and cyclic loading cases are considered. Realistic cyclic uniaxial constitutive laws are adopted for the steel and concrete materials as well as for the shear connection. The finite element response sensitivity analysis is performed according to the Direct Differentiation Method (DDM); its analytical derivation and computer implementation are validated through Forward Finite Difference (FFD) analysis. Sensitivity analysis results are used to gain insight into the effect and relative importance of the various material parameters in regards to the nonlinear monotonic and cyclic response of continuous composite beams, which are commonly used in bridge construction.

• Korean Journal of Computational Design and Engineering
• /
• 제12권3호
• /
• pp.220-232
• /
• 2007

Current geometric modeling and analysis are commonly based on B-Rep modeling and a finite elements method respectively. Furthermore, it is difficult to represent an object whose material property is heterogeneous using the B-Rep method because the B-Rep is basically used for homogeneous models. In addition, meshes are required to analyze a property of a model when the finite elements method is applied. However, the process of generating meshes from B-Rep is cumbersome and sometimes difficult especially when the model is deformed as time goes by because the topology of deforming meshes are changed. To overcome those problems in modeling and analysis including homogeneous and heterogeneous materials, we suggest a unified modeling and analysis method based on implicit representation of the model using R-function which is suggested by Rvachev. For implicit modeling of an object a distance field is approximated and blended for a complex object. Using the implicit function mesh-free analysis is possible where meshes are not necessary. Generally mesh-free analysis requires heavy computational cost compared to a finite elements method. To improve the computing time of function evaluation, we utilize GPU programming. Finally, we give an example of a simple pipe design problem and show modeling and analysis process using our unified modeling and analysis method.