• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1998.04a
• /
• pp.17-17
• /
• 1998

연소반응을 가지는 후류(wake)는 가스터빈 연소실의 flame holder 등에서 발생한다. 후류유동의 안정성 혹은 불안정성은 이러한 유동에 있어서 많은 영향을 끼치므로 상당히 중요하다. 본 연구는 위와 같이 연소 반응에 의한 밀도구배를 가지는 후류유동에 대하여 불안정성 해석을 수행하였다. 교란에 대한 지배방정식은 Navier-Stokes 방정식에서 점성항을 제외한 Euler 방정식을 고려하였다. 충류유동의 압력은 일정하다고 가정하였다. 교란 방정식을 유도하기 위하여 충류 유동이 평행하여 유동 방향에 수직한 방향의 구배만이 존재한다고 가정하였다. 모든 변수들은 충류 유동의 값과 움직이는 파장의 형태를 가지는 작은 교란의 합으로 생각하여 압력에 대한 교란방정식을 구하며, shooting법과 Newton-Raphson법에 근거를 두는 반복법을 사용하여 풀었다. 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 비압축성의 경우 우선 Gaussian Profile 가정함과 동시에 연소반응을 포함하는 유동장을 보다 정확히 구하기 위하여 Navier-Stokes 방정식으로부터 구한 결과를 사용하였다. 연소반응을 포함하는 유동장을 구할 때에는 해석상 편의를 위해 예혼합물질은 이상기체로, 화학반응은 1단계의 비가역반응으로 가정하였으며, 반응열로 인한 부력의 효과는 무시하였다. 위와 같은 유동장을 가지고 불안정성 해석을 수행한 결과 후류유동은 두 개의 변곡점을 가지며 sinuous 모드와 varicose모드의 두 개의 불안정성 모드가 있음을 보였다. 밀도 변화가 있는 경우나 밀도 변화가 없는 경우 모두 sinuous 모드의 가장 불안정한 모드가 varicose 모드의 가장 불안정한 모드보다 더 불안정함을 보여주어 후류 유동은 자유 유동에 가까운 위상 속도를 가지는 sinuous 모드에 의해 지배될 것임을 예측할 수 있다. 연소반응이 완전연소에 가까울수록 그리고 압축성 효과가 클수록 유동내부의 온도가 증가하고 점성 또한 증가하여 후류유동은 안정됨을 알 수 있었다 유동변수들의 contour로부터 유동의 특성을 예측한 결과 baroclinic 항이 dilatational 항보다 상대적으로 크며, 중심선 상하에 생기는 vortex를 더욱 성장시킬 것으로 생각된다.

• Journal of the Korean Data and Information Science Society
• /
• v.25 no.5
• /
• pp.1025-1038
• /
• 2014

Linear mixed model has often been utilized for genetic association analysis with family-based samples. The correlation matrix for family-based samples is constructed with kinship coefficient and assumes that parental phenotypes are independent and the amount of correlations between parent and offspring is same as that of correlations between siblings. However, for instance, there are positive correlations between parental heights, which indicates that the assumption for correlation matrix is often violated. The statistical validity and power are affected by the appropriateness of assumed variance covariance matrix, and in this thesis, we provide the linear mixed model with flexible variance covariance matrix. Our results show that the proposed method is usually more efficient than existing approaches, and its application to genome-wide association study of body mass index illustrates the practical value in real data analysis.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.464-464
• /
• 2012

