• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.399-412
• /
• 2017

The Gauss-Legendre integral method is applied to numerically evaluate the Green function and its derivatives in finite water depth. In this method, the singular point of the function in the traditional integral equation can be avoided. Moreover, based on the improved Gauss-Laguerre integral method proposed in the previous research, a new methodology is developed through the Gauss-Legendre integral. Using this new methodology, the Green function with the field and source points near the water surface can be obtained, which is less mentioned in the previous research. The accuracy and efficiency of this new method is investigated. The numerical results using a Gauss-Legendre integral method show good agreements with other numerical results of direct calculations and series form in the far field. Furthermore, the cases with the field and source points near the water surface are also considered. Considering the computational efficiency, the method using the Gauss-Legendre integral proposed in this paper could obtain the accurate numerical results of the Green function and its derivatives in finite water depth and can be adopted in the near field.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.264-277
• /
• 2015

This article investigates various transcription techniques for the Legendre pseudospectral (PS) method to compare the pros and cons of each approach. Eight combinations from four different types of collocation points and two discretization methods for dynamic constraints, which differentiate Legendre PS transcription techniques, are implemented to solve a carefully selected test set of nonlinear optimal control problems (NOCPs). The convergence property and prediction accuracy are compared to provide a useful guideline for selecting the best combination. The tested NOCPs consist of the minimum time, minimum energy, and problems with state and control constraints. Therefore, the results drawn from this comparative study apply to the solution of similar types of NOCPs and can mitigate much debate about the best combinations. Additionally, important findings from this study can be used to improve the numerical efficiency of the Legendre PS methods. Three PS applications to the aerospace engineering problems are demonstrated to prove this point.

• 호남수학학술지
• /
• 제38권1호
• /
• pp.47-57
• /
• 2016

In the discrete formulation of the bubble stabilized Legendre Galerkin methods, the system of equations includes the artificial viscosity term as the parameter. We investigate the estimation of this parameter to get the least-squares solution which minimizes the sum of the squares of errors at each node points. Some numerical results are reported.

• 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.

• East Asian mathematical journal
• /
• 제16권1호
• /
• pp.97-103
• /
• 2000

In this work, we consider a finite difference operator $L^2_N$ corresponding to $$Lu:=-(u_{xx}+u_{yy})\;in\;{\Omega},\;u=0\;on\;{\partial}{\Omega}$$, in $S_{h^2,1}$. We derive the relation between the absolute value of the bilinear form $l_{N^2}$(u, v) on $S_{h^2,1}{\times}S_{h^2,1}$ and Sobolev $H^1$ norms.

• 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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 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

• 대한수학회보
• /
• 제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.

• International Journal of Aeronautical and Space Sciences
• /
• 제14권3호
• /
• pp.282-295
• /
• 2013

In order to communicate with a ground station, the tracking profile design problem for a directional antenna system is considered. Because the motions of the gimbal angles in the antenna system affect the image quality, the main object is to minimize the motion of the gimbal angles during the satellite's imaging phase. For this goal, parameter optimization problems in the imaging and maneuver phases are formulated separately in the body-frame, and solved sequentially. Also, several mechanical constraints, such as the limitation of the gimbal angle and rate, are considered in the problems. The tracking profiles of the gimbal angles in the maneuver phases are designed with N-th order polynomials, to continuously connect the tracking profiles between two imaging phases. The results confirm that if the vector trace of the desired antenna-pointing vector is within the antenna's beam-width angle, motions of the gimbal angles are not required in the corresponding imaging phase. Also, through numerical examples, it is shown that motion of the gimbal angles in the imaging phase can be minimized by the proposed design process.

• 한국전산구조공학회논문집
• /
• 제15권3호
• /
• pp.491-499
• /
• 2002

직교이방성 적층평판해석을 위해 퇴화 쉘요소에 기초를 둔 p-version 유한요소법이 제안되었다. 이 모델의 비선형 정식화과정에서 기하비선형의 경우 von Karman의 대변형-소변형률 가정을 설명하기 위해 Total Lagrangian 방법이 채택되었으며, 재료비선형의 경우 Huber-Mises의 항복기준과 변형률경화 항복함수에 근거를 둔 Prandtl-Reuss 유동법칙이 사용되었다. 재료모델은 이방성을 표현하는 매개변수에 의해 이방겅재료를 고려할 수 있도록 하였다. 적층평판이론으로는 전단변형 효과를 고려할 수 있는 등가단출이론(ESL Theory)에 기초를 두었기 때문에 두 적층간 계면에서의 전단변형률은 연속이라는 조건을 갖게된다 적분형 르장드르 다항식이 형상함수로 사용되었으며 형상함수의 차수는 1차에서 10차까지 변화시킬 수 있다. 또한, Causs-Lobatto 수치적될법을 사용하기 때문에 기존의 가우스 적분점에서 계산되던 응력값은 이 적분법의 적분점이 절점에 위치하므로 절점에서 바로 응력값이 산출되도록 하였다 극한하중 수렴성, 비선형 효과, 소성역의 형상 등의 비교관점을 통해 p-version 유한요소 모델의 적정성을 보이고자 하였다.