• Geomechanics and Engineering
• /
• v.30 no.4
• /
• pp.363-372
• /
• 2022

One of the most frequent issues in tunnel excavation is the collapse of rock blocks and the dropping of rock fragments from the tunnel face. The tunnel face can be reinforced using a number of techniques. One of the most popular and affordable solutions is the use of face longitudinal dowels, which has benefits including high strength, flexibility, and ease of cutting. In order to examine the reinforced face, this work shows the longitudinal deformation profile and ground response curve for a tunnel face. This approach is based on assumptions made during the analysis phase of problem solving. By knowing the tunnel face response and dowel behavior, the interaction of two elements can be solved. The rock element equation derived from the rock bolt method is combined with the dowel differential equation to solve the reinforced ground response curve (GRC). With a straightforward and accurate analytical equation, the new differential equation produces the reinforced displacement of the tunnel face at each stage of excavation. With simple equations and a less involved computational process, this approach offers quick and accurate solutions. The FLAC3D simulation has been compared with the suggested analytical approach. A logical error is apparent from the discrepancies between the two solutions. Each component of the equation's effect has also been described.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.4
• /
• pp.238-266
• /
• 2013

The Interaction Picture (IP) method is a valuable alternative to Split-step methods for solving certain types of partial differential equations such as the nonlinear Schr$\ddot{o}$dinger equation or the Gross-Pitaevskii equation. Although very similar to the Symmetric Split-step (SS) method in its inner computational structure, the IP method results from a change of unknown and therefore do not involve approximation such as the one resulting from the use of a splitting formula. In its standard form the IP method such as the SS method is used in conjunction with the classical 4th order Runge-Kutta (RK) scheme. However it appears to be relevant to look for RK scheme of higher order so as to improve the accuracy of the IP method. In this paper we investigate 5th order Embedded Runge-Kutta schemes suited to be used in conjunction with the IP method and designed to deliver a local error estimation for adaptive step size control.

• Earthquakes and Structures
• /
• v.10 no.2
• /
• pp.389-408
• /
• 2016

Recently, seismic hazard analysis has become a very significant issue. New systems and available data have been also developed that could help scientists to explain the earthquakes phenomena and its physics. Scientists have begun to accept the role of uncertainty in earthquake issues and seismic hazard analysis. However, handling the existing uncertainty is still an important problem and lack of data causes difficulties in precisely quantifying uncertainty. Ground Motion Prediction Equation (GMPE) values are usually obtained in a statistical method: regression analysis. Each of these GMPEs uses the preliminary data of the selected earthquake. In this paper, a new fuzzy method was proposed to select suitable GMPE at every intensity (earthquake magnitude) and distance (site distance to fault) according to preliminary data aggregation in their area using ${\alpha}$ cut. The results showed that the use of this method as a GMPE could make a significant difference in probabilistic seismic hazard analysis (PSHA) results instead of selecting one equation or using logic tree. Also, a practical example of this new method was described in Iran as one of the world's earthquake-prone areas.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.1
• /
• pp.19-30
• /
• 2019

In this paper, we develop an accurate explicit finite difference method for the two-dimensional Black-Scholes equation with a hybrid boundary condition. In general, the correlation term in multi-asset options is problematic in numerical treatments partially due to cross derivatives and numerical boundary conditions at the far field domain corners. In the proposed hybrid boundary condition, we use a linear boundary condition at the boundaries where at least one asset is zero. After updating the numerical solution by one time step, we reduce the computational domain so that we do not need boundary conditions. To demonstrate the accuracy and efficiency of the proposed algorithm, we calculate option prices and their Greeks for the two-asset European call and cash-or-nothing options. Computational results show that the proposed method is accurate and is very useful for nonlinear boundary conditions.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.10a
• /
• pp.147-158
• /
• 2001

In general, two methods have been used to predict settlement of soft ground. One method is Terzaghi's one dimensional consolidation theory which gives time-settlement relationship using the standard consolidation test results. The other is forecasting method of ground settlement to be occured in the future using in-situ monitoring data. The above both methods have some defects in application manner or in itself especially in very deep and soft clayey ground. In view of the lessons and experiences of soft ground improvement projects, several techniques were proposed for more accurate theorectical calculation of consolidation settlement as follows ; ① Subdivision of soft ground, ② Consideration of secondary compression, ③ Using the modified compression index, etc. And also, revised hyperbolic fitting method was suggested to minimize the error of predicted future settlement. In addition, revised De-Beer equation of immediate settlement of loose sandy soil was proposed to overcome the tendency to show too small settlement calculation results by original De-Deer equation. And also, considering the various effects of settlement delay in the improved ground by vertical drains, time-settlement caculation equation(Onoue method) was revised to match the tendency of settlement delay by using the characteristics of discharge capacity decreases of vertical drain with time elapse by the pattern of hyperbolic equation.