수문분야에 있어서 빈도해석의 목적은 특정 재현기간에 대한 발생 가능한 수문량의 규모를 파악하는데 있으며, 빈도해석의 정확도는 적합한 확률분포모형의 선택과 매개변수 추정방법에 의존하게 된다. 일반적으로 각 확률분포모형의 특성을 대표하는 매개변수를 추정하기 위해서는 모멘트 방법, 확률가중 모멘트 방법, 최대우도법 등을 이용하게 된다. 모멘트 방법에 의한 매개변수 추정은 해를 구하기 위한 과정이 단순한 반면, 비대칭형의 왜곡된 분포를 갖는 자료들에 대해서는 부정확한 결과를 나타내게 된다. 확률가중 모멘트 방법은 표본의 크기가 작거나 왜곡된 자료일 경우에도 비교적 안정적인 결과를 제공하는 반면, 확률 가중치가 정수로만 제한되는 단점을 갖고 있다. 그리고 대수 우도함수를 이용하여 매개변수를 추정하게 되는 최우도법은 가장 효율적인 매개변수 추정치를 얻을 수 있는 것으로 알려져 있으나, 비선형 연립방정식으로 표현되는 해를 구하기 위해서는 Newton-Raphson 방법을 사용하는 등 절차가 복잡하며, 때로는 수렴이 되지 않아 해룰 구하지 못하는 경우가 발생되게 된다. 이에 반해, 최근의 Genetic Algorithm, Ant Colony Optimization 및 Simulated Annealing과 같은 Meta-Heuristic Algorithm들은 복잡합 공학적 최적화 문제 있어서 효율적인 대안으로 주목받고 있으며, Hassanzadeh et al.(2011)에 의해 수문학적 빈도해석을 위한 매개변수 추정에 있어서도 그 적용성이 검증된바 있다. 본 연구의 목적은 연 최대강수 자료의 빈도해석에 적용되는 확률분포모형들의 매개변수 추정을 위해 Meta-Heuristic Algorithm을 적용하고자 함에 있다. 따라서 본 연구에서는 매개변수 추정을 위한 방법으로 Genetic Algorithm 및 Harmony Search를 적용하였고, 그 결과를 최우도법에 의한 결과와 비교하였다. GEV 분포를 이용하여 Simulation Test를 수행한 결과 Genetic Algorithm을 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 비교적 유사한 분포를 나타내었으나 과도한 계산시간이 요구되는 것으로 나타났다. 하지만 Harmony Search를 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 유사한 분포를 나타내었을 뿐만 아니라 계산시간 또한 매우 짧은 것으로 나타났다. 또한 국내 74개소의 강우관측소 자료와 Gamma, Log-normal, GEV 및 Gumbel 분포를 이용한 실증연구에 있어서도 Harmony Search를 이용한 매개변수 추정은 효율적인 매개 변수 추정치를 제공하는 것으로 나타났다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.387-387
• /
• 2015

최근 기후변화에 따른 국지성 집중호우 및 돌발홍수 증가로 도심지역에 많은 침수 피해가 발생하고 있다. 본 연구는 이러한 내외수에 영향을 받는 도시 중소하천의 침수 저감을 위하여 침수 예측모형에 적용 가능한 수리학적 계산모형을 개발하는 것이다. 일반적으로 자연하천은 수지형 수계가 대부분이며, 도시 우수관로와 관계배수 시스템 등은 폐합형 수계에 포함되어 수지형 계산모형으로는 도시 중소하천의 유출 모의를 할 수 없다. 폐합형 수계의 계산모형은 수지형에 비해 복잡하지만, 적용대상 수계가 폐합형이 아닌 경우에도 합류점의 유입량을 처리하기가 편리하고, 역방향의 월류 흐름이 존재하는 감조하천에서의 월류 흐름 모의가 가능한 장점을 갖고 있어 도시 내배수 시스템은 물론 자연하천에도 적용 할 수 있다. 본 모형은 절점, 수로 및 계산점으로 구성되는데 동력학적 방법인 1차원 Saint-Venant의 연속 방정식과 운동량 방정식에 수치해법을 이용하여 구하고자 하는 시간과 지점의 수위와 유량을 계산할 수 있게 구성하였다. 수치해법으로는 가장 보편적으로 사용되는 유한차분법 중 안정성과 정확성이 우수한 것으로 평가된 Preissmann의 4점 음해법으로 차분방정식에 Newton-Raphson 방법을 사용하여 유량과 수위 보정치에 Double Sweep 알고리즘을 적용하였다. 유한차분법은 안정성 문제를 수반할 수 있기에 시간 증분이 작을수록 Courant 조건을 만족할 수 있다. 모형 비교 검증을 위하여 동력학적 방법을 적용한 대표적인 폐합형 수계모형인 EPA SWMM을 지원하는 CHI사의 PC-SWMM 프로그램을 이용하여 가상의 폐합 수계를 구축하였다. 일반적으로 상류지점의 경계조건은 하류로 추적될 입력 자료로써, 상류지점과 합류지점은 유입되는 유량값을 그리고 하류지점은 유출되는 수위값을 경계조건으로 입력하였다. 운동량방정식의 에너지 경사와 마찰경사 항에 포함된 조도계수는 변화량에 중요한 물리적 요소이지만, 고정 상수값인 0.03을 적용하여 검증에 용의하도록 하였다. 구축된 모형과 PC-SWWM을 통해 산출된 계산점별 수위와 유량에 RMS 오차를 비교한 결과 만족할만한 결과를 얻을 수 있었다. 따라서 향후 내외수를 연계에 침수예측모형에 적용 가능할 것으로 사료된다.

