• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.635-647
• /
• 2014

In this paper, a generalized perturbed midpoint rule is established. Various error bounds for this generalization are also obtained.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.11
• /
• pp.1907-1916
• /
• 2003

During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

• Bulletin of the Korean Mathematical Society
• /
• v.42 no.2
• /
• pp.213-230
• /
• 2005

Companions of Ostrowski's integral inequality for absolutely continuous functions and applications for composite quadrature rules and for p.d.J.'s are provided.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.2
• /
• pp.221-229
• /
• 2010

Some new $\check{C}$eby$\check{s}$ev type inequalities have been developed by working on functions whose first derivatives are absolutely continuous and the second derivatives belong to the usual Lebesgue space $L_{\infty}[a,\;b]$. A unified treatment of the special cases is also given.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.99-104
• /
• 2006

A semi-Lagrangian finite element scheme with objective time stepping algorithm for solving viscoelastic flow problem is presented. The convection terms in the momentum and constitutive equations are treated using a quasi-monotone semi-Lagrangian scheme, in which characteristic feet on a regular grid are traced backwards over a single time-step. Concerned with the generalized midpoint rule type of algorithms formulated to exactly preserve objectivity, we use the geometric transformation such as pull-back, push-forward operation. The method is applied to the 4:1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. 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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• The Journal of Churna Manual Medicine for Spine and Nerves
• /
• v.14 no.2
• /
• pp.135-146
• /
• 2019

Objectives : There are many theories about the anatomical relationships between the thoracic spinous and transverse processes. However, these studies were all conducted on cadavers. Thus, there might be differences when applied to living individuals. Therefore, the aim of this study was to determine whether the theories were similar when measured in living individuals. Methods : We conducted studies with 10 individuals aged 20-30 using 3D-CT. Four different measurements were taken between the thoracic spinous and transverse processes. To facilitate the relationship analysis, the spinous process positioned between the traverse process of the segment and the lower segment was set to 0.5. The spinous process located at the height of the traverse process of the lower segment was set to 1. The spinous process located between the traverse process of the lower segment and the two segments below was set to 1.5. Results : Therefore, based on the aforementioned settings, T1-T4 were 1, T5-T7 were 1.5, T8-T10 were 1, and T11 and T12 were 0.5. This indicated that the spinous processes of T1-T4 and T8-T10 were equal in height to the transverse processes of the lower segment, and that T5-T7 were at the midpoint of the height of the transverse processes between the lower segment and two segments below, and T11, T12 corresponded to the midpoint between the traverse process of the segment and the lower segment. Conclusions : Neither 'Geelhoed's rule'nor 'the rule of three'was accurate when applied to living individuals. However, this study had some limitations, so further research is needed.