• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.1
• /
• pp.106-115
• /
• 1988

A numerical procedure for determining the airfoil shape that maximizes the lift is presented. The structure of the flow field is calculated by iteratively coupling potential flow and boundary analysis using the viscous-inviscid interaction method. The potential flow field is obtained by the vortex panel method and boundary layer flow is analyzed by means of integral approximation method which is capable of handling the laminar, transition and turbulent flow regimes. As the result of this study, it is found that the calculated flow regimes have good agreement with the existing experimented data. Davidon-Fletcher-Powell method and Augmented Lagrange Multiplier method are used for the optimal techniques. NACA 23012, NACA 65-3-21, NACA 64-2-415, NACA 64-2-A215 airfoils are used for determining the optimal airfoil shapes as a basic and compensate airfoils. Optimal design showed that the lift coefficients are increased by 17.4% at M$_{0}$=0.2 and 29% at M$_{0}$=0.3, compared with those of basic airfoil.oil.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.25 no.2
• /
• pp.245-250
• /
• 2014

This paper introduces the inverse synthetic aperture radar(ISAR) imaging technique for rear view target of an automobile, which uses both linear frequency modulation-frequency shift keying(LFM-FSK) waveform and monopulse tracking. LFM-FSK waveform consists of two sequential stepped frequency waveforms with some frequency offset, and thus, can be used to generate ISAR images of rear view target of an automobile. However, ISAR images can often be blurred due to non-uniform change rate of relative aspect angle between radar and target. In order to address this problem, one-dimensional(1-D) Lagrange interpolation technique in conjunction with angle information obtained from the monopulse tracking is applied to generate uniform data across the radar's aspect angle. Simulation results show that the proposed method can provide focused ISAR images.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.34 no.7A
• /
• pp.562-567
• /
• 2009

In this paper, two-hop beamforming relay system, with source, relay, and destination nodes, is considered and the transmit- and receive-beamforming vectors and the relay processing matrix are designed for minimizing a mean square error (MMSE) between the transmit and receive signals. Here, to reduce the transmit power of the source or the relay, two local inequality constraints are involved with MMSE problem. By adopting the Lagrange method, closed formed Karush-Kuhn-Tucker (KKT) conditions (equalities) are derived and an iterative algorithm is developed to solve the entangled KKT equalities. Due to the inequality power constraints, the source or the relay can reduce its transmit power when the received signal-to-noise ratios (SNRs) of the first- and the second-hop are different. Meanwhile, the destination can achieve almost identical bit-error-rate performance compared to an optimal beamforming system maximizing the received SNR. This claim is supported by a computer simulation.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.9
• /
• pp.775-780
• /
• 2016

In this paper, a computational algorithm of an improved and versatile structural analysis applicable for large-size flexible nonlinear structures is developed. In more detail, nonlinear finite element based on the co-rotational (CR) framework is developed. Then, a finite element tearing and interconnecting method using local Lagrange multipliers (FETI-local) is combined with the nonlinear CR finite element. The resulting computational algorithm is presented and applied for nonlinear static analyses, i.e., cantilevered beam and multibody structure. Finally, the proposed analysis is evaluated with regard to its parallel computation performance, and it is compared with those obtained by serial computation using the sparse matrix linear solver, PARDISO.

• Journal of the Korean Geosynthetics Society
• /
• v.20 no.2
• /
• pp.35-43
• /
• 2021

In this paper, the sub-structuring technique-applied train-bridge interaction analysis model, which is formulated based on the simplified three-dimensional train-bridge interaction analysis model for high-speed bridge-train interaction analysis, is presented. In the sub-structuring technique, the super-structure and the supporting structure of railway bridges can be modeled as sub-structures, and train-bridge interaction analysis can be efficiently performed. As a train analysis model, two-dimensional train model is used, and the Lagrange equation of motion is applied to derive the equation of motion of two-dimensional train. In the sub-structuring technique, the number of degrees of freedom can be reduced by using the condensation method, thus reducing the time and cost for calculating the eigenvalues and eigenvectors, and the time and cost for the subsequent calculation. In this paper, Guyan reduction method is used as sub-structuring technique. By combining simplified three-dimensional bridge-train interaction analysis and Guyan reduction method, the efficient and accurate bridge-train interaction analysis can be performed.

• Explosives and Blasting
• /
• v.39 no.1
• /
• pp.1-9
• /
• 2021

As the national economic growth and the rapid increase in industrial structures are aging, the demand for removing steel structures is increasing, and research on improving the penetration performance of the linear shape charge explosives. In the study, numerical analyses were performed on the effect of the type of explosive used in the self-made shape charging container and the initiation method on the cutting performance of the steel plate and the effect on the shaped explosive installed close to it. ANSYS LS-DYNA, which can analyze the large deformation problem of materials due to explosion, was used, and an ALE(Arbitrary-Lagrange-Eulerian) model was applied that enables interlocking analysis of gases, liquids, and solid.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.191-197
• /
• 2021

To calculate the induced charge of atoms in molecular dynamics, linear equations for the induced charges need to be solved. As induced charges are determined at each time step, the process involves considerable computational costs. Hence, an efficient method for calculating the induced charge distribution is required when analyzing large systems. This paper introduces the Uzawa method for solving saddle point problems, which occur in linear systems, for the solution of the Lagrange equation with constraints. We apply the accelerated Uzawa algorithm, which reduces computational costs noticeably using the Schur complement and preconditioned conjugate gradient methods, in order to overcome the drawback of the Uzawa parameter, which affects the convergence speed, and increase the efficiency of the matrix operation. Numerical models of molecular dynamics in which two gold nanoparticles are placed under external electric fields reveal that the proposed method provides improved results in terms of both convergence and efficiency. The computational cost was reduced by approximately 1/10 compared to that for the Gaussian elimination method, and fast convergence of the conjugate gradient, as compared to the basic Uzawa method, was verified.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.6A
• /
• pp.989-999
• /
• 2006

In this paper, numerical analysis method to perform linear dynamic analysis of bridge considering the road surface roughness and bridge-vehicle interaction when vehicle is moving on bridge is presented. The vehicle and bridge are modeled as three-dimension where contact length of tire and pitching of tandem spring are considered and single truck with 2-axles and 3- axles, and tractor-trailer with 5-axles are modeled as 7-D.O.F., 8-D.O.F., and 14-D.O.F., respectively. Dynamic equations of vehicle are derived from the Lagrange's equation and solution of the equation is obtained by Newmark-${\beta}$ method. The surface roughness of bridge deck for this analysis is generated from power spectral density (PSD) function. Beam element for the main girder, shell element for concrete deck and rigid link between main girder and concrete deck are used. The equations of the motion of bridges are solved by mode-superposition procedures. The proposed procedure is validated by comparing the results with the experimental data by Whittemore and Fenves.

• International Journal of Precision Engineering and Manufacturing
• /
• v.7 no.1
• /
• pp.24-29
• /
• 2006

In this paper, the effect of open crack on the dynamic behavior of simply supported Timoshenko beam with a moving mass was studied. The influences of the depth and the position of the crack on the beam were studied on the dynamic behavior of the simply supported beam system by numerical methods. The equation of motion is derived by using Lagrange's equation. The crack section is represented by a local flexibility matrix connecting two undamaged beam segments. The crack is modeled as a rotational spring. This flexibility matrix defines the relationship between the displacements and forces on the crack section and is derived by applying fundamental fracture mechanics theory. As the depth of the crack increases, the mid-span deflection of the Timoshenko beam with a moving mass is increased.