• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.7
• /
• pp.989-996
• /
• 2001

High-speed and low-speed flows are simulated numerically by flowfield-dependent mixed explicit-implicit (FDMEI) method. This algorithm depends on implicitness parameters of convection, diffusion, diffusion gradients, and source terms which are calculated from the changes of local Mach, Reynolds, Peclet, and Damkohler numbers between adjacent nodes. Convection phenomena or shock waves are resolved from Mach number-dependent implicitness parameters whereas diffusion or viscous actions are simulated by Reynolds number or Peclet number-dependent implicitness parameters. Fluctuation components of all variables are properly accommodated spatially and temporally in the FDMEI procedure. To illustrate, some benchmark example problems are presented for comparisons of the FDMEI results with other available data. These results appear to be encouraging and point toward the need for further investigations of the FDMEI theory.

• Proceedings of the KSR Conference
• /
• 2007.11a
• /
• pp.1605-1609
• /
• 2007

In the case of Linear induction motor(LIM), numerical analysis method like Finite Element Method(FEM) has been mainly used to analyze the travelling magnetic field problem which includes the velocity-induced electromotive force. If the problem including the velocity-induced electromotive force is analyzed by FEM using Galerkin method, the solution can be oscillated according to Peclet Number, which is determined by conductivity, permeability, moving velocity and size of mesh. Consequently, the accuracy of the solution can be low and the Loops of flux can be occurred at the secondary back-iron. These loops of flux occurred at the secondary back-iron does not exist physically, but it can be occurred in the analysis. In this case, the loops of flux can be generally removed by using Up-Wind method which is impossible to apply a conventional S/W tool(Maxwell 2D). Therefore, in this paper, authors examined the Loops of flux occurred at the secondary back-iron of LIM according to variations of Peclet Number, and analyzed whether these loops of flux affect on the thrust force characteristics of LIM or not.

• Wind and Structures
• /
• v.29 no.6
• /
• pp.441-455
• /
• 2019

In this study, vortex induced vibrations of a cylinder mounted on a flexible rod are analyzed. This simple configuration represents the key element of new conception bladeless wind turbine (Whitlock 2015). In this study the structure oscillations equation coupled to the wake oscillation equation for this configuration are solved using analytical perturbation method, for the first time. An analytical expression that predicts the lock-in phenomena range of wind speed is derived. The discretized equations of motion are also solved using RKF45 numerical method. The equations of motion are discretized by Galerkin method. Free vibration mode shape of the structure taking into account the discontinuity of the cross section are used as comparison function. Numerical results are compared to the analytical results, and they show a satisfying agreement. The effect of system parameters on the oscillations of structure and wake as well as on the lock-in domain are presented. Moreover, it is shown that the values of wind speed triggering the start and the stop of the lock-in phenomenon, for increasing wind speed are different from those values obtained during the reverse process, i.e., when the wind speed decreases.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.24 no.11
• /
• pp.898-906
• /
• 2014

This paper presents a reduction method for a straight-line path error of a flexible beam deploying from a rotating rigid hub. Previous studies discussed about only vibration phenomena of flexible beams deploying from rotating hubs; however, this study investigates a vibration reduction of a rotating beam with variable length. The equation of motion and associated boundary conditions are derived for a flexible beam deploying from a rotating rigid hub, and then they are transformed to a variational equation. By applying the Galerkin method, the discretized equations are obtained from the variational equation. Based on the discretized equations, the dynamic responses of a rotating/deploying beam are analyzed when the beam end has a straight line motion. A reduction method for the trajectory error is proposed, using the average length of a rotating/deploying beam. It is shown that the proposed method is able to reduce the residual vibration of a rotating/deploying beam.

• Steel and Composite Structures
• /
• v.32 no.2
• /
• pp.225-237
• /
• 2019

In this research, the dynamic stability and nonlinear vibration behavior of a smart rotating sandwich cylindrical shell is studied. The core of the structure is a functionally graded material (FGM) which is integrated by functionally graded piezoelectric material (FGPM) layers subjected to electric field. The piezoelectric layers at the inner and outer surfaces used as actuator and sensor, respectively. By applying the energy method and Hamilton's principle, the governing equations of sandwich cylindrical shell derived based on first-order shear deformation theory (FSDT). The Galerkin method is used to discriminate the motion equations and the equations are converted to the form of the ordinary differential equations in terms of time. The perturbation method is employed to find the relation between nonlinear frequency and the amplitude of vibration. The main objective of this research is to determine the nonlinear frequencies and nonlinear vibration control by using sensor and actuator layers. The effects of geometrical parameters, power law index of core, sensor and actuator layers, angular velocity and scale transformation parameter on nonlinear frequency-amplitude response diagram and dynamic stability of sandwich cylindrical shell are investigated. The results of this research can be used to design and vibration control of rotating systems in various industries such as aircraft, biomechanics and automobile manufacturing.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.349-367
• /
• 2022

