• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.859-865
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered bearing and flexible supporting structures by using the finite element method and the mode superposition method. The appropriate finite element equations for polygon mirror are described by rotating annular sector element using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. The rotating components except for the polygon mirror are modeled by Timoshenko beam element including the gyroscopic effect. The flexible supporting structures are modeled by using a 4-node tetrahedron element and 4-node shell element with rotational degrees of freedom. Finite element equations of each component of the polygon mirror scanner motor and the flexible supporting structures are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. The rigid link constraints are also imposed at the interface area between sleeve and sintered bearing to describe the physical motion at this interface. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by solving the associated eigenvalue problem by using the restarted Arnoldi iteration method. Unbalance responses in time and frequency domain are performed by superposing the eigenvalues and eigenvectors from the free vibration analysis. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results. This research also shows that the flexibility of supporting structures plays an important role in determining the unbalance response of the polygon mirror scanner motor.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.5
• /
• pp.575-584
• /
• 2014

In this study, the NREL 5 MW wind turbine tower model was optimized according to the multi-body dynamics and reliability-based design. The mathematical model was defined as a link-joint system including dynamic characteristics derived from Timoshenko's beam theory. For the optimization problem, the sensitivities to variations in the tower thicknesses and inner and outer diameters were acquired and arranged in terms of safety and efficiency according to bending stress and buckling standards. An optimal design was calculated with the advanced first-order second moment method and used to define a finite element model for validation. The finite element model was simulated by static analysis. The relationship between the multi-body dynamic and finite element method throughout the process was investigated, and the optimal model, which had high endurance despite its low mass, was determined.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.6
• /
• pp.447-455
• /
• 2016

Precise propulsion shafting alignment of ships is very important to prevent damage of its support bearings due to excessive reaction forces caused by hull deflection, forces acted on propeller and crankshaft, and so forth. In this paper, a new iterative shafting alignment calculation procedure considering the interaction between shaft deflection and oil film pressure of Sterntube Journal Bearing (SJB) bush with single or multiple slopes is proposed. The procedure is based on a pressure analysis to evaluate distributed equivalent support stiffness of SJB by solving Reynolds equation and a deflection analysis of shafting system by a finite element method based on Timoshenko beam theory. SJB is approximated with multi-point biaxial elastic supports equally distributed to its length. Their initial stiffness values are estimated from dynamic reaction force calculated by assuming SJB as single rigid support. Then, the shaft deflection and the support stiffness of SJB are sequentially and iteratively calculated by applying a criteria on deflection variation between sequential calculation results. To demonstrate validity and applicability of the proposed procedure for optimal slope design of SJB, numerical analysis results for a shafting system are described.

• Structural Monitoring and Maintenance
• /
• v.8 no.4
• /
• pp.345-360
• /
• 2021

This paper presents an analytical approach to calculate the buckling load of the cylindrical ring structures subjected to both hydrostatic and hydrodynamic pressures. Based on the conservative law of energy and Timoshenko beam theory, a theoretical formula, which can be used to evaluate the critical pressure of buckling, is first derived for the simplified cylindrical ring structures. It is assumed that the hydrodynamic pressure can be treated as an equivalent hydrostatic pressure as a cosine function along the perimeter while the thickness ratio is limited to 0.2. Note that this paper limits the deformed shape of the cylindrical ring structures to an elliptical shape. The proposed analytical solutions are then compared with the numerical simulations. The critical pressure is evaluated in this study considering two possible failure modes: ultimate failure and buckling failure. The results show that the proposed analytical solutions can correctly predict the critical pressure for both failure modes. However, it is not recommended to be used when the hydrostatic pressure is low or medium (less than 80% of the critical pressure) as the analytical solutions underestimate the critical pressure especially when the ultimate failure mode occurs. This implies that the proposed solutions can still be used properly when the subsea vehicles are located in the deep parts of the ocean where the hydrostatic pressure is high. The finding will further help improve the geometric design of subsea vehicles against both hydrostatic and hydrodynamic pressures to enhance its strength and stability when it moves underwater. It will also help to control the speed of the subsea vehicles especially they move close to the sea bottom to prevent a catastrophic failure.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.6
• /
• pp.400-406
• /
• 2021

