• Earthquakes and Structures
• /
• v.15 no.2
• /
• pp.113-122
• /
• 2018

This study presents a new beam-column model comprising material nonlinearity and joint flexibility to predict the nonlinear response of reinforced concrete structures. The nonlinear behavior of connections has an outstanding role on the nonlinear response of reinforced concrete structures. In presented research, the joint flexibility is considered applying a rotational spring at each end of the member. To derive the moment-rotation behavior of beam-column connections, the relative rotations produced by the relative slip of flexural reinforcement in the joint and the flexural cracking of the beam end are taken into consideration. Furthermore, the considered spread plasticity model, unlike the previous models that have been developed based on the linear moment distribution subjected to lateral loads includes both lateral and gravity load effects, simultaneously. To confirm the accuracy of the proposed methodology, a simply-supported test beam and three reinforced concrete frames are considered. Pushover and nonlinear dynamic analysis of three numerical examples are performed. In these examples the nonlinear behavior of connections and the material nonlinearity using the proposed methodology and also linear flexibility model with different number of elements for each member and fiber based distributed plasticity model with different number of integration points are simulated. Comparing the results of the proposed methodology with those of the aforementioned models describes that suggested model that only uses one element for each member can appropriately estimate the nonlinear behavior of reinforced concrete structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.3
• /
• pp.279-293
• /
• 2004

The analysis of linear static and free vibration problems of isotropic and laminated composite plates and shells is performed by the improved 9-node shell element with the new strain displacement relationship. In that relationship, the effect of new additional terms between the bending strain and displacement has been investigated in the warping problem. Natural co ordinate based strains, stresses and constitutive equations are used. The assumed natural strain method is used to alleviate both membrane and shear locking behavior from the element. The Lanczos method is employed in the calculation of the eigenvalues of laminated composite structures and the Gauss integration rule is adopted to evaluate the mass matrix. The numerical examples are compared with the analytical solutions to validate the current formulation and the results presented could be useful for the understanding of the behaviour of laminates under free vibration conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.2
• /
• pp.79-87
• /
• 2013

In this paper, an Euler equation solver based on a Cartesian-grid method and non-uniform staggered grid system is applied to predict the ship motion response and added resistance in waves. Water, air, and solid domains are identified by a volume-fraction function for each phase and in each cell. For capturing the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with a weighed line interface calculation (WLIC) method. The volume fraction of solid body embedded in a Cartesian-grid system is calculated by a level-set based algorithm, and the body boundary condition is imposed by volume weighted formula. Added resistance is calculated by direct pressure integration on the ship surface. Numerical simulations for a Wigley III hull and an S175 containership in regular waves have been carried out to validate the newly developed code, and the ship motion responses and added resistances are compared with experimental data. For S175 containership, grid convergence test has been conducted to investigate the sensitivity of grid spacing on the motion responses and added resistances.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.6 no.3
• /
• pp.11-21
• /
• 2002

The evaluation of displacement ductility is performed by direct method through tracking the inelastic hysteretic behavior of RC bridge columns subject to cyclic loading using a flexibility-based fiber element mode. To reasonably track the inelastic behavior until the RC bridge column reaches its ultimate state, the average stress-average strain relations and joint elements, which agree well with experiments, are modified and applied considering the tension stiffening behavior and discontinuous displacement between the column and its base. In addition the evaluation of displacement ductility is performed by a direct method easily applicable to numerical analysis. Locations for the integration points, values for the post-crushing concrete strength and low-cycle fatigue failure of longitudinal reinforcement that affect the calculation of yielding and ultimate displacements are proposed for the application to flexibility-based fiber element model. Since less than 10% of error occurs during the displacement ductility analysis, the yielding and ultimate displacements evaluated by the applied analysis method and model appear to be valid.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.129-136
• /
• 2017

3D seismic data processing methods such as full waveform inversion or reverse-time migration require 3D wave propagation modeling and heavy calculations. We compared efficiency and accuracy of a Xeon Phi coprocessor to those of a high-end server CPU using 3D frequency-domain wave propagation modeling. We adopted the OpenMP parallel programming to the time-domain finite difference algorithm by considering the characteristics of the Xeon Phi coprocessors. We applied the Fourier transform using a running-integration to obtain the frequency-domain wavefield. A numerical test on frequency-domain wavefield modeling was performed using the 3D SEG/EAGE salt velocity model. Consequently, we could obtain an accurate frequency-domain wavefield and attain a 1.44x speedup using the Xeon Phi coprocessor compared to the CPU.

• Nuclear Engineering and Technology
• /
• v.49 no.2
• /
• pp.395-403
• /
• 2017

