• 대한기계학회논문집B
• /
• 제21권11호
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• 한국지진공학회논문집
• /
• 제1권1호
• /
• pp.1-10
• /
• 1997

최근의 철근 콘크리트 구조물의 내진 설계 방식은 비탄성 거대 변형에 의한 에너지 방출에 의존하고 있다. 이러한 구조물의 거동에 대한 비선형 동적 해석은 특히 계산이 여러 번 반복되어 질 때 많은 시간과 비용이 요구된다. 그러므로 효율적이고 한편 정확한 계산 방법의 채택이 중요하게 되었다. 예측 접근 방법(PASM) 이라 불리는 새로운 방법을 제시하는 것이 현 연구의 주목적이다. 일반적인 동적 해석 방법에서는 매 시간 단계 혹은 반복 계산 때마다 수식계산을 위하여 메트릭스 삼각 분해가 요구되어지나, 예측 접근방법에서는 구조물이 정적 반복하중으로 비선형 범위로 변형되어졌을 때의 강성 상태에서 미리 얻어진 한정적 수의 분해된 메트릭스를 동적 해석에서 이용하게 된다. 이곳에서 제시될 접근 방법은 강성치를 매 시각 단계 혹은 반복 계산 단계마다 재산출해야 하는 다른 접근 방법들과 비교할 때 전체적 수치 해석 양을 줄이게 될 것이다.

• Steel and Composite Structures
• /
• 제30권6호
• /
• pp.535-550
• /
• 2019

Prestressed concrete (PC) bridges using corrugated steel webbing have emerged as one of the most promising forms of steel-concrete composite bridge. However, their long-term behavior is not well understood, especially in the case of large-span bridges. In order to study the time-dependent performance, a large three-span PC bridge with corrugated steel webbing was compared to a similar conventional PC bridge to examine their respective time-dependent characteristics. In addition, a three-dimensional finite element method with step-by-step time integration that takes into account cantilever construction procedures was used to predict long-term behaviors such as deflection, stress distribution and prestressing loss. These predictions were based upon four well-established empirical creep prediction models. PC bridges with a corrugated steel web were observed to have a better long-term performance relative to conventional PC bridges. In particular, it is noted that the pre-cambering for PC bridges with a corrugated steel web could be smaller than that of conventional PC bridges. The ratio of side-to-mid span has great influence on the long-term deformation of PC bridges with a corrugated steel web, and it is suggested that the design value should be between 0.4 and 0.6. However, the different creep prediction models still showed a weak homogeneity, thus, the further experimental research and the development of health monitoring systems are required to further progress our understanding of the long-term behavior of PC bridges with corrugated steel webbing.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2005년도 춘계학술발표대회 논문집
• /
• pp.172-175
• /
• 2005

The numerical simulation is performed for the acoustic emission and the wave propagation due to fiber breakage in single fiber composite plates by the finite element transient analysis. The acoustic emission and the following wave motions from a fiber breakage under a static loading is simulated to investigate the applicability of the explicit finite element method and the equivalent volume force model as a simulation tool of wave propagation and a modeling technique of an acoustic emission. For such a simple case of the damage event under static loading, various parameters affecting the wave motion are investigated for reliable simulations of the impact damage event. The high velocity and the small wave length of the acoustic emission require a refined analysis with dense distribution of the finite element and a small time step. In order to fulfill the requirement for capturing the exact wave propagation and to cover the 3-D simulation, we utilize the parallel FE transient analysis code and the parallel computing technology.

• 대한산업공학회지
• /
• 제16권2호
• /
• pp.135-147
• /
• 1990

