• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.59-68
• /
• 2003

An optimization Program is developed to produce new topologies of plane structures under multiload case. A four-node finite element is used in the response analysis to reduce the computation time and to ultimately achieve practical topology optimization. The bilinear finite element is prone to produce chequer-boarding phenomenon and a simple filtering process is therefore adopted. An artificial material model is employed to represent the structural material and the resizing algorithm based on the optimality criteria is adopted to update the material density parameter during optimization process. With newly developed optimization program, the comparison study has been made between single and multiload cases and its results are described in this paper. From numerical results, it appears that multiload case should be considered to achieve the practical topology optimization.

• Archives of Plastic Surgery
• /
• v.46 no.4
• /
• pp.386-389
• /
• 2019

The development of breast implant technology continues to evolve over time, but changes in breast shape after implantation have not been fully elucidated. Thus, we performed computerized finite element analysis in order to better understand the trajectory of changes and stress variation after breast implantation. The finite element analysis of changes in breast shape involved two components: a static analysis of the position where the implant is inserted, and a dynamic analysis of the downward pressure applied in the direction of gravity during physical activity. Through this finite element analysis, in terms of extrinsic changes, it was found that the dimensions of the breast implant and the position of the top-point did not directly correspond to the trajectory of changes in the breast after implantation. In addition, in terms of internal changes, static and dynamic analysis showed that implants with a lower top-point led to an increased amount of stress applied to the lower thorax. The maximum stress values were 1.6 to 2 times larger in the dynamic analysis than in the static analysis. This finding has important implications for plastic surgeons who are concerned with long-term changes or side effects, such as bottoming-out, after anatomic implant placement.

• Transactions of Materials Processing
• /
• v.11 no.4
• /
• pp.302-311
• /
• 2002

A generalized numerical approach based on the finite element method to analysis and design of hydroforming process is proposed in this paper. The special attention is focused on comparison of an implicit and an explicit finite element method widely used for hydroforming simulation. Furthermore, in order to meet the increasing real needs for prediction of forming limit, a ductile fracture criterion combined with finite element method is introduced and then applied to hydroforming process of an automobile lower m Consequently, the numerical analysis and design for hydroforming process presented here will facilitate the development and application of the tube hydrofoniung process to a new level.

• Structural Engineering and Mechanics
• /
• v.29 no.5
• /
• pp.469-488
• /
• 2008

Finite element methods have often been used for structural analyses of various mechanical problems. When finite element analyses are utilized to resolve mechanical systems, numerical uncertainties in the initial data such as structural parameters and loading conditions may result in uncertainties in the structural responses. Therefore the initial data have to be as accurate as possible in order to obtain reliable structural analysis results. The typical finite element method may not properly represent discrete systems when using uncertain data, since all input data of material properties and applied loads are defined by nominal values. An interval finite element analysis, which uses the interval arithmetic as introduced by Moore (1966) is proposed as a non-stochastic method in this study and serves a new numerical tool for evaluating the uncertainties of the initial data in structural analyses. According to this method, the element stiffness matrix includes interval terms of the lower and upper bounds of the structural parameters, and interval change functions are devised. Numerical uncertainties in the initial data are described as a tolerance error and tree graphs of uncertain data are constructed by numerical uncertainty combinations of each parameter. The structural responses calculated by all uncertainty cases can be easily estimated so that structural safety can be included in the design. Numerical applications of truss and frame structures demonstrate the efficiency of the present method with respect to numerical analyses of structural uncertainties.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.2
• /
• pp.423-440
• /
• 2016

The bubble stabilization technique of Chebyshev-Legendre high-order element methods for one dimensional advection-diffusion equation is analyzed for the proposed scheme by Canuto and Puppo in [8]. We also analyze the finite element lower-order preconditioner for the proposed stabilized linear system. Further, the numerical results are provided to support the developed theories for the convergence and preconditioning.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.301-316
• /
• 2002

An important goal of computational wind engineering is to impact the design process with simulations of flow around buildings and bridges. One challenging aspect of this goal is to solve the Navier-Stokes (NS) equations accurately. For the unsteady computations, an adaptive finite element technique may reduce the computer time and storage. The preliminary application of a p-version as well as an h-version adaptive technique to computational wind engineering has been reported in previous paper. The details on the implementation of p-adaptive technique will be discussed in this paper. In this technique, two posteriori error estimations, which are based on the velocity and vorticity, are first presented. Then, the polynomial order of the interpolation function is increased continuously element by element until the estimated error is less than the accepted. The second through sixth orders of hierarchical functions are used as the interpolation polynomials. Unequal order interpolations are used for velocity and pressure. Using the flow around a circular cylinder with Reynolds number of 1000 the two error estimators are compared. The result show that the estimated error based on the velocity is lower than that based on the vorticity.

