• Transactions of the Korean Society of Automotive Engineers
• /
• v.1 no.1
• /
• pp.140-148
• /
• 1993

When a large automobile sheet metal part is formed in a draw die, the binder wrap is first calculated to predict the initial punch contact location for avoiding wrinkles and severe stretching of its thin blank sheet. Since the boundary of a pseudo blank in calculating a binder wrap by means of a geometrically nonlinear finite element method is unknown in advance, an iteration method is generally used. This paper presents an effective iteration method for correction of the pseudo blank in a binder wrap calculation. For the performance test, two examples are adopted. The calculated results for both examples show the good convergence which wasted solutions are obtained in the second iteration step.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1994.10a
• /
• pp.852-859
• /
• 1994

Optmal topology design is to search the optimal layout of the structure which can be used fot the shape of the conceptual design stage. Our objective is to maximize the stiffness of the structure under a material usage constraint. The density of each finite element is the design variable and its relationship with Young's modulus is expressed by quadratic form. The shape is represented by the entire density distribution, the structural analysis is performed by finite element method and the optimization is achieved by feasible direction method. Unlike optimality criteria method,feasible direction method can handle various problems simultaneously, that is, multi- objectives and multi-constraints. Total optimization time can be reduced by the approximation of the material property and fewer design variables than homogenization method. Topology optimization is applied to design the shape of ribs.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.703-715
• /
• 1996

In stochastic analysis, the randomness of the structural parameters is taken into consideration and the response variability is obtained in addition to the conventional (mean) response. In the present paper the structural response variability of plate structure is calculated using the weighted integral method and is compared with the results obtained by different methods. The stochastic field is assumed to be normally distributed and to have the homogeneity. The decomposition of strain-displacement matrix enabled us to extend the formulation to the stochastic analysis with the quadratic elements in the weighted integral method. A new auto-correlation function is derived considering the uncertainty of plate thickness. The results obtained in the numerical examples by two different methods, i.e., weighted integral method and Monte Carlo simulation, are in a close agreement. In the case of the variable plate thickness, the obtained results are in good agreement with those of Lawrence and Monte Carlo simulation.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.3 no.6
• /
• pp.96-111
• /
• 1995

A bellows is a component installed in the automobile exhaust system to reduce the impact from an engine. It's stiffness has a great influence on the natural frequency of the system. Therefore, it must be designed to keep the specified stiffness that requires in the system. This study present the finite element analysis of U-typed bellows using a curved conical frustum element and the shape optimal design with specified stiffness. The finite element analysis is verified by comparing with the experimental results. In the shape optimal design, the weight is considered as the cost function. The specified stiffness from the system design is transformed to equality constraints. The formulation has inequality constraints imposed on the fatigue limit, the natural frequencies, the buckling load and the manufacturing conditions. A procedure for shape optimization adopts a thickness, a corrugation radius, and a length of annular plate as optimal design variables. The external loading conditions include the axial and lateral loads with a boundary condition fixed at an end of the bellows. The recursive quadratic programming algorithm is selected to solve the problem. The result are compared with the existing bellows, and the characteristics of the bellows is investigated through the optimal design process. The optimized shape of the bellows are expected to give quite a good guideline to the practical design.

• Bulletin of the Korean Mathematical Society
• /
• v.48 no.5
• /
• pp.1079-1092
• /
• 2011

We consider the Newton's method for an direct solver of the optimal control problems of the Navier-Stokes equations. We show that the finite element solutions of the optimal control problem for Stoke equations may be chosen as the initial guess for the quadratic convergence of Newton's algorithm applied to the optimal control problem for the Navier-Stokes equations provided there are sufficiently small mesh size h and the moderate Reynold's number.

• Journal of Powder Materials
• /
• v.30 no.1
• /
• pp.29-34
• /
• 2023

A fixed-point iteration is proposed to integrate the stress and state variables in the incremental analysis of plastic deformation. The Conventional Newton-Raphson method requires a second-order derivative of the yield function to generate a complicated code, and the convergence cannot be guaranteed beforehand. The proposed fixed-point iteration does not require a second-order derivative of the yield function, and convergence is ensured for a given strain increment. The fixed-point iteration is easier to implement, and the computational time is shortened compared with the Newton-Raphson method. The plane-stress condition is considered for the biaxial loading conditions to confirm the convergence of the fixed-point iteration. 3-dimensional tensile specimen is considered to compare the computational times in the ABAQUS/explicit finite element analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2148-2158
• /
• 1996

Uncoupled, quasi-static and linear thermoviscoelastic problems are analyzed in time domain by the finite element approximation which is developed using the principle of virtual work and viscoelasticity matrices instead of shear and bulk relaxation functions as in usual formulations. The material is assumed to be isotropic, homegeneous and thermorheologically simple, which means that the temperature-time equivalence postulate is effective. The stress-strain laws are expressed by relaxation-type hereditary integrals. In spatial and time discritizations, isoparametric quadratic quadrilateral finite elements and linear time variations are adopted. For explicit derivations, the viscoelastic material is assumed to behave standard linear solid in shear and elastically in dilatation. Two-dimensional examples are solved under general temperature distributions T = T(x, t), and compared with other opproximate solutions to show the versatility of the presented analysis.

• Transactions of Materials Processing
• /
• v.7 no.6
• /
• pp.594-602
• /
• 1998

The design analysis of axisymmetric, multi-stage deep drawing dies was performed using the rigid-viscoplastic finite element formulation. In the formulation the axisymmetric CFS algorithm was employed. Hill's non-quadratic normal anisotropic yield criterion and isotropic hardening rule were considered. For trial initial displacements and tool contact points. the geometric force equilibrium method was adopted. In order to see the validity of the formulation, the multi-stage deep drawing processes of shell-cylinder front part of hydraulic booster were simulated. The simulation showed good agreements with measurments and PAM-STAMP results.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.855-859
• /
• 2004

We develop a new four-node flat shell element which is accurate, efficient, and suitable to be used on general purpose. The new element has a hybrid Trefftz element with drilling degrees of freedom as a membrane part. We define the two independent displacement field: the internal displacement field that satisfies governing equations in the domain a priori and the boundary displacement field that is usually used as a conventional finite element method. The hybrid Trefftz variational formulation connects these two displacement fields on the boundary of the domain. To add drilling degrees of freedom, we introduce the Allman's quadratic displacement field to the boundary displacement field. As a result, our flat shell element has 6 degrees of freedom per a node. We also use the well-known DKMQ plate bending element for the plate part of the proposed element. The DKMQ element satisfies Mindlin-Reissner‘s plate theory along the edge of the element and gives proper behavior regardless of the thickness. A series of numerical experiments shows that the performance of the new element such as accuracy, rate of convergence, robustness to mesh quality, and so on.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.