• 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.

• Steel and Composite Structures
• /
• v.6 no.6
• /
• pp.459-477
• /
• 2006

A theoretical research project is undertaken to develop integrated analysis and design tools for long span composite beams in modern high-rise buildings, and it aims to develop non-linear finite element models for practical design of composite beams. As the first paper in the series, this paper presents the development study as well as the calibration exercise of the proposed finite element models for simply supported composite beams. Other practical issues such as continuous composite beams, the provision of web openings for passage of building services, the partial continuity offered by the connections to columns as well as the behaviour of both unprotected and protected composite beams under fires will be reported separately. In this paper, details of the finite elements and the material models for both steel and reinforced concrete are first described, and finite element studies of composite beams with full details of test data are then presented. It should be noted that in the proposed finite element models, both steel beams and concrete slabs are modelled with two dimensional plane stress elements whose widths are assigned to be equal to the widths of concrete flanges, and the flange widths and the web thicknesses of steel beams as appropriate. Moreover, each shear connector is modelled with one horizontal spring and one vertical spring to simulate its longitudinal shear and pull-out actions based on measured load-slippage curves of push-out tests of shear connectors. The numerical results are then carefully analyzed and compared with the corresponding test results in terms of load mid-span deflection curves as well as load end-slippage curves. Other deformation characteristics of the composite beams such as stress and strain distributions across the composite cross-sections as well as distributions of shear forces and slippages in shear connectors along the beam spans are also examined in details. It is shown that the numerical results of the composite beams compare well with the test data in terms of various load-deformation characteristics along the entire deformation ranges. Hence, the proposed analysis and design tools are considered to be simple and yet effective for composite beams with practical geometrical dimensions and arrangements. Structural engineers are strongly encouraged to employ the models in their practical work to exploit the full advantages offered by composite construction.

• Journal of the Korea Concrete Institute
• /
• v.14 no.6
• /
• pp.892-899
• /
• 2002

Among the methods for enhancement of load-carrying capacity on flexural concrete member, recently, a concept is being investigated which replaces the steel in a conventional reinforced concrete member with a fiber reinforced polymer(FRP) shell. This study focuses on modeling of the structural behavior of concrete surrounded with FRP shells in flexural bending members. A numerical model of fiber cross-sectional analysis is proposed to predict the stress and deformation state of the FRP shell and concrete. The stress-strain relationship of concrete confined by a FRP shell is formulated to be based on the constitutive law of concrete in multi-axial compressive stress state, in assuming that the compression response is dependent on the radial expansion of the concrete. To describe the FRP shell behavior, equivalent orthotropic properties of in-plane behavior from classical lamination theory are used. The present model is validated to compare with the experiments of 4-point bending tests of FRP shell concrete beam, and has well predicted the moment-curvature relationships of the members, axial and hoop strains in the section, and the enhancement of confinement effect in concrete surrounded by FRP shell.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.463-474
• /
• 2009

Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.24 no.3
• /
• pp.259-267
• /
• 2004

A finite element method is used to simulate interaction of laser-based ultrasounds with surface breaking tracks in elastic media. The laser line source focused on the surface of semi-infinite medium is modeled as a shear dipole in 2-D plane strain finite elements. The shear dipole-finite clement model is found to give correct directivity patterns for generated longitudinal and shear waves. The interaction of surface waves with surface breaking cracks (2-D machined slot) is considered in two ways. Both the source and receiver are fixed with respect to the cracks in the first case, while the source is moving in another case. It is shown that the crack depth tested in the range of 0.3-5.0mm $({\lambda}_R/d=0.21{\sim}3.45)$ can be measured using the corner reflected waves produced by the fixed laser source. The moving laser source is found to cause a large amplitude change of reflected waves near crack, and the crack whose depth is one order lower than the wavelength ran be detected from this change.

• Economic and Environmental Geology
• /
• v.15 no.1
• /
• pp.41-47
• /
• 1982

