• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.2
• /
• pp.117-131
• /
• 2006

A coupled T-H-M(Thermo-Hydro-Mechanical) analysis was carried out for KENTEX (KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System), which is a facility for validating the coupled T-H-M behavior in the engineered barrier system of the Korean reference HLW(high-level waste) disposal system. The changes of temperature, water saturation, and stress were estimated based on the coupled T-H-M analysis, and the influence of the types of mechanical constitutive material laws was investigated by using elastic model, poroelastic model, and poroelastic-plastic model. The analysis was done using ABAQUS, which is a commercial finite element code for general purposes. From the analysis, it was observed that the temperature in the bentonite increased sharply for a couple of days after heating the heater and then slowly increased to a constant value. The temperatures at all locations were nearly at a steady state after about 37.5 days. In the steady state, the temperature was maintained at $90^{\circ}C$ at the interface between the heater and the bentonite and at about $70^{\circ}C$ at the interface between the bentonite and the confining cylinder. The variation of the water saturation with time in bentonite was almost same independent of the material laws used in the coupled T-H-M processes. By comparing the saturation change of T-H-M and that of H-M(Hydro-Mechanical) processes using elastic and poroelastic material mod31 respectively, it was found that the degree of saturation near the heater from T-H-M calculation was higher than that from the coupled H-M calculation mainly because of the thermal flux, which seemed to speed up the saturation. The stresses in three cases with different material laws were increased with time. By comparing the stress change in H-M calculation using poroelasetic and poroelasetic-plastic model, it was possible to conclude that the influence of saturation on the stress change is higher than the influence of temperature. It is, therefore, recommended to use a material law, which can model the elastic-plastic behavior of buffer, since the coupled T-H-M processes in buffer is affected by the variation of void ratio, thermal expansion, as well as swelling pressure.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

• Geotechnical Engineering
• /
• v.12 no.4
• /
• pp.157-178
• /
• 1996

Current design methods for reinforced earth structures take no account of the magnitude of the strains induced in the tensile members as these are invariably manufactured from high modulus materials, such as steel, where straits are unlikely to be significant. With fabrics, however, large strains may frequently be induced and it is important to determine these to enable the stability of the structure to be assessed. In the present paper internal design method of analysis relating to the use of fabric reinforcements in reinforced earth structures for both stress and strain considerations is presented. For the internal stability analysis against rupture and pullout of the fabric reinforcements, a strain compatibility analysis procedure that considers the effects of reinforcement stiffness, relative movement between the soil and reinforcements, and compaction-induced stresses as studied by Ehrlich 8l Mitchell is used. I Bowever, the soil-reinforcement interaction is modeled by relating nonlinear elastic soil behavior to nonlinear response of the reinforcement. The soil constitutive model used is a modified vertsion of the hyperbolic soil model and compaction stress model proposed by Duncan et at., and iterative step-loading approach is used to take nonlinear soil behavior into consideration. The effects of seepage pressures are also dealt with in the proposed method of analy For purposes of assessing the strain behavior oi the fabric reinforcements, nonlinear model of hyperbolic form describing the load-extension relation of fabrics is employed. A procedure for specifying the strength characteristics of paraweb polyester fibre multicord, needle punched non-woven geotHxtile and knitted polyester geogrid is also described which may provide a more convenient procedure for incorporating the fablic properties into the prediction of fabric deformations. An attempt to define improvement in bond-linkage at the interconnecting nodes of the fabric reinforced earth stracture due to the confining stress is further made. The proposed method of analysis has been applied to estimate the maximum tensions, deformations and strains of the fabric reinforcements. The results are then compared with those of finite element analysis and experimental tests, and show in general good agreements indicating the effectiveness of the proposed method of analysis. Analytical parametric studies are also carried out to investigate the effects of relative soil-fabric reinforcement stiffness, locked-in stresses, compaction load and seepage pressures on the magnitude and variation of the fabric deformations.

• Journal of the Korean GEO-environmental Society
• /
• v.19 no.4
• /
• pp.5-16
• /
• 2018

