• Journal of the Korea Concrete Institute
• /
• v.22 no.4
• /
• pp.595-603
• /
• 2010

Monte-Carlo simulation technique is often used in order to predict service life of concrete structure subjected to chloride penetration in marine environment based on probability theory. Monte-Carlo simulation method, however, the method gives different results every time that the simulation is run. On the other hand, moment method, which is frequently used in reliability analysis, needs negligible computational cost compared with simulation technique and gives a constant result for the same problem. Thus, in this study, moment method was applied to the calculation of corrosion-initiation probability. For this purpose, computer programs to calculate failure probabilities are developed using first-order second moment (FOSM) and second-order second moment (SOSM) methods, respectively. From the analysis examples with the developed programs, SOSM was found to give a more accurate result than FOSM does. The sensitivity analysis has shown that the factor affecting the corrosion-initiation probability the most was the cover depth, and the corrosion-initiation probability was influenced more by its coefficient of variation than its mean value.

• Smart Structures and Systems
• /
• v.18 no.2
• /
• pp.355-374
• /
• 2016

This paper presents and compares a one-dimensional (1D) bending theory for piezoelectric thin beam-type structures with resistive-inductive electrodes to ANSYS$^{(R)}$ three-dimensional (3D) finite element (FE) analysis. In particular, the lateral deflections and vibrations of slender piezoelectric beams are considered. The peculiarity of the piezoelectric beam model is the modeling of electrodes in such a manner that is does not fulfill the equipotential area condition. The case of ideal, perfectly conductive electrodes is a special case of our 1D model. Two-coupled partial differential equations are obtained for the lateral deflection and for the voltage distribution along the electrodes: the first one is an extended Bernoulli-Euler beam equation (second-order in time, forth order in space) and the second one the so-called Telegrapher's equation (second-order in time and space). Analytical results of our theory are validated by 3D electromechanically coupled FE simulations with ANSYS$^{(R)}$. A clamped-hinged beam is considered with various types of electrodes for the piezoelectric layers, which can be either resistive and/or inductive. A natural frequency analysis as well as quasi-static and dynamic simulations are performed. A good agreement between the extended beam theory and the FE results is found. Finally, the practical relevance of this type of electrodes is shown. It is found that the damping capability of properly tuned resistive or resistive-inductive electrodes exceeds the damping performance of beams, where the electrodes are simply linked to an optimized impedance.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.1063-1067
• /
• 2005

A second-order accuracy upwind scheme is used to investigate the run-up heights of tsunamis in the East Sea and the predicted results are compared with field observed data and results of a first-order accuracy upwind scheme, In the numerical model, the governing equations solved by the finite difference scheme are the linear shallow-water equations in deep water and nonlinear shallow-water equations in shallow water The target events is 1993 Hokktaido Nansei Oki Tsunami. The predicted results represent reasonably the run-up heights of tsunamis in the East Sea. And, The results of simulation is used to design inundation map.

• Journal of Korea Technical Association of The Pulp and Paper Industry
• /
• v.31 no.3
• /
• pp.35-46
• /
• 1999

Mathematical modeling provided a systematic analysis for the dynamic behavior of stock fluid in a paper machine pressurized headbox. Dynamic responses of liquid level, sheet basis weight and hydraulic pressure were predicted from the simulation model which represents the system. A unit step and asinusoidal wave load were considered as the input forcing functions in the headbox. Results are summarized as follows : 1. The dependence of sheet basis weight on liquid level in the pressurized-headbox was non -linear. 2. Liquid level in the head-box showed first-order lag with a unit step forcing to fluid input rate ; 3 . The amplitude of wave response of liquid level was inversely proportional to the time content for the sinusoidal input changes ; 4.Sheet basis weight showed second-order oscillating underamped responses for the step input load of flow rate ; 5. The damping factor in the second-order system was a function of air-pressure in the headbox ; and, 6. Dead-time existed in the measuring process for the headbox slice pressure.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.46-51
• /
• 1996

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. The fourth -order artificial damping is added to the continuity equation for numerical stability. A H-H type multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. The simulation shows three dimensional vortex-pair just behind body. And the flow separation is also observed the rear of the body. It has concluded that the results of present study properly agree with physical flow phenomena.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.4
• /
• pp.1-8
• /
• 1998

In order to solve the conflict problem between the ride comfort and the road holding, the optimal design of shock absorber that minimizes the r.m.s. of sprung mass vertical acceleration and pitch rate with the understeer characteristics constraints in the high speed stability is proposed. The design of experiments and the nonlinear optimization algorithm are used together to obtain the optimal design of shock absorber. The second order regression models of the input variables(front and rear damping coefficients) and the output variables (ride comfort index and road holding one) are obtained by the central composite design in the design of experiments. Then the optimal design of shock absorber can be systematically adjusted with applying the nonlinear optimization algorithm to the obtained second order regression model. The frequency response analysis of sprung mass acceleration and pitch rate shows the effectiveness of the proposed optimal design of shock absorber in the sprung mass resonance range with the understeer characteristics constraints.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.4
• /
• pp.1-8
• /
• 2000