A nuclear power plant can be viewed as a large complex man-machine system where high system reliability is obtained by ensuring that sub-systems are designed to operate at a very high level of performance. The chance of severe accident involving at least partial core-melt is very low but once it happens the consequence is very catastrophic. The prediction of risk in low probability, high-risk incidents must be examined in the contest of general engineering knowledge and operational experience. Engineering knowledge forms part of the prior information that must be quantified and then updated by statistical evidence gathered from operational experience. Recently, Bayesian procedures have been used to estimate rate of accident and to predict future risks. The Bayesian procedure has advantages in that it efficiently incorporates experts opinions and, if properly applied, it adaptively updates the model parameters such as the rate or probability of accidents. But at the same time it has the disadvantages of computational complexity. The predictive distribution for the time to next incident can not always be expected to end up with a nice closed form even with conjugate priors. Thus we often encounter a numerical integration problem with high dimensions to obtain a predictive distribution, which is practically unsolvable for a model that involves many parameters. In order to circumvent this difficulty, we propose a method of approximation that essentially breaks down a problem involving many integrations into several repetitive steps so that each step involves only a small number of integrations.

• 한국항공우주학회지
• /
• 제40권5호
• /
• pp.405-412
• /
• 2012

높은 효율의 풍력터빈 블레이드을 위해 와류 셀이 장착된 에어포일의 정지상태 및 동실 속 상태에서의 유동제어 특성을 수치적으로 연구하였다. 수치기법은 Roe의 flux-difference-splitting을 사용한 격자점 중심 유한체적법과 이중시간 전진 기법을 사용하는 내재적 시간적분법을 사용하였다. 계산결과 와류 셀을 장착한 경우 셀 내부의 부압으로 인해 양항비증가를 얻을 수 있음을 확인하였다. 동실속의 경우 셀 내부의 와류에 의해 hysterisis 현상을 상당히 감소시킬 수 있음을 확인하였다.

• Structural Engineering and Mechanics
• /
• 제70권5호
• /
• pp.551-562
• /
• 2019

In the present study, a suitable mathematical model considering parabolic transverse shear strains for dynamic analysis of laminated composite skew plates under different types of impulse and spatial loads was presented for the first time. The proposed mathematical model satisfies zero transverse shear strain at the top and bottom of the plate. On the basis of the cubic variation of thickness coordinate in in-plane displacement fields of the present mathematical model, a 2D finite element (FE) model was developed including skew transformations in the mathematical model. No shear correction factor is required in the present formulation and damping effect was also incorporated. This is the first FE implementation considering a cubic variation of thickness coordinate in in-plane displacement fields including skew transformations to solve the forced vibration problem of composite skew plates. The effect of transverse shear and rotary inertia was incorporated in the present model. The Newmark-${\beta}$ scheme was adapted to perform time integration from step to step. The $C^0$ FE formulation was implemented to overcome the problem of $C^1$ continuity associated with the cubic variation of thickness coordinate in in-plane displacement fields. The numerical studies showed that the present 2D FE model predicts the result close to the analytical results. Many new results varying different parameter such as skew angles, boundary conditions, etc. were presented.

• 한국전산구조공학회논문집
• /
• 제20권5호
• /
• pp.521-532
• /
• 2007

원형축이 축방향으로 충격하중을 받으면 외경에서 반사된 파가 축의 중앙으로 집중되어 순간적으로 큰 응력이 발생하게 된다. 본 연구에서는 여러 가지 충격 축하중을 받는 횡등방성 반-무한 원형축을 대상으로 중실축 내의 종방향 응력전파를 축대칭 유한요소법과 Houtolt 시간적분법을 이용하여 프로그램을 작성하고 수치적으로 해석하여 그 결과를 횡등방성 재료의 재료구성비에 따라 자세히 설명한다. 제시된 해법의 타당성은 본 논문 수치 결과와 기 해석된 다른 해법에 의한 수치결과의 비교를 통해 검증된다. 여러 종류의 충격하중들에 따른 파동의 결과를 2차원, 3차원적으로 제시하여 축응력 전파를 이해하는데 기본 자료가 되도록 하였다. 또한 유한요소법을 이용하여 수치해석을 함에 있어 정확한 수치결과를 얻기 위한 무차원 동특성 시간변수에 대해 기술하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2007년도 춘계학술대회논문집
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.