• Ocean Systems Engineering
• /
• 제4권4호
• /
• pp.327-342
• /
• 2014

The tension leg platform (TLP) is one of the compliant structures which are generally used for deep water oil exploration. With respect to the horizontal degrees of freedom, it behaves like a floating structure moored by vertical tethers which are pretension due to the excess buoyancy of the platform, whereas with respect to the vertical degrees of freedom, it is stiff and resembles a fixed structure and is not allowed to float freely. In the current study, a numerical study for square TLP using modified Morison equation was carried out in the time domain with water particle kinematics using Airy's linear wave theory to investigate the effect of changing the tether tension force on the stiffness matrix of TLP's, the dynamic behavior of TLP's; and on the fatigue stresses in the cables. The effect was investigated for different parameters of the hydrodynamic forces such as wave periods, and wave heights. The numerical study takes into consideration the effect of coupling between various degrees of freedom. The stiffness of the TLP was derived from a combination of hydrostatic restoring forces and restoring forces due to cables. Nonlinear equation was solved using Newmark's beta integration method. Only uni-directional waves in the surge direction was considered in the analysis. It was found that for short wave periods (i.e., 10 sec.), the surge response consisted of small amplitude oscillations about a displaced position that is significantly dependent on tether tension force, wave height; whereas for longer wave periods, the surge response showed high amplitude oscillations that is significantly dependent on wave height, and that special attention should be given to tethers fatigue because of their high tensile static and dynamic stress.

• 한국소음진동공학회논문집
• /
• 제23권3호
• /
• pp.240-248
• /
• 2013

자동차 수동변속기 내 래틀과 같은 회전체 진동은 엔진 연소폭발 주기와 연동하여 발생한다. 이러한 충격형 진동 특성은 선형 시간 불변 분석법을 적용하여 살펴볼 수 있다. 기어의 동적 특성을 이해하기 위해, 특정 형태의 앞 바퀴 굴림 방식 수동 변속기에 대한 연구를 수행하였다. 첫째, 동특성 고유값과 주파수 응답함수를 기반으로 자유도를 줄이는 방법을 제시하였는데, 이는 행렬의 크기를 줄이는 효과적인 역할을 한다. 둘째, 단일 질량 플라이휠과 이중 질량 플라이휠의 동적 특성을 비교하였다. 모빌리티 분석을 기반으로 한 이중 질량 플라이휠의 효과를 검토하였는데, 이로부터 진동에 의해 가진되는 충격을 회피하기 위한 기본 개념을 이해 할 수 있다. 마지막으로, 두 단계 클러치 댐퍼로부터 유효 강성 값을 도출하여 선형 불변형 시스템 모델을 연구하였다. 두 단계의 서로 다른 클러치 강성조합을 이용하여 클러치 댐퍼의 동적특성에 대한 관계를 예측할 수 있다.

• 콘크리트학회논문집
• /
• 제25권5호
• /
• pp.529-537
• /
• 2013

이 연구의 목적은 하중기반 유한요소 정식화에 의한 FRP 보강된 철근콘크리트 보 또는 기둥 부재의 비선형 층상화의 등매개 골조 유한요소모델을 개발하는데 있다. 단면에서 콘크리트는 3축 응력-변형률 관계로 모델화하고 FRP 피복층은 2차원의 적층복합재료로 모델화하였다. 하중기반 유한요소의 요소강성행렬은 변위형상함수의 가정이 없고 하중보간함수를 갖고 있다. 횡 하중을 받는 GFRP 시트 보강된 철근콘크리트 기둥의 실험에 대해 개발된 하중기반 유한요소모델에 의한 해석을 수행하였다. 기존 강성도법의 유한요소해석과 비교하여 하중기반 유한요소해석은 전체적인 하중-변위 관계 뿐만 아니라 기둥의 소성힌지영역에서의 비선형 변형 및 손상을 보다 정확히 예측해 주었다.

• Structural Engineering and Mechanics
• /
• 제13권5호
• /
• pp.557-568
• /
• 2002

