• 대한기계학회논문집A
• /
• 제26권5호
• /
• pp.939-951
• /
• 2002

We propose a new multiscale Galerkin method based on interpolation wavelets for two-dimensional Poisson's and plane elasticity problems. The major contributions of the present work are: 1) full multiresolution numerical analysis is carried out, 2) general boundaries are handled by a fictitious domain method without using a penalty term or the Lagrange multiplier, 3) no special integration rule is necessary unlike in the (bi-)orthogonal wavelet-based methods, and 4) an efficient adaptive scheme is easy to incorporate. Several benchmark-type problems are considered to show the effectiveness and the potentials of the present approach. is 1-2m/s and impact deformation of the electrode depends on the strain rate at that velocity, the dynamic behavior of the sinter-forged Cu-Cr is a key to investigate the impact characteristics of the electrodes. The dynamic response of the material at the high strain rate is obtained from the split Hopkinson pressure bar test using disc-type specimens. Experimental results from both quasi-static and dynamic compressive tests are Interpolated to construct the Johnson-Cook model as the constitutive relation that should be applied to simulation of the dynamic behavior of the electrodes. The impact characteristics of a vacuum interrupter are investigated with computer simulations by changing the value of five parameters such as the initial velocity of a movable electrode, the added mass of a movable electrode, the wipe spring constant, initial offset of a wipe spring and the virtual fixed spring constant.

• Wind and Structures
• /
• 제27권4호
• /
• pp.235-245
• /
• 2018

In this study the stress analysis of orthotropic thin plate with arbitrary shapes for different boundary conditionsis investigated. Meshfreemethod is applied to static analysis of thin plates with various geometries based on the Kirchhoff classical plate theory. According to the meshfree method the domain of the plates are expressed through a set of nodes without using mesh. In this method, a set of nodes are defined in a standard rectangular domain, then via a third order map, these nodes are transferred to the main domain of the original geometry; therefore the analysis of the plates can be done. Herein, Meshless local Petrov-Galerkin (MLPG) as a meshfree numerical method is utilized. The MLS function in MLPG does not satisfy essential boundary conditions using Delta Kronecker. In the MLPG method, direct interpolation of the boundary conditions can be applied due to constructing node by node of the system equations. The detailed parametric study is conducted, focusing on the arbitrary geometries of the thin plates. Results show that the meshfree method provides better accuracy rather than finite element method. Also, it is found that trend of the figures have good agreement with relevant published papers.

• 대한기계학회논문집A
• /
• 제24권1호
• /
• pp.190-195
• /
• 2000

It has been established that a crack has an important effect on the dynamic behavior of a structure. This effect depends mainly on the location and depth of the crack. To identify the location and depth of a crack in a structure. Nikolakopoulos et. al. used the intersection point of the superposed contours that correspond to the eigenfrequency caused by the crack presence. However the intersecting point of the superposed contours is not only difficult to find but also incorrect to calculate. A method is presented in this paper which uses optimization technique for the location and depth of the crack. The basic idea is to find parameters which use the structural eigenfrequencies on crack depth and location and optimization algorithm. With finite element model of the structure to calculate eigenfrequencies, it is possible to formulate the inverse problem in optimization format. Method of optimization is augmented lagrange multiplier method and search direction method is BFGS variable metric method and one dimensional search method is polynomial interpolation.

• 대한기계학회논문집
• /
• 제16권6호
• /
• pp.1163-1170
• /
• 1992

본 연구에서는 자동차 전자제어식 공기유량계를 Fig.2와 같이 평행평판 안에 사각 지주가 있다고 단순화하고 공기는 이차원 비 압축성 점성유동으로 가정했다. 지배방정식은 유체 운동량방정식(navier-Stokes equation)을 와도 전달 방정식(vorti- city transport equation)과 유량 함수 방정식(stream function equation)으로 변환하 여 사용하였다. Peacemanrachford ADI 방법으로 수치해석 하였으며, 유량 함수 방정 식의 수렴성을 좋게 하기 위하여 Wachspress parameter를 사용하였다. 벽면의 경계 조건은 Briely의 4th-order Lagrange interpolation 방법을 따랐다. Reynolds 수 200과 500에서의 비정상유동(unsteady flow)을 계산하였으며, 유동이 정상상태(steady state)에 도달하였을 때에 유동을 교란시켜 와동 흘림(vortex shedding)을 구하였다.

• International Journal of CAD/CAM
• /
• 제8권1호
• /
• pp.1-10
• /
• 2009

This paper presents a topology optimization approach using element-free Galerkin method (EFGM) for the optimal design of compliant mechanisms with geometrically non-linearity. Meshless method has an advantage over the finite element method(FEM) because it is more capable of handling large deformation resulted from geometrical nonlinearity. Therefore, in this paper, EFGM is employed to discretize the governing equations and the bulk density field. The sensitivity analysis of the optimization problem is performed by incorporating the adjoint approach with the meshless method. The Lagrange multipliers method adjusted for imposition of both the concentrated and continuous essential boundary conditions in the EFGM is proposed in details. The optimization mathematical formulation is developed to convert the multi-criteria problem to an equivalent single-objective problem. The popularly applied interpolation scheme, solid isotropic material with penalization (SIMP), is used to indicate the dependence of material property upon on pseudo densities discretized to the integration points. A well studied numerical example has been applied to demonstrate the proposed approach works very well and the non-linear EFGM can obtain the better topologies than the linear EFGM to design large-displacement compliant mechanisms.

