• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.1
• /
• pp.39-45
• /
• 1989

Many mechanical systems including ships and/or offshore structures have poor dynamic response characteristics such as undesirable natural frequencies and undesirable mode shapes. It is mandatory to redesign the structure. In this paper a procedure for the dynamic redesign of an undamped structural system is presented. The method which uses a penalty function with a penalty term containing error in equilibrium for a given vibration mode may have a shortcoming. This method includes unconstrained eigenvector degrees of freedom as unknowns. In the work developed here, only constrained mode shape changes are used in the solution procedure, resulting in a reduction of the unnecessary calculations. Among the set of equations which characterizes the redesign of the structural systems, the under constrained problem is discussed here and formulated as an optimization problem, with an optimal criterion such as minimum change or minimum structural weight of the system. Four simple numerical applications illustrate the efficiency of the method. The method can be applied to the vibration problems of ships and/or offshore structures with an implementation of the commercial FE codes.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• Journal of Wetlands Research
• /
• v.7 no.3
• /
• pp.57-66
• /
• 2005

This paper is to present the determination of the optimal loss rate parameters and unit hydrographs from the observed single rainfall-runoff event using optimization model. The linear program models has been formulated to derive the optimal unit hydrographs and loss rate parameters for the site of the Su-Jik Bridge in the HwangGuJichen River; one minimizes the summation of the absolute residual between predicted and observed runoff ordinates. In the perturbation stage of parameters the trial and error method has been adopted to determine the loss rate parameters for Kostiakov's, Philip's, Horton's, and Green-Ampt's equation. The unique unit hydrograph ordinates for a given rainfall-runoff event is exclusively obtained with ${\Phi}$ index, but unit hydrograph ordinates depend upon the parameters for each loss rate equations. In this paper the single rainfall-runoff event observed from the sample watershed is considered to test the proposed method. The optimal unit hydrograph obtained by the optimization model has smaller deviations than the ones by the conventional method.

• Journal of Soil and Groundwater Environment
• /
• v.9 no.1
• /
• pp.79-86
• /
• 2004

Hydraulic conductivity along rock fracture is mainly dependent on fracture geometries such as orientation, aperture, roughness and connectivity. Therefore, it needs to consider fracture geometries sufficiently on a fracture model for a numerical analysis to calculate permeability coefficient in a fracture. This study performed new type of numerical analysis using a homogenization analysis method to calculate permeability coefficient accurately along single fractures with several fracture models that were considered fracture geometries as much as possible. First of all, fracture roughness and aperture variation due to normal stress applied on a fracture were directly measured under a confocal laser scaning microscope (CLSM). The acquired geometric data were used as input data to construct fracture models for the homogenization analysis (HA). Using the constructed fracture models, the homogenization analysis method can compute permeability coefficient with consideration of material properties both in microscale and in macroscale. The HA is a new type of perturbation theory developed to characterize the behavior of a micro inhomogeneous material with a periodic microstructure. It calculates micro scale permeability coefficient at homogeneous microscale, and then, computes a homogenized permeability coefficient (C-permeability coefficient) at macro scale. Therefore, it is possible to analyze accurate characteristics of permeability reflected with local effect of facture geometry. Several computations of the HA were conducted to prove validity of the HA results compared with the empirical equations of permeability in the previous studies using the constructed 2-D fracture models. The model can be classified into a parallel plate model that has fracture roughness and identical aperture along a fracture. According to the computation results, the conventional C-permeability coefficients have values in the range of the same order or difference of one order from the permeability coefficients calculated by an empirical equation. It means that the HA result is valid to calculate permeability coefficient along a fracture. However, it should be noted that C-permeability coefficient is more accurate result than the preexisting equations of permeability calculation, because the HA considers permeability characteristics of locally inhomogeneous fracture geometries and material properties both in microscale and macroscale.

• Journal of Astronomy and Space Sciences
• /
• v.22 no.1
• /
• pp.21-34
• /
• 2005

