• Structural Engineering and Mechanics
• /
• v.48 no.4
• /
• pp.501-518
• /
• 2013

Apart from thinning of cortical layers, the local bone curvature, varying along bone periphery, modulates ultrasound waves as well, which is however often underestimated or overlooked in clinical quantitative ultrasound (QUS). A dedicated three-dimensional finite element modelling technique for cortical bones was established, for quantitatively exploring and calibrating the effect of local curvature of cortical bone on ultrasound. Using a correlation-based mode extraction technique, high-velocity group (HVG) and low-velocity group (LVG) wave modes in a human radius were examined. Experimental verification using acrylic cylinders and in vitro testing using a porcine femur were accomplished. Results coherently unravelled the cortical curvature exerts evident influence on bone-guided ultrasound when RoC/${\lambda}$ <1 for HVG mode and RoC/${\lambda}$ <2 for LVG mode (RoC/${\lambda}$: the ratio of local bone curvature radius to wavelength); the sensitivity of LVG mode to bone curvature is higher than HVG mode. It has also been demonstrated the local group velocity of an HVG or LVG mode at a particular skeletal site is equivalent to the velocity when propagating in a uniform cylinder having an outer radius identical to the radius of curvature at that site. This study provides a rule of thumb to compensate for the effect of bone curvature in QUS.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.311-319
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate, hollow cylinder and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.55-60
• /
• 2011

In the present work, isogeometric analysis in linear elasticity problem is conducted using the basis functions from NURBS. The objectives of isogeometric analysis introduced is to integrate both geometric modeling(CAD) and computational analysis(CAE), and this can be accomplished from direct usage of geometric modeling by NURBS as the computational mesh. The merit of the isogeometry analysis is that NURBS surface are able to represent exact geometry from the control points and knot vectors, and also subsequent refinement is relatively simple relatively. In order to verify the computer codes developed in this study, it has been applied to two structural models of which geometry are simple ; 1) circular cylinder subjected to the constant internal pressure loading, 2) square plate with circular hole at center subjected to uniform tension. The exact solutions of these two models are available. Convergence of the approximate solutions by the present code for the isogeometry analysis are investigated by mesh refinement with inserting knots (h-refinement) and by mesh refinement with order elevation of the basis functions (p-refinement).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.671-680
• /
• 2007

The objective of this work is to reduce drag on a bluff body within a viscous flow by applying suction or injection of fluid along the surface of the body. In addition to minimizing drag, the optimal solution tends to reduce boundary layer separation and flow recirculation. When discretized by finite elements, the optimal control problem can be posed as a large-scale nonlinearly-constrained optimization problem. The constraints correspond to the discretized form of the Navier-Stokes equations. Unfortunately, solving such large-scale problems directly is essentially intractable. We developed several Sequential Quadratic Programming methods that are tailored to the structure of the control problem. Example problems of laminar flow around an infinite cylinder in two dimensions are solved to demonstrate the methodology. We use these optimal control techniques to study the influence of number of suction/injection holes and location of holes on the resulting optimized flow. We compare the proposed SQP methods against one another, as well as against available methods from the literature, from the point of view of efficiency and robustness. The most efficient of the proposed methods is two orders of magnitude faster than existing methods.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.7
• /
• pp.629-636
• /
• 2013

Pneumatic levitation is based upon Bernoulli's principle. However, this method is known to require a large gas flow rate that can lead to an increase in the cost of products. In this case, the gas flow rate should be increased, and the compressible effects of the gas may be of practical importance. In the present study, a computational fluid dynamics method has been used to obtain insights into Bernoulli levitation flows. Three-dimensional compressible Navier-Stokes equations in combination with the SST k-${\omega}$ turbulence model were solved using a fully implicit finite volume scheme. The gas flow rate, workpiece diameter,and clearance gap between the workpiece and the circular cylinder were varied to investigate the flow characteristics inside. It is known that there is an optimal clearance gap for the lifting force and that increasing the supply gas flow rate results in a larger lifting force.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.22 no.3_1spc
• /
• pp.573-579
• /
• 2013

Once-through helically coiled steam generator tubes subjected to external pressure are of interest because of their application to advanced small- and medium-sized integral reactors, in which a primary coolant with a relatively higher pressure flows outside the tubes, while secondary water with a relatively lower pressure flows inside the tubes. Another notable point is that the values of the mean radius to thickness ratio of these steam generator tubes are very small, which means that a thick-walled cylinder is employed for these steam generator tubes. In the present paper, the maximum allowable pressure of helically coiled and thick-walled steam generator tubes with through-wall cracks under external pressure is investigated based on a detailed nonlinear three-dimensional finite element analysis. In terms of the crack orientation, either circumferential or axial through-wall cracks are considered. In particular, in order to quantify the effect of the crack location on the maximum external pressure, these cracks are assumed to be located in the intrados, extrados, and flank of helically coiled cylinders. Moreover, an evaluation is also made of how the maximum external pressure is affected by the ovality, which might be inherently induced during the tube coiling process used to fabricate the helically coiled steam generator tubes.

• Design & Manufacturing
• /
• v.10 no.1
• /
• pp.51-54
• /
• 2016

Metal sheet forming has been commonly used as the core technology in manufacturing parts of automobiles. It guarantees the highest production rate due to the process of mass production employing the press die. For precision of the product, the accuracy of the molds and its mechanic structures are considered as essential factors. One of these is fine blanking, which is utilized for the production of the metal sheet spring, with which clear sheer surfaces can be achieved in one operation from the materials. However, the current designs of press dies perform the forming analysis with the molds of rigid body, so they are focused on weight lightening by a rule of thumb. Therefore, this paper practice structural analysis about developing the semi fine-fine blanking technology. The semi fine-blanking can be run through the combination of the hydraulic cylinders and normal presses, so this paper analyze the amount of deformation according to the oil pressure. In addition, based on the plasticity of 50CrV4, the materials of the mold parts, the structural analysis and life analysis are proceeded, so they are expected to be useful as data for manufacturing the mold.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Journal of Ocean Engineering and Technology
• /
• v.19 no.5 s.66
• /
• pp.16-25
• /
• 2005

The effect of vertical riser support system on the dynamic behaviour of a classical spar platform is investigated. Spar platform generally uses buoyancy-can riser support system, but as water depth gets deeper the alternative riser support system is required due to safety and cost issues. The alternative riser support system is to hang risers off the spar platform using pneumatic cylinders rather than the buoyancy-can. The existing numerical model for hull/mooring/riser coupled dynamics analysis treats riser as an elastic rod truncated at the keel (truncated riser model), thus, in this model, the effect of riser support system can not be modeled correctly. Due to this reason, the truncated riser model tends to overestimate the spar pitch and heave motion. To evaluate more realistic global spar motion, mechanical coupling among risers, guide frames and support cylinders inside of spar moon-pool should be modeled. In the newly developed model, the risers are extended through the moon-pool by using nonlinear finite element methods with realistic boundary condition at multiple guide frames. In the simulation, the vertical tension from pneumatic cylinders is modeled by using ideal-gas equation and the vertical tension from buoyancy-cans is modeled as constant top tension. The different dynamic characteristics between buoyancy-can riser support system and pneumatic riser support system are extensively studied. The alternative riser support system tends to increase spar heave motion and needs damper system to reduce the spar heave motion.