• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.345-352
• /
• 2014

In this paper, the CAD data for the optimal shape design obtained by isogeometric shape optimization is directly used to fabricate the specimen by using 3D printer for the experimental validation. In a conventional finite element method, the geometric approximation inherent in the mesh leads to the accuracy issue in response analysis and design sensitivity analysis. Furthermore, in the finite element based shape optimization, subsequent communication with CAD description is required in the design optimization process, which results in the loss of optimal design information during the communication. Isogeometric analysis method employs the same NURBS basis functions and control points used in CAD systems, which enables to use exact geometrical properties like normal vector and curvature information in the response analysis and design sensitivity analysis procedure. Also, it vastly simplify the design modification of complex geometries without communicating with the CAD description of geometry during design optimization process. Therefore, the information of optimal design and material volume is exactly reflected to fabricate the specimen for experimental validation. Through the design optimization examples of elasticity problem, it is experimentally shown that the optimal design has higher stiffness than the initial design. Also, the experimental results match very well with the numerical results. Using a non-contact optical 3D deformation measuring system for strain distribution, it is shown that the stress concentration is significantly alleviated in the optimal design compared with the initial design.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.10
• /
• pp.1139-1146
• /
• 2012

An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.12
• /
• pp.1401-1409
• /
• 2009

All the forces in the real world act dynamically on structures. Design and analysis should be performed based on the dynamic loads for the safety of structures. Dynamic (transient or vibrational) responses have many peaks in the time domain. Topology optimization, which gives an excellent conceptual design, mainly has been performed with static loads. In topology optimization, the number of design variables is quite large and considering the peaks is fairly costly. Topology optimization in the frequency domain has been performed to consider the dynamic effects; however, it is not sufficient to fully include the dynamic characteristics. In this research, linear dynamic response topology optimization is performed in the time domain. First, the necessity of topology optimization to directly consider the dynamic loads is verified by identifying the relationship between the natural frequency of a structure and the excitation frequency. When the natural frequency of a structure is low, the dynamic characteristics (inertia effect) should be considered. The equivalent static loads (ESLs) method is proposed for linear dynamic response topology optimization. ESLs are made to generate the same response field as that from dynamic loads at each time step of dynamic response analysis. The method was originally developed for size and shape optimizations. The original method is expanded to topology optimization under dynamic loads. At each time step of dynamic analysis, ESLs are calculated and ESLs are used as the external loads in static response topology optimization. The results of topology optimization are used to update the design variables (density of finite elements) and the updated design variables are used in dynamic analysis in a cyclic manner until the convergence criteria are satisfied. The updating rules and convergence criteria in the ESLs method are newly proposed for linear dynamic response topology optimization. The proposed updating rules are the artificial material method and the element elimination method. The artificial material method updates the material property for dynamic analysis at the next cycle using the results of topology optimization. The element elimination method is proposed to remove the element which has low density when static topology optimization is finished. These proposed methods are applied to some examples. The results are discussed in comparison with conventional linear static response topology optimization.

• Journal of Powder Materials
• /
• v.11 no.5
• /
• pp.399-411
• /
• 2004

Crack generation during die compaction and distortion during sintering have been critical problems for the conventional pressing and sintering process. Until now, trial and error approach with engineers' industrial experiences has been only solution to protect the crack generation and distortion. However, with complexity in shape and process it is very difficult to design process conditions without CAE analysis. We developed the exclusive CAE software (PMsolver/Compaction) for die compaction process. The accuracy of PMsolver is verified by comparing the finite element simulation results with experimental results. The simplified procedures to find material properties are proposed and verified with iron based powder and tungsten carbide powder. Based on the accurate simulation result by PMsolver, the optimal process conditions are designed to get uniform density distribution in a powder compact after die compaction process by using a derivative based optimization scheme. In addition, the effect of non-uniform density distribution in a powder compact on distortion during sintering is shown in case of the fabrication of tungsten carbide insert.

• Journal of Ocean Engineering and Technology
• /
• v.20 no.6 s.73
• /
• pp.123-129
• /
• 2006

