• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5A
• /
• pp.457-465
• /
• 2009

This study presents a moving least squares (MLS) finite difference method for solving elastic crack problems with stress singularity at the crack tip. Near-tip functions are intrinsically employed in the MLS approximation to model near-tip field inducing singularity in stress field. employment of the functions does not lose the merit of the MLS Taylor polynomial approximation which approximates the derivatives of a function without actual differentiating process. In the formulation of crack problem, computational efficiency is considerably improved by taking the strong formulation instead of weak formulation involving time consuming numerical quadrature Difference equations are constructed on the nodes distributed in computational domain. Numerical experiments for crack problems show that the intrinsically enriched MLS finite difference method can sharply capture the singular behavior of near-tip stress and accurately evaluate stress intensity factors.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.6
• /
• pp.527-536
• /
• 2013

In order to predict structural vibration and radiated noise of high-voltage transformer in operation, it is necessary to precisely find the excitation force generated by the coils and core. However, finding the excitation force through experiments of high voltage transformer in operation is not possible. Therefore, this paper deals with identifying the excitation force by using the acceleration data measured through experiments and the transfer function estimated through finite element model. A method to predict structural vibration and radiated noise was also proposed. Three-phase windings and the core are the source of high-voltage transformer. The excitation forces were identified using the acceleration data and the transfer function of the surface of the tank. Structural vibration and radiated noise from the surface of the tank was predicted by using the identified excitation force. As a result of the interpretation of the experimental and computational analysis of structural vibration from the surface of the tank and radiated noise from the field point, the interpretation of the computational analysis showed relatively good accordance with the experiment.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.8
• /
• pp.78-86
• /
• 2000

This study identified a linear time-invariant mathematical model of each stand of a five-stand tandem cold mill. Two model identification methods are applied to construct a linear model of each stand of the tandem cold mill. For the model identification the input-output data that have interstand interference property in tandem cold rolling are obtained from a nonlinear simulator of the tandem cold mill. And a linear model of each stand is identified with N4SD(numerical algorithms for subspace state space system identification) method based on a state-space model and Least Square algorithm based on a transfer function. Furthermore a modeling error of the tandem cold mill is quantitatively analyzed from a maximum singular value plot of error function between an identified nominal model and uncertain model. In conclusion the comparison of the output signals between the existing Taylor linearized model the identified linear model and the nonlinear model of the tandem cold mill shows the accuracy and the applicability of the proposed identified model.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.5
• /
• pp.300-308
• /
• 2018

The contribution analysis of various external aeroacoustic noise sources to interior noise is important, enabling to design an automobile with a low interior noise level. With a new technique, the CD (Cholesky Decomposition), it is proposed to decompose an overall interior noise spectrum into multiple spectra, each representing the contribution of a specific noise source to the interior noise. In order to validate this method, three kinds of experiments were conducted. Furthermore, it is proposed to improve the CD-based contribution analysis method to be integrated with existing exterior microphone arrays in the wind tunnel. This method was validated with an experiment with two speakers.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.9
• /
• pp.1643-1653
• /
• 1992

A method, which uses analytical or numerical modal analysis results, e.g. from finite element analysis, to select desirable response measurement and excitation points for an efficient modal testing is introduced. First, points of master degree of freedom(DOP) are determined so as to statistically minimize errors between responses of a full order model and those estimated from the reduced order model. Such master DOF's are selected as the response measurement points. Then a criterion named 'driving point model constant(DPMC)' related to the magnitudes of resonance peaks of the driving point freqency response functions used to select the point of excitation out of the master DOF's. In this work, the method is demonstrated through applications to modal testing on a one dimensional cantilever beam and an aluminum plate and the results are compared with those by another technique. also, the method is applied to a two dimensional structural component of a passenger car.

• Journal of Korean Society for Atmospheric Environment
• /
• v.27 no.6
• /
• pp.650-662
• /
• 2011

