• Korean Journal of Computational Design and Engineering
• /
• v.1 no.3
• /
• pp.224-233
• /
• 1996

An automatic mesh generation scheme with triangular finite elements on three-dimensional surfaces has been developed. The surface triangulation process is performed as follows. To begin, surfaces with key nodes are transformed to two-dimensional planes and the meshes with triangular elements are constructed in these planes. Finally, the constructed meshes are transformed back to the original 3D surfaces. For the mesh generation, an irregular mesh generation scheme is employed in which local mesh densities are assigned by the user along the boundaries of the analysis domain. For this purpose a looping algorithm combined with an advancing front technique using basic operators has been developed, in which the loops are recursively subdivided into subloops with the use of the best split lines and then the basic operators generate elements. Using the split lines, the original boundaries are split recursively until each loop contains a certain number of key nodes, and then using the basic operators such as type-1 and type-2, one or two triangular elements are generated at each operation. After the triangulation process has been completed for each meshing domain, the resulting meshes are finally improved by smoothing process. Sample meshes are presented to demonstrate the versatility of the algorithm.

• Geomechanics and Engineering
• /
• v.9 no.4
• /
• pp.513-530
• /
• 2015

Commonly, the base stability of sheeted excavation pits against seepage failure by heave is evaluated by using two-dimensional groundwater flow models and Terzaghi's failure criterion. The objective of the present study is to investigate the effect of three-dimensional groundwater flow on the heave for sheeted excavation pits with various dimensions. For this purpose, the steady-state groundwater flow analyses are performed by using the finite element program ABAQUS 6.12. It has been shown that, in homogeneous soils depending on the ratio of half of excavation width to embedment depth b/D, the ratio of safety factor obtained from 3D analyses to that obtained from 2D analyses $FS_{(3D)}/FS_{(2D)}$ can reach up to 1.56 and 1.34 for square and circular shaped excavations, respectively. As failure body, both an infinitesimal soil column adjacent to the wall (Baumgart & Davidenkoff's criterion) and a three-dimensional failure body with the width suggested by Terzaghi for two-dimensional cases are used. It has been shown that the ratio of $FS_{(Terzaghi)}/FS_{(Davidenkoff)}$ varies between 0.75 and 0.94 depending on the ratio of b/D. Additionally, the effects of model size, the shape of excavation pit and anisotropic permeability on the heave are studied. Finally, the problem is investigated for excavation pits in stratified soils, and important points are emphasized.

• Proceeding of EDISON Challenge
• /
• 2016.03a
• /
• pp.195-202
• /
• 2016

In this paper, crack detection method using eigen value analysis is presented. Three methods are used: theoretical analysis, finite element method with the cracked beam elements and finite element method with three dimensional continuum elements. Finite element formulation of the cracked beam element is introduced. Additional term about stress intensity factor based on fracture mechanics theory is added to flexibility matrix of original beam to model the crack. As using calculated stiffness matrix of cracked beam element and mass matrix, natural frequencies are calculated by eigen value analysis. In the case of using continuum elements, the natural frequencies could be calculated by using EDISON CASAD solver. Several cases of crack are simulated to obtain natural frequencies corresponding the crack. The surface of natural frequency is plotted as changing with crack location and depth. Inverse analysis method is used to find crack location and depth from the natural frequencies of experimental data, which are referred by another papers. Predicted results are similar with the true crack location and depth.

• Structural Engineering and Mechanics
• /
• v.35 no.2
• /
• pp.241-256
• /
• 2010

The interaction between steel tube and concrete core is the key issue for understanding the behavior of concrete-filled steel tube columns (CFTs). This study investigates the force transfer by natural bond or by mechanical shear connectors and the interaction between the steel tube and the concrete core under three types of loading. Two and three-dimensional nonlinear finite element models are developed to study the force transfer between steel tube and concrete core. The nonlinear finite element program ABAQUS is used. Material and geometric nonlinearities of concrete and steel are considered in the analysis. The damage plasticity model provided by ABAQUS is used to simulate the concrete material behavior. Comparisons between the finite element analyses and own experimental results are made to verify the finite element models. A good agreement is observed between the numerical and experimental results. Parametric studies using the numerical models are performed to investigate the effects of diameterto-thickness ratio, uniaxial compressive strength of concrete, length of shear connectors, and the tensile strength of shear connectors.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3159-3174
• /
• 1994