• Plant Journal
• /
• v.13 no.4
• /
• pp.48-50
• /
• 2017

A program for analyzing the transient behaviors of industrial gas turbines was developed. Each component (compressor, combustor, turbine and ducts)of gas turbine is modeled as a fully module to enhance the expandability of the program. We used object-oriented programing for this purpose. The mass and energy balance equations are solved numerically by Multivariable Newton Raphson method. The characteristic maps for the compressor and turbine were used for predicting the performance of a gas turbine engine. Combustion in the combustor is assumed to be complete combustion. PID control is used to maintain constant rotational speed and turbine exhaust temperature by the control of the fuel flow rate and the changing of the compressor inlet guide vane angle at the same time. It was confirmed that stable control of the gas turbine was possible, even for a rapid load change.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.785-802
• /
• 2018

In this paper, two different computational methods, called Rayleigh-Ritz and collocation are developed to estimate the ultimate strength of composite plates. Progressive damage behavior of moderately thick composite laminated plates is studied under in-plane compressive load and uniform lateral pressure. The formulations of both methods are based on the concept of the principle of minimum potential energy. First order shear deformation theory and the assumption of large deflections are used to develop the equilibrium equations of laminated plates. Therefore, Newton-Raphson technique will be used to solve the obtained system of nonlinear algebraic equations. In Rayleigh-Ritz method, two degradation models called complete and region degradation models are used to estimate the degradation zone around the failure location. In the second method, a new energy based collocation technique is introduced in which the domain of the plate is discretized into the Legendre-Gauss-Lobatto points. In this new method, in addition to the two previous models, the new model named node degradation model will also be used in which the material properties of the area just around the failed node are reduced. To predict the failure location, Hashin failure criteria have been used and the corresponding material properties of the failed zone are reduced instantaneously. Approximation of the displacement fields is performed by suitable harmonic functions in the Rayleigh-Ritz method and by Legendre basis functions (LBFs) in the second method. Finally, the results will be calculated and discussions will be conducted on the methods.

• The Journal of Korea Institute of Information, Electronics, and Communication Technology
• /
• v.15 no.6
• /
• pp.492-505
• /
• 2022

This study is a study on the results of improving the logic to effectively transmit and decode compressed data in real time by improving the encoding and decoding performance of BL-beta codes that can be used for lossless real-time transmission of IoT sensing data. The encoding process of BL-beta code includes log function, exponential function, division and square root operation, etc., which have relatively high computational burden. To improve them, using bit operation, binary number pattern analysis, and initial value setting of Newton-Raphson method using bit pattern, a new regularity that can quickly encode and decode data into BL-beta code was discovered, and by applying this, the encoding speed of the algorithm was improved by an average of 24.8% and the decoding speed by an average of 5.3% compared to previous study.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6D
• /
• pp.895-904
• /
• 2008

CWR (Continuous Welded Rail) may be broken when a temperature drop below the neutral temperature changes in axial force, causing tensile fracture and rail gap, in winter. Rail-breaks may lead to the damage of the rail and wheel by dynamic load, and the reduction of running safety if not detected before the passage of a train. In this study, the track and train coupled model with open gap for dynamic interaction analysis, is proposed. Linear track and train systems is coupled by the nonlinear Herzian contact spring and the complete system matrices of total track-train system is constructed. And the interaction phenomenon considering open gap, was defined by assigning the irregularity functions between the two sides of a gap. Time history analysis, which have an iteration scheme such as $Newmark-{\beta}$ method based on Modified Newton-Raphson methods, was performed to solve the nonlinear equation. Finally, numerical studies are performed to assess the effect of various parameters of system, apply to various speeds, open gap size and the support stiffness of rail.