The statistical tools such as empirical orthogonal function (EOF), and singular value decomposition (SVD) have been applied to analyze the characteristic of air pollutant over southeast Asia as well as to evaluate Zimeke's tropospheric column ozone (ZTO) determined by tropospheric residual method. In this study, we found that the EOF and SVD analyses are useful methods to extract the most significant temporal and spatial pattern from enormous amounts of satellite data. The EOF analyses with OMI $NO_2$ and OMI HCHO over southeast Asia revealed that the spatial pattern showed high correlation with fire count (r=0.8) and the EOF analysis of CO (r=0.7). This suggests that biomass burning influences a major seasonal variability on $NO_2$ and HCHO over this region. The EOF analysis of ZTO has indicated that the location of maximum ZTO was considerably shifted westward from the location of maximum of fire count and maximum month of ZTO occurred a month later than maximum month (March) of $NO_2$, HCHO and CO. For further analyses, we have performed the SVD analyses between ZTO and ozone precursor to examine their correlation and to check temporal and spatial consistency between two variables. The spatial pattern of ZTO showed latitudinal gradient that could result from latitudinal gradient of stratospheric ozone and temporal maximum of ZTO in March appears to be associated with stratospheric ozone variability that shows maximum in March. These results suggest that there are some sources of error in the tropospheric residual method associated with cloud height error, low efficiency of tropospheric ozone, and low accuracy in lower stratospheric ozone.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.10
• /
• pp.3048-3057
• /
• 1996

Analysis model containing two inclined surface cracks on semi-infinite elastic body is established and analyzed on the basis of linear fracture mechanics to examine mutual interference of two surface cracks. Muskhelishvili's complex stress functions are introduced and a set of singular integral equations is obtained for a dislocation density function. The stress intensity factors at crack tip are obtained by using the Gerasoulis'method. When two surface cracks are parallel and have the same length, the values of $K_1$and $\Delta K_11$(variation of $K_11$) for crack 1 and crack 2 decrease by the mutual interference of two surface cracks as the distance between the two surface cracks shortens. The effect of mutual interference is remarkable in high friction coefficient. In case that two surface cracks are parallel, the values of $K_1$and $\Delta K_11$for crack 2 decrease as the length ratio ot crack 2 to crack 1 becomes small. As the crack inclination angle rises, the value of $K_1$ and the mutual interference of $K_1$for crack 2 increase and the value of$\Delta K_11$ for crack 1 becomes smaller than that for crack 2.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.215-222
• /
• 2019

The heat conduction problem with discontinuous material coefficients generally consists of the conservative equation, boundary condition, and interface condition, which should be additionally satisfied in the solution procedure. This feature often makes the development of new numerical schemes difficult as it induces a layered singularity in the solution fields; thus, a special approximation is required to capture the singular behavior. In addition to the approximation, the construction of a total system of equations is challenging. In this study, a wedge function is devised for enriching the approximation, and the interface condition itself is embedded in the moving least squares(MLS) derivative approximation to consistently satisfy the interface condition. The heat conduction problem is then discretized in a strong form using the developed derivative approximation, which is named as the interface immersed MLS difference method. This method is able to efficiently provide a numerical solution for such interface problems avoiding both numerical quadrature as well as extra difference equations related to the interface condition enforcement. Numerical experiments proved that the developed numerical method was highly accurate and computationally efficient at solving the heat conduction problem with interfacial jump as well as the problem with a geometrically induced interfacial singularity.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.38-46
• /
• 1990

A numerical solution method of the boundary integral equation for axisymmetric potential flows is presented. Those are represented by ring source and ring vorticity distribution. Strengths of ring source and ring vorticity are approximated by linear functions of a parameter $\zeta$ on a segment. The geometry of the body is represented by a cubic B-spline. Limiting integral expressions as the field point tends to the surface having ring source and ring vorticity distribution are derived upto the order of ${\zeta}ln{\zeta}$. In numerical calculations, the principal value integrals over the adjacent segments cancel each other exactly. Thus the singular part proportional to $\(\frac{1}{\zeta}\)$ can be subtracted off in the calculation of the induced velocity by singularities. And the terms proportional to $ln{\zeta}$ and ${\zeta}ln{\zeta}$ can be integrated analytically. Thus those are subtracted off in the numerical calculations and the numerical value obtained from the analytic integrations for $ln{\zeta}$ and ${\zeta}ln{\zeta}$ are added to the induced velocity. The four point Gaussian Quadrature formula was used to evaluate the higher order terms than ${\zeta}ln{\zeta}$ in the integration over the adjacent segments to the field points and the integral over the segments off the field points. The root mean square errors, $E_2$, are examined as a function of the number of nodes to determine convergence rates. The convergence rate of this method approaches 2.