The primary effect of the far-field underwater explosion (UNDEX) is the whipping of the ship hull girder. This paper aims to verify why inelastic effects should be considered in the whipping response estimations from the UNDEX simulations. A navy ship was modeled using Timoshenko beam elements over the ship length uniformly keeping the constant midship section modulus. The transient UNDEX pressure was produced using two types of the Geers-Hunter doubly-asymptotic models: compressible and incompressible fluids. Because the UNDEX model based on incompressible fluid assumption provided more increased fluid volume acceleration in the bubble phase, the incompressible fluid-based UNDEX model was adopted for the inelastic whipping response analyses. The non-linear hull girder bending moment-curvature curve was used to embed inelastic effects in the UNDEX analyses where the Smith method was applied to derive the non-linear stiffness. We assumed two stand-off distances to see more apparent inelastic effects: 40.5 m and 35.5 m. In the case of the 35.5 m stand-off distance, there was a statistically significant inelastic effect in terms of the average of peak moments and the average exceeding proportional limit moments. For the conservative design of a naval ship under UNDEX, it is recommended to use incompressible fluid. In the viewpoint of cost-effective naval ship design, the inelastic effects should be taken into account.

• Tribology and Lubricants
• /
• v.40 no.1
• /
• pp.8-16
• /
• 2024

This study conducts numerical modeling and eigen-analysis of a rod-fastened rotor, which is mainly used in aircraft gas turbine engines in which multiple disks are in contact through curvic coupling. Nayak's theory is adopted to calculate surface parameters measured from the tooth profile of the curvic coupling gear. Surface parameters are important design parameters for predicting the stiffness between contact surfaces. Based on the calculated surface parameters, elastoplastic contact analysis is performed according to the interference between two surfaces based on the Greenwood-Williamson model. The equivalent bending stiffness is predicted based on the shape and elastoplastic contact stiffness of the curvic coupling. An equation of motion of the rod-fastened rotor, including the bending stiffness of the curvic coupling, is developed. Methods for applying the bending stiffness of a curvic coupling to the equation of motion and for modeling the equation of motion of a rotor that includes both inner and outer rotors are introduced. Rotordynamic analysis is performed through one-dimensional finite element analysis, and each element is modeled based on Timoshenko beam theory. Changes in bending stiffness and the resultant critical speed change in accordance with the rod fastening force are predicted, and the corresponding mode shapes are analyzed.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.155-155
• /
• 2011

일반적으로 I형 보에 횡하중이 작용하는 경우, 횡 변위와 함께 회전을 동반하는 횡-비틀림 좌굴(Lateral-Torsional Buckling)이 발생하게 된다. 이러한 I형 보의 탄성 및 비탄성 횡-비틀림 좌굴에 대한 해석적 이론적 연구는 이미 많은 연구자들에 의해 수행되었다(Timoshenko 등, 1961; Galambos, 1963; Lindner, 1974; Trahair, 1993). I형 보의 비지지 길이 내 하중이 작용할 때 모멘트 구배계수(Cb)는 하중이 부재 단면에 작용하는 위치에 따라 달라지게 되는데 이를 하중고 효과(Load Height Effects)라고 한다. 탄성 영역 내 비지지길이가 존재하는 I형 보의 하중고 효과를 고려한 모멘트 구배계수 제안식은 Nethercot & Rockey(1971)에 의해 연구된 바 있다. 또한 Helwig 등(1997)은 Nethercot & Rockey(1971)의 제안식을 간략화 하여 탄성 영역 내 비지지길이가 존재하는 I형 보의 하중고 효과를 고려한 모멘트 구배계수식을 제안하였다. 그러나 현재까지 진행 된 하중고 효과에 대한 연구는 탄성 영역 내 비지지 길이가 존재하는 I형 보에 대한 제안식이며 현재까지 비탄성 영역 내 비지지 길이를 갖는 I형 보의 하중고 효과에 대한 연구는 진행된 바 없다. 본 연구는 비탄성 영역 내 비지지 길이가 존재하는 I형 보의 하중고 효과를 고려한 비탄성 횡-비틀림 좌굴강도에 대한 연구를 수행하였다. 하중조건으로는 집중하중 과 등분포 하중을 적용시켰으며, 비선형 횡-비틀림 좌굴 해석을 위해 잔류응력 및 초기변형을 고려하였다. Pi와 Trahair(1995)이 고려한 단순직선분포를 잔류응력으로 가정하였으며, 국내 I형강 표준 치수 허용치(현대제철, 2006)에 근거하여 부재 길이의 0.1%를 초기 최대 횡 변위로 적용하여 초기제작오차로 고려하였다. 유한요소해석결과를 바탕으로 Nethercot & Rockey(1971)와 Helwig 등(1997)의 연구내용을 바탕으로 범용구조해석 프로그램(ABAQUS, 2007)을 이용하여 비탄성 영역 내 존재하는 I형보의 횡-비틀림 좌굴강도를 산정하였다. 유한요소해석결과를 바탕으로 Nethercot & Rockey(1971)및 Helwig 등(1997)의 모멘트구배계수 제안식과 비교 분석 하였고 회기분석프로그램 MINITAB(2006)을 이용하여 비탄성 영역 내 비지지길이가 존재하는 I형보의 하중고 효과를 고려한 모멘트구배계수식을 개발 제안하였다. 본 연구에서 개발된 제안식은 경제적이고 합리적인 휨부재 강도평가에 적극 활용될 수 있으며, 비탄성 영역내 I형보의 횡-비틀림 좌굴강도 및 휨강도 연구에 널리 활용될 것이다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.10
• /
• pp.773-782
• /
• 2020

