• Ocean Systems Engineering
• /
• v.7 no.2
• /
• pp.143-155
• /
• 2017

For complex flexible structures such as nets, the determination of drag forces and its deformation is a challenging task. The accurate prediction of loads on cages is one of the key steps in designing fish farm facilities. The basic physics with a simple cage, can be addressed by the use of experimental studies. However, to design more complex cage system for various environmental conditions, a reliable numerical simulation tool is essential. In this work, the current load on a cage is calculated using a Morison-force model applied at instantaneous positions of equivalent-net modeling. Variations of solidity ratio ($S_n$) of the net and current speed are considered. An equivalent array of cylinders is built to represent the physical netting. Based on the systematic comparisons between the published experimental data for Raschel nets and the current numerical simulations, carried out using the commercial software OrcaFlex, a new formulation for $C_d$ values, used in the equivalent-net model, is presented. The similar approach can also be applied to other netting materials following the same procedure. In case of high solidity ratio and current speed, the hybrid model defines $C_d$ as a function of Re (Reynolds number) and $S_n$ to better represent the corresponding weak diffraction effects. Otherwise, the conventional $C_d$ values depending only on Re can be used with including shielding effects for downstream elements. This new methodology significantly improves the agreement between numerical and experimental data.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• v.1
• /
• pp.181-186
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. The present model of the fluid motion is based on the Navier-Stokes equations incorporating velocity-pressure formulation because of need to model the nonlinear wave interaction with spatially non-uniform current field. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. As an intermediate stage of development, solution procedure and characteristics aspects of the present modeling and numerical method are addressed in detail, and preliminary numerical results prove its accuracy and convergence.

• International Journal of Ocean System Engineering
• /
• v.2 no.2
• /
• pp.116-138
• /
• 2012

As in the other fields of mechanics, analytical methods represent an important analysis tool in marine hydrodynamics. The analytical approach is interesting for different reasons : it gives reference results for numerical codes verification, it gives physical insight into some complicated problems, it can be used as a simplified predesign tool, etc. This approach is of course limited to some simplified geometries (cylinders, spheres, ...), and only the case of one or more cylinders, truncated or not, will be considered here. Presented methods are basically eigenfunction expansions whose complexity depends on the boundary conditions. The hydrodynamic boundary value problem (BVP) is formulated within the usual assumptions of potential flow and is additionally simplified by the perturbation method. By using this approach, the highly nonlinear problem decomposes into its linear part and the higher order (second, third, ...) corrections. Also, periodicity is assumed so that the time dependence can be factorized i.e. the frequency domain formulation is adopted. As far as free surface flows are concerned, only cases without or with small forward speed are sufficiently simple to be solved semi-analytically. The problem of the floating body advancing in waves with arbitrary forward speed is far more complicated. These remarks are also valid for the general numerical methods where the case of arbitrary forward speed, even linearized, is still too difficult from numerical point of view, and "it is fair to say that there exists at present no general practical numerical method for the wave resistance problem" [9], and even less for the general seakeeping problem. We note also that, in the case of bluff bodies like cylinders, the assumptions of the potential flow are justified only if the forward speed is less than the product of wave amplitude with wave frequency.

• Geomechanics and Engineering
• /
• v.2 no.1
• /
• pp.1-18
• /
• 2010

The analysis of structure response and design of buried structures subjected to dynamic destructive loads have been receiving increasing interest due to recent severe damage caused by strong earthquakes and terrorist attacks. For a comprehensive design of buried structures subjected to blast loads to be conducted, the whole system behaviour including simulation of the explosion, propagation of shock waves through the soil medium, the interaction of the soil with the buried structure and the structure response needs to be simulated in a single model. Such a model will enable more realistic simulation of the fundamental physical behaviour. This paper presents a complete model simulating the whole system using the finite element package ABAQUS/Explicit. The Arbitrary Lagrange Euler Coupling formulation is used to model the explosive charge and the soil region near the explosion to eliminate the distortion of the mesh under high deformation, while the conventional finite element method is used to model the rest of the system. The elasto-plastic Drucker-Prager Cap model is used to model the soil behaviour. The explosion process is simulated using the Jones-Wilkens-Lee equation of state. The Concrete Damage Plasticity model is used to simulate the behaviour of concrete with the reinforcement considered as an elasto-plastic material. The contact interface between soil and structure is simulated using the general Mohr-Coulomb friction concept, which allows for sliding, separation and rebound between the buried structure surface and the surrounding soil. The behaviour of the whole system is evaluated using a numerical example which shows that the proposed model is capable of producing a realistic simulation of the physical system behaviour in a smooth numerical process.

• Journal of the Korean Geotechnical Society
• /
• v.39 no.11
• /
• pp.23-31
• /
• 2023