A simple octree encoding algorithm based on a tetrahedron root has been developed to be used for fully automatic generation of three dimensional finite element meshes. This algorithm starts octree decomposition from a tetrahedron root node instead of a hexahedron root node so that the terminal mode has the same topology as the final tetrahedral mesh. As a result, the terminal octant can be used as a tetrahedral finite element without transforming its topology. In this part(I) of the thesis, an efficient algorithm for the tetrahedron-based octree is proposed. For this development, the following problems have been solved, : (1) an efficient data structure for storing the octree and finite elements, (2) an encoding scheme of a tetrahedral octree, (3) a neighbor finding technique for the tetrahedron-based octree.

• Restorative Dentistry and Endodontics
• /
• v.18 no.2
• /
• pp.291-316
• /
• 1993

Restorative procedures can lead to tooth fracture due to the relatively small amount of the remaining tooth structure. It is essential to prevent fractures by having a 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, MO amalgam cavity were prepared on maxillary first premolar and filled with amalgam. Three dimensional, model with 1365 8-node brick elements was made by serial photographic method. In this model, isthmus was varied in width at 1/4, 1/3, 1/2 and 2/3 of intercuspal width and material properties were given for three element groups, i.e., enamel, dentin and amalgam. A load of 500 N was applied vertically on amalgam and enamel. In case of enamel loading, 2 model (with and without amalgam) was compared to consider the possibility of play at the interface between tooth material and amalgam. These models were analyzed with three dimensional finite element method. The results were as follows: 1. The stress was concentrated on the facio-pulpal line angle and distal marginal ridge of the cavity. 2. With the increase of the isthmus width, the stress spread around the facio-pulpal line angle and the area of stress concentration moved toward the proximal box. 3. In case of narrow isthmus width, the initiation point of crack would be in the area of isthmus corner of the cavity, and with the increase of the isthmus width, it would move toward the proximal box and at the same time the possibility of crack increase at the distal marginal ridge. 4. The direction of crack progressed outward and downward from the facio-pulpal line angle, and with the increase of the isthmus width, it approximated vertical direction. At the marginal ridge, it occurred in vertical direction. 5. It would be favorable to make the isthmus width narrower than a third of the intercuspal width, and to cover the cusp if isthmus width were wider than half of the intercuspal width. 6. It is necessary to apply the possibility of play to the finite element analysis.

• Geotechnical Engineering
• /
• v.13 no.3
• /
• pp.33-52
• /
• 1997

A three dimensional nonlinear finite element analysis is used to study the influence of various design decisions for tieback walls. The numerical model simulates the soldier piles and the tendon bonded length of the anchors with beam elements, the unbonded tendon with a spring element, the wood lagging with the shell elements, and the soil with solid 3D nonlinear elements. The soil model used is a modified hyperbolic model with unloading hysteresis. The complete sequence of construction is simulated including the excavation, and the placement and stressing of the anchors. The numerical model is calibrated against a full scale instrumented tieback wall at the National Geotechnical Experimentation Site (NGES) on the Riverside Campus of Texas A&M University. Then a parametric study is conducted. The results give information on the influence of the following factors on the wall behavior : location of the first anchor, length of the tendon unbonded zone, magnitude of the anchor forces, embedment of the soldier piles, stiffness of the wood lagging, and of the piles. The implications in design are discussed.

• Restorative Dentistry and Endodontics
• /
• v.19 no.1
• /
• pp.231-254
• /
• 1994

