• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.439-446
• /
• 2008

We have developed a method to evaluate pumping rates from a system of pumping-well family. For a given system actual pumping rates depend on pump characteristics and the sum of the static head and the dynamic head. The static head is the elevation difference between the natural groundwater level and the outlet of the pipeline that connects all the wells. Major components of the dynamic head are groundwater drawdown in the well and pipeline head loss. The dynamic head and the pump characteristics depend on the pumping rates. Actual pumping rates are determined at the intersections of the system total-head curves and the pump characteristic curves. The Newton-Raphson's method is used to solve the nonlinear simultaneous equations. The method is applied to a hypothetical well family. Impacts of various design and operational parameters on the pumping rates are analyzed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.3
• /
• pp.215-223
• /
• 2009

In this paper, the helicopter aerodynamics is simulated in hovering and forward flight. Also, an overlapped grid technique is applied in this simulation to consider the blade motion and moving effects. The Caradonna & Tung's rotor blade was selected to analyze the unsteady aerodynamics in hovering and non-lift forward flight. Also, the AH-1G rotor blade was selected in forward flight. In forward flight case, the numerical trim was applied to determine the cyclic pitching angles using Newton-Raphson method, and the numerical results were in good agreement with experimental data, especially, the BVI effects were well simulated in advancing side in comparison other numerical results. The governing equation is a three dimensional unsteady Euler equation, and the Riemann invariants condition is used for inflow and outflow at the boundary.

• Journal of Korean Society of Steel Construction
• /
• v.19 no.4
• /
• pp.345-355
• /
• 2007

During winter, the CWR (continuous welded rail) may be broken when a temperature drop below the neutral level changes the axial force, causing tensile fracture and creating a rail gap. The passage of a train on a rail with an open gap may lead to very costly derailments. In this paper, the use of a track-and-train-coupled model whose rail has an open gap is proposed for dynamic interaction analysis. Linear track and train systems were coupled in this study by a nonlinear Herzian contact spring, and the complete system matrices of the total track-train system were constructed. Moreover, the interaction phenomenon considering the presence of an open gap in the rail was toughly defined by assigning the irregularity functions between the two sides of the gap. Time history analysis, which has an iteration scheme such as the Newmark-$\beta$ method (based on the Modified Newton-Raphson methods), was conducted to solve the nonlinear equation. .Finally, numerical studies were conducted to assess the effect of the various parameters of the system when applied to various speeds, open-gap sizes, and support stiffnesses of the rail.

• Advances in nano research
• /
• v.7 no.6
• /
• pp.405-417
• /
• 2019

This research is devoted to study post-buckling analysis of functionally graded carbon nanotubes reinforced composite (FG-CNTRC) micro plate with cut out subjected to magnetic field and resting on elastic medium. The basic formulation of plate is based on first order shear deformation theory (FSDT) and the material properties of FG-CNTRCs are presumed to be changed through the thickness direction, and are assumed based on rule of mixture; moreover, nonlocal Eringen's theory is applied to consider the size-dependent effect. It is considered that the system is embedded in elastic medium and subjected to longitudinal magnetic field. Energy approach, domain decomposition and Rayleigh-Ritz methods in conjunction with Newton-Raphson iterative technique are employed to trace the post-buckling paths of FG-CNTRC micro cut out plate. The influence of some important parameters such as small scale effect, cut out dimension, different types of FG distributions of CNTs, volume fraction of CNTs, aspect ratio of plate, magnitude of magnetic field, elastic medium and biaxial load on the post-buckling behavior of system are calculated. With respect to results, it is concluded that the aspect ratio and length of square cut out have negative effect on post-buckling response of micro composite plate. Furthermore, existence of CNTs in system causes improvement in the post-buckling behavior of plate and different distributions of CNTs in plate have diverse response. Meanwhile, nonlocal parameter and biaxial compression load on the plate has negative effect on post-buckling response. In addition, imposing magnetic field increases the post-buckling load of the microstructure.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.326-326
• /
• 2017

개수로에서 비에너지(specific energy)는 수로바닥을 기준으로 단위무게의 물이 가지는 에너지로 정의되며 흐름의 위치수두와 속도수두의 합으로 표현된다. 비에너지는 수로단면의 변화에 따른 수심의 변화를 해석하기 위하여 사용되는 중요한 개념이다. 사각형 개수로에서의 비에너지 관계식은 3차방정식의 형태이며, 해석적으로 3개의 해(3개의 수심)를 가지나, 물리적인 의미를 가지는 해는 2개이며 나머지 하나의 해는 음수이므로 물리적인 의미를 가지지 않는다. 물리적인 의미를 가지는 2개의 해는 각각 흐름이 상류(subcritical flow)인 경우와 사류(supercritical flow)인 경우에 대한 수심이다. 즉, 일정한 유량이 흐르는 조건에서 동일한 비에너지를 가지는 수심이 상류와 사류에 각각 존재하는데, 이 2개의 수심을 대응수심(alternate depths)이라 정의한다. 이러한 사각형 개수로에 대한 비에너지 관계식은 3차방정식이므로 그 해석해를 구할 수 있어, 수로단면의 변화에 따른 흐름의 변화를 비교적 쉽게 해석할 수 있다. 사각형 개수로가 아닌 경우의 비에너지 관계식을 이론적으로 고찰하는 연구는 찾아보기 힘들다. 이에 본 연구에서는 포물선형 개수로에 대해서 비에너지 관계식을 유도하였다. 유도된 비에너지 관계식은 비선형 음함수의 형태로 해석적으로 해를 구할 수 없다. 유도된 관계식의 해법으로 2차의 정밀도를 가지는 Newton-Raphson방법을 이용하였으며, 계산의 초기치는 상용화된 Excel에서 쉽게 구할 수 있는 회귀식을 이용하여 구하였다. 적용 예를 통해, 단순 회귀식을 이용하는 경우에는 정해와의 상대오차가 2 - 8% 내외였는데, 본 연구에서 제안하는 방법을 사용하는 경우에는 동일한 조건에서 상대오차가 0.25% 내외를 보였다. 즉 본 연구에서 제시하고 있는 양해법을 이용하면, 포물선형 개수로 흐름의 대응수심을 용이하게 그리고 정확도가 매우 높게 산정할 수 있다.

