• Communications of the Korean Mathematical Society
• /
• v.17 no.2
• /
• pp.349-361
• /
• 2002

The ill-posed elliptic wave propagation problems can be transformed into well-posed initial value problems of the reflection and transmission operators characterizing the material structure of the given model by the combination of wave field splitting and invariant imbedding methods. In general, the derived scattering operator equations are of first-order in range, nonlinear, nonlocal, and stiff and oscillatory with a subtle fixed and movable singularity structure. The phase space and path integral analysis reveals that construction and reconstruction algorithms depend crucially on a detailed symbol analysis of the scattering operators. Some information about the singularity structure of the scattering operator symbols is presented and analyzed in the transversely homogeneous limit.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.7
• /
• pp.1840-1850
• /
• 1994

Enhancement of numerical algorithms for low speed compressible flow will be considered. Contemporary time-marching algorithm has been widely accepted and applied as the method of choice for transonic, supersonic and hypersonic flows. In the low Mach number regime, time-marching algorithms do not fare as well. When the velocity is small, eigenvalues of the system of compressible equations differ widely so that the system becomes very stiff and the convergence becomes very slow. This characteristic can lead to difficulties in computations of many practical engineering problems. In the present approach, the time-derivative preconditioning method will be used to control the eigenvalue stiffness and to extend computational capabilities over a wide range of flow conditions (from very low Mach number to supersonic flow). Computational capabilities of the above algorithm will be demonstrated through computation of a variety of practical engineering problems.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.27-30
• /
• 2007

In recent years, scientific community has found renewed interest in hypersonic flight research. These hypersonic vehicles undergo severe aero-thermal environments during their flight regimes. One of the most important topics of research in hypersonic aerodynamics is to find a reasonable way of calculating either the surface temperature or the heat flux to surface when its temperature is held fixed. This requires modeling of physical and chemical processes. Hyperbolic system of equations with stiff relaxation method are being identified in recent literature as a novel method of predicting long time behavior of systems such as gas at high temperatures. In present work, Energy Relaxation Method (ERM) has been considered to simulate the real gas flow over a 2-D cylinder. Present heat flux results over the cylinder compared well with the experiment. Thus, real gas effects in hypersonic flows can be modeled through energy relaxation method.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.59-64
• /
• 1999

The Equations of Chemical kinetics ate very stiff, which forces the use of an implicit scheme. The problem of implicit scheme, however, is that the jacobian must be solved at each time step. In this paper, we examined the approximate chemical jacobian methods such as Gauss-Seidel, Jacobi partial jacobian and diagonalized jacobian that can be stable without full jacobian, We show that Gauss-Seidel jacobian method is stable and accurate as well as full jacobian and that this is more efficient in supersonic combustion problem about $20\%$ than the full jacobian method with same accuracy,

• Proceedings of the Korea Concrete Institute Conference
• /
• 2005.05a
• /
• pp.451-454
• /
• 2005

The ㄱ type perfobond rib shear connector is a ㄱ type flat steel plate with a number of holes punched through. This connector can be effectively used in girder with high shear. The ㄱ type perfobond rib shear connector exhibit very stiff behaviour under service load conditions and also had the characteristic of retaining a significant amount of load after the attainment of ultimate capacity. A regression analysis, which is based on a model that takes into account the contributions of concrete dowels formed by the rib holes, the transverse reinforcement, the strength of concrete in front of the rib, and the ㄱ type plate as well as a nonlinear finite element analysis, is used in the derivation. An empirical equation for the design of ㄱ type perfobond rib shear connector is proposed.

• Journal of KSNVE
• /
• v.9 no.6
• /
• pp.1259-1266
• /
• 1999

The objective of this investigation is to develop a compliant track link model and apply this model to the multi-body dynamic analysis of high mobility tracked vehicles. Two major difficulties encountered in developing the compliant track models. The first one is that the integration step size must be kept small in order to maintain the numerical stability of the solution. This solution deals with high oscillatory signals resulting from the impulsive contact forces and stiff compliant elements to represent the joints between the track links. The second difficulty is due to the large number of the system equations of motion of the three dimensional multibody tracked vehicle model. This problem was sloved by decoupling the equations of motion of the chassis subsystem and the track subsystems. Recursive methods are used to obtain a minimum set of equations for the chassis subsystem. Several simulation scenarios were tested for the high mobility tracked vehicle including accelaeration, high speed cruising, braking, and turning motion in order to demonstrate the effectiveness and validity of the methods proposed in this investigation.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.12 no.5
• /
• pp.23-29
• /
• 2008

