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

In this study, we developed a methodology to determine the reference parameter for an aircraft impact induced risk assessment of nuclear power plant (NPP) using finite element impact analysis of containment building. The target structure used to develop the method of reference parameter selection is one of the typical Korean PWR type containment buildings. We composed a three-dimensional finite element model of the containment building. The concrete damaged plasticity model was used for the concrete material model. The steels in the tendon, rebar, and liner were modeled using the piecewise-linear stress-strain curves. To evaluate the correlations between structural response and each candidate parameter, we developed Riera's aircraft impact force-time history function with respect to the variation of the loading parameters, i.e., impact velocity and mass of the remaining fuel. For each force-time history, the type of aircraft is assumed to be a Boeing 767 model. The variation ranges of the impact velocity and remaining fuel percentage are 50 to 200m/s, and 30 to 90%, respectively. Four parameters, i.e., kinetic energy, total impulse, maximum impulse, and maximum force are proposed for candidates of the reference parameter. The wellness of the correlation between the reference parameter and structural responses was formulated using the coefficient of determination ($R^2$). From the results, we found that the maximum force showed the highest $R^2$ value in most responses in the materials. The simplicity and intuitiveness of the maximum force parameter are also remarkable compared to the other candidate parameters. Therefore, it can be concluded that the maximum force is the most proper candidate for the reference parameter to assess the aircraft impact induced risk of NPPs.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.19 no.4
• /
• pp.257-268
• /
• 2003

The purpose of this study was to assess the loading distributing characteristics of implant prosthesis according to position and direction of load, under vertical and inclined loading using FEA analysis. The finite element model was designed according to standard fixture (4.1mm restorative component x 11.5mm length). The crown for mandibular first molar was made using UCLA abutment. Each three-dimensional finite element model was created with the physical properties of the implant and surrounding bone. This study simulated loads of 200N at the central fossa in a vertical direction (loading condition A), 200N at the outside point of the central fossa with resin filling into screw hole in a vertical direction (loading condition B), 200N at the centric usp in a $15^{\circ}$ inward oblique direction (loading condition C), 200N at the in a $30^{\circ}$ inward oblique direction (loading condition D) or 200N at the centric cusp in a $30^{\circ}$ outward oblique direction (loading condition E) individually. Von Mises stresses were recorded and compared in the supporting bone, fixture, and abutment screw. The following results have been made based on this study: 1. Stresses were concentrated mainly at the ridge crest around implant in both vertical and oblique loading but stresses in the cancellous bone were low in both vertical and oblique loading. 2. Bending moments resulting from non-axial loading of dental implants caused stress concentrations on cortical bone. The magnitude of the stress was greater with the oblique loading than with the vertical loading. 3. An offset of the vertical occlusal force in the buccolingual direction relative to the implant axis gave rise to increased bending of the implant. 4. The relative positions of the resultant line of force from occlusal contact and the center of rotation seems to be more important. 5. The magnitude of the stress in the supporting bone, fixture and abutment screw was greater with the outward oblique loading than with the inward oblique loading and was the greatest under loading at the centric cusp in a $30^{\circ}$ outward oblique direction. Conclusively, this study provides evidence that bending moments resulting from non-axial loading of dental implants caused stress concentrations on cortical bone. But it seems to be more important that how long is the distance from center of rotation of the implant itself to the resultant line of force from occlusal contact(leverage). The goal of improving implants should be to avoid bending of the implant.

• Composites Research
• /
• v.25 no.4
• /
• pp.98-104
• /
• 2012

In recent years, degenerative spinal instability has been effectively treated with a cage. However, little attention is focused on the stiffness of the cage. Recent advances in the medical implant industry have resulted in the use of medical carbon fiber reinforced polymer (CFRP) cages. The biomechanical advantages of using different cage material in terms of stability and stresses in bone graft are not fully understood. A previously validated three-dimensional, nonlinear finite element model of an intact L2-L5 segment was modified to simulate posterior interbody fusion cages made of CFRP and titanium at the L4-L5 disc with pedicle screw, to investigate the effect of cage stiffness on the biomechanics of the fused segment in the lumbar region. From the results, it could be found that the use of a CFRP cage would not only reduce stress shielding, but it might also have led to increased bony fusion.

