• 대한조선학회지
• /
• 제26권1호
• /
• pp.39-45
• /
• 1989

선체구조물이나 해양구조물의 동적응답중 원치 않는 고유진동수와 고유진동형태를 가지게 되는 경우가 있으며, 이러한 구조물은 동적 구조 재설계가 필수적이다. 본 소고에서는 비감쇄 구조물의 고유진동수와 진동형태를 기진력에 의한 특정한 진동수와 공진하지 않도록 또는 구조물의 중요한 부분이 특정 진동형태의 최대치에 오지 않도록 구조물의 질량과 강성을 최적하게 변화시키는 방법에 대해 논의하고 있다, 이 방법은 기존의 방법에서 사용되는 모든 고유진동형태의 수식포함과 달리 구속된 고유진동형태만을 미지수로 수식중에 사용하여 불필요한 계산과정을 줄이고 있다. 동적 구조 재설계중 최적화 문제에 중점을 두었으며 목적함수로는 구조물의 최소의 변화와 또는 최소의 중량을 취하였고, 예제를 통하여 본 방법의 응용과 효율성이 입증되었다. 예제에서는 간단한 구조물을 다루었으나 본 방법은 상용 유한요소코드의 연계이용으로 각종 선체구조물과 해양구조물의 진동문제해결에 응용될 수 있음은 자명한 일이다.

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

• 한국습지학회지
• /
• 제7권3호
• /
• pp.57-66
• /
• 2005

본 연구는 최적화모형을 이용하여 관측된 단일 강우-유출 사상으로부터 최적단위도와 침투율 공식의 매개변수를 결정하였다. 단일 강우-유출사상으로부터 절대오차누계를 최소로 하는 선형계획모형을 정립하고, 황구지천 수직교 지점에서의 최적침투율과 단위도를 결정하였다. 침투율 공식의 섭동 단계에서는 시산법을 적용하여 Kostiakov, Philip, Horton 및 Green-Ampt 공식의 매개변수를 결정하였다. ${\Phi}$지표법에 의한 단위도 종거의 값은 유일하게 결정되지만 침투율 공식의 매개변수에 따라 단위도의 종거는 변한다. 본 연구에서는 관측된 단일 강우-유출사상에 관하여 제안된 모형을 적용하였으며 종래의 방법에 의한 오차보다 작은 단위도를 구할 수 있었다.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제9권1호
• /
• pp.79-86
• /
• 2004

암석 내 균열을 따른 수리전도도는 균열의 기하학적 요소, 즉 방향, 간극, 거칠기 그리고 상호 연결도에 주로 좌우된다. 따라서, 균열 내 투수계수를 정확하게 계산하기 위해서는 이와 같은 기하 요소들을 최대한 계산모델에 반영할 필요가 있다. 이 연구에서는 균열 기하양상을 최대한 정확히 반영한 균열모델에서 기존 수치해석과는 다른 새로운 방법인 균질화 해석법(homogenization analysis method)을 이용하여 균열을 따른 투수계수를 구하기 위해 수치해석을 수행하였다. 먼저, 공초점 레이저 스캔 현미경(Confocal Laser Scanning Microscope)을 이용하여 암석시료의 균열 조도와 균열에 가한 수직압축력의 변화에 따른 간극 변화량을 직접 측정하고, 이와 같이 획득한 자료는 균열모델 재현을 위한 입력자료로 사용되었다. 재현된 균열모델을 토대로 한 균질화 해석법은 미시규모(microscale) 매질특성과 거시규모(macroscale) 매질특성을 동시에 고려하여 투수계수를 계산할 수 있는 것이다. 즉, 균질화 해석법은 주기적 미세구조(microstructure)를 갖는 미소 불균질 물질의 거동특성을 구명하기 위해 개발된 새로운 형태의 섭동(perturbation) 이론이다. 이는 균질한 미시규모에서 미시 투수특성을 계산한 후, 거시규모에서의 균질화 투수계수를 계산하게 된다. 그러므로, 이 방법은 균열 기하양상의 국부적 영향을 고려한 투수특성을 정확히 해석할 수 있다. 균질화법을 이용한 투수계수 산정결과를 기존 연구에서 제안한 경험식과 비교하여 그 타당성을 검증하기 위해 전술한 2차원 균열모델을 이용한 투수계수 계산을 수행하였다. 균열모델은 거칠기(roughness)를 반영하고 동일한 간극을 할당한 평행판 모델을 가정하였다. 계산결과에 의하면, 균질화 해석법에 의해 계산한 C-투수계수는 실내투수시험에 의해 구한 투수계수와 같은 범위의 값을 가지거나 $10^1$ 정도의 차이를 보여, 그 계산결과는 타당하다고 볼 수 있다. 그러나, 균질화 해석법은 국부적으로 불균질한 균열 기하양상과 물질특성이 미시규모와 거시규모에서 모두 고려되므로, 이들 특성을 정확히 알고 있을 경우 기존에 제안된 경험식들에 의한 계산결과 보다 균질화 해석법의 결과가 훨씬 정확함을 주목하여야 한다.

• Journal of Astronomy and Space Sciences
• /
• 제22권1호
• /
• pp.21-34
• /
• 2005

