• Geomechanics and Engineering
• /
• v.27 no.5
• /
• pp.465-479
• /
• 2021

One of the important causes of building and infrastructure failure, such as bridges on pile foundations, is the placement of the piles in liquefiable soil that can become unstable under seismic loads. Therefore, the overarching aim of this study is to investigate the seismic behavior of a soil-pile system in liquefiable soil using three-dimensional numerical FEM analysis, including soil-pile interaction. Effective parameters on concrete pile response, involving the pile diameter, pile length, soil type, and base acceleration, were considered in the framework of finite element non-linear dynamic analysis. The constitutive model of soil was considered as elasto-plastic kinematic-isotropic hardening. First, the finite element model was verified by comparing the variations on the pile response with the measured data from the centrifuge tests, and there was a strong agreement between the numerical and experimental results. Totally 64 non-linear time-history analyses were conducted, and the responses were investigated in terms of the lateral displacement of the pile, the effect of the base acceleration in the pile behavior, the bending moment distribution in the pile body, and the pore pressure. The numerical analysis results demonstrated that the relationship between the pile lateral displacement and the maximum base acceleration is non-linear. Furthermore, increasing the pile diameter results in an increase in the passive pressure of the soil. Also, piles with small and big diameters are subjected to yielding under bending and shear states, respectively. It is concluded that an effective stress-based ground response analysis should be conducted when there is a liquefaction condition in order to determine the maximum bending moment and shear force generated within the pile.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.77-86
• /
• 1999

In the previous $paper^{(1),(2)}$, a geometrically non-linear finite element formulation of spatial cables subjected to self-weights and support motions was presented using multiple noded cable elements and how to determine the initial equililbrium state of cables was addressed. In this paper, in order to perform dynamic non-linear analysis of ocean cables subjected to support motions and earthquakes, a numerical method to calculate Morison forces and incorporate effects of earthquake motions is presented based on the Newmark method. Challenging example problems are presented in order to investigate dynamic non-linear behaviors of ocean cables subjected to support motions and earthquake loadings.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2002.04a
• /
• pp.122-126
• /
• 2002

Contaminant transport in porous media is characterized by solving an advection-dispersion equation(ADE). The ADE can cover equilibrium phenomena of interest, which include sorption, decay, and chemical reactions. Among these phenomena, sorption mechanism is described by several types of sorption isotherm. If we assume the sorption isotherm as linear, the solution of ADE can be easily procured. However, if we consider the sorption isotherm as non-linear isotherm like a Dual Reactive Domain Model (DRDM), the resulting differential equation becomes non-linear. In this case, the solution of ADE cannot be easily acquired by an analytic method. In this paper, we present the numerical analysis of ADE using a DRDM. The results reveal that even if sorption data may be fitted well using linear or non-linear isotherm, the characteristics of contaminant transport of the two cases are different from each other. To be concrete, the retardation of linear isotherm has stronger effect than that of the DRDM. As the non-linearity of sorption isotherm increases, the difference of retardation effects of the two cases becomes larger. For a pulse source, the maximum concentration of the linear model is higher than that of the DRDM, but the plume of the DRDM moves faster than that of the linear model. Behaviors of contaminant transport using the DRDM are consistent with common features of a linear model. For instance, biodegradation effect becomes larger as time goes by The faster the seepage velocity is, the faster the plume of contaminant moves. The plume of the contaminant is distributed evenly over overall domain in the event of high dispersion coefficient.

• Structural Engineering and Mechanics
• /
• v.12 no.1
• /
• pp.35-50
• /
• 2001

A numerical methodology is presented in this paper for the geometrically non-linear analysis of slender uni-dimensional structural elements under unilateral contact constraints. The finite element method together with an updated Lagrangian formulation is used to study the structural system. The unilateral constraints are imposed by tensionless supports or foundations. At each load step, in order to obtain the contact regions, the equilibrium equations are linearized and the contact problem is treated directly as a minimisation problem with inequality constraints, resulting in a linear complementarity problem (LCP). After the resulting LCP is solved by Lemke's pivoting algorithm, the contact regions are identified and the Newton-Raphson method is used together with path following methods to obtain the new contact forces and equilibrium configurations. The proposed methodology is illustrated by two examples and the results are compared with numerical and experimental results found in literature.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.107-114
• /
• 2002

This paper deals with the non-linear analysis of cantilever beams with constant volume. Numerical methods are developed for solving the elastica of cantilever ben subjected to a tip Point load and a tip couple. The linear, parabolic and sinusoidal tapers with the regular polygon cross-section are considered, whose material volume and span length are always held constant. The Runge-Kutta and Regula-Falsi methods, respectively, are used to integrate the governing differential equations and to compute the unknown value of the tip deflection. The numerical results obtained herein are shown in tables and figures. Also the shapes of strongest beams are determined by reading the minimum values form the deflection versus section ratio curves.