In the current study, the engineering behaviour of prebored and precast steel pipe piles was examined from a series of full-scale field measurements by conducting static pile load tests, dynamic pile load tests (EOID and restrike tests) and Class-A and Class-C1 type numerical analysis. The study includes the pile load - settlement relations, allowable pile capacity and shear stress transfer mechanism. Compared to the allowable pile capacity obtained from the static pile load tests, the dynamic pile load tests and the numerical simulation showed surprisingly large variations. Overall among these the restrike tests displayed the best results, however the reliability of the predictions from the numerical analysis was lower than those estimated from the dynamic pile load tests. The allowable pile capacity obtained from the EOID tests and the restrike tests indicated 20.0%-181.0% (avg: 69.3%) and 48.2%-181.1% (avg: 92.1%) of the corresponding measured values from the static pile loading tests, respectively. Furthermore, the computed results from the Class-A type analysis showed the largest scatters (37.1%-210.5%, avg: 121.2%). In the EOID tests, a majority of the external load were carried by the end bearing pile capacity, however, similar skin friction and end bearing capacity in magnitude were mobilised in the restrike tests. The measured end bearing pile capacity from the restrike tests were smaller than was measured from the EOID tests. The present study has revealed that if the impact energy is not sufficient in a restrike test, the end bearing pile capacity most likely will be underestimated. The shear stresses computed from the numerical analysis deviated substantially from the measured pile force distributions. It can be concluded that the engineering behaviour of the pile is heavily affected if a slime layer exists near the pile tip, and that the smaller the stiffness of the slime and the thicker the slime, the greater the settlement of the pile.

• The Journal of Korean Academy of Prosthodontics
• /
• v.48 no.3
• /
• pp.224-231
• /
• 2010

Purpose: This study was to investigate how the crestal module design could affect the level of marginal bone stress around dental implant. Materials and methods: A submerged implant of 4.1 mm in diameter and 10 mm in length was selected as baseline model (Dentis Co., Daegu,Korea).A total of 5 experimental implants of different crestal modules were designed (Type I model : with microthread at the cervical 3 mm, Type II model : the same thread pattern as Type I but with a trans-gingival module, Type III model: the same thread pattern as the control model but with a trans-gingival module, Type IV model: one piece system with concave transgingival part, Type V model: equipped with beveled platform). Stress analysis was conducted with the use of axisy mmetric finite element modeling scheme. A force of 100 N was applied at 30 degrees from the implant axis. Results: Stress analysis has shown no stress concentration around the marginal bone for the control model. As compared to the control model, the stress levels of 0.2 mm areas away from the recorded implant were slightly lower in Type I and Type IV models, but higher in Type II, Type III and Type V models. As compared to 15.09 MPa around for the control model, the stress levels were 14.78 MPa, 18.39 MPa, 21.11 MPa, 14.63 MPa, 17.88 MPa in the cases of Type I, II, III, IV and V models. Conclusion: From these results, the conclusion was drawn that the microthread and the concavity with either crestal or trans-gingival modules maybe used in standard size dental implants to reduce marginal bone stress.

• The Journal of Korean Academy of Prosthodontics
• /
• v.48 no.2
• /
• pp.111-121
• /
• 2010

Purpose: This 3D-FEA study was performed to investigate the influence of marginal bone loss pattern around the implant to the stress distribution. Material and methods: From the right second premolar to the right second molar of the mandible was modeled according to the CT data of a dentate patient. Teeth were removed and an implant ($\Phi\;4.0{\times}10.0mm$) was placed in the first molar area. Twelve bone models were created; Studied bone loss conditions were horizontal bone loss and vertical bone loss, assumed bone loss patterns during biologic width formation, and pathologic vertical bone loss with or without cortification. Axial, buccolingual, and oblique force was applied independently to the center of the implant crown. The Maximum von Mises stress value and stress contour was observed and von Mises stresses at the measuring points were recorded. Results: The stress distribution patterns were similar in the non-resorption and horizontal resorption models, but differed from those in the vertical resorption models. Models assuming biologic width formation showed altered stress distribution, and weak bone to implant at the implant neck area seams accelerates stress generation. In case of vertical bone resorption, contact of cortical bone to the implant may positively affect the stress distribution.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.325-334
• /
• 2008

