• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.619-628
• /
• 2002

This paper deals with the FE analysis for the free vibration of sloshing in horizontal cylindrical tank with baffles. We use Laplace equation based on potential theory as governing equation. This problem is solved by FEM using lineal isoparametric elements. We assume that the tank as well as baffles is rigid body and by separating nodes into two at the baffle location, baffle effect is obtained by separating nodes into two at the baffle location. For the calculation of natural frequencies and mode shapes, we introduce Lanczos transformation and Jacobi iteration methods. Numerical results of the first longitudinal and transverse modes, while comparing with literature cited, are very good. In order for the baffle effects on the free vibration of sloshing, various combinations of baffle parameters, which are location, inner diameter and number, are examined.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• The korean journal of orthodontics
• /
• v.31 no.4 s.87
• /
• pp.425-438
• /
• 2001

The delivery of optimal orthodontic treatment is greatly influenced by clinician's ability to predict and control tooth movement by applying well-known force system to dentition. It is very important to determine the location of the centers of resistance of a tooth or teeth in order to have better understanding the nature of displacement characteristics under various force levels. In this study, three dimensional finite element analysis was used to measure the initial displacement of the consolidated teeth under loading. The purpose of this study was to define the location of the centers of resistance at the upper six anterior segment. To observe the changes of six anterior segment, 200gm, 250gm, 300gm, and 350gm forces at right and left hand side each were imposed toward lingual direction. For this study, two cases, six anterior teeth and six anterior teeth after corticotomy, were reviewed. In addition, it was reviewed the effects of changes on the location of the center of resistance in both cases based on different degree of forces aforementioned. The results were that : 1. The instantaneous center of resistance for the six anterior teeth was vertically located between level 4 and level 5, which is, at 6.76mm, $44.32\%$ apical to the cementoenamel junction level. 2. The instantaneous center of resistance for the six anterior teeth after corticotomy was located vertically between level 4 and level 5, that is, at 7.09mm $46.38\%$ apical to the cementoenamel junction level. 3. Changes of force showed little effect on the location of the center of resistance in each case. 4. It was observed that the location of the instantaneous center of resistance for the six anterior teeth after corticotomy was changed more than the six anterior teeth without corticotomy to the apical part, and the displacement of the consolidated anterior teeth moved further in case of the consolidated teeth after corticotomy.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.189-196
• /
• 2009

This paper presents the subparametric finite element model formulated by partial-linear layerwise theory for the analysis of laminate composites. The proposed model is based on refined approximations of two dimensional plane for orthotropic thick laminate plate as well as thin case. Three dimensional problem can be reduced to two dimensional case by assuming piecewise linear variation of in-plane displacement and a constant value of out-of-plane displacement across the thickness. The integrals of Legendre polynomials are chosen to define displacement fields and Gauss-Lobatto numerical integration is implemented in order to directly obtain maximum values occurred at the nodal points of each layer without other extrapolation techniques. The validity and characteristics of the proposed model have been tested by using orthotropic multilayered plate problem as compared to the values available in the published references. In this study, the convergence test has been carried out to determine the optimal layer model in terms of central deflection and stresses. Also, the distribution of displacements and stresses across the thickness has been investigated as the number of layer is increased.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.458-467
• /
• 2016

In this paper, a strongly coupled method for investigating the interaction between a cable and tow-fish is presented. The nodal position finite element method was utilized to analyze the nonlinear cable dynamics, and 6DOF equations of motion were employed to describe the 3D rigid body motion of the tow-fish. Combining cable and tow-fish systems into a single formulation allowed the two nonlinear systems to be strongly coupled into a unified nonlinear system. This strongly coupled system was numerically integrated in the time domain using a predictor/multi-corrector Newmark algorithm. To demonstrate the validity, efficacy, and applicability of the current approach, two different scenarios (virtual and sea trial) were simulated, and the simulation results were validated using the physical plausibility and the sea trial test.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• The Journal of the Acoustical Society of Korea
• /
• v.16 no.1
• /
• pp.16-23
• /
• 1997

As analysis of the forced dynamic response and vibrational mode for the cantilevered beam is described. Experimental results are compared with the natural frequencies and vibrational modes for the cantilevered beam using the theory of Bernoulli-Euler and finite element method. We have found 1st and 2nd resonance frequency of the cantilevered beam by means of the various external frequencies, $1{\sim}70Hz$, using magnetic transducer. And we have studied the vibrational displacement at obtained resonance frequency of the cantilevered beam. The experimental results for the nodes of cantilevered beam were 0 in 1st mode and 0,0.786 in 2nd mode. close agreement between the theoretically predicted results and experimental result was obtained for the vibrational mode.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2012

In this study, the ground motion was synthesized using the finite fault model by the stochastic green function method, and the difference in the ground motions was evaluated by using various values of the source parameters. An earthquake with a moment magnitude of 6.5 was assumed for the example fault model. The distribution of the slip in the fault plane was calculated using the statistical data of the asperity area. The source parameters considered in this study were the location of the hypocenter in the fault plane and the ratio of the rupture to the shear wave velocity, the rise time, the corner frequency of the source spectrum, and a high frequency filter. The values of the parameters related to the stochastic element source model were adjusted for different tectonic regions, and the others were selected for several possible cases. The response spectra were constructed from the synthesized ground motion time history and compared with the different parameter values. The frequency range affected by each parameter and the differences of the spectral accelerations were evaluated.

• The Journal of Korean Academy of Prosthodontics
• /
• v.47 no.4
• /
• pp.394-405
• /
• 2009

Statement of problem: Crestal bone loss, a common problem associated with dental implant, has been attributed to excessive bone stresses. Design of implant's transgingival (TG) part may affect the crestal bone stresses. Purpose: To investigate if concavely designed geometry at a dental implant's TG part reduces peri-implant bone stresses. Material and methods: A total of five differently configured TG parts were compared. Base model was the ITI one piece implant (Straumann, Waldenburg, Switzerland) characterized by straight TG part. Other 4 experimental models, i.e. Model-1 to Model-4, were designed to have concave TG part. Finite element analyses were carried out using an axisymmetric assumption. A vertical load of 50 N or an oblique load of 50 N acting at $30^{\circ}$ with the implant's long axis was applied. For a systematic stress comparison, a total of 19 reference points were defined on nodal points around the implant. The peak crestal bone stress acting at the intersection of implant and crestal bone was estimated using regression analysis from the stress results obtained at 5 reference points defined along the mid plane of the crestal bone. Results: Base Model with straight configuration at the transgingival part created highest stresses on the crestal bone. Stress level was reduced when concavity was imposed. The greater the concavity and the closer the concavity to the crestal bone level, the less the crestal stresses. Conclusion: The transgingival part of dental implant affect the crestal bone stress. And that concavely designed one may be used to reduce bone stress.