Fracture of cusp, on posterior teeth, especially those carious or restored, is major cause of tooth loss. Inappropriate treatments, such as unnecessarily wide cavity preparations, increase the potential of further trauma and possible fracture of the remaining tooth structures. Fracture potential may be directly related to the stresses exerted upon the tooth during masticatory function. The purpose of this study is to evaluate the fracture resistance of tooth, restored with composite resin inlay. In this study, MOD inlay cavity prepared on maxillary first premolar and restored with composite resin inlay. Three dimensional finite element models with eight nodes isoparametric solid element, developed by serial grinding-photographing technique. These models have various occlusal isthmus and depth of cavity, 1/2, 1/3 and 1/4 of isthmus width and 0.7, 0.85 and 1.0 of depth of cavity. The magnitude of load was 474 N and 172 N as presented to maximal biting force and normal chewing force. These loads applied onto ridges of buccal and lingual cusp. These models analyzed with three dimensional finite element method. The results of this study were as follows : 1. There is no difference of displacement between width of occlusal isthmus and depth of cavity. 2. The stress concentrated at bucco-mesial comer, bucco-disal comer, pulpal line angle and the interface area between internal slopes of cusp and resin inlay. 3. The vector of stress direct to buccal and lingual side from center of cavity, to tooth surface going on to enamel. The magnitude of vector increase from occlusal surface to cervix. 4. The crack of tooth start interface area, between internal slop of buccal cusp and resin inlay. It progresses through buccopulpal line angle to cervix at buccomesial and buccodistal comer. 5. The influence with depth of cavity to fracture of tooth was more than width of isthmus. 6. It would be favorable to make the isthmus width narrower than a third of the intercuspal distance and depth of cavity is below 1 : 0.7.

• Journal of Biomedical Engineering Research
• /
• v.17 no.1
• /
• pp.61-70
• /
• 1996

In cementless total hip arthroplasty(THA), an initial stability of the femoral component is mandatory to achieve bony inyowth and secondary long term fixation. Primary stability of the femoral component can be obtained by minimizing the magnitude of relative micromotions at bone stem interface. An accurate evaluation of interf'ace micromotion and stress/strain fields in the bone-implant system may be relevant for better understanding of clinical situations and improving THA design. Recently finite element method(FEM) was introduced in'orthopaedic research field due to its unique capacity to evaluate stress in structure of complex shape, loading and material behavior. The authors developed the 3-dimensional finite element model of proximal femur with $Multilock^{TM}$ stem of 1179 blick elements to analyse the micromotions and mechanical behaviors at the bone-stem inteface in early post-operative period for the load simulating single leg stance. The results indicates that the values of relative motion for this well fit stem were $150{\mu}m$ in maximum $82{\mu}m$ in minimum and the largest relative motion was developed in medial region of Proximal femur and in anterior-posterior direction. The motion in the proximal bone was much greater than in the distal bone and the stress pattern showed high stress concentration on the cortex near the tip of the stem. These findings indicate that the loading on the hip joint in the early postoperative situation before achieving bony ingrowth could produce large micromotion of $150{\mu}m$ and clinicaly non-cemented THA patient should not be allowed weight bearing strictly early in the postoperative period.

• The Journal of Advanced Prosthodontics
• /
• v.13 no.4
• /
• pp.216-225
• /
• 2021

PURPOSE. The aim of this study is to evaluate the effects of canine guidance occlusion and group function occlusion on the degree of stress to the bone, implants, abutments, and crowns using finite element analysis (FEA). MATERIALS AND METHODS. This study included the implant-prosthesis system of a three-unit bridge made of monolithic zirconia and hybrid abutments. Three-dimensional (3D) models of a bone-level implant system and a titanium base abutment were created using the original implant components. Two titanium implants, measuring 4 × 11 mm each, were selected. The loads were applied in two oblique directions of 15° and 30° under two occlusal movement conditions. In the canine guidance condition, loads (100 N) were applied to the canine crown only. In the group function condition, loads were applied to all three teeth. In this loading, a force of 100 N was applied to the canine, and 200-N forces were applied to each premolar. The stress distribution among all the components of the implant-bridge system was assessed using ANSYS SpaceClaim 2020 R2 software and finite element analysis. RESULTS. Maximum stress was found in the group function occlusion. The maximum stress increased with an increase in the angle of occlusal force. CONCLUSION. The canine guidance occlusion with monolithic zirconia crown materials is promising for implant-supported prostheses in the canine and premolar areas.