• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.1487-1492
• /
• 2007

In this paper, we apply a forging simulator to automatic simulation of a connecting-rod forging process and compare its results with others found in the literature. The process information comes from the distributed examples of DEFORM3D. The process is fully automatically simulated using the tetrahedral element capability of AFDEX3D, developed by the authors. Our results are compared with the results found from the related literature, already simulated using DEFORM3D by other researchers. The comparison shows that our results are relatively excellent especially in terms of mesh quality on which the solution accuracy depends mainly.

• Journal of Power System Engineering
• /
• v.14 no.3
• /
• pp.71-76
• /
• 2010

The safety block which prevents drop of laborers at high altitude was analyzed by finite element method. Elastic analysis was done by Ansys ver. 11.0. and tetrahedral meshing was used. As load applied more vertically at the fixed face of saw tooth, the stress concentration became smaller and the load distributed broader. When load worked at saw tooth and the shape was more straight to the direction of load, most stresses except principal stress became smaller. When the area of the load increased, principal stress and equivalent stress could be decreased simultaneously. A principal stress and other various stresses occurred in 3D shape, therefore revised model which has smaller equivalent stress than other models shows excellence on the stability and the credibility.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.7
• /
• pp.74-79
• /
• 2000

Integrating a 3D solid modeler with a general purpose FEM code, an automatic nonlinear analysis system of the 3D crack problems has been developed. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy knowledge processing. Nodes are generated by the bucketing method, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. The complete finite element(FE) model generated, and a stress analysis is performed. In this system, burden to analysts fur introducing 3D cracks to the FE model as well as fur estimating their fracture mechanics parameters can be dramatically reduced. This paper describes the methodologies to realize such functions, and demonstrates the validity of the present system.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• Proceedings of the KOSOMBE Conference
• /
• v.1995 no.11
• /
• pp.34-38
• /
• 1995

Using a 3-D finite element method (FEM), the biomechanical characteristics of a threaded truncated acetabular component and a porous coated hemispherical acetabular component were studied. The Von-Mises stress/strain patterns in the acetabulum reconstructed with these two different types of cementless acetabular cups were investigated. The geometry and dimensions of human hemi-pelvis used in the present shape modeling for finite element analysis were scanned with a 3-D laser scanner(TDS-9000, Cyberware, USA). The scanned data was numerically handled with a shape modelling software 'Pro-Engineer'. Using 19836, 16853 tetrahedral elements, respectively, the stress and displacement field of the acetabulum reconstructed with the two different types of the acetabular components were computed. While the hemi-sphere component was found to show a relatively similar stress/strain patterns to those in the normal hip, the results with the threaded cup showed a considerably different patterns from those in the normal condition. Several regions in cancellous bone near the threads and the edge of the truncated cup was found to be overstressed, especially in the superior-lateral part of the acetabulum. It was postulated that the excessive reaming-out of subchondral bone layer when the truncated cup was used can cause the presence of these overstressed regions of cancellous bone. This theoretical prediction for the implanted acetabulum appeared to consistent with the pathological observation of proximal/medial migration of the threaded truncated acetabular prostheses in the previous publications.

• The Journal of Advanced Prosthodontics
• /
• v.12 no.2
• /
• pp.67-74
• /
• 2020

PURPOSE. This study evaluated the influence of prosthesis weight and number of implants on the bone tissue microstrain. MATERIALS AND METHODS. Fifteen (15) fixed full-arch implant-supported prosthesis designs were created using a modeling software with different numbers of implants (4, 6, or 8) and prosthesis weights (10, 15, 20, 40, or 60 g). Each solid was imported to the computer aided engineering software and tetrahedral elements formed the mesh. The material properties were assigned to each solid with isotropic and homogeneous behavior. The friction coefficient was set as 0.3 between all the metallic interfaces, 0.65 for the cortical bone-implant interface, and 0.77 for the cancellous bone-implant interface. The standard earth gravity was defined along the Z-axis and the bone was fixed. The resulting equivalent strain was assumed as failure criteria. RESULTS. The prosthesis weight was related to the bone strain. The more implants installed, the less the amount of strain generated in the bone. The most critical situation was the use of a 60 g prosthesis supported by 4 implants with the largest calculated magnitude of 39.9 mm/mm, thereby suggesting that there was no group able to induce bone remodeling simply due to the prosthesis weight. CONCLUSION. Heavier prostheses under the effect of gravity force are related to more strain being generated around the implants. Installing more implants to support the prosthesis enables attenuating the effects observed in the bone. The simulated prostheses were not able to generate harmful values of peri-implant bone strain.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.18 no.1
• /
• pp.122-126
• /
• 2008

