• Journal of Biomedical Engineering Research
• /
• v.28 no.3
• /
• pp.397-409
• /
• 2007

Owing to needs of biomechanical comprehension and analysis to obtain various medical treatment designs which are related with the spine in order to cure and diagnose LBP patients, the FE modeling and nonlinear analysis of lumbosacrum including a partial ilium and iliolumbar ligaments, were carried out. First, we investigated whether the geometrical configuration of vertebrae displayed by DICOM slice files is regular and normal condition. After constructing spinal vertebrae including a partial ilium, a sacrum and five lumbars (from L1 to L5)with anatomical shape reconstructed using softwares such as image modeler and CAD modeler, we added iliolumbar ligaments, lumbar ligaments, discs and facet joints, etc.. And also, we assigned material property and discretized the model using proper finite element types, thus it was completely modeled through the above procedure. For the verification of each segment, average sagittal ROM, average coronal ROM and average transversal ROM under various loading conditions(${\pm}10Nm$), average vertical displacement under compression(400N), ALL(Anterior Longitudinal Ligament) and PLL(Posterior Longitudinal Ligament) force at L12 level, strains of seven ligaments on sagittal plane at L45 level and maximal strain of disc fibers according to various loading conditions at L45 level, etc., they were compared with experimental results. For the verification of multilevel-lumbosacrum spine including partial ilium and iliolumbar ligaments, the cases with and without iliolumbar ligaments were compared with ROM of experiment. The results were obtained from analysis of the verified FE model as follows: I) Iliolumbar ligaments played a stabilizing role as mainly posterior iliolumbar ligaments under flexion and as both posterior and anterior iliolumbar ligaments of one side under lateral bending. 2) The iliolumbar ligaments decreased total ROM of 1-8% in total model according to various motion conditions, which changed facet contact forces of L5S level by approximately 0.8-1.4 times and disc forces of L5S level by approximately 0.8-1.5 times more than casewithout ilioligaments, under various loading conditions. 3) The force of lower discs such as L45 and L5S was bigger than upper discs under flexion, left and right bending and left and right twisting, except extension. 4) It was predicted that strains of posterior ligaments among iliolumbar ligaments would produce the maximum 16% under flexion and the maximum 10% under twisting. 5) It's expected that this present model applies to the development and design of artificial disc, since it was comparatively in agreement with the experimental datum.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary element.

• Korean Journal of Materials Research
• /
• v.3 no.3
• /
• pp.253-260
• /
• 1993

High strength and fracture toughness of Al-Li-Cu alloy(2090 Al alloy) have been achieved by the improvement of melting and casting, extrusion and heat treatment techniques. To establish the sucessful process for semi-industrial scale ingot(20Kg) the following areas have been investigated: (1) Improvement of melting and casting techniques for ingot by introducing atmospheric modifications, vacuum and rotary degassing, and deslagging. (2) The effect of heat treatment on mechanical properties (3) Mechanical characterization by tensile test, fracture toughness test and fatigue crack propagation test. High mechanical properties were found to be intimately related with ingot soundness. Tensile strength of final products varied from 534MPa to 566MPa in peak aged condition while elongation/ductility ranged from 9.0% to 11.9%. From the fracture toughness test with using compact tensile specimen, plane strain fracture toughness($K_{Ic}$) appeared to be 39MPa${\surd}$m in peak aged condition and 23MPa${\surd}$ m in underaged condition. When load ratios of 0.1, 0.3 and 0.5 were given ${\Delta}K_{th}$ was 6.0MPa${\surd}$ m, 5.3MPa${\surd}$ m and 4.3MPa${\surd}$ m respectively.

• Korean Journal of Materials Research
• /
• v.11 no.2
• /
• pp.88-93
• /
• 2001

