• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.1 no.1
• /
• pp.53-68
• /
• 1981

Most clays under sustained load exhibit time-dependent deformation because of creep movement of soil particles and many investigators have attempted to relate their findings to the creep behavior of natural ground and to the long-term stability of slopes. Since the creep behavior of clays may assume a variety of forms depending on such factors as soil plasticity, activity and water content, it is difficult and complicated to analyse the creep behavior of clays. Rheological models composed of linear springs in combination with linear or nonlinear dashpots and sliders, are generally used for the mathematical description of the time-dependent behavior of soils. Most rheological models, however, have been proposed to simulate the behavior of secondary compression for saturated clays and few definitive data exist that can evaluate the behavior of non-saturated clays under the action of sustained stress. The clays change gradually from a solid state through plastic state to a liquid state with increasing water content, therefore, the rheological models also change. On the other hand, creep is time-dependent, and also the effect of thixotropy is time-function. Consequently, there may be certain correlations between creep behavior and the effects of thixotropy in compacted clays. In addition, the states of clay depend on water content and hence the height of the specimen under drained conditions. Futhermore, based on present and past studies, because immediate elastic deformation occurs instantly after the pressure increment without time-delayed behavior, the factor representing immediate elastic deformations in the rheological model is necessary. The investigation described in this paper, based on rheological model, is designed to identify the immediate elastic deformations and the effects of thixotropy and height of clay specimens with varing water content and stress level on creep deformations. For these purposes, the uniaxial drain-type creep tests were performed. Test results and data for three compacted clays have shown that a linear top spring is needed to account for immediate elastic deformations in the rheological model, and at lower water content below the visco-plastic limit, the effects of thixotropy and height of clay specimens can be represented by the proposed rheological model not considering the effects. Therefore, the rheological model does not necessitate the other factors representing these effects. On the other hand, at water content higher than the visco-plastic limit, although the state behavior of clays is visco-plastic or viscous flow at the beginning of the test, the state behavior, in the case of the lower height sample, does not represent the same behavior during the process of the test, because of rapid drainage. In these cases, the rheological model does not coincide with the model in the case of the higher specimens.

• Journal of Korean Society of Steel Construction
• /
• v.14 no.6
• /
• pp.735-746
• /
• 2002

The task of plastic theory is twofold: first, to set up relationships between stress and strain that adequately describe the observed plastic deformation of metals, and second, to develop techniques for using these relationships in studying of the mechanics of metal forming processes, and the anlaysis and design of structures. One of the major problems in the theory of plasticity is to describe the behavior of work-hardening materials in the plastic range for complex loading histories. This can be achieved by formulating constitutive laws either in the integral or differential forms. To adequately predict the response of steel members during cyclic loading, the hardening rule must account for the features of cyclic stress-strain behavior. Neithe of the basic isotropic and kinematic hardening rules is suitable for describing cyclic streess-strain behavior, although a kinematic hardening rule describes the nearly linear portions of the stabilized hystersis loops. There is also a limited expansion of the yield surface as predicted by the isotropic hardening rule. Strong ground motions or wind gusts affect the complex and nonproportional loading histories in the inelastic behavior of structues rather than the proportional loading. Nonproportional loading is defined as externally applied forces on the structure, with variable ratios during the entire loading history. This also includes the rate of time-dependency of the loads. For nonproportional loading histories, unloading may take place along a chord instead of the radius of the load surface. In such cases, the shape of the stress-strain curve has to be determined experimentally for all non-radial loading conditions. The plasticity models including two surface models ae surveyed based on a yield surface and a bound surface that represent a state of maximum stress. This paper is concerned with the improvement of a plasticity models of the two-surface type for structural steel. This is follwed by an overview of plasticity models on structural steel. Finally the need for further research is identified.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.04a
• /
• pp.108-108
• /
• 2009

The main object of this research is to minimize the shock effects which frequently result in fatal damage in wind met mast on impact of barge. The collision between wind met mast and barge is generally a complex problem and it is often not practical to perform rigorous finite element analyses to include all effects and sequences during the collision. LS-dyna generally purpose explicit finite element code, which is a product of ANSYS software, is used to model and analyze the non-linear response of the met mast due to barge collision. A significant part of the collision energy is dissipated as strain energy and except for global deformation modes, the contribution from elastic straining can normally be neglected. On applying impact force of a barge to wind met mast, the maximum acceleration, internal energy and plastic strain were calculated for each load cases using the finite element method and then compare it, varying to the velocity of barge, with one varying to the thickness of rubber fender conditions. Hence, we restrict the present research mainly to the wind met mast and also parametric study has been carried out with various velocities of barge, thickness of wind met mast, thickness and Mooney-Rivlin coefficient of rubber fender with experimental data. The equation of motion of the wind met mast is derived under the assumption that it was ignored vertical movement effect of barge on sea water. Such an analyzing method which was developed so far, make it possible to determine the proper size and material properties of rubber fender and the optimal moving conditions of barge, and finally, application method can be suggested in designing process of rubber fender considering barge impact.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.85-94
• /
• 2011