• Journal of the Korean Society for Precision Engineering
• /
• v.24 no.11
• /
• pp.116-125
• /
• 2007

Human muscle skeletal model was developed for biomechanical study. The human model was consists with 19 bone-skeleton and 122 muscles. Muscle number of upper limb, trunk and lower limb part are 28, 60, 34 respectively. Bone was modeled with 3D beam element and muscle was modeled with spar element. For upper limb muscle modelling, rectus abdominis, trapezius, deltoideus, biceps brachii, triceps brachii muscle and other main muscles were considered. Lower limb muscle was modeled with gastrocenemius, gluteus maximus, gluteus medius and related muscles. The biomechanical stress and strain analysis of human was conducted by proposed finite element analysis model under Kumdo head hitting motion. In this study structural analysis has been performed in order to investigate the human body impact by Kumdo head hitting motion. As the results, the analytical displacement, stress and strain of human body are presented.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.8
• /
• pp.1577-1584
• /
• 2002

This paper develops a finite element model for studying the occupant behavior and injury cofficients of a large-sized cab-over type truck. Since it does not have a room to absorb collision energy and deformation in front of the passenger compartment the deformation is directly transmitted to the passenger compartment. Moreover, since its steering column is attached on the frame, severe deformation of the frame directly affects on the steering wheel's movement. Therefore, if the occupant behavior and injury coefficients analysis is performed using a finite element model developed based on a sled test, it is very difficult to expect acquiring satisfactory results. Thus, the finite element model developing in this paper is based on the frontal crash test in order to overcome the inherent problems of the sled test based model commonly used in the passenger car. The occupant behavior and injury coefficients analysis is performed using PAM-CRASH installed in super-computer SP2. In order to validate the reliability of the developed finite element model, a frontal crash test is carried out according to a test method used fur developing truck occupant's secondary safety system in european community and japan. That is, test vehicle's collision direction is vertical to the rigid barrier and collision velocity is 45kph. Thus, measured vehicle pulses at the lower parts of the left and right B-pilla., dummy chest and head deceleration profiles, HIC(head injury criterial) and CA(chest acceleration) values, and dummy behavior from the frontal crash test are compared to the analysis results to validate reliability of the developed model.

• The Journal of Korean Academy of Prosthodontics
• /
• v.43 no.6
• /
• pp.727-735
• /
• 2005

Statement of problem. Platform switching in implant prosthesis has been used for esthetic and biological purpose. But there are few reports for this concept. Purpose. The purpose of this study is evaluation of platform switching in wide implant by three dimensional finite element analysis. Materials and Methods. The single implant and prosthesis was modeled in accordance with the geometric designs for Osstem implant system. Three-dimensional finite element models were developed for (1) a wide diameter 3i type titanium implant 5 mm in diameter, 13 mm in length with wide cemented abutment, titanium alloy abutment screw, and prosthesis (2) a wide diameter 3i type titanium implant 5 mm in diameter, 13 mm in length with regular cemented abutment, titanium alloy abutment screw and prosthesis(platform switching) was made for finite element analysis. The abutment screws were subjected to a tightening torque of 30 Ncm. The amount of preload was hypothesized to 650N, and round and flat type prostheses were loaded to 200 N. Four loading offset point (0, 2, 4, 6 mm from the center of the implants) were evaluated. Models were processed by the software programs HyperMesh and ANSA. The PAM-CRASH 2G simulation software was used for analysis of stress. The PAM-VIEW and HyperView were used for post processing. Results. The results from experiment were as follows; 1. von Mises stress value is increased in order of bone, abutment, implant and abutment screw. 2. von Mises stress of abutment screw is lower when platform switching. 3. von Mises stress of implant is lower when platform switching until loading offset 4 mm. 4. von Mises stress of abutment is similar between each other. 5. von Mises stress of bone is slightly higher when platform switching. Conclusion. The von Mises stress pattern of implant components is favor when platform switch ing but slightly higher in bone stress distribution than use of wide abutment. The research about stress distribution is essential for investigation of the cortical bone loss.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.406-411
• /
• 2001

The expandable implant employs an inner expansion screw in order to expand several legs of implant. Compressive stresses are produced at the bone tissue surrounding the implant, and the contact area between the implant and the bone tissue is increased, which result in increased resistance to horizontal and vertical pressure loads. The stress distribution in implant is also an important factor. Three types of implant models including an existing one have been investigated by using the Finite Element Method, and an improved design model has been suggested.