• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.10a
• /
• pp.497-500
• /
• 1995

In this paper, a program for real time simulation of a vehicle is developed. The program uses relative coordinates and BEF(Backward Difference Formula) numerical integration method. Numerical tests showed that the proposed implicit method is more stable in carring out the numerical integration for vehicl dynamics than the explicit method. Hardware requirements for real time simulation are suggested. Algorithms of parallel processing is developed with DSP (digital signal processor).

• Elastomers and Composites
• /
• v.53 no.4
• /
• pp.207-212
• /
• 2018

Taylor rod impact tests have been the subject of many theoretical and experimental investigations. This paper discusses the numerical methods for simulating the Taylor impact test, which is widely used to obtain constitutive equations and failure conditions under high-velocity collisions of materials. With this in mind, a particle-based MPM (material point method) for linear viscoelastic solid materials was implemented, and MPM simulations for viscoelastic deformation behavior were numerically verified and confirmed by comparing the MPM and FEM results. In addition, this modeling and numerical approach could be extended to more complex viscoelastic models for basic understanding and to analyze the deformation and fracture behavior of more complicated viscoelastic material systems.

• Journal of KSNVE
• /
• v.9 no.2
• /
• pp.310-316
• /
• 1999

4-Pole parameter method based on an acoustic theory is very popular for the analysis of the acoustic behavior of the car exhaust system. However, this method is applicable only for the simple shape of acoustic elements of the muffler. Numerical methods such as FEM(Finite Element Method) or BEM(Boundary Element Method) can also provide acceptable results for the acoustic analysis of the car exhaust system. Even though these numerical methods have benefits for the analysis of complicated shape of acoustic elements of the muffler, time consuming is another problem during modeling and numerical calculation. Combining benefits of both methods, the new code called the hybrid method for car exhaust system is introduced. And the developed code is utilized for calculation of the transmission loss of a main muffler of an automobile comparing with the experimental results.

• Journal of the Korean Geosynthetics Society
• /
• v.17 no.4
• /
• pp.41-52
• /
• 2018

Since the composite ground design method is easy to apply for calculation or numerical analysis, it is applied to the design of cement mixing methods. However, the comparison studies between analysis and actual results such as a trial test and construction for the cement mixing method are few because the composite ground design method was developed for the compaction pile (SCP, GCP) methods. In this study, the results of various analysis methods, such as the composite ground analysis method (1 case) and the individual pile method (3 cases), were compared with actual measurements through a two-dimensional finite element numerical analysis. In case of the surface settlements, the results of study show that the individual plate method was larger than the actual measurements, while other methods are similar. The settlements at the under ground of the improved area is overestimated in all analysis methods. When comparing numerical analysis results for the horizontal displacement, and ground reaction forces, the individual pile method in equivalent wall concept was found to be the most similar to the measurements. The composite ground method was not able to predict the behavior of stress transfer (Arching effect) and it turned out that the prediction of horizontal displacement was too large.

• Kyungpook Mathematical Journal
• /
• v.60 no.3
• /
• pp.629-645
• /
• 2020

A uniformly convergent numerical method is developed for solving singularly perturbed 1-D parabolic convection-diffusion problems. The developed method applies a non-standard finite difference method for the spatial derivative discretization and uses the implicit Runge-Kutta method for the semi-discrete scheme. The convergence of the method is analyzed, and it is shown to be first order convergent. To validate the applicability of the proposed method two model examples are considered and solved for different perturbation parameters and mesh sizes. The numerical and experimental results agree well with the theoretical findings.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.26 no.1
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.

• Structural Engineering and Mechanics
• /
• v.56 no.4
• /
• pp.589-603
• /
• 2015

The mixed-mode stress intensity factors of 2-D angled cracks are evaluated by Petrov-Galerkin natural element (PG-NE) method in which Voronoi polygon-based Laplace interpolation functions and CS-FE basis functions are used for the trial and test functions respectively. The interaction integral is implemented in a frame of PG-NE method in which the weighting function defined over a crack-tip integral domain is interpolated by Laplace interpolation functions. Two Cartesian coordinate systems are employed and the displacement, strains and stresses which are solved in the grid-oriented coordinate system are transformed to the other coordinate system aligned to the angled crack. The present method is validated through the numerical experiments with the angled edge and center cracks, and the numerical accuracy is examined with respect to the grid density, crack length and angle. Also, the stress intensity factors obtained by the present method are compared with other numerical methods and the exact solution. It is observed from the numerical results that the present method successfully and accurately evaluates the mixed-mode stress intensity factors of 2-D angled cracks for various crack lengths and crack angles.

• Structural Engineering and Mechanics
• /
• v.88 no.3
• /
• pp.239-249
• /
• 2023

The Radial Point Interpolation Method (RPIM) has been proposed to overcome the difficulties associated with the use of the Radial Basis Functions (RBFs). The RPIM has the following properties: Simple implementation in terms of boundary conditions as in the Finite Element Method (FEM). A less expensive CPU time compared to other collocation meshless methods such as the Moving Least Square (MLS) collocation method. In this work, we propose an adaptive high-order numerical algorithm based on RPIM to simulate the thermoviscoplastic behavior of a material mixing observed in the Friction Stir Welding (FSW) process. The proposed adaptive meshfree RPIM algorithm adapts well to the geometric and physical data by choosing a good shape parameter with a good precision. Our numerical approach combines the RPIM and the Asymptotic Numerical Method (ANM). A numerical procedure is also proposed in this work to automatically determine an improved shape parameter for the RBFs. The efficiency of the proposed algorithm is analyzed in comparison with an iterative algorithm.

• Geomechanics and Engineering
• /
• v.15 no.6
• /
• pp.1207-1217
• /
• 2018

Geomechanics parameters are critical to numerical simulation, stability analysis, design and construction of geotechnical engineering. Due to the limitations of laboratory and in situ experiments, back analysis is widely used in geomechancis and geotechnical engineering. In this study, a hybrid back analysis method, that coupling numerical simulation, response surface (RS) and relevance vector machine (RVM), was proposed and applied to identify geomechanics parameters from hydraulic fracturing. RVM was adapted to approximate complex functional relationships between geomechanics parameters and borehole pressure through coupling with response surface method and numerical method. Artificial bee colony (ABC) algorithm was used to search the geomechanics parameters as optimal method in back analysis. The proposed method was verified by a numerical example. Based on the geomechanics parameters identified by hybrid back analysis, the computed borehole pressure agreed closely with the monitored borehole pressure. It showed that RVM presented well the relationship between geomechanics parameters and borehole pressure, and the proposed method can characterized the geomechanics parameters reasonably. Further, the parameters of hybrid back analysis were analyzed and discussed. It showed that the hybrid back analysis is feasible, effective, robust and has a good global searching performance. The proposed method provides a significant way to identify geomechanics parameters from hydraulic fracturing.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.