• Geomechanics and Engineering
• /
• v.29 no.2
• /
• pp.155-170
• /
• 2022

An accurate analysis of structures supported on soft soils and subjected to seismic loading requires the consideration of the soil-foundation-structure interaction. An important aspect of this interaction lies with the energy dissipation due to soil material damping. Unlike advanced constitutive models that can induce energy loss, the use of simple elastoplastic constitutive models requires additional damping. The frequency dependent Rayleigh damping is a formulation that is frequently used in dynamic analysis. The main concern of this formulation is the correct selection of the target damping ratio and the frequency range where the response is frequency independent. The objective of this study is to investigate the effects of the Rayleigh damping parameters in soil-pile-structure and soil-inclusion-platform-structure systems in the presence of soft soil under seismic loading. Three-dimensional analyses of both systems are carried out using the finite difference software Flac3D. Different values of target damping ratios and minimum frequencies are utilized. Several earthquakes are used to study the influence of different excitation frequencies in the systems. The soil response in terms of accelerations, displacements and strains is obtained. For the rigid elements, the results are presented in terms of bending moments and normal forces. The results show that when the frequency of the input motion is close to the minimum (central) frequency in the Rayleigh damping formulation, the overdamping amount is reduced, and the surface spectral acceleration of the analyzed pile and inclusion systems increases. Thus, the bending moments and normal forces throughout the piles and inclusions also increase.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.13 no.3
• /
• pp.237-245
• /
• 2008

We developed a computational method on coupled dynamics of tracked vehicle on seafloor and long flexible pipe. The tracked vehicle is modeled as rigid-body vehicle, and the linked flexible pipe is discretized according to a lumped-parameter model. The equations of motion of the rigid-body vehicle on the soft seafloor are combined with the governing equations of flexible pipe dynamics. Four Euler parameters method is used to express the orientations of the vehicle and the flexible pipe. In order to solve the nonlinear coupled dynamics of vehicle and flexible pipe an incremental-iterative formulation is implemented. For the time-domain integration $Newmark-\beta$ method is adopted. The total Jacobean matrix has been derived based on the incremental-iterative formulation. The interactions between the dynamics of flexible pipe and the mobility of the tracked vehicle on soft seafloor are investigated through numerical simulations in time domain.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

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

A powerful and efficient method for finding approximate solutions to initial-boundary-value problems in the transient coupled thermoelasticity is formulated in time domain using the finite element technique with time-marching strategy. The final system equations can be derived by the Guritin's variational principle using the definition of convolution integral. But, the finite element formulation for the equations of motion is modified by differentiating in time. Numerical results to some test problems are compared with analytical and other sophisticated approximate solutions. Stable responces are observed in all the given examples irrespective of incremental time steps and mesh shapes. In addition, it is shown that good numerical results are obtained even in coarser mesh or larger time step comparing to other numerical methods.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.2
• /
• pp.86-97
• /
• 2005

Accurate estimation of the wave-induced pore water pressure in the seabed is key factor in studying the stability of the seabed in the vicinity of coastal structure. Most of the existing numerical models for wave structure seabed interaction have been linked through applying hybrid numerical technique which is analysis method separating the wave field and seabed regime. Therefore, it is necessary to develope a numerical model f3r simulating accurately wave$\cdot$structure$\cdot$ seabed interaction under wave loadings by the single domain approach for wave field and seabed regime together. In this study, direct numerical simulation is newly proposed. In this model, modeled fluid drag has been used to detect the hydraulic properties according to the varied geometrical shape inside the porous media by considering the turbulence resistance as well as laminar resistance. Contrary to hybrid numerical technique, direct numerical simulation avoids the explicit formulation of the boundary conditions at the fluid/porous media interface. A good agreement has been obtained by the comparison between existed experimental results by hydraulic model test and direct numerical simulation results far wave $\cdot$structure$\cdot$seabed interaction. Therefore, the newly proposed numerical model is a powerful tool for estimating the nonlinear dynamic responses among a structure, its seabed foundation and water waves.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.20 no.3
• /
• pp.191-199
• /
• 2000

In this paper, a numerical analysis using the finite element method (FEM) is presented which models the eddy current testing (ECT) of tubes with 3-dimensional defects. For the description of 3-dimensional eddy current problems, the governing equation is derived from the Maxwell's equations. The 3-dimensional FEM formulation with hexahedral elements is carried out using the Galerkin weighted residual method. The INCONEL 600 steam generator tube with inner and outer diameter defects is adopted for the numerical analysis, and the ECT signal, which is the trajectory of the probe impedance, is calculated. For the verification of the numerical analysis method, results of numerical calculations and experiments are compared and they show good agreements. Based on this verification, several defect signals are predicted and their characteristics are investigated with the variation in the defect depth and the circumferential angle of the defect.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.22-30
• /
• 2008

This paper deals with the nonlinear optimal control approach to helicopter maneuver problems using the indirect method. We apply a penalty function to the deviation from a prescribed trajectory to convert the system optimality to an unconstrained optimal control problem. The resultant two-point boundary value problem has been solved by using the multiple-shooting method. This paper focuses on the effect of the number of shooting nodes and initialization methods on the numerical solution in order to define the minimum number of shooting nodes required for numerical convergence and to provide a method increasing convergence radius of the indirect method. The results of this study can provide an approach to improve numerical stability and convergence of the indirect method when we solve the optimal control problems of an inherently unstable helicopter system.

• Computers and Concrete
• /
• v.19 no.5
• /
• pp.579-588
• /
• 2017

Creep and shrinkage are the main types of volume change with time in concrete. These changes cause deflection, cracking and stresses that affect durability, serviceability, long-term reliability and structural integrity of civil engineering infrastructure. Although laboratory test may be undertaken to determine the deformation properties of concrete, these are time-consuming, often expensive and generally not a practical option. Therefore, relatively simple empirically design code models are relied to predict the creep strain. This paper reviews the accuracy of creep and shrinkage predictions of reinforced concrete (RC) shear walls structures strengthened with carbon fibre reinforced polymer (CFRP) plates, which is characterized by a widthwise varying fibre volume fraction. This review is yielded by three commonly used international "code type" models. The assessed are the: CEB-FIP MC 90 model, ACI 209 model and Bazant & Baweja (B3) model. The time-dependent behavior was investigated to analyze their seismic behavior. In the numerical formulation, the adherents and the adhesives are all modelled as shear wall elements, using the mixed finite element method. Several tests were used to demonstrate the accuracy and effectiveness of the proposed method. Numerical results from the present analysis are presented to illustrate the significance of the time-dependency of the lateral displacements and eigenfrequencies modes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Structural Engineering and Mechanics
• /
• v.75 no.5
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.