The thermal behavior and the crystallographic characteristics of an epitaxial $C49-TiSi_2$ island formed in a Si (001) substrate by $N_2$, treatment were investigated by X-ray diffraction (XRD) and high-resolution transmission electron microscopy (HRTEM). It was found from the analyzed results that the epitaxial $C49-TiSi_2$ was thermally stable even at high temperature of $1000^{\circ}C$ therefore did not transform into the C54-stable phase and did not deform morphologically. HRTEM results clearly showed that the epitaxial $TiSi_2$ phase and Si have the orientation relationship of (060)[001]$TiSi_2$//(002)[110]Si, and the lattice strain energy at the interface was mostly relaxed by the formation of misfit dislocations. Furthermore, the mechanism on the formation of the epitaxial $_C49-TiSi2$ in Si and stacking faults lying on the (020) plane of the C49 Phase were discussed through the analysis of the HRTEM image and the atomic modeling.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.2
• /
• pp.67-75
• /
• 2019

The diaphragm used for CFT columns has a small amount of steel to be used, but has a disadvantage that welding is difficult and openings are required because the steel tube and four sides must be welded. The improved diaphragm to be examined in this study was cut into four corners by cutting the center hole for concrete filling. In the improved diaphragm, the width of the center hole is the same as that of the previous diaphragm, but the width of the diaphragm contacting the steel tube is reduced, thereby reducing the welding length by about 70% compared to the previous diaphragm. The in-plane strain of each specimen was analyzed when the same load was applied to the interior diaphragm through a simple tensile test. Using the general FEM program(ANSYS 19.2), the analysis was performed under the same conditions as the actual simple tensile test, and the load transfer between the improved diaphragm and the previous diaphragm was compared. When the width of the diaphragm is equal to or smaller than the flange width, stress is concentrated from the end of the diaphragm, and when the flange width is larger, stress is concentrated at the center.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.11
• /
• pp.753-764
• /
• 2020

In this study, the plastic zone and internal earth pressure of the tunnel were calculated using the following three methods: metal plasticity to analyze the deformation of metal during plastic processing, Terzaghi's earth pressure theory from the geotechnical perspective and modified Terzaghi's earth pressure theory, and slip line theory using Mohr-Coulomb yield conditions. All three methods are two-dimensional mathematical analysis models for analyzing the plane strain conditions of isotropic materials. Using the theory of metallurgical plastics, the plastic zone and the internal earth pressure of the ground were obtained by assuming that the internal pressure acts on the tunnel, so different results were derived that did not match the actual tunnel site, where only gravity was applied. An analysis of the plasticity zone and earth pressure via the slip-line method showed that a failure line is formed in a log-spiral, which was found to be similar to the real failure line by comparing the results of previous studies. The earth pressure was calculated using a theoretical method. Terzaghi's earth pressure was calculated to be larger than the earth pressure considering the dilatancy effect.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.1
• /
• pp.49-63
• /
• 2014

Bearing capacity factor $N_{\gamma}$ and shape factor were studied for rigid strip and circular footings with a rough base on sand by numerical modelling considering the effect of dilation angle. The numerical model was developed with an explicit finite difference code. Loading procedures and interpretation methods were devised in order to shorten the running time while eliminating the exaggeration of the reaction caused by the explicit scheme. Using the Mohr-Coulomb plasticity model with associated (${\psi}={\phi}$) and nonassociated (${\psi}$ < ${\phi}$) flow-rules, the bearing capacity factor $N_{\gamma}$ was evaluated for various combinations of internal friction angles and dilation angles. Bearing capacity factor decreased as the dilation angle was reduced from the associated condition. An equation applicable to typical sands was proposed to evaluate the relative bearing capacity for the nonassociated condition compared to the associated condition on which most bearing capacity factor equations are based. The shape factor for the circular footing varied substantially when the plane-strain effect was taken into account for the strip footing. The numerical results of this study showed closer trends with the previous experimental results when the internal friction angle was increased for the strip footing. Discussions are made on the reason that previous equations for the shape factor give different results and recommendations are made for the appropriate design shape factor.

• Economic and Environmental Geology
• /
• v.34 no.1
• /
• pp.133-146
• /
• 2001

