• 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 Korean Society of Steel Construction
• /
• v.23 no.4
• /
• pp.429-438
• /
• 2011

In this study, the explicit numerical algorithm was proposed to simulate the stress erection process and ultimate-load analysis of the strarch (stressed arch) system. The strarch system is a unique and innovative structural system and member prestress comprising prefabricated plane truss frames erected through a post-tensioning stress erection procedure. The flexible bottom chord, which has sleeve and gap details, is closed by the reaction force of the prestressing tendon. The prestress imposed on the tendon will enable the strarch system to be erected. This post-tensioning process is called "stress erection process." During this process, plastic rigid-body rotation occurs to the flexible top chord due to the excessive amount of plastic strain, and the structural characteristic is unstable. In this study, the dynamic relaxation method (DRM) was adopted to calculate the nonlinear equilibrium equation of the system, and a displacement-based finite-element-formulated filament beam element was used to simulate the nonlinear behavior of the top chord sections of the strarch system. The section of the filament beam element was composed by the amount of filaments, which can be modeled by various material models. The Ramberg-Osgood and bilinear kinematic elastic plastic material models were formulated for the nonlinear material behaviors of the filaments. The numerical results that were obtained in the present study were compared with the experiment results of the stress erection and with the results of the ultimate-load analysis of the strarch unit frame. The results of the present studies are in good agreement with the previous experiment results, and the explicit DRM enabled the analysis of the post-buckling behaviors of the strarch unit frame.

• The korean journal of orthodontics
• /
• v.34 no.1 s.102
• /
• pp.93-107
• /
• 2004

Enlow's counterpart analysis explains the complex with anatomic and developmental characteristics where craniofacial aspect of Individuals has been developed. Counterpart analysis does not compare individual measurement with the normal value from the average of majority but analyzes by comparison of values that each individual has. In this study we examined surgical changes in skeletal Class III malocclusion patients(male 40, female 40) and compared them with normal occlusion patients using counterpart analysis. The results indicated that : 1. Skeletal anterior-posterior discrepancy was relieved by shortening of the ramus width(B3). 2. The ramus alignment(R3, R4) was displaced posteriorly and the occlusal plane angle(R5) was rotated clockwise. 3. Skeletal Class III pattern was relieved in the post-operative group, but differences in the level of the cranium(R1, R2) was remaining compared to the normal occlusion patients. 4. In the comparison of surgery methods, the two-jaw surgery group presented changes In the maxillary length(A4), ramus alignment(R3, R4) and occlusal plane angle(R5) compared to the one-jaw surgery group, but the differences were not significant. In the past study about Korean skeletal Class m patients, the skeletal characteristics are upward backward rotation of the cranial base, posterior displacement of the maxilla, forward inclination of the ramus and lengthening of the mandibular body, but in this study, skeletal Class m pattern was relieved by shortening of the ramus width and maxillary advancement by orthognathic surgery, because orthognathic surgery is usually performed on limited areas in the maxilla and the mandible.

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

• Wind and Structures
• /
• v.11 no.4
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.2
• /
• pp.23-38
• /
• 1990

In this study, decompositive optimization method of two levels was selected to optimize effectively the geometry of the truss which takes the multi-loading condition, and the allowable stress, bucking stress, displacement and natural frequency constraints into consideration. The algorithm of this study is made up of sectional optimization using the feasible direction method in level 1, and geometrical optimization employing Powell's one-direction search method which menimizes only objictive function in level 2. The results of this study acquired by beenning applied to structural model of the truss are as follows : 1. It is verified that the algorithm of this study effectively converges, independent of the initial geometry of the truss and the applied various constraints. 2. The optimum goemetry of the truss varies more considerably according to the constraints selected. 3. Under the condition of the same design, the weight of the truss can be decreased more considerably by means of optimizing even the geometry of truss than by means of optimizing the section of truss while fixing geometrical configuration of it, even though there might be a little difference according to the initial geometry of the truss and the design condition.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.67-75
• /
• 2000

The numerical simulation of Brazilian fracture toughness test is carried out using PFC code and the influence parameters are analyzed such as shape of loading plane, size of Brazilian disc and unit panicle of model, loading angle and loading rate. The flattened Brazilian disc is adopted for applying uniform load. The range of loading angle(2$\alpha$) necessary to induce the tensile crack at disc center and to obtain the load-displacement curve giving the critical load for the stable crack propagation is shown as 20$^{\circ}$~40$^{\circ}$. In condition that the loading angle is 20$^{\circ}$, the mode-I fracture toughness is evaluated almost constant in the range of particle size less than I mm and loading rate less than 0.01 mm/s. This range of influence parameters seems appropriate condition for the tensile crack initiation at disc center and the control of stable crack propagation, which can give the reliance in evaluation of fracture toughness by Brazilian test.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.320-328
• /
• 2000

The numerical simulation of Brazilian fracture toughness test is carried out using PFC code and the influence parameters are analyzed such as shape of loading plane, size of Brazilian disc and unit particle of model, loading angle and loading rate. The flattened Brazilian disc is adopted for applying uniform load. The range of loading angle(2$\alpha$) necessary to induce the tensile crack at disc center and to obtain the load-displacement curve giving the critical load for the stable crack propagation is shown as 20°∼40°. In condition that the loading angle is 20°, the mode-I fracture toughness is evaluated almost constant in the range of particle size less than 1 mm and loading rate less than 0.01㎜/s. This range of influence parameters seems appropriate condition for the tensile crack initiation at disc center and the control of stable crack propagation, which can give the reliance in evaluation of fracture toughness by Brazilian test.

• 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 earth science society
• /
• v.43 no.1
• /
• pp.151-164
• /
• 2022

According to the records of historical and instrumental earthquakes, the southeastern part of the Korean Peninsula is considered the highest seismic activity area. Owing to recent reports of numerous Quaternary faults along the Yangsan and Ulsan fault zones, paleoseismological studies are being actively conducted in these areas. The study area is located in the central part of the Ulsan fault zone, where the largest number of active faults have been reported. Based on lineament and geomorphic analysis using LiDAR images and aerial photographs, fault-related landforms showing topographic relief were observed and a trench survey was conducted. The trench length 20 m, width 5 m, depth 5 m is located approximately 300 m away to the northeast from the previously reported Malbang fault. From the trench section, we interpreted the geometric and kinematic characteristics of the fault based on the deformed features of the Quaternary sedimentary layers. The attitude of the reverse fault, N26°W/33°NE, is similar to those of the reported faults distributed along the Ulsan fault zone. Although a single apparent displacement of approximately 40 cm has been observed, the true displacement could not be calculated due to the absence of the slickenline on the fault plane. Based on the geochronological results of the cryogenic structure proposed in a previous study, the most recent faulting event has been estimated as being earlier than the late Wurm glaciation. We interpreted the thrust fault system of the study area as an imbrication structure based on the previous studies and the fault geometry obtained in this additional trench. Although several previous investigations including many trench surveys have been conducted, they found limited success in obtaining the information on fault parameters, which could be due to complex characteristics of the reverse fault system. Additional paleoseismic studies will contribute to solving the mentioned problems and the comprehensive fault evolution.