• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.225-245
• /
• 2008

T-splines are recently proposed mathematical tools for geometric modeling, which are generalizations of B-splines. Local refinement can be performed effectively using T-splines while it is not the case when B-splines or NURBS are used. Using T-splines, patches with unmatched boundaries can be combined easily without special techniques. In the present study, an analysis framework using T-splines is proposed. In this framework, T-splines are used both for description of geometries and for approximation of solution spaces. This analysis framework can be a basis of a CAD/CAE integrated approach. In this approach, CAD models are directly imported as the analysis models without additional finite element modeling. Some numerical examples are presented to illustrate the effectiveness of the current analysis framework.

• Structural Engineering and Mechanics
• /
• v.30 no.5
• /
• pp.513-534
• /
• 2008

This paper examines local vibrations in the stay cables of a cable-stayed bridge subjected to wind gusts. The wind loads, including the self-excited load and the buffeting load, are converted into time-domain values using the rational function approximation and the multidimensional autoregressive process, respectively. The global motion of the girder, which is generated by the wind gusts, is analyzed using the modal analysis method. The local vibration of stay cables is calculated using a model in which an inclined cable is subjected to time-varying displacement at one support under global vibration. This model can consider both forced vibration and parametric vibration. The response characteristics of the local vibrations in the stay cables under wind gusts are described using an existing cable-stayed bridge. The results of the numerical analysis show a significant difference between the combined parametric and forced vibrations and the forced vibration.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.88-89
• /
• 2011

풍하중이 작용하는 교량의 응답을 구하기 위하여 RFA(Rational Function Approximation)와 같은 시간 영역해석법이 널리 사용되고 있다. 교량 단면의 공기역학적 특성을 정의하는 플러터 계수는 주파수 영역에서 정의되기 때문에, 시간 영역해석을 위하여 inverse Fourier transform을 통해 얻어진 impulse response function을 이용한 중첩 적분법이 제안되었다. 시간 영역해석을 위해서는 플러터 계수에 상관성이 존재해야 함을 밝히고, 최적화 방법을 이용하여 시간 영역 해석을 위한 플러터 계수 산정법을 제안하고자 한다. B/D=20의 구형 단면에 적용하여 제안한 방법의 타당성을 검증하고자 한다.

• Structural Engineering and Mechanics
• /
• v.41 no.5
• /
• pp.617-632
• /
• 2012

An isogeometric approach is presented for static analysis of thin plate problems of various geometries. Non-Uniform Rational B-Splines (NURBS) basis function is applied for approximation of the thin plate deflection, as for description of the geometry. The governing equation based on Kirchhoff plate theory, is discretized using the standard Galerkin method. The essential boundary conditions are enforced by the Lagrange multiplier method. Several typical examples of thin plate and thin plate on elastic foundation are solved and compared with the theoretical solutions and other numerical methods. The numerical results show the robustness and efficiency of the proposed approach.

• Journal of Drive and Control
• /
• v.15 no.4
• /
• pp.67-73
• /
• 2018

The paper deals with an approach to time domain simulation for closed end at the downstream of pipe, hydraulic lines terminating into a tank and series lines with change of cross sectional area. Time domain simulation of a fluid power systems containing hydraulic lines is very complex and difficult if the transfer functions consist of hyperbolic Bessel functions which is the case for the distributed parameter dissipative model. In this paper, the magnitudes and phases of the complex transfer functions of hydraulic lines are calculated, and the MATLAB Toolbox is used to formulate a rational polynomial approximation for these transfer functions in the frequency domain. The approximated transfer functions are accurate over a designated frequency range, and used to analyze the time domain response. This approach is usefully to simulate fluid power systems with hydraulic lines without to approximate the frequency dependent viscous friction.

• Honam Mathematical Journal
• /
• v.45 no.3
• /
• pp.542-554
• /
• 2023

Given a real number α, the Lagrange number of α is the supremum of all real numbers L > 0 for which the inequality |α - p/q| < (Lq²)^-1 holds for infinitely many rational numbers p/q. All Lagrange numbers less than 3 can be arranged as a set {l_p/q : p/q ∈ ℚ ∩ [0, 1]} using the Farey index. The present paper considers a function C(α) devised from Sturmian words. We demonstrate that the function C(α) contains all information on Lagrange numbers less than 3. More precisely, we prove that for any real number α ∈ (0, 1], the value C(α) - C(0) is equal to the sum of all numbers 3 - l_p/q where the Farey index p/q is less than α.

