• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.216-232
• /
• 1991

This Finite Element Program(TAFEM) has been developed to be able to carry out the structural analsis of tunnel section and simulate the surrounding ground behaviour due to New Austrian Tunnelling Method, of which main support is the surrounding ground, itself. The Elasto-plastic theory has been applied. The used finite elements are 8-noded isoparametric element(rock & shotcrete), 2 or 3-noded rod element(rock bolt) and infinite boundary element. The load incremental method and tangential stiffness method has been used. Associated flow rule was applied to plastic flow and yield criteria inclued not only Mohr-Coulomb but also Drucker-Prager. In this paper, Drucker-Prager yield criterion has been used. The relationship between plastic strain and stress is based on the incremental strain concept and stress-strain equation on the basis of the stress path of each gauss point has been adopted. It may be rational that rock is considered to be no-tension material, so that no-tension analysis has been adopted in accordance with the brittle fracture constitutive equation.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1998.04a
• /
• pp.167-172
• /
• 1998

An algorithm of Thermal-elastohydrodynamic lubrication analysis for the piston ring is developed. This algorithm contains cavitation boundary condition so it automatically satisfies conservation of mass. 1-D Reynolds equation and 2-D energy equation are solved simultaneously by using Gauss-Jordan method and Newton-Raphson method. Minimum film thickness and friction force are calculated for 1 cycle. There is little difference between the results caculated by isothermal rigid and EHL analysis in entire cycle. In the results of THL, shear heating effect and temperature boundary condition affect the minimum film thickness and friction force prediction. The minimum film thickness and the friction force calculated by THL are lower than those caculated using isothermal assumption.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.957-962
• /
• 2004

This paper deals with the free vibration analysis of beam-columns resting on Pasternak foundation by the Differential Quadrature Method. Based on the differential equation subjected to the boundary conditions, adopted from the open literature, which governs the free vibrations of such member, this equation is applied to the Differential Quadrature Method. For computing natural frequencies, the numerical procedures are developed by QR Algorithm, in which the Chebyshev-Gauss-Lobatto method is used for choosing the grid points. The numerical methods developed herein for computing natural frequencies are programmed in FORTRAN code, and all solutions obtained in this study are quite agreed with those in the open literature.

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.101-103
• /
• 1999

This paper presents an analysis of 2-dimensional(2-D) multi-regions problem using boundary integral equation method(BIEM). When compared with finite element method(FEM), there are only a few unknown variables in BIEM because it implements numerical analysis only for the surface or boundary of a model. As a result, a lot of computational memory and time can be saved. Procedure to analyze 2-D multi-regions problem using potentials and its derivatives in a boundary as unknown variables, first, numerical analysis is performed for each of subregions. And then interface continuity condition is applied to the interface between them and Gauss Quadrature Formula are adopted to solve singular integral in a boundary in this paper.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2006.05a
• /
• pp.163-164
• /
• 2006

In this paper, an analytical laser ablation model with Maxwell equation will be addressed by considering relationship between laser ablation and material. The Maxwell equation consists of four equations: two Gauss laws for electric and magnetic fields, Faraday's law, and Ampere's law. This analytical model will be calculated by employing Finite Difference Time Domain (FDTD). This method also makes it possible to simulate the laser beam propagation in a wide range of materials, such as metals, semiconductors, and dielectrics. Therefore, in this study, a numerical model for short pulse laser interaction with materials is developed, focusing on the accurate description of laser beam propagation and ablation process into the material with each pulse.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.2
• /
• pp.487-494
• /
• 2016

Let p be an odd prime and c be a fixed integer with (c, p) = 1. For each integer a with $1{\leq}a{\leq}p-1$, it is clear that there exists one and only one b with $0{\leq}b{\leq}p-1$ such that $ab{\equiv}c$ mod p. Let N(c, p) denote the number of all solutions of the congruence equation $ab{\equiv}c$ mod p for $1{\leq}a$, $b{{\leq}}p-1$ in which a and $\bar{b}$ are of opposite parity, where $\bar{b}$ is defined by the congruence equation $b{\bar{b}}{\equiv}1$ mod p. The main purpose of this paper is using the mean value theorem of Dirichlet L-functions and the properties of Gauss sums to study the computational problem of one kind mean value function related to $E(c,p)=N(c,p)-{\frac{1}{2}}{\phi}(p)$, and give its an exact computational formula.

• Structural Engineering and Mechanics
• /
• v.91 no.5
• /
• pp.459-468
• /
• 2024

This study presents a frictionless receding contact problem for an orthotropic elastic layer. It is assumed that the layer is supported by two rigid quarter planes and the material of the layer is orthotropic. The layer of thickness h is indented by a rigid cylindrical punch of radius R. The problem is modeled by using the singular integral equation method with the help of the Fourier transform technique. Applying the boundary conditions of the problem the system of singular integral equations is obtained. In this system, the unknowns are the contact stresses and contact widths under the punch and between the layer and rigid quarter planes. The Gauss-Chebyshev integration method is applied to the obtained system of singular integral equations of Cauchy type. Five different orthotropic materials are considered during the analysis. Numerical results are presented to interpret the effect of the material property and the other parameters on the contact stress and the contact width.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1072-1078
• /
• 2002

The constrained motion requires the determination of constraint force acting on unconstrained systems for satisfying given constraints. Most of the methods to decide the force depend on numerical approaches such that the Lagrange multiplier method, and the other methods need vector analysis or complicated intermediate process. In 1992, Udwadia and Kalaba presented the generalized inverse method to describe the constrained motion as well as to calculate the constraint force. The generalized inverse method has the advantages which do not require any linearization process for the control of nonlinear systems and can explicitly describe the motion of holonomically and/or nongolonomically constrained systems. In this paper, an explicit equation to describe the constrained motion is derived by minimizing the performance index, which is a function of constraint force vector, with respect to the constraint force. At this time, it is shown that the positive-definite weighting matrix in the performance index must be the inverse of mass matrix on the basis of the Gauss's principle and the derived differential equation coincides with the generalized inverse method. The effectiveness of this method is illustrated by means of two numerical applications.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.2410-2413
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAB, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. In the pre-processor of 3D FEM package, a new 3D mesh generating algorithm can make information on 3D Delaunay tetrahedral mesh elements for analyses of 3D FEM problems. The solver of the 3D FEM package offers three methods for solving the linear algebraic matrix equation, i.e., Gauss-Jordan elimination solver, Band solver, and Skyline solver. The post-processor visualizes the results for 3D FEM problems such as the deformed position and the stress. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.48 no.9
• /
• pp.1081-1087
• /
• 1999

Load Flow calulation methods can usually be divided into Gauss-Seidel method, Newton-Raphson method and decoupled method. Load flow calculation is a basic on-line or off-line process for power system planning. operation, control and state analysis. These days Newton-Raphson method is mainly used since it shows remarkable convergence characteristics. It, however, needs considerable calculation time in construction and calculation of inverse Jacobian matrix. In addition to that, Newton-Raphson method tends to fail to converge when system loading is heavy and system has a large R/X ratio. In this paper, matrix equation is used to make algebraic expression and then to slove load flow equation and to modify above defects. And it preserve P-Q bus part of Jacobian matrix to shorten computing time. Application of mentioned algorithm to 14 bus, 39 bus, 118 bus systems led to identical results and the same numbers of iteration obtained by Newton-Raphson method. The effect of computing time reduction showed about 28% , 30% , at each case of 39 bus, 118 bus system.