• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.10a
• /
• pp.435-442
• /
• 2001

Derivation procedures of exact dynamic element stiffness matrix of shear deformable nonsymmetric thin-walled straight beams are rigorously presented for the spatial free vibration analysis. An exact dynamic element stiffness matrix is established from governing equations for a uniform beam element with nonsymmetric thin-walled cross section. First this numerical technique is accomplished via a generalized linear eigenvalue problem by introducing 14 displacement parameters and a system of linear algebraic equations with complex matrices. Thus, the displacement functions of dispalcement parameters are exactly derived and finally exact stiffness matrices are determined using member force-displacement relationships. The natural frequencies are evaluated and compared with analytic solutions or results of the analysis using ABAQUS' shell elements for the thin-walled straight beam structure in order to demonstrate the validity of this study.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.385-392
• /
• 2002

Derivation procedures of exact elastic element stiffness matrix of thin-walled curved beams are rigorously presented for the static analysis. An exact elastic element stiffness matrix is established from governing equations for a uniform curved beam element with nonsymmetric thin-walled cross section. First this numerical technique is accomplished via a generalized linear eigenvalue problem by introducing 14 displacement parameters and a system of linear algebraic equations with complex matrices. Thus, the displacement functions of displacement parameters are exactly derived and finally exact stiffness matrices are determined using member force-displacement relationships. The displacement and normal stress of the section are evaluated and compared with thin-walled straight and curved beam element or results of the analysis using shell elements for the thin-walled curved beam structure in order to demonstrate the validity of this study.

• Structural Engineering and Mechanics
• /
• v.59 no.1
• /
• pp.101-122
• /
• 2016

Elasticity solutions for bi-directional functionally graded beams subjected to arbitrary lateral loads are conducted, with emphasis on the end effects. The material is considered macroscopically isotropic, with Young's modulus varying exponentially in both axial and thickness directions, while Poisson's ratio remaining constant. In order to obtain an exact analysis of stress and displacement fields, the symplectic analysis based on Hamiltonian state space approach is employed. The capability of the symplectic framework for exact analysis of bi-directional functionally graded beams has been validated by comparing numerical results with corresponding ones in open literature. Numerical results are provided to demonstrate the influences of the material gradations on localized stress distributions. Thus, the material properties of the bi-directional functionally graded beam can be tailored for the potential practical purpose by choosing suitable graded indices.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1997.10a
• /
• pp.609-612
• /
• 1997

This paper deals with in-plane vibration analysis of curved beams. The exact dynamic element method is applied to obtain the dynamic model for curved beams. Numerical examples are provided to validate the proposed modeling and analysis method. The numerical results show that the proposed method is useful for the dynamic analysis of curved beams.

• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.279-296
• /
• 2008

This paper presents an exact integration for the hypersingular boundary integral equation of two-dimensional elastostatics. The boundary is discretized by straight segments and the physical variables are approximated by discontinuous quadratic elements. The integral for the hypersingular boundary integral equation analysis is given in a closed form. It is proven that using the exact integration for discontinuous boundary element, the singular integral in the Cauchy Principal Value and the hypersingular integral in the Hadamard Finite Part can be obtained straightforward without special treatment. Two numerical examples are implemented to verify the correctness of the derived exact integration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.4
• /
• pp.417-424
• /
• 1983

This paper presents two methods of obtaining approximate analytic solutions for the temperature distributions and heat flow to two-dimensional transient heat conduction problems in a finite strip with constant thermal properties using the Heat Balance Integral. The methods introduced in this study are as follows; one using the Heat Balance Integral only, and the other successively using the Heat Balance Integral and an exact analytic method. Both methods are applicable to a large number of the two-dimensional unsteady conduction problems in finite regions such as extended surfaces with uniform thickness, but in this paper only solutions for the unsteady problems in a finite strip with boundary condition at the base expressed in terms of step function are provided as an illustration. Results obtained by both methods are compared with those by the exact two-dimensional transient analysis. It is found that both approximate methods generate small time solutions, which can not be obtained easily by any exact analytic method for small values of Fourier numbers. In the case of applying the successive use of the Heat Balance Integral and Laplace transforms, the analysis shows good agreement with the exact solutions for any Fourier number in the range of Biot numbers less than 0.5.