The standard prestressed concrete I-girder bridge (PSC I-girder bridge) is one of the most prevalent types for small and medium bridges in Korea. When determining the member forces in a section to assess the safety of girder in this type of bridge, the general practice is to use the simplified practical equations or the live load distribution factors proposed in design standards rather than the precise analysis through the finite element method or so. Meanwhile, the live load distribution factors currently used in Korean design practice are just a reflection of overseas research results or design standards without alterations. Therefore, it is necessary to develop an equation of the live load distribution factors fit for the design conditions of Korea, considering the standardized section of standard PSC I-girder bridges and the design strength of concrete. In this study, to develop an equation of the live load distribution factors, a parametric analysis and sensitivity analysis were carried out on the parameters such as width of bridge, span length, girder spacing, width of traffic lane, etc. As a result, the major variables to determine the size of distribution factors were girder spacing, overhang length and span length in case of external girders. For internal adjacent girders, the determinant factors were girder spacing, overhang length, span length and width of bridge. For internal girders, the factors were girder spacing, width of bridge and span length. Then, an equation of live load distribution factors was developed through the multiple linear regression analysis on the results of parametric analysis. When the actual practice engineers design a bridge with the equation of live load distribution factors developed here, they will determine the design of member forces ensuring the appropriate safety rate more easily. Moreover, in the preliminary design, this model is expected to save much time for the repetitive design to improve the structural efficiency of PSC I-girder bridges.

• Journal of the Korea Concrete Institute
• /
• v.24 no.3
• /
• pp.293-303
• /
• 2012

Recently, as construction technology improved, concrete structures not only became larger, taller and longer but were able to perform various functions. However, if extreme loads such as impact, blast, and fire are applied to those structures, it would cause severe property damages and human casualties. Especially, the structural responses from extreme loading are totally different than that from quasi-static loading, because large pressure is applied to structures from mass acceleration effect of impact and blast loads. Therefore, the strain rate effect and damage levels should be considered when concrete structure is designed. In this study, the low velocity impact loading test of steel fiber reinforced concrete (SFRC) slabs including 0%~1.5% (by volume) of steel fibers, and strengthened with two types of FRP sheets was performed to develop an impact resistant structural member. From the test results, the maximum impact load, dissipated energy and the number of drop to failure increased, whereas the maximum displacement and support rotation were reduced by strengthening SFRC slab with FRP sheets in tensile zone. The test results showed that the impact resistance of concrete slab can be substantially improved by externally strengthening using FRP sheets. This result can be used in designing of primary facilities exposed to such extreme loads. The dynamic responses of SFRC slab strengthened with FRP sheets under low velocity impact load were also analyzed using LS-DYNA, a finite element analysis program with an explicit time integration scheme. The comparison of test and analytical results showed that they were within 5% of error with respect to maximum displacements.

• Computational Structural Engineering
• /
• v.7 no.3
• /
• pp.95-103
• /
• 1994

The safety related components in the nuclear power plant should be designed to withstand the seismic load. Among these components the integrity of reactor internals under earthquake load is important in stand points of safety and economics, because these are classified to Seismic Class I components. So far the modelling methods of reactor internals have been investigated by many authors. In this paper, the dynamic behaviour of reactor internals of Yong Gwang 1&2 nuclear power plants under SSE(Safe Shutdown Earthquake) load is analyzed by using of the simpled Global Beam Model. For this, as a first step, the characteristic analysis of reactor internal components are performed by using of the finite element code ANSYS. And the Global Beam Model for reactor internals which includes beam elements, nonlinear impact springs which have gaps in upper and lower positions, and hydrodynamical couplings which simulate the fluid-filled cylinders of reactor vessel and core barrel structures is established. And for the exciting external force the response spectrum which is applied to reactor support is converted to the time history input. With this excitation and the model the dynamic behaviour of reactor internals is obtained. As the results, the structural integrity of reactor internal components under seismic excitation is verified and the input for the detailed duel assembly series model could be obtained. And the simplicity and effectiveness of Global Beam Model and the economics of the explicit Runge-Kutta-Gills algorithm in impact problem of high frequency interface components are confirmed.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.