• 대한기계학회논문집
• /
• 제15권2호
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

• 한국정밀공학회지
• /
• 제21권7호
• /
• pp.76-82
• /
• 2004

The solutions of neighboring optimal control are typically obtained using the sweep method or transition matrices. Due to the numerical integration, however, the gain matrix can become infinite as time go to final one in the transition matrices, and the Riccati solution can become infinite when the final time free. To overcome these disadvantages, this paper proposes the pseudospectral Legendre method which is to first discreteize the linear boundary value problem using the global orthogonal polynomial, then transforms into an algebraic equations. Because this method is not necessary to take any integration of transition matrix or Riccati equation, it can be usefully used in real-time operation. Finally, its performance is verified by the numerical example for the space vehicle's orbit transfer.

• Steel and Composite Structures
• /
• 제29권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.

• International Journal of Aeronautical and Space Sciences
• /
• 제10권1호
• /
• pp.112-121
• /
• 2009

This paper deals with supersonic air-breathing missile system. A supersonic air-breathing missile system has very complicated and incoherent thrust characteristics with respect to outer and inner environment during operation. For this reason, the missile system has many maneuver constraints and is allowed to operate within narrow flight envelope. In this paper, trajectory optimization of the missile is accomplished. The trajectory optimization problem is formulated as a discrete parameter optimization problem. For this formulation, Legendre Pseudo-Spectral method is introduced. This method is based on calculating the state and control variables on Legendre-Gauss-Lobatto (LGL) points. This approach helps to find approximated derivative and integration quantities simply. It is shown that, for this trajectory optimization, trend analysis is performed from thrust characteristics on various conditions so that the trajectory optimization is accomplished with fine initial guess with these results.

• 지구물리와물리탐사
• /
• 제15권2호
• /
• pp.57-65
• /
• 2012

탄성파 수치 모형 계산에 있어서 다양한 방법들이 개발되어 적용되었다. 최근에는 특히 탄성파 수치 모형 계산에 있어 혁신적인 방법인 SEM (Spectral Element Method)가 개발되어 사용되어 왔다. 이 방법은 지형을 자유롭게 표현하는데 있어 유연한 유한요소법의 장점에 정확성을 높인 방법이다. 일반적으로 Weak Formulation 형태의 파동방정식에 육면체 요소와 Gauss-Lobatto-Legendre 적분법을 적용한 방법이 널리 사용된다. 일반적인 SEM에서는 PML (Perfectly Matched Layer)경계조건을 적용하기 어려워 속도-응력 변분식으로 파동방정식을 변경하였다. CFS-PML (Complex frequency Shifted PML)경계조건을 ADE (Auxiliary Differential Equation)방정식으로 변경하여 속도-응력 파동방정식에 적용함으로써 분리할 필요가 없는 PML을 적용한 SEM 수치 모형 계산 알고리듬을 구현하였다. 1차원 수치모형과 3차원 수치모형 실험을 통하여 SEM에 적용한 비분리 CFS-PML이 유한경계에서 인공적으로 반사되는 반사파를 효과적으로 제거하는 것을 확인하였다.

• Structural Engineering and Mechanics
• /
• 제66권5호
• /
• pp.557-568
• /
• 2018

An investigation is made in the present work on the post-buckling and geometrically nonlinear behaviors of moderately thick perfect and imperfect rectangular plates made-up of functionally graded materials. Spectral collocation approach based on Legendre basis functions is developed to analyze the functionally graded plates while they are subjected to end-shortening strain. The material properties in this study are varied through the thickness according to the simple power law distribution. The fundamental equations for moderately thick rectangular plates are derived using first order shear deformation plate theory and taking into account both geometric nonlinearity and initial geometric imperfections. In the current study, the domain of interest is discretized with Legendre-Gauss-Lobatto nodes. The equilibrium equations will be obtained by discretizing the Von-Karman's equilibrium equations and also boundary conditions with finite Legendre basis functions that are substituted into the displacement fields. Due to effect of geometric nonlinearity, the final set of equilibrium equations is nonlinear and therefore the quadratic extrapolation technique is used to solve them. Since the number of equations in this approach will always be more than the number of unknown coefficients, the least squares technique will be used. Finally, the effects of boundary conditions, initial geometric imperfection and material properties are investigated and discussed to demonstrate the validity and capability of proposed method.

• 정보처리학회논문지D
• /
• 제10D권7호
• /
• pp.1145-1148
• /
• 2003

공학문제 해결을 위한 수치적 프로그램에서 원하는 해와 그 해의 변이 값에 대하여 같은 수준의 오차를 유지할 수 있는 기존의 복합유한 요소방법을 소개하고 이에 대한 효과적인 프로그램 재사용을 이용한 Matrix 생성기법을 소개한다. 또한, 원하는 임의의 차수의 기저에 대한 Matrix의 자동 생성기법을 제안한다. 여기서, 자동 생성된 Matrix는 최소한의 nonzero element를 갖고, 이는 Inverse Matriix 형성에 있어서 최소오차와 효율성을 보장한다. 위에서 제안한 MatriBt 생성기법을 최소표면적 문제에 적용하여 본다.

• 대한수학회보
• /
• 제51권2호
• /
• pp.595-612
• /
• 2014

The spectral collocation method for a second order elliptic boundary value problem on a domain ${\Omega}$ with curved boundaries is studied using the Gordon and Hall transformation which enables us to have a transformed elliptic problem and a square domain S = [0, h] ${\times}$ [0, h], h > 0. The preconditioned system of the spectral collocation approximation based on Legendre-Gauss-Lobatto points by the matrix based on piecewise bilinear finite element discretizations is shown to have the high order accuracy of convergence and the efficiency of the finite element preconditioner.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 봄 학술발표회 논문집
• /
• pp.319-326
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted for in the sense of von Karman hypothesis. The material model Is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The Integrals of Legendre Polynomials we used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone

• 한국전산구조공학회논문집
• /
• 제31권2호
• /
• pp.87-95
• /
• 2018

이 논문에서는 스펙트럴 요소법과 외연적 시간적분법을 이용해 SH파의 전파 거동을 계산하는 수치해석 기법을 제시한다. 2차원 영역에서의 탄성파 해석을 위해 해석영역을 유한 영역으로 한정하고 파동이 반사되지 않도록 수치적 파동흡수 경계조건인 perfectly matched layer(PML)를 도입하였다. PML이 포함된 시간영역 파동방정식의 유한요소해법을 위해 스펙트럴 요소법을 적용하였고 Legendre- Gauss-Lobatto 수치적분법을 사용하여 질량행렬을 대각화하였다. 2차 미분방정식 시스템의 파동방정식을 1차 미분방정식 시스템으로 변환하였고 병렬화를 통한 탄성파 해석 성능의 최적화를 위해 외연적 시간적분법인 4차 Runge-Kutta 방법을 이용해 해석영역에서의 변위응답을 계산하였다. 2차원 해석영역에서 SH파의 전파 거동을 계산하는 수치예제를 통해 제시한 외연적 스펙트럴 요소법의 정확성을 검증하였고 PML로 인한 반사파의 감쇠효과를 확인하였다. 외연적 시간적분법을 통한 탄성파 해석 기법은 3차원 영역과 같은 대규모 문제에서의 탄성파 수치해석을 효율적으로 수행할 수 있을 것으로 기대된다.