• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.12
• /
• pp.204-213
• /
• 1999

This paper suggests a numerical method of finding optimal geometric path and minimum-time motion for a manipulator arm. To find the minimum-time motion, the optimal geometric path is searched first, and the minimum-time motion is searched on this optimal path. In the algorithm finding optimal geometric path, the objective function is minimizing the combination of joint velocities, joint-jerks, and actuator forces as well as avoiding several static obstacles, where global search is performed by adjusting the seed points of the obstacle models. In the minimum-time algorithm, the traveling time is expressed by the linear combinations of finite-term quintic B-splines and the coefficients of the splines are obtained by nonlinear programming to minimize the total traveling time subject to the constraints of the velocity-dependent actuator forces. These two search algorithms are basically similar and their convergences are quite stable.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.355-372
• /
• 2011

Based on the fully nonlinear velocity potential theory, the liquid sloshing in a three dimensional tank under random excitation is studied. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing scheme, B-spline curve, is applied to both the longitudinal and transverse directions of the tank to eliminate the possible saw-tooth instabilities. When the tank is undergoing one dimensional regular motion of small amplitude, the calculated results are found to be in very good agreement with linear analytical solution. In the simulation, the normal standing waves, travelling waves and bores are observed. The extensive calculation has been made for the tank undergoing specified random oscillation. The nonlinear effect of random sloshing wave is studied and the effect of peak frequency used for the generation of random oscillation is investigated. It is found that, even as the peak value of spectrum for oscillation becomes smaller, the maximum wave elevation on the side wall becomes bigger when the peak frequency is closer to the natural frequency.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.13 no.4
• /
• pp.58-65
• /
• 2005

A method measuring angular speed and estimating angular acceleration of an automotive wind shield wiper pivot with limited resources has been proposed. Limited resources refer to the fact that processes cannot be operated in real-time with a regular notebook running a Microsoft Windows. Also, they refer to the fact that data acquisition cards have only two general purpose counters as many generic cards do. An optical incremental encoder has been employed for measuring angular motion. To measure the angular speed of the pivot, periods for the encoder's output pulses have been measured as the speed is related to the reciprocal of the period. Since only information acquired from one counter channel is the magnitude of the angular speed, sign correction is necessary. Also the information for the exact time when a pivot passes left and right dead points is also missing and the situation is inherent to the hardware setup. To find out the zero-crossing time of the angular speed, a linear interpolation technique has been employed. Lastly, to overcome the imperfection of the mechanical encoders, the angular speed has been curve fitted to a spline. Angular acceleration can be obtained by a differentiation of the angular speed.

• Current Optics and Photonics
• /
• v.7 no.5
• /
• pp.518-528
• /
• 2023

As an effective means of remotely detecting the spectral information of the object, the spectral calibration for the Savart polarization interference imaging spectrometer (SPIIS) is a basis and prerequisite of information quantification, and its experimental calibration scheme is firstly proposed in this paper. In order to evaluate the accuracy of the spectral information acquisition, the linear interpolation, cubic spline interpolation, and piecewise cubic interpolation algorithms are adopted, and the precision of the quadratic polynomial fitting is the highest, whose fitting error is better than 5.8642 nm in the wavelength range of [500 nm, 820 nm]. Besides, the inversed value of the spectral resolution for the monochromatic light is greater than the theoretical value, and the deviation between them becomes larger with the wavelength increasing, which is mainly caused by the structural design of the SPIIS, together with the rationality of the spectral restoration algorithm and the selection of the maximum optical path difference (OPD). This work demonstrates that the SPIIS has achieved high performance assuring the feasibility of its practical use in various fields.

• Asian Pacific Journal of Cancer Prevention
• /
• v.13 no.5
• /
• pp.1851-1856
• /
• 2012

Immigration to a nation with a stronger anti-smoking environment has been hypothesized to make smoking less common. However, little is known about how environments influence risk of smoking across the lifecourse. Research suggested a linear decline in smoking over the lifecourse but these associations, in fact, might not be linear. This study assessed the possible nonlinear associations between age and smoking and examined how these associations differed by environment through comparing Koreans in Seoul, South Korea and Korean Americans in California, United States. Data were drawn from population based telephone surveys of Korean adults in Seoul (N=500) and California (N=2,830) from 2001-2002. Locally weighted scatterplot smoothing (lowess) was used to approximate the association between age and smoking with multivariable spline logistic regressions, including adjustment for confounds used to draw population inferences. Smoking differed across the lifecourse between Korean and Korean American men. The association between age and smoking peaked around 35 years among Korean and Korean American men. From 18 to 35 the probability of smoking was 57% higher (95%CI, 40 to 71) among Korean men versus 8% (95%CI, 3 to 19) higher among Korean American men. A similar difference in age after 35, from 40 to 57 years of age, was associated with a 2% (95%CI, 0 to 10) and 20% (95%CI, 16 to 25) lower probability of smoking among Korean and Korean American men. A nonlinear pattern was also observed among Korean American women. Social role transitions provide plausible explanations for the decline in smoking after 35. Investigators should be mindful of nonlinearities in age when attempting to understand tobacco use.

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.

