• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.6
• /
• pp.401-406
• /
• 2011

In this study, generation methods of oblique incident wave are newly proposed and examined using the fully non-linear numerical model with non-reflected wave generation system(LES-WASS-3D). In order to verify, free surface elevation and horizontal velocities are compared with $3^{rd}$ -order Stokes wave theory in 3-D oblique incident wave field. As a results, it is revealed that the numerical results by newly proposed technique are in good agreement with the theory.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.7 s.238
• /
• pp.814-819
• /
• 2005

The phase change heat transfer has been applied to the processes of machines as well as of manufacturing. The cycle in a heat exchanger includes the phase change phenomena of coolant for air conditioning, the solidification in casting process makes use of the characteristics of phase change of metal, and the welding also proceeds with melting and solidification. To predict the phase change processes, the experimental and numerical approaches are available. In the case of numerical analysis, the Enthalpy method is most widely applied to the phase change problem, comparing to the other numerical methods, i.e. the Equivalent Specific Heat method and the Temperature Recovery method. It's because that the Enthalpy method is accurate and straightforward. The Enthalpy method does not include any correction step while the correction of final temperature field is inevitable in the Equivalent Specific Heat method and the Temperature Recovery method. When the temperature field is to be used in the calculation, however, there must be converting process from enthalpy to temperature in the calculation scheme of Enthalpy method. In this study, an improved method for the Equivalent Specific Heat method is introduced whose method dose not include the correction steps and takes temperature as an independent variable so that the converting between enthalpy and temperature does not need any more. The improved method is applied to the solidification process of pure metal to see the differences of conventional and improved methods.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.565-570
• /
• 2007

Fiber reinforced polymer(FRP) composite decks are new to bridge applications and hence not much literature exists on their structural mechanical behavior. As there are many differences between numerical displacements through static analysis of the primary model and experimental displacements through static load tests, system identification (SI)techniques such as Neural Networks (NN) and support vector machines (SVM) utilized in the optimization of the FE model. During the process of identification, displacements were used as input while stiffness as outputs. Through the comparison of numerical displacements after SI and experimental displacements, it can note that NN and SVM would be effective SI methods in modeling an FRP deck. Moreover, two methods such as response surface method and iteration were proposed to optimize the estimated stiffness. Finally, the results were compared through the mean square error (MSE) of the differences between numerical displacements and experimental displacements at 6 points.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.177-184
• /
• 1997

In this paper, numerical methods are developed for solving the elastica of simple beams with variable-arc-length subjected to a point loading. The beam model is based on Bernoulli-Euler beam theory. The Runge-Kutta and Regula-Falsi methods, respectively, are used to solve the governing differential equations and to compute the beam's rotation at the left end of the beams. Extensive numerical results of the elastica responses, including deflected shapes, rotations of cross-section and bending moments, are presented in non-dimensional forms. The possible maximum values of the end rotation, deflection and bending moment are determined by analyzing the numerical data obtained in this study.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.691-697
• /
• 2007

Numerical eigen analysis of beam-like structure can be easily and effectively done by various conventional beam theory-based methods. However, in case of the structures composed of composite-sectioned beams, the application of conventional numerical methods requires one to derive both equivalent material and geometry properties. In the present paper, these equivalent properties are derived by the transformed section method and the test FEM program is coded. The numerical accuracy of the proposed method is verified through the comparison with the ANSYS 3-D model.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Journal of Electrical Engineering and Technology
• /
• v.13 no.3
• /
• pp.1110-1122
• /
• 2018

This paper proposes an alternative method to evaluate the effect of wind power to the power system stability with small disturbance. Alternatively, available techniques for stability analysis of a power system based on deterministic methods are less accurate for high penetration of wind power. Numerical simulations of random behaviors are computationally expensive. A stochastic stability index (SSI) is proposed for the power system stability evaluation based on the theory of stochastic stability and energy function, specifically the stochastic derivative of the relative well-defined energy function and the critical energy. The SSI is implemented on the modified nine-bus system including wind turbines under different conditions. A doubly-fed induction generator (DFIG) wind turbine is characterized and modeled using measured wind data from several sites in Thailand. Each of the obtained wind power data is analyzed. The wind power effect is modeled considering the aggregated effect of wind turbines. With the proposed method, the system behavior is properly predicted and the stability is quantitatively evaluated with less computational effort compared with conventional numerical simulation methods.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2084-2093
• /
• 2022

Probabilistic safety assessment is widely used to quantify the risks of nuclear power plants and their uncertainties. When the lognormal distribution describes the uncertainties of basic events, the uncertainty of the top event in a fault tree is approximated with the sum of lognormal random variables after minimal cutsets are obtained, and rare-event approximation is applied. As handling complicated analytic expressions for the sum of lognormal random variables is challenging, several approximation methods, especially Monte Carlo simulation, are widely used in practice for uncertainty analysis. In this study, a theoretical approach for analyzing the sum of lognormal random variables using an efficient numerical integration method is proposed for uncertainty analysis in probability safety assessments. The change of variables from correlated random variables with a complicated region of integration to independent random variables with a unit hypercube region of integration is applied to obtain an efficient numerical integration. The theoretical advantages of the proposed method over other approximation methods are shown through a benchmark problem. The proposed method provides an accurate and efficient approach to calculate the uncertainty of the top event in probabilistic safety assessment when the uncertainties of basic events are described with lognormal random variables.

• Structural Engineering and Mechanics
• /
• v.88 no.5
• /
• pp.493-500
• /
• 2023

The multiparameter eigenvalue method can be used to solve the damped finite element model updating problems. This method transforms the original problems into multiparameter eigenvalue problems. Comparing with the numerical methods based on various optimization methods, a big advantage of this method is that it can provide all possible choices of physical parameters. However, when solving the transformed singular multiparameter eigenvalue problem, the proposed method based on the generalised inverse of a singular matrix has some computational challenges and may fail. In this paper, more details on the transformation from the dynamic model updating problem to the multiparameter eigenvalue problem are presented and the structure of the transformed problem is also exposed. Based on this structure, the rigorous mathematical deduction gives the upper bound of the number of possible choices of the physical parameters, which confirms the singularity of the transformed multiparameter eigenvalue problem. More importantly, we present a row and column compression method to overcome the defect of the proposed numerical method based on the generalised inverse of a singular matrix. Also, two numerical experiments are presented to validate the feasibility and effectiveness of our method.