• Journal of the Korean Geotechnical Society
• /
• v.29 no.3
• /
• pp.29-40
• /
• 2013

Seismic response of single degree of freedom system supported by shallow foundation was analyzed by using nonlinear explicit finite difference element code. Numerical analysis results were verified with dynamic centrifuge test results of the same soil profile and structural dimensions with the numerical analysis model at a centrifugal acceleration of 20 g. Differences between the analysis and the test results induced by the boundary conditions of control points can be reduced by adding additional local damping to the natural born cyclic hysteretic damping of the soil strata. The analysis results show good agreement with the test results in terms of both time histories and response spectra. Thus, it can be concluded that the nonlinear explicit finite difference element code will be a useful technique for estimating seismic residual displacement, earthpressure etc. which are difficult to measure during laboratory tests and real earthquake.

• Journal of Periodontal and Implant Science
• /
• v.36 no.2
• /
• pp.473-488
• /
• 2006

Adequate bone quality and stress distribution to the bone are of decisive importance for implant success. Even though the success rates of dental implants have been high, implant failures do occur. Overloading has been identified as a primary factor behind dental implant failure. The purpose of this study was to theoretically investigate the effect of two types of implants on the stress distribution in poor bone quality. Employing the finite element method, the study modeled a 4.1 mm diameter, 12.0 mm length implant placed in cortical or spongeous bone. A static loading of lOON was applied at the occlusal surface at 0, 30 degrees angle to the vertical axis of the implant. von Mises stresses concentrations in the supporting bone were analyzed with finite element analysis program. The results were as follows; 1. The stresses at the marginal bone were higher under buccal oblique load(30 degrees off of the long axis) than under vertical load. 2. Under buccal oblique load, the stresses were higher at the lingual marginal bone than at the buccal marginal bone, and the differences were almost the same. 3, Under vertical and oblique load, the stress was the highest at the marginal bone and lowest at the bone around apical portions of implant in cortical bone. 4, Under vertical load, Model 1 showed more effective stress distribution than Model 2 irrespective of bone types. On the other hand, Model 2 showed lower stress concentration than Model 1 under buccal oblique load.

• The Journal of Korean Academy of Prosthodontics
• /
• v.44 no.1
• /
• pp.85-102
• /
• 2006

Statement of problem: Currently, there are some 20 different geometric variations in implant/abutment interface available. The geometry is important because it is one of the primary determinants of joint strength, joint stability, locational and rotational stability. Purpose: As the effects of the various implant-abutment connections and the prosthesis height variation on stress distribution are not yet examined this study is to focus on the different types of implant-abutment connection and the prosthesis height using three dimensional finite element analysis. Material and method. The models were constructed with ITI, 3i TG, Bicon, Frialit-2 fixtures and solid abutment, TG post, Bicon post, EstheticBase abutment respectively. And the super structures were constructed as mandibular second premolar shapes with 8.5 mm, 11 mm, 13.5 mm of crown height. In each model, 244 N of vertical load and 244 N of $30^{\circ}$ oblique load were placed on the central pit of an occlusal surface. von Mises stresses were recorded and compared in the crowns, abutments, fixtures. Results: 1. Under the oblique loading, von Mises stresses were larger in the crown, abutment, fixture compared to the vertical loading condition. 2. The stresses were increased proportionally to the crown height under oblique loading but showed little differences with three different crown heights under vertical loading. 3. In the crown, the highest stress areas were loading points under vertical loading, and the finish lines under oblique loading. 4. Under the oblique loading, the higher stresses were located in the fixture/abutment interface of the Bicon and Frialit-2 systems compared to the ITI and TG systems. Conclusions: The stress distribution patterns of each implant-abutment system had difference among them and adequate crown height/implant ratio was important to reduce the stresses around the implants.

• Computers and Concrete
• /
• v.24 no.4
• /
• pp.379-397
• /
• 2019

Steel-concrete composite beams are widely employed in constructions and their performance at the serviceability stage is of concern among practitioners and design regulations. In this context, an accurate evaluation of long-term deflections via various rheological concrete models is needed. In this work, the performance and predict capability of some concrete creep and shrinkage models ACI, CEB, B3, FIB and GL2000 are ascertained, and compared by using statistical bias indicators. Ten steel-concrete composite beams with existing experimental and numerical results are then modeled for this purpose. The proposed modeling technique uses the finite element method, where the concrete slab and steel beam are modeled with shell finite elements. Concrete is considered as an aging viscoelastic material and cracking is treated with the common smeared approach. The results show that when the experimental ultimate shrinkage strain is used for calibration, all studied rheological models predict nearly similar deflections, which agree with the experimental data. In contrast, significance differences are encountered for some models, when none calibration is made prior to. A value between twenty and thirty times the cracking strain is recommended for the ultimate tensile strain in the tension stiffening model. Also, increasing the relative humidity and decreasing the ambient temperature can lead to a substantial reduction of slab cracking for beams under negative flexure. Finally, there is not a unique rheological model that clearly excels in all scenarios.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.3
• /
• pp.335-342
• /
• 1999