• Geomechanics and Engineering
• /
• v.1 no.2
• /
• pp.121-142
• /
• 2009

The paper deals with the applications of spectral finite element method to the dynamic analysis of framed foundations supporting high speed machines. Comparative performance of approximate dynamic stiffness methods formulated using static stiffness and lumped or consistent or average mass matrices with the exact spectral finite element for a three dimensional Euler-Bernoulli beam element is presented. The convergence of response computed using mode superposition method with the appropriate dynamic stiffness method as the number of modes increase is illustrated. Frequency proportional discretisation level required for mode superposition and approximate dynamic stiffness methods is outlined. It is reiterated that the results of exact dynamic stiffness method are invariant with reference to the discretisation level. The Eigen-frequencies of the system are evaluated using William-Wittrick algorithm and Sturm number generation in the $LDL^T$ decomposition of the real part of the dynamic stiffness matrix, as they cannot be explicitly evaluated. Major's method for dynamic analysis of machine supporting structures is modified and the plane frames are replaced with springs of exact dynamic stiffness and dynamically flexible longitudinal frames. Results of the analysis are compared with exact values. The possible simplifications that could be introduced for a typical machine induced excitation on a framed structure are illustrated and the developed program is modified to account for dynamic constraint equations with a master slave degree of freedom (DOF) option.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.4
• /
• pp.1303-1315
• /
• 2018

For principal component analysis (PCA) to efficiently analyze large scale matrices, it is crucial to find a few singular vectors in cheaper computational cost and under lower memory requirement. To compute those in a fast and robust way, we propose a new stochastic method. Especially, we adopt the stochastic variance reduced gradient (SVRG) method [11] to avoid asymptotically slow convergence in stochastic gradient descent methods. For that purpose, we reformulate the PCA problem as a unconstrained optimization problem using a quadratic penalty. In general, increasing the penalty parameter to infinity is needed for the equivalence of the two problems. However, in this case, exact penalization is guaranteed by applying the analysis in [24]. We establish the convergence rate of the proposed method to a stationary point and numerical experiments illustrate the validity and efficiency of the proposed method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.1 s.256
• /
• pp.105-112
• /
• 2007

The linkage framework of geometric modeling based on NURBS(Non-Uniform Rational B-Spline) surface and shell finite analysis is developed in the present study. For this purpose, geometrically exact shell finite element is implemented. NURBS technology is employed to obtain the exact geometric quantities for the analysis. Especially, because NURBS is the most powerful and wide-spread method to represent general surfaces in the field of computer graphics and CAD(Computer Aided Design) industry, the direct computation of surface geometric quantities from the NURBS surface equation without approximation shows great potential for the integration between geometrically exact shell finite element and geometric modeling in the CAD systems. Some numerical examples are given to verify the performance and accuracy of the developed linkage framework. In additions, trimmed surfaces with some cutouts are considered for more practical applications.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.16 no.2
• /
• pp.21-27
• /
• 2008

An inverse method is introduced to construct the problem for the error analysis of the numerical solution of initial value problem. These problems constructed through this method have a known exact solution, even though analytical solutions are generally not obtainable. The process leading to the exact solution makes use of an initially available approximate numerical solution. A smooth interpolation of the approximate solution is forced to exactly satisfy the differential equation by analytically deriving a small forcing function to absorb all of the errors in the interpolated approximate solution. Using this special case exact solution, it is possible to investigate the relationship between global errors of a candidate numerical solution process and the associated tuning parameters for a given problem. Under the assumption that the original differential equation is well-posed with respect to the small perturbations, we thereby obtain valuable information about the optimal choice of the tuning parameters and the achievable accuracy of the numerical solution.