• Journal of information and communication convergence engineering
• /
• v.4 no.1
• /
• pp.28-35
• /
• 2006

In this paper, we propose to synchronization method for mutual cooperative control method that have unstable limit cycles in a chaos trajectory surface in the chaotic UAVs. We assume all obstacles in the chaos trajectory surface have a Van der Pol equation with an unstable limit cycle. We also show computer simulation results of Arnold equation, Chua's equation trajectories with one or more Van der Pol as a obstacles. We proposed and verified the results of the method to make the embedding chaotic UAV to synchronization with the chaotic trajectory in any plane.

• Journal of the Korean Mathematical Society
• /
• v.61 no.4
• /
• pp.621-657
• /
• 2024

In this study, the numerical solution of a space-fractional Burgers' equation with initial and boundary conditions is considered. This equation is the simplest nonlinear model for diffusive waves in fluid dynamics. It occurs in a variety of physical phenomena, including viscous sound waves, waves in fluid-filled viscous elastic pipes, magneto-hydrodynamic waves in a medium with finite electrical conductivity, and one-dimensional turbulence. The proposed QBS/CNG technique consists of the Galerkin method with a function basis of quadratic B-splines for the spatial discretization of the space-fractional Burgers' equation. This is then followed by the Crank-Nicolson approach for time-stepping. A linearized scheme is fully constructed to reduce computational costs. Stability analysis, error estimates, and convergence rates are studied. Finally, some test problems are used to confirm the theoretical results and the proposed method's effectiveness, with the results displayed in tables, 2D, and 3D graphs.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1284-1296
• /
• 2004

In order to investigate effects of anisotropic fiber packing on stresses in composites, a Volume Integral Equation Method is applied to calculate the elastostatic field in an unbounded isotropic elastic medium containing multiple orthotropic inclusions subject to remote loading, and a Mixed Volume and Boundary Integral Equation Method is introduced for the solution of elastostatic problems in unbounded isotropic materials containing multiple anisotropic inclusions as well as one void under uniform remote loading. A detailed analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out for square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively. Also, an analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out, when it is assumed that a void is replaced with one inclusion adjacent to the central inclusion of square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively, due to manufacturing and/or service induced defects. The effects of random orthotropic fiber packing on stresses at the interface between the isotropic matrix and the central orthotropic inclusion are compared with the influences of square and hexagonal orthotropic fiber packing on stresses. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with multiple orthotropic inclusions and one void, it will be established that these new methods are very accurate and effective for investigating effects of general anisotropic fiber packing on stresses in composites.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.709-714
• /
• 2000

A splitting method for the direct numerical simulation of solid-liquid mixtures is presented, where a symmetric pressure equation is newly proposed. Through numerical experiment, it is found that the newly proposed splitting method works well with a matrix-free formulation fer some bench mark problems avoiding an erroneous pressure field which appears when using the conventional pressure equation of a splitting method. When deriving a typical pressure equation of a splitting method, the motion of a solid particle has to be approximated by the 'intermediate velocity' instead of treating it as unknowns since it is necessary as a boundary condition. Therefore, the motion of a solid particle is treated in such an explicit way that a particle moves by the known form drag (pressure drag) that is calculated from the pressure equation in the previous step. From the numerical experiment, it was shown that this method gives an erroneous pressure field even for the very small time step size as a particle velocity increases. In this paper, coupling the unknowns of particle velocities in the pressure equation is proposed, where the resulting matrix is reduced to the symmetric one by applying the projector of the combined formulation. It has been tested over some bench mark problems and gives reasonable pressure fields.

• Journal of the Korean Mathematical Society
• /
• v.48 no.3
• /
• pp.487-497
• /
• 2011

In this paper we propose a simple iterative method for finding a root of a nonlinear equation. It is shown that the new method, which does not require any derivatives, has a quadratic convergence order. In addition, one can find that a hybrid method combined with the non-iterative method can further improve the convergence rate. To show the efficiency of the presented method we give some numerical examples.