• Steel and Composite Structures
• /
• v.37 no.6
• /
• pp.663-678
• /
• 2020

This paper proposes a novel topology optimization method generating multiple materials for external linear plane crack structures based on the combination of IsoGeometric Analysis (IGA) and eXtended Finite Element Method (X-FEM). A so-called eXtended IsoGeometric Analysis (X-IGA) is derived for a mechanical description of a strong discontinuity state's continuous boundaries through the inherited special properties of X-FEM. In X-IGA, control points and patches play the same role with nodes and sub-domains in the finite element method. While being similar to X-FEM, enrichment functions are added to finite element approximation without any mesh generation. The geometry of structures based on basic functions of Non-Uniform Rational B-Splines (NURBS) provides accurate and reliable results. Moreover, the basis function to define the geometry becomes a systematic p-refinement to control the field approximation order without altering the geometry or its parameterization. The accuracy of analytical solutions of X-IGA for the crack problem, which is superior to a conventional X-FEM, guarantees the reliability of the optimal multi-material retrofitting against external cracks through using topology optimization. Topology optimization is applied to the minimal compliance design of two-dimensional plane linear cracked structures retrofitted by multiple distinct materials to prevent the propagation of the present crack pattern. The alternating active-phase algorithm with optimality criteria-based algorithms is employed to update design variables of element densities. Numerical results under different lengths, positions, and angles of given cracks verify the proposed method's efficiency and feasibility in using X-IGA compared to a conventional X-FEM.

• The Korean Journal of Applied Statistics
• /
• v.32 no.4
• /
• pp.607-618
• /
• 2019

A new and interesting task in statistics is to effectively analyze functional data that frequently comes from advances in modern science and technology in areas such as meteorology and biomedical sciences. Functional linear regression with scalar response is a popular functional data analysis technique and it is often a common problem to determine a functional association if a functional predictor variable affects the scalar response in the models. Recently, Kong et al. (Journal of Nonparametric Statistics, 28, 813-838, 2016) established classical testing methods for this based on functional principal component analysis (of the functional predictor), that is, the resulting eigenfunctions (as a basis). However, the eigenbasis functions are not generally suitable for regression purpose because they are only concerned with the variability of the functional predictor, not the functional association of interest in testing problems. Additionally, eigenfunctions are to be estimated from data so that estimation errors might be involved in the performance of testing procedures. To circumvent these issues, we propose a testing method based on fixed basis such as B-splines and show that it works well via simulations. It is also illustrated via simulated and real data examples that the proposed testing method provides more effective and intuitive results due to the localization properties of B-splines.

• Journal of Preventive Medicine and Public Health
• /
• v.56 no.3
• /
• pp.282-290
• /
• 2023

Objectives: This study aimed to evaluate the potential interaction between kidney function and the non-linear association between serum calcium levels and cardiovascular disease (CVD) mortality. Methods: This study included 8927 participants enrolled in the Dong-gu Study. Albumin-corrected calcium levels were used and categorized into 6 percentile categories: <2.5th, 2.5-25.0th, 25.0-50.0th, 50.0-75.0th, 75.0-97.5th, and >97.5th. Restricted cubic spline analysis was used to examine the non-linear association between calcium levels and CVD mortality. Cox proportional hazard regression was used to estimate hazard ratios (HRs) for CVD mortality according to serum calcium categories. All survival analyses were stratified by the estimated glomerular filtration rate. Results: Over a follow-up period of 11.9±2.8 years, 1757 participants died, of whom 219 died from CVD. A U-shaped association between serum calcium and CVD mortality was found, and the association was more evident in the low kidney function group. Compared to the 25.0-50.0th percentile group for serum calcium levels, both low and high serum calcium tended to be associated with CVD mortality (<2.5th: HR, 6.23; 95% confidence interval [CI], 1.16 to 33.56; >97.5th: HR, 2.56; 95% CI, 0.76 to 8.66) in the low kidney function group. In the normal kidney function group, a similar association was found between serum calcium levels and CVD mortality (<2.5th: HR, 1.37; 95% CI, 0.58 to 3.27; >97.5th: HR, 1.65; 95% CI, 0.70 to 3.93). Conclusions: We found a non-linear association between serum calcium levels and CVD mortality, suggesting that calcium dyshomeostasis may contribute to CVD mortality, and kidney function may modify the association.