• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2681-2692
• /
• 1996

The finite-volume method for radiation in a three-dimensional non-orthogonal gas turbine combustion chamber with absorbing, emitting and anisotropically scattering medium is presented. The governing radiative transfer equation and its discretization equation using the step scheme are examined, while geometric relations which transform the Cartesian coordinate to a general body-fitted coordinate are provided to close the finite-volume formulation. The scattering phase function is modeled by a Legendre polynomial series. After a benchmark solution for three-dimensional rectangular combustor is obtained to validate the present formulation, a problem in three-dimensional non-orthogonal gas turbine combustor is investigated by changing such parameters as scattering albedo, scattering phase function and optical thickness. Heat flux in case of isotropic scattering is the same as that of non-scattering with specified heat generation in the medium. Forward scattering is found to produce higher radiative heat flux at hot and cold wall than backward scattering and optical thickness is also shown to play an important role in the problem. Results show that finite-volume method for radiation works well in orthogonal and non-orthogonal systems.

• Structural Engineering and Mechanics
• /
• v.73 no.3
• /
• pp.259-269
• /
• 2020

In this paper, bending-stretching analysis of thick functionally graded piezoelectric rectangular plates is studied using the higher-order shear and normal deformable plate theory. On the basis of this theory, Legendre polynomials are used for approximating the components of displacement field. Also, the effects of both normal and shear deformations are encountered in the theory. The governing equations are derived using the principle of virtual work and variational approach. It is assumed that plate is made of piezoelectric materials with functionally graded distribution of material properties. Hence, exponential function is used to modify mechanical and electrical properties through the thickness of the plate. Finally, the effect of material properties, electrical boundary conditions and dimensions are investigated on the static response of plate. Also, it is shown that results of the presented model are close to the three dimensional elasticity solutions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.10a
• /
• pp.38-45
• /
• 1994

P-version finite element model for the computation of stress intensity factors in two dimensional cracked panels by J-integral method is presented. The proposed model is based on high order theory and hierarchical shape function. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three part such as basic mode, side mode, integral mode. The stress intensity factors are computed by J-integral method. The example models for validating the proposed p-version model are centrally cracked panel, single and double edged crack in a rectangular panel under pure Mode I. And the analysis results are compared with those by the h-version of FEM and empirical solutions in literatures. Very good agreement with the existing solution are shown.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.4 s.35
• /
• pp.82-97
• /
• 1999

A deep sea vehicle must be designed to ensure its safety under ultra-high pressure circumstances. If a pressure housing of a deepsea vehicle is collapsed by ultra-high pressure, the deepsea vehicle may be lost. The objective of this paper is to introduce a design collapse pressure for the deep sea pressure vessel which is composed of one cylinder and two hemispheres. Especially the collapse pressure of hemispherical shell with a hole at top is analyzed by a variational approach (weighted residual method). And for the purpose of design, the salty factor of collapse pressure is presented which is analyzed by interpolation method.

• Honam Mathematical Journal
• /
• v.35 no.4
• /
• pp.701-706
• /
• 2013

Summation theorems for hypergeometric series $_2F_1$ and generalized hypergeometric series $_pF_q$ play important roles in themselves and their diverse applications. Some summation theorems for $_2F_1$ and $_pF_q$ have been established in several or many ways. Here we give a proof of Watson's classical summation theorem for the series $_3F_2$(1) by following the same lines used by Rakha [7] except for the last step in which we applied an integral formula introduced by Choi et al. [3].

• Structural Engineering and Mechanics
• /
• v.2 no.3
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 1993.11a
• /
• pp.99-102
• /
• 1993

A consistent general order nodal method for solving the three-dimensional neutron diffusion equation in (x-y-z) geometry has been derived by using a weighted integral technique and expanding the spatial variable by the Legendre orthogonal series function. The equation set derived can be converted into any order nodal schemes. It forms a compact system for general order of nodal moments. The method utilizes fewer unknown variables in the schemes for iterative-convergence solution than other nodal methods listed in the literatures, and because the method utilizes the analytic solutions of the transverse-integrated one dimensional equations and a consistent approximation for a given spatial variable through all the solution procedures, which renders the use of an approximation for the transverse leakages no longer necessary, we can expect extremely accurate solutions and the solution would converge exactly when the mesh width is decreased or the approximation order is increased.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.04a
• /
• pp.62-69
• /
• 1996

In this study, an equivalent domain integral (EDI) method is presented to estimate the track-till integral parameter, J-value, for two dimensional cracked elastic bodies which may quantify the severity of the crack-tit) stress fields. The conventional J-integral method based on line integral has been converted to equivalent area or domain integrals by using the divergence theorem. It is noted that the EDI method is very attractive because all the quantities necessary for computation of the domain integrals are readily available in a finite element analysis. The details and its implementation are extened to both h-version finite element model with 8-node isoparametric element and p-version finite element model with high order hierarchic element using Legendre type shape fuctions. The variations with respect to the different path of domain integrals from the crack-tip front and the choice of 5-function have been tested by several examples.

• Communications for Statistical Applications and Methods
• /
• v.17 no.6
• /
• pp.811-827
• /
• 2010

This paper deals with a density estimation method in binary choice models that can be regarded as a statistical inverse problem. We use an orthogonal basis to estimate density function and consider the choice of an appropriate truncation parameter to reflect the model complexity and the prediction accuracy. We propose a data-dependent rule to choose the truncation parameter in the context of binary choice models. A numerical simulation is provided to illustrate the performance of the proposed method.