• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.367-371
• /
• 2001

In the shape optimization based on the finite element method, the accuracy of finite element analysis of a given structure is important to determine the final shape. In case of a bending dominant problem, finite element solutions by the full integration scheme are not reliable because of the locking phenomenon. Furthermore, in the process of shape optimization, the mesh distortion is large due to the change of the structure outline: therefore, we cannot guarantee the accurate result unless the finite element itself is accurate. We approach to more accurate shape optimization to diminish these inaccuracies by improving the existing finite element. The shape optimization using the modified finite element is applied to a two-dimensional simple beam. Results show that the modified finite element have improved the optimization results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.99-104
• /
• 2006

A semi-Lagrangian finite element scheme with objective time stepping algorithm for solving viscoelastic flow problem is presented. The convection terms in the momentum and constitutive equations are treated using a quasi-monotone semi-Lagrangian scheme, in which characteristic feet on a regular grid are traced backwards over a single time-step. Concerned with the generalized midpoint rule type of algorithms formulated to exactly preserve objectivity, we use the geometric transformation such as pull-back, push-forward operation. The method is applied to the 4:1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

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

This paper is concerned with development of an enhanced quadratic Mindlin plate bending element. The behavior of the proposed plate element is further improved by the coupled use of non-conforming displacement modes, the selectively reduced integration scheme, and the assumed shear strain fields. The improvement may be attributable to the fact that the merits of these improvement techniques are merged in the formation of the new element in a complementary manner. The proposed quadratic finite element passes the patch tests, does not show spurious mechanism, and does not produce shear locking phenomena even with distorted meshes. It is shown that the element produces reliable solutions through numerical tests for standard benchmark problems. It is also noted that the element is applicable to transient dynamic analysis of Mindlin plates.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.485-500
• /
• 2000

In this paper the adaptive crack propagation analysis based on the estimated local and global error in the element-free Galerkin (EFG) method is presented. It is possible to keep consistency and accuracy of analysis in each propagation step by adaptive analysis. The adaptivity analysis in crack propagation is achieved by adding and removing the node along the background integration cell that are refined or recovered as estimated error. These errors are obtained by calculating the difference between the values of the projected stresses and original EFG stresses. To evaluate the performance of proposed adaptive procedure, the convergence behavior is investigated lot several examples. The results of these examples show the efficiency of proposed scheme in crack propagation analysis.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.379-395
• /
• 2021

For the complicated structural details in ships and offshore structures, the traditional hotspot stress approaches are known to be sensitive to the element variables of element topologies, sizes, and integration schemes. This motivated to develop a new approach for predicting reasonable hotspot stresses, which is less sensitive to the element variables and easy to be implemented the real marine structures. The three-point bending tests were conducted for the longitudinal attachments with the round and rectangular weld toes. The tests were reproduced in the numerical simulations using the solid and shell element models, and the simulation technique was validated by comparing the experimental stresses with the simulated ones. This paper considered three hotspot stress approaches: the ESM method based on surface stress extrapolation, the Dong's method based on nodal forces along a weld toe, and the proposed method based on nodal forces perpendicular to an imaginary vertical plane at a weld toe. In order to study the element sensitivities of each method, 16 solid element models and 8 shell element models were generated under the bending and tension loads, respectively. The element sensitivity was analyzed in terms of Stress Concentration Factors (SCFs) in viewpoints of two statistical quantities of mean and bias with respect to the reference SCFs. The average SCFs predicted by the proposed method were remarkably in good agreement with the reference SCFs based on the experiments and the ship rules. Negligibly small Coefficients of Variation (CVs) of the SCFs, which is measure of statistical bias, were drawn by the proposed method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.05a
• /
• pp.258-261
• /
• 2003

In this study, finite element analysis for multi-stage deep drawing process of rectangular configuration with extreme aspect ratio is carried out especially for the blank design. The analysis of rectangular deep drawing process with extreme aspect ratio is likewise very difficult with respect to the design process parameters including the intermediate die profile. In order to solve the difficulties, numerical approach using finite element method is performed in the present analysis and design. A series of experiments for multi-stage rectangular deep drawing process are conducted and the deformed configuration is investigated by comparing with the results of the finite element analysis. Additionally, to minimize amount of removal material after trimming process, finite element simulation is applied for the blank modification. The analysis incorporates brick elements for a rigid-plastic finite element method with an explicit time integration scheme using LS-DYNA3D.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.4
• /
• pp.602-609
• /
• 2001

This paper is concerned with the orthotropic problem of diaphragm-type air springs which consist of rubber linings, nylon reinforced rubber composite and bead ring. The analysis is carried out with a finite element method developed to consider the orthotropic properties, geometric nonlinearity using four-node degenerated shell element with reduced integration. Physical stabilization scheme is used to control the zeroenergy mode of the element. The analysis includes an inflation analysis and a lateral analysis of an air spring for the deformed shape and the spring load with respect to the vertical and l ateral deflection. Numerical results demonstrate the variation of the outer diameter, the fold height, the vertical force and the lateral force with respect to the inflation pressure and the lateral deflection.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.274-279
• /
• 2000

This paper is concerned with the orthotropic problem of diaphragm-type air springs which consist of rubber linings, nylon reinforced rubber composite and bead ring. The analysis is carried out with a finite element method developed to consider the orthotropic properties, geometric nonlinearity using four-node degenerated shell element with reduced integration. Physical stabilization scheme is used to control the zero-energy modes of the element. An inflation analysis and a lateral deformation analysis of an air spring are carried out. Numerical analysis results demonstrate the variation of the outer diameter, the fold height, the vertical force and the lateral force.