This study investigates the dynamic instability of strength-limited bilinear single degree of freedom (SDF) systems under seismic excitation. The strength-limited bilinear hysteretic model best replicates the hysteretic behavior of the steel moment resisting frames. To estimate the dynamic instability of SDF systems, the collapse strength ratio is used, which is the yield-strength reduction factor when collapse occurs. Statistical studies are carried out to estimate median collapse strength ratios and those dispersions of strength-limited bilinear SDF systems with given natural periods, hardening stiffness ratios, post-capping stiffness ratios, ductility and damping ratios ranging from 2 to 20% subjected to 240 earthquake ground motions recorded on stiff soil sites. Equations to calculate median and standard deviation of collapse strength ratios in strength-limited bilinear SDF systems are obtained through nonlinear regression analysis. By using the proposed equations, this study estimated the probabilistic distribution of collapse strength ratios, and compared this with the exact values from which the accuracy of the proposed equations was verified.

• Structural Engineering and Mechanics
• /
• v.82 no.1
• /
• pp.121-131
• /
• 2022

The main goal of this study is to prepare a program for analyzing High Strength Steel Fibrous Reinforced Concrete (HSSFRC) slabs and predict the response and strength of the slab instead of preparing a prototype and testing it in the laboratory. For this purpose, new equations are proposed to represent the material properties of High Strength Steel Fibrous Reinforced Concrete. The proposed equations obtained from performing regression analysis on many experimental results using statistical programs. The finite element method is adopted for non-linear analysis of the slabs. The eight-node "Serendipity element" (3 DoF) is chosen to represent the concrete. The layered approach is adopted for concrete elements and the steel reinforcement is represented by a smeared layer. The compression properties of the concrete are modeled by a work hardening plasticity approach and the yield condition is determined depending on the first two stress invariants. A tensile strength criterion is adopted in order to estimate the cracks propagation. many experimental results for testing slabs are compared with the numerical results of the present study and a good agreement is achieved regarding load-deflection curves and crack pattern. The response of the load deflection curve is slightly stiff at the beginning because the creep effect is not considered in this study and for assuming perfect bond between the steel reinforcement and the concrete, however, a great agreement is achieved between the ultimate load from the present study and experimental results. For the models of the tension stiffening and cracked shear modulus, the value of Bg and Bt (Where Bg and Bt are the curvature factor for the cracked shear modulus and tension stiffening models respectively) equal to 0.005 give good results compared with experimental result.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.2
• /
• pp.243-250
• /
• 2003

The extinction behavior and the unsteady response of augmented reduced mechanism(ARM) have been investigated by adopting an OPPDIF code and a numerical solver for the flamelet equations. By comparing the performance of the ARM based on Miller and Bowman's mechanism(MB-ARM) with that of the ARM based on GRI-Mech 3.0(GRI-3.0-ARM), it is identified that the MB-ARM is more suitable for the unsteady calculation because it is relatively less stiff than GRI-3.0-ARM during an ignition process. The steady results using the MB-ARM, which is modified to predict reasonably the extinction point of experiment, are in excellent agreement with those from full mechanism. Under the sinusoidal transient disturbances of scalar dissipation rate, the unsteady responses of the flame temperature and species concentrations using a modified MB-ARM show in very close agreement with those from full mechanism. It is presumed that above modified MB-ARM is very suitable for the unsteady simulation of turbulent flames because it gives not only a low computational cost but also a good prediction performance for flame structure, extinction point and unsteady response.

• Steel and Composite Structures
• /
• v.15 no.4
• /
• pp.379-397
• /
• 2013

Conoidal shells are doubly curved stiff surfaces which are easy to cast and fabricate due to their singly ruled property. Application of laminated composites in fabrication of conoidal shells reduces gravity forces and mass induced forces compared to the isotropic constructions due to the high strength to weight ratio of the material. These light weight shells are preferred in the industry to cover large column free open spaces. To ensure design reliability under service conditions, detailed knowledge about different behavioral aspects of conoidal shell is necessary. Hence, in this paper, static bending, free and forced vibration responses of composite conoidal shells are studied. Lagrange's equation of motion is used in conjunction with Hamilton's principle to derive governing equations of the shell. A finite element code using eight noded curved quadratic isoparametric elements is developed to get the solutions. Uniformly distributed load for static bending analysis and three different load time histories for solution of forced vibration problems are considered. Eight different stacking sequences of graphite-epoxy composite and two different boundary conditions are taken up in the present study. The study shows that relative performances of different shell combinations in terms of static behaviour cannot provide an idea about how they will relatively behave under dynamic loads and also the fact that the points of occurrence of maximum static and dynamic displacement may not be same on a shell surface.