For the evaluation of the capability of a tubular member of an offshore structure to absorb the collision energy, a simple method can be employed for the collision analysis without performing the detailed analysis. The most common simple method is the rigid-plastic method. However, in this method any characteristics for horizontal movement and rotation at the ends of the corresponding tubular member are not included. In a real structural system of an offshore structure, tubular members sustain a certain degree of elastic support from the adjacent structure. End fixity has influences in the behaviors of a tubular member. Three-dimensional FEM analysis can include the effect of end fixity fully, however in viewpoints of the inherent computational complexities of the 3-D approach, this is not the recommendable analysis at the initial design stage. In this paper, influence of end fixity on the behaviors of a tubular member is investigated, through a new approach and other approaches. A new analysis approach that includes the flexibility of the boundary points of the member is developed here. The flexibility at the ends of a tubular element is extracted using the rational reduction of the modeling characteristics. The property reduction is based on the static condensation of the related global stiffness matrix of a model to end nodal points of the tubular element. The load-displacement relation at the collision point of the tubular member with and without the end flexibility is obtained and compared. The new method lies between the rigid-plastic method and the 3-demensional analysis. It is self-evident that the rigid-plastic method gives high strengthening membrane effect of the member during global deformation, resulting in a steeper slope than the present method. On the while, full 3-D analysis gives less strengthening membrane effect on the member, resulting in a slow going load-displacement curve. Comparison of the load-displacement curves by the new approach with those by conventional methods gives the figures of the influence of end fixity on post-yielding behaviors of the relevant tubular member. One of the main contributions of this investigation is the development of an analytical rational procedure to figure out the post-yielding behaviors of a tubular member in offshore structures.

• 대한조선학회지
• /
• 제22권1호
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• 대한조선학회논문집
• /
• 제29권4호
• /
• pp.202-213
• /
• 1992

선체(船體)의 횡강도(橫强度) 부재(部材)를 설계(設計)하기 위해서는 필수적(必須的)으로 횡강도(橫强度) 해석(解析)을 수반(隨伴)하여야 하며 이에 따라 많은 작업(作業) 시간(時間) 및 계산(計算) 시간(時間)이 필요(必要)하게 된다. 선체(船體)의 횡강도(橫强度) 해석(解析)을 위해 종래(從來)에는 경사(傾斜) 처짐법을 이용(利用)한 해석적(解析的) 방법(方法)이 사용(使用)되어 왔지만 부재(部材)의 신축(伸縮)을 무시(無視)함에 따라 해(解)의 정도(精度)가 상당히 낮을 뿐아니라 층방정식(層方程式)을 표현(表現)하기가 까다로워 프로그램 구성(構成)이 어렵다. 또한 최근(最近) computer의 발달(發達)과 함께 급속도(急速度)로 발전(發展)하고 있는 유한요소법(有限要素法)을 이용(利用)하여 선체(船體)의 횡강도(橫强度) 해석(解析)을 수행(遂行)하고 있지만 아직도 구조(構造) 설계(設計) 및 최적(最適) 구조(構造) 설계(設計)를 수행(遂行)하기에는 계산(計算) 시간(時間)의 극복(克服)이 어려운 실정(實情)이다. 본(本) 연구(硏究)에서는 선체(船體)의 구조(構造) 해석(解析) 및 설계(設計)를 위해 bracket이 붙은 부분을 span point 개념(槪念)을 도입(導入)하여 처리(處理)하고, 기존(旣存)의 경사(傾斜) 처짐법에서 무시(無視)하였던 부재(部材)의 신축(伸縮)에 따른 축방향(軸方向)의 변위(變位)를 고려(考慮)하여 각 절점(節点)에서의 평형방정식(平衡方程式)으로만 해(解)를 구할 수 있도록 하고 matrix method와 결합(結合)하여 2차원(次元) 및 3차원(次元)에 대한 일반화(一般化) 경사(傾斜) 처짐법을 유도(誘導)하였으며 기존(旣存)의 경사(傾斜) 처짐법 및 유한요소법(有限要素法)과 계산(計算) 시간(時間) 및 정도(精度)를 비교하여 본(本) 해석(解析) 방법(方法)의 우수성(優秀性)을 입증(立證)하였다.

• Structural Engineering and Mechanics
• /
• 제81권4호
• /
• pp.411-428
• /
• 2022