• Maxillofacial Plastic and Reconstructive Surgery
• /
• v.30 no.4
• /
• pp.345-352
• /
• 2008

Purpose: This study was performed to evaluate the stress distribution in the bone and the displacement distribution of the miniscrew under orthopedic force with two different types of miniplate design as skeletal anchorage for orthopedic treatment. Materials and methods: Finite element models were made for 6-hole miniplate (0.8mm in thickness), which were designed in two different shapes-one is curvilinear shaped (C plate, Jeil Medical Co., Korea) and another, Y shaped (Y plate), fixed with 3 pieces of miniscrew 2mm-diameter and 6mm-long respectively. A traction force of 4 N was applied in $0^{\circ}$, $30^{\circ}$ and $60^{\circ}$ to imaginary axis connecting two unfixed distalmost holes of the miniplate. Results: The maximum von Mises stress in the bone was much greater in the cortical portion rather than in the cancellous portion. C plate showed greater maximum von Mises stress in the cortical bone than Y plate. The maximum displacement of the miniscrew was greater in C plate than Y plate. The more increased the angle of the applied orthopedic force, the greater maximum von Mises stress in the bone and maximum displacement of the miniscrew. It was observed that in C plate, the von Mises stress in the bone and displacement of the miniscrew were distributed around the distalmost screw-fixed area. Conclusions: The results suggest that Y plate should have the advantage over C plate and in the placement of the miniplate, its imaginary axis should be placed as parallel as possible to the direction of orthopedic force to obtain its primary stability.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.4
• /
• pp.5-15
• /
• 2017

Three-hinge failure occurs in a jointed rock slope with a joint set parallel with slope and a conjugate joint set. Limit Equilibrium Method (LEM) and Finite Element Method (FEM) which are commonly used for slope design, are not suitable for evaluating stability against three-hinge failure, and this study performed parametric study to analyze the failure mechanism and to find influence factors causing three-hinge failure using UDEC which is a commercial two-dimensional DEM based numerical program. Numerical analyses were performed for various joint structural conditions and joint properties as well as ground water conditions. It was found that pore water pressure is the main factor triggering the three-hinge failure and the mode of failure depends on friction angle of basal joint and bedding joint set. The results obtained from this study can be used for adequate and economic footwall slope reinforcement design and construction.

• Restorative Dentistry and Endodontics
• /
• v.18 no.2
• /
• pp.395-411
• /
• 1993

The purpose of this study was to examine the clear concept of the designs for cavity preparations. Among the several parameters in cavity designs, profound understanding of isthmus width factor would facilitate selection of the appropriate cavity preparation for a specific clinical situation. In this study, the cavities were prepared on maxillary first premolar and filled with gold inaly. A two - dimensional model was composed of 1037 - node triangle elements. In this model, isthmus was varied in width at 1/4, 1/3 and 1/2 of intercuspal width and material properties were given for four element groups, i.e., enamel, dentin, pulp and gold. The 500N occlusal load varied in direction and it was examined using three types of load : concentrated load, divided load and distributed load. The models were also examined with empty cavities using the devided load and distributed load. These models were analyzed the displacement and strees distribution by the two - dimensional Finite Element Method. The results were as follows : 1. All experimental models which filled with gold inlay after cavity preparation were similar direction of displacement with control model under same load type. But in the models with empty cavities, as isthmus width was wider, the degree of displacement was increased at same load type. 2. Among the experimental models which were filled with gold inaly after cavity preparation, the model II showed the least stress concentration under concentrated load and divided load. But in the models with empty cavities, the model III showed the largest stress concentration and tooth fracture is expected regardless isthmus width. 3. All experimental models showed similar displacement pattern beneath restorative material under a concentrated load. In the models with empty cavities, a divided load resulted in a lingual displacement of the lingual cusp, but a distributed load resulted in a buccal displacement of the lingual cusp. In regard to the above results, the restored models were stronger than empty models in respect to the bending moment and tensile stress. The empty models are expected to fracture regardless isthmus width. The safest isthmus width was 1/3 of intercuspal distance, which showed the least stress concentration in respect to the effect of stress distribution.