The initial conditions that generate bounded motion in eccentric reference orbit are determined for satellite formation flying. Because Hill's equations cannot describe the relative motion between two satellites in eccentric orbit, a new relative dynamics utilizing the nonlinearity and eccentricity correction for Hill's initial conditions is implemented. The constraint that matches angular rates of chief and deputy satellites is used to obtain the bounded motion between them. The constraint can be applied to satellite formation motions in eccentric orbit, since it implicates J2 perturbation due to the central body's aspherical gravitational forces. The periodic bounded motions are analyzed for the orbit with the eccentricity of less than 0.05 and about 0.5 km relative distance between chief and deputy satellites. It is mainly illustrated that the satellite formations in small eccentric orbits can have hounded motions; consequently, the formation can be kept by matching angular rates of the satellites. These results demonstrate an useful method that reduces the cost for operating satellites by providing effective initial conditions for satellite formation flying in eccentric orbit.

• KSTLE International Journal
• /
• v.7 no.2
• /
• pp.27-34
• /
• 2006

Knudsen number is the ratio of molecular mean free path versus mm thickness and the criterion to determine the flow form. When its value is lower than 0.01, the flow can be assumed to has no slip boundary condition. And in the case that the value is between 0.01 and 10, then the flow has slip boundary condition at both the adjacent walls. The condition of the air flow between the rotating journal and top foil in the air foil bearing is determined by the rotating speed and load, and the Knudsen number is also varied by those values. Because the molecular mean free path is variable to the pressure and temperature, more exact formulation is necessary to understand and analyze the flow regime. In this study, the analysis considering Knudsen number formulated with those variables (pressure, temperature and mm thickness) was executed. The approximate value was examined using the equation to confirm whether the flow has the slip or no-slip boundary condition. From the analytic investigation, it was decided to range approximately 0.01 to 1.0 and the flow can be supposed to have the slip boundary condition. Under the condition of the slip flow, the static characteristics of the air foil bearing were examined using modified Reynolds equations. The results were compared with those considering no slip condition. It shows that the slip condition makes the flow decelerates and the load carrying capacity decreases compared with no slip condition. And as the bearing number and eccentricity ratio increase, the load carrying capacity also increased at both the cases. From this result, it can be supposed that the bearing torque also increases. In the analysis of the dynamic characteristics, the perturbed Knudsen number was taken into consideration. Because the Knudsen number is expressed as the terms of each variable, the perturbed equation can be simply derived. The results of both cases considering and not considering Knudsen number were compared each other. In the case of the direct terms of the stiffness and damping coefficients, the difference between both cases was little and increased as the bearing number and eccentricity ratio increased. And the cross terms have less or more differences.

• Korean Journal of Optics and Photonics
• /
• v.16 no.3
• /
• pp.182-190
• /
• 2005

The fiber Raman amplifier(FRA) is a distinctly advantageous technology. Due to its wider, flexible gain bandwidth, and intrinsically lower noise characteristics, FRA has become an indispensable technology of today. Various FRA modeling methods, with different levels of convergence speed and accuracy, have been proposed in order to gain valuable insights for the FRA dynamics and optimum design before real implementation. Still, all these approaches share the common platform of coupled ordinary differential equations(ODE) for the Raman equation set that must be solved along the long length of fiber propagation axis. The ODE platform has classically set the bar for achievable convergence speed, resulting exhaustive calculation efforts. In this work, we propose an alternative, highly efficient framework for FRA analysis. In treating the Raman gain as the perturbation factor in an adiabatic process, we achieved implementation of the algorithm by deriving a recursive relation for the integrals of power inside fiber with the effective length and by constructing a matrix formalism for the solution of the given FRA problem. Finally, by adiabatically turning on the Raman process in the fiber as increasing the order of iterations, the FRA solution can be obtained along the iteration axis for the whole length of fiber rather than along the fiber propagation axis, enabling faster convergence speed, at the equivalent accuracy achievable with the methods based on coupled ODEs. Performance comparison in all co-, counter-, bi-directionally pumped multi-channel FRA shows more than 102 times faster with the convergence speed of the Average power method at the same level of accuracy(relative deviation < 0.03dB).