• Journal of Navigation and Port Research
• /
• v.43 no.5
• /
• pp.328-334
• /
• 2019

Underwater architecture in providing a comfortable living space underwater is mandated to survive prevailing environmental loads, especially hydrostatic ambient water pressure exerted on the structure of individual habitat hulls at depth and hydrodynamic fluctuation of external forces that perturb the postural equilibrium and mooring stability of the underwater housing system, for which the design including the hull shape and mooring system constraint the responses. In this study, the postural stability of a proposed underwater floating housing system with three vertically connected ellipsoidal-shape habitat hulls of different sizes are theorized and calculated for hydrostatic stability, using MATLAB in the volumetric integration of a hull and the weight of operational loads under assumed scenarios. The assumptions made in the numerical method to estimate the stability of the habitat system include the fixed weight of the hulls, and their adjustable loads within operational limits for the set meteorological oceanic conditions. The purpose of this study was to numerically manipulate a) The buoyancy and b) The adjusted center of mass of the system within the range of designed external and internal load changes, by which the effective mooring system capability and postural equilibrium requirements were argued with the quantitative analysis.

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

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.9
• /
• pp.3119-3124
• /
• 2010

The purpose of this paper is to analyze dynamic behaviors of the oil film thickness and engine bearings in both aligned and misaligned operation conditions of a crankshaft using computer simulation techniques. A crankshaft as an elastic body is modeled for a misaligned crankshaft which is very important design parameter of the film thickness and engine bearings. In this analysis, a dynamic characteristic of a minimum oil film is analyzed based on the elastohydrodynamic lubrication theory. The boundary conditions for analyzing the film behaviors include non-linear constraint forces and bending moments in engine bearings. The more expedient model of an engine bearing is extended to consider the effect of crankshaft misalignment. The computed results indicate that the minimum oil film thickness that causes a major influence on the performance of engine bearings has showed a decrease of 16% to 24% for the misaligned crankshaft compared with an aligned crankshaft. The computed results show that the misalignment of a crankshaft inevitably brings the reduction of minimum oil film thickness and this may increase the failure of a bearing. These results as design parameters are very useful for a bearing designer as a firm reference data of an automotive engine.

• Journal of Korean Society of Steel Construction
• /
• v.16 no.1 s.68
• /
• pp.123-134
• /
• 2004

In this research, a multilevel decomposition technique to enhance the efficiency of the configuration optimization of truss structures was proposed. On the first level, the nonlinear programming problem was formulated considering cross-sectional areas as design variables, weight, or volume as objective function and behavior under multiloading condition as design constraint. Said nonlinear programming problem was transformed into a sequential linear programming problem. which was effective in calculation through the approximation of member forces using behavior space approach. Such approach has proven to be efficient in sensitivity analysis and different form existing shape optimization studies. The modified method of feasible direction (MMFD) was used for the optimization process. On the second level, by treating only shape design variables, the optimum problem was transformed into and unconstrained optimal design problem. A unidirectional search technique was used. As numerical examples, some truss structures were applied to illustrate the applicability. and validity of the formulated algorithm.

• Computers and Concrete
• /
• v.14 no.2
• /
• pp.127-143
• /
• 2014

LOCA (Loss Of Coolant Accident) is one of the most important utmost accidents for Prestressed Concrete Containment Vessel (PCCV) due to its coupled effect of high temperature and inner pressure. In this paper, heat conduction analysis is used to obtain the LOCA temperature distribution of PCCV. Then the elastic internal force of PCCV under LOCA temperature is analyzed by using both simplified theoretical method and FEM (finite element methods) method. Considering the coupled effect of LOCA temperature, a nonlinear elasto-plasitic analysis is conducted for PCCV under utmost internal pressure considering three failure criteria. Results show that the LOCA temperature distribution is strongly nonlinear along the shell thickness at the early time; the moment result of simplified analysis is well coincident with the one of numerical analysis at weak constraint area; while in the strong constrained area, the value of moments and membrane forces fluctuate dramatically; the simplified and numerical analysis both show that the maximum moment occurs at 6hrs after LOCA.; the strain of PCCV under LOCA temperature is larger than the one of no temperature under elasto-plastic analysis; the LOCA temperature of 6hrs has the greatest influence on the ultimate bearing capacity with 8.43% decrease for failure criteria 1 and 2.65% decrease for failure criteria 3.

• Structural Engineering and Mechanics
• /
• v.59 no.6
• /
• pp.1139-1153
• /
• 2016

