• 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.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.634-639
• /
• 2006

Active control method is applied to a flexible beam excited by a shock impulse in order to reduce the residual vibrations after the shock event. It is assumed that the shock input can be measured and is always occurred on the same point of the beam. If the system is well identified and the corresponding inverse system is designed reliably, it has shown that a very simple feed-forward active control method may be applied to suppress the residual vibrations without using error sensors and adaptive algorithm. Both numerical simulations and experimental results show a promising possibility of applying to a practical problem. Also, the performance of the method is examined by considering various practical aspects : shock duration, shock magnitude, and control point.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.6 s.111
• /
• pp.651-657
• /
• 2006

Active control method is applied to a flexible beam excited by a shock impulse in order to reduce the residual vibrations after the shock event. It is assumed that the shock input can be measured and is always occurred on the same point of the beam. If the system is well identified and the corresponding inverse system is designed reliably, it has shown that a very simple feed-forward active control method may be applied to suppress the residual vibrations without using error sensors and adaptive algorithm. Both numerical simulations and experimental results show a promising Possibility of applying to a practical problem. Also, the performance of the method is examined by considering various practical aspects : shock duration, shock magnitude, and control point.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.26 no.3
• /
• pp.292-298
• /
• 2017

Numerical analysis for various design parameters should be preceded by optimal design of composite materials. Numerous studies have been conducted on the bolting of interconnecting beams. In this study, the response surface method was applied to optimize the design of bolted joints connected by laminated wood composite beams. The response surface was created by combining the FEA code for composite analysis and the algorithm for forming the response surface. Optimization on this response surface was performed with a genetic algorithm to derive the results. The determination of the optimum bolt-hole position for the connection of composite beams is an optimization problem. Tsai-Wu composite failure index, maximum deflection, and simple von Mises stress are set as the objective functions. It has been proved that the design results of the optimized bolt-hole are superior to the design performance of the existing conventional bolt-hole position.

• Structural Engineering and Mechanics
• /
• v.8 no.5
• /
• pp.439-452
• /
• 1999

A simple numerical scheme suitable for tracing the fracture propagation path for structures idealized by means of Hillerborg's classical cohesive crack model is presented. A direct collocation, multidomain boundary element method is adopted for the required space discretization. The algorithm proposed is necessarily iterative in nature since the crack itinerary is a priori unknown. The fracture process is assumed to be governed by a path-dependent generally nonlinear softening law. The potentialities of the method are illustrated through two examples.

• Journal of the Chungcheong Mathematical Society
• /
• v.20 no.1
• /
• pp.11-21
• /
• 2007

Given a nonlinear function f : $\mathbb{R}{\rightarrow}\mathbb{R}$ that has a simple real zero ${\alpha}$, a new numerical method to be called k-fold pseudo- Halley's method is proposed and it's error analysis is under investigation to confirm the convergence behavior near ${\alpha}$. Under the assumption that f is sufficiently smooth in a small neighborhood of ${\alpha}$, the order of convergence is found to be at least k+3. In addition, the corresponding asymptotic error constant is explicitly expressed in terms of k, ${\alpha}$ and f as well as the derivatives of f. A zero-finding algorithm is written and has been successfully implemented for numerous examples with Mathematica.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.

• Journal of Advanced Marine Engineering and Technology
• /
• v.21 no.3
• /
• pp.236-245
• /
• 1997

The flow characteristics of developing turbulent steady flow are investigated numerically and experimentally in the entrance region of a square duct ($40 mm{\times}40 mm$ and 4, 000 mm). The numerical anaysis are incorporated by finite- volume discretization with staggered grid system and SIMPLE algorithm. The numerical solution are compared with experimental results of mean velocity profiles, turbulence intensity and entrance length. For turbulent steady flow, the turbulent components in the velocity waveforms increase as the dimensionless transverse position approaches the wall. Thrbulence intensity increases as the dimensionless transverse position increases from the center to the wall of the duct for the developing turbulent steady flows. The entrance length of the turbulent steady flow is about 40 times as large as the hydraulic diameter under the present experimental condition.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.3
• /
• pp.396-405
• /
• 1995

Numerical calculations are carried out for the swirling turbulent flow in a pipe. Calculations are made for the flow with swirl parameter of 2.25 and the Reynolds number of 24,300. The turbulence closure models used in these calculations are two different types of Reynolds stress model, and the results are compared with those of $k-{\varepsilon}$ model and the experimental data. The finite volume method is used for the discretization, and the power-law scheme is employed as a numerical scheme. The SIMPLE algorithm is used for velocity-pressure correction. The computational results show that GL model gives the results better than those of SSG model in the predictions of velocity and stress components.