• Plant Journal
• /
• v.17 no.3
• /
• pp.42-46
• /
• 2021

The equation of motion of the drill string along the excavation trajectory was analyzed using the Lagrangian approach together with the finite element method (FEM). A drill string of circular cross section is constructed by combining a plurality of circular axes each having 12 degrees of freedom (DOF). FEM analysis can observe the vibration and dynamic changes of the entire drill string, and it is easy to apply comprehensive boundary conditions to reproduce the simulation of a realistic drill string. In this study, the constructed FEM motel was simulated. In order to apply the FEM program to the actual drill trajectory, the dynamic analysis of the curved beam was verified by comparison with the actual values. The dynamic change over time was observed.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.9-16
• /
• 2014

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The "Hamilton-Jacobi(H-J)" equation and computationally robust numerical technique of "up-wind scheme" lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H-J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes are not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.1
• /
• pp.239-246
• /
• 2017

There are various transversal structures (small dams or drop structures) in median and small streams in Korea. Most of them are concrete structures and it is so hard to exclude low-level water. Unless drainage valves and/or gates would not be installed near bottom of bed, sediment from upstream should be deposited and also contaminants attached to the sediments would devastatingly threaten the water quality and ecosystem. One of countermeasures for such problem is the bypass pipe installed underneath the transversal structure. However, there is still issued whether it would be workable if the gravels and/or stones would roll into and be not excluded. Therefore, in this study, the conditions to exclude the rip stone which enter into the bypass pipe was reviewed. Based on sediment transport phenomenon, the behavior of stones was investigated with the concepts from the critical shear stress of sediment and d'Alembert principle. As final results, the basis condition (${\tau}_c{^*}$) was derived using the Lagrangian description since the stones are in the moving state, not in the stationary state. From hydraulic experiments the relative velocity could be obtained. In order to minimize the scale effect, the extra wide channel of 5.0 m wide and 1.0 m high was constructed and the experimental stones were fully spherical ones. Experimental results showed that the ratio of flow velocity to spherical particle velocity was measured between 0.5 and 0.7, and this result was substituted into the suggested equation to identify the critical condition wether the stones were excluded. Regimes about the exclusion of stone in bypass pipe were divided into three types according to particle Reynolds number ($Re_p$) and dimensionless critical shear force (${\tau}_c{^*}$) - exclusion section, probabilistic exclusion section, no exclusion section. Results from this study would be useful and essential information for bypass pipe design in transveral structures.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.33 no.1
• /
• pp.9-26
• /
• 1997

Stokes drift(SD) and Lagrangian discharge(LD) are important factors for analysis of flushing time, tidal exchange, solute transport and pollutant dispersion. The factors should be calculated using the approached method to flow phenomena. The aim of this paper re-examines the previous procedures for computing the SD and LD, and is to propose the new method approached to stratified flow field in the cross-section of coastal region, e.g. Masan Bay. The intensity of velocity near the bottom boundary layer(BBL) depends on the sea-bed irregularity in the coastal estuaries. So we calculated the depth mean velocity(DMV) considering that of BBL omitted in Kjerfve's calculation method. It revealed that BBL effect resulting in application of the bay acts largely on DMV in half more among 1l stations. The new expression of SD and LD per unit width in the cross-section using the developed DMV and proposed decomposition procedure of current were derived as follow : $$Q=u_0+\frac{1}{2}H_1{U_1cos(\varphi_h-\varphi_u)+U_3cos(\varphi_h-\varphi{ud})} LD ED SD$(Q_{skim}+Q_{sk2}) The third term, $Q_{sk2}$, on the right-hand of the equation is showed newly and arise from vertical oscillatory shear. According to the results applied in 3 cross-sections including 11 stations of the bay, the volume difference between proposed and previous SD was founded to be almost 2 times more at some stations. But their mean transport volumes over all stations are 18% less than the previous SD. Among two terms of SD, the flux of second term, $Q_{skim}$, is larger than third term, $Q_{sk2}$, in the main channel of cross-section, so that $Q_{skim}$ has a strong dependence on the tidal pumping, whereas third term is larger than second in the marginal channel. It means that $Q_{sk2}$ has trapping or shear effect more than tidal pumping phenomena. Maximum range of the fluctuation in LD is 40% as compared with the previous equations, but mean range of it is showed 11% at all stations, namely, small change. It mean that two components of SD interact as compensating flow. Therefore, the computation of SD and LD depend on decomposition procedure of velocity component in obtaining the volume transport of temporal and spacial flow through channels. The calculation of SD and LD proposed here can separate the shear effect from the previous SD component, so can be applied to non-uniform flow condition of cross-section, namely, baroclinic flow field.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.39 no.4
• /
• pp.106-113
• /
• 1997