• Journal of the Korean GEO-environmental Society
• /
• v.22 no.4
• /
• pp.15-23
• /
• 2021

Various micropiles have been developed through research related to micropiles, which have been carried out with the increased use of micropiles. Among the micropile construction methods being developed, the triaxial micropile (tmp), which is recently developed for the purpose of increasing the horizontal bearing capacity (seismic resistance), is representative. The three-axis micropile has the advantage of a method that can resist horizontal load more effectively because three micropiles installed inclined on each axis resist horizontal load. However, there is a problem in effectively using this pile method due to insufficient research on the support characteristics of the triaxial group micropile. In order to effectively utilize the triaxial group micropile (tmp), it is required to evaluate the bearing capacity for the factors that affect the horizontal bearing capacity of the pile. Therefore, in this study, field horizontal loading Tests were performed for each load direction, field loading Tests were verified through three-dimensional finite element analysis, behavioral characteristics of triaxial micropiles were evaluated, and appropriate horizontal bearing capacity was analyzed in consideration of horizontal load directions.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.411-428
• /
• 2022

This paper presents the Campbell diagram analysis of the rotordynamic system using the full order model (FOM) and the reduced order model (ROM) techniques to determine the critical speeds, identify the stability and reduce the computational time. Due to the spin-speed-dependent matrices (e.g., centrifugal stiffening matrix), several model order reduction (MOR) techniques may be considered, such as the modal superposition (MS) method and the Krylov subspace-based MOR techniques (e.g., Ritz vector (RV), quasi-static Ritz vector (QSRV), multifrequency quasi-static Ritz vector (MQSRV), multifrequency/ multi-spin-speed quasi-static Ritz vector (MMQSRV) and the combined Ritz vector & modal superposition (RV+MS) methods). The proposed MMQSRV method in this study is extended from the MQSRV method by incorporating the rotational-speed-dependent stiffness matrices into the Krylov subspace during the MOR process. Thus, the objective of this note is to respond to the question of whether to use the MS method or the Krylov subspace-based MOR technique in establishing the Campbell diagram of the shaft-disc-blade assembly systems in three-dimensional (3D) finite element analysis (FEA). The Campbell diagrams produced by the FOM and various MOR methods are presented and discussed thoroughly by computing the norm of relative errors (ER). It is found that the RV and the MS methods are dominant at low and high rotating speeds, respectively. More precisely, as the spinning velocity becomes large, the calculated ER produced by the RV method is significantly increased; in contrast, the ER produced by the MS method is smaller and more consistent. From a computational point of view, the MORs have substantially reduced the time computing considerably compared to the FOM. Additionally, the verification of the 3D FE rotordynamic model is also provided and found to be in close agreement with the existing solutions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.39-46
• /
• 1993

The need to develop method which can improve the shape design efficiency using high quality approximation is being brought up. In this study, to perform shape optimal design of three-dimensional continuum structures an efficient approximation method for stress constraints is proposed, based on expanding the nodal forces in Taylor series with respect to shape variables. Numerical examples are performed using the 3-D cantilever beam and fixed-fixed beam and compared with other method to demonstrate the efficiency and convergence rate of the Force Approximation method. It is shown that by taking advantage of this high quality approximation, the total number of finite element analysis required for shape optimization of 3-D continuum structures can be reduced significantly, resulting to the same level of efficiency achieved previously in sizing optimization problems. Also, shape representation by super curve technique applied to obtain optimal shape finds useful method.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.