• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 춘계학술대회논문집
• /
• 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.

• 대한기계학회논문집A
• /
• 제38권5호
• /
• pp.575-584
• /
• 2014

본 연구에서는 NREL 5MW 풍력발전기 타워 설계 모델을 다물체 동역학과 신뢰성 기반 최적 설계를 이용하여 최적화하는 연구를 수행하였다. 타워 모델은 티모센코 빔 이론을 이용하여 얻은 동특성을 내포한 링크와 조인트로 이루어진 수학적 모델로 표현하였다. 최적화 문제에서는 높이가 일정한 타워에서 두께, 내 외곽 지름이 변할 때 나타나는 민감도 변화를 비교하여 결과를 도출하였으며, 비교 기준으로 굽힘 응력과 좌굴 안정성을 사용하였다. 일계 이차 모멘트법을 이용한 최적화 알고리즘에서 얻은 최종 모델은 유한요소법을 이용한 정하중 해석에서 최대 응력 분포를 이용한 안전성을 고려를 통해 유효성을 검증하였다. 본 연구 방법을 통해 동역학적 모델과 유한요소 모델간의 연계성을 확인하고, 낮은 타워 설치 비용으로 더 강건한 시스템을 구축할 수 있는 설계 방향을 제시하였다.

• 대한조선학회논문집
• /
• 제53권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
• /
• 제8권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.

• 대한조선학회논문집
• /
• 제58권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
• /
• 제40권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년도 정기 학술발표대회
• /
• 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형보의 횡-비틀림 좌굴강도 및 휨강도 연구에 널리 활용될 것이다.

• 한국항공우주학회지
• /
• 제48권10호
• /
• pp.773-782
• /
• 2020

심층신경망 기반 하이브리드 유한요소해석을 위한 조인트 모델 방법 구축을 소개한다. 트렉터의 앞차축에서 다양한 체결 조건에 의해 유발되는 복잡한 거동 상태를 가지는 볼트와 베어링의 재료 모델을 심층신경망으로 대체했다. 볼트는 6자유도를 갖는 1차원 티모센코 빔 요소를 이용했고, 베어링은 3차원 솔리드 요소를 이용했다. 다양한 하중 조건을 바탕으로 유한요소해석을 한 뒤, 모든 요소에서 응력-변형률 데이터를 추출하고 텐서플로를 이용하여 학습시켰다. 신경망 기반 유한요소해석을 할 때 추출된 데이터를 바탕으로 학습된 심층신경망은 ABAQUS 서브루틴 안에 포함되어 현재 해석 증분의 응력을 예측하고 접선강도행렬을 계산할 수 있게 했다. 학습된 심층신경망 조인트 모델의 일반화 성능은 훈련에 사용되지 않은 새로운 하중 조건에서 해석하여 검증하였다. 최종적으로 이 방법을 이용하여 심층신경망 기반 앞차축 해석을 진행하고 응력장 분포를 검증했다. 또한, 실제 트렉터의 3점 굽힘 실험 결과와 비교하여 심층신경망 기반 해석의 타당성을 검토했다.

• 대한조선학회지
• /
• 제22권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.