• Korean Journal of Soil Science and Fertilizer
• /
• v.34 no.3
• /
• pp.213-236
• /
• 2001

The ionic composition of soil solution is related to a nutrient uptake by plant. Many models for estimating ionic composition of solution have been developed, and most of them have been used for calculating a content of mineral and ionic species in a geochemical point of view. An approximation model considering both cation and anion in soil solution was developed. Variables such as pH, Eh, EC, cations(K, Ca, Mg. Na, Fe, Mn, Al, $NH_4{^+}$), anions(Si, S, P, CY, $NO_3{^-}$, $HCO_3{^-}$ and chemical equilibria of ionic species in soil solution were input into Excel sheet. The activities of soluble ion, ionpairs and complexes of input element were estimated by Newton-Raphson method using conditional equilibrium constant calculated by Davies equation and special models. Equilibrium contents of insoluble minerals and complexes were also calculated.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.2 no.2
• /
• pp.87-94
• /
• 1998

Behavioral characteristics of shallow circular arches with dynamic loading and different end conditions are analysed. Geometric nonlinearity is modelled using Lagrangian description of the motion. The finite element analysis procedure is used to solve the dynamic equation of motion, and the Newmark method is adopted in the approximation of time integration. The behavior of arches is analysed using the buckling criterion and non-dimensional time, load and shape parameters which Humphreys suggested. But a new deflection-ratio formula including the effect of horizontal displacement plus vertical displacement is presented to apply for the non-symmetric buckling problems. Through the model analysis, it's confirmed that fix-ended arches have higher buckling stability than hinge-ended arches, and arches with the same shape parameter have the same deflection ratio at the same time parameter when loaded with the same parametric load.

• Journal of rehabilitation welfare engineering & assistive technology
• /
• v.7 no.1
• /
• pp.21-27
• /
• 2013

Manual wheelchair propulsion can lead to pain and injuries of users due to mechanical inefficiency of wheelchair propulsion motion. The kinetic analysis of the upper limbs during manual wheelchair propulsion needs to be studied. A two dimensional inverse dynamic model of upper limbs was developed to compute the joint torque during manual wheelchair propulsion. The model was composed of three segments corresponding to upper arm, lower arm and hand. These segments connected in series by revolute joints constitute open chain mechanism in sagittal plane. The inverse dynamic method is based on Newton-Euler formalism. The model was applied to data collected in experiments. Kinematic data of upper limbs during wheelchair propulsion were obtained from three dimensional trajectories of markers collected by a motion capture system. Kinetic data as external forces applied on the hand were obtained from a dynamometer. The joint rotation angles and joint torques were computed using the inverse dynamic model. The developed model is for upper limbs biomechanics and can easily be extended to three dimensional dynamic model.

• Tribology and Lubricants
• /
• v.30 no.6
• /
• pp.337-342
• /
• 2014

The reduction of drag torque is an important research issue in terms of improving transmission efficiency. Drag torque in a wet clutch occurs because of the viscous drag generated by the transmission fluid in a narrow gap (clearance) between the friction plate and a separate plate. The objective of this paper is to observe the effects of the friction plate area and the clearance on the drag torque using finite element simulation. The two-phase flow of air and oil fluid is considered and modeled for the simulation. The simulation analysis reveals that as the rotational speed increases, the drag torque generally increases to a critical point and then decreases sharply at a high speed regime. The clearance between the two plates plays an important role in controlling drag torque peak. An increase in the clearance causes a decrease in shear stress; thus, the drag torque also decreases according to Newton's law of viscosity. An observation of the effect of the area of contact between transmission fluid and friction plate shows that the drag torque increases with the contact area. The flow vectors inside the flow channel present clear evidence that the velocity of the fluid flows is faster with a larger friction plate, that is, in the case of a larger contact area. Therefore, the optimum size of the friction plate should be determined carefully, considering both the clutch performance and drag reduction. It is expected that the results from this study can be very useful as a database for clutch design and to predict the drag torque for the initial design with respect to various clutch parameters.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.661-669
• /
• 2007

Eyer since the Prandtl's experiment in 1934 and X-21 airjet test in 1950 both attempting to reduce drag, it was found that controlling the velocities of surface for extremely fast-moving object in the air through suction or injection was highly effective and active method. To obtain the right amount of suction or injection, however, repetitive trial-and error parameter test has been still used up to now. This study started from an attempt to decide optimal amount of suction and injection of incompressible Navier-Stokes by employing optimization techniques. However, optimization with traditional methods are very limited, especially when Reynolds number gets high and many unexpected variables emerges. In earlier study, we have proposed an algorithm to solve this problem by using step by step method in analysis and introducing SQP method in optimization. In this study, we propose more effective and robust algorithm and techniques in solving flow optimization problem.

• Archives of design research
• /
• v.15 no.3
• /
• pp.5-14
• /
• 2002

The amusement design of a large amusement park, including theme park, largely consists of attraction design, the main facilities, and auxiliary facilities design. The purpose of this study was to discuss the physical motion characteristics of major ride, a recently highlighted representative ride attraction, modern people's thirst for stronger stimuli, and some of the right future directions for amusement design. In theme park, there is an growing trend for rides, especially major rides, to be more thrilling, and Newton's Law motion is basically applied. The rides we can often see provides a thrilling experience that is different from that offered by vehicles we use everyday, due to a variety of speed change, namely of acceleration. Every entertainment facilities, including ride attractions, are externally cutting edge with the use of new technology. but internally, poor theme or concept not only makes it difficult to function properly as attraction but poses a severe threat to theme park survival itself. Theme park will be able to exist as a comprehensive entertainment entity, when attractions are designed to have a perfect harmony among story, method and environment with cultural universality.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.3
• /
• pp.263-272
• /
• 2005

For the Communication, Ocean and Meteorological Satellite (COMS) which will be launched in 2008, an algorithm for finding the precise location of the sun-glint point on the ocean surface is studied. The precise locations of the sun-glint are estimated by considering azimuth and elevation angles of Sun-satellite-Earth geometric position and the law of reflection. The obtained nonlinear equations are solved by using the Newton-Raphson method. As a result, when COMS is located at $116.2^{\circ}E$ or $128.2^{\circ}E$ longitude, the sun-glint covers region of ${\pm}10^{\circ}(N-S)$ latitude and $80-150^{\circ}(E-W)$ longitude. The diurnal path of the sun-glint in the southern hemisphere is curved towards the North Pole, and the path in the northern hemisphere is forwards the south pole. The algorithm presented in this paper can be applied to predict the precise location of sun-glint region in any other geostationary satellites.

• Journal of Biosystems Engineering
• /
• v.38 no.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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.