• Steel and Composite Structures
• /
• 제38권5호
• /
• pp.497-521
• /
• 2021

In this study, an effective numerical method is introduced for nonlinear inelastic analyses of rectangular concrete-filled steel tubular (CFST) frames for the first time. A steel-concrete composite fiber beam-column element model is developed that considers material, and geometric nonlinearities, and residual stresses. This is achieved by using stability functions combined with integration points along the element length to capture the spread of plasticity over the composite cross-section along the element length. Additionally, a multi-spring element with a zero-length is employed to model the nonlinear semi-rigid beam-to-column connections in CFST frame models. To solve the nonlinear equilibrium equations, the generalized displacement control algorithm is adopted. The accuracy of the proposed method is firstly verified by a large number of experiments of CFST members subjected to various loading conditions. Subsequently, the proposed method is applied to investigate the nonlinear inelastic behavior of rectangular CFST frames with fully rigid, semi-rigid, and hinged connections. The accuracy of the predicted results and the efficiency pertaining to the computation time of the proposed method are demonstrated in comparison with the ABAQUS software. The proposed numerical method may be efficiently utilized in practical designs for advanced analysis of the rectangular CFST structures.

• Korean Journal of Mathematics
• /
• 제22권2호
• /
• pp.355-365
• /
• 2014

We show the existence of nontrivial solutions for some fourth order elliptic boundary value problem with fully nonlinear term. We obtain this result by approaching the variational method and using a linking theorem. We also get a uniqueness result.

• Ocean Systems Engineering
• /
• 제9권3호
• /
• pp.241-259
• /
• 2019

This paper presents a fully three-dimensional numerical approach for analyzing deepwater drilling riser-conductor system vortex-induced vibrations (VIV) including nonlinear soil-structure interactions (SSI). The drilling riser-conductor system is modeled as a tensioned beam with linearly distributed tension and is solved by a fully implicit discretization scheme. The fluid field around the riser-conductor system is obtained by Finite-Analytic Navier-Stokes (FANS) code, which numerically solves the unsteady Navier-Stokes equations. The SSI is considered by modeling the lateral soil resistance force according to nonlinear p-y curves. Overset grid method is adopted to mesh the fluid domain. A partitioned fluid-structure interaction (FSI) method is achieved by communication between the fluid solver and riser motion solver. A riser-conductor system VIV simulation without SSI is firstly presented and served as a benchmark case for the subsequent simulations. Two SSI models based on a nonlinear p-y curve are then applied to the VIV simulations. Also, the effects of two key soil properties on the VIV simulations of riser-conductor systems are studied.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제5권4호
• /
• pp.513-528
• /
• 2013

A floating Oscillating Water Column (OWC) wave energy converter, a Backward Bent Duct Buoy (BBDB), was simulated using a state-of-the-art, two-dimensional, fully-nonlinear Numerical Wave Tank (NWT) technique. The hydrodynamic performance of the floating OWC device was evaluated in the time domain. The acceleration potential method, with a full-updated kernel matrix calculation associated with a mode decomposition scheme, was implemented to obtain accurate estimates of the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT was based on the potential theory and the boundary element method with constant panels on the boundaries. The mixed Eulerian-Lagrangian (MEL) approach was employed to capture the nonlinear free surfaces inside the chamber that interacted with a pneumatic pressure, induced by the time-varying airflow velocity at the air duct. A special viscous damping was applied to the chamber free surface to represent the viscous energy loss due to the BBDB's shape and motions. The viscous damping coefficient was properly selected using a comparison of the experimental data. The calculated surface elevation, inside and outside the chamber, with a tuned viscous damping correlated reasonably well with the experimental data for various incident wave conditions. The conservation of the total wave energy in the computational domain was confirmed over the entire range of wave frequencies.

• 한국해양공학회지
• /
• 제21권5호
• /
• pp.1-8
• /
• 2007

A long wave induced by a Gaussian-shape submarine landslide is simulated by a 2D fully nonlinear numerical wave tank (NWT). The NWT is based on the boundary element method and the mixed Eulerian/Lagrangian approach. Using the NWT, physical characteristics of land-slide tsunami, including wave generation, propagation, particle kinematics, hydrodynamic pressure, run-up and depression, are simulated for the early stage of long wave generation and propagation. Various sliding mass heights are applied to the developed model for a systematic sensitivity analysis. In particular, the fully nonlinear NWT results are compared with linear results (exact body-boundary conditions with linear free-surface conditions) to identify the nonlinear effects in the respective cases.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2000년도 춘계학술대회 논문집
• /
• pp.140-145
• /
• 2000

The spectral method which is composed of an eigenfunction expansion of free modes in the wave number domain is used to produce two dimensional unsteady inviscid wave simulation such as progressive waves in a numerical pneumatic wave tank. A spatial and time dependent free surface elevation and the potential are calculated by integrating ODE derived from fully nonlinear kinematic and dynamic free surface boundary condition at each time step. The nonlinear characteristics in the waves by this method were notable as increasing wave steepness. This method is very useful and powerful in terms of saving computational time caused by rapid convergence exponentially with increasing number of nodes, even preserving accurate numerical results. Moreover, it will given us many possibilities to apply to naval and ocean engineering fields.

• 한국해안해양공학회지
• /
• 제13권4호
• /
• pp.296-308
• /
• 2001

Wei et al.의 완전비선형 Boussinesq방정식을 4차의 Adams predictorcorrector기법을 사용하여 차분하고 면 내부조파기법과 스폰지 경계충을 사용하였으며 쇄파구조를 추가하였다. 면 내부조파기법을 사용해 목적파를 잘 재현할 수 있었다. 비선형성이 부각되는 고립파의 천수실험을 통해 완전비선형 모형의 정화성을 확인할 수 있었고 해저평붕으로 인한 규칙파의 변형을 모의해 보았다. 쇄파 수치실험을 통해 적용된 쇄파구조의 특성을 확인해 보았고 수중천퇴에 대한 이차원 파랑전파 수치실험을 통해 비선형 모형의 중요성을 확인하였다.

• 대한수학회보
• /
• 제57권1호
• /
• pp.219-244
• /
• 2020

In this paper, we approximate the solution of the coupled nonlinear Schrödinger equations by using a fully discrete finite element scheme based on the standard Galerkin method in space and implicit midpoint discretization in time. The proposed scheme guarantees the conservation of the total mass and the energy. First, a priori error estimates for the fully discrete Galerkin method is derived. Second, the existence of the approximated solution is proved by virtue of the Brouwer fixed point theorem. Moreover, the uniqueness of the solution is shown. Finally, convergence orders of the fully discrete Crank-Nicolson scheme are discussed. The end of the paper is devoted to some numerical experiments.

• 한국해양공학회지
• /
• 제23권6호
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• Nuclear Engineering and Technology
• /
• 제41권7호
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.