This paper presents the Campbell diagram analysis of the rotordynamic system using the full order model (FOM) and the reduced order model (ROM) techniques to determine the critical speeds, identify the stability and reduce the computational time. Due to the spin-speed-dependent matrices (e.g., centrifugal stiffening matrix), several model order reduction (MOR) techniques may be considered, such as the modal superposition (MS) method and the Krylov subspace-based MOR techniques (e.g., Ritz vector (RV), quasi-static Ritz vector (QSRV), multifrequency quasi-static Ritz vector (MQSRV), multifrequency/ multi-spin-speed quasi-static Ritz vector (MMQSRV) and the combined Ritz vector & modal superposition (RV+MS) methods). The proposed MMQSRV method in this study is extended from the MQSRV method by incorporating the rotational-speed-dependent stiffness matrices into the Krylov subspace during the MOR process. Thus, the objective of this note is to respond to the question of whether to use the MS method or the Krylov subspace-based MOR technique in establishing the Campbell diagram of the shaft-disc-blade assembly systems in three-dimensional (3D) finite element analysis (FEA). The Campbell diagrams produced by the FOM and various MOR methods are presented and discussed thoroughly by computing the norm of relative errors (ER). It is found that the RV and the MS methods are dominant at low and high rotating speeds, respectively. More precisely, as the spinning velocity becomes large, the calculated ER produced by the RV method is significantly increased; in contrast, the ER produced by the MS method is smaller and more consistent. From a computational point of view, the MORs have substantially reduced the time computing considerably compared to the FOM. Additionally, the verification of the 3D FE rotordynamic model is also provided and found to be in close agreement with the existing solutions.

• Structural Engineering and Mechanics
• /
• 제12권6호
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.

• 대한토목학회논문집
• /
• 제14권6호
• /
• pp.1299-1308
• /
• 1994

본 연구는 사장교의 차량하중에 의한 동적거동을 파악하고자 수치해석상 비교적 간단한 일본의 풍리(豊里)(Toyosato)대교(大橋)의 자료를 근거로하여 수치해석 대상모델을 구조형식별로, 여러가지 설계변수-즉, (1) 경간비, (2) 중앙경간장과 주탑높이와의 비, (3) 거어더의 강성, (4) 주탑의 강성, (5) 케이블의 강성-을 변화시켜 수치해석을 수행하여 동적거동을 파악하고, 그 결과를 가지고 설계변수의 영향 및 충격계수의 변화에 대하여 비교 분석하였다. 이때 변위 및 단면력의 영향선을 구하기 위한 해석은 전달행렬법을 이용하였으며, 동적해석에 있어서는 평면구조계의 집중질량계로 모델을 가정하여 차량과 교량의 운동방정식을 유도한 후 모드중첩법을 사용하여 각 질점에 대한 변위 및 단면력의 동적시간이력을 구하였다.

• 한국전산구조공학회논문집
• /
• 제23권5호
• /
• pp.475-482
• /
• 2010

본 논문에서는 대규모 자유도 시스템의 병렬처리를 위하여 2단계로 이루어진 영역분할법(Domain Decomposition Method) 기반의 병렬 알고리즘을 제안하였다. 분할된 영역의 내부 및 외부 경계를 상위영역문제로 정의하고 국부영역문제는 변위 경계조건이 모두 주어지는 분할영역에서의 Dirichlet 문제로 구성한다. 상위영역에서는 전체 상위영역에 대한 강성 행렬의 어셈블이 필요없는 반복법을 통하여 변위를 구하고, 이를 바탕으로 국부영역에서 Multi-Frontal Sparse Solver (MFSS)를 이용하여 변위를 계산한다. 상위영역문제의 연산에서 프로세서 간의 데이터 교환을 최소화하여 계산효율을 유지하며, 동시에 해석 가능한 자유도를 증대시키는 병렬 PCG(Preconditioned Conjugate Gradient)법 기반의 알고리즘을 개발하였다. 제안된 알고리즘을 적용하여 수치해석을 수행한 결과, 프로세서 수가 증가할수록 계산성능의 손실없이 해석 가능한 자유도가 비례하여 증가하는 선형 확장성을 관찰할 수 있었으며, 대규모 자유도 문제에 효과적으로 사용 가능함을 확인하였다.