• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.95-102
• /
• 2014

In this paper, a coupling scheme for applying finite element analysis(FEA) programs, such as, LS-DYNA and MIDAS/Civil, to a nonlinear soil structure interaction analysis by the boundary reaction method(BRM) is presented. With the FEA programs, the structure and soil media are discretized by linear or nonlinear finite elements. To absorb the outgoing elastic waves to unbounded soil region as much as possible, the PML elements and viscous-spring elements are used at the outer FE boundary, in the LS-DYNA model and in MIDAS/Civil model, respectively. It is also assumed that all the nonlinear elements in the problem are limited to structural region. In this study, the boundary reaction forces for the use in the BRM are calculated using the KIESSI-3D program by solving soil-foundation interaction problem subjected to incident seismic waves. The effectiveness of the proposed approach is demonstrated with a linear SSI seismic analysis problem by comparing the BRM solution with the conventional SSI solution. Numerical comparison indicates that the BRM can effectively be applied to a nonlinear soil-structure analysis if motions at the foundation obtained by the BRM for a linear SSI problem excluding the nonlinear structure is conservative.

• Geomechanics and Engineering
• /
• v.8 no.4
• /
• pp.559-594
• /
• 2015

A nonlinear creep function incorporated into the elastic visco-plastic model may describe the long-term soil deformation more accurately. However, by applying the conventional procedure, there are challenges to determine the model parameters due to limitation of suitable data points. This paper presents a numerical solution to obtain several parameters simultaneously for a nonlinear elastic visco-plastic (EVP) model using the available consolidation data. The finite difference scheme using the Crank-Nicolson procedure is applied to solve a set of coupled partial differential equations of the time dependent strain and pore water pressure dissipation. The model parameters are determined by applying the algorithm of trust-region reflective optimisation in conjunction with the finite difference solution. The proposed method utilises all available consolidation data during dissipation of the excess pore water pressure to determine the required model parameters. Moreover, the reference time in the elastic visco-plastic model can readily be adopted as a unit of time; denoting creep is included in the numerical predictions explicitly from the very first time steps. In this paper, the settlement predictions of thick soft clay layers are presented and discussed to evaluate and compare the accuracy and reliability of the proposed method against the graphical procedure to obtain the model parameters. In addition, comparison of the available experimental results to the numerical predictions confirms the accuracy of the numerical procedure.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.4
• /
• pp.45-54
• /
• 2007

In the context of auto-landing containers from a container ship to a truck or automatic guided vehicle and vice versa, this research investigates three schemes, one in Part I and two in Part II, for measuring the absolute position of a container. Coordinate transformations between the reference-coordinate, sensor-coordinate, and body-coordinate systems are briefly discussed. The scheme explored in Part I aims the use of three laser-slit sensors, which are relatively inexpensive. In this case, nine nonlinear equations are formulated for six unknown variables (three for orientation and three for position), so a closed-form solution is not available. Instead, an approximate solution through linearization was derived. An advantage of the method in Part I is its ability to measure an absolute position in 3D space, while a disadvantage is the computation time required to obtain pseudo-inverses and the approximate nature of the obtained solution. Numerical examples are provided.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1975-1981
• /
• 1991

A new numerical solution of the elastohydrodynamic lubrication(EHL) problem of an axially profiled cylindrical roller is presented. A finite difference method and the Newton-Raphson method are used to solve the nonlinear system equations. A non-uniform grid system is adopted to reduce the number of grid points and to obtain accurate solution. For two different types of profiles which have similar elastostatic pressure distribution, the EHL results show large differences. Especially the difference in film shape is larger than in pressure distribution. Therefore, the magnitude of the minimum film thickness should be a major criteria to design the axial profile of the roller. Variations of the minimum film thickness with dimensionless parameters show considerably different behavior from those of infinite solution and show a good agreement with the experimental data in literatures. Present numerical scheme can be used generally in the analysis of three-dimensional EHL problem.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4
• /
• pp.23-31
• /
• 1992

In this study, Behavior of Eccentrically Loaded Restrained Concrete Slender Columns are investigated the continuation. The paper is : nonlinear constitutive relations of reinforced concrete columns are formulated, the multivariable nonlinear equation formulated in this study is solved by the incremental iterative method, the solution scheme is suitable for computer application, the proposed method bas been compared with the numerical examples and other research results and is shown to be in good agreement Accordingly, the objectives of this investigation are to use the method developed to conduct a parametric study of the different variable (Rotational Restraints,Ratio of End Restraints. End Eccentricity, Ratio of End Eccentricities,Lateral Restraint, Concrete Strength, Ratio of Reinforcement. Yield Point of Reinforcement) influencing the behavior and the ultimate strength for columns. Chosen values for different parameters and their variations are based on values encountered in practice and on results of other analylitical studies.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.659-669
• /
• 2018

The Dynamic equations in the polar coordinates are drawn out using the MLPG method for the non-symmetric FG cylindrical shell. To simulate the mechanical properties of FGM, the nonlinear volume fractions for radial direction are used. The shape function applied in this paper is a form of the radial basis functions, by using this function all the requirements for an effective and suitable shape function are established. Hence in this study, the multiquadrics (MQ) radial basis functions are exploited as the shape function governing the problem. The MLPG method is combined with the the Newmark time approximation scheme to solve dynamic equations in the time domain. The obtained results by the MLPG method to be verified are compared with the analytical solution and the FEM. The obtained results through the MLPG method show a good agreement in comparison to other results and the MLPG method has high accuracy for dynamic analysis of the non-symmetric FG cylindrical shell. To demonstrate the capability of the present method to dynamic analysis of the non-symmetric FG cylindrical shell, it is analyzed dynamically with different volume fraction exponents under harmonic and rectangular shock loading. The present method shows high accuracy, efficiency and capability to dynamic analysis of the non-symmetric FG cylindrical shell with nonlinear grading patterns.

• Journal of Advanced Navigation Technology
• /
• v.10 no.1
• /
• pp.20-25
• /
• 2006

In these days, OFDM(Orthogonal Frequency Division Multiplexing) is adopted to support high-speed data communication based on multi-path RF channel, but it has some weak point. One of those is that it has a higher PAR(peak-to-average power ratio) compared with single-carrier method. If some PAR of the transmitted signal is high, nonlinear amplitude distortion has occurred when it pass through the HPA(high power amplifier). There is a solution to prevent nonlinear distortion using higher peak power HPA, but it makes inefficiency and a cost problem. In this paper, we choose the SLM(Selected Mapping) scheme, which transmit the lowest PAR signal after OFDM symbol mapping, in various schemes reducing PAR for OFDM system. And we derived the performances of SLM method in fading channel through computer simulations.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].