• The KIPS Transactions:PartA
• /
• v.14A no.1 s.105
• /
• pp.73-76
• /
• 2007

The determination of the conformal maps from the unit disk onto a Jordan region has been completed by solving the Theodorsen equation which is an nonlinear equation for the boundary correspondence function. Wegmann's method has been well known for the efficient mothed among the many suggestions for the Theodorsen equation. We proposed an improved method for convergence by applying a low-frequency pass filter to the Wegmann's method and theoretically proved convergence of improved iteration[1, 2]. And we proposed an effective method which makes it possible to estimate an error even if the real value is nut acquired[3]. In this paper, we propose an automatic algorithm for numerical conformal mapping bared on this error analysis in our early study. By this algorithm numerical conformal mapping is determined automatically according to the given domain of problem and the required accuracy. The discrete numbers and parameters of the low-frequency filter were acquired only by experience. This algorithm, however, is able to determine the discrete numbers and parameters of the low-frequency filter automatically in accordance with the given region This results from analyzing the function, which may decide the shape of the given domain under the assumption that the degree of the problem depends of the transformation of a given domain, as seen in the Fourier Transform. This proposed algorithm is also ploved by numerical experience.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.6
• /
• pp.561-567
• /
• 2011

The inverse radiation analysis of a two-dimensional irregular configuration using unstructured triangular meshes is presented. In this study, an enclosure filled with an absorbing, emitting and scattering medium with diffusely emitting and reflecting opaque boundaries is considered. The finite volume method is applied to solve the radiative transfer equation in order to simulate the measured incident radiation values which are used as input data for the inverse analysis. The conjugate gradient method is adopted for the estimation of wall emissivities by minimizing the objective function at each iteration step. To verify the performance of the unstructured grid system, we compare the results with those using a structured grid system for the two-dimensional lopsided shape. The effect of measurement errors on the estimation accuracy is also investigated.

• 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.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.15-26
• /
• 1990

This paper describes a potential-baoed panel method formulated for the analysis cf a supercavitating two-dimensional symmetri strut. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type, With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lifting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged.

• Structural Engineering and Mechanics
• /
• v.87 no.3
• /
• pp.255-272
• /
• 2023

Spatial cable systems can provide more transverse stiffness and torsional stiffness without sacrificing the vertical bearing capacity compared with conventional vertical cable systems, which is quite lucrative for long-span earth-anchored suspension bridges' development. Higher economy highlights the importance of refined form-finding analysis. Meanwhile, the internal connection between the lateral and vertical sags has not yet been specified. Given this, an analytic algorithm of form-finding for the earth-anchored suspension bridge with spatial cables is proposed in this paper. Through the geometric compatibility condition and mechanical equilibrium condition, the expressions for cable segment, the recurrence relationship between catenary parameters and control equations of spatial cable are established. Additionally, the nonlinear general reduced gradient method is introduced into fast and high-precision numerical analysis. Furthermore, the analytic expression of the lateral and vertical sags is deduced and discussed. This is very significant for the space design above the bridge deck and the optimization of the sag-to-span ratio in the preliminary design stage of the bridge. Finally, the proposed method is verified with the aid of two examples, one being an operational self-anchored suspension bridge (with spatial cables and a 260 m main span), and the other being an earth-anchored suspension bridge under design (with spatial cables and a 500 m main span). The necessity of an iterative calculation for hanger tensions on earth-anchored suspension bridges is confirmed. It is further concluded that the main cable and their connected hangers are in very close inclined planes.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.3 s.96
• /
• pp.251-258
• /
• 2005

The free vibration of a spinning flexible disk-spindle system supported by hydro dynamic bearings (HDB) in an HDD is analyzed by FEM. The spinning flexible disk is described using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. It is discretized by annular sector element. The rotating spindle which includes the clamp, hub, permanent magnet and yoke, is modeled by Timoshenko beam including the gyroscopic effect. The flexible supporting structure with a complex shape which includes stator core, housing, base plate, sleeve and thrust pad is modeled by using a 4-node tetrahedron element with rotational degrees of freedom to satisfy the geometric compatibility. The dynamic coefficients of HDB are calculated from the HDB analysis program, which solves the perturbed Reynolds equation using FEM. Introducing the virtual nodes and the rigid link constraints defined in the center of HDB, beam elements of the shaft are connected to the solid elements of the sleeve and thrust pad through the spring and damper element. The global matrix equation obtained by assembling the finite element equations of each substructure is transformed to the state-space matrix-vector equation, and the associated eigen value problem is solved by using the restarted Arnoldi iteration method. The validity of this research is verified by comparing the numerical results of the natural frequencies with the experimental ones. Also the effect of supporting structures to the natural modes of the total HDD system is rigorously analyzed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.11a
• /
• pp.572-578
• /
• 2003