The time simulation of motion responses of dolphin-moored BMP in waves is presented. The equation of motion based on Cummin's theory of impulse responses are employed, and solved in time domain by using the Newmark $\beta$ method. The hydrodynamic coefficient and first order wave exciting forces involved in the equations are obtained from a three-dimensional panel method in the frequency domain. The second order wave drift forces and mooring for dolphin system are taken into account. As for numerical example, time domain analysis are carried out for a BMP in irregular wave condition.

• Journal of the Korean Society of Safety
• /
• v.35 no.3
• /
• pp.49-57
• /
• 2020

In comparison with the existing static reliability analysis methods, the dynamic reliability analysis(DyRA) method is more suitable for estimating the failure probability of a structure subjected to earthquake excitations because it can take into account the frequency characteristics and damping capacity of the structure. However, the DyRA is known to have an issue of numerical stability due to the uncertainty in random sampling of the earthquake excitations. In order to solve this numerical stability issue in the DyRA approach, this study proposed two earthquake-scale factors. The first factor is defined as the ratio of the first earthquake excitation over the maximum value of the remaining excitations, and the second factor is defined as the condition number of the matrix consisting of the earthquake excitations. Then, we have performed parametric studies of two factors on numerical stability of the DyRA method. In illustrative example, it was clearly confirmed that the two factors can be used to verify the numerical stability of the proposed DyRA method. However, there exists a difference between the two factors. The first factor showed some overlapping region between the stable results and the unstable results so that it requires some additional reliability analysis to guarantee the stability of the DyRA method. On the contrary, the second factor clearly distinguished the stable and unstable results of the DyRA method without any overlapping region. Therefore, the second factor can be said to be better than the first factor as the criterion to determine whether or not the proposed DyRA method guarantees its numerical stability. In addition, the accuracy of the numerical analysis results of the proposed DyRA has been verified in comparison with those of the existing first-order reliability method(FORM), Monte Carlo simulation(MCS) method and subset simulation method(SSM). The comparative results confirmed that the proposed DyRA method can provide accurate and reliable estimation of the structural failure probability while maintaining the superior numerical efficiency over the existing methods.

• Earthquakes and Structures
• /
• v.9 no.3
• /
• pp.481-501
• /
• 2015

Current design codes and technical recommendations often provide rough indications on how to assess effective stiffness of Reinforced Concrete (R/C) frames subjected to seismic loads, which is a key factor when a linear analysis is performed. The Italian design code (NTC-2008), Eurocode 8 and ACI 318 do not take into account all the structural parameters affecting the effective stiffness and this may not be on the safe side when second-order $P-{\Delta}$ effects may occur. This paper presents a study on the factors influencing the effective stiffness of R/C beams, columns and walls under seismic forces. Five different approaches are adopted and analyzed in order to evaluate the effective stiffness of R/C members, in accordance with the scientific literature and the international design codes. Furthermore, the paper discusses the outcomes of a parametric analysis performed on an actual R/C building and analyses the main variables, namely reinforcement ratio, axial load ratio, concrete compressive strength, and type of shallow beams. The second-order effects are investigated and the resulting displacements related to the Damage Limit State (DLS) under seismic loads are discussed. Although the effective stiffness increases with steel ratio, the analytical results show that the limit of 50% of the initial stiffness turns out to be the upper bound for small values of axial-load ratio, rather than a lower bound as indicated by both Italian NTC-2008 and EC8. As a result, in some cases the current Italian and European provisions tend to underestimate second-order $P-{\Delta}$ effects, when the DLS is investigated under seismic loading.

• Structural Engineering and Mechanics
• /
• v.11 no.4
• /
• pp.423-442
• /
• 2001

This paper presents a method of elastic-plastic analysis for planar steel frames that provides the accuracy of distributed plasticity methods with the computational efficiency that is greater than that of distributed plasticity methods but less than that of plastic-hinge based methods. This method accounts for the effect of spread of plasticity accurately without discretization through the cross-section of a beam-column element, which is achieved by the following procedures. First, nonlinear equations describing the relationships between generalized stresses and strains of the cross-section are derived analytically. Next, nonlinear force-deformation relationships for the beam-column element are obtained through lengthwise integration of the generalized strains. Elastic-plastic flexibility coefficients are then calculated by differentiating the above element force-deformation relationships. Finally, an elastic-plastic stiffness matrix is obtained by making use of the flexibility-stiffness transformation. Adding the conventional geometric stiffness matrix to the elastic-plastic stiffness matrix results in the tangent stiffness matrix, which can readily be used to evaluate the load carrying capacity of steel frames following standard nonlinear analysis procedures. The accuracy of the proposed method is verified by several examples that are sensitive to the effect of spread of plasticity.