In this research, an approach combining a semi-analytical method and an analytical method is presented to investigate the static and dynamic post-buckling behavior of the sandwich functionally graded (FG) porous cylindrical shells exposed to external pressure. The sandwich cylindrical shell considered is composed of a viscoelastic core and two FG porous (FGP) face layers. The viscoelastic core is made of Kelvin-Voigt-type material. The material properties of the FG porous face layer are considered continuous through each face thickness according to a porosity coefficient and a volume fraction index. Two types of sandwich FG porous viscoelastic cylindrical shells named Type A and Type B are considered in the research. Type A shell has the porosity evenly distributed across the thickness direction, and Type B has the porosity unevenly distributes across the thickness direction. The FG face layers are considered in two cases: outside metal surface, inside ceramic surface (OMS-ICS), and inside metal surface, outside ceramic surface (IMS-OCS). According to Donnell shell theory, von-Karman equation, and Galerkin's method, a discretized nonlinear governing equation is derived for analyzing the behavior of the shells. The explicit expressions for static and dynamic critical buckling loading are thus developed. To study the dynamic buckling of the shells, the governing equation is examined via a numerical approach implementing the fourth-order Runge-Kutta method. With a procedure presented by Budiansky-Roth, the critical load for dynamic post-buckling is obtained. The effects of various parameters, such as material and geometrical parameters, on the post-buckling behaviors are investigated.

• Structural Engineering and Mechanics
• /
• v.91 no.5
• /
• pp.503-516
• /
• 2024

The present research examines the primary resonance (PR) behaviors of oblique stiffened multilayer functionally graded (OSMFG) shallow shells featuring an FG porous (FGP) core under an external excitation. The research considers two distinct types of FGP cores: one characterized by uniform porosity distribution (UPD) and the other by non-uniform porosity distribution (NPD) along the thickness direction. Furthermore, the study explores two types of shallow shells: one with external oblique stiffeners and one with internal oblique stiffeners, which might have angles that are similar or different from each other. Using the stress function alongside the first-order shear deformation theory (FSDT), the research establishes a nonlinear model for OSMFG shallow shells. The strain-displacement relationships are obtained utilizing FSDT and von-Kármán's geometric assumptions. The Galerkin approach is utilized to discretize the nonlinear governing equations, allowing for the analysis of stiffeners at varied angles. To validate the obtained results, a comparison is made not only with the findings of previous research but also with the response of PR obtained theoretically with the method of multiple scales, using the P-T method. Renowned for its superior accuracy and reliability, the P-T method is deemed an apt selection within this framework. Additionally, the study investigates how differences in material characteristics and stiffener angles affect the system's PR behaviors. The results of this study can be used as standards by engineers and researchers working in this area, and they can offer important information for the design and evaluation of the shell systems under consideration.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.13-23
• /
• 2008

This paper is an extension of previous study[1] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite (serendipity) cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires four degrees-of-freedom at each element corners. Those degrees-of-freedom are the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational basis functions from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[2].

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1897-1907
• /
• 1991

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies, .omega.$_{n}$. Three mode interactions, .omega.$_{2}$=3.omega.$_{1}$, and .omega.$_{3}$=.omega.$_{1}$+2.omega.$_{2}$, are considered and their influence on the response is studied. The case of two mode interaction, .omega.$_{2}$=3.omega.$_{1}$, is also considered in order to compare it with the case of three mode interactions. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous nonlinear ordinary differential equations. The method of multiple scales is applied to obtain steady-state responses of the system. Results of numerical investions show that there exists no significant difference between both modal interactions.

• Journal of the Korean Society for Precision Engineering
• /
• v.13 no.7
• /
• pp.148-157
• /
• 1996

This paper deals with the applicability of linear elastic fracture mechanics under centrifugal force. Stress intensity factors K are calculated as a function of the inclination crack of length 2a, the position at different angular velocities 1200rpm, 2400rpm and at different values of the inclination crack angle .phi. ( .phi. = 0 .deg. , 15 .deg. , 30 .deg. , 45 .deg. , 60 .deg. , 75 .deg. , 90 .deg. ) and are measured in models of rotation disks using a boundary element method. Especially, stress intensity factors $K_{l}$ and $K_{ll}$ obtained separately from the crack tip of the mixed mode, were used to further investigate the influence of $K_{l}$ and $K_{ll}$ on fracture in rotating disks. With the increase in the speed of rotation, the effect of K/ sub l/became larger where as that of $K_{ll}$ became small. For the increase in the inclination crack angle .phi. , a decrease in $K_{l}$ and an increase in $K_{ll}$ were observed.