• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1679-1688
• /
• 1990

The method of nonlinear parameter estimation(NPE), which is a statistical and an inverse method, is used to estimate the thermal diffusivity of the porous insulation material. In order to apply the NPE method for measuring the thermal diffusivity, and algorithm for programing suitable to IBM personal computer is established, and is studied the statistical treatment of experimental data and theory of estimation. The experimental data obtained by discrete measurement using a constant heat flux technique are used to find the boundary conditions, initial conditions, and the thermal diffusivity, and then the final values are compared with the values obtained by some different methods. The results are presented as follows:(1) NPE method is used to establish the estimation of the thermal diffusivity and compared results with experimental output shows, that this method can be applicable to define the thermal diffusivity without considering hear flux types. (2) Because of all of the temperatures obtained by the discrete measurement on each steps of time are used to estimate the thermal diffusivity. Although some error in the temperature measurements of temperature are included in estimating process, its influences on the final value are minimzed in NPE method. (3) NPE method can reduce the experimental time including the time of data collecting in a few minutes and can take smaller specimen compared with steady state method. If the tube-type furnace is used, also the adjusting time of surrounding temperature can be reduced.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Structural Engineering and Mechanics
• /
• v.25 no.4
• /
• pp.445-466
• /
• 2007

As the LNG (Liquefied Natural Gas) tank contains cryogenic liquid, realistic thermal analyses are of a primary importance for a successful design. The structural details of the LNG tank are so complicated that some strategies are necessary to reasonably predict its temperature distribution. The proposed heat transfer model can consider the beneficial effects of insulation layers and a suspended deck on temperature distribution of the outer concrete tank against cryogenic conditions simply by the boundary conditions of the outer tank model. To this aim, the equilibrium condition or heat balance in a steady state is utilized in a various way, and some aspects of heat transfer via conduction, convection and radiation are implemented as necessary. Overall thermal analysis procedures for the LNG tank are revisited to examine some unjustifiable assumptions of conventional analyses. Concrete and insulation properties under cryogenic condition and a reasonable conversion procedure of the temperature-induced nonlinear stress into the section forces are discussed. Numerical examples are presented to verify the proposed schemes in predicting the actual temperature and stress distributions of the tank as affected by the cryogenic LNG for the cases of normal operation and leakage from the inner steel tank. It is expected that the proposed schemes enable a designer to readily detect the effects of insulation layers and a suspended deck and, therefore, can be employed as a useful and consistent tool to evaluate the thermal effect in a design stage of an LNG tank as well as in a detailed analysis.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.147-147
• /
• 2019

The Numerical Weather Prediction (NWP) models determine the future state of the weather by forcing current weather conditions into the atmospheric models. The NWP models approximate mathematically the physical dynamics by nonlinear differential equations; however these approximations include uncertainties. The errors of the NWP estimations can be related to the initial and boundary conditions and model parameterization. Development in the meteorological forecast models did not solve the issues related to the inevitable biases. In spite of the efforts to incorporate all sources of uncertainty into the forecast, and regardless of the methodologies applied to generate the forecast ensembles, they are still subject to errors and systematic biases. The statistical post-processing increases the accuracy of the forecast data by decreasing the errors. Error prediction of the NWP models which is updating the NWP model outputs or model output statistics is one of the ways to improve the model forecast. The regression methods (including linear, polynomial and scaling regression) are applied to the present study to improve the real time forecast skill. Such post-processing consists of two main steps. Firstly, regression is built between forecast and measurement, available during a certain training period, and secondly, the regression is applied to new forecasts. In this study, the WRF real-time forecast data, in comparison with the observed data, had systematic biases; the errors related to the NWP model forecasts were reflected in the underestimation of the meteorological data forecast by the WRF model. The promising results will indicate that the post-processing techniques applied in this study improved the meteorological forecast data provided by WRF model. A comparison between various bias correction methods will show the strength and weakness of the each methods.

• Geomechanics and Engineering
• /
• v.25 no.3
• /
• pp.195-205
• /
• 2021

Time effect on the deformation and strength characteristics of geogrid reinforced sand retaining wall has become an important issue in geotechnical and transportation engineering. Three physical model tests on geogrid reinforced sand retaining walls performed under various loading conditions were simulated to study their rate-dependent behaviors, using the presented nonlinear finite element method (FEM) analysis procedure. This FEM was based on the dynamic relaxation method and return mapping scheme, in which the combined effects of the rate-dependent behaviors of both the backfill soil and the geosynthetic reinforcement have been included. The rate-dependent behaviors of sands and geogrids should be attributed to the viscous property of materials, which can be described by the unified three-component elasto-viscoplastic constitutive model. By comparing the FEM simulations and the test results, it can be found that the present FEM was able to be successfully extended to the boundary value problems of geosynthetic reinforced soil retaining walls. The deformation and strength characteristics of the geogrid reinforced sand retaining walls can be well reproduced. Loading rate effect, the trends of jump in footing pressure upon the step-changes in the loading rate, occurred not only on sands and geogrids but also on geogrid reinforced sands retaining walls. The lateral earth pressure distributions against the back of retaining wall, the local tensile force in the geogrid arranged in the retaining wall and the local stresses beneath the footing under various loading conditions can also be predicted well in the FEM simulations.

• Proceedings of the KWS Conference
• /
• 2009.11a
• /
• pp.97-97
• /
• 2009