The escalating settlements observed in concrete slab tracks pose a significant challenge in Korea, raising concerns about their adverse impact on the safe operation of high-speed railways and the substantial costs involved in restoration. A primary contributor to these settlements is identified as the utilization of rock materials sourced from tunnel construction, incorporated into the lower subgrade without the requisite soil mixing to achieve an appropriate particle size distribution. This study employs numerical analysis to evaluate the efficacy of partial reinforcement in reducing settlements in rock-filled lower subgrades. Column-shaped reinforcement areas strategically positioned at regular intervals in the lower subgrade induce soil arching in the upper subgrade, leading to a concentration of soil loads on the reinforced areas and consequent settlement reduction. The analysis employs finite element methods to investigate the influence of the size, stiffness, and spacing of the reinforced areas on settlement reduction in the lower subgrade. The numerical results guide the formulation of an optimal design approach, proposing a method to determine the minimum spacing required for reinforcements to effectively limit settlements within acceptable bounds. This research contributes valuable insights into addressing the challenges associated with settlement in concrete slab tracks, offering a basis for informed decision-making in railway infrastructure management.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.181-197
• /
• 2024

A numerical method is presented in this paper, for buckling analysis of thin arbitrary stiffened composite cylindrical shells under axial compression. The stiffeners can be placed inside and outside of the shell. The shell and stiffeners are operated as discrete elements, and their interactions are taking place through the compatibility conditions along their intersecting lines. The governing equations of motion are obtained based on Koiter's theory and solved by utilizing the principle of the minimum potential energy. Then, the buckling load coefficient and the critical buckling load are computed by solving characteristic equations. In this formulation, the elastic and geometric stiffness matrices of a single curved strip of the shell and stiffeners can be located anywhere within the shell element and in any direction are provided. Moreover, five stiffened composite shell specimens are made and tested under axial compression loading. The reliability of the presented method is validated by comparing its numerical results with those of commercial software, experiments, and other published numerical results. In addition, by using the ANSYS code, a 3-D finite element model that takes the exact geometric arrangement and the properties of the stiffeners and the shell into consideration is built. Finally, the effects of Poisson's ratio, shell length-to-radius ratio, shell thickness, cross-sectional area, angle, eccentricity, torsional stiffness, numbers and geometric configuration of stiffeners on the buckling of stiffened composite shells with various end conditions are computed. The results gained can be used as a meaningful benchmark for researchers to validate their analytical and numerical methods.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2002.10b
• /
• pp.123-124
• /
• 2002

A press shoe is an element of a machine for squeezing water from wood pulp in the field of manufacturing paper. This is used to compress the pulp enveloped by felt sheet with a large roller. The squeezing force is made by hydraulic pressure. The press shoe has a mechanism similar to a partial hydrostatic bearing. The pressure profile between press shoe and roller affects their squeezing ability, and partial peak pressure can tear the wet pulp. The curvature of the surface of press shoe varies to reduce the peak pressure and increase the mean pressure simultaneously. Therefore, the prediction of pressure distribution considering partially changed curvature of hydrostatic bearing is very important for designing the press shoe. In this study, the difference formulation of Reynolds' equation for partial hydrostatic bearing is by direct numerical method and a computer program to calculate the pressure distribution is developed. We investigate the effect of partially changed curvature of bearing surface on the pressure distribution. Other design parameter for hydrostatic bearing such as depth of pocket and relative velocity are also studied.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.4
• /
• pp.244-254
• /
• 1991

The numerical analysis on tile behaviour of mooring system in the offshore guyed tower is presented. The governing equilibrium equations are derived by the principle of virtual work. The drag and inertia effects of fluid are included in a Morrison type equation. The finite element method is used in the computation. Geometric nonlinearities for the analysis of the mooring line are considered in which both modified Newton-Raphson method and Newmark-$\beta$ method are employed. Numerical experiments show the validity and the capability of the developed mathematical formulation.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.123-132
• /
• 2001

An improved version of the Element-free Galerkin method (EFGM) is presented here for addressing the problem of transverse shear locking in shear-deformable beams with a high length over thickness ratio. Based upon Timoshenko's theory of thick beams, it has been recognized that shear locking will be completely eliminated if the rotation field is constructed to match the field of slope, given by the first derivative of displacement. This criterion is applied directly to the most commonly implemented version of EFGM. However in the numerical process to integrate strain energy, the second derivative of the standard Moving Least Square (MLS) shape functions must be evaluated, thus requiring at least a $C^1$ continuity of MLS shape functions instead of $C^0$ continuity in the conventional EFGM. Yet this hindrance is overcome effortlessly by only using at least a $C^1$ weight function. One-dimensional quartic spline weight function with $C^2$ continuity is therefore adopted for this purpose. Various numerical results in this work indicate that the modified version of the EFGM does not exhibit transverse shear locking, reduces stress oscillations, produces fast convergence, and provides a surprisingly high degree of accuracy even with coarse domain discretizations.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.7
• /
• pp.1120-1131
• /
• 2003

A recently developed numerical method based on a mixed volume and boundary integral equation method is applied to calculate the accurate stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks subject to external loads. Firstly, it should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. Secondly, this method takes full advantage of the capabilities developed in FEM and BIEM. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and volume integral equation method. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.