The purpose of present research is to develop and efficient numerical method for the calculation of potential flow and predict the wave-making resistance for the application to ship design of tuna purse seiner. Havelock was considered the wave resistance of a post extending vertically downwards through the water from the surface, its section by a horizontal plane being the same at all depths and having its breadth small compared with its length. This enables us to elucidate certain points of interest in ship resistance. However, the ship has not infinite draft. So, the problem which is investigated ind detail in this paper is the wave resistance of a mathematical quadratic model in a uniform stream. The paper deals with the numerical calculation of potential flow around the series 60 with forward velocity by the new slender ship theory. This new slender ship theory is based on the asymptotic expression of the Kelvin-source, distributed over the small matrix at each transverse section so as to satisfy the approximate hull boundary condition due to the assumption of slender body. The numerical results using the panel shift method and finite difference method are compared with the experimental results for wigley mono hull. There are no differences in the wave resistance. However, it costs much time to compute not only wave resistance but also wave pattern over some range of Froude numbers. More improvements are strongly desired in the numerical procedure.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1975-1981
• /
• 1991

A new numerical solution of the elastohydrodynamic lubrication(EHL) problem of an axially profiled cylindrical roller is presented. A finite difference method and the Newton-Raphson method are used to solve the nonlinear system equations. A non-uniform grid system is adopted to reduce the number of grid points and to obtain accurate solution. For two different types of profiles which have similar elastostatic pressure distribution, the EHL results show large differences. Especially the difference in film shape is larger than in pressure distribution. Therefore, the magnitude of the minimum film thickness should be a major criteria to design the axial profile of the roller. Variations of the minimum film thickness with dimensionless parameters show considerably different behavior from those of infinite solution and show a good agreement with the experimental data in literatures. Present numerical scheme can be used generally in the analysis of three-dimensional EHL problem.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.9
• /
• pp.935-942
• /
• 2009

In this paper, the purpose is to investigate the natural frequency of a cracked Timoshenko cantilever beams by FEM(finite element method) and experiment. In addition, a method for detection of crack in a cantilever beams is presented based on natural frequency measurements. The governing differential equations of a Timoshenko beam are derived via Hamilton's principle. The two coupled governing differential equations are reduced to one fourth order ordinary differential equation in terms of the flexural displacement. The crack is assumed to be in the first mode of fracture and to be always opened during the vibrations. The detection method of a crack location in a beam based on the frequency measurements is extended here to Timoshenko beams, taking the effects of both the shear deformation and the rotational inertia into account. The differences between the actual and predicted crack positions and sizes are less than 6 % and 23 % respectively.

• Wind and Structures
• /
• v.9 no.6
• /
• pp.473-492
• /
• 2006

This paper presents numerical results of pedestrian-level wind environment around the base of a row of tall buildings by CFD. Four configurations of building arrangement are computed including a single square tall building. Computed results of pedestrian-level wind flow patterns and wind speeds are compared to previous wind tunnel measurement data to enable an assessment of CFD predictions. The CFD model uses the finite-volume method with RNG $k-{\varepsilon}$ model for turbulence closure. It is found that the numerical results can reproduce key features of pedestrian-level wind environment such as corner streams around corners of upwind building, sheltered zones behind buildings and channeled high-speed flow through a building gap. However, there are some differences between CFD results and wind tunnel data in the wind speed distribution and locations of highest wind speeds inside the corner streams. In locations of high ground-level wind speeds, CFD values match wind tunnel data within ${\pm}10%$.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.1 no.1
• /
• pp.120-130
• /
• 1993

In the analysis of automobile exhaust system, the exciting forces from the engine determine the dynamic behavior of the system and the dynamic characteristics influence the riding quality. Therefore, the identification of the force in numeric value is quite important for the vibrational reduction. However, the value is difficult to obtain by experiments due to harsh conditions around the engine. In this research, an optimization technology is adopted to evaluate the exciting forces. An experimental method is conducted for the verification of the finite element modeling. Displacements on the end of the exhaust system are measured under the idling environment. cost function is set up to minimize the differences between the displacements of the numerical simulation and the experiment. Design variables are the components of the exciting forces. That is, optimization is utilized to estimate the forces with existing data. Excellent estimations have been calculated efficiently and the information is used again for the forced vibration of the exhaust system.