The free vibration of a spinning flexible disk-spindle system supported by hydro dynamic bearings in a HDD is analyzed by FEM. The spinning flexible disk is described using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. It is discretized by annular sector element. The rotating spindle which includes the clamp, hub, permanent magnet and yoke, is modeled by Timoshenko beam including the gyroscopic effect. The flexible supporting structure with a complex shape which includes stator core, housing, base plate, sleeve and thrust pad is modeled by using a 4-node tetrahedron element with rotational degrees of freedom to satisfy the geometric compatibility. The dynamic coefficients of HDB are calculated from the HDB analysis program, which solves the perturbed Raynolds equation using FEM. Introducing the virtual nodes and the rigid link constraints defined in the center of HDB, beam elements of the shaft are connected to the solid elements of the sleeve and thrust pad through the spring and damper element. The global matrix equation obtained by assembling the finite element equations of each substructure is transformed to the state-space matrix-vector equation, and the associated eigenvalue problem is solved by using the restarted Arnoldi iteration method. The validity of this research is verified by comparing the numerical results of the natural frequencies with the experimental ones. Also the effect of supporting structures to the natural modes of the total HDD system is rigorously analyzed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.3
• /
• pp.49-59
• /
• 1985

Formulation of the geometric optimization for truss structures based on the elasticity theory turn out to be the nonlinear programming problem which has to deal with the Cross sectional area of the member and the coordinates of its nodes simultaneously. A few techniques have been proposed and adopted for the analysis of this nonlinear programming problem for the time being. These techniques, however, bear some limitations on truss shapes loading conditions and design criteria for the practical application to real structures. A generalized algorithm for the geometric optimization of the truss structures which can eliminate the above mentioned limitations, is developed in this study. The algorithm developed utilizes the two-phases technique. In the first phase, the cross sectional area of the truss member is optimized by transforming the nonlinear problem into SUMT, and solving SUMT utilizing the modified Newton-Raphson method. In the second phase, the geometric shape is optimized utilizing the unidirctional search technique of the Rosenbrock method which make it possible to minimize only the objective function. The algorithm developed in this study is numerically tested for several truss structures with various shapes, loading conditions and design criteria, and compared with the results of the other algorithms to examme its applicability and stability. The numerical comparisons show that the two-phases algorithm developed in this study is safely applicable to any design criteria, and the convergency rate is very fast and stable compared with other iteration methods for the geometric optimization of truss structures.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.15 no.6
• /
• pp.3448-3454
• /
• 2014

Generally existing structure analysis was applied to elastic analysis basically in practice. Considering the nonlinear material and the nonlinear geometric to be a more precise analysis, for this reason, The necessity for a structual analysis have been constantly required. Therefore, after optimization is performed, designed a simple model which is applied the principle of nonlinear in this study, a structural analysis of existing experienced users, have a aims at presenting theory and a method in order to perform anyone the analysis easily. In this study, the proposed model applied to die ribs, Regarding the shear load, less strain and stress was generated but strength was sufficient. The initial strain and stress was reconfigured to fit the size and shape, A hyperstudy in conjunction with Abaqus with nonlinear structural analysis, revealed an acceptable maximum and minimum range of stress and under the conditions of minimum strain, the plate made with a constant increment. In the experimental models, the plate thickness was given a power of 40 Newton, according to the thickness of the press die through an iterative process. When the stress and strain was applied to the die thickness, 7-8mm thickness could be obtained by optimizing.

• Journal of Korea Water Resources Association
• /
• v.49 no.1
• /
• pp.61-72
• /
• 2016

In Korea, the aquifers at the coastal areas are mostly shallow alluvial unconfined aquifers. To simulate the flow and dispersion in unconfined aquifer, a FDM model has been developed to solve the nonlinear Boussinesq equation. Related analysis and verification have been executed. The iteration method is used to solve the nonlinearity, and the model shows 3-D shape because it is a 2-D y model that consider the undulation of water table and bottom. For the verification of the model, the output of flow module is compared to the 1-D analytic solution of Lee (1989) which have the drawdown or uplift boundary condition, and the two results show almost the same value. and the mass balance of dispersion module shows about 10% error. The developed model can be used for the analysis and design of the flow and dispersion in the unconfined aquifers. The model has been applied to the estuary area of Ssangcheon watershed, and the parameters have been deduced as a result : hydraulic conductivity is 90 m/day, and longitudinal dispersivity is 15 m. And the analysis with these parameters shows that the wells are situated in the influence circle of each others except for No. 7 well. Groundwater discharge to sea is $3700m^3/day$. And the chlorine ion ($cl^-$) concentration at the pumping wells increase at least 1000 mg/L if groundwater dam is not exist, so the groundwater dam plays an important role for the prevention of sea water intrusion.