• Coupled systems mechanics
• /
• v.10 no.1
• /
• pp.1-19
• /
• 2021

Composite nature of the masonry structures in general causes complex and non-linear behaviour, especially in intense vibration conditions. The presence of different types and forms of structural elements and different materials is a major problem for the analysis of these type of structures. For this reason, the analysis of the behaviour of masonry structures requires a combination of experimental tests and non-linear mathematical modelling. The famous UNESCO Heritage Old Bridge in Mostar was selected as an example for the analysis of the global behaviour of reinforced stone arch masonry bridges. As part of the experimental research, a model of the Old Bridge was constructed in a scale of 1:9 and tested on a shaking table platform for different levels of seismic excitation. Non-linear mathematical modelling was performed using a combined finite-discrete element method (FDEM), including the effect of connection elements. The paper presents the horizontal displacement of the top of the arch and the failure mechanism of the Old Bridge model for the experimental and the numerical phase, as well as the comparison of the results. This research provided a clearer insight into the global behaviour of stone arch masonry structures reinforced with steel clamps and steel dowels, which is significant for the structures classified as world cultural heritage.

• Composites Research
• /
• v.18 no.2
• /
• pp.30-37
• /
• 2005

In this paper, geometrically non-linear finite element analyses were performed to study the mechanical behavior of the material system of the envelope of stratospheric airships. The microstructure of the load-bearing plain weave layer was identified and modeled. The Updated Lagrangian formulation was employed to consider the geometric non-linearity as well as the induced structural non-linearity for the fiber tows. The stress-strain behavior was predicted and the effective elastic modulus was calculated by numerical experiments. It was found the non-linear stress-strain curves were largely different from those by linear analysis. And non-linear elastic moduli were much higher than linear elastic moduli. The difference was more distinguishable when the tow waviness ratio was smaller.

• Journal of Advanced Marine Engineering and Technology
• /
• v.15 no.4
• /
• pp.46-53
• /
• 1991

In the present study, Navier-Stokes equation is numerically solved by use of a Finite analytic method to obtain the 2-dimensional flow field in the square cavity. The basic idea of F.A.M. is the incorporation of local analytic solutions in the numerical solution of linear or non-linear partial differential equations. In the F.A.M., the total problem is subdivided into a number of all elements. The local analytic solution is obtained for the small element in which the governing equation, if non-linear, to be linearized. The local analytic solutions are then expressed in algebraic form and are overlapped to cover the entire region of the problem. The assembly of these local analytic solutions, which still preserve the overall nonlinearity of the governing equations, results in a system of linear algebraic equations. The system of algebraic equations is then solved to provide the numerical solutions of the total problem. The computed flow field shows the same characteristics to physical concept of flow phenomena.

• Structural Engineering and Mechanics
• /
• v.85 no.4
• /
• pp.485-500
• /
• 2023

The non-linear static analysis of reinforced concrete (RC) structures using the three-dimensional (3D) finite element method is a time-consuming and challenging task. Moreover, this type of analysis encounters numerical problems such as the lack of convergence of results in the stages of growth and propagation of cracks in the structure. The time integration analysis along with the mass scaling (MS) technique is usually used to overcome these limitations. Despite the use of this method in the 3D finite element analysis of RC structures, a comprehensive study has not been conducted so far to assess the effects of the MS method on the accuracy of results. This study aims to evaluate the accuracy of the MS method in the non-linear quasi-static finite element analysis of RC structures. To this aim, different types of RC structures were simulated using the finite element approach based on the implicit time integration method and the mass scaling technique. The influences of effective parameters of the MS method (i.e., the allowable values of increase in the mass of the RC structure, the relationship between the duration of the applied load and fundamental vibration period of the RC structure, and the pattern of applied loads) on the accuracy of the simulated results were investigated. The accuracy of numerical simulation results has been evaluated through comparison with existing experimental data. The results of this study show that the achievement of accurate structural responses in the implicit time integration analyses using the MS method involves the appropriate selection of the effective parameters of the MS method.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.5
• /
• pp.1-11
• /
• 2001

The purpose of this study is to develop the main program and pre/post processor for the geometric and plastic non-linear analysis of steel jacket structures subjected to wave and earthquake loadings. In this paper, steel jackets are modelled using geometric non-linear space frames and wave loadings re evaluated based on Morrison equation using the linear Airy theory and the fifth Stokes theory. Random wave is generated using JONSWAP spectrum. For earthquake analysis, dynamic analysis is performed using artificial earthquake time history. Also the plastic hinge method is presented for limit analysis of steel jacket. In the companion paper, the pre/post processor is developed and the numerical examples are presented for linear and non-linear dynamic analysis of steel jackets.