An efficient evaluation method for the probability of a tornado missile strike without using the Monte Carlo method is proposed in this paper. A major part of the proposed probability evaluation is based on numerical results computed using an in-house code, Tornado-borne missile analysis code, which enables us to evaluate the liftoff and flight behaviors of unconstrained objects on the ground driven by a tornado. Using the Tornado-borne missile analysis code, we can obtain a stochastic correlation between local wind speed and flight distance of each object, and this stochastic correlation is used to evaluate the conditional strike probability, $Q_V(r)$, of a missile located at position r, where the local wind speed is V. In contrast, the annual exceedance probability of local wind speed, which can be computed using a tornado hazard analysis code, is used to derive the probability density function, p(V). Then, we finally obtain the annual probability of tornado missile strike on a structure with the convolutional integration of product of $Q_V(r)$ and p(V) over V. The evaluation method is applied to a simple problem to qualitatively confirm the validity, and to quantitatively verify the results for two extreme cases in which an object is located just in the vicinity of or far away from the structure.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.37-46
• /
• 1991

A method is developed for solving the natural frequencies and mode shapes of linearly variable tapered beams. The governing differential equation for the tapered beam is derived. Three kinds of cross sectional shape are considered in differential equation. The Runge-Kutta method and the determinant search method are used to perform the integration of the differential equation and to determine the natural frequencies, respectively. The hinged-hinged, hinged-clamped, damped-clamped and free-damped end constraints are investigated in numerical examples. The lowest four nondimensional natural frequencies are obtained as functions of $d_b/d_a$. ratio. The effects of end constraints and cross sectional shapes on frequencies are analyzed and typical mode shapes are also presented.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.9-21
• /
• 1992

The effectiveness of the isolation system installed at the base of the primary structure and at the support of the substructure mounted on the primary structure is evaluated for reducing of structural responses under different earthquakes in this paper. The structural responses are analyzed to identify its behavior due to the input motion characteristics such as various peak acceleration and frequency content. Three analytical models are used to evaluate the effectiveness of the isolation system in this study as follows: fixed-base primary structure with support-fixed substructure, base-isolated primary structure with support-fixed substructure, and fixed-base primary structure with support-isolated substruciure. A computer code (KBISAP) is used for numerical integration of equation of motion considering the interaction between the primary structure and the secondary structure. The matrix condensation technique and constant average acceleration method are utilized in this program. And also, the effective stiffness of the base-isolator on reducing the structural response are evaluated for various earthquakes through the relationship of the acceleration - displacement.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.413-426
• /
• 2017

Probabilistic methods are used in engineering where a computational model contains random variables. The proposed method under development: Direct Optimized Probabilistic Calculation (DOProC) is highly efficient in terms of computation time and solution accuracy and is mostly faster than in case of other standard probabilistic methods. The novelty of the DOProC lies in an optimized numerical integration that easily handles both correlated and statistically independent random variables and does not require any simulation or approximation technique. DOProC is demonstrated by a collection of deliberately selected simple examples (i) to illustrate the efficiency of individual optimization levels and (ii) to verify it against other highly regarded probabilistic methods (e.g., Monte Carlo). Efficiency and other benefits of the proposed method are grounded on a comparative case study carried out using both the DOProC and MC techniques. The algorithm has been implemented in mentioned software applications, and has been used effectively several times in solving probabilistic tasks and in probabilistic reliability assessment of structures. The article summarizes the principles of this method and demonstrates its basic possibilities on simple examples. The paper presents unpublished details of probabilistic computations based on this method, including a reliability assessment, which provides the user with the probability of failure affected by statistically dependent input random variables. The study also mentions the potential of the optimization procedures under development, including an analysis of their effectiveness on the example of the reliability assessment of a slender column.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.2
• /
• pp.73-78
• /
• 2010

Wave energy is dissipated mainly by friction on the seabed until the waves reach the surf zone. Many researchers have investigated the mechanism of wave friction and the bottom shear stress induced by wave motion at a certain point is now well estimated by introducing the wave friction factor related to the near bed velocity given by linear wave theory. The variation of wave energy or wave height over a long distance can be, however, estimated by an iteration process when the propagation of waves is strongly influenced by bed friction. In the present study simple semi-theoretical equation has been developed to compute the variation of wave height for the condition of wave propagation on a constant beach slope. The ratio of wave height is determined by the product of shoalng factor and wave height friction factor (frictional wave energy dissipation factor). The wave height estimated by the new equation is compared with the wave height estimated by the solution of numerical integration for the condition that the waves propagate on a constant slope.