• Journal of the Society of Disaster Information
• /
• v.16 no.4
• /
• pp.832-841
• /
• 2020

Purpose: This study attempts at analyzing mechanical response of functionally graded material (FGM) plates in bending. An accurate and effective numerical approach based on isogeometric analysis (IGA) combined with higher-order shear deformation plate theory to predict the nonlinear flexural behavior is developed. Method: A higher-order shear deformation theory(HSDT) which accounts for the geometric nonlinearity in the von Karman sense is presented and used to derive the equilibrium and governing equations for FGM plate in bending. The nonlinear equations are solved by the modified Newton-Raphson iterative technique. Result: The volume fraction, plate length-to-thickness ratio and boundary condition have signifiant effects on the nonlinear flexural behavior of FGM plates. Conclusion: The proposed IGA method can be used as an accurate and effective numerical tool for analyzing the mechanical responses of FGM plates in flexure.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.151-157
• /
• 2021

This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

• Journal of the Society of Disaster Information
• /
• v.17 no.4
• /
• pp.839-847
• /
• 2021

Purpose: This study investigates mechanical behavior of functionally graded (FG) carbon nanotube-reinforced composite (CNTRC) plate in flexure. Isogeometric analysis (IGA) method coupled with shear deformable theory of higher-order (HSDT) to analyze the nonlinear bending response is presented. Method: Shear deformable plate theory into which a polynomial shear shape function and the von Karman type geometric nonlinearity are incorporated is used to derive the nonlinear equations of equilibrium for FG-CNTRC plate in bending. The modified Newton-Raphson iteration is adopted to solve the system equations. Result: The dispersion pattern of carbon nanotubes, plate geometric parameter and boundary condition have significant effects on the nonlinear flexural behavior of FG-CNTRC plate. Conclusion: The proposed IGA method coupled with the HSDT can successfully predict the flexural behavior of FG-CNTRC plate.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2B
• /
• pp.199-209
• /
• 2010

A finite element model for analyzing surface water flow was developed. Shallow water equation was discretized and solved by Galerkin and Newton-Raphson method. Triangular or rectangular elements can be mixed together to construct meshes. The algebraic equation was solved by frontal method which is very efficient in finite element problem. The developed model was applied to rectangular meandering channel with two bends and transverse velocities and water depth distributions were examined. High velocity was located near the inner bank at the apexes of the bends and velocity distribution was symmetrical about the centerline at the midsection of two bend and super elevation also occurred. Simulation results showed very good agreement with measured data. Another numerical simulation was carried out in mild, steep, adverse and abrupt bottom change slope and channels with weir. 12 water surface profiles of gradually varied flow were correct in terms of hydraulic interpretation.

• The KIPS Transactions:PartA
• /
• v.12A no.5 s.95
• /
• pp.413-420
• /
• 2005

The Newton-Raphson iterative algorithm for finding a floating point reciprocal square mot calculates it by performing a fixed number of multiplications. In this paper, a variable latency Newton-Raphson's reciprocal square root algorithm is proposed that performs multiplications a variable number of times until the error becomes smaller than a given value. To find the rediprocal square root of a floating point number F, the algorithm repeats the following operations: '$X_{i+1}=\frac{{X_i}(3-e_r-{FX_i}^2)}{2}$, $i\in{0,1,2,{\ldots}n-1}$' with the initial value is '$X_0=\frac{1}{\sqrt{F}}{\pm}e_0$'. The bits to the right of p fractional bits in intermediate multiplication results are truncated and this truncation error is less than '$e_r=2^{-p}$'. The value of p is 28 for the single precision floating point, and 58 for the double precision floating point. Let '$X_i=\frac{1}{\sqrt{F}}{\pm}e_i$, there is '$X_{i+1}=\frac{1}{\sqrt{F}}-e_{i+1}$, where '$e_{i+1}{<}\frac{3{\sqrt{F}}{{e_i}^2}}{2}{\mp}\frac{{Fe_i}^3}{2}+2e_r$'. If '$|\frac{\sqrt{3-e_r-{FX_i}^2}}{2}-1|<2^{\frac{\sqrt{-p}{2}}}$' is true, '$e_{i+1}<8e_r$' is less than the smallest number which is representable by floating point number. So, $X_{i+1}$ is approximate to '$\frac{1}{\sqrt{F}}$. Since the number of multiplications performed by the proposed algorithm is dependent on the input values, the average number of multiplications Per an operation is derived from many reciprocal square root tables ($X_0=\frac{1}{\sqrt{F}}{\pm}e_0$) with varying sizes. The superiority of this algorithm is proved by comparing this average number with the fixed number of multiplications of the conventional algorithm. Since the proposed algorithm only performs the multiplications until the error gets smaller than a given value, it can be used to improve the performance of a reciprocal square root unit. Also, it can be used to construct optimized approximate reciprocal square root tables. The results of this paper can be applied to many areas that utilize floating point numbers, such as digital signal processing, computer graphics, multimedia, scientific computing, etc.