Typical stratigraphic sequences of Paleozoic form the Baegunsan Syncline (equal to Hambaeg Syncline) and Taebaegsan Group is the basement in the southern part of Kangweon Province. Deformation pattern depends on their stratigraphic site and their internal properties. In the biotite schists and meta-sandstones of Pre-Cambrian sequences, flow folds and ptygmatic folds are developed, and high strain deformation is pervasive. Deformational patterns of Cambro-Ordovician are variable because of their various formations. Fracturing is pervasive in the Jangsan Quartzite and Daegi and Maggol Limestone. Welldeveloped slaty cleavage and angular folds of kink and chevron types are the prevailing structures of pelitic rocks. The most characteristic feature of limestone alternated with argillite is the "Compositional cleavage (Author's proposal)" known as the "Worm-eatern" structures. It was known that this structures have a sedimentary origin. But their preferred orienation of long axis of erosional lime holes, originated from tectonic deformation of folding. And this structures have the same character as the axial plane cleavages. Fracturing and long wave concentric folds are dominant characters of the Upper Paleozic sequences. In this area, two folds are superimposed. Field studying and analysis of structures show that the cross-folds of NS-direction are pervasive and were deformed by the EW-directional Baegunsan Syncline. The cross-folding is the major and penetrative deformation and prior to Baegunsan Syncline.

• Geotechnical Engineering
• /
• v.10 no.2
• /
• pp.41-56
• /
• 1994

Under the groundwater level, the tunnel face is affected by the seepage force and the groundwater flow may cause a trouble to the tunnel support systems. The appropriate methods of analysis and design in the tunnel face and the lining, considering groundwater flow according to tunnel drainage condition are presented in this thesis. First, the effect of seepage on the stability of tunnel face was studied. Seepage force was estimated by the 3-D finite element analysis and the stability of tunnel face was checked by analytical method. Furthermore, using the finite difference method the stress and displacement on the face were computed for either case, where the seepage force is or is not considered, and the effect of seepage on the tunnel face stability was evaluated. Second, the effect of seepage force on the tunnel lining when construction is finished and steady state seepage flow occurs was studied and a design methodology considering seepage effect was made. Consequently, in case where the groundwater level remains almost unchanged and the steady state groundwater flow occurs, the proper countermeasures for face staility are required according to the condition of groundwater flow. Moreover, the tunnel lining should be designed and constructed considering the seepage force occuring by the groundwater flow toward the tunnel linings.

• Journal of the Korean Geosynthetics Society
• /
• v.16 no.4
• /
• pp.83-92
• /
• 2017

This paper deals with numerical analysis of behavior of curved mechanically stabilized earth(MSE) walls with geosynthetics reinforcement. Unlike typical concrete retaining walls, MSE wall enables securing stability of higher walls without being constrained by backfill height and is currently and widely used to create spaces for industrial and residential complexes. The design of MSE walls is carried out by checking external stability, similarly to the external checks of conventional retaining wall. In addition, internal stability check is mandatory. Typical stability check based on numerical analysis is done assuming 2-dimensional condition (plane strain condition). However, according to the former studies of 3-dimensional MSE wall, the most weakest part of a curved geosynthetic MSE wall is reported as the convex location, which is also identified from the studies of the laboratory model tests and field monitoring. In order to understand the behaviour of the convex location of the MSE wall, 2-dimensional analysis clearly reveals its limitation. Furthermore, laboratory model tests and field monitoring also have restriction in recognizing their behaviour and failure mechanism. In this study, 3-dimensional numerical analysis was performed to figure out the behaviour of the curved part of the geosynthetic reinforced wall, and the results of the straight-line and curved part in the numerical analysis were compared and analysed. In addition, the behaviour characteristics at each condition were compared by considering the overburden load and relative density of backfill.

• 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.

• Journal of the Korean Chemical Society
• /
• v.40 no.3
• /
• pp.161-166
• /
• 1996

6-ethyl-5,6-dihydrouracil($C_6H_10N_2O_2$) is monoclinic, space group $$P2_{1}c}$$ with a=10.302(2), b=10.419(3), $c=7.095(1)\AA$, $\beta=106.6(0)$, Z=4, $V=729.7(3)\AA$^3$$, $D_c=1.29 g/cm^3,\;{\lambda}(MoK\alpha)=0.71073\AA$, $\mu=0.010cm^{-1}$, F(000)=304, and R=0.054 for 1070 unique observed reflection with F>4.0 $\sigma(F).$ The structure was solved by direct methods and refined by full-matrix least-squares refinement with the fixed C-H bond length at $0.96\AA.$ The hydrouracil molecule makes an envelope conformation with the ethyl substituent oriented to an axial position attainable to a varying degree of steric strain. There are two intermolecular hydrogen-bondings via N-H---O interactions, being nearly parallel to the 100 plane. The shortest distance between molecules is $3.187\AA$ of C(4) and O(8) (-x,-y, 1-z).