Dynamic loading of structures often causes excursions of stresses well into the inelastic range, and the influence of the geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. A structure in a nuclear power plant is a structure of importance which puts emphasis on safety. A nuclear container is a pressure vessel subject to internal pressure and this structure is constructed by a reinforced concrete or a pre-stressed concrete. In this study, the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion is numerically solved by a central difference scheme. The constitutive relation of concrete is modeled according to a Drucker-Prager yield criterion in compression. The reinforcing bars are modeled by a smeared layer at the location of reinforcements, and the steel layer model under Von Mises yield criteria is adopted to represent an elastic-plastic behavior. To investigate the dynamic response of a nuclear reinforced concrete containment structure, the steel-ratios of 0, 3, 5 and 10 percent, are considered. The results obtained from the analysis of an example were summarized as follows 1. As the steel-ratio increases, the amplitude and the period of the vertical displacements in apex of dome decreased. The Dynamic Magnification Factor(DMF) was some larger than that of the structure without steel. However, the regular trend was not found in the values of DMF. 2. The dynamic response of the vertical displacement and the radial displacement in the dome-wall junction were shown that the period of displacement in initial step decreased with the steel-ratio increases. Especially, the effect of the steel on the dynamic response of radial displacement disapeared almost. The values of DMF were 1.94, 2.5, 2.62 and 2.66, and the values increased with the steel-ratio. 3. The characteristics of the dynamic response of radial displacement in the mid-wall were similar to that of dome-wall junction. The values of DMF were 1.91, 2.11, 2.13 and 2.18, and the values increased with the steel-ratio. 4. The amplitude and the period of the hoop-stresses in the dome, the dome-wall junction, and the mid-wall were shown the decreased trend with the steel-ratio. The values of DMF were some larger than those of the structure without steel. However, the regular trend was not found in the values of DMF.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.17 no.1
• /
• pp.59-70
• /
• 2022

In this paper, we proposed an active fault tolerant control (AFTC) method for the position control of a quadrotor with complete loss of effectiveness of one motor. We obtained the dynamics of a quadrotor using Lagrangian equation without small angle assumption. For detecting the fault on a motor, we designed a fault detection module, which consists of the fault detection and diagnosis (FDD) module and the fault detection and isolation (FDI) module. For the FDD module, we designed a nonlinear observer that observes the states of a quadrotor based on the obtained dynamics. Using the observed states of a quadrotor, we designed residual signals and set the appropriate threshold values of residual signals to detect the fault. Also, we designed an FDI module to identify the fault location using the designed additional conditions. To make a quadrotor track the desired path after detecting the fault of a motor, we designed a fault tolerant controller based on the multiple sliding surface control (MSSC) technique. Finally, through simulations, we verified the effectiveness of the proposed AFTC method for a quadrotor with complete loss of effectiveness of one motor.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.180-193
• /
• 2021

The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.547-558
• /
• 2008

With the advent of super cargo ship due to the explosive increase in the amount of cargo shipped via seas, some mega ports are under construction in South Korea, to accommodate the super cargo ship, and some of them already enter their final phase. To sustain the harbor tranquility, mega ports usually comprise huge vertical type breakwaters which are intrinsically vulnerable to the attack of long waves. In this rationale, we present the chamber type breakwater with a circular curtain wall - Eco-breaker 2, to alleviate the reflection of long waves and numerically investigate the hydraulic characteristics of Eco-breaker 2. As a wave driver, we use the Navier-Stokes eq., the most robust wave driver, using SPH (Smoothed Particle Hydrodynamics) and LES (Large Eddy Simulation). For the verification of numerical results, we also carried out hydraulic model test. It is shown that Eco-breaker 2 can effectively alleviate the reflection of long waves with its inherited large organized eddies encompassing the water chamber and some region off the curtain wall of varying size. It is also shown that the scope and strength of large organized eddies strongly depends on the incident wave period, and the reflection coefficient can be lowered to 0.18 by tuning the size of water chamber such that resident time at the chamber is just short of the half period of incident waves. Based on these results, we present the specification of Eco-breaker 2 to boost its use on the development of water environment friendly harbor worldwide.