Mixtures of sand and rubber particles ($D_{sand}/D_{rubber}=1$) are investigated to explore their characteristics under different stain level. Mixtures are prepared with different volumetric sand fractions ($sf=V_{sand}/V_{total}$). Experimental data are gathered from a resonant column, an instrumented oedometer, and a direct shear tests. Results show that sand and rubber differently control the behavior of the whole mixture with strain level. Non-linear degradation of small strain stiffness is observed for the mixtures with $sf{\geq}0.4$, while the mixtures with low sand fraction ($sf{\leq}0.2$) show significantly high elastic threshold strain. Vertical stress-deformation increases dramatically when the rubber particle works as a member of force chain. The strength of the mixtures increases as the content of rubber particle decreases, and contractive behavior is observed in the mixtures with $sf{\leq}0.8$. Rubber particle plays different roles with strain level in the mixture: it increases a coordination number and controls a plasticity of the mixture in small strain; it prevents a buckling of force chain in intermediate strain; it leads a contractive behavior in large strain.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.523-530
• /
• 2017

In this study, a method of probabilistic evaluation of the performance point of the structure obtained by capacity spectrum method (CSM) is presented. The performance point of the 4-story and 1-bay steel structure was determined by using CSM according to ATC-40. In order to analyze whether the demand spectrums exceed the performance limit of the structure, the limit displacements are derived for the performance limit of the structure defined from the plastic deformation angle of the structural member. In addition, by selecting a total of 30 artificial seismic wave having the response spectrum similar to the design response spectrum, the fragility curves were derived by examining whether the response spectrum obtained from the artificial seismic wave were exceeded each performance limit according to the spectral acceleration. The maximum likelihood method was used to derive the fragility curve using observed excess probabilities. It has been confirmed that there exists a probability that the response acceleration value of the design response spectrum corresponding to each performance limit exceeds the performance limit. This method has a merit that the stochastic evaluation can be performed considering the uncertainty of the seismic waves with respect to the performance point of the structure, and the analysis time can be shortened because the incremental dynamic analysis (IDA) is not necessary.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.6
• /
• pp.4267-4275
• /
• 2015

Prefabricated houses are fabricated at the factory for approximately 60 to 80% of the entire construction process, and assembled in the field. In the process of transporting and lifting, internal and external finishes of the unit module are concerned about damages. The purpose of this study is to improve the fixing equipment by analyzing load behavior. The improved fixing equipment would minimize the deformation of internal and external finishes. In order to develop the improved fixing equipment, transporting load on the fixing equipment is analyzed using Monte Carlo simulations, and structural performance is verified by the non-linear finite element analysis. Statistical analysis shows load distribution of unit module is similar with extreme value distribution. Based on the statistical analysis and Monte Carlo simulation, the maximum transporting load is 28.9kN and 95% confidence interval of transporting load is -1.22kN to 9.5kN. The nonlinear structural analysis shows improved fixing equipment is not destructed to the limit load of 35.3kN and withstands the load-bearing in the 95% confidence interval of transporting load.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

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

Considering the non-linear behavior of structure and soil when evaluating a nuclear power plant's seismic safety under a beyond-design basis earthquake is essential. In order to obtain the nonlinear response of a nuclear power plant structure, a time-domain SSI analysis method that considers the nonlinearity of soil and structure and the nonlinear Soil-Structure Interaction (SSI) effect is necessary. The Boundary Reaction Method (BRM) is a time-domain SSI analysis method. The BRM can be applied effectively with a Perfectly Matched Layer (PML), which is an effective energy absorbing boundary condition. The BRM has a characteristic that the magnitude of the response in far-field soil increases as the boundary interface of the effective seismic load moves outward. In addition, the PML has poor absorption performance of low-frequency waves. For this reason, the accuracy of the low-frequency response may be degraded when analyzing the combination of the BRM and the PML. In this study, the accuracy of the analysis response was improved by adjusting the PML input parameters to improve this problem. The accuracy of the response was evaluated by using the analysis response using KIESSI-3D, a frequency domain SSI analysis program, as a reference solution. As a result of the analysis applying the optimal PML parameter, the average error rate of the acceleration response spectrum for 9 degrees of freedom of the structure was 3.40%, which was highly similar to the reference result. In addition, time-domain nonlinear SSI analysis was performed with the soil's nonlinearity to show this study's applicability. As a result of nonlinear SSI analysis, plastic deformation was concentrated in the soil around the foundation. The analysis results found that the analysis method combining BRM and PML can be effectively applied to the seismic response analysis of nuclear power plant structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.4A
• /
• pp.281-290
• /
• 2009

In this study, an analysis technique of hybrid girder considering nonlinearity of steel-concrete contact surface is presented. Steel-concrete hybrid girder shows partial-interaction behavior due to the deformation of shear connectors, slip and detachment at the interface, and cracks under the applied loads. Therefore, the partial-interaction approach becomes more reasonable. Contact surface is modeled by interface element and analyzed nonlinearly because of cost of time and effort to detailed model and analysis. Steel and Concrete are modeled considering non-linearity of materials. Material property of contact surface is obtained from push-out test and input to interface element. For the constitutive models, Drucker-Prager and smeared cracking model are used for concrete in compression and tension, respectively, and a von-Mises model is used for steel. This analysis technique is verified by comparing it with test results. Using verified analysis technique, various analyses are performed with different parameters such as nonlinear material property of interface element and prestress. The results are compared with linear analysis result and analysis result with the assumption of full-interaction.