Elastic constraints are usually simplified as "spring forces" exerted on beam ends without considering the "spring deformation". The partial differential equation governing the free vibrations of a cantilever Bernoulli-Euler beam considering the deformation of elastic constraints is firstly established, and is nondimensionalized to obtain two dimensionless factors, $k_v$ and $k_r$, describing the effects of elastically vertical and rotational end constraints, respectively. Then the frequency equation for the above Bernoulli-Euler beam model is derived using the method of separation of variables. A numerical analysis method is proposed to solve the transcendental frequency equation for the continuous change of the frequency with $k_v$ and $k_r$. Then the mode shape functions are given. Finally, effects of $k_v$ and $k_r$ on free vibration characteristics of the beam with different slenderness ratios are calculated and analyzed. The results indicate that the effects of $k_v$ are larger on higher-order free vibration characteristics than on lower-order ones, and the impact strength decreases with slenderness ratio. Under a relatively larger slenderness ratio, the effects of $k_v$ can be neglected for the fundamental frequency characteristics, while cannot for higher-order ones. However, the effects of $k_r$ are large on both higher- and lower-order free vibration characteristics, and cannot be neglected no matter the slenderness ratio is large or small.

• Transactions of Materials Processing
• /
• v.23 no.8
• /
• pp.501-506
• /
• 2014

Progressive shearing with blanking dies is commonly employed to produce large quantities of tiny sheet metal electronic parts. Sheet metal pins, which are narrow and long, that are sheared with a progressive die set are often twisted. The twist in the sheet metal pins, which usually occurs in the final shearing operation, generally decreases with increasing blank holding force. The blank holding forces in all shearing operations are not the same because of different shearing positions and areas. In the current study, the optimal layout of the springs in a progressive die set to minimize the twist of the sheet metal pin is proposed. In order to find the holding force acting on the tiny narrow blanks produced with the proposed springs during the shearing process, the equivalent area method is used in the structural analysis. The shearing of the sheet-metal pin was simulated to compute the twist angle associated with the blank holding force. The constraint condition satisfying the pre-set blank holding force from the previous shearing operations was imposed. A design of experiments (DOE) was numerically implemented by analyzing the progressive die structure and by simulating the shearing process. From the meta-model created from the experimental results and by using a quadratic response surface method (PQRSM), the optimal layout of the springs was determined. The twist of sheet metal pin associated with the optimal layout of the springs found in the current study was compared with that of an existing progressive die to obtain a minimal amount of twist.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.6
• /
• pp.37-46
• /
• 2001

In this research, a controller design method based on optimization is proposed that can satisfy constraints on maximum responses of building structures subject to around excitation modeled by partially stochastic process. The class of controllers to be optimized is restricted to LQR. Weighting matrix on controlled outputs is used as design variable. Objective function, constraint functions and their gradients are computed by the parameterization of control gain with Riccati matrix. Full state feedback controllers designed by proposed optimization method satisfy various design objectives and their necessary maximum control forces are computed for the production of actuator. LQG controllers composed of Kalman filter and LQR designed by proposed method perform well with little deterioration. So it is possible to design output feedback controllers satisfying constraints on various maximum responses of structures.

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.39 no.1
• /
• pp.98-104
• /
• 2016

With the development of modern science and technology, weapon systems such as tanks, submarines, combat planes, radar are also dramatically advanced. Among these weapon systems, the ballistic missile, one of the asymmetric forces, could be considered as a very economical means to attack the core facilities of the other country in order to achieve the strategic goals of the country during the war. Because of the current ballistic missile threat from the North Korea, establishing a missile defense (MD) system becomes one of the major national defense issues. This study focused on the optimization of air defense artillery units' deployment for effective ballistic missile defense. To optimize the deployment of the units, firstly this study examined the possibility of defense, according to the presence of orbital coordinates of ballistic missiles in the limited defense range of air defense artillery units. This constraint on the defense range is originated from the characteristics of anti-ballistic missiles (ABMs) such as PATRIOT. Secondly, this study proposed the optimized mathematical model considering the total covering problem of binary integer programming, as an optimal deployment of air defense artillery units for defending every core defense facility with the least number of such units. Finally, numerical experiments were conducted to show how the suggested approach works. Assuming the current state of the Korean peninsula, the study arbitrarily set ballistic missile bases of the North Korea and core defense facilities of the South Korea. Under these conditions, numerical experiments were executed by utilizing MATLAB R2010a of the MathWorks, Inc.

• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.