• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.5
• /
• pp.1084-1094
• /
• 1993

An efficient method for determining forced responses of a general linear structural system in time domain using subtructure modes and Lagrange multipliers is presented. Compared with the conventional mode synthesis methods, the suggested method does not construct the equations of motion of the combined whole structure and thus the modal parameters of the whole structure are not required. Only modal parameters of each substructure and geometric compatibility conditions are needed. Both the loaded interface free-free modes and free interface modes can be employed as the modal bases of each substructure. Recurrence discrete-time state equations based upon state transition matrix are formulated for the transient analysis of a parameter-changing system. It is shown form numerical examples that the suggested method is very accurate and efficient to calculate transient responses compares with the direct numerical integration method.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.423-434
• /
• 1995

It was used boundary element method(BEM) and analysed axisymmetric problem to solve hydraulic cylinder for large vessel acting uniform internal pressure(25N/m super(2)) within elastic limit. This paper was utilized the carbon steel tubes for machine structural purposed model, inner radius was 150mm and outer radius was 250mm, axial length was semi-infinite and the isoparametric element was used. The important results obtained in this study were summarized as follows. Radial, tangential and shearing stress occured the maximum stresses(48, -20 and 34MPa) at the inner radius and the minimum stresses(32, -4 and 18MPa) at the outer radius of the hydraulic cylinder for large vessel. But negative signs have meaning compressive stress and stress diminution ratio was about 0.15MPa/mm. The use of isoparametric element raised accuracy and the increment of input data lessened the error in internal point but computer run-time was increased. The double node was improved the internal solutions to settle discontinuity at corner and the double exponential formula lessened error of stress value at boundary neighborhood. And then coincidence between the analytical and exact results is found to be fairly good, showing that the proposed analytical by BEM is reliable.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.55-63
• /
• 1995

This paper describes the procedure to analyze the performance of the end-plated propeller(EPP) by a boundary integral method. The screw blade(SB) and end-plate(EP) are represented by a set of quadrilateral panels, where the source and normal dipole of uniform strength are distributed. The perturbation velocity potential, being the only unknown via the potential-based formulation, is determined by satisfying the flow tangency condition on the blade and the end-plate at the same time. The Kutta condition is satisfied through an iterative process by requiring the null pressure jump across the upper and lower sides of the trailing edges of both the SH and the EP. Sample calculations indicate that the EP increases the loading near the tip of the SB while spreading the trailing vortices along the trailing edge of the EP, thus avoiding the strong tip-vortex formation. Predicted performance of the EPP shows good correlations with the experimental results. The method is therefore considered applicable in designing and analyzing the EPP which may be an alternative for energy-saving propulsive devices.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.9 no.5
• /
• pp.47-57
• /
• 2000

In this study, a fundamental approach to make clear the mechanism of the mutual interference and coalescence of stress fields in the vicinity of two crack tips on the process of their slow growth, using boundary element method. Automatic generation of quadratic discontinuous elements along both of the crack boundaries which can be defined by an arbitrary piece-wise straight geometry. The direction of the crack-extension increment is predicted by the maximum principal stress criterion, corrected to account for the discreteness of the crack extension. Along the computed direction, the crack is extended one increment. Automatic incremental crack-extension analysis with no remeshing, computation of the stress intensity factors by J-integral. Numerical stress intensity factors for two growing cracks in plane-homogeneous regions were determined.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.4
• /
• pp.20-27
• /
• 2001

Recently numerical approaches for wave loads acting on the vertical caisson of breakwater, and resulting wave reflection and transmission coefficients have been performed. Although the numerical studies by Sulisz's(1997) and Kim et al.(2000) are suggested representatively, theoretical formulation for nonlinear wave pressure is not developed yet. And experimental results of Sulisz(1997) revealed that nonlinear uplift pressure on the caisson may be produced largely on the case of caisson founded on the high rubble mound. From the results of this study, the nonlinear theory for the uplift wave pressure acting on the caisson by applying boundary integral method of Green theorem is formulated, and also the characteristics of nonlinear uplift pressure and run-up height on the caisson are evaluated numerically, according to the variations of hydraulic properties of the rubble mound.

• Proceedings of the KIEE Conference
• /
• 1998.11a
• /
• pp.88-90
• /
• 1998

In this paper, the indirect BIEM(boundary integral equation method) is adopted to analyze 3-D eddy currents in cover plate of power transformer. In indirect BIEM, the equivalent magnetic surface charge density and the eqivalent magnanetic surface current density are the unknowns. Using triangular constant elements, the integral equations are discretized into boundary element equations of minimum order. Eddy currents are obtained in terms of euqivalent magnetic surface sources. And the locad overheating can be predicted using the eddy currents distribution in cover plate of power transformer.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.803-810
• /
• 1992

With increases in the rotational speed of hydraulic machine, studies on the hydrodynamic characteristics of supercavitating cascade are important on the view of flow analysis and design of fluid machinery. In the present paper, the complex functions of nonlinear theory corresponding to the flow of supercavitating cascade can be obtained by distributing singulary singulary points such as sources, vortexes and doublets on hydrofoil and free streamline. The numerical calculations on the closed wake model and semi-closed wake model are carried out in order to show the flow characteristics around the supecavitating cascade with finite with finite cavity length. As the result of this study, the flow characteristics such as lift, drag and cavitation coefficients are predicted by the flow conditions of supercavitating cascade in the fluid machinery.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.521-528
• /
• 2017

In this paper, the numerical analysis for heat conduction of silica/phenolic composite material, which is used for solid rocket nozzle liner or insulator, was conducted. 1-D Finite Difference Method for the analysis of silica/phenolic during the firing of solid rocket motor was used to calculate the heat conduction considering the surface ablation and the thermal decomposition. The boundary condition at the nozzle wall took into account the convective heat transfer, which was obtained by integration equation. The numerical results of the surface ablation and char depth were compared with the results of test motor that is TPEM-10. It was found that the result of calculation is favorably agreed with the thermal response of test motor.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.4
• /
• pp.76-84
• /
• 2018

In this paper, the numerical analysis for heat conduction of silica/phenolic composite material, used for solid rocket nozzle liners or insulators, is conducted. A 1-dimensional finite difference method for the analysis of silica/phenolic during the firing of a solid rocket motor is used to calculate heat conduction, considering surface ablation and thermal decomposition. The boundary condition at the nozzle wall, considering the convective heat transfer, is obtained via integration equations. The numerical results of the surface ablation and char depth are compared with the results of a TPEM-10 test motor, finding that the result of calculation agrees with the thermal response of the test motor.