• International Journal of CAD/CAM
• /
• v.4 no.1
• /
• pp.33-45
• /
• 2004

In this work, a rational procedure has been formulated for the selection of points approximating slice contours cut in LOM (Laminated Object manufacturing) with first order approximation. It is suggested that the number of points representing a slice contour can be 'minimised' or 'optmised' by equating the horizontal chordal deviation (HCD) to the user-defined surface form tolerance. It has been shown that such optimization leads to substantial reduction in slice height calculations and NC codes file size for cutting out the slices. Due to optimization, the number of contour points varies from layer to layer, so that points on successive layer contours have to be matched by four sided ruled surface patches and triangular patches. The technological problems associated with the cutting out of triangular patches have been addressed. A robust algorithm has been developed for the determination of slice height for optimum and arbitrary numbers of contour points with different strategies for error calculations. It has been shown that optimisation may even lead to detection and appropriate representation of elusive surface features. An index of optimisation has been defined and calculations of the same have been tabulated.

• Nuclear Engineering and Technology
• /
• v.52 no.9
• /
• pp.1907-1916
• /
• 2020

In this work, the depletion capability implemented in Monte Carlo code MCS is investigated to predict the isotopic compositions of spent nuclear fuel (SNF). By comparison of MCS calculation results to post irradiation examination (PIE) data obtained from one pressurized water reactor (PWR), the validation of this capability is conducted. The depletion analysis is performed with the ENDF/B-VII.1 library and a fuel assembly model. The transmutation equation is solved by the Chebyshev Rational Approximation Method (CRAM) with a depletion chain of 3820 isotopes. 18 actinides and 19 fission products are analyzed in 14 SNF samples. The effect of statistical uncertainties on the calculated number densities is discussed. On average, most of the actinides and fission products analyzed are predicted within ±6% of the experiment. MCS depletion results are also compared to other depletion codes based on publicly reported information in literature. The code-to-code analysis shows comparable accuracy. Overall, it is demonstrated that the depletion capability in MCS can be reliably applied in the prediction of SNF isotopic inventory.

• Korean Journal of Computational Design and Engineering
• /
• v.3 no.3
• /
• pp.183-191
• /
• 1998

There are many available algorithms based on the different approaches to solve the intersection problems between two curves. Among them, the implicitization method is frequently used since it computes precise solutions fast and is robust in lower degrees. However, once the degrees of curves to be intersected are higher than cubics, its computation time increases rapidly and the numerical stability gets worse. From this observation, it is natural to transform the original problem into a set of easier ones. Therefore, curves are subdivided appropriately depending on their geometric behavior and approximated by a set of rational quadratic Bezier cures. Then, the implicitization method is applied to compute the intersections between approximated ones. Since the solutions of the implicitization method are intersections between approximated curves, a numerical process such as Newton-Raphson iteration should be employed to find true intersection points. As the seeds of numerical process are close to a true solution through the mix-and-match process, the experimental results illustrates that the proposed algorithm is superior to other algorithms.

• Structural Engineering and Mechanics
• /
• v.45 no.6
• /
• pp.759-777
• /
• 2013

In this paper, a topology optimization method based on the element independent nodal density (EIND) is developed for continuum solids with multiple load cases and multiple constraints. The optimization problem is formulated ad minimizing the volume subject to displacement constraints. Nodal densities of the finite element mesh are used a the design variable. The nodal densities are interpolated into any point in the design domain by the Shepard interpolation scheme and the Heaviside function. Without using additional constraints (such ad the filtering technique), mesh-independent, checkerboard-free, distinct optimal topology can be obtained. Adopting the rational approximation for material properties (RAMP), the topology optimization procedure is implemented using a solid isotropic material with penalization (SIMP) method and a dual programming optimization algorithm. The computational efficiency is greatly improved by multithread parallel computing with OpenMP to run parallel programs for the shared-memory model of parallel computation. Finally, several examples are presented to demonstrate the effectiveness of the developed techniques.