The development of joint FE models for deep learning neural network (DLNN)-based hybrid FEA is presented. Material models of bolts and bearings in the front axle of tractor, showing complex behavior induced by various tightening conditions, were replaced with DLNN models. Bolts are modeled as one-dimensional Timoshenko beam elements with six degrees of freedom, and bearings as three-dimensional solid elements. Stress-strain data were extracted from all elements after finite element analysis subjected to various load conditions, and DLNN for bolts and bearing were trained with Tensorflow. The DLNN-based joint models were implemented in the ABAQUS user subroutines where stresses from the next increment are updated and the algorithmic tangent stiffness matrix is calculated. Generalization of the trained DLNN in the FE model was verified by subjecting it to a new loading condition. Finally, the DLNN-based FEA for the front axle of the tractor was conducted and the feasibility was verified by comparing with results of a static structural experiment of the actual tractor.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.45-53
• /
• 1985

For the analysis of vertical vibrations of a ship's hull, the Timoshenko beam analogy is accepted up to seven or eight-node modes provided that the system parameters are properly calculated. As to the shear coefficient, it has been a common practice to apply the strain energy method or the projected area method. The theoretical objection to the former is that it ignores lateral contraction due to Poisson's ratio, and the latter is of extreme simplifications. Recently, Cowper's and Stephen's shear coefficient formulas have drawn ship vibration analysts' attentions because these formulas, derivation of which are based on an integrations of the equations of three-dimensional elasticity, take Poisson's ratio into account. Providing computer programs for calculation of the shear coefficient of ship sections modeled as thin-walked multicell sections by each of the forementioned methods, the authors calculated natural vibration characteristics of a bulk carrier and of a container ship by the transfer matrix method using shear coefficients obtained by each of the methods, and discussed the results in comparision. The major conclusions resulted from this investigation are as follows: (1) The shear coefficients taking account of the effects of Poisson's ratio, Cowper's $K_c$ and Stephen's $K_s$, result in higher values of about 10% in maximum as compared with the shear coefficient $K_o$ based on the conventional strain energy methods; (a) $K_c/K_o{\cong}1.05\;and\;K_s/K_o{\cong}1.10$ for ships having single skin side-shell such as a bulk carrier. (b) $K_c/K_o{\cong}1.02\;and\;K_s/K_o{\cong}1.05$ for ships having longitudinally through bulkheads and/or double side-shells in the portion of the cargo hod such as a container carrier. (2) The distributions of the effective shear area along the ship's hull based on each of $K_o,\;K_c\;and\;K_s$ are similar each another except the both end portions. (3) Natural frequencies and mode shapes of the hull based on each of $K_c\;and\;K_s$ are of small differences as compared each other. (4) In cases of using $K_c\;or\;K_s$ in ship vibration analysis, it is also desirable to have the bending rigidity be corrected according to the effective breadth concept. And then, natural frequencies and mode shapes calculated with the bending rigidity corrected in the above and with each of $K_o,\;K_c\;and\;K_s$ result in small differences as compared each another. (5) Referring to those mentioned in the above (3) and (4) and to the full-scale experimental results reported by Asmussen et al.[17], and considering laboursome to prepare the computer input data, the following suggestions can safely be made; (a) Use of $K_o$ in ship vibration analysis is appropriate in practical senses. (b) Use of $K_c$ is appropriate even for detailed vibration analysis of a ship's hull. (6) The effective shear area based on the projected area method is acceptable for the two-node mode.