Integrating a 3D solid modeler with a general purpose FEM code, an automatic stress intensity factor analysis system of the 3D crack problems has been developed. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy knowledge processing. Nodes are generated and quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. Finally, the complete finite element(FE) model generated, and a stress analysis is performed. This paper describes the methodologies to realize such functions, and demonstrates the validity of the present system.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.1
• /
• pp.1-6
• /
• 2013

Three Surface cracks are among the more common flaws in aircraft and pressure vessel components. Accurate stress intensity analyses and crack growth rate data of surface-cracked components are needed for reliable prediction of their fatigue life and fracture strengths. Three Dimensional finite element method (FEM) was used to obtain the stress intensity factor for surface cracks existing in structures. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Nodes are generated by bucket method, and quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. To examine accuracy and efficiency of the present system, the stress intensity factor for a semi-elliptical surface crack in cylindrical structures subjected to pressure is calculated. Analysis results by present system showed good agreement with those by ASME equation and Raju-Newman's equation.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Economic and Environmental Geology
• /
• v.28 no.4
• /
• pp.351-367
• /
• 1995

The elemental geochemistry and oxygen isotopes of illite/smectite (I/S) have been studied in relationship to the mineralogical trend in the Reindeer D-27 well, Beaufort-Mackenzie Basin. The increase in concentrations of $K_2O$, Rb and rare earth elements (REE), the decrease in concentrations of tetrahedral elements such as Mg, Ti, Sc, Zn and Zr, and the increase in concentrations of tetrahedral elements such as Be and V can be related to I/S compositions that vary systematically with depth. Layer formulae of S- and I-layers are estimated as $[Al_{1.57}Fe_{.19}Mg_{.31}Ti_{.07}][Si_{3.84}Al_{.16}]O_{10}(OH)_2$ and $[Al_{1.84}Mg_{.16}][Si_{3.33}Al_{.67}]O_{10}(OH)_2$, respectively. The mobilization of REE appears to occur during illitization. The increase in concentrations of REE, especially La and Ce, with depth is probably linked to incorporation of ions with high valency (e.g. $V^{5+}$) in tetrahedral sites. The excess valency due to V is partly counter-balanced by ions with low valency (e.g. $Be^{2+}$) and, in turn, the local valency deficiency caused by $Be^{2+}$ could be compensated by high-charge interlayer cations such as REE (＋3). ${\delta}^{18}O$ values of I/S range from 2.91 to 15.72‰ (SMOW), and increase with depth, contrasting to trends observed in the Gulf Coast and elsewhere. The increase in ${\delta}^{18}O$ of I/S results from the rapid increase in ${\delta}^{18}O$ of pore water that overcomes the decrease in temperature-dependent fractionation values with increasing burial depth (${\delta}^{18}O_{pore\;water}>-d{\Delta}/_{I/S-water};\;d{\delta}^{18}O_{I/S}>0$). Calculated ${\delta}^{18}O$ values of pore water in equilibrium with I/S suggest that the original water was probably meteoric water. The stratification of pore water is postulated from the presence of an isotopically light interval, about 450m thick. The depth range of the isotopically light zone overlaps, but does not coincide with the interval of lowered I-content and $K_2O$ concentrations, suggesting that oxygens may have been exchanged independently of mineralogical and geochemical reactions.