We report the results of structural field observation and measurement of anisotropy of magnetic susceptibility (AMS) of the diamictitic Hwanggangri Formation distributed in Chungju-Suanbo area of the Okchon Belt, Korea. The outcrops of the Hwanggangri Formation show two types of cleavage in general: slaty cleavage (SI) and crenulation cleavage (5z). 5] cleavage is, however, well observable only in the notheastem (NE) part of study area, while overwhelmed by 52 cleavage in the southwestern (5W) part, indicating stronger later deformation in 5W part of the study area. This partitioning of the study area is corroborated by both IRM and AMS parameters: NE part of the study area is characterized by higher IRM intensity, higher bulk magnetic susceptibility, higher AM5 degree, and by oblate shape of magnetic susceptibility ellipsoid. Their values become drastically lowered toward southwest, and reach to a stable minimum in the whole 5W part of the study area. In addition, degree of both metamorphism and deformation tends to increase gradually from northeast toward southwest and also from northwest toward southeast in the study area. Based on the distribution pattern of the principal axes ( $k_1, k_2, k_3$ axes) of magnetic anisotropy ellipsoids revealed in the NE part of the study area, three episodes of deformation ( $D_1, D_2, D_3$ ) are recognized: D_1$ deformation produced $S_2$ cleavage with NE-5W trend, which is caused by a strong NW-SE tlattening of a coaxial pure shear. $D_2$ deformation produced 5z cleavage characterized by a non-coaxial deformation. It was caused by a ductile or semi-ductile thrusting toward NW and concurrent sinistral shearing along $S_2$ cleavage plane. Lastly, $D_3$ deformation produced tlexural folding of all previous structures with a nearly horizontal NE fold axis. Distribution pattern of the principal axes of magnetic anisotropy ellipsoid from the SW part of the study area, on the other hand, does not show any coherency among sites or samples. We interpret that this dispersed pattern of $k_1, k_2, k_3$ axes together with lower anisotropy strength indicates that magnetic fabrics in the SW part have been disturbed either by a superposition of strong deformation/metamorphism or by a kind of reciprocal strain due to an overlapping of $D_1$ and $D_2$ or by both processes.

• Korean Journal of Remote Sensing
• /
• v.35 no.6_1
• /
• pp.933-944
• /
• 2019

This study analyzes the tide deformation of land boundary regions on the east (Region A) and west (Region B) sides of the Ross Ice Shelf in Antarctica using Double-Differential Interferometric Synthetic Aperture Radar (DDInSAR). A total of seven Sentinel-1A SAR images acquired in 2015-2016 were used to estimate the accuracy of tide prediction model and Young's modulus of ice shelf. First, we compared the Ross Sea Height-based Tidal Inverse (Ross_Inv) model, which is a representative tide prediction model for the Antarctic Ross Sea, with the tide deformation of the ice shelf extracted from the DDInSAR image. The accuracy was analyzed as 3.86 cm in the east region of Ross Ice Shelf and it was confirmed that the inverse barometric pressure effect must be corrected in the tide model. However, in the east, it is confirmed that the tide model may be inaccurate because a large error occurs even after correction of the atmospheric effect. In addition, the Young's modulus of the ice was calculated on the basis of the one-dimensional elastic beam model showing the correlation between the width of the hinge zone where the tide strain occurs and the ice thickness. For this purpose, the grounding line is defined as the line where the displacement caused by the tide appears in the DDInSAR image, and the hinge line is defined as the line to have the local maximum/minimum deformation, and the hinge zone as the area between the two lines. According to the one-dimensional elastic beam model assuming a semi-infinite plane, the width of the hinge region is directly proportional to the 0.75 power of the ice thickness. The width of the hinge zone was measured in the area where the ground line and the hinge line were close to the straight line shown in DDInSAR. The linear regression analysis with the 0.75 power of BEDMAP2 ice thickness estimated the Young's modulus of 1.77±0.73 GPa in the east and west of the Ross Ice Shelf. In this way, more accurate Young's modulus can be estimated by accumulating Sentinel-1 images in the future.