• Journal of Positioning, Navigation, and Timing
• /
• 제3권4호
• /
• pp.155-161
• /
• 2014

The purpose of this study is to develop precise point positioning (PPP) algorithms based on GLONASS code-pseudorange, verify their performance and present their utility. As the basic correction models of PPP, we applied Inter Frequency Bias (IFB), relativistic effect, satellite antenna phase center offset, and satellite orbit and satellite clock errors, ionospheric errors, and tropospheric errors that must be provided on a real-time basis. The satellite orbit and satellite clock errors provided by Information-Analytical Centre (IAC) are interpolated at each observation epoch by applying the Lagrange polynomial method and linear interpolation method. We applied Global Ionosphere Maps (GIM) provided by International GNSS Service (IGS) for ionospheric errors, and increased the positioning accuracy by applying the true value calculated with GIPSY for tropospheric errors. As a result of testing the developed GLONASS PPP algorithms for four days, the horizontal error was approximately 1.4 ~ 1.5 m and the vertical error was approximately 2.5 ~ 2.8 m, showing that the accuracy is similar to that of GPS PPP.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Journal of Biosystems Engineering
• /
• 제38권2호
• /
• pp.138-148
• /
• 2013

Purpose: Determining an appropriate path is a top priority in order for a robot to maneuver in a dynamically efficient way especially in a pick-and-place task. In a non-standardized work environment, current robot arm executes its motion based on the kinematic displacements of joint variables, though resulting motion is not dynamically optimal. In this research we suggest analyzing and applying motion patterns of the human arm as an alternative to perform near optimum motion trajectory for arbitrary pick-and-place tasks. Methods: Since the motion of a human arm is very complicated and diverse, it was simplified into two links: one from the shoulder to the elbow, and the other from the elbow to the hand. Motion patterns were then divided into horizontal and vertical components and further analyzed using kinematic and dynamic methods. The kinematic analysis was performed based on the D-H parameters and the dynamic analysis was carried out to calculate various parameters such as velocity, acceleration, torque, and energy using the Newton-Euler equation of motion and Lagrange's equation. In an attempt to assess the efficacy of the analyzed human motion pattern it was compared to the virtual motion pattern created by the joint interpolation method. Results: To demonstrate the efficacy of the human arm motion mechanical and dynamical analyses were performed, followed by the comparison with the virtual robot motion path that was created by the joint interpolation method. Consequently, the human arm was observed to be in motion while the elbow was bent. In return this contributed to the increase of the manipulability and decrease of gravity and torque being exerted on the elbow. In addition, the energy required for the motion decreased. Such phenomenon was more apparent under vertical motion than horizontal motion patterns, and in shorter paths than in longer ones. Thus, one can minimize the abrasion of joints by lowering the stress applied to the bones, muscles, and joints. From the perspectives of energy and durability, the robot arm will be able to utilize its motor most effectively by adopting the motion pattern of human arm. Conclusions: By applying the motion pattern of human arm to the robot arm motion, increase in efficiency and durability is expected, which will eventually produce robots capable of moving in an energy-efficient manner.

• Nuclear Engineering and Technology
• /
• 제53권7호
• /
• pp.2104-2125
• /
• 2021

This paper presents the verification and validation (V&V) of a calculation module for isotope inventory prediction to control the back-end cycle of spent nuclear fuel (SNF). The calculation method presented herein was implemented in a two-step code system of a lattice code STREAM and a nodal diffusion code RAST-K. STREAM generates a cross section and provides the number density information using branch/history depletion branch calculations, whereas RAST-K supplies the power history and three history indices (boron concentration, moderator temperature, and fuel temperature). As its primary feature, this method can directly consider three-dimensional core simulation conditions using history indices of the operating conditions. Therefore, this method reduces the computation time by avoiding a recalculation of the fuel depletion. The module for isotope inventory calculates the number densities using the Lagrange interpolation method and power history correction factors, which are applied to correct the effects of the decay and fission products generated at different power levels. To assess the reliability of the developed code system for back-end cycle analysis, validation study was performed with 58 measured samples of pressurized water reactor (PWR) SNF, and code-to-code comparison was conducted with STREAM-SNF, HELIOS-1.6 and SCALE 5.1. The V&V results presented that the developed code system can provide reasonable results with comparable confidence intervals. As a result, this paper successfully demonstrates that the isotope inventory prediction code system can be used for spent nuclear fuel analysis.

• Nuclear Engineering and Technology
• /
• 제53권1호
• /
• pp.44-60
• /
• 2021

In this paper, we validate the decay heat calculation capability via a two-step method to analyze spent nuclear fuel (SNF) discharged from pressurized water reactors (PWRs). The calculation method is implemented with a lattice code STREAM and a nodal diffusion code RAST-K. One of the features of this method is the direct consideration of three-dimensional (3D) core simulation conditions with the advantage of a short simulation time. Other features include the prediction of the isotope inventory by Lagrange non-linear interpolation and the use of power history correction factors. The validation is performed with 58 decay heat measurements of 48 fuel assemblies (FAs) discharged from five PWRs operated in Sweden and the United States. These realistic benchmarks cover the discharge burnup range up to 51 GWd/MTU, 23.2 years of cooling time, and spanning an initial uranium enrichment range of 2.100-4.005 wt percent. The SNF analysis capability of STREAM is also employed in the code-to-code comparison. Compared to the measurements, the validation results of the FA calculation with RAST-K are within ±4%, and the pin-wise results are within ±4.3%. This paper successfully demonstrates that the developed decay heat calculation method can perform SNF back-end cycle analyses.