The mechanical properties of bellows, such as the extensibility and the strength can be changed depending on the shape. For the shipbuilding material, it is desirable that the fatigue life is long due to the elastic property and the reduction of thermal stress in piping system. Nowadays, the domestic production and design of bellows are based on the E.J.M.A. Code. Therefore, the design standard is in need because of much errors and lack of detailed analysis. In this study, it is attempted to find out the optimal shape of U-type bellows using the finite element analysis. The design factors, mountain height, length, thickness, and the number of convolutions are considered and the proper values are chosen for the simulation. The results shaw that as the number of convolutions reduces, the volume decreases while the stress increases. However, as the number of convolutions increases, the volume increases above the standard volume and the stress obviously increases. In addition, the effect of the thickness of bellows on the stress is very large. Both of the mass and stress are decreasing at a certain lower value region. Also, we investigated shape optimization with considering maximum stress distribution tendency.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.3
• /
• pp.173-184
• /
• 2013

A multi-level design optimization framework for aerodynamic design of rotary wing such as propeller and helicopter rotor blades is presented in this study. Strategy of the proposed framework is to enhance aerodynamic performance by sequentially applying the planform and sectional design optimization. In the first level of a planform design, we used a genetic algorithm and blade element momentum theory (BEMT) based on two-dimensional aerodynamic database to find optimal planform variables. After an initial planform design, local flow conditions of blade sections are analyzed using high-fidelity CFD methods. During the next level, a sectional design optimization is conducted using two dimensional Navier-Stokes analysis and a gradient based optimization algorithm. When optimal airfoil shape is determined at the several spanwise locations, a planform design is performed again. Through this iterative design process, not only an optimal flow condition but also an optimal shape of an EAV propeller blade is obtained. To validate the optimized propeller-blade design, it is tested in wind-tunnel facility with different flow conditions. An efficiency, which is slightly less than the expected improvement of 7% predicted by our proposed design framework but is still satisfactory to enhance the aerodynamic performance of EAV system.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.4
• /
• pp.372-378
• /
• 2013

An acoustical shape optimization problem is formulated for optimal design of a perforated reactive muffler with offset inlet/outlet. The mean transmission loss value in a target frequency range is maximized for an allowed pressure drop value between an inlet and an outlet. Partitions in the chamber are divided into several sub-partitions, whose lengths are selected as design variables. Each sub-partition has the same number of holes, whose sizes are equal. A finite element model is employed for acoustical and flow analyses. A gradient-based optimization algorithm is used to obtain an optimal muffler. The acoustical and fluidic characteristics of the optimal muffler are compared with those of a reference muffler. Validation experiment is carried out to support the effectiveness of our suggested method.

• Journal of the Korean Society for Precision Engineering
• /
• v.31 no.5
• /
• pp.411-416
• /
• 2014

In this research we present design optimization methods for a vehicle-mounted satellite antenna positioner. Our initial antenna positioner was conservatively designed to satisfy a worst case scenario where wind blew across the positioner at the speed of 120 km/h. Investigating stresses and safety based on Finite Element Methods (FEM), we find reflector support frames can be optimized to significantly reduce the weight of the positioner system. Thus, we optimize the reflector support frame from the given initial design while considering weight, maximum stress, maximum allowable deflection, cross section, and thickness. As a result, Shape C and the thickness of 2 mm are determined for the cross section of the reflector support frame. Applying this result, the weight of the new antenna positioner is 57.343 kg, which is decreased by 10.74% compared to the initial conservative design.

• Transactions of Materials Processing
• /
• v.31 no.5
• /
• pp.302-308
• /
• 2022

In this study, the weight of the platen was reduced using the structural strength analysis and topology optimization design of the extruder by finite element analysis. The main components of the extruder such as the stem and billet, were modeled, and the maximum stress and safety factor were verified through structural strength analysis. Based on the results of the structural strength analysis, the optimal phase that satisfies the limitation given to the design area of the structure and maximizes or minimizes the objective function was obtained through a numerical method. The platen was redesigned with a phase-optimal shape, the weight was reduced by 40% (from the initial weight of 11.1 tons to 6.6 tons), and the maximum stress was 147.49 MPa safety factor of 1.86.

• Architectural research
• /
• v.15 no.4
• /
• pp.207-214
• /
• 2013

A form-finding technique is proposed for three-dimensional spatial structures. Two-step discrete finite element (FE) mesh generator based on computer aided geometric design (CAGD) is introduced and used to control the shape of three-dimensional spatial structures. Mathematical programming technique is adopted to search new forms (or shapes) of spatial structures. For this purpose, the strain energy is introduced as the objective function to be minimized and the initial volume (or the initial weight) is considered as constraint function. Numerical examples are carried out to test the capability of the proposed form-finding techniques and provided as benchmark tests.