• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.317-330
• /
• 2008

Stress intensity factors for a planar crack parallel to a bimaterial interface are considered. The formulation leads to a system of hypersingular integral equations whose unknowns are three modes of crack opening displacements. In the numerical analysis, the unknown displacement discontinuities are approximated by the products of the fundamental density functions and polynomials. The numerical results show that the present method yields smooth variations of stress intensity factors along the crack front accurately. The mixed mode stress intensity factors are indicated in tables and figures with varying the shape of crack, distance from the interface, and elastic constants. It is found that the maximum stress intensity factors normalized by root area are always insensitive to the crack aspect ratio. They are given in a form of formula useful for engineering applications.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.53-59
• /
• 2007

An axisymmetric supersonic screeching jet is numerically simulated to examine the length scales of screech frequency as well as screech tone generation mechanism. The axisymmetric Reynolds-averaged Navier-Stokes equations in conjuction with a modified Spalart-Allmaras turbulence model are employed. It is demonstrated that the axisymmetric jet screech tones can be simulated correctly and the numerical results are in good agreement with the experimental data. Instability waves, shock-cell structures and the phenomena of shock motion are investigated in detail to identify the screech tone generation mechanism. Shock spacings and standing wave length are analyzed to determine the dominent length scale crucial to the screech frequency formulation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.175-182
• /
• 2005

The present study is concerned with the application of dynamic relaxation method in the investigation of the large deflection behavior of spatial structures. This numerical algorithm do not require the computation or formulation of any tangent stiffness matrix. The convergence to the solution is achieved by using only vectorial quantities and no stiffness matrix is required in its overall assembled form. In an effort to evaluate the merits of the methods, extensive numerical studies were carried out on a number of selected structural systems. The advantages of using dynamic relaxation methods, in tracing the post-buckling behavior of spatial structures, are demonstrated.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.03a
• /
• pp.125-132
• /
• 2002

In this study, a numerical method is developed for nonlinear analysis for soil-pile-structure interaction system in time domain. Finite elements considering material nonlinearity are used for the near field and boundary elements for the far field. In the near field, frame elements are used for modeling a pile and plane-strain elements for surrounding soil and superstructure. In. the far field, boundary element formulation using the dynamic fundamental solution is adopted and coupled with the near field. Transformation of stiffness matrices of boundary elements into time domain is performed by inverse FFT. Stiffness matrices in the near field and far field are coupled. Newmark direct time integration method is applied. Developed soil-pile-structure interaction analysis method is verified with available literature and commercial code. Also, parametric studies by developed numerical method are performed. And seismic response analysis is performed using actual earthquake records.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.03a
• /
• pp.117-124
• /
• 2002

The major purpose of this study is to determine the dynamic behavior of soil-pile-structure interaction system considering the underground cavity. For the analysis, a numerical method fur ground response analysis using FE-BE coupling method is developed. The total system is divided into two parts so called far field and near field. The far field is modeled by boundary element formulation using the multi-layered dynamic fundamental solution that satisfied radiational condition of wave. And this is coupled with near field modeled by finite elements. For the verification of dynamic analysis in the frequency domain, both forced vibration analysis and free-field response analysis are performed. The behavior of soil non-linearity is considered using the equivalent linear approximation method. As a result, it is shown that the developed method can be an efficient numerical method to solve the seismic response analysis considering the underground cavity in 2D problem.

• Steel and Composite Structures
• /
• v.18 no.4
• /
• pp.889-907
• /
• 2015

This paper aims to present an alternative analytical method for transient vibration analysis of doubly-curved laminated shells subjected to dynamic loads. In the method proposed, the governing differential equations of laminated shell are derived using the dynamic version of the principle of virtual displacements. The governing equations of first order shear deformation laminated shell are obtained by Navier solution procedure. Time-dependent equations are transformed to the Laplace domain and then Laplace parameter dependent equations are solved numerically. The results obtained in the Laplace domain are transformed to the time domain with the help of modified Durbin's numerical inverse Laplace transform method. Verification of the presented method is carried out by comparing the results with those obtained by Newmark method and ANSYS finite element software. Also effects of number of laminates, different material properties and shell geometries are discussed. The numerical results have proved that the presented procedure is a highly accurate and efficient solution method.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.165-176
• /
• 2015

A numerical analysis was conducted to investigate the inlet buzz and combustion oscillation in an axisymmetric ramjet engine with wedge-type flame holders. The physical model of concern includes the entire engine flow path, extending from the leading edge of the inlet center-body through the exhaust nozzle. The theoretical formulation is based on the Farve-averaged conservation equations of mass, momentum, energy, and species concentration, and accommodates finite-rate chemical kinetics and variable thermo-physical properties. Turbulence closure is achieved using a combined scheme comprising of a low-Reynolds number k-${\varepsilon}$ two-equation model and Sarkar's compressible turbulence model. Detailed flow phenomena such as inlet flow aerodynamics, flame evolution, and acoustic excitation as well as their interactions, are investigated. Mechanisms responsible for driving the inlet buzz are identified and quantified for the engine operating at subcritical conditions.

• Steel and Composite Structures
• /
• v.27 no.1
• /
• pp.27-33
• /
• 2018

In this study, a mixed-interpolated tensorial component 4 nodes method (MITC4) is treated as a numerical analysis model for topology optimization using multiple materials assigned within Reissner-Mindlin plates. Multi-material optimal topology and shape are produced as alternative plate retrofit designs to provide reasonable material assignments based on stress distributions. Element density distribution contours of mixing multiple material densities are linked to Solid Isotropic Material with Penalization (SIMP) as a design model. Mathematical formulation of multi-material topology optimization problem solving minimum compliance is an alternating active-phase algorithm with the Gauss-Seidel version as an optimization model of optimality criteria. Numerical examples illustrate the reliability and accuracy of the present design method for multi-material topology optimization with Reissner-Mindlin plates using MITC4 elements and steel materials.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.350-357
• /
• 1997

The transient incompressible flow behind the axisymmetric bluff body is numerically simulated using the random vortex method(RVM). Based on the vorticity formulation of the unsteady Navier-Stokes equations, the Lagrangian approach with a stochastic simulation of diffusion using random walk technique is employed to account for the transport processes of the vortex elements. The numerical solutions for 2-dimensional recirculating flow behind a backward-facing step in the laminar range of Reynolds number are compared with experimental data. The present simulation focuses on the transitional flow regime where the recirculation zone behind the bluff body becomes highly unsteady and large-scale vortex eddies are shed from the bluff body wake due to intrinsic shear layer instabilities. The unsteady vertical flow structures and the mixing characteristics behind the bluff body are discussed in detail.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.193-200
• /
• 2009

Recently, there have been efforts to construct hybrids among the existing methodologies for multiphase flow such as VOF, Level Set, and Front Tracking with the intention of facilitating simulations of general three-dimensional problems. As one of the hybrid method, we have developed the Level Contour Reconstruction Method (LCRM) for general three-dimensional multiphase flows including phase change. The main idea was focused on simplicity and a robust algorithm especially for the three-dimensional case. It combines characteristics of both Front Tracking and Level Set methods. While retaining an explicitly tracked interface using interfacial elements, the calculation of a vector distance function plays a crucial role in the periodic reconstruction of the interface elements in the LCRM method to maintain excellent mass conservation and interface fidelity. In addition, compact curvature formulation is incorporated for the calculation of the surface tension force thereby reducing parasitic currents to a negligible level.