The use of thin plate increases due to the need for light weight in large ship. Thin plate is easily distorted and has residual stress by welding heat. Therefore, the thin plate should be carefully joined to minimize the welding deformation which costs time and money for repair. For one effort to reduce welding deformation, it is very useful to predict welding deformation before welding execution. There are two methods to analyze welding deformation. One is simple linear analysis. The other is nonlinear analysis. The simple linear analysis is elastic analysis using the equivalent load method or inherent strain method from welding experiments. The nonlinear analysis is thermo-elastic analysis which gives consideration to the nonlinearity of material dependent on temperature and time, welding current, voltage, speed, sequence and constraint. In this study, the welding deformation is analyzed by using thermo-elastic method for PCTC(Pure Car and Truck Carrier) which carries cars and trucks. PCTC uses thin plates of 6mm thickness which is susceptible to welding heat. The analysis dimension is 19,200mm(length) * 13,825mm(width) * 376mm(height). MARC and MENTAT are used as pre and post processor and solver. The boundary conditions are based on the real situation in shipyard. The simulations contain convection and gravity. The material of the thin block is mild steel with $235N/mm^2$ yield strength. Its nonlinearity of conductivity, specific heat, Young's modulus and yield strength is applied in simulations. Welding is done in two pass. First pass lasts 2,100 second, then it rests for 900 second, then second pass lasts 2,100 second and then it rests for 20,000 second. The displacement at 0 sec is caused by its own weight. It is maximum 19mm at the free side. The welding line expands, shrinks during welding and finally experiences shrinkage. It results in angular distortion of thin block. Final maximum displacement, 17mm occurs around welding line. The maximum residual stress happens at the welding line, where the stress is above the yield strength. Also, the maximum equivalent plastic strain occurs at the welding line. The plastic strain of first pass is more than that of second pass. The flatness of plate in longitudinal direction is calculated in parallel with the direction of girder and compared with deformation standard of ${\pm}15mm$. Calculated value is within the standard range. The flatness of plate in transverse direction is calculated in perpendicular to the direction of girder and compared with deformation standard of ${\pm}6mm$. It satisfies the standard. Buckle of plate is calculated between each longitudinal and compared with the deformation standard. All buckle value is within the standard range of ${\pm}6mm$.

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

Fixed connection is generally used for beam and column connections of concrete structures, but significant damages at the connection due to severe earthquakes have been reported. In order to reduce damages of the connection and improve seismic performance of the connection, several innovative connections have been suggested. One newly proposed connection type allows a rotation of the connection for applications in rotating or rocking beams, columns, and shear walls. Such structural elements would provide a nonlinear lateral force-displacement response since their contact depth developed during rotation is gradually reduced and the stress across the sections of the elements is non-linearly distributed around a contact area, which is called an elastic hinge region in the present study. The purpose of the present study is to define the elastic hinge region or length for the rocking columns, through investigating the cross-sectional stress distribution during their lateral behavior. Performing a finite element analysis (FEA), several parameters are considered including axial load levels (5% and 10% of nominal strength), different boundary conditions (confined-ends and cantilever types), and slenderness ratios (length/depth = 5, 7, 10). The FEA results showed that the elastic hinge length does not directly depend on the parameters considered, but it is governed by a contact depth only. The elastic hinge length started to develop after an opening state and increased non-linearly until a rocking point(pre-rocking). However, the length did not increase any more after the rocking point (post-rocking) and remained as a constant value. Half space model predicting the elastic hinge length is adapted and the results are compared with the numerical results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• The Journal of the Convergence on Culture Technology
• /
• v.7 no.2
• /
• pp.405-410
• /
• 2021

In this study, a finite element analysis was performed applying a nonlinear material model and fatigue load conditions to evaluate the service life and spring stiffness of the resilient pad for rail fastening system. As a result of the fatigue analysis, the rate of change in spring stiffness compared to the initial condition was about 16%, indicating that fatigue hardening occurred. As for the stress generated in the longitudinal direction of the resilient pad, the difference between the stress generated at the center and the edge was about 10 times or more. In addition, it was analyzed that the equivalent stress of the outer boundary was more than twice as large as that of the central part. Therefore, it was analyzed that the damage and deformation of the resilient pad are the corners of the resilient pad under actual service conditions. The fatigue life diagram of the resilient pad (S-N curve) was derived using the equivalent stress of the resilient pad according to the fatigue cycles. Using the fatigue life diagram of the resilient pad derived in this study, it is considered that it can be used to predict the fatigue life under the relevant conditions by calculating the equivalent stress of the resilient pad under various load conditions.

• Journal of the Korea Institute of Building Construction
• /
• v.13 no.5
• /
• pp.491-497
• /
• 2013

The objective of the present study is to establish a straightforward mixture proportioning procedure for structural lightweight aggregate concrete (LWAC), and evaluate the selection range of the targeted dry density of concrete against the designed concrete compressive strength. In developing this procedure, mathematical models were formulated based on a nonlinear regression analysis over 347 data sets and two boundary conditions of the absolute volume and dry density of concrete. The proposed procedure demonstrated the appropriate water-to-cement ratio and dry density of concrete to achieve the designed strength decrease with the increase in volumetric ratio of coarse aggregates. This trend was more significant in all-LWAC than in sand-LWAC. Overall, the selection range of the dry density of LWAC exists within a certain range according to the designed strength, which can be obtained using the proposed procedure.