본 연구에서는 타원궤도상에서 위성의 편대비행을 유지하기 위하여 필요한 포기조건을 결정하고자 한다. 타원궤도일 경우 Hill 방정식으로는 위성간의 상대운동을 기술할 수 없기 때문에, Hill 방정식의 초기조건에 비선형성과 이심률에 대한 보정을 하여 얻은 새로운 운동방정식을 사용했다. 편대비행에서 상대적 거리를 유지하기 위하여 주위성과 부위성의 평균각속도를 일치시키는 구속조건을 이용했다. 이 구속조건은 J2 섭동항을 고려한 것이므로, 이 구속조건을 만족하는 편대비행의 초기조건은 타원궤도에서의 위성편대비행을 유지하는데 잘 적용될 수 있다. 타원궤도에서의 상대운동방정식 초기조건에 J2 섭동을 고려한 구속조건을 적용할 때, 이심률이 0.05 이하이고 위성간의 상대거리가 0.5km 정도인 경우만이 주기적으로 일정하게 간격이 유지되는 결과를 얻을 수 있다. 따라서 이심률이 크지 않은 타원궤도에서는 평균각속도 일치의 구속조건을 사용하여 위성간의 상대거리를 유지할 수 있었다. 이러한 결과를 이용하여 타원궤도에서의 위성편대비행을 위한 효율적인 초기조건을 제공할 수 있고, 위성편대비행의 운용에 있어서 비용을 절감할 수 있는 방법을 제시할 수 있다.

• KSTLE International Journal
• /
• 제7권2호
• /
• pp.27-34
• /
• 2006

Knudsen number is the ratio of molecular mean free path versus mm thickness and the criterion to determine the flow form. When its value is lower than 0.01, the flow can be assumed to has no slip boundary condition. And in the case that the value is between 0.01 and 10, then the flow has slip boundary condition at both the adjacent walls. The condition of the air flow between the rotating journal and top foil in the air foil bearing is determined by the rotating speed and load, and the Knudsen number is also varied by those values. Because the molecular mean free path is variable to the pressure and temperature, more exact formulation is necessary to understand and analyze the flow regime. In this study, the analysis considering Knudsen number formulated with those variables (pressure, temperature and mm thickness) was executed. The approximate value was examined using the equation to confirm whether the flow has the slip or no-slip boundary condition. From the analytic investigation, it was decided to range approximately 0.01 to 1.0 and the flow can be supposed to have the slip boundary condition. Under the condition of the slip flow, the static characteristics of the air foil bearing were examined using modified Reynolds equations. The results were compared with those considering no slip condition. It shows that the slip condition makes the flow decelerates and the load carrying capacity decreases compared with no slip condition. And as the bearing number and eccentricity ratio increase, the load carrying capacity also increased at both the cases. From this result, it can be supposed that the bearing torque also increases. In the analysis of the dynamic characteristics, the perturbed Knudsen number was taken into consideration. Because the Knudsen number is expressed as the terms of each variable, the perturbed equation can be simply derived. The results of both cases considering and not considering Knudsen number were compared each other. In the case of the direct terms of the stiffness and damping coefficients, the difference between both cases was little and increased as the bearing number and eccentricity ratio increased. And the cross terms have less or more differences.

• 한국광학회지
• /
• 제16권3호
• /
• pp.182-190
• /
• 2005

본 논문은 효율적으로 라만 증폭기의 이득을 예측하고 잡음 특성을 분석하기 위해서 라만 증폭기 상관된 상미분방정식 (Ordinary Differential Equation: ODE)을 유효거리(Effective Length) 기반의 상관된 적분형(closed integral form)방정식 및 매트릭스로 전개하고 이 전개를 활용한 라만 증폭기 모델링 및 수치해석 알고리즘을 기술한다. 광섬유 라만 증폭기는 유연하고 넓은 이득 대역폭 및 낮은 잡음 등의 장점 때문에 최근 광통신 시스템에서 핵심기술로 각광받고 있으며, 특히 멀티채널 펌핑구조에서 성능예측을 위해 많은 라만 증폭기 모델링 방법들이 연구되어 왔다. 그러나 기존의 많은 연구들은 라만 증폭기 상관된 상미분방정식의 해를 "fiber propagation axis"를 기반으로 구해왔기 때문에 광섬유 길이에 의존적이고 복잡한 계산으로 상당히 많은 시간이 필요했으며, 실제 전송 시스템에서의 활용이 어려웠다. 이러한 문제점을 해결하기 위해서 본 논문에서는 기존의 상관된 상미분방정식을 접근이 용이한 "유효거리" 기반의 적분형 방정식으로 전개하고 매트릭스 및 벡터 형태로 알고리즘을 공식화하여 빠른 라만 이득 계산과 "iteration axis"를 이용한 해의 도출을 통해 새로운 라만 증폭기 모델링 방법과 수치해석 알고리즘을 제시하였다. 제안한 수치해석 알고리즘을 전방, 후방 및 양방향 펌핑구조의 라만 증폭기가 도입된 시스템에 적용하는 컴퓨터 모의실험을 수행한 결과 기존의 "Average power method와 비교하여 라만 이득과 광선로 내의 펌프 및 신호 광의 진행정도를 18배 이상의 정밀도에서 0.03 dB 이내의 매우 작은 오차범위 및 100배 이상 단축된 짧은 시간으로 정확히 예측하였다. 또한, 수치해석 알고리즘을 통해서 얻은 신호 광 파워의 분포를 바탕으로 Amplified Spontaneous Emission(ASE), 후방 ASE의 레일레이 산란 및 신호의 이중 레일레이 산란과 같은 라만 증폭기 잡음 요소들을 분석하였다. 새롭게 제안한 수치해석 알고리즘은 실제 광통신 시스템에 적용되어 신속하고 효율적으로 라만 증폭기 성능을 예측하고 분석할 수 있을 것으로 기대된다.